変分法と汎関数の基礎（解析力学）

関数は「数を入力すると数を返す」ものですが、変分法では「関数を入力すると数を返す」対象を扱います。この「関数の関数」を汎関数と呼び、解析力学の数学的基盤を支える重要な概念です。

汎関数とは

普通の関数 $f (x)$ は、数 $x$ を受け取って数 $f (x)$ を返します。一方、汎関数は関数 $y (x)$ を丸ごと受け取って、1 つの数を返します。汎関数は $J [y]$ のように角括弧で表記するのが慣例です。

たとえば、2 点間を結ぶ曲線 $y (x)$ の長さは次の積分で与えられます。

$J [y] = \int_{a}^{b} 1 + (\frac{d y}{d x})^{2} d x$

この $J [y]$ は汎関数です。曲線の形 $y (x)$ を変えると $J$ の値も変わります。直線なら短く、くねくねした曲線なら長くなるわけです。

関数

数 $x$ を受け取って数 $f (x)$ を返す。 $f : R \to R$

汎関数

関数 $y (x)$ を受け取って数 $J [y]$ を返す。 $J : {関数} \to R$

身近な汎関数の例

汎関数は抽象的に見えますが、物理や幾何の問題では自然に登場します。

曲線の長さ

2 点 $A, B$ を結ぶ曲線 $y (x)$ の長さ。どんな経路をとるかによって値が変わり、直線のとき最短になります。

最速降下線

重力下で点 $A$ から点 $B$ まで滑り落ちるとき、所要時間が最短になる曲線の形を求める問題。答えはサイクロイドになることが知られています。

物理の作用

ラグランジアン $L$ を時間で積分した量 $S [q] = \int_{t_{1}}^{t_{2}} L d t$ は、解析力学における最も重要な汎関数です。

いずれの問題でも「どの関数を選べば汎関数の値が最大または最小になるか」を求めたいという共通点があります。この問いに答えるのが変分法です。

変分とは

関数の微分では $x$ を少しずらして $x + d x$ とします。変分法ではこれに対応して、関数 $y (x)$ を少しずらして $y (x) + δy (x)$ とします。この $δy (x)$ を変分と呼びます。

ここで重要なのは、 $δy (x)$ は $x$ のあらゆる点で独立に値をとれる「関数のずれ」だという点です。ただし端点は固定するのが通常の設定で、次の条件を課します。

$δy (a) = 0, δy (b) = 0$

汎関数 $J [y]$ に対して、 $y$ を $y + δy$ に置き換えたときの変化の 1 次の部分を $δ J$ と書きます。

$J [y + δy] = J [y] + δ J + (高次の項)$

普通の関数で $f^{'} (x) = 0$ となる点が極値を与えるのと同様に、汎関数が極値をとる条件は $δ J = 0$ です。変分法の基本的な流れをまとめると次のようになります。

関数 $y (x)$ を $y (x) + δy (x)$ にずらす

汎関数の変化 $δ J$ を計算する

$δ J = 0$ を要求して極値条件を得る

オイラー＝ラグランジュ方程式の導出

具体的に計算してみましょう。次の形の汎関数を考えます。

$J [y] = \int_{a}^{b} F (x, y, y^{'}) d x$

ここで $y^{'} = d y / d x$ です。 $y$ を $y + δy$ に置き換えると、被積分関数の変化は 1 次まで展開すると次のようになります。

$δ F = \frac{\partial F}{\partial y} δy + \frac{\partial F}{\partial y ^{'}} δ y^{'}$

$δ y^{'} = d (δy) / d x$ であることに注意して、第 2 項を部分積分します。

$\int_{a}^{b} \frac{\partial F}{\partial y ^{'}} δ y^{'} d x = [\frac{\partial F}{\partial y ^{'}} δy]_{a}^{b} - \int_{a}^{b} \frac{d}{d x} \frac{\partial F}{\partial y ^{'}} δy d x$

端点条件 $δy (a) = δy (b) = 0$ より境界項は消えます。したがって $δ J = 0$ の条件は次の形に集約されます。

$\int_{a}^{b} (\frac{\partial F}{\partial y} - \frac{d}{d x} \frac{\partial F}{\partial y ^{'}}) δy d x = 0$

$δy$ は任意の関数なので、括弧内がすべての $x$ でゼロでなければなりません。こうして得られるのがオイラー＝ラグランジュ方程式です。

$\frac{\partial F}{\partial y} - \frac{d}{d x} \frac{\partial F}{\partial y ^{'}} = 0$

この方程式は「汎関数を極値にする関数 $y (x)$ が満たすべき微分方程式」を与えてくれます。変分法の中心的な結果であり、解析力学のすべての議論がここから始まります。

具体例：最短経路問題

2 点 $(a, y_{a})$ と $(b, y_{b})$ を結ぶ最短曲線を求めてみましょう。曲線の長さの汎関数は次の通りです。

$J [y] = \int_{a}^{b} 1 + y^{'2} d x$

ここで $F (x, y, y^{'}) = 1 + y^{'2}$ とおきます。 $F$ は $y$ を陽に含まないため $\partial F / \partial y = 0$ となり、オイラー＝ラグランジュ方程式は次のように簡略化されます。

$\frac{d}{d x} \frac{\partial F}{\partial y ^{'}} = 0$

$\partial F / \partial y^{'}$ を計算すると $y^{'} / 1 + y^{'2}$ なので、これが定数であることから $y^{'} = const$ が得られます。つまり最短経路は直線です。直感的に当然の結果ですが、変分法を使えばこれを厳密に証明できるわけです。

解析力学への接続

変分法の枠組みを力学に応用するとき、被積分関数 $F$ に当たるものがラグランジアン $L (q, \overset{q}{˙}, t)$ になります。

汎関数 $J [y]$ に対応するのが作用と呼ばれる量で、 $S [q] = \int_{t_{1}}^{t_{2}} L (q, \overset{q}{˙}, t) d t$ と定義されます。

系の運動経路全体を評価する量。実際の運動は作用を停留させる経路に一致する（最小作用の原理）。

オイラー＝ラグランジュ方程式の変数を力学の記法に読み替えると、対応は明快です。

変分法の記法

$\frac{\partial F}{\partial y} - \frac{d}{d x} \frac{\partial F}{\partial y ^{'}} = 0$

解析力学の記法

$\frac{\partial L}{\partial q} - \frac{d}{d t} \frac{\partial L}{\partial q ˙} = 0$

変数の名前が変わっただけで、数学的構造は完全に同じです。変分法という数学的道具が、座標の選び方によらない普遍的な運動方程式を導く力を持っています。

前回の記事で導入した一般化座標 $q$ とラグランジアン $L$ をこの方程式に入れれば、任意の力学系の運動方程式が得られます。変分法と汎関数は、解析力学の論理構造を支える数学的な柱といえるでしょう。