tecktom.com

第一積分:方程式の階数を1つ下げる

2026年10月15日9分で読めますAI生成教科書『虎の巻』第1部 第4章 Euler–Lagrange 方程式と境界条件(数学だけ)
目次
  1. がある座標を含まないとき
  2. が時刻を陽に含まないとき
  3. 全微分を足しても方程式は変わらない
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第4章 §1.4.2 ★★★ LL が qjq_j を含まなければ ∂L/∂q˙j\partial L/\partial\dot q_j が、tt を含まなければ h=∑iq˙i∂L/∂q˙i−Lh=\sum_i\dot q_i\partial L/\partial\dot q_i-L が一定。全微分を足しても方程式は変わらない

2階の微分方程式を解くには、普通は2回積分する必要がある。解に沿って値が一定になる q,q˙,tq,\dot q,t の関数が1つ見つかれば、それを「定数 == その関数」とおくことで、方程式を1階に下げられる。この一定になる関数を第一積分と呼ぶ。Euler–Lagrange 方程式では、LL の形を見るだけで第一積分が2種類見つかる。

LL がある座標を含まないとき

LL が qjq_j そのものを含まず、q˙j\dot q_j だけを含むとき、∂L/∂qj=0\partial L/\partial q_j=0 なので、jj 番目の Euler–Lagrange 方程式は

ddt∂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}

になる。LL に現れない座標 qjq_j を循環座標と呼ぶ。∂L/∂q˙j\partial L/\partial\dot q_j は後の章で正準運動量 pjp_j と名付ける量で、循環座標に対応する運動量は保存する。「LL が qjq_j を含まない」とは「qjq_j を一斉に定数だけずらしても LL が変わらない」ということで、変換で変わらないことと保存量の関係は第1部第10章(Noether の定理)で一般の形にする。

LL が時刻を陽に含まないとき

命題 1.4.2(LL が tt を陽に含まないときの第一積分 ★★★)

h(q,q˙,t)≔∑i=1fq˙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}

とおくと、Euler–Lagrange 方程式の解に沿って

dhdt=−∂L∂t(1.4.10)\frac{\mathrm{d} h}{\mathrm{d} t}=-\frac{\partial L}{\partial t}\tag{1.4.10}

が成り立つ。特に LL が tt を陽に含まなければ、hh は解に沿って一定である。

証明

LL を経路に沿って tt で微分すると、連鎖律で dLdt=∑i(∂L∂qiq˙i+∂L∂q˙iq¨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} である。一方、積の微分でddt(∑iq˙i∂L∂q˙i)=∑i(q¨i∂L∂q˙i+q˙iddt∂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) である。差をとると q¨i\ddot q_i の項は打ち消し合い、

dhdt=∑iq˙i(ddt∂L∂q˙i−∂L∂qi)−∂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}

が残る。解の上では括弧の中が Euler–Lagrange 方程式により 00 なので、式 (1.4.10) を得る。

hh をエネルギー関数と呼ぶ。L=12mq˙2−V(q)L=\frac12m\dot q^2-V(q) なら h=mq˙2−L=12mq˙2+V(q)h=m\dot q^2-L=\frac12m\dot q^2+V(q) で、力学的エネルギーになる。hh は q,q˙q,\dot q の関数で、これを qq と p=∂L/∂q˙p=\partial L/\partial\dot q の関数に書き直したものが 第1部第8章 の Hamiltonian HH である(書き直しの操作は 第1部第7章 の Legendre 変換)。1変数の場合、h=q˙ ∂L/∂q˙−L=Ch=\dot q\,\partial L/\partial\dot q-L=C は qq と q˙\dot q だけの式(1階の微分方程式)で、変分法の教科書では Beltrami の恒等式とも呼ばれる。1.4.4節と1.4.5節の2つの問題は、どちらも独立変数を陽に含まない LL なので、この第一積分から解く。

全微分を足しても方程式は変わらない

命題 1.4.3(全微分の不定性 ★★)

Λ(q,t)\Lambda(q,t) を任意の滑らかな関数とし、

L′(q,q˙,t)=L(q,q˙,t)+ddtΛ(q,t)=L+∑i∂Λ∂qiq˙i+∂Λ∂tL'(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′L' と LL の Euler–Lagrange 方程式は同じである。

証明

作用の差は S′[q]−S[q]=∫t1t2dΛdtdt=Λ(q(t2),t2)−Λ(q(t1),t1)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′=δS\delta S'=\delta S であり、汎関数微分も停留の条件も一致する。

直接の計算でも確かめられる。足した項 Λqq˙+Λt\Lambda_q\dot q+\Lambda_t(1変数で書き、Λq=∂Λ/∂q\Lambda_q=\partial\Lambda/\partial q などとした)の Euler–Lagrange の左辺はddtΛq−(Λqqq˙+Λqt)=(Λqqq˙+Λqt)−(Λqqq˙+Λqt)=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 である。1つの方程式を導く LL は1つに決まらず、全微分だけの不定性がある。電磁場中の荷電粒子でゲージ変換が LL に全微分を足す操作になること(第2部第2章)と、Noether の定理で「LL が全微分だけ変わる変換」にも保存量が対応すること(第1部第10章)は、この性質を使う。

例題 1.4.2 全微分の項と、第一積分で階数を下げる ★★

オリジナル

  1. (1) LA=12q˙2+q˙sin⁡qL_A=\frac12\dot q^2+\dot q\sin q の Euler–Lagrange 方程式を直接計算し、L=12q˙2L=\frac12\dot q^2 のものと同じであることを確かめよ。足した項は何の全微分か。
  2. (2) LB=12(1+q2)q˙2L_B=\frac12(1+q^2)\dot q^2 の Euler–Lagrange 方程式を書け。
  3. (3) LBL_B のエネルギー関数 hh を求め、初期条件 q(0)=0q(0)=0, q˙(0)=1\dot q(0)=1 の解について、qq が 11 に達する時刻を求めよ。
解答

(1) ∂LA∂q=q˙cos⁡q\frac{\partial L_A}{\partial q}=\dot q\cos q, ∂LA∂q˙=q˙+sin⁡q\frac{\partial L_A}{\partial \dot q}=\dot q+\sin q なので

ddt(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 .

Euler–Lagrange 方程式は q¨=0\ddot q=0 で、L=12q˙2L=\frac12\dot q^2 と同じである。足した項は q˙sin⁡q=ddt(−cos⁡q)\dot q\sin q=\frac{\mathrm{d} }{\mathrm{d} t}(-\cos q) で、Λ=−cos⁡q\Lambda=-\cos q の全微分である。

(2) ∂LB∂q=qq˙2\frac{\partial L_B}{\partial q}=q\dot q^2, ∂LB∂q˙=(1+q2)q˙\frac{\partial L_B}{\partial \dot q}=(1+q^2)\dot q で、ddt[(1+q2)q˙]=(1+q2)q¨+2qq˙2\frac{\mathrm{d} }{\mathrm{d} t}\bigl[(1+q^2)\dot q\bigr]=(1+q^2)\ddot q+2q\dot q^2 なので

(1+q2)q¨+qq˙2=0.(1+q^2)\ddot q+q\dot q^2=0 .

q¨\ddot q について解けるが、このまま2回積分するのは面倒である。

(3) LBL_B は tt を陽に含まないので h=q˙ ∂LB∂q˙−LB=(1+q2)q˙2−12(1+q2)q˙2=12(1+q2)q˙2h=\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=12h=\frac12 なので (1+q2)q˙2=1(1+q^2)\dot q^2=1。q˙(0)>0\dot q(0)>0 で、q˙=1/1+q2\dot q=1/\sqrt{1+q^2} は 00 にならないので q˙>0\dot q>0 が続き、

1+q2 dq=dt,t=∫0q1+u2 du=12(q1+q2+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).

q=1q=1 を入れると t=12(2+arsinh⁡1)=12(2+ln⁡(1+2))≈1.148t=\frac12\bigl(\sqrt2+\operatorname{arsinh}1\bigr)=\frac12\bigl(\sqrt2+\ln(1+\sqrt2)\bigr)\approx1.148 である。2階の方程式を2回積分する代わりに、第一積分で1階に下げたので、変数分離の1回の積分で済んだ。(数値の確かめ:(2) の方程式を数値積分すると、q=1q=1 に達する時刻は 1.14781.1478 で一致する。)

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ