Ode45

Fermé
ecohydrologie Messages postés 1 Date d'inscription lundi 24 mars 2008 Statut Membre Dernière intervention 3 avril 2008 - 3 avril 2008 à 14:13
Sacabouffe Messages postés 9427 Date d'inscription dimanche 19 août 2007 Statut Membre Dernière intervention 29 mai 2009 - 10 avril 2008 à 18:03
Bonjour,
jé un probléme pour simuler une equation différentiel non autonome en matlab en faite je veux utilisé ode45pour la simulation mais jé des probléme

merci davance
voici l equation


******************************************

function dydt = myode1(t,y)


mu1 = 1/12;
mu2 = 1/12;
mu3 = 1/12;
gamaG = 1/4;
gamaR = 1/4;
K = 130;
%f = interp1(ft,f,t); % Interpolate the data set (ft,f) at time t
%g = interp1(gt,g,t); % Interpolate the data set (gt,g) at time t

dy = zeros(3,1); % a column vector
dy(1) = g(t) * y(3) - (mu1 + gamaG) * y(1);
dy(2) = f(t)* (y(3) + y(2)) * (1 - (y(3) + y (2)) / K ) - (mu2 + gamaR) * y(2) ;
dy(3) = gamaG * y(1) + gamaR * y(2) - mu3 * y(3);

************************************************

1 réponse

Sacabouffe Messages postés 9427 Date d'inscription dimanche 19 août 2007 Statut Membre Dernière intervention 29 mai 2009 1 835
10 avril 2008 à 18:03
Salut
Il y a moult erreurs et imprécisions, je peux rien faire avec ça pour le moment.
Sur quel intervalle de temps tu veux résoudre ton équa diff?
Qu'est-ce que f et g?
Il fat pas les commenter si tu t'en sers ensuite. A moins que tu veuilles les utiliser comme elles sont définies précédemment. Sauf que là elles le sont pas...
Je pense qu'ensuite c'est dydt le nom, et non dy tout court. A moins que tu n'aies pas mis ton programme en entier et qu'il y ait quelque chose après.
A plus
0