Discussion :

[Mat32]

Hello,

J'aimerais résoudre un système d'équations différentielles, cependant il y a une équation algébrique. Est-il possible de l'inclure ?

Voici le système d'équations :

Code :

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

1
2
3
4
5
6
 
dy(1) = F;
dy(2) = -p.k1*y(2)*y(3) + p.k2*y(4)*y(5) - y(2)/y(1)*dy(1);
dy(3) = -p.k1*y(2)*y(3) + p.k2*y(4)*y(5) - y(3)/y(1)*dy(1);
dy(4) = +p.k1*y(2)*y(3) - p.k2*y(4)*y(5) - y(4)/y(1)*dy(1);
dy(5) = +p.k1*y(2)*y(3) - p.k2*y(4)*y(5) - y(5)/y(1)*dy(1);

L’équation algébrique à ajouter est la suivante, cependant le code qui suit ne fonctionne pas :

Code :

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

 y(6) = y(2)*p.k3

Merci

[FR119492]

Salut!
Quelques questions et remarques pour clarifier le problème:

1. Tes dy(1) à dy(5) semblent être des dérivées, mais par rapport à quoi? En d'autres termes, quelle est la variable indépendante?
2. Les grandeurs F et p sont-elles des variables ou des constantes?
3. La variable y(6) n'intervient pas dans ton système différentiel. Elle doit donc être calculée à part, à partir de y(2), à chaque pas d'intégration.

Jean-Marc Blanc

[Mat32]

Quelle est la variable indépendante?

C'est le temps : t

Les grandeurs F et p sont-elles des variables ou des constantes?

La grandeur F est variable mais pas les grandeur p, mais celles-ci sont gérés par une autre partie du code

La variable y(6) n'intervient pas dans ton système différentiel. Elle doit donc être calculée à part, à partir de y(2), à chaque pas d'intégration.

C'est exactement ma question, comment calculer y(6) à chaque pas car il faudrait, dans l'idéal, utiliser la variable y(6) dans le système d'équations différentielles.

[Invité]

C'est exactement ma question, comment calculer y(6) à chaque pas car il faudrait, dans l'idéal, utiliser la variable y(6) dans le système d'équations différentielles.

Et Jean Marc y a répondu
y(6) n'intervenant pas dans ton système différentiel, tu le résous sans, puis tu le calcules à partir des valeurs de y(2) obtenues.

[Mat32]

y(6) n'intervenant pas dans ton système différentiel, tu le résous sans

y(6) n'intervient pas car je ne sais pas comment faire mais cette variable devrait intervenir...

Merci

[Invité]

y(6) n'intervient pas car je ne sais pas comment faire mais cette variable devrait intervenir...

Dans ce cas, comment intervient-elle ?
Quel est ton système d'équations ?

[Mat32]

Voici le système :

Code :

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

1
2
3
4
5
6
7
dy(1) = F;
dy(2) = -p.k1*y(2)*y(3) + p.k2*y(4)*y(5) - y(2)/y(1)*dy(1);
dy(3) = -p.k1*y(2)*y(3) + p.k2*y(4)*y(5) - y(3)/y(1)*dy(1);
dy(4) = +p.k1*y(2)*y(3) - p.k2*y(4)*y(5) - y(4)/y(1)*dy(1);
dy(5) = +p.k1*y(2)*y(3) - p.k2*y(4)*y(5) - y(5)/y(1)*dy(1) + p.k4*y(6);
 
 y(6) = p.k3*y(2)

Comment pourrais-je calculer y(6) à chaque pas ?
Comment pourrais-je résoudre le système d'équations sans substituer y(6) par p.k3*y(2) dans le code ?

Merci

[Invité]

Tu possèdes y(6) = p.k3*y(2)... il n'est pas difficile de d'obtenir dy(6) à partir de là

[Mat32]

Mais il est impossible de calculer y(6) (ou une autre valeur) à chaque pas ?

[membreComplexe12]

Mais il est impossible de calculer y(6) (ou une autre valeur) à chaque pas ?

si il est possibe mais il faut ecrire tous le systeme sous forme différentielle et il faut donc que tu remplace dans ton systeme ta ligne y(6) par une ligne dy(6) qui sera la dérivée par rapport au temps de y(6)

Ensuite ODE45 va te résoudre ce systeme diffférentielle et te donnerai l'évolution en fonction du temps de y(1) y(2) y(3)...y(6)

[FR119492]

Bonjour à tous!
Cette discussion est une magnifique illustration du débat en cours sur ce site:

Quelque chose ne va vraiment pas avec les développeurs "modernes", un développeur à "l'ancienne" critique la multiplication des bibliothèques

Mat32 essaie d'utiliser une fonction toute faite, à savoir ODE45; mais connait-il l'algorithme utilisé? On ne lui a probablement jamais enseigné comment programmer la méthode de Runge-Kutta. Alors que, s'il savait le faire, et s'il utilisait un langage de bas niveau, comme le C ou le Fortran, il n'aurait aucune peine à calculer y(6) à chaque pas. Je pense que William Edwards a raison et que, même si Matlab est un outil remarquable, il n'est pas fait pour les débutants.
Jean-Marc Blanc