Erhalten Sie die Lagrange-Funktion aus dem System der gekoppelten Gleichung [geschlossen]

In diesem speziellen Papier

"Interaktion zwischen einem sich bewegenden Spiegel und Strahlungsdruck: Eine Hamiltonsche Formulierung" von CKLaw, PhysRevA.51.2537

(2.6) Q ¨ k = ω k 2 Q k + 2 Q ˙ Q J G k J Q ˙ J + Q ¨ Q Q ˙ 2 Q 2 J G k J Q J + Q ˙ 2 Q 2 J G J k G J Q

(2.7) M Q ¨ = v ( Q ) Q + 1 Q k , J ( 1 ) k + J ω k ω J Q k Q J
Wo k , J , N + { 1 , 2 , 3 , } .

Hier die positionsabhängigen Frequenzen ω k werden von gegeben

(2.8) ω k ( Q ) = k π Q
und die dimensionslosen Koeffizienten G k J werden von gegeben

(2.9-2.10) G k J = { ( 1 ) k + J 2 k J J 2 k 2 k J 0 k = J

(3.1) L ( Q , Q ˙ , Q k , Q ˙ k ) = 1 2 k [ Q ˙ k 2 ω k 2 ( Q ) Q k 2 ] + 1 2 M Q ˙ 2 v ( Q ) Q ˙ Q J , k G k J Q ˙ k Q J + Q ˙ 2 2 Q 2 J , k , G k J G k Q Q J

wie man die in Gleichung 3.1 angegebene Lagrange-Funktion aus dem System der gekoppelten Gleichungen 2.6 und 2.7 erhält

Das Problem, mit dem ich konfrontiert bin, besteht darin, den kanonischen Impuls in einer solchen Gleichung zu identifizieren und auch nicht in Euler-Lagrange-Form formulieren zu können, um den Lagrange zu erhalten.

Die kanonischen Impulse P k , P konjugieren zu Q k , Q sind in der Veröffentlichung jeweils durch die folgenden Gleichungen gegeben

(3.3) P k = Q ˙ k Q ˙ Q J G k J Q J (3.4) P = M Q ˙ 1 Q J k G k J P k Q J

Ich habe den Back-Step-Prozess ausprobiert, bei dem ich versucht hatte, das System von eqns zu erhalten. aus der gegebenen Lagrange-Funktion ohne Erfolg. Vielleicht brauche ich einen anderen Ansatz

Ich bin sehr dankbar für jede Art von Hilfe in dieser Angelegenheit.

Die Bezeichnungen k, j, l scheinen sich über das Ganze zu ziehen N . Haben Sie in den Gleichungen (2.9-2.10). G k J = 0 für k = ± J oder sind k , J 0 ?

Antworten (1)

HAUPTTEIL: Der Lagrange

Lassen Sie die Bewegungsgleichungen und die Euler-Lagrange-Gleichungen mit null rechten Seiten ausdrücken

(01a) Q ¨ k + ω k 2 Q k 2 Q ˙ Q J G k J Q ˙ J Q ¨ Q Q ˙ 2 Q 2 J G k J Q J Q ˙ 2 Q 2 J G J k G J Q = 0

(01b) M Q ¨ + v ( Q ) Q 1 Q k , J ( 1 ) k + J ω k ω J Q k Q J = 0

(02a) D D T ( L Q ˙ k ) L Q k = 0

(02b) D D T ( L Q ˙ ) L Q = 0

Wo L ( Q , Q ˙ , Q k , Q ˙ k ) der Lagrange.

Wir gehen zu den folgenden Definitionen über, um die große Menge an Variablen und Indizes mittels komprimierter vereinfachter Ausdrücke zu handhaben:

(03) Q def [ Q 1 Q 2 Q k ] Q ˙ def [ Q ˙ 1 Q ˙ 2 Q ˙ k ] Q ¨ def [ Q ¨ 1 Q ¨ 2 Q ¨ k ]

(04) G def [ 0 G 12 G 13 G 1 k G 21 0 G 23 G 2 k G k 1 G k 2 G k 3 0 ] = G T

(05) Ω ( Q ) def [ ω 1 0 0 0 ω 2 0 0 0 ω k ] = π Q [ 1 0 0 0 2 0 0 0 k ] = Ω T ( Q )

(06) ϕ ( Q , Q ˙ ) def Q ˙ Q

Wir definieren unten auch den reellen Skalar, so etwas wie das innere Produkt von reellen Vektoren

(07) < Q , P > def k Q k P k

Unter diesen Definitionen und unter Verwendung der Gleichungen (A-01), siehe HILFSABSCHNITT, haben wir die folgenden Ausdrücke (08) anstelle der Bewegungsgleichungen (01) und (09) anstelle von (02):

(08a) Q ¨ + Ω 2 ( Q ) Q 2 ϕ ( Q , Q ˙ ) G Q ˙ ϕ ˙ ( Q , Q ˙ ) G Q + ϕ 2 ( Q , Q ˙ ) G 2 Q = 0

(08b) M Q ¨ + v ( Q ) Q 1 Q < Ω 2 ( Q ) Q , Q > + 1 Q < Ω 2 ( Q ) Q , G Q > = 0
(09a) D D T ( L Q ˙ ) L Q = 0
(09b) D D T ( L Q ˙ ) L Q = 0
während die Lagrange-Funktion des Systems, siehe Gleichung (3.1) in der Frage, unter Verwendung der Gleichungen (A-02) ausgedrückt wird als
L ( Q , Q ˙ , Q , Q ˙ ) = 1 2 < Q ˙ , Q ˙ > 1 2 < Ω 2 ( Q ) Q , Q > + 1 2 M Q ˙ 2 v ( Q ) ϕ < G Q , Q ˙ > 1 2 ϕ 2 < G 2 Q , Q >
(10) -------------------------------------------------- ---------
und in noch kompakterer Form
(10 ' ) L ( Q , Q ˙ , Q , Q ˙ ) = 1 2 ( ϕ G Q Q ˙ 2 Ω Q 2 ) + 1 2 M Q ˙ 2 v
Wir werden versuchen, die Lagrange-Funktion Schritt für Schritt durch ein Trial-and-Error-Verfahren zu erstellen.

Wir erwarten also, dass der 1. Term von Gleichung (08a) aus einem Lagrange-Teil stammt L 1 ( Q ˙ ) so dass durch (09a)

(11) D D T ( L 1 Q ˙ ) = Q ¨ L 1 Q ˙ = Q ˙
Aus der Regel (A-3.d), siehe HILFSABSCHNITT, L 1 Ist
(12) L 1 ( Q ˙ ) = 1 2 < Q ˙ , Q ˙ >
seit
(13) ( < Q ˙ , Q ˙ > ) Q ˙ = 2 Q ˙

Für den 2. Term von Gleichung (08a) erwarten wir einen Lagrange-Anteil L 2 ( Q , Q ) so dass durch (09a)

(14) L 2 Q = Ω 2 ( Q ) Q
So
(15) L 2 ( Q , Q ) = 1 2 < Ω 2 ( Q ) Q , Q >
da aus der Regel (A-3.c) und der symmetrischen (genauer: diagonalen) Matrix Ω 2 = ( Ω 2 ) T
(16) ( < Ω 2 Q , Q > ) Q = [ Ω 2 + ( Ω 2 ) T ] Q = 2 Ω 2 Q
Aber als Lagrange-Teil L 2 ( Q , Q ) ist eine Funktion von Q außerdem erzeugt es Elemente in den Bewegungsgleichungen, wenn es in den 2. Term von (09b) eingefügt wird:
(17) L 2 Q = + 1 2 ( < Ω 2 Q , Q > ) Q = + < Ω Ω Q Q , Q > = 1 Q < Ω 2 ( Q ) Q , Q >
das ist genau der 3. Term in Gleichung (08b).

Andererseits sind die ersten beiden Terme von (08b) die eines Teilchens, das sich in einem Potential bewegt, also stammen sie aus einem Lagrange-Teil L 3 ( Q , Q ˙ ) :

(18) L 3 ( Q , Q ˙ ) = 1 2 M Q ˙ 2 v ( Q )
Dieser Teil L 3 ( Q , Q ˙ ) wenn in (9a) eingesetzt, ergibt sich nichts (kein Term in Bewegungsgleichungen). Nun, in (08a) geben die Hälfte des 3. Terms und des 4. Terms an
(19) ϕ ( Q , Q ˙ ) G Q ˙ ϕ ˙ ( Q , Q ˙ ) G Q = D D T ( ϕ G Q )
wir erwarten also einen Lagrange-Anteil L 4 ( Q , Q ˙ , Q , Q ˙ ) so dass durch (09a)
(20) L 4 Q ˙ = ϕ ( Q , Q ˙ ) G Q
das ist
(21) L 4 ( Q , Q ˙ , Q , Q ˙ ) = ϕ ( Q , Q ˙ ) < G Q , Q ˙ >
Aber wegen der Antisymmetrie von G , dieser Teil kann auch ausgedrückt werden als
(22) L 4 ( Q , Q ˙ , Q , Q ˙ ) = + ϕ ( Q , Q ˙ ) < G Q ˙ , Q >
also einfügen in den 2. Term von (09a)
(23) L 4 Q = ϕ ( Q , Q ˙ ) G Q ˙
das ist die andere Hälfte des 3. Terms in (08a). Das bedeutet, dass L 4 , wenn in (09a) eingefügt, erzeugt den 3. und 4. Term von (08a)
(24) D D T ( L 4 Q ˙ ) L 4 Q = 2 ϕ ( Q , Q ˙ ) G Q ˙ ϕ ˙ ( Q , Q ˙ ) G Q
Die Ausgabe der Einfügung von L 4 in (09b) würde später zusammen mit untersucht L 5 . Der 5. Term von (08a) kann aus einem Lagrange-Teil stammen L 5 ( Q , Q ˙ , Q , Q ˙ ) so dass durch (09a)
(25) L 5 Q = ϕ 2 ( Q , Q ˙ ) G 2 Q
So
(26) L 5 ( Q , Q ˙ , Q , Q ˙ ) = 1 2 ϕ 2 < G 2 Q , Q >
da von (A-03.c) und die Symmetrie von G 2
(27) ( < G 2 Q , Q > ) Q = ( G 2 + ( G 2 ) T ) Q = 2 G 2 Q
Es kann bewiesen werden, siehe A PROOF ABSCHNITT, dass die Summe L 45 = L 4 + L 5
(28) L 45 ( Q , Q ˙ , Q , Q ˙ ) = L 4 + L 5 = 1 2 ϕ 2 < G 2 Q , Q > ϕ ( Q , Q ˙ ) < G Q , Q ˙ >
wenn in (09b) eingefügt, entsteht der 4. Term von (08b)
(29) D D T ( L 45 Q ˙ ) L 45 Q = + 1 Q < Ω 2 ( Q ) Q , G Q >

In Gleichung (30) unten summieren wir die gefundenen Lagrange-Teile und die endgültige Lagrange-Funktion ist

L ( Q , Q ˙ , Q , Q ˙ ) = 1 2 < Q ˙ , Q ˙ > L 1 1 2 < Ω 2 ( Q ) Q , Q > L 2 + 1 2 M Q ˙ 2 v ( Q ) L 3 ϕ < G Q , Q ˙ > L 4 1 2 ϕ 2 < G 2 Q , Q > L 5
(30) -------------------------------------------------- ---------
identisch mit der in der Veröffentlichung angegebenen Gleichung (10).

Gleichungen (31) sind die Bewegungsgleichungen (08) mit geschweiften Klammern unter Punkten, die von den Lagrange-Termen angegeben sind L M Diese Artikel stammen von:

(31a) Q ¨ L 1 + Ω 2 ( Q ) Q L 2 2 ϕ ( Q , Q ˙ ) G Q ˙ ϕ ˙ ( Q , Q ˙ ) G Q L 4 + ϕ 2 ( Q , Q ˙ ) G 2 Q L 5 = 0
(31b) M Q ¨ + v ( Q ) Q L 3 1 Q < Ω 2 ( Q ) Q , Q > L 2 + 1 Q < Ω 2 ( Q ) Q , G Q > L 4 + L 5 = 0

Beachten Sie, dass die kanonischen Impulse P , P konjugieren zu Q , Q jeweils sind
(32a) P = L Q ˙ = Q ˙ Q ˙ Q G Q (32b) P = L Q ˙ = M Q ˙ 1 Q < G Q , P >
wo für den Beweis von (32b)
P = L Q ˙ = M Q ˙ 1 Q < G Q , Q ˙ > Q ˙ Q 2 < G 2 Q , Q > = M Q ˙ 1 Q < G Q , Q ˙ > + Q ˙ Q 2 < G Q , G Q > (32b ' ) = M Q ˙ 1 Q < G Q , Q ˙ Q ˙ Q G Q P > = M Q ˙ 1 Q < G Q , P >
Die Gleichungen (32a) und (32b) sind identisch mit (3.3) bzw. (3.4) der unten angegebenen Veröffentlichung
(3.3) P k = Q ˙ k Q ˙ Q J G k J Q J (3.4) P = M Q ˙ 1 Q J k G k J P k Q J


HILFSTEIL : Komprimierte vereinfachte Ausdrücke und partielle Differenzierungsregeln

Gleichungen (A-01) sind nützlich für die Umwandlung der Bewegungsgleichungen von der Form (01) in die Form (08) :

(A-01.a) ω k 2 Q k = [ Ω 2 ( Q ) Q ] k (A-01.b) J G k J Q J = [ G Q ] k (A-01.c) J G k J Q ˙ J = [ G Q ˙ ] k (A-01.d) J G J k G J Q = ( J G k J G J ) Q = ( G 2 ) k Q = ( G 2 Q ) k (A-01.e) Q ¨ Q Q ˙ 2 Q 2 = D D T ( Q ˙ Q ) = D ϕ ( Q , Q ˙ ) D T = ϕ ˙ ( Q , Q ˙ ) (A-01.f) k , J ( 1 ) k + J ω k ω J Q k Q J = < Ω 2 ( Q ) Q , Q > < Ω 2 ( Q ) Q , G Q >
Der Beweis von (A-01.f) läuft wie folgt ab
(A-01.f ' ) k , J ( 1 ) k + J ω k ω J Q k Q J = k ω k 2 Q k 2 < Ω 2 ( Q ) Q , Q > + k , J k ( 1 ) k + J ω k ω J Q k Q J < Ω 2 ( Q ) Q , G Q >
seit
k , J k ( 1 ) k + J ω k ω J Q k Q J = ( π Q ) 2 k , J k ( 1 ) k + J k J Q k Q J = ( π Q ) 2 k , J k ( 1 ) k + J 2 k J J 2 k 2 G k J J 2 k 2 2 Q k Q J = 1 2 k , J G k J ( ω J 2 ω k 2 ) Q k Q J = 1 2 J ( ω J 2 Q J ) [ Ω 2 ( Q ) Q ] J k G J k Q k [ G Q ] J 1 2 k ( ω k 2 Q k ) [ Ω 2 ( Q ) Q ] k J G k J Q J [ G Q ] k = < Ω 2 ( Q ) Q , G Q >
(A-01.f ' ' ) -------------------------------------------------- ---------

Die Gleichungen (A-02) und (A-03) sind nützlich für die Umwandlung der Lagrange-Funktion von der Form (3.1), siehe fragliche Gleichung, in die Form (10) und für die schrittweise Konstruktion dieser Lagrange-Funktion aus der Bewegungsgleichungen (08) :

(A-02.a) k Q ˙ k 2 = < Q ˙ , Q ˙ > = Q ˙ 2 (A-02.b) k ω k 2 ( Q ) Q k 2 = < Ω 2 Q , Q > = < Ω Q , Ω T Q > = < Ω Q , Ω Q > = Ω Q 2 (A-02.c) J , k G k J Q ˙ k Q J = < G Q , Q ˙ > = < G Q ˙ , Q > (A-02.d) J , k , G k J G k Q Q J = < G 2 Q , Q > = < G Q , G Q > = G Q 2
Die Gleichungen (A-02.c) und (A-02.d) werden jeweils wie folgt bewiesen
J , k G k J Q ˙ k Q J = k ( J G k J Q J ) Q ˙ k = k [ G Q ] k [ Q ˙ ] k = (A-02.c ' ) < G Q , Q ˙ > = < Q , G T Q ˙ > = < Q , G Q ˙ > = < G Q ˙ , Q >
J , k , G k J G k Q Q J = k ( J G k J Q J ) ( G k Q ) = k [ G Q ] k [ G Q ] k (A-02.d ' ) = < G Q , G Q > = < G T G Q , Q > = < G 2 Q , Q >

Die nachstehenden Gleichungen (A-03) sind gewissermaßen partielle Differentiationsregeln einer Skalarfunktion einer Vektorvariablen S in Bezug auf diese Variable. Die Skalarfunktionen sind normalerweise Skalarprodukte und der variable Vektor ist es S = Q oder Q ˙ . Im Folgenden A , R sind Vektoren und F lineare Transformation alle unabhängig vom variablen Vektor S . Normalerweise F = Ω , Ω 2 , G , G 2 :

(A-03.a) ( < A , S > ) S = ( < S , A > ) S = A (A-03.b) ( < R , F S > ) S = ( < F T R , S > ) S = F T R (A-03.c) ( < F S , S > ) S = ( F + F T ) S (A-03.d) ( < S , S > ) S = 2 S
(A-03.b) ist ein Sonderfall von (A-03.a) mit A = F T R und (A-03.d) ist ein Spezialfall von (A-03.c) mit F = ICH .

Eine im folgenden Abschnitt nützliche Identität ist

(A-04) G T = G < G S , S > = 0 , für jeden reellen Vektor  S
seit
(A-04 ' ) < G S , S > = < S , G T S > = < S , ( G ) S > = < G S , S >


A BEWEISABSCHNITT: Beweis von Gleichung (29) bei gegebener Gleichung (28).

Wir werden Gleichung (29) aus (28) beweisen, wobei die beiden Gleichungen hier der Einfachheit halber wiederholt werden

(29) D D T ( L 45 Q ˙ ) L 45 Q = + 1 Q < Ω 2 ( Q ) Q , G Q >
Wo
(28) L 45 ( Q , Q ˙ , Q , Q ˙ ) def 1 2 ϕ 2 ( Q , Q ˙ ) < G 2 Q , Q > ϕ ( Q , Q ˙ ) < G Q , Q ˙ >

L 45 Q = ϕ ϕ Q < G 2 Q , Q > + ϕ Q < G Q , Q ˙ > = ( Q ˙ Q ) ( Q ˙ Q ) Q < G 2 Q , Q > + ( Q ˙ Q ) Q < G Q , Q ˙ >
So
(B-01) L 45 Q = ( Q ˙ 2 Q 3 ) < G 2 Q , Q > + ( Q ˙ Q 2 ) < G Q , Q ˙ >
Jetzt
L 45 Q ˙ = ϕ ϕ Q ˙ < G 2 Q , Q > ϕ Q ˙ < G Q , Q ˙ > (B-02) L 45 Q ˙ = ( Q ˙ Q 2 ) < G 2 Q , Q > + ( 1 Q ) < G Q , Q ˙ >
Differenzieren (B-02) bzgl T
D D T ( L 45 Q ˙ ) = ( Q ¨ Q 2 Q ˙ 2 Q 3 ) < G 2 Q , Q > + ( Q ˙ Q 2 ) < G 2 Q ˙ , Q > + ( Q ˙ Q 2 ) < G 2 Q , Q ˙ > + ( Q ˙ Q 2 ) < G Q , Q ˙ > + ( 1 Q ) < G Q ˙ , Q ˙ > = 0 , siehe (A-04) (B-03) + ( 1 Q ) < G Q , Q ¨ >
Hinzufügen von (B-01) und (B-03)
D D T ( L 45 Q ˙ ) L 45 Q = ( Q ¨ Q Q ˙ 2 Q 3 ) < G 2 Q , Q > + ( 2 Q ˙ Q 2 ) < G 2 Q , Q ˙ > + ( 1 Q ) < G Q , Q ¨ > = + 1 Q < ( Q ¨ Q ˙ 2 Q 2 ) G Q + 2 Q ˙ Q G Q ˙ Q ¨ , G Q > = + 1 Q < ϕ ˙ G Q + 2 ϕ G Q ˙ Q ¨ = Ω 2 ( Q ) Q + ϕ 2 G 2 Q ,  siehe (08a) , G Q > + 1 Q < Ω 2 ( Q ) Q + ϕ 2 G 2 Q , G Q > = + 1 Q < Ω 2 ( Q ) Q , G Q > + ϕ 2 Q < G 2 Q , G Q > = 0 , siehe (A-04)
So
(B-04) D D T ( L 45 Q ˙ ) L 45 Q = + 1 Q < Ω 2 ( Q ) Q , G Q >
QED.