ОСНОВНОЙ РАЗДЕЛ: Лагранжиан
Выразим уравнения движения и уравнения Эйлера-Лагранжа с нулевыми правыми частями
Вопрос¨к+ю2кВопроск− 2д˙д∑Джгк ДжВопрос˙Дж−д¨д−д˙2д2∑Джгк ДжВопросДж−д˙2д2∑Дж ℓгj kгДж ℓВопросℓ= 0(01а)
мд¨+∂В( q)∂д−1д∑к , дж( − 1)к + джюкюДжВопроскВопросДж= 0(01б)
ггт(∂л∂Вопрос˙к) —∂л∂Вопроск= 0(02а)
ггт(∂л∂д˙) —∂л∂д= 0(02б)
гдел ( q,д˙,Вопроск,Вопрос˙к)
лагранжиан.
Мы переходим к следующим определениям, чтобы обрабатывать большое количество переменных и индексов с помощью сжатых упрощенных выражений:
Вопрос≡деф⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢Вопрос1Вопрос2⋮Вопроск⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥Вопрос˙≡деф⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢Вопрос˙1Вопрос˙2⋮Вопрос˙к⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥Вопрос¨≡деф⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢Вопрос¨1Вопрос¨2⋮Вопрос¨к⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥(03)
г≡деф⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢0г21⋮гк 1⋮г120⋮гк 2⋮г13г23⋮гк 3⋮⋯⋯⋮⋯⋮г1 тыс.г2 к⋮0⋮⋯⋯⋮⋯⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥= -гТ(04)
Ом ( q)≡деф⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢ю10⋮0⋮0ю2⋮0⋮⋯⋯⋮⋯⋮00⋮юк⋮⋯⋯⋮⋯⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥"="πд⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢10⋮0⋮02⋮0⋮⋯⋯⋮⋯⋮00⋮к⋮⋯⋯⋮⋯⋮⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥"="ОмТ( q)(05)
ϕ ( q,д˙)≡дефд˙д(06)
Мы также определяем действительный скаляр ниже, что-то вроде внутреннего произведения действительных векторов.
< Q , Р >≡деф∑кВопроскпк(07)
При этих определениях и использовании уравнений (A-01), см. ВСПОМОГАТЕЛЬНЫЙ РАЗДЕЛ, мы имеем следующие выражения (08) вместо уравнений движения (01) и (09) вместо (02):
Вопрос¨+Ом2( q) Q − 2 ϕ ( q,д˙) ГВопрос˙−ф˙( q,д˙) Г К +ф2( q,д˙)г2Q = 0(08а)
мд¨+∂В( q)∂д−1д<Ом2( q) Q , Q > +1д<Ом2( q) Q , G Q > = 0(08б)
ггт(∂л∂Вопрос˙) —∂л∂Вопрос= 0(09а)
ггт(∂л∂д˙) —∂л∂д= 0(09б)
в то время как лагранжиан системы, см. уравнение (3.1) в вопросе, с использованием уравнений (A-02) выражается как
л ( q,д˙, Q ,Вопрос˙) =12<Вопрос˙,Вопрос˙> −12<Ом2( q) Q , Q > +12мд˙2− В( q) − ϕ < G Q ,Вопрос˙> −12ф2<г2Q , Q >
-------------------------------------------------- ---------(10)
и в еще более компактном виде
л ( q,д˙, Q ,Вопрос˙) =12(∥∥ϕ G Q -Вопрос˙∥∥2−∥ Ом Q ∥2) +12мд˙2− В(10′)
Попробуем построить лагранжиан шаг за шагом методом проб и ошибок.
Итак, мы ожидаем, что 1-й член уравнения (08a) будет исходить из лагранжевой частил1(Вопрос˙)
такое, что в силу (09а)
ггт(∂л1∂Вопрос˙) =Вопрос¨⟹∂л1∂Вопрос˙"="Вопрос˙(11)
Из правила (A-3.d) см. ВСПОМОГАТЕЛЬНЫЙ РАЗДЕЛ,
л1
является
л1(Вопрос˙) =12<Вопрос˙,Вопрос˙>(12)
с
∂( <Вопрос˙,Вопрос˙> )∂Вопрос˙= 2Вопрос˙(13)
Для 2-го члена уравнения (08а) мы ожидаем, что лагранжева частьл2( q, В )
такое, что в силу (09а)
−∂л2∂Вопрос"="Ом2( q) В(14)
так
л2( q, Q ) = -12<Ом2( q) Q , Q >(15)
поскольку из правила (A-3.c) и симметричной (точнее: диагональной) матрицы
Ом2"="(Ом2)Т
∂( <Ом2Q , Q > )∂Вопрос= [Ом2+(Ом2)Т] Q =2Ом2Вопрос(16)
Но так как лагранжева часть
л2( q, В )
является функцией
д
кроме того, он производит элементы в уравнениях движения, если вставить их во 2-й член (09b):
−∂л2∂д= +12∂( <Ом2Q , Q > )∂Вопрос= + < Ом∂Ом∂дQ , Q >=-1д<Ом2( q) Q , Q >(17)
это точно 3-й член в уравнении (08b).
С другой стороны, первые два члена (08b) относятся к частице, движущейся в потенциале, поэтому они происходят из лагранжевой частил3( q,д˙)
:
л3( q,д˙) =12мд˙2− В( q)(18)
Эта часть
л3( q,д˙)
если его вставить в (9а), он ничего не даст (ни одного члена в уравнениях движения). Теперь в (08a) половина 3-го члена и 4-й член дают
− ϕ ( q,д˙) ГВопрос˙−ф˙( q,д˙) Г Q =ггт( - ϕ G Q )(19)
поэтому мы ожидаем лагранжевой части
л4( q,д˙, Q ,Вопрос˙)
такое, что в силу (09а)
∂л4∂Вопрос˙= − ϕ ( q,д˙) Г К(20)
то есть
л4( q,д˙, Q ,Вопрос˙) =−ϕ ( q,д˙) < G Q ,Вопрос˙>(21)
Но из-за антисимметрии
г
, эта часть может быть выражена также как
л4( q,д˙, Q ,Вопрос˙) =+ϕ ( q,д˙) < ГВопрос˙, Q >(22)
поэтому вставив это во 2-й член (09a)
−∂л4∂Вопрос= − ϕ ( q,д˙) ГВопрос˙(23)
что является другой половиной третьего члена в (08a). Это значит, что
л4
, если вставить в (09a), дает 3-й и 4-й члены (08a)
ггт(∂л4∂Вопрос˙) —∂л4∂Вопрос= − 2 ϕ ( q,д˙) ГВопрос˙−ф˙( q,д˙) Г К(24)
Результат вставки
л4
в (09b) будет рассмотрен позже вместе с
л5
. 5-й член (08a) может быть получен из лагранжевой части
л5( q,д˙, Q ,Вопрос˙)
такое, что в силу (09а)
−∂л5∂Вопрос"="ф2( q,д˙)г2Вопрос(25)
так
л5( q,д˙, Q ,Вопрос˙) =-12ф2<г2Q , Q >(26)
поскольку из (A-03.c) и симметрии
г2
∂( <г2Q , Q > )∂Вопрос= (г2+(г2)Т) Q =2г2Вопрос(27)
Можно доказать (см. раздел ДОКАЗАТЕЛЬСТВО), что сумма
л45"="л4+л5
л45( q,д˙, Q ,Вопрос˙) =л4+л5= -12ф2<г2Q , Q >−ϕ ( q,д˙) < G Q ,Вопрос˙>(28)
если вставить в (09b), получится 4-й член (08b)
ггт(∂л45∂д˙) —∂л45∂д= +1д<Ом2( q) Q , G Q >(29)
В уравнении (30) ниже мы суммируем найденные части лагранжиана, и окончательный лагранжиан равен
л ( q,д˙, Q ,Вопрос˙) =12<Вопрос˙,Вопрос˙>л1−12<Ом2( q) Q , Q >л2+12мд˙2− В( q)л3− ϕ < G Q ,Вопрос˙>л4−12ф2<г2Q , Q >л5
-------------------------------------------------- ---------(30)
идентичное приведенному в статье уравнению (10).
Уравнения (31) представляют собой уравнения движения (08) со скобками под указанными пунктами, из которых лагранжевы членылм
эти предметы происходят из:
Вопрос¨л1+Ом2( q) Вл2− 2 ϕ ( q,д˙) ГВопрос˙−ф˙( q,д˙) Г Кл4+ф2( q,д˙)г2Вопросл5= 0(31а)
мд¨+∂В( q)∂дл3−1д<Ом2( q) Q , Q >л2+1д<Ом2( q) Q , G Q >л4+л5= 0(31б)
Заметим, что канонические импульсы
п ,р
сопряжено с
Q ,q
соответственно
пп"="∂л∂Вопрос˙"="Вопрос˙−д˙дGQ _"="∂л∂д˙= мд˙−1д< GQ , Р > _(32а)(32б)
где для доказательства (32b)
п"="∂л∂д˙= мд˙−1д< Q , _Вопрос˙> −д˙д2<г2Q , Q >= мд˙−1д< Q , _Вопрос˙> +д˙д2< GQ , GQ > _ _= мд˙−1д< Q , _Вопрос˙−д˙дGQ _п> = мд˙−1д< GQ , Р > _(32б′)
Уравнения (32a) и (32b) идентичны уравнениям (3.3) и (3.4) статьи соответственно, приведенным ниже.
пкп"="Вопрос˙к−д˙д∑Джгк ДжВопросДж= мд˙−1д∑j kгк ДжпкВопросДж(3.3)(3.4)
ВСПОМОГАТЕЛЬНЫЙ РАЗДЕЛ: Сжатые упрощенные выражения и правила частичного дифференцирования
Уравнения (A-01) полезны для преобразования уравнений движения из формы (01) в форму (08):
ю2кВопроск"="[Ом2( q) В ]к∑Джгк ДжВопросДж"="[ Г В ]к∑Джгк ДжВопрос˙Дж"="[ ГВопрос˙]к∑Дж ℓгj kгДж ℓВопросℓ= -∑ℓ(∑Джгк ДжгДж ℓ)Вопросℓ= -∑ℓ(г2)к лВопросℓ= -(г2В )кд¨д−д˙2д2"="ггт(д˙д) =гϕ ( q,д˙)гт"="ф˙( q,д˙)∑к , дж( − 1)к + джюкюДжВопроскВопросДж= <Ом2( q) Q , Q > − <Ом2( q) Q , G Q >(А-01.а)(А-01.б)(А-01.с)(А-01.г)(А-01.е)(А-01.f)
Доказательство (A-01.f) выполняется следующим образом.
∑к , дж( − 1)к + джюкюДжВопроскВопросДж"="∑кю2кВопрос2к<Ом2( q) Q , Q >+∑к , дж ≠ к( − 1)к + джюкюДжВопроскВопросДж− <Ом2( q) Q , G Q >(А-01.f′)
с
∑к , дж ≠ к( − 1)к + джюкюДжВопроскВопросДж"="(πд)2∑к , дж ≠ к( − 1)к + джк ДжВопроскВопросДж"="(πд)2∑к , дж ≠ к( − 1)к + дж2 кДж _Дж2−к2гк ДжДж2−к22ВопроскВопросДж"="12∑к , джгк Дж(ю2Дж−ю2к)ВопроскВопросДж"="−12∑Дж(ю2ДжВопросДж)[Ом2( q) В ]Дж∑кгj kВопроск[ Г В ]Дж−12∑к(ю2кВопроск)[Ом2( q) В ]к∑Джгк ДжВопросДж[ Г В ]к= - <Ом2( q) Q , G Q >
-------------------------------------------------- ---------(А-01.f′ ′)
Уравнения (A-02) и (A-03) полезны для преобразования лагранжиана из формы (3.1), см. соответствующее уравнение, в форму (10) и для пошагового построения этого лагранжиана из уравнения движения (08) :
∑кВопрос˙2к= <Вопрос˙,Вопрос˙> =∥∥Вопрос˙∥∥2∑кю2к( q)Вопрос2к= <Ом2Q , Q >=<Ω Q ,ОмТQ >=<Ом Q ,Ом Q >=∥ Ом Q ∥2∑дж , кгк ДжВопрос˙кВопросДж= < G Q ,Вопрос˙> = - < ГВопрос˙, Q >∑Дж , К , лгк Джгк лВопросℓВопросДж= - <г2Q , Q >=< G Q , G Q >=∥ G Q ∥2(А-02.а)(А-02.б)(А-02.с)(А-02.г)
Уравнения (A-02.c) и (A-02.d) доказываются соответственно следующим образом.
∑дж , кгк ДжВопрос˙кВопросДж"="∑к(∑Джгк ДжВопросДж)Вопрос˙к"="∑к[ Г В ]к[Вопрос˙]к"="< Q , _Вопрос˙> = < Q ,гТВопрос˙> = < Q , − GВопрос˙> = - < ГВопрос˙, Q >(А-02.с′)
∑Дж , К , лгк Джгк лВопросℓВопросДж"="∑к(∑Джгк ДжВопросДж) (∑ℓгк лВопросℓ) =∑к[ Г В ]к[ Г В ]кзнак равно < G Q , G Q > = <гТG Q , Q >=-<г2Q , Q >(А-02.д′)
Уравнения (A-03) ниже в некотором смысле являются правилами частичного дифференцирования скалярной функции векторной переменной.С
относительно этой переменной. Скалярные функции обычно являются внутренними произведениями, а переменный вектор равенС = QилиВопрос˙
. В следующихА , Р
являются векторами иФ
линейное преобразование, все они не зависят от переменного вектораС
. ОбычноF =Ом,Ом2, г ,г2
:
∂( < А , S > )∂С"="∂( < S , А > )∂С= А∂( < Р , Ф С > )∂С"="∂( <ФТР , С > )∂С"="ФТр∂( < Ф С , С > )∂С= ( Ф +ФТ) С∂( < С , С > )∂С= 2 с(А-03.а)(А-03.б)(А-03.с)(А-03.г)
(A-03.b) является частным случаем (A-03.a) с
А =ФТр
и (A-03.d) является частным случаем (A-03.c) с
Ф = Я
.
Идентичность, полезная в следующем разделе,
гТ= - Г⟹< G S , S > = 0 ,для любого действительного вектора С(А-04)
с
< G S , S > знак равно < S ,гТS >знак равно< S , ( - г ) S >знак равно-< г S , S >(А-04′)
РАЗДЕЛ ДОКАЗАТЕЛЬСТВА: Доказательство уравнения (29) с учетом уравнения (28).
Мы докажем уравнение (29) из (28), два уравнения повторены здесь для удобства.
ггт(∂л45∂д˙) —∂л45∂д= +1д<Ом2( q) Q , G Q >(29)
где
л45( q,д˙, Q ,Вопрос˙)≡деф−12ф2( q,д˙) <г2Q , Q >−ϕ ( q,д˙) < G Q ,Вопрос˙>(28)
−∂л45∂д= ф∂ф∂д<г2Q , Q >+∂ф∂д< Q , _Вопрос˙>= (д˙д)∂(д˙д)∂д<г2Q , Q >+∂(д˙д)∂д< Q , _Вопрос˙>
так
−∂л45∂д= ( -д˙2д3) <г2Q , Q >+ ( -д˙д2) < G Q ,Вопрос˙>(Б-01)
Сейчас
∂л45∂д˙∂л45∂д˙= - ϕ∂ф∂д˙<г2Q , Q >-∂ф∂д˙< Q , _Вопрос˙>⟹= ( -д˙д2) <г2Q , Q >+ ( -1д) < G Q ,Вопрос˙>(Б-02)
Дифференциация (B-02) по
т
ггт(∂л45∂д˙) = ( -д¨д− 2д˙2д3) <г2Q , Q >+ ( -д˙д2) <г2Вопрос˙, Q >+ ( -д˙д2) <г2Q ,Вопрос˙> + (д˙д2) < G Q ,Вопрос˙> + ( -1д)< ГВопрос˙,Вопрос˙>= 0 , см. (А-04)+ ( -1д) < G Q ,Вопрос¨>(Б-03)
Добавление (B-01) и (B-03)
ггт(∂л45∂д˙) —∂л45∂д"="( -д¨д−д˙2д3) <г2Q , Q >+ ( -2д˙д2) <г2Q ,Вопрос˙> + ( -1д) < G Q ,Вопрос¨> =+1д< (д¨−д˙2д2) Г К +2д˙дгВопрос˙−Вопрос¨, Q > = + _1д<ф˙G Q +2ϕ GВопрос˙−Вопрос¨"="Ом2( q) К +ф2г2Q , см. (08а), Г К >+1д<Ом2( q) К +ф2г2Q , G Q >=+1д<Ом2( q) Q , G Q > +ф2д<г2Q , G Q >= 0 , см. (А-04)
так
ггт(∂л45∂д˙) —∂л45∂д= +1д<Ом2( q) Q , G Q >(Б-04)
КЭД.
удрв