Ode45
ecohydrologie
Messages postés
1
Statut
Membre
-
Sacabouffe Messages postés 10427 Statut Membre -
Sacabouffe Messages postés 10427 Statut Membre -
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);
************************************************
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
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
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
