tecktom.com

この章のまとめ

2026年10月15日8分で読めますAI生成教科書『虎の巻』第1部 第3章 汎関数と変分(数学だけ)

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第3章 §1.3.7 ★★★ 汎関数・第1変分・汎関数微分・基本補題・2次の変分の5つをそろえ、次の章で Euler–Lagrange 方程式を導く

汎関数は関数を入れて数を返すもので、角括弧 F[q]F[q] で書く。経路を分点の値で表すと汎関数は多変数関数になり、偏微分を刻み幅 Δt\Delta t で割ったものの極限が汎関数微分になる。変分は、経路の族 q+ϵηq+\epsilon\eta を作って実数 ϵ\epsilon で微分する操作として定義した。向き η\eta は端で 00 になる関数に限り、こうして決めた δ\delta は時刻 tt を止めた微分なので、δq˙=ddtδq\delta\dot q=\frac{\mathrm{d} }{\mathrm{d} t}\delta q と δ∫=∫δ\delta\int=\int\delta が成り立つ。第1変分を ∫(δF/δq(t)) δq(t) dt\int(\delta F/\delta q(t))\,\delta q(t)\,\mathrm{d} t の形に書いたときの係数が汎関数微分で、∫V(q) dt\int V(q)\,\mathrm{d} t なら V′(q)V'(q)、∫12q˙2 dt\int\frac12\dot q^2\,\mathrm{d} t なら部分積分で −q¨-\ddot q、1点の値 q(t′)q(t') ならデルタ関数 δ(t−t′)\delta(t-t') になる。変分法の基本補題により、すべての向きで δF=0\delta F=0 となる停留の条件は「汎関数微分 =0=0」と同じで、汎関数の停留の条件は微分方程式になる。停留は最小を意味せず、∫(12q˙2−12ω2q2) dt\int(\frac12\dot q^2-\frac12\omega^2q^2)\,\mathrm{d} t では区間の長さが π/ω\pi/\omega を超えると停留点が最小点でなくなる。次の章(第1部第4章)では、一般の L(q,q˙,t)L(q,\dot q,t) について汎関数微分を計算して Euler–Lagrange 方程式を出し、端を固定しない場合の境界条件、懸垂線と最速降下線の問題を扱う。

この章の道具

  • 第1変分:δF=ϵ ddϵF[q+ϵη]∣ϵ=0\delta F=\epsilon\,\frac{\mathrm{d} }{\mathrm{d} \epsilon}F[q+\epsilon\eta]\big|_{\epsilon=0}、δq=ϵη\delta q=\epsilon\eta、η(t1)=η(t2)=0\eta(t_1)=\eta(t_2)=0。
  • 汎関数微分:δF=∫δFδq(t)δq(t) dt\delta F=\int\frac{\delta F}{\delta q(t)}\delta q(t)\,\mathrm{d} t。離散化では 1Δt∂SN∂qk\frac{1}{\Delta t}\frac{\partial S_N}{\partial q_k} の極限。
  • 基本補題:すべての η\eta で ∫gη dt=0\int g\eta\,\mathrm{d} t=0 なら g≡0g\equiv0。よって停留 ⇔\Leftrightarrow δF/δq=0\delta F/\delta q=0。
  • 2次の変分 δ2F\delta^2F の符号で最小かどうかを判定する。停留は最小の必要条件にすぎない。

総合問題 端を決めた2次の汎関数を4通りに調べる ★★★

オリジナル0<a<π0<a<\pi を定数とし、[0,1][0,1] 上の汎関数

F[q]=∫01(12q˙(t)2−a22q(t)2)dtF[q]=\int_0^1\Bigl(\frac12\dot q(t)^2-\frac{a^2}{2}q(t)^2\Bigr)\mathrm{d} t

を、q(0)=0q(0)=0, q(1)=1q(1)=1 を満たす経路の中で考える。

  1. (1) (使う道具:1.3.2節/近似の数理 第4章)[0,1][0,1] を N=2N=2 等分し、式 (1.3.4) の形の和 S2(q1)S_2(q_1) を書け(q0=0q_0=0, q2=1q_2=1)。その停留点 q1q_1 を aa で表し、a=1a=1 での値を求めよ。
  2. (2) (使う道具:1.3.3節)許される向き η\eta について第1変分 δF\delta F を、η\eta と η˙\dot\eta を含む積分で表せ。
  3. (3) (使う道具:1.3.4節)汎関数微分 δF/δq(t)\delta F/\delta q(t) を求めよ。
  4. (4) (使う道具:1.3.5節の基本補題)停留点 q∗q_* を求め、q∗(1/2)q_*(1/2) の a=1a=1 での値を (1) と比べよ。
  5. (5) (使う道具:1.3.5節の2次の変分と Wirtinger の不等式)F[q∗+η]−F[q∗]F[q_*+\eta]-F[q_*] を正確に求め、q∗q_* が最小点であることを示せ。a>πa>\pi では何が起こるか。
  6. (6) (使う道具:1.3.4節の部分積分)F[q∗]=a2cot⁡aF[q_*]=\dfrac a2\cot a を示し、a=1a=1 での値を求めよ。
解答

(1) Δt=12\Delta t=\frac12。和の2つの項は k=0k=0 が [12(q1/Δt)2−a22⋅02]Δt=q12\bigl[\frac12(q_1/\Delta t)^2-\frac{a^2}{2}\cdot0^2\bigr]\Delta t=q_1^2、k=1k=1 が [12((1−q1)/Δt)2−a22q12]Δt=(1−q1)2−a24q12\bigl[\frac12\bigl((1-q_1)/\Delta t\bigr)^2-\frac{a^2}{2}q_1^2\bigr]\Delta t=(1-q_1)^2-\frac{a^2}{4}q_1^2 なので

S2(q1)=q12+(1−q1)2−a24q12.S_2(q_1)=q_1^2+(1-q_1)^2-\frac{a^2}{4}q_1^2 .

dS2dq1=2q1−2(1−q1)−a22q1=0\frac{\mathrm{d} S_2}{\mathrm{d} q_1}=2q_1-2(1-q_1)-\frac{a^2}{2}q_1=0 から q1=48−a2q_1=\dfrac{4}{8-a^2}。a=1a=1 では q1=47≈0.571q_1=\frac47\approx0.571。これは折れ線の空間(内側の分点が1つ)での近似で、近似の数理 第4章の言葉では部分空間を狭めた A 型の近似である。

(2) ddϵ[12(q˙+ϵη˙)2−a22(q+ϵη)2]ϵ=0=q˙η˙−a2qη\frac{\mathrm{d} }{\mathrm{d} \epsilon}\Bigl[\frac12(\dot q+\epsilon\dot\eta)^2-\frac{a^2}{2}(q+\epsilon\eta)^2\Bigr]_{\epsilon=0}=\dot q\dot\eta-a^2q\eta なので

δF=ϵ∫01(q˙ η˙−a2q η)dt.\delta F=\epsilon\int_0^1\bigl(\dot q\,\dot\eta-a^2q\,\eta\bigr)\mathrm{d} t .

(3) 第1項を部分積分し、端の項 [q˙η]01[\dot q\eta]_0^1 が η(0)=η(1)=0\eta(0)=\eta(1)=0 で消えるので δF=ϵ∫01(−q¨−a2q)η dt\delta F=\epsilon\int_0^1(-\ddot q-a^2q)\eta\,\mathrm{d} t。よってδFδq(t)=−q¨(t)−a2q(t)\dfrac{\delta F}{\delta q(t)}=-\ddot q(t)-a^2q(t)。

(4) 基本補題から、停留は q¨+a2q=0\ddot q+a^2q=0 と同じで、一般解は q=Asin⁡at+Bcos⁡atq=A\sin at+B\cos at。q(0)=0q(0)=0 から B=0B=0、q(1)=1q(1)=1 から A=1/sin⁡aA=1/\sin a(0<a<π0<a<\pi なので sin⁡a≠0\sin a\neq0)。よって

q∗(t)=sin⁡atsin⁡a.q_*(t)=\frac{\sin at}{\sin a}.

a=1a=1 では q∗(1/2)=sin⁡12/sin⁡1=0.4794/0.8415≈0.570q_*(1/2)=\sin\frac12/\sin1=0.4794/0.8415\approx0.570 で、(1) の 0.5710.571 と3桁目まで近い。分点を増やすと離散化の停留点は q∗q_* に近づく(N=400N=400 で最大の差は 4×10−84\times10^{-8} 程度)。

(5) 代入して展開すると、FF が2次式なので

F[q∗+η]−F[q∗]=∫01(q˙∗η˙−a2q∗η)dt+12∫01(η˙2−a2η2)dt.F[q_*+\eta]-F[q_*]=\int_0^1\bigl(\dot q_*\dot\eta-a^2q_*\eta\bigr)\mathrm{d} t+\frac12\int_0^1\bigl(\dot\eta^2-a^2\eta^2\bigr)\mathrm{d} t .

第1項は (3) と同じ部分積分で ∫(−q¨∗−a2q∗)η dt=0\int(-\ddot q_*-a^2q_*)\eta\,\mathrm{d} t=0。第2項に Wirtinger の不等式 (1.3.19)(T=1T=1)を使うと

F[q∗+η]−F[q∗]=12∫01(η˙2−a2η2)dt ≥ π2−a22∫01η2 dt>0(η≠0).F[q_*+\eta]-F[q_*]=\frac12\int_0^1\bigl(\dot\eta^2-a^2\eta^2\bigr)\mathrm{d} t\ \ge\ \frac{\pi^2-a^2}{2}\int_0^1\eta^2\,\mathrm{d} t>0\qquad(\eta\neq0).

よって 0<a<π0<a<\pi で q∗q_* は唯一の最小点である。a>πa>\pi では η=sin⁡πt\eta=\sin\pi t を入れると右辺の第2項が π2−a24<0\frac{\pi^2-a^2}{4}<0 となり、停留点があっても最小点ではない(aa が π\pi の整数倍でなければ停留点 q∗q_* は存在するが、その向きに動かすと FF は減る)。これは1.3.5節の式 (1.3.18) と同じ現象で、ωT\omega T が aa にあたる。

(6) q˙2=ddt(qq˙)−qq¨\dot q^2=\frac{\mathrm{d} }{\mathrm{d} t}(q\dot q)-q\ddot q を使って第1項を書き換えると

F[q∗]=12[q∗q˙∗]01−12∫01q∗(q¨∗+a2q∗)dt=12q∗(1)q˙∗(1)=12⋅1⋅acos⁡asin⁡a=a2cot⁡a.F[q_*]=\frac12\bigl[q_*\dot q_*\bigr]_0^1-\frac12\int_0^1q_*\bigl(\ddot q_*+a^2q_*\bigr)\mathrm{d} t=\frac12q_*(1)\dot q_*(1)=\frac12\cdot1\cdot\frac{a\cos a}{\sin a}=\frac a2\cot a .

停留点では被積分関数の部分が消え、値は端の項だけで決まる。a=1a=1 では 12cot⁡1≈0.321\frac12\cot1\approx0.321。(離散化の値は N=400N=400 で 0.32170.3217 で、NN を増やすと 0.3210.321 に近づく。)

木下の勉強ステータス

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

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

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