jeudi 12 mars 2020

Using a condition to plot

I solve diff system with euler's method. Before solve system, i find a constant:

for a2 = 0.1:0.1:1.5
    for b1 = 0.1:0.1:1.5
        for b2 = 0.1:0.1:1.5
            for T1 = 0:0.1:1
                for T2 = 0:0.1:1
                    %euler's method
                    for i = 1:n - h

Then i use condition for system's order(or plot order?) and write const to txt file

if ((x1(end) - x1(1) <= 10^3))&&((mean(x1(end - 9:end))) < eps)
                        fprintf(fid, '%2.5f\t', a2);

In other script, for plot read txt file

A = load('param3.txt');
a2 = A(169,1);
b1 = A(169,2);
b2 = A(169,3);
T1 = A(169,4);
T2 = A(169,5);

Solve the diff system again and draw the plot:

plot(t,x1, t,x2,  t, a1); hold on
yline(x_w_st1, '--');
hold off

My plot with order x10^203

My plot with order x10^203

But i get wrong plot. Order is x10^203 and i need no more than x10^3. So, my condition if ((x1(end) - x1(1) <= 10^3)) doesn't work. What condition i need for order no more than x10^3?

P.S. I try set less step for system, this doesn't work too.

Aucun commentaire:

Enregistrer un commentaire