Discussion :

[touttty15]

Bonjour,
Je suis débutante avec MATLAB et je voudrais bien que vous m'aidez.
Je travaille sur la modélisation cinétique. J'ai un système d'équations différentielles de premier ordre avec des paramètres cinétiques km1 et km2 définie dans le fichier alpham.m suivant:

Code :

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

1
2
3
4
5
6
7
8
function dym= alpham (t,ym);
%réponse du modèle
dym(1)= -km1.*ym(1).*ym(2)-km2.*ym(1).*ym(3);
dym(2)=-km1.*ym(1).*ym(2);
dym(3)=km1.*ym(1).*ym(2)-km2.*ym(1).*ym(3);
dym(4)=km2.*ym(1).*ym(3);
dym=dym';
end

j'ai crée le script qui permet de résoudre ce système dans le fichier simalpha.m suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
%résolution du modèle alpha
grid
title ('modèle alpha')
xlabel ('temps(s)')
ylabel ('concentration')
hold on
horizon = [0:100:10000];
ym0 = [1  1  0  0 ];
[t,ym]=ode15s ('alpham' ,horizon, ym0,km1,km2)
plot (t,ym)
hold off

Ma question est comment puis-je faire un programme principal dont les arguments d'entrées sont les valeurs des paramètres km1 et km2 et qui fait appel à alpha.m et simalpha.m? (cad jun programme qui me permet de résoudre ce système à n'importe quelles valeurs de km1 et km2 et qui m'affiche le résultat ym en fonction du temps)

merci

[Invité]

Bonjour,

Une solution: les fonction imbriquées