Voir la source de Équation différentielle linéaire d'ordre deux

Une '''équation différentielle linéaire d'ordre 2''' est une [[équation différentielle]] de la forme :
:<math>a\,y''+b\,y'+c\,y=d</math>
où <math>y</math> est une [[fonction (mathématiques)|fonction]] inconnue (d'une [[variable (mathématiques)|variable]]) et <math>a</math>, <math>b</math>, <math>c</math> et <math>d</math> quatre fonctions données (de la même variable). La variable et les fonctions peuvent être [[Fonction réelle d'une variable réelle|réelles]] ou [[Fonction d'une variable complexe différentiable au sens réel|complexes]].

Ces équations différentielles ne peuvent pas toutes être résolues explicitement mais plusieurs méthodes existent pour résoudre celles qui peuvent l'être, ou à défaut pour en étudier les solutions qualitativement. Les plus fréquentes et les plus simples à résoudre sont les équations à coefficients constants (où <math>a</math>, <math>b</math> et <math>c</math> sont des constantes).

La [[linéarité]] de l'équation a pour conséquence qu'il est possible d'appliquer des procédés de superposition de solutions, et d'exploiter des résultats d'[[algèbre linéaire]]. Un rôle particulier est dévolu aux '''équations différentielles homogènes''' (où {{math|''d'' {{=}} 0}}). Il existe une théorie générale des [[équation différentielle linéaire|équations différentielles linéaires]] (vectorielles), mais celles étudiées dans cet article comptent parmi les plus simples et les plus fréquemment rencontrées, notamment en physique.

== Équation différentielle homogène ==

L'équation est homogène lorsque {{math|''d'' {{=}} 0}}{{sfn|Deschamps|Warusfel|1999|p=459}}. Dans ce cas, une somme de deux solutions de l'équation est encore solution, ainsi que le produit d'une solution par une constante. L'ensemble des solutions est donc un [[espace vectoriel]]{{sfn|Deschamps|Warusfel|1999|p=459}} et contient notamment une solution évidente, la [[fonction nulle]].

=== À coefficients constants ===
{{Voir autre projet|wikiversity=Équation différentielle/Équation différentielle linéaire du deuxième ordre à coefficients constants#Résolution|wikiversity titre=Équation différentielle linéaire homogène d'ordre 2 à coefficients constants}}

Elles sont de la forme <math>ay'' + by' + cy = 0</math> où {{mvar|a}}, {{mvar|b}} et {{mvar|c}} sont réels ou complexes, {{mvar|a}} non nul (on peut toujours, en divisant par {{mvar|a}}, se ramener au cas {{math|''a'' {{=}} 1}}){{sfn|Deschamps|Warusfel|1999|p=459}}.

On les rencontre, entre autres, dans la modélisation de mouvement avec force de rappel (type [[ressort]]), avec ou sans amortissement<ref>{{Note autre projet|Wikiversité|Oscillateurs}}</ref> ou pendule idéal{{sfn|id=Petite Encyclopédie|texte= Petite Encyclopédie|p=555}} ou encore dans les [[Circuit électrique|circuits électriques]] comportant une [[inductance]] et un [[condensateur (électricité)|condensateur]]{{sfn|Deschamps|Warusfel|1999|p=457}}.

On cherche des solutions sous forme exponentielle, c'est-à-dire telles que {{math|''f''(''x'') {{=}} e{{exp|''λx''}}}}. Une telle fonction sera solution de l'équation différentielle si et seulement si {{mvar|λ}} est solution de
:<math>a\lambda^2 + b\lambda + c = 0.~</math>
Cette équation est appelée '''[[équation caractéristique d'une équation différentielle linéaire à coefficients constants|équation caractéristique]]''' de l'équation différentielle{{sfn|Deschamps|Warusfel|1999|p=467}}.
==== Cas réel ====
Il s'agit du cas où les coefficients sont des réels et où les solutions sont cherchées parmi les fonctions à valeurs réelles.
Comme pour toute [[équation du second degré]], trois cas se présentent selon le signe du [[discriminant]] {{math|Δ}} de l'équation caractéristique{{sfn|Deschamps|Warusfel|1999|p=469}}.

===== Si {{math|Δ > 0}} =====
L'équation caractéristique possède deux solutions réelles {{math|''λ''{{ind|1}}}} et {{math|''λ''{{ind|2}}}}.

Les solutions sont alors les fonctions définies sur ℝ par <math>f(x) =A{\rm e}^{\lambda_1x}+B{\rm e}^{\lambda_2x}</math> où {{mvar|A}} et {{mvar|B}} sont deux réels{{sfn|Deschamps|Warusfel|1999|p=469}}.

Pour déterminer ces deux constantes, il est naturel de donner deux informations sur la fonction.
*Cela se fait en général en donnant des '''conditions initiales''' en un point {{math|''x''{{ind|0}}}}, c'est-à-dire en précisant les valeurs {{math|''y''{{ind|0}}}} et {{math|''y{{'}}''{{ind|0}}}} de {{mvar|y}} et {{mvar|y'}} à cet instant. Dans ce cas l'existence et l'unicité de la solution vérifiant ces conditions initiales sont garanties{{sfn|Deschamps|Warusfel|1999|p=477}}.
*Pour de nombreux problèmes physiques, il est fréquent de donner des '''conditions aux limites''' en précisant les valeurs {{math|''y''{{ind|1}}}} et {{math|''y''{{ind|2}}}} aux instants {{math|''x''{{ind|1}}}} et {{math|''x''{{ind|2}}}}<ref>{{lien web|auteur=Gjon Koliqi|titre=Équations différentielles 3MOS|site=Gymnsase de Beaulieu|url=https://www.gymnasedebeaulieu.ch/wp-content/uploads/2019/12/equa-diff-koliqi-19.12.03.pdf |année=2019|p=16}}</ref>. Il y a alors fréquemment existence et unicité des solutions, mais ce n'est pas toujours vrai.

===== Si {{math|Δ {{=}} 0}} =====
L'équation caractéristique a une solution double {{mvar|λ}}.

Les solutions sont alors les fonctions {{mvar|f}} définies sur ℝ par <math>f(x)=(Ax + B){\rm e}^{\lambda x}</math> où {{mvar|A}} et {{mvar|B}} sont des réels quelconques{{sfn|Deschamps|Warusfel|1999|p=469}}.
 
Pour déterminer {{mvar|A}} et {{mvar|B}}, il faut, comme dans le cas précédent posséder deux informations sur {{mvar|f}}.

===== Si {{math|Δ < 0}} =====
L'équation caractéristique a deux solutions complexes conjuguées distinctes, {{math|''u'' + i''v''}} et {{math|''u'' - i''v''}}.

Les solutions sont alors les fonctions {{mvar|f}} définies sur ℝ par
: <math>f(x)={\rm e}^{ux}(A\cos(|v|x)+B\sin(|v|x))</math> où ''A'' et ''B'' sont deux réels quelconques{{sfn|Deschamps|Warusfel|1999|p=469}}.

Remarque : on peut écrire cette solution sous la forme : <math>f(x)= q{\rm e}^{ux}\cos(|v|x+r)</math>, {{mvar|q}} et {{mvar|r}} sont deux réels quelconques (cette forme est parfois plus pratique). On peut également changer la fonction cosinus en fonction sinus, la solution restant correcte.

La détermination de {{mvar|A}} et {{mvar|B}} (ou {{mvar|q}} et {{mvar|r}}) se fait, comme dans les cas précédents, par la donnée de deux informations sur {{mvar|f}}.

==== Cas complexe ====
Dans le cas où les coefficients sont complexes et où les solutions sont cherchées parmi les fonctions de la variable réelle à valeurs complexes, il ne reste alors que deux cas à étudier{{sfn|Deschamps|Warusfel|1999|p=467}} :
* Si le discriminant est non nul, l'équation caractéristique possède deux racines distinctes : {{math|''λ''{{ind|1}}}} et {{math|''λ''{{ind|2}}}} et les fonctions de ℝ dans ℂ solutions de l'équation différentielle sont les fonctions définies par <math>f(x)=C_1{\rm e}^{\lambda_1x}+C_2{\rm e}^{\lambda_2x}</math> où {{math|''C''{{ind|1}}}} et {{math|''C''{{ind|2}}}} sont deux complexes quelconques.
* Si le discriminant est nul, l'équation caractéristique ne possède qu'une seule solution {{mvar|λ}} et les solutions sont les fonctions {{mvar|f}} définies sur ℝ par <math>f(x)=(Ax + B){\rm e}^{\lambda x}</math> où {{mvar|A}} et {{mvar|B}} sont des complexes quelconques.

La détermination des constantes se fait, comme dans les cas précédents, par la donnée de deux informations sur {{mvar|f}}.

=== Équation à coefficients non constants ===

Il s'agit d'une équation de la forme <math>ay''+by'+cy=0</math>, où cette fois {{mvar|a}}, {{mvar|b}} et {{mvar|c}} sont des fonctions numériques, supposées continues sur l'intervalle d'étude ''I'', la fonction {{mvar|a}} ne s'annulant en aucun point de {{mvar|I}} (on peut toujours, en divisant par {{mvar|a}}, se ramener au cas {{math|''a'' {{=}} 1}}).

Il n'existe pas d'expression générale des solutions d'une telle équation{{sfn|id=Petite Encyclopédie|texte= Petite Encyclopédie|p=555}}. C'est pour cette raison qu'au {{s-|XIX}} furent introduites de nombreuses [[Fonction spéciale|fonctions spéciales]] ([[fonction de Bessel|fonctions de Bessel]], [[fonction d'Airy]]…) définies comme solutions d'équations qu'il est impossible de résoudre explicitement. Toutefois, dès lors qu'une solution particulière (non nulle) de l'équation est connue, il est possible de la résoudre complètement.

==== Le théorème de structure ====

Le [[théorème de Cauchy-Lipschitz]] affirme que l'ensemble ''S'' des solutions de l'équation constitue un [[Plan vectoriel|espace vectoriel de dimension deux]]{{sfn|id=Petite Encyclopédie|texte= Petite Encyclopédie|p=554}}. Dès lors, résoudre l'équation différentielle revient à exhiber deux fonctions solutions non proportionnelles : elles formeront une [[Base (algèbre linéaire)|base]] de l'espace ''S'' de solutions. Une telle base est appelée '''système fondamental de solutions'''{{sfn|id=Petite Encyclopédie|texte= Petite Encyclopédie|p=554}}.

En outre, pour tout point {{math|''x''{{ind|0}}}} de ''I'', l'application de conditions initiales en {{math|''x''{{ind|0}}}}
:<math>\begin{matrix}S & \longrightarrow & \R^2 \\ y & \longmapsto &(y(x_0),y'(x_0))\end{matrix}</math>
constitue un [[isomorphisme d'espaces vectoriels]].

==== Le wronskien ====

Étant donné deux solutions <math>y_1,y_2</math> de l'équation, leur [[wronskien]] est défini comme la fonction
:<math>W(x)=\begin{vmatrix} y_1(x) & y_2(x)\\ y'_1(x) & y'_2(x) \end{vmatrix}=y_1(x)y'_2(x)-y_2(x)y'_1(x)</math>
En utilisant les propriétés d'annulation du [[déterminant (mathématiques)|déterminant]], l'isomorphisme de conditions initiales peut être réécrit en termes de wronskien : les deux solutions forment un système fondamental si et seulement si leur wronskien est non nul en un point{{sfn|id=Petite Encyclopédie|texte= Petite Encyclopédie|p=555}}. Et dans ce cas, le wronskien ne s'annule en aucun point.

Un calcul direct montre que le wronskien de deux solutions vérifie la relation
:<math>W'(x)=-\frac{b(x)}{a(x)}W(x)</math>
qui est un cas particulier de la [[formule de Liouville]].

Il s'agit d'une [[équation différentielle linéaire d'ordre 1]] (sans second membre) qu'on sait résoudre. Donc la valeur du wronskien est connue, à une constante multiplicative près en notant {{mvar|A}} une primitive de la fonction <math>-\frac b a</math>, on obtient <math>W(x)=W_0{\rm e}^{A(x)}</math> pour une certaine constante {{math|''W''{{ind|0}}}}.

==== Application à la résolution ====

Une application fondamentale de cette propriété est la possibilité de résoudre l'équation si une solution {{math|''y''{{ind|1}}}} non nulle est connue. En effet la relation <math>W(x)=y_1(x)y'_2(x)-y_2(x)y'_1(x)</math> peut maintenant se lire comme une nouvelle équation différentielle linéaire d'ordre 1, d'inconnue {{math|''y''{{ind|2}}}}.

;Exemple
:Soit l'équation <math>x^2y''-xy'+y=0</math> à résoudre sur <math>]0,+\infty[.</math>
:Il s'agit bien d'une équation vérifiant les conditions demandées avec les fonctions {{mvar|a}}, {{mvar|b}}, {{mvar|c}} continues sur l'intervalle d'étude et {{mvar|a}} jamais nulle. La fonction définie par {{math|''y''{{ind|1}}(''x'') {{=}} ''x''}} est solution évidente. Le wronskien vérifie l'équation <math>W'(x)=\frac1xW(x)</math> ; il est de la forme{{Retrait|<math>W(x)=W_0 x = xy'_2(x)-y_2(x)</math>}}
:C'est bien une équation différentielle linéaire d'ordre 1, avec la fonction {{mvar|x}} non nulle sur l'intervalle. Elle a une solution de la forme{{Retrait|<math>y_2(x)=W_0 x\ln x +A x</math>}}
:où {{mvar|A}} est une nouvelle constante arbitraire. Nous venons d'établir que toutes les solutions de l'équation d'ordre 2 initialement considérée sont de cette forme, on y reconnaît bien la description d'un espace vectoriel de dimension 2 de solutions, un système fondamental étant donné par les solutions <math>y_1:x\mapsto x</math> et <math>x\mapsto x\ln x</math>.

== Équation avec second membre ==

=== Application du principe de superposition ===

Lorsque l'équation différentielle possède un second membre ({{mvar|d}} est une fonction non nulle), il reste possible d'exploiter ce qui précède. L'équation obtenue en remplaçant {{mvar|d}} par la fonction nulle est appelée '''équation homogène associée''' à l'équation différentielle ; on la suppose résolue.

Il suffit alors de trouver une solution {{math|''y''{{ind|0}}}} de l'équation avec second membre, pour les connaître toutes. En effet, les solutions de l'équation différentielle sont les fonctions {{math|''y''{{ind|0}} + ''g''}} où {{mvar|g}} est une solution générale de l'équation homogène associée{{sfn|id=Petite Encyclopédie|texte= Petite Encyclopédie|p=556}}.

Si le second membre {{mvar|d}} est la somme de deux fonctions {{math|''d''{{ind|1}}}} et {{math|''d''{{ind|2}}}} : <math>ay''+by'+cy=d_1+d_2</math>, on peut chercher une solution particulière {{math|''s''{{ind|1}}}} de l'équation différentielle de second membre {{math|''d''{{ind|1}}}} : <math>ay''+by'+cy=d_1</math>, puis une solution particulière {{math|''s''{{ind|2}}}} de l'équation différentielle de second membre {{math|''d''{{ind|2}}}} : <math>ay''+by'+cy=d_2</math>.
La somme {{math|''s {{=}} s''{{ind|1}} + ''s''{{ind|2}}}} de ces deux solutions particulières est solution particulière de l'équation de départ{{sfn|Deschamps|Warusfel|1999|p=475}}.
Ce principe se généralise facilement au cas où {{mvar|d}} est la somme de plus de deux fonctions.

On remarquera que les constantes d'une solution doivent être déterminées en toute fin de résolution, lorsque l'on a additionné la solution générale de l'équation homogène et la (ou les) solution(s) particulière(s) de l'équation (ou des équations) inhomogènes.

=== Recherche de solutions ayant une forme particulière ===

'''''Si {{mvar|a}}, {{mvar|b}} et {{mvar|c}} sont des constantes''''' (''a'' non nul) et si ''d'' est produit d'une [[fonction polynomiale]], d'une [[fonction trigonométrique]] et d'une [[Exponentielle de base a|fonction exponentielle]], on cherche une solution particulière de la même forme. Plus précisément :
Si ''d'' est de la forme<center><math>d:x\mapsto {\rm e}^{mx}(A(x)\cos(\omega x)+B(x)\sin(\omega x))</math></center>où ''A'' et ''B'' sont des polynômes de degré inférieur ou égal à ''n'' et où ''m'' et ω sont des réels, il existe une {{refnec|unique}} solution de la forme{{sfn|Deschamps|Warusfel|1999|p=475}}<center><math>y:x\mapsto x^p{\rm e}^{mx}(\alpha(x)\cos(\omega x)+\beta(x)\sin(\omega x))</math></center>où α et β sont des polynômes de degré inférieur ou égal à ''n'' et ''p'' peut prendre trois valeurs :
{| class="wikitable"
|-
! Par rapport à l'équation caractéristique !! Valeur de ''p''
|-
| ''m'' + [[unité imaginaire|{{math|i}}ω]] n'est pas racine de l'équation caractéristique || ''p'' = 0
|-
| ''m'' + {{math|i}}ω est racine simple de l'équation caractéristique || ''p'' = 1
|-
| ''m'' + {{math|i}}ω est racine double de l'équation caractéristique || ''p'' = 2
|}

Ceci inclut le cas où les « polynômes » ''A'' et ''B'' sont simplement des constantes (''n'' = 0), dont l'une peut être nulle, et/ou ''d'' ne contient pas d'exponentielle (''m'' = 0) et/ou pas de fonction trigonométrique (ω = 0).

{{Démonstration/début|titre=Exemples}}

* Prenons l'équation différentielle suivante :
<center><math>\ddot{y}(x)+\omega_0^2y(x)=\frac gm\qquad     (E)</math></center>
où ω{{ind|0}} est une constante non nulle. On est dans le cas ''p'' = 0. On cherche donc une solution particulière constante :

<center><math>y(x)=K.</math></center>

On calcule les dérivées par rapport à {{mvar|x}} :
<center><math>y(x)=K</math></center>
<center><math>\dot y(x)=\ddot y(x)=0</math></center>

L'équation (''E'') devient alors :

<center><math>K\omega_0^2=\frac gm</math></center>

et l'on trouve alors une solution de l'équation inhomogène :

<center><math>y(x)=K=\frac g{m\omega_0^2}.</math></center>

* Si on a une équation du type

<center><math>\ddot y(x)+b\dot y(x)=d,</math></center>
avec {{mvar|b}} non nul, on est dans le cas ''p'' = 1. On cherche donc une solution particulière de la forme ''Kt''.

* Soit une équation différentielle de la forme
<center><math>\ddot y(x)+\Gamma\dot y(x)+\omega_0^2y(x)=A_0 \cos(\omega_e x)\qquad (E)</math></center>
(c'est l'équation différentielle régissant un oscillateur harmonique périodique qui subit des frottements en présence d'une force extérieure).

On sait qu'il existe dans ce cas une solution de la forme :

<center><math>y(x)=B \cos(\omega_ex+\phi)</math></center>

On peut utiliser la notation de l'exponentielle complexe en prenant <math>Y(x)=B{\rm e}^{{\rm i}(\omega_ex+\phi)}</math>. En posant <math>B{\rm e}^{{\rm i}\phi}=\Delta</math>, on obtient <math>\Delta {\rm e}^{{\rm i}\omega_e t}</math> telle que la partie réelle de {{math|''Y''(''t'')}} soit la solution de l'équation différentielle considérée. On calcule alors les dérivées en temps :

:<math>Y(x)=\Delta {\rm e}^{{\rm i}\omega_ex}</math>
:<math>\dot Y(x)={\rm i}\omega_e\Delta {\rm e}^{{\rm i}\omega_ex}={\rm i}\omega_eY(x)</math>
:<math>\ddot Y(x)=-\omega_e^2\Delta {\rm e}^{{\rm i}\omega_ex}=-\omega_e^2Y(x).</math>

On injecte ces trois équations dans l'équation différentielle et on trouve :

<center><math>(-\omega_e^2 +\Gamma {\rm i}\omega_e +\omega_0^2)\Delta {\rm e}^{{\rm i}\omega_ex} = A_0{\rm e}^{{\rm i}\omega_ex}.</math></center>

On simplifie par les deux exponentielles (qui sont non nulles) :

<center><math>[(\omega_0^2-\omega_e^2) +{\rm i}(\omega_e\Gamma)]\Delta= A_0</math></center>

et l'on trouve alors que

<center><math>\Delta=\frac{A_0}{(\omega_0^2-\omega_e^2) +{\rm i}(\omega_e\Gamma)}.</math></center>

Ainsi, la constante {{mvar|B}} est égale au module de ce nombre complexe et la phase {{math|Φ}} à son argument :

<center><math>|\Delta|=B=\frac{A_0}{\sqrt{(\omega_0^2-\omega_e^2)^2 +(\omega_e\Gamma)^2}}</math></center>

et

<center><math>\arg(\Delta)=\phi=\arctan\left(\frac{\operatorname{Im}(\Delta)}{\operatorname{Re}(\Delta)}\right) = \arctan\left(\frac{\omega_e\Gamma}{\omega_e^2-\omega_0^2}\right).</math></center>

Ainsi, une solution particulière de l'équation inhomogène (''E'') est

<center><math>y(x)=\frac{A_0}{\sqrt{(\omega_0^2-\omega_e^2)^2 +(\omega_e\Gamma)^2}}\cos\left(\omega_ex+\arctan\left(\frac{\omega_e\Gamma}{\omega_e^2-\omega_0^2}\right)\right)</math></center>

Il est à noter que cette « méthode de l'exponentielle complexe » ne peut pas être appliquée pour une équation de la forme
<center><math>\ddot y(x)+\omega_0^2y(x)=\cos(z(x)x+\phi).</math></center>

Dans ce cas, il faut linéariser l'expression en utilisant les développements limités.
{{Démonstration/fin}}

Des considérations physiques ([[Système oscillant à un degré de liberté#Oscillations forcées|oscillations forcées]], [[résonance]]) peuvent aussi guider vers une forme particulière de solution.

=== Méthode générale ===
{{article détaillé|Méthode de variation des constantes}}

Il existe un procédé systématique de recherche des solutions, connu sous le nom de '''méthode de variation des constantes'''. Elle peut être justifiée par la théorie générale des [[équation différentielle linéaire|équations différentielles linéaires]].

Soit l'équation <math>y'' + a(x) \cdot y' + b(x) \cdot y = d(x)</math>, soient deux solutions <math>y_1, y_2</math> indépendantes de l'équation homogène (donc un système fondamental de solutions). Alors les solutions de l'équation avec second membre sont les fonctions de la forme <math>y = \lambda \cdot y_1+\mu \cdot y_2\,</math>, où les fonctions <math>\lambda, \mu</math> sont de classe <math>{\mathcal C}^1</math> et données par le système{{sfn|id=Petite Encyclopédie|texte= Petite Encyclopédie|p=556}}
:<math>\forall x \in I , \qquad \begin{cases}\lambda'(x) \cdot y_1(x)+\mu'(x) \cdot y_2(x)&=0\\
\lambda'(x) \cdot y'_1(x)+\mu'(x) \cdot y'_2(x)&=d(x)\end{cases} </math>

Dans la pratique, on écrira donc ce système, qui admet une solution pour chaque ''x''. Les fonctions solutions peuvent être [[primitive|primitivées]] et on obtient non seulement une mais toutes les solutions de l'équation avec second membre (si l'on prend en compte les constantes d'intégration dans ce système).

'''Remarque.''' Si une fonction multiplie <math>y''</math> , on pourra diviser toute l'équation par cette fonction afin de revenir au cas étudié ici.

== Les équations différentielles linéaires d'ordre deux en physique ==

Les équations d'ordre deux sont sans doute les équations différentielles les plus utilisées dans les domaines les plus variés. En particulier, les problèmes de dynamique basés sur la [[Lois du mouvement de Newton|deuxième loi de Newton]] aboutissent à une équation du deuxième ordre (ou plusieurs) dans laquelle le produit de la masse d'un corps par son accélération est égal à la somme des forces appliquées{{sfn|Deschamps|Warusfel|1999|p=457}}.

{{refsou|Sauf éventuellement des équations de la physique fondamentale}}, aucune équation représentant un phénomène réel n'est vraiment linéaire. Néanmoins l'hypothèse des petits mouvements permet souvent de négliger des petits termes non linéaires et les méthodes évoquées ci-dessus s'appliquent pour donner une solution « exacte ». C'est un problème classique dans divers domaines techniques (voir [[Systèmes oscillants à un degré de liberté]]).

Dans les cas où les mouvements ne peuvent plus être considérés comme infiniment petits, il arrive fréquemment que le système soit soumis à des forces bien individualisées : une force de rappel qui ne dépend que du mouvement, une force d'amortissement qui ne dépend que de la vitesse et une force extérieure qui ne dépend ni de l'un ni de l'autre. Dans ces conditions il existe des techniques de [[linéarisation équivalente|linéarisation]] qui permettent de calculer des coefficients adaptés à l'amplitude de l'excitation.

Enfin, on rencontre parfois des cas beaucoup plus compliqués dans lesquels on ne peut séparer les différentes forces. Par exemple, la force extérieure dépend de la position du système au même instant, éventuellement selon une loi très compliquée. Dans ces conditions les solutions analytiques doivent faire place à des solutions numériques. Parmi celles-ci, la [[méthode de Newmark]] présente une mise en œuvre très simple, une grande souplesse et une convergence sans trop de difficultés.

=== Exemple ===
En physique, on utilise souvent l'équation différentielle
:<math>\frac{d^2x}{dt^2}+2\lambda\frac{dx}{dt}+\omega_0^2 x = y</math>
L'équation différentielle homogène associée
:<math>\frac{d^2x}{dt^2}+2\lambda\frac{dx}{dt}+\omega_0^2 x = 0</math>
possède selon le signe de <math>\lambda^2-\omega_0^2</math> les solutions suivantes :

* <math>\lambda^2 > \omega_0^2 : x(t)={\rm e}^{-\lambda t}\left(A {\rm e}^{\alpha t} + B e^{-\alpha t}\right)</math>, avec <math>\alpha = \sqrt{\lambda^2-\omega_0^2}</math>, dit '''[[régime apériodique]]''',
* <math>\lambda^2 = \omega_0^2 : x(t)={\rm e}^{-\lambda t}\left(A t + B\right)</math>, dit '''[[régime critique]]''',
* <math>\lambda^2 < \omega_0^2 : x(t)={\rm e}^{-\lambda t}\left(A \cos (\omega t) + B \sin (\omega t)\right)</math>, avec <math>\omega = \sqrt{\omega_0^2-\lambda^2}</math>, dit '''[[régime pseudo-périodique]]'''.

On note aussi cette équation différentielle <math>\ddot x + 2\lambda \dot x + \omega_0^2 x = y</math> (en fonction du temps).

On peut également avoir une équation différentielle de la forme <math>a\ddot x + b\dot x + c x = d</math>. Dans ce cas, <math>2\lambda = \frac{b}{a}, \omega_0^2 = \frac{c}{a}</math> et <math>y = \frac{d}{a}</math>.

== Théorie de Sturm-Liouville ==

{{Article détaillé|Théorie de Sturm-Liouville}}

La [[théorie de Sturm-Liouville]] étudie le cas particulier des équations différentielles linéaires de la forme
:<math> -{d\over dx}\left[p(x){dy\over dx}\right]+q(x)y=\lambda w(x)y, </math>
dans laquelle le paramètre <math>\lambda</math> fait partie comme la fonction ''y'' des inconnues. Cette équation est fréquemment posée sur un segment {{math|[''a'' , ''b'']}} et accompagnée de conditions « au bord » reliant les valeurs <math>y(a),y'(a),y(b),y'(b)</math>. Les solutions {{mvar|λ}} et {{mvar|y}} du problème apparaissent alors comme [[Valeur propre (synthèse)|valeur propre et vecteur propre]] d'un certain opérateur [[Espace euclidien#Adjoint d'un endomorphisme|autoadjoint]] dans un [[espace de Hilbert]]. Le résultat principal de la théorie est l'existence d'une [[base de Hilbert|base hilbertienne]] de vecteurs propres associés à des valeurs propres formant une suite strictement croissante ([[Théorème de Riesz-Fredholm]]).

Par multiplication par un [[facteur intégrant]] convenable, toute équation différentielle linéaire d'ordre deux peut être mise sous la forme d'une équation de Sturm-Liouville.

== Références ==
{{Références}}
== Bibliographie==
* {{Demailly1}}
* {{ouvrage|auteur1=[[Claude Deschamps]]|auteur2=[[André Warusfel]]|titre=Mathématiques 1ère année MPSI, PCSI, PTSI|collection =J'intègre|éditeur=Dunod|année=1999}}
* {{ouvrage
 | prénom1 = W |nom1 = Gellert
 | prénom2 = H. | nom2 = Küstner
 | prénom3 = M. | nom3 = Hellwich
 | prénom4 = H. | nom4 = Kästner
 | titre = Petite encyclopédie des mathématiques
 | éditeur = Didier
 | année = 1980
 | id = Petite Encyclopédie 
}}

== Voir aussi ==
{{Autres projets
  |v=Équation différentielle/Équation différentielle linéaire du deuxième ordre à coefficients variables
}}

* [[Équation différentielle]]
* [[Équation différentielle linéaire d'ordre un]]
* [[Histoire des équations]]

{{Portail|analyse}}

[[Catégorie:Équation différentielle|Lineaire d'ordre deux]]