=====================================
変分法2
=====================================
変分法1_ では、関数 $f(x,y,y')$ によって決まる汎関数 $I[y]\equiv \int _{a}^{b}f(x,y,y')dx$ の停留値問題を解くためには、オイラー方程式を解けば良いのだということを勉強しました。
\displaystyle {\partial f\over \partial y}-{d\over dx}\Big({\partial f\over \partial y'}\Big)=0 \tag{0}
ここでは、オイラー方程式の、さらに発展的な形について勉強します。もし、なんだかちょっと難しいぞ、と感じる人は、もう一度 変分法1_ を復習して下さい。
色々な関数形が出てきますが、その全てを理解する必要はありません。全部読むのは長くて大変ですから、必要なときに必要な部分だけ、この記事を参考書代わりに読んで頂ければ十分だと思います。次のセクションが目次を兼ねています。
色々な関数形
-----------------------------------------
汎関数 $I[y]$ を決める関数として、 変分法1_ では $f(x,y,y')$ という形のものしか考えませんでした。しかし、言うまでもなく世の中には、もっと色々な形の関数があります。これよりももっと簡単な形の関数もあれば、もっと複雑なものもあるわけです。次に挙げるような関数に対して、オイラー方程式( 変分法1_ で導いたオイラー方程式と区別するために、これらを変形オイラー方程式と呼びます)を導くのがこの記事の目的です。
.. contents::
1.〜3.は、 $f$ が、特に簡単な関数形をしている場合です。4.は、 $y$ 以外にも変数を含む場合です。 $z$ も $x$ の関数であることに注意して下さい。これは一般に $n$ 個の変数を含む場合に容易に拡張できますので、本文の方でその場合にも触れます。5.は $f$ が $y$ の二次導関数まで含む場合です。5.の場合を更に一般の高次導関数まで拡張したものが6.です。5.と6.は本質的に同じ問題ですが、少し計算が長くなるので別にしました。1.〜6.に出てくる $y$ や $z$ は $x$ の関数でしたから、基本変数は $x$ だけでした。7.として、独立に動ける基本変数が $x_{1},x_{2}$ のように複数個あり、それらの関数である $u=u(x_{1},x_{2})$ と、その偏導関数 $u_{x_{1}}$ . $u_{x_{2}}$ が $f$ に含まれる場合を考えます。
以下に、この7つのタイプの変分問題について、順番に変形オイラー方程式を導いていきます。本文では一応導出の過程も全て述べますが、これらを覚える必要は全くありません。最後に幾つかの例題を載せておきますので、出来れば答えを見ないでチャレンジしてみてください。
色々な変形オイラー方程式
---------------------------------------------------------------------
1. $f=f(x,y')$
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
これは関数 $f$ に $y$ が含まれていない簡単な形です。 $\frac{\partial f}{\partial y}$ が零になりますので、オイラー方程式(0)は次式のようになります。
\frac{d}{dx}\Big( {\partial f\over \partial y'} \Big) = 0
両辺を $x$ で積分して次のような形になります。これが1番目の形の変形オイラー方程式です。
\displaystyle {\partial f\over \partial y'}
=const. \tag{1}
とても計算が簡単だったので、なにもわざわざ勿体ぶって公式などと呼ぶ必要はないと感じた読者の方が多いと思います。しかし、実際に変分法の問題を解くときには少しでも簡単な形から始めるということが、とても重要なのです。オイラー方程式は二階微分方程式ですが、式(1)は一階微分方程式の形をしています。 $f$ が少し複雑な形をしている場合、無駄な微分の手間を一回省くというのは非常に大きな節約になります。
2. $f=f(y,y')$
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
まず $f$ の全微分を考えます。
df(y,y')={\partial f\over \partial y}dy+{\partial f\over \partial y'}dy'
少し式変形をして、次のような形にまとめます。
{\partial f\over \partial y}y'={df \over dx}-{\partial f\over \partial y'}{\partial y' \over \partial x} \tag{2}
一方、オイラー方程式(0)の両辺に $y'$ を掛けたものを考えます。
y' \Big({\partial f\over \partial y}-\frac{d}{dx}{\partial f\over \partial y'} \Big)= y'{\partial f\over \partial y}-y'\frac{d}{dx}{\partial f\over \partial y'}=0 \tag{3}
式 (2)(3)を連立して ${\partial f\over \partial y}y'$ を消去します。
f-y'\displaystyle {\partial f\over \partial y'}
=const. \tag{4}
この公式は1868年にベルトラミによって導かれたので、特にベルトラミの公式という名前がついています。この公式が使える問題の中で有名なものに、 最速降下曲線_ を求めるというものがあります。
3. $f=f(x,y)$
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
関数 $f$ に $y'$ が含まれていない場合です。 ${\partial f\over \partial y'}$ の項を落とすと、オイラー方程式(0)は次のような形になります。
{\partial f\over \partial y}
=0 \tag{5}
これは直ちに積分できて、 $f=const.$ が言えるでしょう。変分問題としては、全然面白くない形の関数でした。
4. $f=f(x,y,y',z,z')$
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
関数 $f$ が、 $y$ 以外にも変数を含む場合です。汎関数は次式(6)のように表されます。
\displaystyle I[y_{1},y_{2}]\equiv \int _{a}^{b}f(x,y,y',z,z')dx \tag{6}
ここで、 $y$ も $z$ も $x$ の関数であることに注意しましょう。関数形の微小な変形を、任意関数 $\eta _{1}(x)$ , $\eta _{2}(x)$ と任意の微小量 $\varepsilon$ を使って、次のように表現することにします。
y(x)=y_{0}(x)+\varepsilon \cdot \eta _{1}(x) \tag{7}
z(x)=z_{0}(x)+\varepsilon \cdot \eta _{2}(x) \tag{8}
式(7)には $z$ が含まれていませんし、式(8)には $y$ が含まれていません。そこで、式(7)(8)をそれぞれ別の変分問題だと考えることができます。これは 変分法1_ で求めた形のオイラー方程式に帰着する形ですから、結果が次の式(9)(10)のように書けます。オイラー方程式が二本出てきたということです。これらを連立することで変分問題が解けます。
\displaystyle {\partial f\over \partial y}
-{d\over dx}
\Big({\partial f\over \partial y'}
\Big)=0 \tag{9}
\displaystyle {\partial f\over \partial z}
-{d\over dx}
\Big({\partial f\over \partial z'}
\Big)=0 \tag{10}
変数がたくさんある場合には、変数の数だけオイラー方程式が出てきます。一般に、関数 $f$ が $n$ 個の変数を含み、 $f=f(x,y_{1},y_{1}',...,y_{n},y_{n}')$ と表わされる場合、この $f$ に関する変分問題は次のように、 $n$ 次連立オイラー方程式を解くという問題に帰着します。
\displaystyle {\partial f\over \partial y_{i}}
-{d\over dx}
\Big({\partial f\over \partial y_{i}'}
\Big)=0 \ \ \ \ (i=1,...,n) \tag{11}
5. $f=f(x,y,y',y'')$
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
関数 $f$ が $y$ の二次導関数までを含む場合です。基本的には 変分法1_ で行った計算と同じですので、そのときの手順も思い出しながら計算を追ってみて下さい。
汎関数の形は式(12)のように書けることでしょう。また, $y$ の関数形の微小な変化は、任意の関数 $\eta(x)$ と任意の微小量 $\varepsilon$ を用いて式(13)のように書けます。
\displaystyle I[y]\equiv \int _{a}^{b}f(x,y,y',y'')dx \tag{12}
y(x)=y_{0}(x)+\varepsilon \cdot \eta (x) \tag{13}
式(12)で表される汎関数 $I$ に対し、 ${dI\over d\varepsilon }$ が次式のようになるというは大丈夫でしょうか。一行目では $I$ の全微分を考えています。一行目に出てきた ${\partial y\over \partial \varepsilon }$ などを式(13)を用いて変形したのが二行目です。
\displaystyle {dI\over d\varepsilon }
&=\int _{a}^{b}\Big({\partial f\over \partial y}
\cdot {\partial y\over \partial \varepsilon }
+{\partial f\over \partial y'}
\cdot {\partial y'\over \partial \varepsilon }
+{\partial f\over \partial y''}
\cdot {\partial y''\over \partial \varepsilon }
\Big)dx \\
&=\int _{a}^{b}\Big({\partial f\over \partial y}
\cdot \eta (x)+{\partial f\over \partial y'}
\cdot \eta '(x)+{\partial f\over \partial y''}
\cdot \eta ''(x)\Big)dx \tag{14}
次に、式(14)の二行目の第三項を 部分積分_ します。
{dI\over d\varepsilon }
&=\int _{a}^{b} \Big({\partial f\over \partial y}
\eta (x)+{\partial f\over \partial y'}
\eta '(x)+{\partial f\over \partial y''}
\eta ''(x)\Big)dx \\
&=\bigg[{\partial f\over \partial y''}
\eta '(x)\bigg]_{a}^{b}+\int _{a}^{b}\bigg({\partial f\over \partial y}
\eta (x)
+{\partial f\over \partial y'}
\eta '(x)-{d\over dx}
\Big({\partial f\over \partial y'}
\Big)\eta '(x)
\bigg) dx \\
&=0
積分の中から $\eta ''(x)$ が外に出せました。もう一回頑張って、積分の中の第二項・第三項を部分積分しましょう。今度は $\eta '(x)$ をなくすのが目標です。
&\bigg[{\partial f\over \partial y''}
\eta '(x)\bigg]_{a}^{b}+\int _{a}^{b}\bigg({\partial f\over \partial y}
\eta (x)+{\partial f\over \partial y'}
\eta '(x)-{d\over dx}
\Big( {\partial f\over \partial y'}
\Big) \eta '(x)\bigg) dx \\
&=\bigg[ {\partial f\over \partial y''}
\eta '(x)+{\partial f\over \partial y'}
\eta (x)-{d\over dx}
\Big({\partial f\over \partial y'}
\Big)\eta (x)\bigg]_{a}^{b}+\int _{a}^{b}\eta (x)\biggl\{ {\partial f\over \partial y}
-{d\over dx}
\big( {\partial f\over \partial y'}
\big) +{d^{2}\over dx^{2}}
\big( {\partial f\over \partial y'}
\Big) \biggl\} dx \\
&=0
境界条件より、最初の括弧は全て零になります。 $\eta(x)$ は任意の関数だということでしたから、結局、上式が成り立つためには、次式(15)が成り立つことが必要十分条件となります。これが求める変形オイラー方程式です。オイラー方程式(0)と比べてみてください。
\displaystyle {\partial f\over \partial y}
-{d\over dx}
\bigg({\partial f\over \partial y'}
\bigg)+{d^{2}\over dx^{2}}
\bigg({\partial f\over \partial y'}
\bigg)=0 \tag{15}
6. $f=f(x,y,y',...,y^{(n)})$
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
5.の場合の拡張として、 $f$ が $y$ の $n$ 次導関数まで含む場合を考えましょう。5.で行った計算と同じようにして、 $n$ 回部分積分を繰り返しと、次のような変形オイラー方程式が得られます。
\displaystyle &{\partial f\over \partial y}
-{d\over dx}
\bigg({\partial f\over \partial y'}
\bigg)+{d^{2}\over dx^{2}}
\bigg({\partial f\over \partial y''}
\bigg)-
{d^{3}\over dx^{3}}
\bigg({\partial f\over \partial y^{(3)}}
\bigg)+....... \\
&=\sum \limits _{n=0}^{k}(-1)^{n}{d^{n}\over dx^{n}}
\bigg({\partial f\over \partial y^{(n)}}
\bigg)=0 \tag{16}
7. $f=f(x_{1},x_{2},u(x_{1},x_{2}),u_{x_{1}}(x_{1},x_{2},x_{3})),u_{x_{2}}(x_{1},x_{2},x_{3}))$
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
今まで見てきた関数では、 $y$ や $z$ は $x$ の関数でしたから、結局、一番基本的な変数は $x$ だけでした。ここでは、そのような基本変数が二個ある場合を考えます。このとき汎関数は次式のように書けます。重積分になっているところがポイントです。
\displaystyle I\equiv \int _{a_{2}}^{b_{2}} \int _{a_{1}}^{b_{1}} f(x_{1},x_{2},u,u_{x_{1}},u_{x_{2}})dx_{1}dx_{2}
式中、 $u_{x_{1}}$ というのは、 $u$ の $x_{1}$ による偏微分という意味だとします。例によって、 $u$ の微小な変化を、任意の微小量 $\varepsilon$ と、任意の関数 $\eta (x_{1}, x_{2})$ によって、次のように表しましょう。
u(x_{1},x_{2})=u_{0}(x_{1},x_{2})+\varepsilon \cdot \eta (x_{1},x_{2})\
次に、 $\frac{dI}{d\varepsilon }$ を考えます。 $\eta$ の偏導関数が二種類出てくるのがポイントです。
\displaystyle {\partial I\over \partial \varepsilon}
=\int _{a_{2}}^{b_{2}} \limits \int _{a_{1}}^{b_{1}} \limits \bigg({\partial f\over \partial u}
\eta +{\partial f\over \partial u_{x_{1}}}
\eta _{x_{1}}+{\partial f\over \partial u_{x_{2}}}
\eta _{x_{2}}\bigg)dx_{1} dx_{2}=0
右辺の括弧内の、第二項、第三項をそれぞれ $x_{1}$ と $x_{2}$ で部分積分します。
\displaystyle
\Big[
\frac{\partial f}{\partial u_{x_{1}}}\frac{\partial^2 \eta}{\partial {x_{1}}^2}
\Big]_{a_{1}}^{b_{1}}
+\Big[
\frac{\partial f}{\partial u_{x_{2}}}\frac{\partial^2 \eta}{\partial {x_{2}}^2}
\Big]_{a_{2}}^{b_{2}}
+
\int _{a_{2}}^{b_{2}} \limits \int _{a_{1}}^{b_{1}} \limits \biggl\{ {\partial f\over \partial u}
-{d\over dx_{1}}
\bigg({\partial f\over \partial u_{x_{1}}}
\bigg)-{d\over dx_{2}}
\bigg({\partial f\over \partial u_{x_{2}}}
\bigg)\biggr\} \eta dx_{1}dx_{2}=0
部分積分した前へ出した部分は、境界条件によって零になります。また $\eta$ は任意の関数だということでしたので、結局、上式が成り立つための必要十分条件は、次式(17)が成り立つことです。
\displaystyle {\partial f\over \partial u}
-{d\over dx_{1}}
\bigg({\partial f\over \partial u_{x_{1}}}
\bigg)-{d\over dx_{2}}
\bigg({\partial f\over \partial u_{x_{2}}}
\bigg)=0 \tag{17}
オイラー方程式(0)と見比べてみて下さい。
例題
----------------------------------------------------------------------
1.[直線距離]
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
二点 $P(0,0)$ , $Q(a,b)$ を結ぶ最短経路を求めて下さい。
解答
この答えが直線であることは直感的に明らかですが、変分法を使って解くのが練習です。二点を結ぶ任意の曲線の長さ $l$ は、次のように書けますね。
l=\int _{P}^{Q} ds = \int _{P}^{Q} \sqrt{dx^2 + dy^2} = \int _{0}^{a}\sqrt {1+\Big({dy\over dx}\Big)^{2}}dx= \int _{0}^{a}\sqrt {1+y'^{2}}dx
求めたいのは $l$ の最小値ですから、 $l$ を汎関数にとり、さっそく変分法で解きましょう。
\displaystyle I[y]\equiv \int _{0}^{a}f(x,y,y')dx=\int _{a}^{b} \sqrt{1+y'^2}dx
この問題では $f=\sqrt {1+y'^{2}}$ となっていますね。これは $y$ を含まない形ですから、変形オイラー方程式(2)が使える形です。求める変形オイラー方程式は次のようになります。
\displaystyle {y'\over \sqrt {1+y'^{2}}}
=const.
これを一回積分すれば、 $y'=C_{1}$ となります。もう一度積分して次の答えが得られます。
y=C_{1}x+C_{2}
積分定数 $C_{1}$ , $C_{2}$ は、この直線が二点 $P$ , $Q$ を通るという条件から決まります。予想通り、 $P$ , $Q$ を通る直線が、最短経路だということになりました。
2.[最速降下曲線]
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
高さの異なる二点を坂道で結び、玉を転がすとき、玉が一番速く転がる坂の形(曲線)を考えてください。重力は一様だとします。
この問題の解答は、 最速降下曲線_ をご覧下さい。
(ヒント:式(4)が使える問題です。)
3.[懸垂曲線]
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
密度が一様な紐を吊り下げたとき、ポテンシャルエネルギーが最小になるような紐の垂れ方(曲線)を考えて下さい。重力は一様だとします。
解答
紐の微小部分が持つ位置エネルギーを全て合わせたものを、紐全体のポテンシャルエネルギー $U$ だとします。紐の微小部分 $ds$ が持つ位置エネルギーは、その微小部分の高さを $y$ だとして $mgyds=mgy\sqrt{1+y'}dx$ と表わされますので、 $U$ は次のように、積分で表すことができます。
\displaystyle U=mg\int _{P}^{Q}y\sqrt {1+y'^{2}}dx
ここで $U$ の最小値を求めたいのですから、 $I=U$ として、変分問題を解きましょう。
\displaystyle I=\int _{P}^{Q} fdx=mg\int _{P}^{Q}y\sqrt {1+y'^{2}}dx
この問題では $f=y\sqrt {1+y'^{2}}$ となっていて、 $f$ が $x$ を含まない形ですから、変形オイラー方程式(4)が使える形です。解くべき式は次のようになります。
f-y'\displaystyle {\partial f\over \partial y'}
&=y\sqrt {1+y'^{2}}-y'{yy'\over \sqrt {1+y'^{2}}} \\
&={y\over \sqrt {1+y'^{2}}} \\
&=const.
これを積分して、次の答えを得ます。
y=C_{1}\cosh C_{2}x
これは双曲線関数(ハイパーボリックコサイン)と呼ばれる曲線です。 [*]_
.. [*] 余談ですが、西洋の大きな建築物は、たいてい石で出来ています。石を削るのはなかなか大変なことですから、自然と、予め地上で石を設計図通りに切っておいて、上で積み上げるという方式が取られることになります。設計図として簡単に書けるのは、定規とコンパスだけで作図できる図形ですから、古い西洋建築に見られる曲線というのは、大方、円弧なのです。ヨーロッパに旅行に行く機会があれば、どこの街にも古い教会などがありますから、是非窓枠や屋根の曲線を眺めて見てください。ほとんど全て円弧の一部のはずです。さて、たった今求めた双曲線関数ですが、これは $\cosh x = \frac{e^{x}+e^{-x}}{2}$ という関数ですので、ちょっと定規とコンパスだけで作図するのは無理そうです。西洋で双曲線関数を始めて建築に取り入れたのは、スペインのアントニ・ガウディで、19世紀末の話です。ところが、わが日本では、神社の屋根の微妙な反り返りなどがほとんど全て双曲線関数になっているのです!日本では設計図を描いてから木を削るのではなく、大工さんが鎖を垂らして曲線を決めていくので、自然と双曲線関数になってしまうのでした。こんど神社に行くようなことがあれば、是非屋根をながめて見て下さい。西洋で、ようやくここ100年そこそこになって使われ始めた曲線が、日本では何百年も伝統的に使われているとは、なんだか嬉しい話です。
.. image:: Joh-Shrine.jpg
.. _ここ: http://www.york.ac.uk/depts/maths/histstat/people/bernoulli_tree.htm
.. _変分法1:
.. _部分積分:
.. _最速降下曲線:
@@author:Joh@@
@@accept: 2005-03-31@@