オイラー・ラグランジュ方程式汎関数の極値を求めることは、関数の最大値と最小値を求めることと似ています。関数の最大値と最小値は、その導関数がゼロになる点(つまり、導関数がゼロになる点)を見つけることで求められます。汎関数の極値は、汎関数の導関数 がゼロになる関数を見つけることで得られます。これは、関連するオイラー・ラグランジュ方程式 を解くことにつながります。[ g ]
関数を考える
J [ y ] = ∫ x 1 x 2 L ( x 、 y ( x ) 、 y ′ ( x ) ) d x 、 {\displaystyle J[y]=\int _{x_{1}}^{x_{2}}L\left(x,y(x),y'(x)\right)\,dx,}
どこ
x 1 、 x 2 {\displaystyle x_{1},x_{2}} 定数 です。y ( x ) {\displaystyle y(x)} 2回連続微分可能であり、y ′ ( x ) = d y d x 、 {\displaystyle y'(x)={\frac {dy}{dx}},} L ( x 、 y ( x ) 、 y ′ ( x ) ) {\displaystyle L\left(x,y(x),y'(x)\right)} 引数に関して2回連続微分可能であるx 、 y 、 {\displaystyle x,y,} そしてy ′ 。 {\displaystyle y'.} 機能的J [ y ] {\displaystyle J[y]} で局所最小値 に達するf 、 {\displaystyle f,} そしてη ( x ) {\displaystyle \eta (x)} は、少なくとも1つの導関数を持ち、両端でゼロになる任意の関数である。x 1 {\displaystyle x_{1}} そしてx 2 、 {\displaystyle x_{2},} すると任意の数に対してε {\displaystyle \varepsilon } 0に近い、
J [ f ] ≤ J [ f + ε η ] 。 {\displaystyle J[f]\leq J[f+\varepsilon \eta ]\,.}
用語ε η {\displaystyle \varepsilon \eta } これは関数の変形 と呼ばれますf {\displaystyle f} そして、で表されます。δ f 。 {\displaystyle \delta f.} [ 1 ] [ h ]
代入f + ε η {\displaystyle f+\varepsilon \eta } のためにy {\displaystyle y} 機能的にJ [ y ] 、 {\displaystyle J[y],} 結果は関数であるε 、 {\displaystyle \varepsilon ,}
Φ ( ε ) = J [ f + ε η ] 。 {\displaystyle \Phi (\varepsilon )=J[f+\varepsilon \eta ]\,.}
機能的であるためJ [ y ] {\displaystyle J[y]} 最低限のy = f {\displaystyle y=f} 機能Φ ( ε ) {\displaystyle \Phi (\varepsilon )} 最低値はε = 0 {\displaystyle \varepsilon =0} したがって、[ i ]
Φ ′ ( 0 ) ≡ d Φ d ε | ε = 0 = ∫ x 1 x 2 d L d ε | ε = 0 d x = 0 。 {\displaystyle \Phi '(0)\equiv \left.{\frac {d\Phi }{d\varepsilon }}\right|_{\varepsilon =0}=\int _{x_{1}}^{x_{2}}\left.{\frac {dL}{d\varepsilon }}\right|_{\varepsilon =0}dx=0\,.}
の全微分 をとるとL [ x 、 y 、 y ′ ] 、 {\displaystyle L\left[x,y,y'\right],} どこy = f + ε η {\displaystyle y=f+\varepsilon \eta } そしてy ′ = f ′ + ε η ′ {\displaystyle y'=f'+\varepsilon \eta '} の関数として考えられますε {\displaystyle \varepsilon } それよりもx 、 {\displaystyle x,} 収量
d L d ε = ∂ L ∂ y d y d ε + ∂ L ∂ y ′ d y ′ d ε ${\displaystyle {\frac {dL}{d\varepsilon }}={\frac {\partial L}{\partial y}}{\frac {dy}{d\varepsilon }}+{\frac {\partial L}{\partial y'}}{\frac {dy'}{d\varepsilon }}}$
そしてなぜならd y d ε = η {\displaystyle {\frac {dy}{d\varepsilon }}=\eta } そしてd y ′ d ε = η ′ 、 {\displaystyle {\frac {dy'}{d\varepsilon }}=\eta ',}
d L d ε = ∂ L ∂ y η + ∂ L ∂ y ′ η ′ 。 \displaystyle {\frac {dL}{d\varepsilon }}={\frac {\partial L}{\partial y}}\eta +{\frac {\partial L}{\partial y'}}\eta '.}
したがって、
∫ x 1 x 2 d L d ε | ε = 0 d x = ∫ x 1 x 2 ( ∂ L ∂ f η + ∂ L ∂ f ′ η ′ ) d x = ∫ x 1 x 2 ∂ L ∂ f η d x + ∂ L ∂ f ′ η | x 1 x 2 − ∫ x 1 x 2 η d d x ∂ L ∂ f ′ d x = ∫ x 1 x 2 ( ∂ L ∂ f η − η d d x ∂ L ∂ f ′ ) d x {\displaystyle {\begin{aligned}\int _{x_{1}}^{x_{2}}\left.{\frac {dL}{d\varepsilon }}\right|_{\varepsilon =0}dx&=\int _{x_{1}}^{x_{2}}\left({\frac {\partial L}{\partial f}}\eta +{\frac {\partial L}{\partial f'}}\eta '\right)\,dx\\&=\int _{x_{1}}^{x_{2}}{\frac {\partial L}{\partial f}}\eta \,dx+\left.{\frac {\partial L}{\partial f'}}\eta \right|_{x_{1}}^{x_{2}}-\int _{x_{1}}^{x_{2}}\eta {\frac {d}{dx}}{\frac {\partial L}{\partial f'}}\,dx\\&=\int _{x_{1}}^{x_{2}}\left({\frac {\partial L}{\partial f}}\eta -\eta {\frac {d}{dx}}{\frac {\partial L}{\partial f'}}\right)\,dx\\\end{aligned}}}
どこL [ x 、 y 、 y ′ ] → L [ x 、 f 、 f ′ ] {\displaystyle L\left[x,y,y'\right]\to L\left[x,f,f'\right]} いつε = 0 {\displaystyle \varepsilon =0} そして、第2項には部分積分 を用いました。第2行の第2項は、η = 0 {\displaystyle \eta =0} でx 1 {\displaystyle x_{1}} そしてx 2 {\displaystyle x_{2}} 定義により。また、前述のように、等式の左辺はゼロなので、
∫ x 1 x 2 η ( x ) ( ∂ L ∂ f − d d x ∂ L ∂ f ′ ) d x = 0 。 {\displaystyle \int _{x_{1}}^{x_{2}}\eta (x)\left({\frac {\partial L}{\partial f}}-{\frac {d}{dx}}{\frac {\partial L}{\partial f'}}\right)\,dx=0\,.}
変分法の基本補題 によれば、この方程式が任意の選択に対して成り立つという事実は、η {\displaystyle \eta } これは、括弧内の被積分関数の部分がゼロであることを意味します。
∂ L ∂ f − d d x ∂ L ∂ f ′ = 0 {\displaystyle {\frac {\partial L}{\partial f}}-{\frac {d}{dx}}{\frac {\partial L}{\partial f'}}=0}
これはオイラー・ラグランジュ方程式 と呼ばれる。この方程式の左辺は、の関数微分と呼ばれる。 J [ f ] {\displaystyle J[f]} そして、と表記されるδ J {\displaystyle \delta J} またはδ f ( x ) 。 {\displaystyle \delta f(x).}
一般に、これは2階常微分方程式 を与え、これを解くことで極値関数を求めることができる。f ( x ) 。 {\displaystyle f(x).} オイラー・ラグランジュ方程式は極値を求めるための必要 条件ではあるが、十分条件ではない。 J [ f ] 。 {\displaystyle J[f].} 最小値の十分条件は、「最小値の変動と十分条件」の セクションで説明されています。
例 このプロセスを説明するために、極値関数を求める問題を考えてみましょう。y = f ( x ) 、 {\displaystyle y=f(x),} これは2点を結ぶ最短の曲線です( x 1 、 y 1 ) {\displaystyle \left(x_{1},y_{1}\right)} そして( x 2 、 y 2 ) 。 {\displaystyle \left(x_{2},y_{2}\right).} 曲線の弧長は次のように表される。
A [ y ] = ∫ x 1 x 2 1 + [ y ′ ( x ) ] 2 d x 、 {\displaystyle A[y]=\int _{x_{1}}^{x_{2}}{\sqrt {1+[y'(x)]^{2}}}\,dx\,,}
と
y ′ ( x ) = d y d x 、 y 1 = f ( x 1 ) 、 y 2 = f ( x 2 ) 。 {\displaystyle y'(x)={\frac {dy}{dx}}\,,\ \ y_{1}=f(x_{1})\,,\ \ y_{2}=f(x_{2})\,.}
yが x の関数であると仮定すると一般性が失われることに注意してください。理想的には、両方とも他の何らかのパラメータの関数であるべきです。このアプローチは、あくまでも教育的な目的にのみ適しています。
オイラー・ラグランジュ方程式を用いて極値関数を求める。f ( x ) {\displaystyle f(x)} 機能を最小限に抑えるA [ y ] 。 {\displaystyle A[y].}
∂ L ∂ f − d d x ∂ L ∂ f ′ = 0 {\displaystyle {\frac {\partial L}{\partial f}}-{\frac {d}{dx}}{\frac {\partial L}{\partial f'}}=0}
と
L = 1 + [ f ′ ( x ) ] 2 。 {\displaystyle L={\sqrt {1+[f'(x)]^{2}}}\,.}
以来f {\displaystyle f} 明示的には現れないL 、 {\displaystyle L,} オイラー・ラグランジュ方程式の第1項は、すべての場合において消滅する。f ( x ) {\displaystyle f(x)} したがって、
d d x ∂ L ∂ f ′ = 0 。 {\displaystyle {\frac {d}{dx}}{\frac {\partial L}{\partial f'}}=0\,.}
代わりにL {\displaystyle L} そして微分すると、
d d x f ′ ( x ) 1 + [ f ′ ( x ) ] 2 = 0 。 {\displaystyle {\frac {d}{dx}}\ {\frac {f'(x)}{\sqrt {1+[f'(x)]^{2}}}}\ =0\,.}
したがって
f ′ ( x ) 1 + [ f ′ ( x ) ] 2 = c 、 {\displaystyle {\frac {f'(x)}{\sqrt {1+[f'(x)]^{2}}}}=c\,,}
ある定数に対してc {\displaystyle c} 。 それから
[ f ′ ( x ) ] 2 1 + [ f ′ ( x ) ] 2 = c 2 、 {\displaystyle {\frac {[f'(x)]^{2}}{1+[f'(x)]^{2}}}=c^{2}\,,}
どこ
0 ≤ c 2 < 1. {\displaystyle 0\leq c^{2}<1.}
解くと、
[ f ′ ( x ) ] 2 = c 2 1 − c 2 {\displaystyle [f'(x)]^{2}={\frac {c^{2}}{1-c^{2}}}}
これはつまり
f ′ ( x ) = m {\displaystyle f'(x)=m}
は定数であり、したがって2点を結ぶ最短曲線は( x 1 、 y 1 ) {\displaystyle \left(x_{1},y_{1}\right)} そして( x 2 、 y 2 ) {\displaystyle \left(x_{2},y_{2}\right)} は
f ( x ) = m x + b と m = y 2 − y 1 x 2 − x 1 そして b = x 2 y 1 − x 1 y 2 x 2 − x 1 {\displaystyle f(x)=mx+b\qquad {\text{with}}\ \ m={\frac {y_{2}-y_{1}}{x_{2}-x_{1}}}\quad {\text{and}}\quad b={\frac {x_{2}y_{1}-x_{1}y_{2}}{x_{2}-x_{1}}}}
こうして極値関数が見つかった。f ( x ) {\displaystyle f(x)} 機能を最小限に抑えるA [ y ] {\displaystyle A[y]} となることによってA [ f ] {\displaystyle A[f]} 最小値です。直線の式はy = m x + b 。 {\displaystyle y=mx+b.} つまり、2点間の最短距離は直線である。[ j ]
オイラー・ポアソン方程式もしS {\displaystyle S} の高階微分に依存するy ( x ) {\displaystyle y(x)} つまり、
S = ∫ 1 b f ( x 、 y ( x ) 、 y ′ ( x ) 、 … 、 y ( n ) ( x ) ) d x 、 {\displaystyle S=\int _{a}^{b}f(x,y(x),y'(x),\dots ,y^{(n)}(x))dx,}
それからy {\displaystyle y} オイラー・ポアソン 方程式を満たさなければならない。[ 22 ]
∂ f ∂ y − d d x ( ∂ f ∂ y ′ ) + ⋯ + ( − 1 ) n d n d x n [ ∂ f ∂ y ( n ) ] = 0. {\displaystyle {\frac {\partial f}{\partial y}}-{\frac {d}{dx}}\left({\frac {\partial f}{\partial y'}}\right)+\dots +(-1)^{n}{\frac {d^{n}}{dx^{n}}}\left[{\frac {\partial f}{\partial y^{(n)}}}\right]=0.}
デュ・ボワ=レイモンの定理これまでの議論では、極値関数は2つの連続導関数を持つと仮定してきたが、積分の存在はJ {\displaystyle J} 試行関数の1階微分のみを必要とする。極値で1階変分がゼロになるという条件は、オイラー・ラグランジュ方程式の弱形式 とみなすことができる。デュ・ボワ・レイモンの定理は、この弱形式が強形式を意味すると主張している。L {\displaystyle L} すべての引数に関して連続な1階および2階導関数を持ち、
∂ 2 L ∂ f ′ 2 ≠ 0 、 {\displaystyle {\frac {\partial ^{2}L}{\partial f'^{2}}}\neq 0,}
それからf {\displaystyle f} 2つの連続微分を持ち、オイラー・ラグランジュ方程式を満たす。
ラヴレンティエフ現象 ヒルベルトは、オイラー・ラグランジュ方程式が定常解を与えるための良好な条件を初めて示した人物である。凸領域と正の3回微分可能なラグランジアンの範囲内では、解は境界に沿って進むか、内部でオイラー・ラグランジュ方程式を満たす、可算個の断面から構成される。
しかし、 1926年にラヴレンティエフは 、最適解が存在しない場合でも、セクションの数を増やすことで任意に近づけることができる状況があることを示しました。ラヴレンティエフ現象は、許容関数の異なるクラス間で最小化問題の下限に違いがあることを示しています。たとえば、1934年にマニアによって提示された次の問題: [ 23 ]
L [ x ] = ∫ 0 1 ( x 3 − t ) 2 x ′ 6 、 {\displaystyle L[x]=\int _{0}^{1}(x^{3}-t)^{2}x'^{6},}
A = { x ∈ W 1 、 1 ( 0 、 1 ) : x ( 0 ) = 0 、 x ( 1 ) = 1 } 。 {\displaystyle {A}=\{x\in W^{1,1}(0,1):x(0)=0,\ x(1)=1\}.}
明らかに、x ( t ) = t 1 3 {\displaystyle x(t)=t^{\frac {1}{3}}} 関数を最小化しますが、任意の関数を見つけますx ∈ W 1 、 ∞ {\displaystyle x\in W^{1,\infty }} 下限値から離れた範囲の値を与える。
(一次元での)例は、従来はW 1 、 1 {\displaystyle W^{1,1}} そしてW 1 、 ∞ 、 {\displaystyle W^{1,\infty },} しかし、BallとMizel [ 24 ] は、Lavrentiev現象を広範囲にわたって示す最初の関数を取得しました。W 1 、 p {\displaystyle W^{1,p}} そしてW 1 、 q {\displaystyle W^{1,q}} のために1 ≤ p < q < ∞ 。 {\displaystyle 1\leq p<q<\infty .} この現象が発生しない基準を示す結果はいくつかあります。例えば、「標準的な成長」、第2変数に依存しないラグランジアン、またはチェザリの条件(D)を満たす近似シーケンスなどです。しかし、結果はしばしば特殊であり、少数の汎関数クラスにのみ適用可能です。
ラヴレンティエフ現象と関連しているのは反発特性です。ラヴレンティエフ現象を示す関数はすべて弱い反発特性を示します。[ 25 ]
複数の変数の関数 例えば、φ ( x 、 y ) {\displaystyle \varphi (x,y)} ドメイン上の膜の変位を示すD {\displaystyle D} でx 、 y {\displaystyle x,y} 平面の場合、その位置エネルギーは表面積に比例する。
U [ φ ] = ∬ D 1 + ∇ φ ⋅ ∇ φ d x d y 。 {\displaystyle U[\varphi ]=\iint _{D}{\sqrt {1+\nabla \varphi \cdot \nabla \varphi }}\,dx\,dy.}
プラトーの問題は、 境界上で所定の値を仮定しながら表面積を最小化する関数を見つけることである。D {\displaystyle D} 解は極小曲面 と呼ばれる。この問題に対するオイラー・ラグランジュ方程式は非線形である。
φ x x ( 1 + φ y 2 ) + φ y y ( 1 + φ x 2 ) − 2 φ x φ y φ x y = 0. {\displaystyle \varphi _{xx}(1+\varphi _{y}^{2})+\varphi _{yy}(1+\varphi _{x}^{2})-2\varphi _{x}\varphi _{y}\varphi _{xy}=0.}
詳細はクーラント(1950)を参照のこと。
ディリクレの原理膜の小さな変位のみを考慮すれば十分な場合が多く、変位がない場合とのエネルギー差は次のように近似される。
V [ φ ] = 1 2 ∬ D ∇ φ ⋅ ∇ φ d x d y 。 {\displaystyle V[\varphi ]={\frac {1}{2}}\iint _{D}\nabla \varphi \cdot \nabla \varphi \,dx\,dy.}
機能的V {\displaystyle V} すべての試行関数の中で最小化されるべきであるφ {\displaystyle \varphi } 境界上で規定値をとるD {\displaystyle D} 。 もしu {\displaystyle u} は最小化関数であり、v {\displaystyle v} は、境界上でゼロになる任意の滑らかな関数です。D {\displaystyle D} 、そして最初のバリエーションV [ u + ε v ] {\displaystyle V[u+\varepsilon v]} 消え去らなければならない:
d d ε V [ u + ε v ] | ε = 0 = ∬ D ∇ u ⋅ ∇ v d x d y = 0. {\displaystyle \left.{\frac {d}{d\varepsilon }}V[u+\varepsilon v]\right|_{\varepsilon =0}=\iint _{D}\nabla u\cdot \nabla v\,dx\,dy=0.}
ただし、u {\displaystyle u} 2つの導関数を持つので、発散定理を適用して
∬ D ∇ ⋅ ( v ∇ u ) d x d y = ∬ D ∇ u ⋅ ∇ v + v ∇ ⋅ ∇ u d x d y = ∫ C v ∂ u ∂ n d s 、 {\displaystyle \iint _{D}\nabla \cdot (v\nabla u)\,dx\,dy=\iint _{D}\nabla u\cdot \nabla v+v\nabla \cdot \nabla u\,dx\,dy=\int _{C}v{\frac {\partial u}{\partial n}}\,ds,}
どこC {\displaystyle C} 境界はD 、 {\displaystyle D,} s {\displaystyle s} 弧長はC {\displaystyle C} そして∂ u / ∂ n {\displaystyle \partial u/\partial n} は、u {\displaystyle u} の上C 。 {\displaystyle C.} 以来v {\displaystyle v} 消えるC {\displaystyle C} そして最初の変分が消えると、結果は
∬ D v ∇ ⋅ ∇ u d x d y = 0 {\displaystyle \iint _{D}v\nabla \cdot \nabla u\,dx\,dy=0}
すべての滑らかな関数についてv {\displaystyle v} 境界で消えるD {\displaystyle D} 1次元積分の場合の証明をこの場合にも適用すると、次のことが示される。
∇ ⋅ ∇ u = 0 {\displaystyle \nabla \cdot \nabla u=0} でD 。 {\displaystyle D.}
この推論の難点は、最小化関数がu {\displaystyle u} 2つの導関数を持つ必要がある。リーマンは、滑らかな最小化関数の存在は物理的問題との関連によって保証されると主張した。膜は実際に最小ポテンシャルエネルギーを持つ配置をとる。リーマンはこの考えを、師であるペーター・グスタフ・ルジューヌ・ディリクレにちなんで ディリクレ原理 と名付けた。しかし、ワイエルシュトラスは、解のない変分問題の例を挙げた。最小化
W [ φ ] = ∫ − 1 1 ( x φ ′ ) 2 d x {\displaystyle W[\varphi ]=\int _{-1}^{1}(x\varphi ')^{2}\,dx}
すべての機能の中でφ {\displaystyle \varphi } 満足するφ ( − 1 ) = − 1 {\displaystyle \varphi (-1)=-1} そしてφ ( 1 ) = 1. {\displaystyle \varphi (1)=1.} W {\displaystyle W} は、原点の小さな近傍で −1 から 1 の間を遷移する区分的線形関数を選択することによって、任意に小さくすることができます。しかし、 となるような関数は存在しません。W = 0. {\displaystyle W=0.} [ k ] 最終的にディリクレの原理が有効であることが示されましたが、楕円型偏微分方程式 。Jost と Li–Jost (1998) を参照してください。
他の境界値問題への一般化 膜のポテンシャルエネルギーのより一般的な表現は次のとおりである。
V [ φ ] = ∬ D [ 1 2 ∇ φ ⋅ ∇ φ + f ( x 、 y ) φ ] d x d y + ∫ C [ 1 2 σ ( s ) φ 2 + g ( s ) φ ] d s 。 {\displaystyle V[\varphi ]=\iint _{D}\left[{\frac {1}{2}}\nabla \varphi \cdot \nabla \varphi +f(x,y)\varphi \right]\,dx\,dy\,+\int _{C}\left[{\frac {1}{2}}\sigma (s)\varphi ^{2}+g(s)\varphi \right]\,ds.}
これは外部力密度に相当するf ( x 、 y ) {\displaystyle f(x,y)} でD 、 {\displaystyle D,} 外部の力g ( s ) {\displaystyle g(s)} 境界線上C 、 {\displaystyle C,} 弾性力σ ( s ) {\displaystyle \sigma (s)} 行動するC {\displaystyle C} 境界値に制限を設けずに ポテンシャルエネルギーを最小化する関数は、次のように表される。u {\displaystyle u} ただし、f {\displaystyle f} そしてg {\displaystyle g} 連続であり、正則性理論は最小化関数を意味する。u {\displaystyle u} 2つの導関数を持つ。最初の変分をとる場合、増分に境界条件を課す必要はない。v {\displaystyle v} .最初のバリエーションV [ u + ε v ] {\displaystyle V[u+\varepsilon v]} は
∬ D [ ∇ u ⋅ ∇ v + f v ] d x d y + ∫ C [ σ u v + g v ] d s = 0. {\displaystyle \iint _{D}\left[\nabla u\cdot \nabla v+fv\right]\,dx\,dy+\int _{C}\left[\sigma uv+gv\right]\,ds=0.}
発散定理を適用すると、結果は
∬ D [ − v ∇ ⋅ ∇ u + v f ] d x d y + ∫ C v [ ∂ u ∂ n + σ u + g ] d s = 0. {\displaystyle \iint _{D}\left[-v\nabla \cdot \nabla u+vf\right]\,dx\,dy+\int _{C}v\left[{\frac {\partial u}{\partial n}}+\sigma u+g\right]\,ds=0.}
まず最初に設定するとv = 0 {\displaystyle v=0} の上C 、 {\displaystyle C,} 境界積分はゼロになり、前と同様に次の結論に至る。
− ∇ ⋅ ∇ u + f = 0 {\displaystyle -\nabla \cdot \nabla u+f=0}
でD {\displaystyle D} それで、もし私たちが許可すればv {\displaystyle v} 任意の境界値を仮定すると、これは次のことを意味します。u {\displaystyle u} 境界条件を満たさなければならない
∂ u ∂ n + σ u + g = 0 、 {\displaystyle {\frac {\partial u}{\partial n}}+\sigma u+g=0,}
の上C {\displaystyle C} この境界条件は、最小化特性の結果である。u {\displaystyle u} : それは事前に課せられるものではありません。このような条件は自然境界条件 と呼ばれます。
上記の推論は、以下の場合には妥当ではない。σ {\displaystyle \sigma } 同じように消えるC 。 {\displaystyle C.} このような場合、試用機能を許可することができます。φ ≡ c {\displaystyle \varphi \equiv c} 、 どこc {\displaystyle c} は定数です。このような試行関数では、
V [ c ] = c [ ∬ D f d x d y + ∫ C g d s ] 。 {\displaystyle V[c]=c\left[\iint _{D}f\,dx\,dy+\int _{C}g\,ds\right].}
適切な選択によりc {\displaystyle c} 、V {\displaystyle V} 括弧内の量がゼロにならない限り、任意の値をとることができます。したがって、変分問題は、次の条件を満たさない限り無意味です。
∬ D f d x d y + ∫ C g d s = 0. {\displaystyle \iint _{D}f\,dx\,dy+\int _{C}g\,ds=0.}
この条件は、系にかかる正味の外力が平衡状態にあることを意味します。これらの力が平衡状態にある場合、変分問題には解が存在しますが、任意の定数が加えられる可能性があるため、解は一意ではありません。詳細と例については、Courant and Hilbert (1953) を参照してください。
固有値問題 一次元および多次元の固有値問題は 、いずれも変分問題として定式化できる。
シュトゥルム=リウヴィル問題シュトゥルム・リウヴィル固有値問題は、 一般的な二次形式を含む。
Q [ y ] = ∫ x 1 x 2 [ p ( x ) y ′ ( x ) 2 + q ( x ) y ( x ) 2 ] d x 、 {\displaystyle Q[y]=\int _{x_{1}}^{x_{2}}\left[p(x)y'(x)^{2}+q(x)y(x)^{2}\right]\,dx,}
どこy {\displaystyle y} 境界条件を満たす関数に限定される
y ( x 1 ) = 0 、 y ( x 2 ) = 0. {\displaystyle y(x_{1})=0,\quad y(x_{2})=0.}
させてR {\displaystyle R} 正規化積分
R [ y ] = ∫ x 1 x 2 r ( x ) y ( x ) 2 d x 。 {\displaystyle R[y]=\int _{x_{1}}^{x_{2}}r(x)y(x)^{2}\,dx.}
機能p ( x ) {\displaystyle p(x)} そしてr ( x ) {\displaystyle r(x)} は、あらゆる場所で正であり、ゼロから離れた範囲で有界である必要があります。主な変分問題は、比率を最小化することです。Q / R {\displaystyle Q/R} すべての中でy {\displaystyle y} 終点条件を満たすこと、これは最小化することと同等であるQ [ y ] {\displaystyle Q[y]} 制約条件の下でR [ y ] {\displaystyle R[y]} は定数である。以下に示すように、最小化のためのオイラー・ラグランジュ方程式はu {\displaystyle u} は
− ( p u ′ ) ′ + q u − λ r u = 0 、 {\displaystyle -(pu')'+qu-\lambda ru=0,}
どこλ {\displaystyle \lambda } 商は
λ = Q [ u ] R [ u ] 。 {\displaystyle \lambda ={\frac {Q[u]}{R[u]}}.}
最小化は、(GelfandとFomin 1963を参照)示すことができる。u {\displaystyle u} 2つの導関数を持ち、オイラー・ラグランジュ方程式を満たす。関連するλ {\displaystyle \lambda } は、λ 1 {\displaystyle \lambda _{1}} これは、この方程式と境界条件に対する最小の固有値です。関連する最小化関数は、で表されます。u 1 ( x ) {\displaystyle u_{1}(x)} 固有値のこの変分特性はレイリー・リッツ法 につながる。近似値を選択するu {\displaystyle u} 基底関数(例えば三角関数)の線形結合として表し、そのような線形結合の中から有限次元の最小化を行う。この方法は、しばしば驚くほど高い精度を示す。
次に小さい固有値と固有関数は、最小化することによって得られます。Q {\displaystyle Q} 追加の制約の下で
∫ x 1 x 2 r ( x ) u 1 ( x ) y ( x ) d x = 0. {\displaystyle \int _{x_{1}}^{x_{2}}r(x)u_{1}(x)y(x)\,dx=0.}
この手順を拡張することで、問題の固有値と固有関数の完全なシーケンスを取得できます。
変分問題は、より一般的な境界条件にも適用されます。y {\displaystyle y} 端点で消滅し、端点で条件を課すことができないため、
Q [ y ] = ∫ x 1 x 2 [ p ( x ) y ′ ( x ) 2 + q ( x ) y ( x ) 2 ] d x + 1 1 y ( x 1 ) 2 + 1 2 y ( x 2 ) 2 、 {\displaystyle Q[y]=\int _{x_{1}}^{x_{2}}\left[p(x)y'(x)^{2}+q(x)y(x)^{2}\right]\,dx+a_{1}y(x_{1})^{2}+a_{2}y(x_{2})^{2},}
どこ1 1 {\displaystyle a_{1}} そして1 2 {\displaystyle a_{2}} は任意です。y = u + ε v {\displaystyle y=u+\varepsilon v} 比率の最初のバリエーションQ / R {\displaystyle Q/R} は
V 1 = 2 R [ u ] ( ∫ x 1 x 2 [ p ( x ) u ′ ( x ) v ′ ( x ) + q ( x ) u ( x ) v ( x ) − λ r ( x ) u ( x ) v ( x ) ] d x + 1 1 u ( x 1 ) v ( x 1 ) + 1 2 u ( x 2 ) v ( x 2 ) ) 、 {\displaystyle V_{1}={\frac {2}{R[u]}}\left(\int _{x_{1}}^{x_{2}}\left[p(x)u'(x)v'(x)+q(x)u(x)v(x)-\lambda r(x)u(x)v(x)\right]\,dx+a_{1}u(x_{1})v(x_{1})+a_{2}u(x_{2})v(x_{2})\right),}
どこλ {\displaystyle \lambda } 比率によって与えられるQ [ u ] / R [ u ] {\displaystyle Q[u]/R[u]} 前と同様。部分的な統合後、
R [ u ] 2 V 1 = ∫ x 1 x 2 v ( x ) [ − ( p u ′ ) ′ + q u − λ r u ] d x + v ( x 1 ) [ − p ( x 1 ) u ′ ( x 1 ) + 1 1 u ( x 1 ) ] + v ( x 2 ) [ p ( x 2 ) u ′ ( x 2 ) + 1 2 u ( x 2 ) ] 。 {\displaystyle {\frac {R[u]}{2}}V_{1}=\int _{x_{1}}^{x_{2}}v(x)\left[-(pu')'+qu-\lambda ru\right]\,dx+v(x_{1})[-p(x_{1})u'(x_{1})+a_{1}u(x_{1})]+v(x_{2})[p(x_{2})u'(x_{2})+a_{2}u(x_{2})].}
まず最初に要求するのはv {\displaystyle v} 端点で消滅し、最初の変分はすべてのそのようなものに対して消滅しますv {\displaystyle v} のみ
− ( p u ′ ) ′ + q u − λ r u = 0 のために x 1 < x < x 2 。 {\displaystyle -(pu')'+qu-\lambda ru=0\quad {\hbox{for}}\quad x_{1}<x<x_{2}.}
もしu {\displaystyle u} この条件を満たす場合、最初の変分は任意のに対して消滅する。v {\displaystyle v} のみ
− p ( x 1 ) u ′ ( x 1 ) + 1 1 u ( x 1 ) = 0 、 そして p ( x 2 ) u ′ ( x 2 ) + 1 2 u ( x 2 ) = 0. {\displaystyle -p(x_{1})u'(x_{1})+a_{1}u(x_{1})=0,\quad {\hbox{and}}\quad p(x_{2})u'(x_{2})+a_{2}u(x_{2})=0.}
これらの条件は、この問題における自然な境界条件 である。なぜなら、それらは最小化のための試行関数に課されるものではなく、最小化の結果として生じるものだからである。
多次元における固有値問題 高次元における固有値問題は、一次元の場合と同様に定義される。例えば、領域が与えられた場合D {\displaystyle D} 境界付きB {\displaystyle B} 3次元では次のように定義できます
Q [ φ ] = ∭ D p ( X ) ∇ φ ⋅ ∇ φ + q ( X ) φ 2 d x d y d z + ∬ B σ ( S ) φ 2 d S 、 {\displaystyle Q[\varphi ]=\iiint _{D}p(X)\nabla \varphi \cdot \nabla \varphi +q(X)\varphi ^{2}\,dx\,dy\,dz+\iint _{B}\sigma (S)\varphi ^{2}\,dS,}
そして
R [ φ ] = ∭ D r ( X ) φ ( X ) 2 d x d y d z 。 {\displaystyle R[\varphi ]=\iiint _{D}r(X)\varphi (X)^{2}\,dx\,dy\,dz.}
させてu {\displaystyle u} 商を最小化する関数とするQ [ φ ] / R [ φ ] {\displaystyle Q[\varphi ]/R[\varphi ]} 境界に関する条件は規定されていないB 。 {\displaystyle B.} オイラー・ラグランジュ方程式を満たすu {\displaystyle u} は
− ∇ ⋅ ( p ( X ) ∇ u ) + q ( x ) u − λ r ( x ) u = 0 、 {\displaystyle -\nabla \cdot (p(X)\nabla u)+q(x)u-\lambda r(x)u=0,}
どこ
λ = Q [ u ] R [ u ] 。 {\displaystyle \lambda ={\frac {Q[u]}{R[u]}}.}
最小化u {\displaystyle u} 自然境界条件も満たさなければならない
p ( S ) ∂ u ∂ n + σ ( S ) u = 0 、 {\displaystyle p(S){\frac {\partial u}{\partial n}}+\sigma (S)u=0,}
境界線上B 。 {\displaystyle B.} この結果は、楕円型偏微分方程式の正則性理論に基づいています。詳細は、Jost および Li–Jost (1998) を参照してください。完全性に関する結果、固有値の漸近的性質、固有関数の節点に関する結果など、多くの拡張は Courant および Hilbert (1953) に記載されています。
アプリケーション
光学 フェルマーの原理 によれば、光は(局所的に)その両端間の光路長を最小にする経路を通ります。x {\displaystyle x} -座標がパスに沿ったパラメータとして選択され、y = f ( x ) {\displaystyle y=f(x)} 経路に沿って、光路長は次のように与えられる。
A [ f ] = ∫ x 0 x 1 n ( x 、 f ( x ) ) 1 + f ′ ( x ) 2 d x 、 {\displaystyle A[f]=\int _{x_{0}}^{x_{1}}n(x,f(x)){\sqrt {1+f'(x)^{2}}}dx,}
屈折率n ( x 、 y ) {\displaystyle n(x,y)} 素材によります。f ( x ) = f 0 ( x ) + ε f 1 ( x ) {\displaystyle f(x)=f_{0}(x)+\varepsilon f_{1}(x)} 次に最初のバリエーション A {\displaystyle A} (A {\displaystyle A} に関してε {\displaystyle \varepsilon } ) は
δ A [ f 0 、 f 1 ] = ∫ x 0 x 1 [ n ( x 、 f 0 ) f 0 ′ ( x ) f 1 ′ ( x ) 1 + f 0 ′ ( x ) 2 + n y ( x 、 f 0 ) f 1 1 + f 0 ′ ( x ) 2 ] d x 。 {\displaystyle \delta A[f_{0},f_{1}]=\int _{x_{0}}^{x_{1}}\left[{\frac {n(x,f_{0})f_{0}'(x)f_{1}'(x)}{\sqrt {1+f_{0}'(x)^{2}}}}+n_{y}(x,f_{0})f_{1}{\sqrt {1+f_{0}'(x)^{2}}}\right]dx.}
括弧内の第1項を部分積分すると、オイラー・ラグランジュ方程式が得られる。
− d d x [ n ( x 、 f 0 ) f 0 ′ 1 + f 0 ′ 2 ] + n y ( x 、 f 0 ) 1 + f 0 ′ ( x ) 2 = 0. {\displaystyle -{\frac {d}{dx}}\left[{\frac {n(x,f_{0})f_{0}'}{\sqrt {1+f_{0}'^{2}}}}\right]+n_{y}(x,f_{0}){\sqrt {1+f_{0}'(x)^{2}}}=0.}
光線は、この方程式を積分することによって求めることができる。この形式は、ラグランジュ光学 およびハミルトン光学 の分野で用いられる。
3次元におけるフェルマーの原理ベクトル表記を用いるのが都合が良い。X = ( x 1 、 x 2 、 x 3 ) 、 {\displaystyle X=(x_{1},x_{2},x_{3}),} させてt {\displaystyle t} をパラメータとする。X ( t ) {\displaystyle X(t)} 曲線のパラメトリック表現であるC 、 {\displaystyle C,} そしてX ˙ ( t ) {\displaystyle {\dot {X}}(t)} をその接線ベクトルとする。曲線の光学的長さは次式で与えられる。
A [ C ] = ∫ t 0 t 1 n ( X ) X ˙ ⋅ X ˙ d t 。 {\displaystyle A[C]=\int _{t_{0}}^{t_{1}}n(X){\sqrt {{\dot {X}}\cdot {\dot {X}}}}\,dt.}
この積分は、パラメトリック表現の変化に対して不変であることに注意してください。C 。 {\displaystyle C.} 最小化曲線に対するオイラー・ラグランジュ方程式は対称形である
d d t P = X ˙ ⋅ X ˙ ∇ n 、 {\displaystyle {\frac {d}{dt}}P={\sqrt {{\dot {X}}\cdot {\dot {X}}}}\,\nabla n,}
どこ
P = n ( X ) X ˙ X ˙ ⋅ X ˙ 。 {\displaystyle P={\frac {n(X){\dot {X}}}{\sqrt {{\dot {X}}\cdot {\dot {X}}}}}.}
定義から次のことが導かれる。P {\displaystyle P} 満たす
P ⋅ P = n ( X ) 2 。 {\displaystyle P\cdot P=n(X)^{2}.}
したがって、積分は次のように書くこともできます。
A [ C ] = ∫ t 0 t 1 P ⋅ X ˙ d t 。 {\displaystyle A[C]=\int _{t_{0}}^{t_{1}}P\cdot {\dot {X}}\,dt.}
この形式は、関数を見つけることができればψ {\displaystyle \psi } その勾配は次のように与えられる。P 、 {\displaystyle P,} すると積分A {\displaystyle A} は、ψ {\displaystyle \psi } 積分区間の端点において。したがって、積分を停留させる曲線の研究の問題は、等位面の研究と関連付けることができる。ψ {\displaystyle \psi } このような関数を見つけるために、光の伝播を支配する波動方程式に着目します。この形式は、ラグランジュ光学 とハミルトン光学の 文脈で使用されます。
波動方程式との関連性 不均質媒体の波動方程式は
u t t = c 2 ∇ ⋅ ∇ u 、 {\displaystyle u_{tt}=c^{2}\nabla \cdot \nabla u,}
どこc {\displaystyle c} は速度であり、一般的にはX {\displaystyle X} 光の波面は、この偏微分方程式の特性面であり、次の式を満たす。
φ t 2 = c ( X ) 2 ∇ φ ⋅ ∇ φ 。 {\displaystyle \varphi _{t}^{2}=c(X)^{2}\,\nabla \varphi \cdot \nabla \varphi .}
私たちは次のような形で解決策を探すかもしれません
φ ( t 、 X ) = t − ψ ( X ) 。 {\displaystyle \varphi (t,X)=t-\psi (X).}
その場合、ψ {\displaystyle \psi } 満たす
∇ ψ ⋅ ∇ ψ = n 2 、 {\displaystyle \nabla \psi \cdot \nabla \psi =n^{2},}
どこn = 1 / c {\displaystyle n=1/c} 1階偏微分方程式 の理論によれば、P = ∇ ψ 、 {\displaystyle P=\nabla \psi ,} それからP {\displaystyle P} 満たす
d P d s = n ∇ n 、 {\displaystyle {\frac {dP}{ds}}=n\,\nabla n,}
曲線系(光線 )に沿って、
d X d s = P 。 {\displaystyle {\frac {dX}{ds}}=P.}
1階偏微分方程式の解を求めるこれらの式は、次のことを行えばオイラー・ラグランジュ方程式と同一になる。
d s d t = X ˙ ⋅ X ˙ n 。 {\displaystyle {\frac {ds}{dt}}={\frac {\sqrt {{\dot {X}}\cdot {\dot {X}}}}{n}}.}
関数はψ {\displaystyle \psi } 最小化積分の値A {\displaystyle A} 上端点の関数として。つまり、最小化曲線の族を構築すると、光路長の値は波動方程式に対応する特性方程式を満たす。したがって、関連する1階偏微分方程式を解くことは、変分問題の解の族を見つけることと同等である。これがハミルトン・ヤコビ理論 の本質的な内容であり、より一般的な変分問題にも適用できる。
力学 古典力学では、作用、S 、 {\displaystyle S,} ラグランジアンの時間積分として定義され、L {\displaystyle L} ラグランジアンはエネルギーの差であり、
L = T − U 、 {\displaystyle L=T-U,}
どこT {\displaystyle T} は機械システムの運動エネルギー であり、U {\displaystyle U} その位置エネルギー 。ハミルトンの原理 (または作用原理)は、保存的なホロノミック(積分可能な拘束)力学系の運動は、作用積分が
S = ∫ t 0 t 1 L ( x 、 x ˙ 、 t ) d t {\displaystyle S=\int _{t_{0}}^{t_{1}}L(x,{\dot {x}},t)\,dt}
経路の変化に対して静止しているx ( t ) {\displaystyle x(t)} この系に対するオイラー・ラグランジュ方程式は、ラグランジュ方程式として知られています。
d d t ∂ L ∂ x ˙ = ∂ L ∂ x 、 {\displaystyle {\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {x}}}}={\frac {\partial L}{\partial x}},}
そしてそれらは(そのような系における)ニュートンの運動方程式と等価である。
共役運動量P {\displaystyle P} 定義される
p = ∂ L ∂ x ˙ 。 {\displaystyle p={\frac {\partial L}{\partial {\dot {x}}}}.}
例えば、
T = 1 2 m x ˙ 2 、 {\displaystyle T={\frac {1}{2}}m{\dot {x}}^{2},}
それからp = m x ˙ 。 {\displaystyle p=m{\dot {x}}.}
ハミルトン力学の 結果は、共役運動量が代わりに導入された場合にx ˙ {\displaystyle {\dot {x}}} ラグランジアンのルジャンドル変換によってL {\displaystyle L} ハミルトニアンへH {\displaystyle H} 定義される
H ( x 、 p 、 t ) = p x ˙ − L ( x 、 x ˙ 、 t ) 。 {\displaystyle H(x,p,t)=p\,{\dot {x}}-L(x,{\dot {x}},t).}
ハミルトニアンは系の全エネルギーを表す。H = T + U {\displaystyle H=T+U} フェルマーの原理との類推から、ラグランジュ方程式の解(粒子の軌跡)は、ある関数の等位面によって記述できることが示唆される。X {\displaystyle X} この関数はハミルトン・ヤコビ方程式 の解です。
∂ ψ ∂ t + H ( x 、 ∂ ψ ∂ x 、 t ) = 0. {\displaystyle {\frac {\partial \psi }{\partial t}}+H\left(x,{\frac {\partial \psi }{\partial x}},t\right)=0.}
さらなる応用 変分法のその他の応用例としては、以下のようなものがある。
参考文献 1 2 クーラント& ヒルベルト 1953 、p. 1841 2 ゴールドスタイン、ハーマン H. (1980). 17 世紀から 19 世紀までの変分法の歴史 . スプリンガー ニューヨーク. pp. 7–21 . ISBN 978-1-4613-8106-8 。 1 2 Ferguson, James (2004), A Brief Survey of the History of the Calculus of Variations and its Applications , arXiv : math/0402357 , Bibcode : 2004math......2357F 1 2 ローランズ、ピーター (2017). ニュートンと偉大な世界システム . ワールドサイエンティフィック出版 . pp. 36–39 . doi : 10.1142/q0108 . ISBN 978-1-78634-372-7 。↑ Torres, Delfim FM (2021-07-29). "On a Non-Newtonian Calculus of Variations" . Axioms . 10 (3): 171. arXiv : 2107.14152 . doi : 10.3390/axioms10030171 . ISSN 2075-1680 . ↑ ゲルファンド、IM ; フォミン、SV (2000)。シルバーマン、リチャード A. (編)。 変分法 (完全復刻版)。ミネオラ 、 ニューヨーク:ドーバー出版。p. 3。ISBN 978-0-486-41448-5 。1 2 Thiele, Rüdiger (2007). "Euler and the Calculus of Variations" .Bradley, Robert E.; Sandifer, C. Edward (eds.). Leonhard Euler: Life, Work and Legacy . Elsevier. p. 249. ISBN 978-0-08-047129-7 。↑ ゴールドスタイン、ハーマン H. (2012). 17 世紀から 19 世紀までの変分法の歴史 . Springer Science & Business Media. p. 110. ISBN 978-1-4613-8106-8 。1 2 3 van Brunt, Bruce (2004). 変分法 . Springer. ISBN 978-0-387-40247-5 。1 2 Ferguson, James (2004). "変分法の歴史とその応用に関する簡単な概観". arXiv : math/0402357 . ↑ ディミトリ・ベルツェカス 。『動的計画法と最適制御』アテナ・サイエンティフィック、2005年。↑ Bellman, Richard E. (1954). "動的計画法と変分法の新しい形式体系" . Proc. Natl. Acad. Sci . 40 (4): 231– 235. Bibcode : 1954PNAS...40..231B . doi : 10.1073/pnas.40.4.231 . PMC 527981 . PMID 16589462 . ↑ 「リチャード・E・ベルマン制御遺産賞」 。 アメリカ自動制御評議会 。2004年。 2018年10月1日に オリジナル からアーカイブ 。 2013年7月28日 に取得。 ↑ クーラント、R ; ヒルベルト、D (1953). 数理物理学の方法 . 第 I 巻 (初版 ). ニューヨーク: インターサイエンス出版社, Inc. p. 169. ISBN 978-0-471-50447-4 。↑ ゲルファンド& フォミン 2000、12 ~13 ページ ↑ ゲルファンド& フォミン 2000 、 p.13 ↑ ゲルファンド& フォミン 2000、14 ~15 ページ ↑ クーラント、R. ; ヒルベルト、D. (1953). 数理物理学の方法 . 第 I 巻 (初版 ). ニューヨーク: インターサイエンス・パブリッシャーズ, Inc. ISBN 978-0-471-50447-4 。↑ ケランド、フィリップ (1843)。 証明数学の原理に関する講義 。p. 58 – Google ブックス経由。 ↑ ワイスタイン、エリック W. 「オイラー・ラグランジュ微分方程式」 。 mathworld.wolfram.com 。ヴォルフラム。方程式 (5)。 ↑ Montesinos, Merced; Gonzalez, D.; Meza, V. (2026). "対称性とヘルムホルツ条件を組み合わせてラグランジアンを構築する" . Advances in Mathematical Physics . 2026 (1) 9534805. arXiv : 2603.07760 . doi : 10.1155/admp/9534805 . ↑ Kot, Mark (2014). 「第4章:基本的な一般化」. 変分法入門 . アメリカ数学会. ISBN 978-1-4704-1495-5 。↑ マニア、バーナード (1934)。 「ソプラ・ウン・エセンピオ・ディ・ラヴレンチェフ」。 ボレンティーノ・デル・ウニオーネ・マテマティカ・イタリアーナ 。 13 : 147–153 . ↑ Ball & Mizel (1985). "最小化解がオイラー・ラグランジュ方程式を満たさない一次元変分問題". Archive for Rational Mechanics and Analysis . 90 (4): 325–388 . Bibcode : 1985ArRMA..90..325B . doi : 10.1007/BF00276295 . S2CID 55005550 . ↑ フェリエロ、アレッサンドロ (2007)。 「弱反発特性」。 Journal de Mathématiques Pures et Appliquées 。 88 (4): 378–388 。 土井 : 10.1016/j.matpur.2007.06.002 。 ↑ ターンブル。 「リーマン伝記」 。英国:セント・アンドリュー大学。 ↑ ゲルファンド& フォミン 2000、11-12 ページ 、99ページ ↑ ゲルファンド& フォミン 2000、12 ページ 、脚注6↑ ゲルファンド& フォミン 2000 、p.8↑ ゲルファンド& フォミン 2000 、p.6↑ ゲルファンド& フォミン 2000、11 ~12 ページ ↑ ゲルファンド& フォミン 2000、97 ~98ページ↑ ゲルファンド& フォミン 2000、99 ページ ↑ ゲルファンド& フォミン 2000 、 p.100 ↑ ゲルファント& フォミン 2000 、p.100 、定理2
さらに読む Benesova, B. および Kruzik, M.: 「積分汎関数の弱い下半連続性と応用」 SIAM Review 59(4) (2017), 703–766. ボルツァ、O .:変分法の講義。チェルシー出版、1904年、デジタル数学ライブラリで入手可能。第2版は1961年に再版、ペーパーバックは2005年に出版、ISBN 978-1-4181-8201-4 。Cassel, Kevin W.: Variational Methods with Applications in Science and Engineering、ケンブリッジ大学出版局、2013年。 Clegg, JC:変分法、Interscience Publishers Inc.、1968年。 クーラント、R .:ディリクレの原理、等角写像、極小曲面。インターサイエンス、1950年。ダコローニャ、ベルナール :「序論」変分法入門 、第3版、2014年、ワールド・サイエンティフィック・パブリッシング、ISBN 978-1-78326-551-0 。Elsgolc, LE:変分法、Pergamon Press Ltd.、1962年。 Forsyth, AR:変分法、ドーバー出版、1960年。 Fox, Charles: 『変分法入門』、ドーバー出版、1987年。 ジャキンタ、マリアーノ。ステファン、ヒルデブラント: 変分計算 I および II、シュプリンガー・フェルラーク、ISBN 978-3-662-03278-7 ISBN 978-3-662-06201-2 Jost, J. および X. Li-Jost:変分法. Cambridge University Press, 1998.ランチョス、コルネリウス:力学の変分原理(アルバート・アインシュタインに捧ぐ)、トロント大学出版局、ISBN 0-8020-1743-6 1962年、1966年、1970年版がそれに続いた。ISBN 0-486-65067-7 Lebedev, LP および Cloud, MJ:変分法と最適制御および力学への応用を伴う関数解析、World Scientific、2003 年、1-98 ページ。 ローガン、J.デイビッド:応用数学、第3版。ワイリー・インターサイエンス、2006年 パイク、ラルフ W. 「第 8 章:変分法」。『エンジニアリング システムのための最適化 』 。ルイジアナ州立大学 。2007年 7 月 5 日のオリジナルからアーカイブ済み。 Roubicek, T.: 「変分法」。 『物理学者のための数学ツール』 第17章。(M. Grinfeld編)J. Wiley、Weinheim、2014年、ISBN 978-3-527-41188-7 、551~588ページ 。 ハンス・サガン著『変分法入門』、ドーバー出版、1992年。 ワインストック、ロバート:物理学と工学への応用を伴う変分法、ドーバー出版、1974年(1952年版の復刻版)。
外部リンク 変分法。数学百科事典 。 変分法。PlanetMath 。 変分法 。MathWorld 。 変分法。例題。 数学 - 変分法と積分方程式。YouTube での講義 。 測地線場に関する選集。パートI、パートII。