Discussion :

[Mat32]

Hello,

Quelqu'un pourrait-il me renseigner sur la fonction ode45 ?
La réponse que me donne ode45 pour le système d'équations différentielles suivant est différente selon l'ordre des équations dy(...) = ... :

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
function dy = dABPQ(t,y,p)
dy = zeros(5,1);
 
if t>p.tdos;
   F=0;
else
   F=p.F;
end
 
dy(1) = F; 
dy(2) = -p.k1*y(2)*y(3) -p.k2*y(2)*y(3) -y(2)/y(1)*dy(1)                  ;
dy(3) = -p.k1*y(2)*y(3) -p.k2*y(2)*y(3) -y(3)/y(1)*dy(1) +p.Bini/y(1)*F   ;
dy(4) =  p.k1*y(2)*y(3)                 -y(4)/y(1)*dy(1)                  ;
dy(5) =                  p.k2*y(2)*y(3) -y(5)/y(1)*dy(1)                  ;
 
 
end

la réponse si dy(1) = F est déplacé à la fin, et donc devient dy(5), est différente

Pour le contexte il s'agit de deux réactions parallèles et dy(1), y(1) sont respectivement la variation de volume et le volume.

[salseropom]

salut, non

tu peux très bien écrire dy(1) = F et dy(5) = qqch ou bien dy(5) = qqch puis dy(1) = F. L'ordre n'a AUCUNE importance. Ce qui est TRES important, ce sont les indices (1 et 5 dans l'exemple)

[Mat32]

Alors essaie, je te promet que l'ordre change la réponse.

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
function dy = dABPQ(t,y,p)
dy = zeros(5,1);
 
if t > p.tdos;
   F = 0;
else
   F = p.F;
end
 
 
dy(2) = -p.k1*y(2)*y(3)^p.nb1 -p.k2*y(2)*y(3)^p.nb2 -y(2)/y(1)*dy(1)                  ;
dy(3) = -p.k1*y(2)*y(3)^p.nb1 -p.k2*y(2)*y(3)^p.nb2 -y(3)/y(1)*dy(1) +p.Bini/y(1)*F   ;
dy(4) =  p.k1*y(2)*y(3)^p.nb1                       -y(4)/y(1)*dy(1)                  ;
dy(5) =                        p.k2*y(2)*y(3)^p.nb2 -y(5)/y(1)*dy(1)                  ;
dy(1) = F     ;                        
 
end

Ce code ne donne pas la même réponse.

[salseropom]

Je mientiens ce que je dis

maintenant, dans ton code c'est normal :

Code :

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

1
2
3
 
dy(5) =                        p.k2*y(2)*y(3)^p.nb2 -y(5)/y(1)*dy(1)                  ;
dy(1) = F

ici, dans le dy(5), ton dy(1) vaut 0

ce qui explique pourquoi tu n'as pas la même réponse