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
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