跳到主要内容

東京大学 新領域創成科学研究科 先端エネルギー工学専攻 2017年度 第1問(数学)

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

区間 [x0,x1][x_0,x_1] で定義され、g(x0)=g(x1)=0g(x_0)=g(x_1)=0 である連続な関数 g(x)g(x) がある。この区間における曲線 y=g(x)y=g(x) の長さを LL、

S=∫x0x1g(x) dxS=\int_{x_0}^{x_1}g(x)\,dx

とするとき、定数 λ\lambda を用いて I=S+λLI=S+\lambda L と定義する。以下の問に答えよ。

問1 LL を定積分の形で表せ。

問2

I=∫x0x1f(g(x),g′(x)) dxI=\int_{x_0}^{x_1}f\bigl(g(x),g'(x)\bigr)\,dx

と表した場合における関数 f(g(x),g′(x))f(g(x),g'(x)) を定義する。ただし、X′=dX/dxX'=dX/dx とする。

(a) 微小変化 g(x)+η(x)g(x)+\eta(x) を考える。ただし、η(x0)=η(x1)=0\eta(x_0)=\eta(x_1)=0 とする。

f(g+η,g′+η′)≃f+∂f∂gη+∂f∂g′η′f(g+\eta,g'+\eta') \simeq f+\frac{\partial f}{\partial g}\eta +\frac{\partial f}{\partial g'}\eta'

と近似できるとして、

δI=∫x0x1f(g+η,g′+η′) dx−∫x0x1f(g,g′) dx\delta I=\int_{x_0}^{x_1}f(g+\eta,g'+\eta')\,dx -\int_{x_0}^{x_1}f(g,g')\,dx

を表す式を求めよ。

(b) ddx(∂f∂g′η)\dfrac{d}{dx}\left(\dfrac{\partial f}{\partial g'}\eta\right) を計算し、任意の η(x)\eta(x) に対して δI=0\delta I=0 となるための必要条件を求めよ。

(c) (b) の条件のもとで ddx(f−g′∂f∂g′)\dfrac{d}{dx}\left(f-g'\dfrac{\partial f}{\partial g'}\right) を表す式を求め、f−g′∂f∂g′f-g'\dfrac{\partial f}{\partial g'} が任意の xx に対して満たすべき式を示せ。

問3 問2(c) で求めた式が満たされるとき、g(x)g(x) を求め、区間 [x0,x1][x_0,x_1] における曲線 y=g(x)y=g(x) の概形を描け。

题目描述​

连续函数 g(x)g(x) 定义在区间 [x0,x1][x_0,x_1] 上,且 g(x0)=g(x1)=0g(x_0)=g(x_1)=0。记曲线 y=g(x)y=g(x) 在该区间内的长度为 LL,并令

S=∫x0x1g(x) dx,I=S+λL,S=\int_{x_0}^{x_1}g(x)\,dx,\qquad I=S+\lambda L,

其中 λ\lambda 为常数。

问1 将 LL 表示为定积分。

问2 将 II 写为 ∫x0x1f(g,g′) dx\int_{x_0}^{x_1}f(g,g')\,dx,其中撇号表示对 xx 求导。

(a) 考虑微小变化 g+ηg+\eta,其中 η(x0)=η(x1)=0\eta(x_0)=\eta(x_1)=0。利用一阶近似 f(g+η,g′+η′)≃f+fgη+fg′η′f(g+\eta,g'+\eta')\simeq f+f_g\eta+f_{g'}\eta',求变化量 δI=∫f(g+η,g′+η′) dx−∫f(g,g′) dx\delta I=\int f(g+\eta,g'+\eta')\,dx-\int f(g,g')\,dx 的表达式,两积分的上下限均为 x0,x1x_0,x_1。

(b) 计算 ddx(fg′η)\dfrac{d}{dx}(f_{g'}\eta),求对任意满足端点条件的 η\eta 都有 δI=0\delta I=0 的必要条件。

(c) 在 (b) 的条件下计算 ddx(f−g′fg′)\dfrac{d}{dx}(f-g'f_{g'}),写出 f−g′fg′f-g'f_{g'} 应满足的关系式。

问3 在问2(c) 所得关系成立时求 g(x)g(x),并画出曲线在 [x0,x1][x_0,x_1] 上的概形。

Kai​

x0<x1x_0<x_1 とし、以下では変分計算に必要な微分可能性を仮定する。

問1​

弧長要素は ds=dx2+dy2=1+g′(x)2 dxds=\sqrt{dx^2+dy^2}=\sqrt{1+g'(x)^2}\,dx であるから、

L=∫x0x11+g′(x)2 dx.\boxed{L=\int_{x_0}^{x_1}\sqrt{1+g'(x)^2}\,dx}.

問2​

被積分関数は

f(g,g′)=g+λ1+g′2f(g,g')=g+\lambda\sqrt{1+g'^2}

である。

(a)​

η\eta の一次の項まで取ると、

δI≃∫x0x1(∂f∂gη+∂f∂g′η′) dx.\boxed{\delta I\simeq \int_{x_0}^{x_1}\left( \frac{\partial f}{\partial g}\eta +\frac{\partial f}{\partial g'}\eta' \right)\,dx}.

以下、δI\delta I はこの一次変分を表すものとする。

(b)​

積の微分より、

ddx(∂f∂g′η)=ddx(∂f∂g′)η+∂f∂g′η′.\frac{d}{dx}\left(\frac{\partial f}{\partial g'}\eta\right) =\frac{d}{dx}\left(\frac{\partial f}{\partial g'}\right)\eta +\frac{\partial f}{\partial g'}\eta'.

したがって、部分積分と η(x0)=η(x1)=0\eta(x_0)=\eta(x_1)=0 により、

δI=[∂f∂g′η]x0x1+∫x0x1{∂f∂g−ddx(∂f∂g′)}η dx=∫x0x1{∂f∂g−ddx(∂f∂g′)}η dx.\begin{aligned} \delta I &=\left[\frac{\partial f}{\partial g'}\eta\right]_{x_0}^{x_1} +\int_{x_0}^{x_1}\left\{ \frac{\partial f}{\partial g} -\frac{d}{dx}\left(\frac{\partial f}{\partial g'}\right) \right\}\eta\,dx\\ &=\int_{x_0}^{x_1}\left\{ \frac{\partial f}{\partial g} -\frac{d}{dx}\left(\frac{\partial f}{\partial g'}\right) \right\}\eta\,dx. \end{aligned}

η\eta が任意であるから、必要条件はオイラー・ラグランジュ方程式

∂f∂g−ddx(∂f∂g′)=0\boxed{\frac{\partial f}{\partial g} -\frac{d}{dx}\left(\frac{\partial f}{\partial g'}\right)=0}

である。

(c)​

f=f(g,g′)f=f(g,g') は xx に陽に依存しないので、

ddx(f−g′∂f∂g′)=∂f∂gg′+∂f∂g′g′′−g′′∂f∂g′−g′ddx(∂f∂g′)=g′{∂f∂g−ddx(∂f∂g′)}=0.\begin{aligned} \frac{d}{dx}\left(f-g'\frac{\partial f}{\partial g'}\right) &=\frac{\partial f}{\partial g}g' +\frac{\partial f}{\partial g'}g'' -g''\frac{\partial f}{\partial g'} -g'\frac{d}{dx}\left(\frac{\partial f}{\partial g'}\right)\\ &=g'\left\{\frac{\partial f}{\partial g} -\frac{d}{dx}\left(\frac{\partial f}{\partial g'}\right)\right\}=0. \end{aligned}

よって、ある定数 CC に対して

f−g′∂f∂g′=C\boxed{f-g'\frac{\partial f}{\partial g'}=C}

が成り立つ。これはベルトラミの恒等式である。

問3​

問2(c) に f=g+λ1+g′2f=g+\lambda\sqrt{1+g'^2} を代入すると、

g+λ1+g′2=C.(1)g+\frac{\lambda}{\sqrt{1+g'^2}}=C. \tag{1}

ここで停留曲線は問2(b) の条件も満たすので、

1−λddx(g′1+g′2)=0.(2)1-\lambda\frac{d}{dx}\left(\frac{g'}{\sqrt{1+g'^2}}\right)=0. \tag{2}

λ=0\lambda=0 では (2) は成立せず、停留曲線は存在しない。λ≠0\lambda\ne0 のとき、(2) を積分すると、定数 aa を用いて

g′1+g′2=x−aλ.\frac{g'}{\sqrt{1+g'^2}}=\frac{x-a}{\lambda}.

R=∣λ∣R=|\lambda|、σ=sgn⁡λ\sigma=\operatorname{sgn}\lambda とおけば、∣x−a∣<R|x-a|<R において

g′=σx−aR2−(x−a)2,g(x)=C−σR2−(x−a)2.g'=\sigma\frac{x-a}{\sqrt{R^2-(x-a)^2}}, \qquad g(x)=C-\sigma\sqrt{R^2-(x-a)^2}.

したがって、曲線は円

(x−a)2+(g−C)2=R2(x-a)^2+(g-C)^2=R^2

の一部である。端点条件から (x0−a)2=(x1−a)2(x_0-a)^2=(x_1-a)^2 となるため、

m=x0+x12,d=x1−x02m=\frac{x_0+x_1}{2},\qquad d=\frac{x_1-x_0}{2}

とおくと、a=ma=m、C=σR2−d2C=\sigma\sqrt{R^2-d^2} である。よって、

g(x)=sgn⁡(λ){λ2−d2−λ2−(x−m)2}.\boxed{ g(x)=\operatorname{sgn}(\lambda) \left\{\sqrt{\lambda^2-d^2} -\sqrt{\lambda^2-(x-m)^2}\right\} }.

閉区間上で有限な傾きをもつ解には ∣λ∣>d|\lambda|>d が必要である。∣λ∣=d|\lambda|=d は、両端の接線が鉛直になる半円の極限に対応する。∣λ∣<d|\lambda|<d では、両端を結ぶこのような停留曲線は存在しない。

λ>0\lambda>0 では g′′>0g''>0、g(x)≤0g(x)\le0 の下側の円弧であり、λ<0\lambda<0 では g′′<0g''<0、g(x)≥0g(x)\ge0 の上側の円弧である。いずれも直線 x=mx=m に関して対称で、半径は ∣λ∣|\lambda| である。

両端を固定した停留曲線の概形。負のλでは上側、正のλでは下側の円弧となる。

なお、(1) だけなら g≡0g\equiv0 も C=λC=\lambda として満たすが、これは (2) を満たさず、停留曲線ではない。

Reference​