Exemples d'équations différentielles

Cet article présente quelques exemples d'équations différentielles.

Les équations différentielles sont utilisées pour résoudre des problèmes en physique, en ingénierie et dans plusieurs autres sciences.

Équations du premier ordre

Une équation différentielle linéaire à variables séparables

Les équations différentielles les plus simples sont les équations linéaires homogènes du premier ordre. Par exemple, équation différentielle linéaire à variables séparables, de la forme :

${\frac {\mathrm {d} y}{\mathrm {d} t}}+f(t)y=0$

où f est une fonction connue admettant des primitives. Une façon directe de la résoudre est de considérer, pour f(t) non nul, la forme avec les variables séparées :

${\frac {\mathrm {d} y}{y}}=-f(t)\mathrm {d} t.$

Par intégration, il vient alors :

$\ln(y)+C=\int {\frac {\mathrm {d} y}{y}}=-\int f(t)\mathrm {d} t\Longrightarrow y=A\exp \left(-\int f(t)\mathrm {d} t\right)\,$

où A = e^C est une constante arbitraire.

Une équation différentielle linéaire à variables non séparables

Le cas plus général inclut des formes où les variables ne sont pas séparables, comme :

${\frac {\mathrm {d} y}{\mathrm {d} x}}+p(x)y=q(x).$

où p et q sont des fonctions connues admettant des primitives. Une façon directe de la résoudre est de la réécrire avec un facteur intégrant :

$\mu (x)=\exp \left(\int _{x_{0}}^{x}p(t)\mathrm {d} t\right)$

Ce qui donne, en remultipliant par le facteur μ

$\mu (x){\frac {\mathrm {d} y}{\mathrm {d} x}}+\mu (x)p(x)y=\mu (x){\frac {\mathrm {d} y}{\mathrm {d} x}}+{\frac {\mathrm {d} \mu }{\mathrm {d} x}}y={\frac {\mathrm {d} (\mu (x)y)}{\mathrm {d} x}}=\mu (x)q(x).$

dont on déduit la forme générale de la solution :

$y={\frac {\int \mu (x)q(x)\,\mathrm {d} x+C}{\mu (x)}}.$

Équations du second ordre

Une oscillation simple non amortie

Les mouvements périodiques dont on néglige les effets de frottement qui vont le ralentir (comme l'allongement du ressort x(t) à un temps t) peuvent être modélisés par l'équation différentielle suivante :

${\frac {d^{2}x}{dt^{2}}}=-\omega ^{2}x\,$ où ω est un réel positif

Dont les solutions sont :

$x(t)=A\cos(\omega t)+B\sin(\omega t)\,$

Pour déterminer les constantes A et B, il faut utiliser les conditions initiales qui permettent de décrire l'état du système à un instant donné (correspondant en général à t = 0).

Par exemple si nous supposons qu'à l'instant t = 0, l'extension du ressort est d'une unité de longueur (x = 1), et la masse est immobile (dx/dt = 0). Nous pouvons en déduire

$x(0)=A\cos 0+B\sin 0=A=1\,$,

d'où l'on déduit A = 1.

$x'(0)=-\omega A\sin 0+\omega B\cos 0=\omega B=0\,$,

et donc B = 0.

En conséquence $x(t)=\cos(\omega t)\,$ est solution de l'équation différentielle étudiée.

Plus souvent en physique pour les oscillations simples non amorties, on utilise une solution de la forme :

$x(t)=A\cos(\omega t+\phi ),$

avec A étant l'amplitude et ϕ la phase.

Pour l'exemple cité on procède :

$x(0)=A\cos(\phi )=1$

$x'(0)=-A\omega \sin(\phi )=0$

Ce qui donne ϕ = 0 et par conséquent A = 1.

D'où le résultat $x(t)=\cos(\omega t)$

La solution la plus générale en fonction de conditions initiales quelconques x₀ et ${\dot {x}}_{0}$ est donnée par l'équation :

$x(t)=x_{0}\cos {(\omega t)}+{\dfrac {{\dot {x}}_{0}}{\omega }}\sin {(\omega t)}$

Prise en compte de l'amortissement

Oscillation amortie

Le modèle précédent négligeait les forces de frottement. De ce fait l'oscillation libre pouvait durer indéfiniment, ce qui n'est jamais observé en réalité.

Les frottements sont en général une force proportionnelle à la vitesse (dx/dt) et opposée au mouvement. En rajoutant ce terme notre équation différentielle devient :

${\frac {d^{2}x}{dt^{2}}}=-c{\frac {dx}{dt}}-\omega ^{2}x\,$ où c > 0 est le coefficient de frottement.

Ceci est une équation différentielle linéaire à coefficients constants, homogène et du second ordre, que nous pouvons résoudre.

En cherchant une solution de la forme particulière A e_kt, nous constatons que k doit vérifier l'équation caractéristique suivante :

$k^{2}+ck+\omega ^{2}=0\,$.

On revient ainsi à l'étude d'une équation du second degré. Si c < 2ω nous avons deux racines complexes conjuguées (de la forme a ± i b), et la solution (avec les conditions initiales identiques au cas précédent) a la forme suivante :

$x(t)={\rm {e}}^{-{\frac {c}{2}}t}\left(\cos(t{\sqrt {4\omega ^{2}-c^{2}}})+{\frac {-c}{2{\sqrt {4\omega ^{2}-c^{2}}}}}\sin(t{\sqrt {4\omega ^{2}-c^{2}}})\right)\,$

Nous pouvons démontrer que a < 0.

Le système étudié (le pendule pesant dans le référentiel terrestre supposé galiléen) est le siège d'oscillations libres amorties.

Ce sont les positions du centre d'inertie de la masse, en fonction du temps, avec x = 0 correspondant à une position d'équilibre.

NB : la courbe présente une allure proche d'un régime critique : la position d'équilibre est à peine franchie, et on ne compte guère plus d'une pseudo-période d'oscillations.

Articles connexes

• Intégration
• Analyse vectorielle