Équation d'Euler-Lagrange

Équation d’Euler-Lagrande en application

L’équation d'Euler-Lagrange (en anglais, Euler–Lagrange equation ou ELE) est un résultat mathématique qui joue un rôle fondamental dans le calcul des variations. On retrouve cette équation dans de nombreux problèmes réels de minimisation de longueur d'arc, tels que le problème brachistochrone ou bien encore les problèmes géodésiques. Elle est nommée d'après Leonhard Euler et Joseph-Louis Lagrange.

Notations

E désignera un espace vectoriel normé, [t₀ , t₁] un intervalle réel, et ${\mathcal {G}}$ l'espace affine des fonctions x : [t₀ , t₁] → E de classe C¹ telles que $x\left(t_{i}\right)=x_{i}$, où x₀ , x₁ sont deux vecteurs fixés de E.

Le vecteur dérivé d'une fonction $x\in {\mathcal {G}}$ en un point t ∈ [t₀ , t₁] est noté ${\dot {x}}(t)$.

On se donne par ailleurs une fonction ${\mathcal {L}}:\left[t_{0},t_{1}\right]\times E^{2}\to \mathbb {R}$ de classe C¹.

Ses trois variables étant notées $t,x,{\dot {x}}$ (ce qui risque de prêter à confusion avec la notation précédente mais est d'un usage courant), ses trois applications différentielles partielles sont notées

• ${\frac {\partial {\mathcal {L}}}{\partial t}}$ (de $\mathbb {R} \times E^{2}$ dans $\mathbb {R}$) et
• ${\frac {\partial {\mathcal {L}}}{\partial x}},{\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}$ (de $\mathbb {R} \times E^{2}$ dans E', le dual de E).

Lorsqu'on les compose par la fonction $\left[t_{0},t_{1}\right]\to \mathbb {R} \times E^{2},\;t\mapsto \left(t,x(t),{\dot {x}}(t)\right)$ pour une fonction donnée $x\in {\mathcal {G}}$, on obtient trois fonctions définies sur [t₀ , t₁] (encore à valeurs respectivement dans $\mathbb {R}$, E' et E'), que l'on note usuellement de la même façon (bien que, à nouveau, cela prête à confusion), ce qui donne en particulier un sens aux deux fonctions

${\frac {\partial {\mathcal {L}}}{\partial x}}{\text{ et }}{\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}:\left[t_{0},t_{1}\right]\to E'$.

Énoncé

Soit J la fonctionnelle définie sur ${\mathcal {G}}$ par :

$J(x)=\int _{t_{0}}^{t_{1}}{\mathcal {L}}\left(t,x(t),{\dot {x}}(t)\right)\,\mathrm {d} t$.

Pour toute fonction $x\in {\mathcal {G}}$ stationnaire pour J, $t\mapsto {\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}$ est dérivable et

${\frac {\partial {\mathcal {L}}}{\partial x}}-{\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}\right)=0$.

Démonstration partielle

La démonstration qui suit est annoncée comme « partielle » parce qu'elle suppose que ${\dot {x}}$ et ${\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}$ sont de classe C¹ (auquel cas la dérivabilité de $t\mapsto {\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}\left(t,x(t),{\dot {x}}(t)\right)$ est assurée d'emblée). Pour une démonstration supposant seulement que $x$ et ${\mathcal {L}}$ sont de classe C¹, voir l'application du lemme de Du Bois-Reymond au calcul des variations.

L'expression « stationnaire », dans l'énoncé, signifie : vérifiant la condition d'Euler, qui est une condition nécessaire pour que la fonction $x$ rende extrémale la fonctionnelle $J$ (restreinte dans cette preuve aux fonctions de ${\mathcal {G}}$ de classe C²).

Cette condition d'Euler s'écrit : ${\frac {\mathrm {d} J(x+\varepsilon h)}{\mathrm {d} \varepsilon }}_{|\varepsilon =0}=0$, pour toute fonction h : [t₀ , t₁] → E (de classe C²) nulle en t₀ et t₁. Or

${\frac {\mathrm {d} J(x+\varepsilon h)}{\mathrm {d} \varepsilon }}_{|\varepsilon =0}=\int _{t_{0}}^{t_{1}}\left(\left\langle {\frac {\partial {\mathcal {L}}}{\partial x}},h\right\rangle +\left\langle {\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}},{\dot {h}}\right\rangle \right)\,\mathrm {d} t$ (où $\langle ~,~\rangle :E'\times E\to \mathbb {R}$ est le crochet de dualité)

et le second terme de l'intégrale s'exprime, grâce à une intégration par parties (permise par les hypothèses supplémentaires de régularité), sous la forme

$\int _{t_{0}}^{t_{1}}\left\langle {\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}},{\dot {h}}\right\rangle \,\mathrm {d} t=\left[\left\langle {\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}},h\right\rangle \right]_{t_{0}}^{t_{1}}-\int _{t_{0}}^{t_{1}}\left\langle {\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}\right),h\right\rangle \,\mathrm {d} t$.

Le crochet étant nul puisque h(t₀) = h(t₁) = 0, la condition d'Euler s'écrit donc :

$0={\frac {\mathrm {d} J(x+\varepsilon h)}{\mathrm {d} \varepsilon }}_{|\varepsilon =0}=\int _{t_{0}}^{t_{1}}\left\langle {\frac {\partial {\mathcal {L}}}{\partial x}}-{\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}\right),h\right\rangle \,\mathrm {d} t$.

En appliquant le lemme fondamental du calcul des variations, on en déduit :

${\frac {\partial {\mathcal {L}}}{\partial x}}-{\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}\right)=0$.

Exemple

Un exemple est une application du principe de Fermat. L'objectif est de déterminer un chemin optique plan, dont les coordonnées sont notées horizontalement t et verticalement x, pour respecter les notations de l'énoncé ci-dessus. Le rayon lumineux traverse le vide, à l'exception de la zone correspondant aux valeurs de t situées entre –1 et 1. Sur cette bande, on suppose que l'indice n_t n'est plus égal à 1 mais à 1/|t|. Entre les deux bandes, le chemin optique a pour longueur :

$\mathrm {L} =\int _{-1}^{1}f\left(t,x(t),{\dot {x}}(t)\right)\,\mathrm {d} t\quad {\text{avec}}\quad f(t,x,y)=n_{t}{\sqrt {1+y^{2}}}={\frac {\sqrt {1+y^{2}}}{|t|}}$.

Puisqu'ici ${\frac {\partial f}{\partial x}}=0$, l'équation d'Euler-Lagrange stipule que la dérivée partielle de f par rapport à sa troisième variable est une constante, notée ici C, si elle est appliquée aux variables t, x et sa dérivée. On obtient :

$C={\frac {\partial f}{\partial {\dot {x}}}}={\frac {\dot {x}}{|t|{\sqrt {1+{\dot {x}}^{2}}}}}\quad {\text{donc}}\quad {\dot {x}}^{2}=C^{2}t^{2}(1+{\dot {x}}^{2})$.

Ce résultat s'écrit encore, en posant u = C|t| :

${\dot {x}}={\frac {u}{\sqrt {1-u^{2}}}}$.

On reconnait l'équation d'une portion de cycloïde.

Identité de Beltrami

Un cas particulier fréquent est celui où la fonction ${\mathcal {L}}$ est indépendante de t. Un corollaire de l'équation d'Euler-Lagrange est alors l'identité de Beltrami :

${\mathcal {L}}-{\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}{\dot {x}}=C$.

La lettre C désigne une constante réelle, qui est aussi la transformée de Legendre de la fonction ${\mathcal {L}}$ par rapport à la variable ${\textstyle {\dot {x}}}$.

Démonstration

En supposant $x$ deux fois dérivable, dérivons le membre de gauche de l'identité de Beltrami :

${\frac {\mathrm {d} }{\mathrm {d} t}}\left({\mathcal {L}}-{\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}{\dot {x}}\right)={\frac {\partial {\mathcal {L}}}{\partial x}}{\dot {x}}+{\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}{\ddot {x}}-\left({\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}\right){\dot {x}}+{\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}{\ddot {x}}\right)=\left({\frac {\partial {\mathcal {L}}}{\partial x}}-{\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\partial {\mathcal {L}}}{\partial {\dot {x}}}}\right)\right){\dot {x}}=0$

Un exemple historique célèbre est la courbe brachistochrone. La question posée revient à trouver la courbe reliant un point A à un point B, situé à une altitude plus faible, tel qu'un point matériel partant du point A sans vitesse initiale et glissant sans frottement sur la courbe rejoigne le plus rapidement possible le point B.

Lorsque ${\mathcal {L}}$ est une fonction homogène de la variable ${\dot {x}}$, le théorème d'Euler appliqué à l'identité de Beltrami implique ${\mathcal {L}}=C$.

Liens externes

• (en) Eric W. Weisstein, « Euler-Lagrange Differential Equation », sur MathWorld