Équation différentielle linéaire/Exemples et intérêt

Introduction

Nous avons développé un arsenal mathématique pour résoudre ces équations, espérons que c’est utile ! Fort heureusement, ça l'est. En effet, pour beaucoup, la physique regorge d'équations différentielles, et on peut parfois les approcher par des équations linéaires.

On présente ici quelques exemples et rappelons les avantages de ces méthodes.

Rappels utiles

Ces quelques résultats nous seront utiles dans les exemples, bien que probablement connus du lecteur. Tout d’abord sur l’existence de racines :

En particulier, le cas des polynômes d'ordre deux est simplissime : Modèle:Théorème

Rappelons enfin une propriété immédiate pour ces polynômes :

Le fait d’utiliser la notation complexe, pour un problème qui ne traite que de quantité physiques réelles, permet de restreindre les solutions complexes. Démontrons un résultat utile :

Enfin, un théorème d'équivalence, parfois utile :

Exemple : décroissance radioactive

L'équation régissant la décroissance radioactive d'un ensemble de N particules est :

$\dot{N}=-\lambda N$

D'après le chapitre sur les équations différentielles ordinaires linéaires d'ordre un, la solution est :

$N\left(t\right)=N_0{e}^{-\lambda t}$

Avec N₀ le nombre de particules à l'instant t = 0.

Exemple : ressort libre

Soit une masse m, assimilée à un point, astreinte à se déplacer selon un axe x. Elle est retenue par un ressort de raideur k, de longueur à vide nulle. Mettons cela en équation à partir des lois de Newton :

$m\ddot{x}=-kx$

Réécrivons-la sous une forme similaire à celle étudiée dans le chapitre précédent :

$\ddot{x}=-\frac{k}{m}x$

Remarquons qu’il s'agit d'une équation homogène. Le formalisme introduit au chapitre 5 simplifiera l'étude de ce cas, mais on peut dors et déjà résoudre complètement cette équation avec les outils développés jusqu'ici.

Méthode 1 « astucieuse »

Supposons que la solution est une simple exponentielle : $x(t)=e^{\beta t}$. Alors :

$\ddot{x}(t)={\beta }^{2}x(t)$

Ce qui est solution si $\beta =\pm i\sqrt{\frac{k}{m}}=\pm i{\omega }_0$. Les deux solutions qu'on en tire étant linéairement indépendantes, toute solution à l'équation différentielle est de la forme :

$x(t)=A{e}^{i{\omega }_{0}t}+B{e}^{-i{\omega }_{0}t}$

Pour des raisons physiques, la position doit être un nombre réel. On a donc, en fin de compte, une fonction sinusoïdale :

$x(t)=𝒜\cos ({\omega }_{0}t+\varphi )$

Méthode 2 « matricielle »

Posons, toujours comme dans le chapitre précédent, le vecteur $\left(\begin{array}{c}\dot{x}\\ x\end{array}\right)$

On a alors l'équation différentielle :

$\left(\begin{array}{c}\ddot{x}\\ \dot{x}\end{array}\right)=\left(\begin{array}{cc}0& -\frac{k}{m}\\ 1& 0\end{array}\right)\cdot \left(\begin{array}{c}\dot{x}\\ x\end{array}\right)$

Les valeurs propres de la matrice carrée sont :

• ${\lambda }_1=+i\sqrt{\frac{k}{m}}=+i{\omega }_0$ ;
• ${\lambda }_2=-i\sqrt{\frac{k}{m}}=-i{\omega }_0$.

La solution finale est :

$\left(\begin{array}{c}\dot{x}\\ x\end{array}\right)={\mu }_1{e}^{i{\omega }_{0}t}{\overrightarrow{\mathrm{\Lambda }}}_1+{\mu }_2{e}^{-i{\omega }_{0}t}{\overrightarrow{\mathrm{\Lambda }}}_2$

Le même argument physique que dans la première méthode impose des solutions réelles :

$\left(\begin{array}{c}\dot{x}\\ x\end{array}\right)=𝒜\left(\begin{array}{c}{\omega }_0\sin ({\omega }_{0}t+\varphi )\\ \cos ({\omega }_{0}t+\varphi )\end{array}\right)$

Modèle:Bas de page