跳到主要内容

名古屋大学 情報学研究科 複雑系科学専攻 2016年8月実施 微分方程

Author

思齐塾, 祭音Myyura

Description

非線形常微分方程式

d2ydt2c(dydt)2+g=0(t0)\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) の微分方程式を解け。 tt \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)tyt-y 平面上に図示せよ。 tt \to \infty で解が漸近する直線とその傾き, 切片の値を明記せよ。

题目描述

考察非线性常微分方程

d2ydt2c(dydt)2+g=0,t0,\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()=limtv(t);v(\infty)=\lim_{t\to\infty}v(t);
  2. 在初始条件 y(0)=hy(0)=hh>0h>0)下求解 y(t)y(t)

  3. ttyy 平面上画出解 y(t)y(t),并明确写出 tt\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
dvdtcv2+g=0\frac{dv}{dt} - cv^2 + g = 0
dvdt=cv2g\frac{dv}{dt} = cv^2 - g

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

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

である。実際、

v=gsech2(cgt)=cv2gv'=-g\,\operatorname{sech}^2(\sqrt{cg}\,t) =cv^2-g

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

limtv(t)=gc.\boxed{\lim_{t\to\infty}v(t)=-\sqrt{\frac gc}}.
  1. v=yv=y' を積分すると
y(t)=hgc0ttanh(cgs)ds=h1clogcosh(cgt).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. t0t\ge0y(t)0y'(t)\le0y(t)=gsech2(cgt)<0y''(t)=-g\,\operatorname{sech}^2(\sqrt{cg}\,t)<0 なので、グラフは (0,h)(0,h) から水平な接線で出発し、上に凸のまま減少する。また
logcosh(cgt)=cgtlog2+o(1)\log\cosh(\sqrt{cg}\,t) =\sqrt{cg}\,t-\log2+o(1)

より、漸近直線は

y=gct+h+log2c.\boxed{ y=-\sqrt{\frac gc}\,t+h+\frac{\log2}{c} }.

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