Voir la source de Résolution d'équations différentielles simples/Exemples concrets

{{Chapitre
  | niveau    = 14
  | idfaculté = physique
  | numéro    = 3
  | précédent = [[../Équations sans second membre/]]
  | suivant   = 
}}


== Équation du premier ordre avec second membre constant  ==
On prend l'exemple d’un parachutiste tombant à la verticale, qui, à un instant ''t = 0'' ouvre son parachute alors qu’il avait une vitesse ''v(0)''. On suppose, pour simplifier le problème, que son parachute lui fait subir une force de frottement opposée à sa vitesse. On peut alors montrer en mécanique que l'équation différentielle que vérifie la vitesse ''v'' est :
<math>\frac{dv}{dt} + a v(t)= g</math> où ''g'' (accélération de la pesanteur) et ''a'' (coefficient de frottement) sont des constantes.
On cherche d’abord une ''solution particulière'', notée ''v<sub>p</sub>''. Pour cela, on peut remarquer qu'en la supposant constante, le terme faisant intervenir sa dérivée est forcément nul (la dérivée d’une constante est nulle). Ainsi, il resterait ''a v<sub>p</sub> = g'' ou encore ''v<sub>p</sub> = g/a''. On peut vérifier, en introduisant cette relation dans l'équation différentielle, que c’est effectivement une solution particulière : on la retient pour la suite.

Il nous faut maintenant calculer les solutions ''v₀'' de l'équation sans second membre qui s'écrit :
<math>\frac{dv_0}{dt} + a v_0(t) = 0</math>
Le paragraphe sur l'équation différentielle du premier ordre sans second membre nous donne directement la solution à cette équation : 
<math>v_0(t)=C\exp(-at)</math>
où ''C'' est une constante qu’il nous faudra déterminer.
[[Fichier:Exponential limit.PNG|thumb|right|250px|Évolution de la vitesse du parachutiste en fonction du temps.]]
On sait alors que toute solution générale de l'équation de départ est la somme de la solution particulière et de la solution de l'équation sans second membre :
<math>v(t) = v_p + v_0(t) = \frac{g}{a} + C\exp(-at)</math>
Il nous reste à déterminer la constante C. Pour cela, on évalue l’expression précédente en ''t = 0'' :
<math>v(0) = \frac{g}{a} + C \exp(0) = \frac{g}{a} + C</math>
On a donc ''C = v(0) - g/a''. D'où :
<math>v(t) = \frac{g}{a} + \left( v(0)-\frac{g}{a}\right) \exp(-at)</math>
Finalement, on s'aperçoit que la vitesse évolue de façon exponentielle depuis la vitesse initiale ''v(0)'' à une vitesse finale ''g/a''. Donc au bout d’un moment, le parachutiste va être maintenu à une vitesse maximale, appelée ''vitesse limite''.

== Équation du premier ordre avec second membre oscillant  ==
On considère un circuit électrique contenant une résistance ''R'', un condensateur ''C'' et un « générateur basse fréquence » (amplitude ''E'' et fréquence <math>\omega </math>) associés en série. On obtient alors en électricité l'équation différentielle vérifiée par la charge ''q'' du condensateur :
<math>\frac{dq}{dt} + \frac{1}{RC} q = E \cos(\omega t)</math>

Premièrement, on cherche la solution particulière ''q<sub>p</sub>'' sous la forme <math>q_p = A \cos(\omega t)+B\sin(\omega t)</math>. Pour déterminer ''A'' et ''B'', on insère l’expression de ''q<sub>p</sub>'' dans l'équation différentielle, ce qui donne :
<math>\left( -A\omega+\frac{B}{RC}\right) \sin(\omega t) + \left( B\omega +\frac{A}{RC}\right) \cos(\omega t) = \frac{E}{R} \cos(\omega t)</math>
Or cela est vrai quel que soit ''t''. On choisit donc successivement deux valeurs différentes de ''t''. En l’occurrence, les valeurs ''0'' et <math>\pi/2\omega</math> nous donnent respectivement :
<math>B\omega+\frac{A}{RC}=\frac{E}{R}</math> et <math>-A\omega+\frac{B}{RC}=0</math>
On en déduit, après quelques calculs, que
<math>A=\frac{B}{RC\omega}=\frac{CE}{1+R^2C^2\omega^2}</math>
Nous avons donc déterminé la solution particulière.
[[Fichier:Exponential plus cosine.PNG|thumb|right|250px|Évolution de la charge du condensateur en fonction du temps.]]
Deuxièmement on résout l'équation sans second membre qui s'écrit de la manière suivante :
<math>\frac{dq_0}{dt} + \frac{1}{RC} q_0 = 0</math>
Le paragraphe sur les équations différentielles du premier ordre sans second membre nous donne directement la solution :
<math>q_0(t)=\alpha \exp\left(-\frac{t}{RC}\right)</math>
On obtient finalement la solution de l'équation de départ en sommant ''q₀'' et ''q<sub>p</sub>'' :
<math>q(t)=\frac{CE}{1+R^2C^2\omega^2}\left( \cos(\omega t) + RC\omega \sin(\omega t) + \alpha\exp\left(-\frac{t}{RC}\right)\right)</math>
Il ne reste plus qu’à déterminer la constante <math>\alpha</math>. Pour cela, on suppose connue la charge initiale ''q(0)''. Cela donne l’expression finale de ''q(t)'' :
<math>q(t)=\frac{CE}{1+R^2C^2\omega^2}\left( \cos(\omega t) + RC\omega \sin(\omega t)-1\right) + q(0)\exp\left(-\frac{t}{RC}\right)</math>
Ainsi on remarque que la charge du condensateur va suivre, au bout d’un temps de l’ordre de ''R C'', l'oscillation de la tension, avec toutefois un certain décalage.

== Équation du second ordre avec second membre nul  ==
On considère un atome formé d’un noyau supposé ponctuel et fixe, et d’un électron gravitant autour à une distance ''x(t)''. Cet électron est soumis à l'attraction du noyau (force de rappel élastique) et à un ralentissement (car il émet de la lumière en bougeant). On peut alors montrer que l'équation différentielle que vérifie ''x(t)'' est :
<math>\frac{d^2x}{dt^2}+a\frac{dx}{dt}+\frac{b^2}{4}x = 0</math>
où ''a'' et ''b'' sont des constantes positives caractéristiques du ralentissement et de l'attraction respectivement.

Pour résoudre cette équation, on écrit l'''équation caractéristique'' : <math>r^2 + a r + b^2/4 = 0</math>. Ses solutions sont :
<math>r_1=\frac{a}{2}\left( -1+\sqrt{1-b^2/a^2}\right) </math> et <math>r_2=\frac{a}{2}\left( -1-\sqrt{1-b^2/a^2}\right) </math>
Les solutions de l'équation différentielle sont alors :
<math>x(t)=C_1\exp(r_1t)+C_2\exp(r_2t)</math>
[[Fichier:Oscillation amortie.png|thumb|right|Évolution de la position de l'électron dans le régime d'oscillation amortie.]]
On ne va pas chercher à calculer les constantes ''C₁'' et ''C₂'' car le calcul est fastidieux et peu intéressant. À ce stade, on ne sait pas si ''r₁'' et ''r₂'' sont des nombres réels ou complexes. On distingue donc différents cas :
* Si <math>b\le a</math>, alors ''r₁'' et ''r₂'' sont réels, et d’après leur expression, ils sont tous deux négatifs. Chaque exponentielle tend donc vers 0. Cette situation correspond au cas où l'électron se rapproche indéfiniment du noyau de façon monotone.
* Si ''b > a'', alors ''r₁'' et ''r₂'' sont complexes. En séparant leur partie réelle et leur partie imaginaire, on obtient :
<math>x(t)=\exp\left(-\frac{at}{2}\right) \left( C_1\exp\left(i\sqrt{b^2-a^2}\frac{t}{2}\right) +C_2\exp\left(-i\sqrt{b^2-a^2}\frac{t}{2}\right) \right)</math>
:On se place dans le cas où l'électron a, au départ, une vitesse nulle. On peut alors montrer en dérivant l'équation précédente que ''C₁ = - C₂''. D'où la simplification :
<math>x(t)=2 C_1 \exp\left(-\frac{at}{2}\right) \cos\left(i\sqrt{b^2-a^2}\frac{t}{2}\right) </math>
:Le terme avec le cosinus indique que l'électron effectue des oscillations. Le terme avec l'exponentielle impose au mouvement de tendre vers 0. Le mouvement global de l'électron est donc une oscillation amortie.
On retrouve ce type de mouvement dans de nombreux problèmes de physique. Il est important de se souvenir qu’il existe deux régimes : un amortissement simple, ou un amortissement avec oscillation.

{{Bas de page
  | idfaculté = physique
  | précédent = [[../Équations sans second membre/]]
}}