Fonctions > Résolution et optimisation > Solveurs d'équations différentielles > Exemple : Oscillation harmonique amortie dans des circuits
Exemple : Oscillation harmonique amortie dans des circuits
Etant donné un circuit électrique simple contenant une résistance R, une bobine d'induction L et un condensateur C.
Utilisez des équations différentielles pour modéliser la charge Q sur le condensateur C, puis utilisez les fonctions de solveur d'EDO pour trouver d'autres solutions approximatives. Enfin, comparez les résultats à la solution exacte de Q.
Utilisation d'équations différentielles
1. Ecrivez une équation différentielle pour la tension des trois composants dont le total doit être égal à zéro.
Cliquez pour copier cette expression
2. Ecrivez l'équation différentielle correspondant au changement instantané de charge Q sur le condensateur C :
Cliquez pour copier cette expression
3. Définissez la charge Q au moment zéro :
Cliquez pour copier cette expression
4. Partez des valeurs suivantes pour les paramètres d'entrée :
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
5. Transformez l'équation différentielle sous la forme standard :
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
6. Résolvez l'équation pour tous les cas de a et de b :
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
7. Tracez les deux solutions possibles pour Q :
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
La charge Q1 tombe rapidement à zéro lorsque a=b, et oscille pendant un bon moment avant de passer à zéro quand a et b ne sont pas identiques.
Avant d'utiliser les solveurs d'EDO
Définissez les paramètres à transmettre aux fonctions du solveur d'EDO :
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Les solveurs d'EDO sont de deux types : ceux des systèmes raides et ceux des systèmes non raides. Un système d'EDO écrit sous forme matricielle, par exemple y' = Ax, est qualifié de raide si la matrice A est quasi singulière. Sinon, le système est non raide.
La distinction entre systèmes raides et non raides peut être associée aux échelles de dynamique inhérentes de la matrice A telles que définies par ses valeurs propres. Une matrice dotée de valeurs propres disparates (ou de valeurs allant du trop petit au trop grand) constitue généralement un système raide.
Selon les paramètres système spécifiés, ce simple exemple d'oscillation harmonique amortie représente un système non raide.
Les fonctions rkfixed, Rkadapt et Bulstoer utilisent la valeur TOL du document.
Les fonctions Adams, AdamsBDF, BDF et Radau peuvent prendre une valeur spécifique de TOL ou utiliser leur propre valeur TOL par défaut de 10-5.
Utilisation de solveurs d'EDO pour systèmes non raides
Utilisez les solveurs d'EDO pour systèmes non raides pour trouver les solutions approximatives, puis comparez les résultats à la solution exacte de Q.
1. Adams
Cliquez pour copier cette expression
Cliquez pour copier cette expression
2. rkfixed
Cliquez pour copier cette expression
Cliquez pour copier cette expression
3. Rkadapt
Cliquez pour copier cette expression
Cliquez pour copier cette expression
4. Bulstoer
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Utilisation de solveurs d'EDO pour systèmes raides
Utilisez les solveurs d'EDO pour systèmes raides pour trouver les solutions approximatives, puis comparez les résultats à la solution exacte de Q.
1. BDF
Cliquez pour copier cette expression
Cliquez pour copier cette expression
2. Radau
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Utilisation du solveur hybride d'EDO
Utilisez le solveur hybride d'EDO AdamsBDF pour déterminer s'il s'agit d'un système raide ou non raide, puis appelez Adams ou BDF en conséquence pour trouver les solutions approximatives et comparer les résultats à la solution exacte de Q.
1. AdamsBDF
Cliquez pour copier cette expression
Cliquez pour copier cette expression
Conclusion
Dans les solveurs d'EDO pour systèmes non raides, les fonctions Adams et rkfixed renvoient les solutions avec, respectivement, l'erreur la plus petite et l'erreur la plus grande.
Dans les solveurs d'EDO pour systèmes raides, les fonctions Radau et BDF renvoient les solutions avec, respectivement, l'erreur la plus petite et la plus grande des deux.
La fonction hybride AdamsBDF renvoie une valeur qui est la plus petite des deux valeurs renvoyées par Adams ou BDF.
En règle générale, les fonctions Radau et rkfixed renvoient les solutions avec, respectivement, l'erreur la plus petite et l'erreur la plus grande.
Comparaison des résultats
1. Tracez les solutions renvoyées par les solveurs d'EDO pour les systèmes non raides avec les erreurs la plus petite (Adams, G0nss) et la plus grande (rkfixed, G1nss) :
Cliquez pour copier cette expression
Ces deux solutions présentent des différences évidentes par rapport à Q.
2. Tracez les solutions renvoyées par les serveurs d'EDO pour les systèmes raides avec les erreurs la plus petite (Radau, G1ss) et la plus grande (BDF, G0ss) :
Cliquez pour copier cette expression
Les deux solutions semblent identiques par rapport à Q.
3. Tracez les solutions qui ont donné les erreurs la plus petite (Radau, G1ss) et la plus grande (rkfixed, G1nss), renvoyées par les solveurs d'EDO pour les systèmes raides et non raides :
Cliquez pour copier cette expression
La fonction Radau, solveur de système raide, renvoie la meilleure approximation à la solution de Q.
Est-ce que cela a été utile ?