見出し画像

Eulerの微分方程式

$${I[y]}$$は関数$${F=F(x,y(x),y'(x))}$$の積分を含む汎関数とする。
 $${I[y]=\displaystyle\int_a^bF(x,y,y')dx}$$ (1)
$${x=a,  x=b}$$の両端の$${y(x)}$$を固定した区間内で考える。

ここで $${y(x)}$$がある関数形の時に汎関数$${I[y]}$$が停滞値(極大,極小、変曲)をとるものとし、$${y(x)}$$の関数形を少し変えても$${I[y(x)]}$$の値が変わらない関数形$${y(x)+\delta y(x)}$$を考える。
  黒線:$${I[y]}$$が最小値をとる関数形$${y(x)}$$
  赤線:$${I[y]}$$の値が変わらない関数形$${y(x)+\delta y(x)}$$
  また $${\delta y(a)=\delta y(b)=0}$$ (2)
このとき$${I[y(x)+\delta y(x)]}$$と元の$${I[y(x)]}$$との差は$${\delta I[y(x)]}$$は0である。
      (以下$${(x)}$$を省略)
 $${\delta I[y]=I[y+\delta y]-I[y]=0}$$ (3)
 $${\delta I[y]=I[y+\delta y]-I[y]}$$ ← (1)を代入
    $${=\displaystyle\int_a^bF(x,  y+\delta y,  y'+\delta y')dx-\displaystyle\int_a^bF(x,  y,  y')dx}$$
    $${=\displaystyle\int_a^b\Big\{F(x,  y+\delta y,  y'+\delta y')-F(x,  y,  y')\Big\}dx}$$
    $${=\displaystyle\int_a^b\delta F(x,y,y')dx}$$
    $${=\displaystyle\int_a^b\Big({\dfrac{\partial F}{\partial x}\delta x}+\dfrac{\partial F}{\partial y}\delta y+\dfrac{\partial F}{\partial y'}\delta y'\Big)dx}$$
    $${=\displaystyle\int_a^b\Big(\dfrac{\partial F}{\partial y}\delta y+\dfrac{\partial F}{\partial y'}\delta y'\Big)dx}$$ ($${\dfrac{\partial F}{\partial x}\delta x=0}$$、$${x}$$は変化していない )
    $${=\displaystyle\int_a^b\dfrac{\partial F}{\partial y}\delta y  dx+\displaystyle\int_a^b\dfrac{\partial F}{\partial y'}\delta y'dx}$$
    $${=\displaystyle\int_a^b\dfrac{\partial F}{\partial y}\delta y  dx+\displaystyle\int_a^b\dfrac{\partial F}{\partial y'}\delta\Big(\dfrac{dy}{dx}\Big)dx}$$
    $${=\displaystyle\int_a^b\dfrac{\partial F}{\partial y}\delta y  dx+\displaystyle\int_a^b\dfrac{\partial F}{\partial y'}\dfrac{d}{dx}\Big(\delta y\Big)dx}$$ (第2項を部分積分)
    $${=\displaystyle\int_a^b\dfrac{\partial F}{\partial y}\delta y  dx+\Bigg[\dfrac{\partial F}{\partial y'}\delta y\Bigg]_a^b-\displaystyle\int_a^b\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)\delta y  dx}$$
            ((2)$${\delta y(a)=\delta y(a)=0}$$より  $${\Bigg[\dfrac{\partial F}{\partial y'}\delta y\Bigg]_a^b=0}$$)
    $${=\displaystyle\int_a^b\dfrac{\partial F}{\partial y}\delta y  dx-\displaystyle\int_a^b\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)\delta y  dx}$$
    $${=\displaystyle\int_a^b\Big\{\dfrac{\partial F}{\partial y}-\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)\Big\}\delta y  dx}$$
(3)と併せて
 $${\delta I=\displaystyle\int_a^b\Big\{\dfrac{\partial F}{\partial y}-\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)\Big\}\delta y  dx=0}$$
この式が任意の$${\delta y}$$で成り立つためには
 $${\dfrac{\partial F}{\partial y}-\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)=0}$$
よって
 $${\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)-\dfrac{\partial F}{\partial y}=0}$$ (Eulerの微分方程式)(4)
この式を満たす関数$${y(x)}$$が汎関数$${I[y]=\displaystyle\int_a^bF(x,y,y')dx}$$が停滞値(極大、極小、変曲)を持つための条件である。
Lagrangeの運動方程式と同じ形だが、それについては別に考えることにする。(→ https://note.com/greedyoldman/n/nad19e5d1dbb4
 $${\dfrac{d}{dx}\Big(\dfrac{\partial \mathscr{L}}{\partial y'}\Big)-\dfrac{\partial \mathscr{L}}{\partial y}=0}$$ (Lagrangeの運動方程式)

追記
$${F=F(y,  y')}$$で$${x}$$を含まないとき
$${F(y,  y')}$$の全微分は
 $${dF=\dfrac{\partial F}{\partial y}dy+\dfrac{\partial F}{\partial y'}dy'}$$
両辺÷$${dx}$$
 $${\dfrac{dF}{dx}=\dfrac{\partial F}{\partial y}\dfrac{dy}{dx}+\dfrac{\partial F}{\partial y'}\dfrac{dy'}{dx}=\dfrac{\partial F}{\partial y}y'+\dfrac{\partial F}{\partial y'}\dfrac{dy'}{dx}}$$
よって
 $${\dfrac{\partial F}{\partial y}=\dfrac{1}{y'}\Big(\dfrac{dF}{dx}-\dfrac{\partial F}{\partial y'}\dfrac{dy'}{dx}\Big)}$$
これを Euler の微分方程式に代入
 $${\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)-\underline{\dfrac{\partial F}{\partial y}}=0}$$ (4)
 $${\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)-\dfrac{1}{y'}\Big(\dfrac{dF}{dx}-\dfrac{\partial F}{\partial y'}\dfrac{dy'}{dx}\Big)=0}$$
$${× y'}$$
 $${\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}\Big)y'-\dfrac{dF}{dx}+\dfrac{\partial F}{\partial y'}\dfrac{dy'}{dx}=0}$$
 $${\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}  y'\Big)-\dfrac{dF}{dx}=0}$$
 $${\dfrac{d}{dx}\Big(\dfrac{\partial F}{\partial y'}  y'-F\Big)=0}$$
積分して整理すると
 $${F-y'  \dfrac{\partial F}{\partial y'}=const}$$ (Beltrami の恒等式)
が得られる。

この記事が気に入ったらサポートをしてみませんか?