解析力学から様々な理論分野を組み立てる 第1部 第4章 §1.4.2 ★★★ L L L が q j q_j q j を含まなければ ∂ L / ∂ q ˙ j \partial L/\partial\dot q_j ∂ L / ∂ q ˙ j が、t t t を含まなければ h = ∑ i q ˙ i ∂ L / ∂ q ˙ i − L h=\sum_i\dot q_i\partial L/\partial\dot q_i-L h = ∑ i q ˙ i ∂ L / ∂ q ˙ i − L が一定。全微分を足しても方程式は変わらない
2階の微分方程式を解くには、普通は2回積分する必要がある。解に沿って値が一定になる q , q ˙ , t q,\dot q,t q , q ˙ , t の関数が1つ見つかれば、それを「定数 = = = その関数」とおくことで、方程式を1階に下げられる。この一定になる関数を第一積分 と呼ぶ。Euler–Lagrange 方程式では、L L L の形を見るだけで第一積分が2種類見つかる。
L L L がある座標を含まないとき
L L L が q j q_j q j そのものを含まず、q ˙ j \dot q_j q ˙ j だけを含むとき、∂ L / ∂ q j = 0 \partial L/\partial q_j=0 ∂ L / ∂ q j = 0 なので、j j j 番目の Euler–Lagrange 方程式は
d d t ∂ L ∂ q ˙ j = 0 , すなわち ∂ L ∂ q ˙ j = 一定 (1.4.8) \frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \dot q_j}=0,\qquad\text{すなわち}\qquad \frac{\partial L}{\partial \dot q_j}=\text{一定}\tag{1.4.8} d t d ∂ q ˙ j ∂ L = 0 , すなわち ∂ q ˙ j ∂ L = 一定 ( 1.4.8 )
になる。L L L に現れない座標 q j q_j q j を循環座標 と呼ぶ。∂ L / ∂ q ˙ j \partial L/\partial\dot q_j ∂ L / ∂ q ˙ j は後の章で正準運動量 p j p_j p j と名付ける量で、循環座標に対応する運動量は保存する。「L L L が q j q_j q j を含まない」とは「q j q_j q j を一斉に定数だけずらしても L L L が変わらない」ということで、変換で変わらないことと保存量の関係は第1部第10章 (Noether の定理)で一般の形にする。
L L L が時刻を陽に含まないとき
命題 1.4.2(L L L が t t t を陽に含まないときの第一積分 ★★★)
h ( q , q ˙ , t ) ≔ ∑ i = 1 f q ˙ i ∂ L ∂ q ˙ i − L (1.4.9) h(q,\dot q,t)\coloneqq\sum_{i=1}^f\dot q_i\frac{\partial L}{\partial \dot q_i}-L\tag{1.4.9} h ( q , q ˙ , t ) : = i = 1 ∑ f q ˙ i ∂ q ˙ i ∂ L − L ( 1.4.9 )
とおくと、Euler–Lagrange 方程式の解に沿って
d h d t = − ∂ L ∂ t (1.4.10) \frac{\mathrm{d} h}{\mathrm{d} t}=-\frac{\partial L}{\partial t}\tag{1.4.10} d t d h = − ∂ t ∂ L ( 1.4.10 )
が成り立つ。特に L L L が t t t を陽に含まなければ、h h h は解に沿って一定である。
証明 L L L を経路に沿って t t t で微分すると、連鎖律で d L d t = ∑ i ( ∂ L ∂ q i q ˙ i + ∂ L ∂ q ˙ i q ¨ i ) + ∂ L ∂ t \frac{\mathrm{d} L}{\mathrm{d} t}=\sum_i\bigl(\frac{\partial L}{\partial q_i}\dot q_i+\frac{\partial L}{\partial \dot q_i}\ddot q_i\bigr)+\frac{\partial L}{\partial t} d t d L = ∑ i ( ∂ q i ∂ L q ˙ i + ∂ q ˙ i ∂ L q ¨ i ) + ∂ t ∂ L である。一方、積の微分でd d t ( ∑ i q ˙ i ∂ L ∂ q ˙ i ) = ∑ i ( q ¨ i ∂ L ∂ q ˙ i + q ˙ i d d t ∂ L ∂ q ˙ i ) \frac{\mathrm{d} }{\mathrm{d} t}\bigl(\sum_i\dot q_i\frac{\partial L}{\partial \dot q_i}\bigr)=\sum_i\bigl(\ddot q_i\frac{\partial L}{\partial \dot q_i}+\dot q_i\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \dot q_i}\bigr) d t d ( ∑ i q ˙ i ∂ q ˙ i ∂ L ) = ∑ i ( q ¨ i ∂ q ˙ i ∂ L + q ˙ i d t d ∂ q ˙ i ∂ L ) である。差をとると q ¨ i \ddot q_i q ¨ i の項は打ち消し合い、
d h d t = ∑ i q ˙ i ( d d t ∂ L ∂ q ˙ i − ∂ L ∂ q i ) − ∂ L ∂ t \frac{\mathrm{d} h}{\mathrm{d} t}=\sum_i\dot q_i\Bigl(\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \dot q_i}-\frac{\partial L}{\partial q_i}\Bigr)-\frac{\partial L}{\partial t} d t d h = i ∑ q ˙ i ( d t d ∂ q ˙ i ∂ L − ∂ q i ∂ L ) − ∂ t ∂ L
が残る。解の上では括弧の中が Euler–Lagrange 方程式により 0 0 0 なので、式 (1.4.10) を得る。
h h h をエネルギー関数 と呼ぶ。L = 1 2 m q ˙ 2 − V ( q ) L=\frac12m\dot q^2-V(q) L = 2 1 m q ˙ 2 − V ( q ) なら h = m q ˙ 2 − L = 1 2 m q ˙ 2 + V ( q ) h=m\dot q^2-L=\frac12m\dot q^2+V(q) h = m q ˙ 2 − L = 2 1 m q ˙ 2 + V ( q ) で、力学的エネルギーになる。h h h は q , q ˙ q,\dot q q , q ˙ の関数で、これを q q q と p = ∂ L / ∂ q ˙ p=\partial L/\partial\dot q p = ∂ L / ∂ q ˙ の関数に書き直したものが 第1部第8章 の Hamiltonian H H H である(書き直しの操作は 第1部第7章 の Legendre 変換)。1変数の場合、h = q ˙ ∂ L / ∂ q ˙ − L = C h=\dot q\,\partial L/\partial\dot q-L=C h = q ˙ ∂ L / ∂ q ˙ − L = C は q q q と q ˙ \dot q q ˙ だけの式(1階の微分方程式)で、変分法の教科書では Beltrami の恒等式とも呼ばれる。1.4.4 節と1.4.5 節の2つの問題は、どちらも独立変数を陽に含まない L L L なので、この第一積分から解く。
全微分を足しても方程式は変わらない
命題 1.4.3(全微分の不定性 ★★)
Λ ( q , t ) \Lambda(q,t) Λ ( q , t ) を任意の滑らかな関数とし、
L ′ ( q , q ˙ , t ) = L ( q , q ˙ , t ) + d d t Λ ( q , t ) = L + ∑ i ∂ Λ ∂ q i q ˙ i + ∂ Λ ∂ t L'(q,\dot q,t)=L(q,\dot q,t)+\frac{\mathrm{d} }{\mathrm{d} t}\Lambda(q,t)=L+\sum_i\frac{\partial \Lambda}{\partial q_i}\dot q_i+\frac{\partial \Lambda}{\partial t} L ′ ( q , q ˙ , t ) = L ( q , q ˙ , t ) + d t d Λ ( q , t ) = L + i ∑ ∂ q i ∂ Λ q ˙ i + ∂ t ∂ Λ
とおく。L ′ L' L ′ と L L L の Euler–Lagrange 方程式は同じである。
証明 作用の差は S ′ [ q ] − S [ q ] = ∫ t 1 t 2 d Λ d t d t = Λ ( q ( t 2 ) , t 2 ) − Λ ( q ( t 1 ) , t 1 ) S'[q]-S[q]=\int_{t_1}^{t_2}\frac{\mathrm{d} \Lambda}{\mathrm{d} t}\mathrm{d} t=\Lambda(q(t_2),t_2)-\Lambda(q(t_1),t_1) S ′ [ q ] − S [ q ] = ∫ t 1 t 2 d t d Λ d t = Λ ( q ( t 2 ) , t 2 ) − Λ ( q ( t 1 ) , t 1 ) で、端の値だけで決まる。端を固定した変分ではこの差は動かないので、δ S ′ = δ S \delta S'=\delta S δ S ′ = δ S であり、汎関数微分も停留の条件も一致する。
直接の計算でも確かめられる。足した項 Λ q q ˙ + Λ t \Lambda_q\dot q+\Lambda_t Λ q q ˙ + Λ t (1変数で書き、Λ q = ∂ Λ / ∂ q \Lambda_q=\partial\Lambda/\partial q Λ q = ∂ Λ/ ∂ q などとした)の Euler–Lagrange の左辺はd d t Λ q − ( Λ q q q ˙ + Λ q t ) = ( Λ q q q ˙ + Λ q t ) − ( Λ q q q ˙ + Λ q t ) = 0 \frac{\mathrm{d} }{\mathrm{d} t}\Lambda_q-(\Lambda_{qq}\dot q+\Lambda_{qt})=(\Lambda_{qq}\dot q+\Lambda_{qt})-(\Lambda_{qq}\dot q+\Lambda_{qt})=0 d t d Λ q − ( Λ q q q ˙ + Λ q t ) = ( Λ q q q ˙ + Λ q t ) − ( Λ q q q ˙ + Λ q t ) = 0 である。1つの方程式を導く L L L は1つに決まらず、全微分だけの不定性がある。電磁場中の荷電粒子でゲージ変換が L L L に全微分を足す操作になること(第2部第2章 )と、Noether の定理で「L L L が全微分だけ変わる変換」にも保存量が対応すること(第1部第10章 )は、この性質を使う。
例題 1.4.2 全微分の項と、第一積分で階数を下げる ★★
オリジナル
(1) L A = 1 2 q ˙ 2 + q ˙ sin q L_A=\frac12\dot q^2+\dot q\sin q L A = 2 1 q ˙ 2 + q ˙ sin q の Euler–Lagrange 方程式を直接計算し、L = 1 2 q ˙ 2 L=\frac12\dot q^2 L = 2 1 q ˙ 2 のものと同じであることを確かめよ。足した項は何の全微分か。
(2) L B = 1 2 ( 1 + q 2 ) q ˙ 2 L_B=\frac12(1+q^2)\dot q^2 L B = 2 1 ( 1 + q 2 ) q ˙ 2 の Euler–Lagrange 方程式を書け。
(3) L B L_B L B のエネルギー関数 h h h を求め、初期条件 q ( 0 ) = 0 q(0)=0 q ( 0 ) = 0 , q ˙ ( 0 ) = 1 \dot q(0)=1 q ˙ ( 0 ) = 1 の解について、q q q が 1 1 1 に達する時刻を求めよ。
解答 (1) ∂ L A ∂ q = q ˙ cos q \frac{\partial L_A}{\partial q}=\dot q\cos q ∂ q ∂ L A = q ˙ cos q , ∂ L A ∂ q ˙ = q ˙ + sin q \frac{\partial L_A}{\partial \dot q}=\dot q+\sin q ∂ q ˙ ∂ L A = q ˙ + sin q なので
d d t ( q ˙ + sin q ) − q ˙ cos q = q ¨ + q ˙ cos q − q ˙ cos q = q ¨ . \frac{\mathrm{d} }{\mathrm{d} t}(\dot q+\sin q)-\dot q\cos q=\ddot q+\dot q\cos q-\dot q\cos q=\ddot q . d t d ( q ˙ + sin q ) − q ˙ cos q = q ¨ + q ˙ cos q − q ˙ cos q = q ¨ .
Euler–Lagrange 方程式は q ¨ = 0 \ddot q=0 q ¨ = 0 で、L = 1 2 q ˙ 2 L=\frac12\dot q^2 L = 2 1 q ˙ 2 と同じである。足した項は q ˙ sin q = d d t ( − cos q ) \dot q\sin q=\frac{\mathrm{d} }{\mathrm{d} t}(-\cos q) q ˙ sin q = d t d ( − cos q ) で、Λ = − cos q \Lambda=-\cos q Λ = − cos q の全微分である。
(2) ∂ L B ∂ q = q q ˙ 2 \frac{\partial L_B}{\partial q}=q\dot q^2 ∂ q ∂ L B = q q ˙ 2 , ∂ L B ∂ q ˙ = ( 1 + q 2 ) q ˙ \frac{\partial L_B}{\partial \dot q}=(1+q^2)\dot q ∂ q ˙ ∂ L B = ( 1 + q 2 ) q ˙ で、d d t [ ( 1 + q 2 ) q ˙ ] = ( 1 + q 2 ) q ¨ + 2 q q ˙ 2 \frac{\mathrm{d} }{\mathrm{d} t}\bigl[(1+q^2)\dot q\bigr]=(1+q^2)\ddot q+2q\dot q^2 d t d [ ( 1 + q 2 ) q ˙ ] = ( 1 + q 2 ) q ¨ + 2 q q ˙ 2 なので
( 1 + q 2 ) q ¨ + q q ˙ 2 = 0. (1+q^2)\ddot q+q\dot q^2=0 . ( 1 + q 2 ) q ¨ + q q ˙ 2 = 0.
q ¨ \ddot q q ¨ について解けるが、このまま2回積分するのは面倒である。
(3) L B L_B L B は t t t を陽に含まないので h = q ˙ ∂ L B ∂ q ˙ − L B = ( 1 + q 2 ) q ˙ 2 − 1 2 ( 1 + q 2 ) q ˙ 2 = 1 2 ( 1 + q 2 ) q ˙ 2 h=\dot q\,\frac{\partial L_B}{\partial \dot q}-L_B=(1+q^2)\dot q^2-\frac12(1+q^2)\dot q^2=\frac12(1+q^2)\dot q^2 h = q ˙ ∂ q ˙ ∂ L B − L B = ( 1 + q 2 ) q ˙ 2 − 2 1 ( 1 + q 2 ) q ˙ 2 = 2 1 ( 1 + q 2 ) q ˙ 2 が一定である。初期条件から h = 1 2 h=\frac12 h = 2 1 なので ( 1 + q 2 ) q ˙ 2 = 1 (1+q^2)\dot q^2=1 ( 1 + q 2 ) q ˙ 2 = 1 。q ˙ ( 0 ) > 0 \dot q(0)>0 q ˙ ( 0 ) > 0 で、q ˙ = 1 / 1 + q 2 \dot q=1/\sqrt{1+q^2} q ˙ = 1/ 1 + q 2 は 0 0 0 にならないので q ˙ > 0 \dot q>0 q ˙ > 0 が続き、
1 + q 2 d q = d t , t = ∫ 0 q 1 + u 2 d u = 1 2 ( q 1 + q 2 + arsinh q ) . \sqrt{1+q^2}\,\mathrm{d} q=\mathrm{d} t,\qquad t=\int_0^q\sqrt{1+u^2}\,\mathrm{d} u=\frac12\Bigl(q\sqrt{1+q^2}+\operatorname{arsinh}q\Bigr). 1 + q 2 d q = d t , t = ∫ 0 q 1 + u 2 d u = 2 1 ( q 1 + q 2 + arsinh q ) .
q = 1 q=1 q = 1 を入れると t = 1 2 ( 2 + arsinh 1 ) = 1 2 ( 2 + ln ( 1 + 2 ) ) ≈ 1.148 t=\frac12\bigl(\sqrt2+\operatorname{arsinh}1\bigr)=\frac12\bigl(\sqrt2+\ln(1+\sqrt2)\bigr)\approx1.148 t = 2 1 ( 2 + arsinh 1 ) = 2 1 ( 2 + ln ( 1 + 2 ) ) ≈ 1.148 である。2階の方程式を2回積分する代わりに、第一積分で1階に下げたので、変数分離の1回の積分で済んだ。(数値の確かめ:(2) の方程式を数値積分すると、q = 1 q=1 q = 1 に達する時刻は 1.1478 1.1478 1.1478 で一致する。)
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com