Passer de ode45 MATLAB à ode SCILAB
Bonjour à tous.
Je suis actuellement en train de modéliser l'évolution de la pression, de la température et de la masse d'air dans le cylindre d'un moteur à combustion interne. J'ai déjà effectué ce travail sous MATLAB mais je dois à présent le faire sous Scilab (faute de licence). Seulement une chose me pose problème, la syntaxe de l'ode. L'équation différentielle à résoudre est du style : M(t,y)*dy/dt= f(t,y).
Sous Matlab j'utilisais les lignes suivantes :
Code:
1 2 3 4
|
tspan =[0 : 4*pi/719: 4*pi];
options =odeset('Mass',@Masse,'RelTol',1e-6);
[t,y]=ode45(@Force1,tspan,[Padm Tadm masse_ini],options); |
Mon problème se trouve au niveau de la fonction Masse (M(t,y)), je n'arrive pas à trouver comment la mettre dans l'ode sous Scilab.
J'espère avoir été clair dans mes explications :)
Merci
Edit : J'ai oublié de préciser que M est une matrice 3*3 et il y a donc trois équations à résoudre dans la fonction Force1.