Théorie des perturbations 528333 31441118 2008-07-09T16:16:06Z Jean-Luc W 53231 /* Convergence de la série perturbative ? */ lien D'un point de vue [[heuristique]], la '''théorie des perturbations''' est une méthode [[Mathématiques|mathématique]] générale qui permet de trouver une ''solution approchée'' d'une équation mathématique <math>(E_{\lambda})</math> dépendante d'un paramètre <math>\lambda</math> lorsque la solution de l'équation <math>(E_{0})</math>, correspondant à la valeur <math>\lambda = 0</math>, est connue exactement. L'équation mathématique <math>(E_{\lambda})</math> peut être une équation algébrique, une équation différentielle, une équation aux valeurs propres, ... La méthode consiste à chercher la solution approchée de l'équation <math>(E_{\lambda})</math> sous la forme d'un développement en [[série (mathématiques)|série]] des puissances du paramètre <math>\lambda</math>, cette solution approchée étant supposé être une approximation d'autant meilleure de la solution exacte, mais inconnue, que la valeur absolue du paramètre <math>\lambda</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 <math> N</math> [[problème à N corps|corps]] en interaction [[gravitation|gravitationnelle]] 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 \ge 3</math>. Cet aspect de la théorie des perturbations a été synthétisé à la fin du {{XIXe siècle}} dans les ouvrages classiques de 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 XIXe siècle, en 4 volumes. Niveau second cycle universitaire. Cet ouvrage est disponible en fac-similé sur [http://gallica.bnf.fr/document?O=N077586 ''Gallica''].</ref>, 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 XIXe siècle, en 4 volumes. Niveau second cycle universitaire. Cet ouvrage est disponible en fac-similé sur [http://gallica.bnf.fr/document?O=N029031 ''Gallica''].</ref> et Poincaré<ref>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 : [http://gallica.bnf.fr/document?O=N095010 ''Tome I''], [http://gallica.bnf.fr/document?O=N095011 ''Tome II''], [http://gallica.bnf.fr/document?O=N095012 ''Tome III''].</ref>, avant de connaitre 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 : [[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_{\lambda})</math> sous la forme d'un développement en série des puissances du paramètre <math>\lambda</math> ; la question de la ''convergence'' de cette série se pose alors. Ce problème a été réglé pour l'astronomie par [[Henri Poincaré|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 : <blockquote> « ''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.'' </blockquote> <blockquote> ''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 999 999 ème ; mais ils regarderont la seconde comme divergente, parce que le terme général peut croître au delà de toute limite.'' </blockquote> <blockquote> ''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.'' </blockquote> <blockquote> ''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.'' » </blockquote> 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> : <blockquote> « 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. » </blockquote> 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 donnent 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 du premier ordre suivante : <center> {| border=1 cellpadding=10 style="border-collapse:collapse" |----- | bgcolor="#fff8ff" |<math> (E_{\lambda}) \ : \qquad \frac{dx(t)}{dt} \ + \ \frac{1}{\tau} \ x(t) \ + \ \lambda \ \frac{x^2(t)}{\tau \, L_0} \ = \ 0 </math> |}</center> Dans cette équation, <math>t</math> représente le temps, <math>\tau</math> un paramètre fixé [[analyse dimensionnelle|homogène]] à un temps, <math>L_0</math> un paramètre fixé homogène à une longueur, et <math>\lambda</math> le paramètre de perturbation, [[analyse dimensionnelle|sans dimensions]]. On cherche à déterminer la fonction <math>x(t)</math> inconnue, homogène à une longueur, et vérifiant la ''condition initiale'' : à l'instant <math>t = 0</math>, on a : <math>x(0) = X_0</math>. === Théorie des perturbations au premier ordre === Le problème de départ de la théorie des perturbations est l'équation différentielle <math>(E_{0})</math> correspondant à la valeur <math>\lambda = 0</math> : {| align="center" border="0" |<math> (E_{0}) \ : \qquad \frac{dx(t)}{dt} \ + \ \frac{1}{\tau} \ x(t) \ = \ 0 </math> |} dont la solution analytique exacte est bien connue : {| align="center" border="0" |<math> x_{0}(t) \ = \ A \ e^{- \, t/\tau}</math> |} où <math>A</math> est une constante, pour l'instant inconnue. Illustrons la méthode de perturbations en nous limitant pour simplifier au ''premier ordre'' dans le développement en série des puissances du paramètre <math>\lambda</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_{1}(t)</math> est une fonction inconnue, à déterminer. On injecte cette expression dans l'équation différentielle exacte <math>(E_{\lambda})</math>. ''En se limitant aux termes du premier ordre inclus'' et en utilisant le fait que <math>x_{0}(t)</math> est la solution exacte de <math>(E_0)</math>, on obtient la solution physique approchée au premier ordre : <center> {| border=1 cellpadding=10 style="border-collapse:collapse" |----- | bgcolor="#fff8ff" |<math> x_{\lambda}(t) \ = \ X_0 \ \left( \, 1 \, - \, \lambda \, \frac{X_0}{L_0} \, \right) \ e^{- \, t/\tau} \ + \ \lambda \ \frac{X_0^2}{L_0} \ e^{- \, 2t/\tau} \ + \ O(\lambda^2)</math> |}</center> {{boîte déroulante|align=left|titre=Démonstration|contenu= On injecte le développement dans l'équation différentielle exacte <math>(E_{\lambda})</math>. En se limitant au termes du premier ordre inclus et en utilisant le fait que <math>x_{0}(t)</math> est la solution exacte de <math>(E_0)</math>, on obtient l'équation suivante pour la fonction <math>x_{1}(t)</math> : <center> <math>\lambda \ \left[ \ \frac{dx_1(t)}{dt} \ + \ \frac{1}{\tau} \ x_1(t) \ \right] \ + \ \lambda \ \frac{x^2_0(t)}{\tau \, L_0} \ = \ 0 </math> </center> On a omis le terme <math>O(\lambda^2)</math> de reste de Landau pour simplifier l'écriture. Cette équation se réécrit explicitement : <center> <math> \frac{dx_1(t)}{dt} \ + \ \frac{1}{\tau} \ x_1(t) \ = \ - \ \frac{A^2}{\tau \, L_0} \ e^{- \, 2t/\tau}</math> </center> Cette équation différentielle admet une solution analytique exacte de la forme : <center> <math> x_{1}(t) \ = \ B \ e^{- \, 2t/\tau}</math> </center> où <math>B</math> est une constante, qu'on détermine en introduisant l'expression de <math>x_1(t)</math> dans l'équation différentielle. On obtient explicitement : <center> <math> B \ = \ \frac{A^2}{L_0} \quad \Longrightarrow \quad x_{1}(t) \ = \ \frac{A^2}{L_0} \ e^{- \, 2t/\tau}</math> </center> d'où la solution générale de <math>(E_{\lambda})</math> approchée au premier ordre. }} === Comparaison avec la solution exacte === On peut démontrer ici que l'équation différentielle <math> (E_{\lambda})</math> vérifiant la condition initiale : <math>x(0) = X_0</math> admet pour toute les valeurs du paramètre <math>\lambda</math> la solution exacte suivante : <center> {| border=1 cellpadding=10 style="border-collapse:collapse" |----- | bgcolor="#fff8ff" |<math> x_{\lambda}(t) \ = \ \frac{X_0 \ e^{- \, t/\tau}}{1 \ + \ \lambda \ \frac{X_0}{L_0} \ \left( \, 1 \, - \, \ e^{- \, t/\tau} \ \right)}</math> |}</center> Un développement limité de cette expression au premier ordre en <math>\lambda</math> donne explicitement la solution approchée déterminée au paragraphe précédent par la théorie des perturbations au premier ordre : {| align="center" border="0" |<math> x_{\lambda}(t) \ = \ X_0 \ \left( \, 1 \, - \, \lambda \, \frac{X_0}{L_0} \, \right) \ e^{- \, t/\tau} \ + \ \lambda \ \frac{X_0^2}{L_0} \ e^{- \, 2t/\tau} \ + \ O(\lambda^2)</math> |} 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>\lambda</math> allant de 0.1 à 0.5, en prenant : <math>X_0 = L_0 = 1</math> m, <math>\tau = 1</math> s. * en bleu, la solution exacte. * en rouge, la solution approchée au premier ordre. [[Image:Perturb1.jpg]] [[Image:Perturb2.jpg]] [[Image:Perturb3.jpg]] [[Image:Perturb4.jpg]] [[Image:Perturb5.jpg]] == 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> {| border=1 cellpadding=10 style="border-collapse:collapse" |----- | bgcolor="#fff8ff" |<math> (E_{\lambda}) \ : \qquad \frac{d^2x(t)}{dt^2} \ + \ \omega_0^2 \ x(t) \ + \ \lambda \ \frac{\omega_0^2}{L_0^2} \ x^3(t) \ = \ 0 </math> |}</center> Dans cette équation, <math>t</math> représente le temps, <math>\omega_0</math> un paramètre fixé [[analyse dimensionnelle|homogène]] à une pulsation, c’est-à-dire l'inverse d'un temps. <math>L_0</math> est un paramètre fixé homogène à une longueur, et <math>\lambda</math> le paramètre de perturbation, [[analyse dimensionnelle|sans dimensions]]. On cherche à déterminer la fonction <math>x(t)</math> inconnue, homogène à une longueur, et vérifiant les ''conditions initiales'' : à l'instant <math>t = 0</math>, on a : <math>x(0) = X_0</math> 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 <math>m</math> soumise à une force dérivant d'une [[énergie potentielle]] <math>V(x)</math> : {| align="center" border="0" |<math> m \ \frac{d^2x(t)}{dt^2} \ = \ - \ \frac{dV(x(t))}{dx}</math> |} où le potentiel <math>V(x)</math> quartique s'écrit : <center> {| border=1 cellpadding=10 style="border-collapse:collapse" |----- | bgcolor="#fff8ff" |<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>\lambda</math> ''positives'' ou nulles, <math>V(x)</math> représente un ''puits'' de potentiel. La conservation de l'énergie mécanique totale <math>E</math> de la particule : {| align="center" border="0" |<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> |} entraine alors que le mouvement est ''borné'' dans un intervalle <math>[ \, x_1, \, x_2 \, ]</math>, où les points tournants <math>x_1</math> et <math>x_2</math> sont les deux solutions réelles de l'équation : {| align="center" border="0" |<math>E \ = \ V(x)</math> |} ==== 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_{0})</math> correspondant à la valeur <math>\lambda = 0</math> : {| align="center" border="0" |<math> (E_{0}) \ : \qquad \frac{d^2x(t)}{dt^2} \ + \ \omega_0^2 \ x(t) \ = \ 0 </math> |} Cette équation est par définition un [[oscillateur harmonique]] de pulsation <math>\omega_0</math>, dont la solution analytique exacte est bien connue : {| align="center" border="0" |<math> x_{0}(t) \ = \ A \ \cos \left( \, \omega_0 \, t \, + \, \varphi \, \right) </math> |} où <math>A</math> et <math>\varphi</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 : {| align="center" border="0" |<math> x_{\lambda}(t) \ = \ x_{0}(t) \ + \ \lambda \ x_{1}(t) \ + \ O(\lambda^2)</math> |} où <math> x_{1}(t)</math> est une fonction inconnue, à déterminer. On injecte cette expression dans l'équation différentielle exacte <math>(E_{\lambda})</math>. ''En se limitant au termes du premier ordre inclus'' et en utilisant le fait que <math>x_{0}(t)</math> est la solution exacte de <math>(E_0)</math>, on obtient l'expression au premier ordre de la théorie de perturbation : <center> {| border=1 cellpadding=10 style="border-collapse:collapse" |----- | bgcolor="#fff8ff" |<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> {{boîte déroulante|align=left|titre=Démonstration|contenu= On injecte le développement dans l'équation différentielle exacte <math>(E_{\lambda})</math>. ''En se limitant au termes du premier ordre inclus'' et en utilisant le fait que <math>x_{0}(t)</math> est la solution exacte de <math>(E_0)</math>, on obtient l'équation suivante pour la fonction <math> x_{1}(t)</math> : <center> <math> \frac{d^2x_1(t)}{dt^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_{1}(t)</math> : <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</math>, on a : <math>x(0) = X_0</math> 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 \varphi = 0 </math>, donc que : <math>\varphi = 0 </math>. 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 pertubation contient un terme ''proportionnel au temps'' : {| align="center" border="0" |<math>\lambda \ \omega_0 t \ \sin \left( \omega_0 t \right)</math> |} 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> \lambda</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_{\lambda}(t)</math> 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êchent 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 équation différentielles, permet d'éliminer ces termes séculaire<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 [http://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 [http://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 [http://gallica.bnf.fr/document?O=N003070 ''Gallica''].</ref>{{,}}<ref>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 : [http://gallica.bnf.fr/document?O=N095010 ''Tome I''], [http://gallica.bnf.fr/document?O=N095011 ''Tome II''], [http://gallica.bnf.fr/document?O=N095012 ''Tome III''].</ref>. 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[ \, \left( 1 + \lambda \right) \, \omega_0 t \, \right]</math> |} Cette expression clairement ''bornée'' développée au premier ordre en <math>\lambda</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 apparait 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} \ = \ \left( 1 + \lambda \right) \, \omega_0 </math> |} qui est légèrement différente de la pulsation initiale <math>\omega_0</math> 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{d^2x(t)}{dt^2} \ + \ \omega_0^2 \ x(t) \ = \ \lambda \ f \left(x(t), \, \frac{dx(t)}{dt}, \, t, \, \lambda \right) </math> |}</center> où <math>f</math> est une fonction ''paire'' de <math>x</math> et ''impaire'' de <math>\dot{x}</math> qui est de plus soit ''périodique'' en <math>t</math>, soit indépendante de <math>t</math>. La méthode consiste à faire un ''changement d'échelle de temps'' en introduisant une nouvelle variable <math>s</math> [[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>\omega_k, \ k \ge 1</math> devront être choisies afin de faire disparaitre les termes séculaires de la série perturbative 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{d^2x(t)}{dt^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 \to s</math> et définissons la nouvelle fonction inconnue <math>y</math> de la variable <math>s</math>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 chaine donne pour la dérivée première : {| align="center" border="0" |<math>\frac{dx}{dt} \ = \ \frac{dy}{ds} \ \frac{ds}{dt} </math> |} et pour la dérivée seconde : {| align="center" border="0" |<math>\frac{d^2x}{dt^2} \ = \ \frac{d~~}{ds} \left( \frac{dy}{ds} \ \frac{ds}{dt} \right) \ \frac{ds}{dt} \ = \ \frac{d^2y}{ds^2} \ \left( \frac{ds}{dt} \right)^2 </math> |} Comme on a au premier ordre : {| align="center" border="0" |<math>\frac{ds}{dt} \ = \ \omega_0 \ + \ \lambda \omega_1 \ + \ O(\lambda^2)</math> |} on obtient pour les dérivées : {| align="center" border="0" |<math>\frac{dx}{dt} \ = \ ( \omega_0 \ + \ \lambda \omega_1 ) \ \frac{dy}{ds} \ + \ O(\lambda^2) </math> |} {| align="center" border="0" |<math>\frac{d^2x}{dt^2} \ = \ \left( \omega_0^2 \ + \ 2 \lambda \omega_0 \omega_1 \right) \ \frac{d^2y}{ds^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{d^2y(s)}{ds^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>\lambda</math> : {| align="center" border="0" |<math> \omega_0^2 \left(\frac{d^2y_0}{ds^2} + y_{0} \right) \ + \ \lambda \left[ \omega_0^2 \frac{d^2y_1}{ds^2} + 2 \omega_0 \omega_1 \frac{d^2y_0}{ds^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érentielle : {| align="center" border="0" |<math> \frac{d^2y_0}{ds^2} + y_{0} \ = \ 0</math> |} {| align="center" border="0" |<math>\omega_0^2 \frac{d^2y_1}{ds^2} + 2 \omega_0 \omega_1 \frac{d^2y_0}{ds^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_0(s) = A \, \cos (s + \varphi)</math> où <math>A</math> et <math>\varphi</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_1(s)</math> : {| align="center" border="0" |<math>\frac{d^2y_1(s)}{ds^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_{1}(s)</math> : {| align="center" border="0" |<math>\frac{d^2y_1(s)}{ds^2} \ + \ y_1(s) \ = \ \left[ \ \frac{2 \, \omega_1}{\omega_0} \ - \ \frac{3 \, A^3}{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 + \varphi) </math> en posant : <center> {| border=1 cellpadding=10 style="border-collapse:collapse" |----- | bgcolor="#fff8ff" |<math>\frac{2 \, \omega_1}{\omega_0} \ - \ \frac{3 \, A^3}{4L_0^2} \ = \ 0 \quad \Longrightarrow \quad \omega_1 \ = \ \frac{3 \, A^3}{4L_0^2} \, \omega_0 </math> |}</center> On obtient alors l'équation différentielle finale pour la fonction <math> y_{1}(s)</math> : {| align="center" border="0" |<math>\frac{d^2y_1(s)}{ds^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 : A suivre ... == Perturbation singulière == {{ébauche|mathématiques|physique}} == 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'', Editions Jacques Gabay (1990). Réédition d'un ouvrage classique de la fin du XIXe siècle, en 4 volumes. Niveau second cycle universitaire. Cet ouvrage est disponible en fac-similé sur [http://gallica.bnf.fr/document?O=N077586 ''Gallica''] {{Gallica}}. * 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 XIXe siècle, en 4 volumes. Niveau second cycle universitaire. Cet ouvrage est disponible en fac-similé sur [http://gallica.bnf.fr/document?O=N029031 ''Gallica'']. * 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 : [http://gallica.bnf.fr/document?O=N095010 ''Tome I''], [http://gallica.bnf.fr/document?O=N095011 ''Tome II''], [http://gallica.bnf.fr/document?O=N095012 ''Tome III'']. * Anders Lindstedt ; Abh. K. Akad. Wiss. St. Petersburg 31 (1882), 4. * 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 [http://gallica.bnf.fr/document?O=N003053 ''Gallica'']. * 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 [http://gallica.bnf.fr/document?O=N003064 ''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 114 (1892) 1305-1309. Cette note est disponible en fac-similé sur [http://gallica.bnf.fr/document?O=N003070 ''Gallica'']. ===== Théorie post-KAM ===== * Florin Diacu & Philip Holmes ; ''Celestial Encounters - The Origin of Chaos & Stability'', Princeton University Press (1996). Voir également le site internet : [http://www.math.uvic.ca/faculty/diacu/ ''Diacu'']. * V.I. Arnold, V.V. Kozlov & A.I. Neishtadt ; ''Mathematical Aspects of Classical & Celestial Mechanics'', Springer-Verlag (2{{e}} é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/> {{Portail|astronomie|physique|mathématiques}} [[Catégorie:Équation]] [[Catégorie:Méthode mathématique de la physique]] [[Catégorie:Mécanique céleste]] [[Catégorie:Mécanique quantique]]