跳到主要内容

名古屋大学 情報学研究科 複雑系科学専攻 2017年8月実施 数2 [2]

Author​

思齐塾, 祭音Myyura

Description​

出典:名古屋大学公式問題

非線形常微分方程式

d2ydt2−c(dydt)2+g=0(t≥0)\frac{d^2y}{dt^2} - c\left(\frac{dy}{dt}\right)^2 + g = 0 \quad (t \geq 0)

を考える。ただし、 cc および gg は正定数である。

  1. v(t)=dydtv(t) = \frac{dy}{dt} として、 t=0t=0 のときの条件 v(0)=0v(0) = 0 のもとで v(t)v(t) の微分方程式を解け。 t→∞t \to \infty での v(∞)v(\infty) の値を求めよ。

  2. t=0t=0 のときの条件 y(0)=h(>0)y(0) = h(>0) のもとで y(t)y(t) の微分方程式を解け。

  3. 解 y(t)y(t) を t−yt-y 平面上に図示せよ。 t→∞t \to \infty で解が漸近する直線とその傾き, 切片の値を明記せよ。

题目描述​

考察非线性常微分方程

d2ydt2−c(dydt)2+g=0,t≥0,\frac{d^2y}{dt^2} -c\left(\frac{dy}{dt}\right)^2+g=0, \qquad t\ge0,

其中 c,gc,g 均为正常数。

  1. 令 v(t)=dydtv(t)=\dfrac{dy}{dt},在初始条件 v(0)=0v(0)=0 下求解 v(t)v(t) 所满足的微分方程,并求

    v(∞)=lim⁡t→∞v(t);v(\infty)=\lim_{t\to\infty}v(t);
  2. 在初始条件 y(0)=hy(0)=h(h>0h>0)下求解 y(t)y(t);

  3. 在 tt–yy 平面上画出解 y(t)y(t),并明确写出 t→∞t\to\infty 时解所趋近的渐近直线,以及该直线的斜率和截距。

Kai​

  1. Let v(t)=dydtv(t) = \frac{dy}{dt} . Then d2ydt2=dvdt\frac{d^2y}{dt^2} = \frac{dv}{dt} . The equation becomes
dvdt−cv2+g=0\frac{dv}{dt} - cv^2 + g = 0
dvdt=cv2−g\frac{dv}{dt} = cv^2 - g

初期値 v(0)=0v(0)=0 に合う枝は

v(t)=−gc tanh⁡(cg t)\boxed{v(t)=-\sqrt{\frac gc}\, \tanh(\sqrt{cg}\,t)}

である。実際、

v′=−g sech⁡2(cg t)=cv2−gv'=-g\,\operatorname{sech}^2(\sqrt{cg}\,t) =cv^2-g

かつ v(0)=0v(0)=0 である。したがって

lim⁡t→∞v(t)=−gc.\boxed{\lim_{t\to\infty}v(t)=-\sqrt{\frac gc}}.
  1. v=y′v=y' を積分すると
y(t)=h−gc∫0ttanh⁡(cg s) ds=h−1clog⁡cosh⁡(cg t).y(t) =h-\sqrt{\frac gc}\int_0^t \tanh(\sqrt{cg}\,s)\,ds =\boxed{h-\frac1c\log\cosh(\sqrt{cg}\,t)}.

この式は y(0)=hy(0)=h を満たし、二回微分すれば元の方程式へ直接戻る。

  1. t≥0t\ge0 で y′(t)≤0y'(t)\le0 、 y′′(t)=−g sech⁡2(cg t)<0y''(t)=-g\,\operatorname{sech}^2(\sqrt{cg}\,t)<0 なので、グラフは (0,h)(0,h) から水平な接線で出発し、上に凸のまま減少する。また
log⁡cosh⁡(cg t)=cg t−log⁡2+o(1)\log\cosh(\sqrt{cg}\,t) =\sqrt{cg}\,t-\log2+o(1)

より、漸近直線は

y=−gc t+h+log⁡2c.\boxed{ y=-\sqrt{\frac gc}\,t+h+\frac{\log2}{c} }.

その傾きは −g/c-\sqrt{g/c} 、切片は h+(log⁡2)/ch+(\log2)/c である。

下図は τ=cg t\tau=\sqrt{cg}\,t、Y=c(y−h)Y=c(y-h) による表示である。横軸を 1/cg1/\sqrt{cg} 倍し、縦軸を 1/c1/c 倍して hh だけ上へ移すと、tt–yy 平面のグラフを得る。

解 Y=-log cosh τ と漸近線 Y=-τ+log 2