tecktom.com

汎関数微分

2026年10月15日10分で読めますAI生成教科書『虎の巻』第1部 第3章 汎関数と変分(数学だけ)
目次
  1. 定義
  2. 基本的な例
  3. (a) 掛けて積分する汎関数
  4. (b) 被積分関数が だけの関数
  5. (c) 1点の値
  6. (d) 導関数を含む汎関数
  7. 積と合成の規則
  8. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第3章 §1.3.4 ★★★ 第1変分を ∫(⋯ )δq dt\int(\cdots)\delta q\,\mathrm{d} t の形に書いたときの (⋯ )(\cdots) が汎関数微分。偏微分の添字 ii を連続の tt にしたもの

定義

多変数の方向微分 (1.3.1) は、向き hh との内積 ∑i(∂f/∂xi)hi\sum_i(\partial f/\partial x_i)h_i の形をしていた。添字 ii についての和を tt についての積分に置き換えると、次の定義になる。

定義 1.3.3(汎関数微分)

汎関数 FF の第1変分が、すべての許される向き η\eta について

δF=∫t1t2δFδq(t) δq(t) dt(δq=ϵη)(1.3.11)\delta F=\int_{t_1}^{t_2}\frac{\delta F}{\delta q(t)}\,\delta q(t)\,\mathrm{d} t\qquad(\delta q=\epsilon\eta)\tag{1.3.11}

の形に書けるとき、η\eta によらない tt の関数 δF/δq(t)\delta F/\delta q(t) を FF の qq における汎関数微分と呼ぶ。

1.3.2節の離散化で言えば、δF=∑k(∂SN/∂qk)δqk=∑k(1Δt∂SN/∂qk)δqk Δt\delta F=\sum_k(\partial S_N/\partial q_k)\delta q_k=\sum_k\bigl(\frac{1}{\Delta t}\partial S_N/\partial q_k\bigr)\delta q_k\,\Delta t の括弧の中が δF/δq(tk)\delta F/\delta q(t_k) にあたり、和に Δt\Delta t を掛けたものが積分になる。式 (1.3.6) で偏微分を Δt\Delta t で割ったのはこのためである。δF/δq(t)\delta F/\delta q(t) は時刻 tt ごとに値を持つ関数で、「時刻 tt の近くで経路を単位の面積だけ押し上げたときの FF の変化率」と読める。

η\eta について線型な δF\delta F を ∫g η dt\int g\,\eta\,\mathrm{d} t の形に書く方法が1通りしかないこと(gg が一意に決まること)は、1.3.5節の基本補題から従う。したがって汎関数微分は、存在すれば1つに決まる。

基本的な例

(a) 掛けて積分する汎関数

F[q]=∫g(t)q(t) dtF[q]=\int g(t)q(t)\,\mathrm{d} t(gg は与えられた関数)なら δF=∫g δq dt\delta F=\int g\,\delta q\,\mathrm{d} t なので、δF/δq(t)=g(t)\delta F/\delta q(t)=g(t)。多変数の f=∑igixif=\sum_ig_ix_i で ∂f/∂xi=gi\partial f/\partial x_i=g_i となるのと同じである。

(b) 被積分関数が qq だけの関数

F[q]=∫V(q(t)) dtF[q]=\int V(q(t))\,\mathrm{d} t なら、積分の中で普通に微分して δF=∫V′(q) δq dt\delta F=\int V'(q)\,\delta q\,\mathrm{d} t。よって

δδq(t)∫V(q(t′)) dt′=V′(q(t)).(1.3.12)\frac{\delta}{\delta q(t)}\int V(q(t'))\,\mathrm{d} t'=V'\bigl(q(t)\bigr).\tag{1.3.12}

F1[q]=∫q2 dtF_1[q]=\int q^2\,\mathrm{d} t なら δF1/δq(t)=2q(t)\delta F_1/\delta q(t)=2q(t) である。

(c) 1点の値

F2[q]=q(t0)F_2[q]=q(t_0) なら δF2=ϵη(t0)\delta F_2=\epsilon\eta(t_0)。これを ∫(⋯ )ϵη(t) dt\int(\cdots)\epsilon\eta(t)\,\mathrm{d} t の形に書ける普通の関数は無い。Dirac のデルタ関数 δ(t−t0)\delta(t-t_0)(∫δ(t−t0)η(t) dt=η(t0)\int\delta(t-t_0)\eta(t)\,\mathrm{d} t=\eta(t_0) で定義される超関数)を許すとδF2/δq(t)=δ(t−t0)\delta F_2/\delta q(t)=\delta(t-t_0) と書ける。t0t_0 を t′t' と書き直すと

δq(t′)δq(t)=δ(t−t′)(1.3.13)\frac{\delta q(t')}{\delta q(t)}=\delta(t-t')\tag{1.3.13}

となる。多変数の ∂xj/∂xi=δij\partial x_j/\partial x_i=\delta_{ij}(Kronecker のデルタ)の連続版で、場の理論の章(第2部第5章)でこの形をよく使う。分子の δ\delta は変分、右辺の δ(t−t′)\delta(t-t') はデルタ関数で、同じ字だが引数の有無で区別する。

(d) 導関数を含む汎関数

F[q]=∫t1t212q˙2 dtF[q]=\int_{t_1}^{t_2}\frac12\dot q^2\,\mathrm{d} t なら、規則 (1.3.10) を使って δF=∫q˙ δq˙ dt=∫q˙ ddtδq dt\delta F=\int\dot q\,\delta\dot q\,\mathrm{d} t=\int\dot q\,\frac{\mathrm{d} }{\mathrm{d} t}\delta q\,\mathrm{d} t。このままでは δq\delta q でなく ddtδq\frac{\mathrm{d} }{\mathrm{d} t}\delta q が掛かっているので、定義 (1.3.11) の形になっていない。部分積分で導関数を q˙\dot q の側へ移すと

δF=[q˙ δq]t1t2−∫t1t2q¨ δq dt=−∫t1t2q¨ δq dt,δδq(t)∫12q˙2 dt′=−q¨(t).(1.3.14)\delta F=\Bigl[\dot q\,\delta q\Bigr]_{t_1}^{t_2}-\int_{t_1}^{t_2}\ddot q\,\delta q\,\mathrm{d} t=-\int_{t_1}^{t_2}\ddot q\,\delta q\,\mathrm{d} t, \qquad \frac{\delta}{\delta q(t)}\int\frac12\dot q^2\,\mathrm{d} t'=-\ddot q(t).\tag{1.3.14}

許される向きが端で 00 になる条件 (1.3.7) により、端の項は消える。端を固定しない問題では端の項が残り、それが境界条件を決める(第1部第4章)。この結果は1.3.2節の極限 (1.3.6) と、例題 1.3.2 の −q¨-\ddot q に一致する。例題 1.3.3 の F=∫q˙2 dtF=\int\dot q^2\,\mathrm{d} t(係数が2倍)では δF/δq=−2q¨\delta F/\delta q=-2\ddot q で、q=tq=t なら 00、q=t2q=t^2 なら −4-4 である。後者で ∫01(−4)sin⁡πt dt=−8/π\int_0^1(-4)\sin\pi t\,\mathrm{d} t=-8/\pi となり、例題の1次の係数と合う。q=tq=t では汎関数微分が恒等的に 00 なので、どの向きの第1変分も 00 になる。

被積分関数が qq と q˙\dot q の両方を含む一般の L(q,q˙,t)L(q,\dot q,t) でも、(b) と (d) を組み合わせれば汎関数微分は求まる。その結果の式と、それを 00 とおいた方程式(Euler–Lagrange 方程式)は次の章で導く。

積と合成の規則

汎関数微分は、普通の微分と同じ線型性・積の規則・合成の規則(連鎖律)を満たす。どれも、ϵ\epsilon についての普通の微分の規則を第1変分 (1.3.9) に当てはめると出る。たとえば2つの汎関数の積 F[q]G[q]F[q]G[q] なら

δ(FG)δq(t)=δFδq(t)G[q]+F[q]δGδq(t),\frac{\delta(FG)}{\delta q(t)}=\frac{\delta F}{\delta q(t)}G[q]+F[q]\frac{\delta G}{\delta q(t)},

数の関数 Φ\Phi との合成 Φ(F[q])\Phi(F[q]) なら δΦ(F)/δq(t)=Φ′(F[q]) δF/δq(t)\delta\Phi(F)/\delta q(t)=\Phi'(F[q])\,\delta F/\delta q(t) となる。右辺の G[q]G[q]・F[q]F[q]・Φ′(F[q])\Phi'(F[q]) は tt によらない数で、tt への依存はすべて汎関数微分の側にある。

例題 1.3.4 曲線の長さの汎関数微分は曲率 ★★★

オリジナル平面上の曲線 y(x)y(x)(x1≤x≤x2x_1\le x\le x_2、両端の値を固定)の長さ ℓ[y]=∫x1x21+y′2 dx\ell[y]=\int_{x_1}^{x_2}\sqrt{1+y'^2}\,\mathrm{d} x を考える。

  1. (1) 汎関数微分 δℓ/δy(x)\delta\ell/\delta y(x) を求め、−y′′(1+y′2)3/2-\dfrac{y''}{(1+y'^2)^{3/2}} になることを示せ。
  2. (2) 半径 RR の円の上半分 y=R2−x2y=\sqrt{R^2-x^2}(∣x∣<R\lvert x\rvert <R)で、δℓ/δy(x)\delta\ell/\delta y(x) を求めよ。結果を図形の言葉で説明せよ。
  3. (3) 2点を結ぶ曲線のうち ℓ\ell を停留させるものは線分に限ることを示せ(停留と汎関数微分の関係は次の節で示すので、ここでは「汎関数微分が恒等的に 00」を停留の条件として使ってよい)。
解答

(1) δ1+y′2=y′1+y′2 δy′\delta\sqrt{1+y'^2}=\dfrac{y'}{\sqrt{1+y'^2}}\,\delta y' で、δy′=ddxδy\delta y'=\frac{\mathrm{d} }{\mathrm{d} x}\delta y である。部分積分すると端の項が消えて

δℓ=∫x1x2y′1+y′2ddxδy dx=−∫x1x2ddx(y′1+y′2)δy dx.\delta\ell=\int_{x_1}^{x_2}\frac{y'}{\sqrt{1+y'^2}}\frac{\mathrm{d} }{\mathrm{d} x}\delta y\,\mathrm{d} x=-\int_{x_1}^{x_2}\frac{\mathrm{d} }{\mathrm{d} x}\Bigl(\frac{y'}{\sqrt{1+y'^2}}\Bigr)\delta y\,\mathrm{d} x .

商の微分で ddxy′1+y′2=y′′(1+y′2)−y′⋅y′y′′(1+y′2)3/2=y′′(1+y′2)3/2\frac{\mathrm{d} }{\mathrm{d} x}\dfrac{y'}{\sqrt{1+y'^2}}=\dfrac{y''(1+y'^2)-y'\cdot y'y''}{(1+y'^2)^{3/2}}=\dfrac{y''}{(1+y'^2)^{3/2}} なので

δℓδy(x)=−y′′(1+y′2)3/2.\frac{\delta\ell}{\delta y(x)}=-\frac{y''}{(1+y'^2)^{3/2}} .

右辺の y′′/(1+y′2)3/2y''/(1+y'^2)^{3/2} は曲線 y(x)y(x) の(符号つきの)曲率 κ\kappa である。長さの汎関数微分は −κ-\kappa になる。

(2) y′=−x/R2−x2y'=-x/\sqrt{R^2-x^2}, 1+y′2=R2/(R2−x2)1+y'^2=R^2/(R^2-x^2), y′′=−R2/(R2−x2)3/2y''=-R^2/(R^2-x^2)^{3/2} なので

δℓδy(x)=R2(R2−x2)3/2⋅(R2−x2)3/2R3=1R.\frac{\delta\ell}{\delta y(x)}=\frac{R^2}{(R^2-x^2)^{3/2}}\cdot\frac{(R^2-x^2)^{3/2}}{R^3}=\frac1R .

xx によらず正の定数になる。上に凸な円弧をどこで上へ(外側へ)押しても長さは増え、増え方は単位の面積あたり 1/R1/R で、曲がりの強い(RR の小さい)円ほど大きい。

(3) 停留の条件は y′′/(1+y′2)3/2=0y''/(1+y'^2)^{3/2}=0、すなわち y′′=0y''=0 である。解は y=ax+by=ax+b で、両端の値から a,ba,b が1つに決まる。よって停留する曲線は2点を結ぶ線分に限る。線分が実際に最短であることは、次の節の2次の変分で確かめられる(線分のまわりの2次の変分は 12∫δy′2/(1+y′2)3/2 dx≥0\frac12\int\delta y'^2/(1+y'^2)^{3/2}\,\mathrm{d} x\ge0 である)。

(数値の確かめ:y=x2y=x^2(0≤x≤10\le x\le1)と η=sin⁡πx\eta=\sin\pi x で ddϵℓ[y+ϵη]∣0\frac{\mathrm{d} }{\mathrm{d} \epsilon}\ell[y+\epsilon\eta]\big|_0 を数値微分すると −0.52767-0.52767 で、∫01(−2/(1+4x2)3/2)sin⁡πx dx\int_0^1\bigl(-2/(1+4x^2)^{3/2}\bigr)\sin\pi x\,\mathrm{d} x と小数点以下8桁まで一致する。)

木下の勉強ステータス

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

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

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