Discussion :

[armand.favrot]

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 :



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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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?