Voir la source de Théorie des perturbations

La '''théorie des perturbations''' est un domaine des [[mathématiques]], qui consiste à étudier les contextes où il est possible de trouver une ''solution approchée'' à une [[équation]] en partant de la solution d'un problème plus simple. Plus précisément, on cherche une solution approchée à une équation {{math|(''E''{{ind|λ}})}} (dépendante d'un paramètre {{math|λ}}), sachant que la solution de l'équation {{math|(''E''{{ind|0}})}} (correspondant à la valeur {{math|λ{{=}}0}}) est connue exactement. L'équation mathématique {{math|(''E''{{ind|λ}})}} peut être par exemple une [[équation algébrique]] ou une [[équation différentielle]]. La méthode consiste à chercher la solution approchée de l'équation {{math|(''E''{{ind|λ}})}} sous la forme d'un développement en [[série (mathématiques)|série]] des puissances du paramètre {{math|λ}}, cette solution approchée étant supposée être une approximation d'autant meilleure de la solution exacte, mais inconnue, que la [[valeur absolue]] du paramètre {{math|λ}} est plus « petite ».

== Généralités ==
=== Histoire ===

Dès le début du {{XVIIIe siècle}}, la théorie des perturbations a été utilisée par les astronomes pour les besoins de la [[mécanique céleste]] : en effet, les [[équation différentielle|équations différentielles]] décrivant un système de {{mvar|N}} [[problème à N corps|corps]] en interaction [[gravitation]]nelle n'a pas de solution exacte générale<ref>Il existe en fait une solution exacte du problème des trois corps découverte par Sundman (1909). Cette solution exacte, sous la forme d'un développement en série formelle, n'est en pratique pas exploitable, car si la série converge bien « au sens des géomètres », elle le fait si lentement que cela rend son pouvoir prédictif quasi nul. Lire : Malte Henkel ; ''Sur la solution de Sundman du problème des trois corps'', Philosophia Scientiae 5(2) (2001) pp. 161-184. Texte complet disponible sur l'ArXiv : [http://fr.arxiv.org/abs/physics/0203001 ''physics/0203001''].</ref> pour {{math|''N'' ≥ 3}}. Cet aspect de la théorie des perturbations a été synthétisé à la fin du {{XIXe siècle}} dans les ouvrages classiques de [[Pierre-Simon Laplace|Laplace]]<ref>Pierre-Simon Laplace ; ''Traité de mécanique céleste'', Editions Jacques Gabay (1990). Réédition d'un ouvrage classique de la fin du {{s-|XIX|e}}, en 4 volumes. Niveau second cycle universitaire. Cet ouvrage est disponible en fac-similé sur [https://gallica.bnf.fr/document?O=N077586 ''Gallica''].</ref>, [[Félix Tisserand|Tisserand]]<ref>François-Félix Tisserand ; ''Traité de mécanique céleste'', Editions Jacques Gabay (1990). Réédition d'un ouvrage classique de la fin du {{s-|XIX|e}}, en 4 volumes. Niveau second cycle universitaire. Cet ouvrage est disponible en fac-similé sur [https://gallica.bnf.fr/document?O=N029031 ''Gallica''].</ref> et [[Henri Poincaré|Poincaré]]<ref name="Poincaré" />, avant de connaître de nouveaux développements dans la seconde moitié du {{XXe siècle}} avec l'avènement en 1954 de la « [[théorème KAM|théorie KAM]] », du nom de ses trois concepteurs : [[Andreï Kolmogorov|Kolmogorov]], [[Vladimir Arnold|Arnold]] et [[Jürgen Moser|Moser]].

La méthode a par ailleurs été abondamment utilisée au {{XXe siècle}} pour les besoins de la [[physique quantique]], d'abord en [[mécanique quantique]] non relativiste, puis en [[théorie quantique des champs]] perturbative.

=== Convergence de la série perturbative ? ===

On a vu qu'on cherchait ici la solution approchée de l'équation {{math|(''E''{{ind|λ}})}} sous la forme d'un développement en série des puissances du paramètre {{math|λ}} ; la question de la ''convergence'' de cette série se pose alors. Ce problème a été réglé pour l'astronomie par Poincaré en 1892 : la « série » de perturbation doit être comprise mathématiquement comme un [[développement asymptotique]] au voisinage de zéro, et non comme une [[série (mathématiques)|série]] ordinaire convergente uniformément. Le chapitre VIII de la ''Mécanique céleste'' de Poincaré<ref>Henri Poincaré ; ''Les méthodes nouvelles de la mécanique céleste'', Gauthier-Villars (1892).</ref> commence par le commentaire suivant :

{{Citation bloc|II y a entre les géomètres et les astronomes une sorte de malentendu au sujet de la signification du mot convergence. Les géomètres préoccupés de la parfaite rigueur et souvent trop indifférents à la longueur de calculs inextricables dont ils conçoivent la possibilité, sans songer à les entreprendre effectivement, disent qu'une série est convergente quand la somme des termes tend vers une limite déterminée, quand même les premiers termes diminueraient très lentement. Les astronomes, au contraire, ont coutume de dire qu'une série converge quand les 20 premiers termes, par exemple, diminuent très rapidement, quand même les termes suivants devraient croître indéfiniment.

Ainsi, pour prendre un exemple simple, considérons les deux séries qui ont pour terme général :

<center>
<math>
\frac{1000^n}{n !} \quad \mathrm{et : } \quad \frac{n !}{1000^n}
</math>
</center>

Les géomètres diront que la première converge, et même qu'elle converge rapidement, parce que le millionième terme est beaucoup plus petit que le 999999{{e}} ; mais ils regarderont la seconde comme divergente, parce que le terme général peut croître au-delà de toute limite.

Les astronomes, au contraire, regarderont la première série comme divergente, parce que les 1000 premiers termes vont en croissant ; et la seconde comme convergente, parce que les 1000 premiers termes vont en décroissant et que cette décroissance est d'abord très rapide.

Les deux règles sont légitimes : la première, dans les recherches théoriques ;  la seconde, dans les applications numériques. Toutes deux doivent régner, mais dans deux domaines séparés et dont il importe de bien connaître les frontières.}}

Pour conclure cette discussion qualitative sur la convergence, le mathématicien [[Jean-Pierre Ramis]] précise<ref>Jean-Pierre Ramis, [http://www.math.polytechnique.fr/xups/xups91.pdf ''Séries divergentes et théories asymptotiques''], Journées X-UPS (1991).</ref> :

{{Citation bloc|On peut ainsi parler de séries convergentes « au sens des géomètres » ou « au sens des astronomes ». Notons que pratiquement, dans les applications, on constate que, presque toujours, les séries convergentes au sens des astronomes ont un terme général qui croît très vite après avoir d'abord diminué. Ainsi, ce que Poincaré envisageait comme possibilité est en fait la règle.}}

C'est d'ailleurs ce qui fait l'efficacité pratique de la théorie des perturbations en [[physique théorique]] : il suffit le plus souvent de calculer les quelques premiers termes du développement asymptotique - ceux qui semblent commencer par converger - pour obtenir une très bonne approximation du résultat exact inconnu. Ainsi, dans le cadre de l'[[électrodynamique quantique]], [[Freeman Dyson|Dyson]] a montré en 1948 que la [[série perturbative]] était ''divergente'', alors que la prise en compte des trois ou quatre premiers termes seulement donne des prédictions théoriques en accord remarquable avec les résultats expérimentaux.

Notons qu'il existe certaines procédures de « sommation » qui permettent de donner un sens à certaines [[série divergente|séries divergentes]], comme par exemple la [[sommation de Borel]] ou l'[[approximant de Padé]].

== Un premier exemple élémentaire ==

=== Position du problème ===

Considérons à titre d'exemple l'[[Équation différentielle d'ordre un à variables séparées|équation différentielle d'ordre un]] suivante :

<center>
<math> (E_{\lambda}) \ : \qquad \frac{\mathrm{d}x(t)}{\mathrm{d}t} + \frac{1}{\tau} x(t) + \frac{\lambda}{\tau L_0} x^2(t) = 0 </math> .
</center>

Dans cette équation, {{mvar|t}} représente le temps, {{math|τ}} un paramètre fixé [[analyse dimensionnelle|homogène]] à un temps, {{math|''L''{{ind|0}}}} un paramètre fixé homogène à une longueur, et {{math|λ}} le paramètre de perturbation, [[analyse dimensionnelle|sans dimensions]]. On cherche à déterminer la fonction {{math|''x''(''t'')}} inconnue, homogène à une longueur, et vérifiant la ''condition initiale'' : à l'instant ''t'' = 0, on a {{math|''x''(0){{=}}''X''{{ind|0}}}}.

=== Théorie des perturbations au premier ordre ===
{{Section à vérifier|date=mars 2019}}
Le problème de départ de la théorie des perturbations est l'équation différentielle {{math|(''E''{{ind|0}})}} correspondant à la valeur {{math|λ{{=}}0}} :

{| align="center" border="0"
|<math> (E_{0}) \ : \qquad \frac{\mathrm{d}x(t)}{\mathrm{d}t} \ + \ \frac{1}{\tau} \ x(t) \ = \ 0 </math>
|}

dont la solution analytique exacte est bien connue :

{| align="center" border="0"
|<math> x_{0}(t) = A \mathrm{e}^{- \, t/\tau}</math>
|}

où {{mvar|A}} est une constante, pour l'instant inconnue. On illustre la méthode des perturbations en nous limitant pour simplifier au ''premier ordre'' dans le développement en série des puissances du paramètre {{math|λ}} ; on cherche donc la solution approchée sous la forme :

{| align="center" border="0"
|<math> x_{\lambda}(t) = x_{0}(t) + \lambda x_{1}(t) + O(\lambda^2)</math>
|}

où {{math|''x''{{ind|1}}(''t'')}} est une fonction inconnue, à déterminer. On injecte cette expression dans l'équation différentielle exacte {{math|(''E''{{ind|λ}})}}. En se limitant aux termes du premier ordre inclus et en utilisant le fait que {{math|''x''{{ind|0}}(''t'')}} est la solution exacte de {{math|(''E''{{ind|0}})}}, on obtient la solution physique approchée au premier ordre :

<center>
<math> x_{\lambda}(t) = X_0 \left( 1 - \lambda \frac{X_0}{L_0} \right) \mathrm{e}^{- t/\tau} + \lambda \frac{X_0^2}{L_0} \mathrm{e}^{- 2t/\tau} + O(\lambda^2)</math>
</center>

{{Démonstration|contenu=
On injecte le développement dans l'équation différentielle exacte {{math|(''E''{{ind|λ}})}}. En se limitant aux termes du premier ordre inclus et en utilisant le fait que {{math|''x''{{ind|0}}}} est la solution exacte de {{math|(''E''{{ind|0}})}}, on obtient l'équation suivante pour la fonction {{math|''x''{{ind|1}}}} :
<center>
<math>\lambda \ \left[ \ \frac{\mathrm{d}x_1(t)}{\mathrm{d}t} +\frac{1}{\tau} x_1(t) \right] +\lambda \frac{x^2_0(t)}{\tau L_0}  = 0,  </math>
</center>
On omettra dans cette démonstration les restes <math>O(\lambda^2)</math> pour simplifier les écritures. Cette équation se réécrit explicitement :
<center>
<math> \frac{\mathrm{d}x_1(t)}{\mathrm{d}t} + \frac{1}{\tau} x_1(t) = - \frac{A^2}{\tau L_0} \mathrm{e}^{- 2t/\tau}. </math>
</center>
Les solutions de l'équation homogène sont de la forme :
<center>
<math> t \mapsto B\mathrm{e}^{-t/\tau}, </math>
</center>
où {{mvar|B}} parcourt l'ensemble des réels. Pour déterminer une solution particulière de l'équation complète, on peut utiliser une méthode de variation de la constante.
On cherche donc une solution particulière {{math|''x''{{ind|part}}}} sous la forme :
<center>
<math> x_{\textrm{part}}(t) = C(t) \mathrm{e}^{- t/\tau}, </math>
</center>
où {{mvar|C}} est une fonction qui dépend du temps. En injectant cette solution particulière dans l'équation vérifiée par {{math|x{{ind|1}}}}, on obtient que {{mvar|C}} doit vérifier l'équation :
<center>
<math> C'(t) = -\frac{A^2}{\tau L_0} \mathrm{e}^{-t/\tau}, </math>
</center>
d'où :
<center>
<math> C(t) = \frac{A^2}{L_0} \mathrm{e}^{-t/\tau}, </math>
</center>
convient, et alors :
<center>
<math> x_{\textrm{part}}(t) = C(t) \mathrm{e}^{-t/\tau} = \frac{A^2}{L_0} \mathrm{e}^{-2t/\tau} </math>
</center>
est une solution particulière de l'équation vérifiée par {{math|''x''{{ind|1}}}}.
Ainsi, on a :
<center>
<math> x_{1}(t) = B_0 \mathrm{e}^{- t/\tau} + \frac{A^2}{L_0} \mathrm{e}^{-2t/\tau}, </math>
</center>
où {{math|''B''{{ind|0}}}} est une constante qui dépend de la condition initiale implicite et que l'on va déterminer plus loin. 
Il est maintenant possible d'exprimer <math> x_\lambda = x_0 + \lambda x_1 </math> (on omet encore les restes dans l'écriture) :
<center>
<math> x_\lambda (t) = ( A+\lambda B_0) \mathrm{e}^{-t/\tau} + \lambda \frac{A^2}{L_0} \mathrm{e}^{-2t/\tau}. </math>
</center>
En se rappelant que {{math|''x''{{ind|λ}}}} vérifie la condition initiale {{math|''x''{{ind|λ}}(0) {{=}} ''X''{{ind|0}}}}, on obtient, en faisant {{math|''t'' {{=}} 0}} dans l'expression ci-dessus :
<center>
<math> A + \lambda B_0 + \lambda \frac{A^2}{L_0} = X_0. </math>
</center>
Pour {{math|λ}} proche de 0, on fait l'approximation {{math|''A'' {{=}} ''x''{{ind|0}}}} et donc {{math|''B''{{ind|0}} {{=}} -''X''{{indexp|0|2}}/''L''{{ind|0}}}}.
Ainsi, on trouve que :
<center>
<math> x_{\lambda}(t) = \left(X_0 - \lambda \frac{X_0^2}{L_0}\right) \mathrm{e}^{-t/\tau} + \lambda \frac{X_0^2}{L_0} \mathrm{e}^{-2t/\tau}, </math>
</center>
ce qui peut se réécrire, en mettant {{math|''X''{{ind|0}}}} en facteur dans le premier terme :
<center>
<math> x_{\lambda}(t) = X_0\left(1 - \lambda \frac{X_0}{L_0}\right) \mathrm{e}^{-t/\tau} + \lambda \frac{X_0^2}{L_0} \mathrm{e}^{-2t/\tau}, </math>
</center>
C.Q.F.D.
}}

=== Comparaison avec la solution exacte ===

On peut démontrer ici que l'équation différentielle {{math|(''E''{{ind|λ}})}} vérifiant la [[condition initiale]] : {{math|''x''(0){{=}}''X''{{ind|0}}}} admet pour toutes les valeurs du paramètre {{math|λ}} la solution exacte suivante :

<center>
<math> x_{\lambda}(t) = \frac{X_0 \mathrm{e}^{- t/\tau}}{1 + \lambda \frac{X_0}{L_0}\left( 1 - \mathrm{e}^{- t/\tau} \right)}</math>
</center>

Un [[développement limité]] de cette expression au premier ordre en {{math|λ}} donne explicitement la solution approchée déterminée au paragraphe précédent par la théorie des perturbations au premier ordre :

<center>
<math> x_{\lambda}(t) = X_0 \left(  1 - \lambda \frac{X_0}{L_0} \right) \mathrm{e}^{- t/\tau} + \lambda \frac{X_0^2}{L_0} \mathrm{e}^{- 2t/\tau} + O(\lambda^2)</math>
</center>

Pour visualiser l'écart entre la solution approchée et la solution exacte, on trace ci-dessous les graphes des deux fonctions pour une série de valeurs de {{math|λ}} allant de 0,1 à 0,5, en prenant : {{math|''X''{{ind|0}} {{=}} ''L''{{ind|0}} {{=}} 1 m, τ {{=}} 1 s|§=''X''{{ind|0}} {{=}} ''L''{{ind|0}} {{=}} 1 m, τ {{=}} 1 s}}.
* en bleu, la solution exacte.
* en rouge, la solution approchée au premier ordre.
<gallery>
Fichier:Perturb1.jpg|{{math|λ{{=}}0,1}}
Fichier:Perturb2.jpg|{{math|λ{{=}}0,2}}
Fichier:Perturb3.jpg|{{math|λ{{=}}0,3}}
Fichier:Perturb4.jpg|{{math|λ{{=}}0,4}}
Fichier:Perturb5.jpg|{{math|λ{{=}}0,5}}
</gallery>

== Un deuxième exemple : l'oscillateur de Duffing ==

=== Définition et propriétés ===

==== Définition ====

L'[[oscillateur de Duffing]] satisfait à l'équation différentielle du second ordre suivante :

<center>
<math> (E_{\lambda}) \ : \qquad \frac{\mathrm{d}^2x(t)}{\mathrm{d}t^2} + \omega_0^2 \ x(t) +  \frac{\lambda \omega_0^2}{L_0^2} x^3(t) = 0 </math>
</center>

Dans cette équation, {{mvar|t}} représente le temps, {{math|ω{{ind|0}}}} un paramètre fixé [[analyse dimensionnelle|homogène]] à une pulsation, c’est-à-dire l'inverse d'un temps. {{math|''L''{{ind|0}}}} est un paramètre fixé homogène à une longueur, et {{math|λ}} le paramètre de perturbation, [[analyse dimensionnelle|sans dimensions]]. On cherche à déterminer la fonction {{math|''x''(''t'')}} inconnue, homogène à une longueur, et vérifiant les ''conditions initiales'' : à l'instant {{math|''t''{{=}}0}}, on a {{math|''x''(0){{=}}''X''{{ind|0}}}} et <math>\dot{x}(0) = 0</math>.

==== Interprétation physique ====

On peut interpréter cette équation différentielle comme la loi de la dynamique de Newton d'une particule de masse {{mvar|m}} soumise à une force dérivant d'une [[énergie potentielle]] {{math|''V''(''x'')}} :

<center>
<math> m \ \frac{\mathrm{d}^2x(t)}{\mathrm{d}t^2} = - \frac{\mathrm{d}V(x(t))}{\mathrm{d}x}</math>
</center>

où le potentiel {{math|''V''(''x'')}} quartique s'écrit :

<center>
<math> V(x) = m \, \omega_0^2 \ \left( \ \frac{x^2}{2} + \frac{\lambda \ x^4}{4 L_0^2} \ \right)</math>
</center>

==== Caractère borné du mouvement ====

Pour toutes les valeurs de {{math|λ}} ''positives'' ou nulles, {{math|''V''(''x'')}} représente un [[puits de potentiel]]. La [[conservation de l'énergie]] mécanique totale {{mvar|E}} de la particule :

<center>
<math>E = \frac{1}{2}m v^2 + V(x) \quad \Longrightarrow \quad \frac{1}{2}m v^2 = E - V(x) \ \ge \ 0</math>
</center>

entraîne alors que le mouvement est ''borné'' dans un intervalle {{math|[''x''{{ind|1}}, ''x''{{ind|2}}]}}, où les points tournants {{math|''x''{{ind|1}}}} et {{math|''x''{{ind|2}}}} sont les deux solutions réelles de l'équation {{math|''E''{{=}}''V''(''x'')}}

==== Ordre zéro : l'oscillateur harmonique ====

Le problème de départ de la théorie des perturbations est l'équation différentielle {{math|(''E''{{ind|0}})}} correspondant à la valeur {{math|λ{{=}}0}} :

<center>
<math> (E_{0}) \ : \qquad \frac{\mathrm{d}^2x(t)}{\mathrm{d}t^2} + \omega_0^2 x(t) = 0 </math>
</center>

Cette équation est par définition un [[oscillateur harmonique]] de pulsation {{math|ω{{ind|0}}}}, dont la solution analytique exacte est bien connue :

<center>
<math> x_{0}(t) = A \ \cos \left( \omega_0 t + \varphi \right) </math>
</center>

où {{mvar|A}} et {{math|φ}} sont deux constantes, pour l'instant inconnues.

=== Théorie de perturbation naïve au premier ordre ===

On cherche la solution approchée sous la forme :

<center>
<math> x_{\lambda}(t) = x_{0}(t) + \lambda x_{1}(t) + O(\lambda^2)</math>
</center>

où {{math|''x''{{ind|1}}(''t'')}} est une fonction inconnue, à déterminer. On injecte cette expression dans l'équation différentielle exacte {{math|(''E''{{ind|λ}})}}. En se limitant aux termes du premier ordre inclus et en utilisant le fait que {{math|''x''{{ind|0}}(''t'')}} est la solution exacte de {{math|(''E''{{ind|0}})}}, on obtient l'expression au premier ordre de la théorie de perturbation :

<center>
<math> x_{\lambda}(t) = X_0 \ \cos \left( \omega_0  t \right) + \lambda \ \left[ \frac{X_0^3}{32 \, L_0^2} \left[ \cos \left( 3 \omega_0 t \right) - \cos \left( \omega_0 t \right) \right] + \frac{3 X_0^3}{8 L_0^2} \omega_0 t \sin \left( \omega_0 t \right) \right] + O(\lambda^2) </math>
</center>

{{Démonstration|contenu=

On injecte le développement dans l'équation différentielle exacte {{math|(''E''{{ind|λ}})}}. En se limitant aux termes du premier ordre inclus et en utilisant le fait que {{math|''x''{{ind|0}}(''t'')}} est la solution exacte de {{math|(''E''{{ind|0}})}}, on obtient l'équation suivante pour la fonction {{math|''x''{{ind|1}}(''t'')}} :

<center>
<math> \frac{\mathrm{d}^2x_1(t)}{\mathrm{d}t^2} + \omega_0^2 x_1(t) = - \frac{\omega_0^2 A^3}{L_0^2} \cos^3 \left( \omega_0 t + \varphi \right)</math>
</center>

On utilise l'identité trigonométrique suivante :

<center>
<math> \cos^3 \alpha = \frac{1}{4} \ \left[ 3 \cos \alpha + \cos (3 \alpha ) \right]</math>
</center>

d'où l'équation différentielle pour la fonction {{math|''x''{{ind|1}}(''t'')}} :

<center>
<math> \frac{d^2x_1(t)}{dt^2} + \omega_0^2 x_1(t) = - \frac{\omega_0^2 A^3}{4 L_0^2} \left[ 3 \cos \left( \omega_0 t + \varphi \right) + \cos \left( 3 \omega_0 t + 3 \varphi \right) \right] </math>
</center>

On peut démontrer que cette équation différentielle possède la solution générale exacte suivante :

<center>
<math> x_1(t) = + \frac{3 A^3}{8 L_0^2} \omega_0 t \sin \left( \omega_0 t + \varphi \right) + \frac{A^3}{32 L_0^2} \cos \left( 3 \omega_0 t + 3 \varphi \right) </math>
</center>

On a donc pour la fonction inconnue :

<center>
<math> x_{\lambda}(t) = A \cos \left( \omega_0 t + \varphi  \right) + \lambda \left[ \frac{3 A^3}{8 L_0^2} \omega_0 t \sin \left( \omega_0 t + \varphi \right) + \frac{A^3}{32 L_0^2} \cos \left( 3 \omega_0 t + 3 \varphi \right) \right] + O(\lambda^2)</math>
</center>

L'application des conditions initiales : à l'instant {{math|''t''{{=}}0}}, on a : {{math|''x''(0){{=}}''X''{{ind|0}}}} et <math>\dot{x}(0) = 0</math> conduisent au système de deux équations à deux inconnues :

<center>
<math>A \cos \varphi + \lambda \frac{A^3}{32 L_0^2} \cos (3 \varphi) = X_0 \quad (1)</math>
</center>

<center>
<math>- \ A \ \omega_0 \ \sin (\varphi) \ + \ \lambda \ \left[ \ \frac{3 \, \omega_0 \, A^3}{8 \, L_0^2} \ \sin \varphi  \ - \ \frac{3 \omega_0 \, A^3}{32 \, L_0^2} \ \  \sin(3 \varphi) \ \right] \ = \ 0 \quad (2)</math>
</center>

En utilisant la formule trigonométrique :

<center>
<math>\sin(3 \varphi) \ = \ 3 \, \cos^2 \varphi \, \sin \varphi \ - \ \sin^3 \varphi</math>
</center>

dans l'équation (2), on montre que {{math|sin(φ){{=}}0}}, donc que {{math|φ{{=}}0}}. On reporte ensuite ce résultat dans l'équation (1), ce qui donne :

<center>
<math>A \  + \ \lambda \ \frac{A^3}{32 \, L_0^2} \ = \ X_0 \quad \Longrightarrow \quad A \ = \ X_0 \ \left[ \, 1 \, - \, \frac{\lambda \, X_0^2}{32 \, L_0^2} \, \right] \ + \ O(\lambda^2)</math>
</center>

On en déduit l'expression au premier ordre de la théorie de perturbation.
}}

=== Apparition d'un terme séculaire ===

On constate que la perturbation contient un terme ''proportionnel au temps'' :

<center>
<math>\lambda \ \omega_0 t \ \sin \left(  \omega_0  t \right)</math>
</center>
Ce terme ''non borné'' est appelé ''terme séculaire'', du mot latin ''saeculum'' qui signifie siècle. En effet, pour les temps <math> t \sim 1/\omega_0</math>, la perturbation est bien d'ordre {{math|λ}}, c’est-à-dire petite. En revanche, pour des temps plus longs de l'ordre de <math> t \sim 1/(\lambda \omega_0)</math>, la perturbation devient d'ordre 1 et n'est plus petite ;  le problème devient encore pire pour des temps encore plus longs : <math> t \gg 1/(\lambda \omega_0)</math>. Or nous savons que le mouvement réel est borné, donc que {{math|''x''{{ind|λ}}(''t'')}} ne peut pas croître indéfiniment : notre théorie des perturbations « naïve » n'est donc plus valide.

Dans le cadre de l'astronomie, la présence de ces termes séculaires empêche d'étudier le futur à long terme des trajectoires planétaires, l'unité de temps caractéristique du problème étant le siècle.

== Méthode de Lindstedt ==

[[Anders Lindstedt|Lindstedt]] a proposé en 1882 une méthode qui, pour certaines équations différentielles, permet d'éliminer ces termes séculaires<ref>Anders Lindstedt ; Abh. K. Akad. Wiss. St. Petersburg 31 (1882), 4.</ref>{{,}}<ref>Anders Lindstedt ; ''Sur la forme des expressions des distances mutuelles dans le problème des trois corps'', Comptes-Rendus de l'Académie des Sciences 97 (1883) 1276 et 1353. Ces deux notes sont disponibles en fac-similé sur [https://gallica.bnf.fr/document?O=N003053 ''Gallica''].</ref>. On l'appelle aussi ''méthode de Lindstedt-Poincaré'', Poincaré ayant démontré que les séries introduites par Lindstedt devaient être interprétées comme des expressions ''asymptotiques''<ref>Henri Poincaré ; ''Sur les séries de M. Lindstedt'', Comptes-Rendus de l'Académie des Sciences 108 (1889) 21-24. Cette note est disponible en fac-similé sur [https://gallica.bnf.fr/document?O=N003064 ''Gallica''].</ref>{{,}}<ref>Henri Poincaré ; ''Sur l'application de la méthode de M. Lindstedt au problème des trois corps'', Comptes-Rendus de l'Académie des Sciences 114 (1892) 1305-1309. Cette note est disponible en fac-similé sur [https://gallica.bnf.fr/document?O=N003070 ''Gallica''].</ref>{{,}}<ref name="Poincaré" />.

L'idée de Lindstedt est la suivante : dans certains cas, les termes séculaires peuvent être dus au fait que l'on développe incorrectement les expressions. Par exemple, supposons que le résultat exact soit :

{| align="center" border="0"
|<math>x_{\lambda}(t) = X_0 \cos \left[ ( 1 + \lambda ) \omega_0  t \right]</math>
|}

Cette expression clairement ''bornée'' développée au premier ordre en {{math|λ}} donne :

{| align="center" border="0"
|<math>x_{\lambda}(t) = X_0 \ \cos \left( \omega_0  t \right) - \lambda X_0 \omega_0 t \sin \left( \omega_0  t \right) + O(\lambda^2)</math>
|}

et il apparaît un terme séculaire non borné. On voit que la solution exacte est en fait une fonction de la pulsation :

{| align="center" border="0"
|<math>\omega_{\lambda} = ( 1 +  \lambda ) \omega_0  </math>
|}

qui est légèrement différente de la pulsation initiale {{math|ω{{ind|0}}}} du problème. Lindstedt va utiliser cette remarque de façon systématique.

=== Principe de la méthode de Lindstedt ===

La méthode de Lindstedt ne s'applique que pour les équations différentielles du type suivant :

<center>
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-
| bgcolor="#fff8ff" |<math> \frac{\mathrm{d}^2x(t)}{\mathrm{d}t^2} + \omega_0^2 \ x(t) \ = \ \lambda \ f \left(x(t), \, \frac{\mathrm{d}x(t)}{\mathrm{d}t}, \, t, \, \lambda \right) </math>
|}</center>

où {{mvar|f}} est une fonction ''paire'' de {{mvar|x}} et ''impaire'' de <math>\dot{x}</math> qui est de plus soit ''périodique'' en {{mvar|t}}, soit indépendante de {{mvar|t}}. La méthode consiste à faire un ''changement d'échelle de temps'' en introduisant une nouvelle variable {{mvar|s}} [[analyse dimensionnelle|sans dimensions]] définie par le développement en série :

<center>
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-
| bgcolor="#fff8ff" |<math> s = \left( \omega_0  + \lambda \omega_1 + \lambda^2 \omega_2 + \dots \ \right) \ t</math>
|}</center>

Dans cette expression, les valeurs numériques des constantes inconnues {{math|ω{{ind|''k''}}, ''k'' ≥ 0}} devront être choisies afin de faire disparaître les termes séculaires de la série perturbatrice de la solution approchée à l'ordre désiré.

Illustrons la méthode dans le paragraphe suivant avec l'oscillateur de Duffing.

=== Exemple : l'oscillateur de Duffing au premier ordre ===

On a vu plus haut que l'équation différentielle de l'oscillateur de Duffing s'écrivait :

{| align="center" border="0"
|<math> (E_{\lambda}) \ : \qquad \frac{\mathrm{d}^2x(t)}{\mathrm{d}t^2} + \omega_0^2 x(t) = - \lambda \frac{\omega_0^2}{L_0^2} x^3(t) </math>
|}

Faisons le changement d'échelle de temps {{math|''t'' → ''s''}} et définissons la nouvelle fonction inconnue {{mvar|y}} de la variable {{mvar|s}} par :

{| align="center" border="0"
|<math>y(s) = x(t(s)) = (x \circ t)(s)</math>
|}

La règle de Leibniz de dérivation en chaîne donne pour la dérivée première :

{| align="center" border="0"
|<math>\frac{\mathrm{d}x}{\mathrm{d}t} = \frac{\mathrm{d}y}{\mathrm{d}s} \ \frac{\mathrm{d}s}{\mathrm{d}t} </math>
|}

et pour la [[dérivée seconde]] :

{| align="center" border="0"
|<math>\frac{\mathrm{d}^2x}{\mathrm{d}t^2} = \frac{\mathrm{d}}{\mathrm{d}s} \left( \frac{\mathrm{d}y}{\mathrm{d}s} \ \frac{\mathrm{d}s}{\mathrm{d}t} \right) \ \frac{\mathrm{d}s}{\mathrm{d}t} = \frac{\mathrm{d}^2y}{\mathrm{d}s^2} \ \left( \frac{\mathrm{d}s}{\mathrm{d}t} \right)^2 </math>
|}

Comme on a au premier ordre :

{| align="center" border="0"
|<math>\frac{\mathrm{d}s}{\mathrm{d}t} = \omega_0 + \lambda \omega_1 + O(\lambda^2)</math>
|}

on obtient pour les dérivées :

{| align="center" border="0"
|<math>\frac{\mathrm{d}x}{\mathrm{d}t} = ( \omega_0 + \lambda \omega_1 ) \frac{\mathrm{d}y}{\mathrm{d}s} + O(\lambda^2) </math>
|}

{| align="center" border="0"
|<math>\frac{\mathrm{d}^2x}{\mathrm{d}t^2} = \left( \omega_0^2 + 2 \lambda \omega_0 \omega_1 \right) \frac{\mathrm{d}^2y}{\mathrm{d}s^2} + O(\lambda^2)</math>
|}

L'équation différentielle de Duffing devient au premier ordre :

{| align="center" border="0"
|<math> \left( \omega_0^2 + 2 \lambda \omega_0 \omega_1 \right) \frac{\mathrm{d}^2y(s)}{\mathrm{d}s^2} + \omega_0^2 y(s) = - \lambda \frac{\omega_0^2}{L_0^2} y^3(s) + O(\lambda^2)  </math>
|}

Introduisons maintenant dans cette équation différentielle le développement au premier ordre de la solution :

{| align="center" border="0"
|<math> y_{\lambda}(s) = y_{0}(s) + \lambda y_{1}(s) + O(\lambda^2)</math>
|}

Il vient en développant selon les puissances de {{math|λ}} :

{| align="center" border="0"
|<math> \omega_0^2 \left(\frac{\mathrm{d}^2y_0}{\mathrm{d}s^2} + y_{0}  \right) + \lambda \left[ \omega_0^2 \frac{\mathrm{d}^2y_1}{\mathrm{d}s^2} + 2 \omega_0 \omega_1  \frac{\mathrm{d}^2y_0}{\mathrm{d}s^2} + \omega_0^2  y_1  +  \frac{\omega_0^2}{L_0^2} \ y_0^3 \right] = 0  +  O(\lambda^2)  </math>
|}

On a donc le système de deux équations différentielles :

{| align="center" border="0"
|<math> \frac{\mathrm{d}^2y_0}{\mathrm{d}s^2} + y_{0} = 0</math>
|}

{| align="center" border="0"
|<math>\omega_0^2 \frac{\mathrm{d}^2y_1}{\mathrm{d}s^2} + 2 \omega_0 \omega_1  \frac{\mathrm{d}^2y_0}{\mathrm{d}s^2}  + \omega_0^2  y_1  + \frac{\omega_0^2}{L_0^2} y_0^3 = 0 </math>
|}

La première a pour solution générale : {{math|''y''{{ind|0}}(''s''){{=}}''A'' cos (''s'' + φ)}} où {{mvar|A}} et {{math|φ}} sont deux constantes. On reporte alors cette expression dans la seconde équation différentielle, et on obtient pour la première correction {{math|''y''{{ind|1}}(''s'')}} :

{| align="center" border="0"
|<math>\frac{\mathrm{d}^2y_1(s)}{\mathrm{d}s^2} + y_1(s) = \frac{2 \omega_1 A}{\omega_0} \cos (s + \varphi)  - \frac{A^3}{L_0^2} \cos^3 (s + \varphi) </math>
|}

On réutilise la formule trigonométrique :

{| align="center" border="0"
|<math> \cos^3 \alpha = \frac{1}{4} \ \left[ 3 \cos \alpha + \cos (3 \alpha ) \right]</math>
|}

d'où l'équation différentielle pour la fonction {{math|''y''{{ind|1}}(''s'')}} :

{| align="center" border="0"
|<math>\frac{\mathrm{d}^2y_1(s)}{\mathrm{d}s^2} + y_1(s) = A\left[ \frac{2 \omega_1}{\omega_0} - \frac{3 A^2}{4L_0^2} \right] \cos (s + \varphi) - \frac{A^3}{4L_0^2} \cos (3s + 3 \varphi ) </math>
|}

Il suffit alors d'annuler le coefficient devant le terme en {{math|cos (''s'' + φ)}} en posant :

<center>
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-
| bgcolor="#fff8ff" |<math>\frac{2 \omega_1}{\omega_0} - \frac{3 A^2}{4L_0^2} = 0 \quad \Longrightarrow \quad \omega_1 = \frac{3 A^2}{8L_0^2} \omega_0 </math>
|}</center>

On obtient alors l'équation différentielle finale pour la fonction {{math|''y''{{ind|1}}(''s'')}} :

{| align="center" border="0"
|<math>\frac{\mathrm{d}^2y_1(s)}{\mathrm{d}s^2} +  y_1(s) = - \frac{A^3}{4L_0^2} \cos (3s + 3 \varphi ) </math>
|}
<!---
Cette équation différentielle admet une solution générale ''bornée'' sans terme séculaire :
--->

== Perturbation singulière ==

{{Article détaillé|Perturbation singulière}}

Il arrive que le terme dépendant d'un paramètre de petite valeur ne puisse pas être supprimé en fixant le paramètre à 0 sans changer la nature du problème et avoir des conséquences sur les propriétés des solutions. On parle alors de perturbation singulière. Généralement, les perturbations singulières apparaissent quand le paramètre multiplie le terme de plus haut degré (dérivée d'ordre le plus élevé, plus forte puissance...).

== Problèmes multi-échelles ==
{{…}}

== Annexes ==
=== Bibliographie ===
==== Bibliothèque virtuelle ====
* Nils Berglund ; ''Perturbation theory of dynamical systems'' (2001). Disponible sur l'ArXiv : [http://fr.arxiv.org/abs/math?papernum=0111178 math.HO/0111178].

==== Ouvrages de référence ====
* Ali H. Nayfeh ; ''Perturbation Methods'', [[John Wiley & Sons]] (New York-1973), réédité dans la collection : Wiley Classics Library (2000), {{ISBN|0-471-39917-5}}.
* Ali H. Nayfeh ; ''Introduction to Perturbation Techniques'', John Wiley & Sons (New York-1981), {{ISBN|0-471-31013-1}}.
* E. John Hinch ; ''Perturbation Methods'', Cambridge Texts in Applied Mathematics, [[Cambridge University Press]] (1991), {{ISBN|0521378974}}.
* Donald R. Smith ; ''Singular-Perturbation Theory: An Introduction with Applications'', Cambridge University Press (1985), {{ISBN|0-521-30042-8}}.

==== Articles de revue ====
* Donald R. Smith ; ''The Multivariable Method in Singular Perturbation Analysis'', SIAM Review '''17 (2)''' (1975), 221-273.
* A. B. Vasilieva ; ''On the Development of Singular Perturbation Theory at Moscow State University and Elsewhere'', SIAM Review '''36 (3)''' (1994), 440-452.

==== Applications à la mécanique céleste ====
La [[mécanique hamiltonienne]] a connu une avancée spectaculaire en 1954 avec l'avènement de la « [[Théorème KAM|théorie KAM]] ». Les références ci-dessous sont donc classées en fonction de cet évènement.

===== Théorie pré-KAM =====
* Pierre-Simon Laplace, ''Traité de mécanique céleste'', [[Éditions Jacques Gabay]], 1990. Réédition d'un ouvrage classique de la fin du {{s-|XIX|e}}, en 4 volumes, [http://sites.mathdoc.fr/cgi-bin/oeitem?id=OE_LAPLACE__4_1_0 disponible] sur ''[[Gallica]]''.{{Commentaire biblio SRL|Niveau second cycle universitaire.}}
* François-Félix Tisserand, ''Traité de mécanique céleste'', Éditions Jacques Gabay, 1990. Réédition d'un ouvrage classique de la fin du {{s-|XIX|e}}, en 4 volumes, disponible : [https://gallica.bnf.fr/ark:/12148/bpt6k29031t tome 1], [https://gallica.bnf.fr/ark:/12148/bpt6k290325 tome 2], [http://visualiseur.bnf.fr/ark:/12148/bpt6k29033h tome 3], sur ''Gallica'', [https://archive.org/details/traitdemcani04tissuoft tome 4] sur [[Internet Archive]].{{Commentaire biblio SRL|Niveau second cycle universitaire.}}
* Henri Poincaré, ''Leçons de mécanique céleste'', 3 tomes, 1905-1910, réédité par Jacques Gabay, Paris, 2003, disponible en fac-similé sur ''Gallica'' : [https://gallica.bnf.fr/ark:/12148/bpt6k95010b.table ''Tome I''], ''Tome II'' [https://gallica.bnf.fr/ark:/12148/bpt6k95011p partie 1] et [https://gallica.bnf.fr/ark:/12148/bpt6k6227628j/f5.image partie 2], [https://gallica.bnf.fr/ark:/12148/bpt6k62331958 ''Tome III''].{{Commentaire biblio SRL|Une somme de référence, par le grand mathématicien qui a tant contribué au sujet. Niveau second cycle universitaire.}}
* Anders Lindstedt, ''Abh. K. Akad. Wiss. St. Petersburg'', vol. 31, {{n°|4}}, 1882{{référence incomplète|date=janvier 2016}}.
* Anders Lindstedt, « Sur la forme des expressions des distances mutuelles dans le problème des trois corps », ''[[Comptes rendus de l'Académie des sciences]]'', vol. 97, 1883, {{p.|1276}} et 1353, {{Lire sur Wikisource|Fichier:Comptes_rendus_hebdomadaires_des_séances_de_l’Académie_des_sciences,_tome_097,_1883.djvu}}.
* Henri Poincaré, « Sur les séries de M. Lindstedt », ''Comptes rendus de l'Académie des sciences'', vol. 108, 1889, {{p.|21-24}}, [https://gallica.bnf.fr/ark:/12148/bpt6k3064g/f21n4.capture disponible] sur ''Gallica''.
* Henri Poincaré, « Sur l'application de la méthode de M. Lindstedt au problème des trois corps », ''Comptes rendus de l'Académie des sciences'', vol. 114, 1892, {{p.|1305-1309}}, [https://gallica.bnf.fr/ark:/12148/bpt6k3070h/f1305n5.capture disponible] sur ''Gallica''.

===== Théorie post-KAM =====
* Florin Diacu & Philip Holmes ; ''Celestial Encounters - The Origin of Chaos & Stability'', [[Princeton University Press]] (1996).].
* V.I. Arnold, V.V. Kozlov & A.I. Neishtadt ; ''Mathematical Aspects of Classical & Celestial Mechanics'', Springer-Verlag ({{2e}} édition-1997)
* Bibliographie de l'article [[Mécanique céleste]]

=== Articles connexes ===
* [[Notations de Landau]]
* [[Sommabilité de Borel]]
* [[Théorie de la perturbation (mécanique quantique)]]

=== Références ===

<references>
<ref name="Poincaré">Henri Poincaré ; ''Leçons de mécanique céleste'', 3 tomes, (1905-1910), réédité par Jacques Gabay, Paris (2003). Une somme de référence, par le grand mathématicien qui a tant contribué au sujet. Niveau second cycle universitaire. Cet ouvrage est disponible en fac-similé sur Gallica : [https://gallica.bnf.fr/document?O=N095010 ''Tome I''], [https://gallica.bnf.fr/document?O=N095011 ''Tome II''], [https://gallica.bnf.fr/document?O=N095012 ''Tome III''].</ref>
</references>

{{Palette|Équations différentielles}}
{{Portail|astronomie|physique|mathématiques}}

{{DEFAULTSORT:Theorie des perturbations}}
[[Catégorie:Équation]]
[[Catégorie:Méthode mathématique de la physique]]
[[Catégorie:Mécanique céleste]]
[[Catégorie:Mécanique quantique]]