跳到主要内容

東京大学 新領域創成科学研究科 人間環境学専攻 2019年8月実施 微分方程式

Author​

思齐塾, 祭音Myyura

Description​

下記の微分方程式の解 y(x)y(x) を求めよ.

  1. y′=xy+x+2y+2y' = xy + x + 2y + 2

  2. y′′+y=sin⁡xy'' + y = \sin x y(0)=0,y′(0)=0y(0) = 0, y'(0) = 0

  3. y′′′+4y′′+3y′=6xy''' + 4y'' + 3y' = 6x y(0)=3,y′(0)=−203,y′′(0)=12y(0) = 3, y'(0) = -\frac{20}{3}, y''(0) = 12

ただし、

L[f(n)(x)]=snF(s)−sn−1f(0)−⋯−sf(n−2)(0)−f(n−1)(0)\mathcal{L}[f^{(n)}(x)] = s^n F(s) - s^{n-1}f(0) - \cdots - sf^{(n-2)}(0) - f^{(n-1)}(0)

また、

f∗g=∫0xf(x−u)g(u)du,x>0f * g = \int_0^x f(x-u)g(u) du, x > 0

について、

L[f∗g]=L[f]L[g]\mathcal{L}[f * g] = \mathcal{L}[f] \mathcal{L}[g]

が成立することを用いてもよい。

题目描述​

分别求下列微分方程或初值问题的解 y(x)y(x):

  1. 求解:

    y′=xy+x+2y+2.y'=xy+x+2y+2.
  2. 求解初值问题:

    y′′+y=sin⁡x,y(0)=0,y′(0)=0.y''+y=\sin x,\qquad y(0)=0,\quad y'(0)=0.
  3. 求解:

    y′′′+4y′′+3y′=6x,y'''+4y''+3y'=6x,
    y(0)=3,y′(0)=−203,y′′(0)=12.y(0)=3,\qquad y'(0)=-\frac{20}{3},\qquad y''(0)=12.

解题时可以使用导数的拉普拉斯变换公式

L[f(n)(x)]=snF(s)−sn−1f(0)−⋯−sf(n−2)(0)−f(n−1)(0),\mathcal L[f^{(n)}(x)] =s^nF(s)-s^{n-1}f(0)-\cdots-sf^{(n-2)}(0)-f^{(n-1)}(0),

以及卷积

(f∗g)(x)=∫0xf(x−u)g(u) du(x>0)(f*g)(x)=\int_0^x f(x-u)g(u)\,du\qquad(x>0)

所满足的卷积定理

L[f∗g]=L[f] L[g].\mathcal L[f*g]=\mathcal L[f]\,\mathcal L[g].

Kai​

  1. y′=xy+x+2y+2y' = xy + x + 2y + 2 可以写成 dydx=(x+2)(y+1)\frac{dy}{dx} = (x+2)(y+1) .
dyy+1=(x+2)dx\frac{dy}{y+1} = (x+2)dx

两边积分得到:

∫dyy+1=∫(x+2)dx\int \frac{dy}{y+1} = \int (x+2) dx
ln⁡∣y+1∣=x22+2x+C\ln|y+1| = \frac{x^2}{2} + 2x + C

符号を任意定数に含めると

y+1=Aex22+2x(A∈R)y+1 = Ae^{\frac{x^2}{2} + 2x}\qquad(A\in\mathbb{R})
y=Aex22+2x−1y = Ae^{\frac{x^2}{2} + 2x} - 1

A=0A=0 は定数解 y=−1y=-1 を含む。

  1. y′′+y=sin⁡x,y(0)=0,y′(0)=0y'' + y = \sin x, y(0) = 0, y'(0) = 0 特征方程 r2+1=0r^2 + 1 = 0 , r=±ir = \pm i , 所以 yh=c1cos⁡x+c2sin⁡xy_h = c_1 \cos x + c_2 \sin x . 令 yp=Axcos⁡x+Bxsin⁡xy_p = Ax \cos x + Bx \sin x , 则 yp′=Acos⁡x−Axsin⁡x+Bsin⁡x+Bxcos⁡xy_p' = A \cos x - Ax \sin x + B \sin x + Bx \cos x , yp′′=−Asin⁡x−Asin⁡x−Axcos⁡x+Bcos⁡x+Bcos⁡x−Bxsin⁡x=−2Asin⁡x−Axcos⁡x+2Bcos⁡x−Bxsin⁡xy_p'' = -A \sin x - A \sin x - Ax \cos x + B \cos x + B \cos x - Bx \sin x = -2A \sin x - Ax \cos x + 2B \cos x - Bx \sin x . 代入原方程,得到: (−2Asin⁡x−Axcos⁡x+2Bcos⁡x−Bxsin⁡x)+(Axcos⁡x+Bxsin⁡x)=sin⁡x(-2A \sin x - Ax \cos x + 2B \cos x - Bx \sin x) + (Ax \cos x + Bx \sin x) = \sin x −2Asin⁡x+2Bcos⁡x=sin⁡x-2A \sin x + 2B \cos x = \sin x , 所以 A=−12,B=0A = -\frac{1}{2}, B = 0 . 因此 yp=−12xcos⁡xy_p = -\frac{1}{2} x \cos x . 因此 y=c1cos⁡x+c2sin⁡x−12xcos⁡xy = c_1 \cos x + c_2 \sin x - \frac{1}{2} x \cos x . y(0)=c1=0y(0) = c_1 = 0 . y′=−c1sin⁡x+c2cos⁡x−12cos⁡x+12xsin⁡xy' = -c_1 \sin x + c_2 \cos x - \frac{1}{2} \cos x + \frac{1}{2} x \sin x . y′(0)=c2−12=0y'(0) = c_2 - \frac{1}{2} = 0 , 所以 c2=12c_2 = \frac{1}{2} . 因此 y=12sin⁡x−12xcos⁡xy = \frac{1}{2} \sin x - \frac{1}{2} x \cos x .

  2. v=y′v=y' とおくと

v′′+4v′+3v=6x,v(0)=−203,v′(0)=12.v''+4v'+3v=6x,\qquad v(0)=-\frac{20}{3},\quad v'(0)=12.

その一般解を

v=C1e−x+C2e−3x+2x−83v=C_1e^{-x}+C_2e^{-3x}+2x-\frac83

と書ける。初期条件から

C1+C2=−4,C1+3C2=−10,C_1+C_2=-4,\qquad C_1+3C_2=-10,

したがって C1=−1, C2=−3C_1=-1,\ C_2=-3 である。もう一度積分すると

y=e−x+e−3x+x2−83x+C3.y=e^{-x}+e^{-3x}+x^2-\frac83x+C_3.

y(0)=3y(0)=3 より C3=1C_3=1 なので、

y(x)=1+e−x+e−3x+x2−83x.\boxed{y(x)=1+e^{-x}+e^{-3x}+x^2-\frac83x}.