Ecuación de Euler-Lagrange
La ecuación de Euler-Lagrange (en español, ecuación de Euler-Lagrange o ELE ) es un resultado matemático que juega un papel fundamental en el cálculo de las variaciones . Esta ecuación se encuentra en muchos problemas reales de minimización de la longitud del arco , como el problema de la braquistocrona o incluso los problemas geodésicos . Lleva el nombre de Leonhard Euler y Joseph-Louis Lagrange .
Notaciones
E denotará un espacio vectorial normalizado , [ t 0 , t 1 ] un intervalo real , y el espacio afín de funciones x : [ t 0 , t 1 ] → E de clase C 1 tal que , donde x 0 , x 1 son dos vectores conjunto de E .
GRAMO{\ Displaystyle {\ mathcal {G}}}
X(tI)=XI{\ Displaystyle x \ left (t_ {i} \ right) = x_ {i}}![{\ Displaystyle x \ left (t_ {i} \ right) = x_ {i}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/033e0c2e1a593df0a190c889b2be800906e6b181)
Se denota el vector derivado de una función en un punto t ∈ [ t 0 , t 1 ] .
X∈GRAMO{\ Displaystyle x \ in {\ mathcal {G}}}
X˙(t){\ Displaystyle {\ dot {x}} (t)}![{\ Displaystyle {\ dot {x}} (t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/867f9767dea7667cf8deac86ec6bbe3162c8d4b1)
También nos damos una función de clase C 1 .
L:[t0,t1]×mi2→R{\ Displaystyle {\ mathcal {L}}: \ left [t_ {0}, t_ {1} \ right] \ times E ^ {2} \ to \ mathbb {R}}![{\ Displaystyle {\ mathcal {L}}: \ left [t_ {0}, t_ {1} \ right] \ times E ^ {2} \ to \ mathbb {R}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92b307f9ce03c0f76f8744aea7e3dcc683ff34ed)
Se anotan sus tres variables (lo que probablemente dé lugar a confusión con la notación anterior pero es de uso común), se anotan
sus tres aplicaciones diferenciales parcialest,X,X˙{\ Displaystyle t, x, {\ dot {x}}}
-
∂L∂t{\ Displaystyle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial t}}}
(desde adentro ) yR×mi2{\ Displaystyle \ mathbb {R} \ times E ^ {2}}
R{\ Displaystyle \ mathbb {R}}![\ mathbb {R}](https://wikimedia.org/api/rest_v1/media/math/render/svg/786849c765da7a84dbc3cce43e96aad58a5868dc)
-
∂L∂X,∂L∂X˙{\ estilo de visualización {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}}, {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ punto {x}}}} }
(de en E ' , el dual de E ).R×mi2{\ Displaystyle \ mathbb {R} \ times E ^ {2}}![{\ Displaystyle \ mathbb {R} \ times E ^ {2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/88bbe5815a8bab05cfffd74c4c8b3da011f35763)
Cuando las componimos por la función para una función dada , obtenemos tres funciones definidas en [ t 0 , t 1 ] (nuevamente con valores respectivamente en , E ' y E' ), que usualmente denotamos de la misma manera ( aunque, de nuevo, esto es confuso), lo que en particular da significado a las dos funciones
[t0,t1]→R×mi2,t↦(t,X(t),X˙(t)){\ Displaystyle \ left [t_ {0}, t_ {1} \ right] \ to \ mathbb {R} \ times E ^ {2}, \; t \ mapsto \ left (t, x (t), {\ punto {x}} (t) \ right)}
X∈GRAMO{\ Displaystyle x \ in {\ mathcal {G}}}
R{\ Displaystyle \ mathbb {R}}![\ mathbb {R}](https://wikimedia.org/api/rest_v1/media/math/render/svg/786849c765da7a84dbc3cce43e96aad58a5868dc)
∂L∂X y ∂L∂X˙:[t0,t1]→mi′{\ estilo de visualización {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} {\ texto {et}} {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ punto {x}}}}: \ left [t_ {0}, t_ {1} \ right] \ to E '}![{\ estilo de visualización {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} {\ texto {et}} {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ punto {x}}}}: \ left [t_ {0}, t_ {1} \ right] \ to E '}](https://wikimedia.org/api/rest_v1/media/math/render/svg/41e56fd0ac377c538adadf53659b7a443690d4e3)
.
Estados
Sea J el funcional definido por:
GRAMO{\ Displaystyle {\ mathcal {G}}}![{\ mathcal {G}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8a980c59d42c003fd07fdf3646e1fb95ff82f99)
J(X)=∫t0t1L(t,X(t),X˙(t))Dt{\ Displaystyle J (x) = \ int _ {t_ {0}} ^ {t_ {1}} {\ mathcal {L}} \ left (t, x (t), {\ dot {x}} (t ) \ derecha) \, \ mathrm {d} t}![{\ Displaystyle J (x) = \ int _ {t_ {0}} ^ {t_ {1}} {\ mathcal {L}} \ left (t, x (t), {\ dot {x}} (t ) \ derecha) \, \ mathrm {d} t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff876876b0d4344c4b8d2a36a669b3f25e73d25d)
.
Para cualquier función estacionaria para J , es diferenciable y
X∈GRAMO{\ Displaystyle x \ in {\ mathcal {G}}}
t↦∂L∂X˙{\ Displaystyle t \ mapsto {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}}}![{\ Displaystyle t \ mapsto {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7326b12673354f51530b7462aab7186c523b87f7)
∂L∂X-DDt(∂L∂X˙)=0{\ Displaystyle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} - {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ left ({\ frac { \ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} \ right) = 0}![{\ Displaystyle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} - {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ left ({\ frac { \ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} \ right) = 0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/327511c45e8add369b4035a644a5ecbfd95f7f37)
.
Demostración parcial
La prueba que sigue se anuncia como "parcial" porque supone que y son de clase C 1 (en cuyo caso la diferenciabilidad de está asegurada desde el principio). Para una demostración asumiendo solo que y son de clase C 1 , vea la aplicación del lema de Du Bois-Reymond al cálculo de variaciones .
X˙{\ Displaystyle {\ dot {x}}}
∂L∂X˙{\ Displaystyle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}}}
t↦∂L∂X˙(t,X(t),X˙(t)){\ Displaystyle t \ mapsto {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} \ left (t, x (t), {\ dot {x}} ( t) \ derecha)}
X{\ Displaystyle x}
L{\ Displaystyle {\ mathcal {L}}}![{\ mathcal {L}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9027196ecb178d598958555ea01c43157d83597c)
La expresión "estacionario", en el enunciado, significa: satisfacer la condición de Euler , que es una condición necesaria para que la función haga el funcional extremal (restringido en esta prueba a funciones de clase C 2 ).
X{\ Displaystyle x}
J{\ Displaystyle J}
GRAMO{\ Displaystyle {\ mathcal {G}}}![{\ mathcal {G}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8a980c59d42c003fd07fdf3646e1fb95ff82f99)
Esta condición euleriana se escribe :, para cualquier función h : [ t 0 , t 1 ] → E (de clase C 2 ) cero en t 0 y t 1 . Oro
DJ(X+εh)Dε|ε=0=0{\ displaystyle {\ frac {\ mathrm {d} J (x + \ varepsilon h)} {\ mathrm {d} \ varepsilon}} _ {| \ varepsilon = 0} = 0}![{\ displaystyle {\ frac {\ mathrm {d} J (x + \ varepsilon h)} {\ mathrm {d} \ varepsilon}} _ {| \ varepsilon = 0} = 0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b33414c72e2c8327919b49ea068835f6e4c6340c)
DJ(X+εh)Dε|ε=0=∫t0t1(⟨∂L∂X,h⟩+⟨∂L∂X˙,h˙⟩)Dt{\ displaystyle {\ frac {\ mathrm {d} J (x + \ varepsilon h)} {\ mathrm {d} \ varepsilon}} _ {| \ varepsilon = 0} = \ int _ {t_ {0}} ^ {t_ {1}} \ izquierda (\ izquierda \ langle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}}, h \ derecha \ rangle + \ izquierda \ langle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}}, {\ dot {h}} \ right \ rangle \ right) \, \ mathrm {d} t}![{\ displaystyle {\ frac {\ mathrm {d} J (x + \ varepsilon h)} {\ mathrm {d} \ varepsilon}} _ {| \ varepsilon = 0} = \ int _ {t_ {0}} ^ {t_ {1}} \ izquierda (\ izquierda \ langle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}}, h \ derecha \ rangle + \ izquierda \ langle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}}, {\ dot {h}} \ right \ rangle \ right) \, \ mathrm {d} t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f6d73f80851f13c208654b0ed843d6bbbfa5464)
(donde está el
corchete de dualidad )
⟨ , ⟩:mi′×mi→R{\ Displaystyle \ langle ~, ~ \ rangle: E '\ times E \ to \ mathbb {R}}![{\ Displaystyle \ langle ~, ~ \ rangle: E '\ times E \ to \ mathbb {R}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b191e3a67ed2a627a5ea698efd1c487bef0cc30d)
y el segundo término de la integral se expresa, gracias a una integración por partes (permitida por los supuestos adicionales de regularidad), en la forma
∫t0t1⟨∂L∂X˙,h˙⟩Dt=[⟨∂L∂X˙,h⟩]t0t1-∫t0t1⟨DDt(∂L∂X˙),h⟩Dt{\ Displaystyle \ 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 {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x} }}}, h \ right \ rangle \ right] _ {t_ {0}} ^ {t_ {1}} - \ int _ {t_ {0}} ^ {t_ {1}} \ left \ langle {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ izquierda ({\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ punto {x}}}} \ derecha), h \ right \ rangle \, \ mathrm {d} t}![{\ Displaystyle \ 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 {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x} }}}, h \ right \ rangle \ right] _ {t_ {0}} ^ {t_ {1}} - \ int _ {t_ {0}} ^ {t_ {1}} \ left \ langle {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ izquierda ({\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ punto {x}}}} \ derecha), h \ right \ rangle \, \ mathrm {d} t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/23f0808f35ddf5b941b66330efeb91d4976ecb73)
.
El gancho es cero ya que h ( t 0 ) = h ( t 1 ) = 0 , la condición de Euler se escribe por lo tanto:
0=DJ(X+εh)Dε|ε=0=∫t0t1⟨∂L∂X-DDt(∂L∂X˙),h⟩Dt{\ displaystyle 0 = {\ frac {\ mathrm {d} J (x + \ varepsilon h)} {\ mathrm {d} \ varepsilon}} _ {| \ varepsilon = 0} = \ int _ {t_ {0} } ^ {t_ {1}} \ izquierda \ langle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} - {\ frac {\ mathrm {d}} {\ mathrm {d} t }} \ izquierda ({\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} \ derecha), h \ derecha \ rangle \, \ mathrm {d} t}![{\ displaystyle 0 = {\ frac {\ mathrm {d} J (x + \ varepsilon h)} {\ mathrm {d} \ varepsilon}} _ {| \ varepsilon = 0} = \ int _ {t_ {0} } ^ {t_ {1}} \ izquierda \ langle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} - {\ frac {\ mathrm {d}} {\ mathrm {d} t }} \ izquierda ({\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} \ derecha), h \ derecha \ rangle \, \ mathrm {d} t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/45490e94e3c406d200d3fd616c672545544e79a8)
.
Aplicando el lema fundamental del cálculo de variaciones , deducimos:
∂L∂X-DDt(∂L∂X˙)=0{\ Displaystyle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} - {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ left ({\ frac { \ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} \ right) = 0}![{\ Displaystyle {\ frac {\ parcial {\ mathcal {L}}} {\ parcial x}} - {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ left ({\ frac { \ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} \ right) = 0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/327511c45e8add369b4035a644a5ecbfd95f7f37)
.
Ejemplo
Un ejemplo es una aplicación del principio de Fermat . El objetivo es determinar una trayectoria óptica plana, cuyas coordenadas se indican horizontalmente ty verticalmente x , para cumplir con las notaciones del enunciado anterior. El rayo de luz atraviesa el vacío, a excepción de la zona correspondiente a los valores de t situados entre –1 y 1. En esta banda, se supone que el índice n t ya no es igual a 1 sino a 1 / | t |. Entre las dos bandas, el camino óptico tiene la longitud :
L=∫-11F(t,X(t),X˙(t))DtconF(t,X,y)=not1+y2=1+y2|t|{\ Displaystyle \ mathrm {L} = \ int _ {- 1} ^ {1} f \ left (t, x (t), {\ dot {x}} (t) \ right) \, \ mathrm {d } t \ quad {\ text {con}} \ quad f (t, x, y) = n_ {t} {\ sqrt {1 + y ^ {2}}} = {\ frac {\ sqrt {1 + y ^ {2}}} {| t |}}}![{\ Displaystyle \ mathrm {L} = \ int _ {- 1} ^ {1} f \ left (t, x (t), {\ dot {x}} (t) \ right) \, \ mathrm {d } t \ quad {\ text {con}} \ quad f (t, x, y) = n_ {t} {\ sqrt {1 + y ^ {2}}} = {\ frac {\ sqrt {1 + y ^ {2}}} {| t |}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/026474d494d893b46e63f1d072f37747829a0e69)
.
Ya que aquí , la ecuación de Euler-Lagrange establece que la derivada parcial de f con respecto a su tercera variable es una constante, indicada aquí C , si se aplica a las variables t , x y su derivada. Obtenemos :
∂F∂X=0{\ Displaystyle {\ frac {\ parcial f} {\ parcial x}} = 0}![{\ Displaystyle {\ frac {\ parcial f} {\ parcial x}} = 0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8888fba457b0ecc3d0614b8826635536d39eecfc)
VS=∂F∂X˙=X˙|t|1+X˙2EntoncesX˙2=VS2t2(1+X˙2){\ Displaystyle C = {\ frac {\ parcial f} {\ parcial {\ dot {x}}}} = {\ frac {\ dot {x}} {| t | {\ sqrt {1 + {\ dot { x}} ^ {2}}}}} \ quad {\ text {por lo tanto}} \ quad {\ dot {x}} ^ {2} = C ^ {2} t ^ {2} (1 + {\ dot {x}} ^ {2})}![{\ Displaystyle C = {\ frac {\ parcial f} {\ parcial {\ dot {x}}}} = {\ frac {\ dot {x}} {| t | {\ sqrt {1 + {\ dot { x}} ^ {2}}}}} \ quad {\ text {por lo tanto}} \ quad {\ dot {x}} ^ {2} = C ^ {2} t ^ {2} (1 + {\ dot {x}} ^ {2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7241d49e944257e273d81f38a97b2c45790d1636)
.
Este resultado se escribe de nuevo, estableciendo u = C | t | :
X˙=tu1-tu2{\ Displaystyle {\ dot {x}} = {\ frac {u} {\ sqrt {1-u ^ {2}}}}}![\ dot x = \ frac {u} {\ sqrt {1 - u ^ 2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83ba141948f5be2904f9b525c5ba804bd29f4790)
.
Reconocemos la ecuación de una parte de una cicloide .
Identidad Beltrami
Un caso especial frecuente es aquel en el que la función es independiente de t . Un corolario de la ecuación de Euler-Lagrange es entonces la identidad de Beltrami :
L{\ Displaystyle {\ mathcal {L}}}![{\ mathcal {L}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9027196ecb178d598958555ea01c43157d83597c)
L-∂L∂X˙X˙=VS{\ Displaystyle {\ mathcal {L}} - {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} {\ dot {x}} = C}![{\ Displaystyle {\ mathcal {L}} - {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} {\ dot {x}} = C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19ebcef8bfbcf3875774e2123a62038f9d05ac8e)
.
La letra C designa una constante real, que también es la transformada de Legendre de la función f con respecto a la variable .
X˙{\ textstyle {\ dot {x}}}![{\ textstyle \ dot {x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/539bc058f8b19747a57099fdbaf556ad99239fb5)
Demostración
Suponiendo que sean dos veces diferenciables, derivemos el lado izquierdo de la identidad de Beltrami:
X{\ Displaystyle x}![X](https://wikimedia.org/api/rest_v1/media/math/render/svg/87f9e315fd7e2ba406057a97300593c4802b53e4)
DDt(L-∂L∂X˙X˙)=∂L∂XX˙+∂L∂X˙X¨-(DDt(∂L∂X˙)X˙+∂L∂X˙X¨)=(∂L∂X-DDt(∂L∂X˙))X˙=0{\ Displaystyle {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ left ({\ mathcal {L}} - {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} {\ dot {x}} \ right) = {\ frac {\ partial {\ mathcal {L}}} {\ partial x}} {\ dot {x}} + {\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x}}}} {\ ddot {x}} - \ left ({\ frac {\ mathrm {d}} {\ mathrm {d } t}} \ izquierda ({\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ punto {x}}}} \ derecha) {\ punto {x}} + {\ frac {\ parcial {\ mathcal {L}}} {\ partial {\ dot {x}}}} {\ ddot {x}} \ right) = \ left ({\ frac {\ partial {\ mathcal {L}}} {\ parcial x}} - {\ frac {\ mathrm {d}} {\ mathrm {d} t}} \ izquierda ({\ frac {\ parcial {\ mathcal {L}}} {\ parcial {\ dot {x} }}} \ right) \ right) {\ dot {x}} = 0}
Un ejemplo histórico famoso es la curva braquistócrona . La pregunta planteada equivale a encontrar la curva que conecta un punto A con un punto B, ubicado a menor altitud, como un punto material que parte del punto A sin velocidad inicial y deslizándose sin fricción en la curva une lo más rápidamente posible el punto B .
Cuando es una función homogénea de la variable , el teorema de Euler aplicado a la identidad de Beltrami implica .
L{\ Displaystyle {\ mathcal {L}}}
X˙{\ Displaystyle {\ dot {x}}}
L=VS{\ Displaystyle {\ mathcal {L}} = C}![{\ Displaystyle {\ mathcal {L}} = C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5624b3acc5b380ffb5718c9ef55621309b086a36)
enlaces externos
<img src="https://fr.wikipedia.org/wiki/Special:CentralAutoLogin/start?type=1x1" alt="" title="" width="1" height="1" style="border: none; position: absolute;">