1 pièce(s) jointe(s)
Portrait de phase en scilab
Bonjour,
Mon problème est de tracer un jolie portrait de phase d'un système dynamique. J'ai beaucoup cherché y compris avec d'autres logiciel (Python, R) mais je n'ai
pas vraiment de bonne solution. la meilleur que j'ai est avec scilab mais c'est loin d'être parfait, ça fait des gros patés :
Pièce jointe 349952
Sans compé que ma figure a été retourné mais ça c'est un autre problème : j'ai posté un autre message pour ce souci.
Voila le code que j'utilise :
Code:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
|
clear;
funcprot(0);
xdel(winsid());
function [out]=d(t,y) //Dilution différente
out(1)=0.1*(2-y(1))-(y(1)/(0.2+y(1)))*y(2)
out(2)=(y(1)/(0.2+y(1))-0.7)*y(2)
endfunction
scf(2)
t1=linspace(0,200,n_ref);
sol= ode([0.05;0.1], 0, t1, d);
plot2d4(sol(1,:),sol(2,:))
sol= ode([1.9;0.01], 0, t1, d);
plot2d4(sol(1,:),sol(2,:))
sol= ode([1.9;0.7], 0, t1, d);
plot2d4(sol(1,:),sol(2,:))
xlabel("s");
ylabel("x"); |
Il y en a qui se passent des flèches mais je ne trouve pas ça satisfaisant.. Quelqu'un a-t-il une solution?