跳到主要内容

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

Author

思齐塾, 祭音Myyura

Description

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

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

  2. y+y=sinxy'' + 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)sn1f(0)sf(n2)(0)f(n1)(0)\mathcal{L}[f^{(n)}(x)] = s^n F(s) - s^{n-1}f(0) - \cdots - sf^{(n-2)}(0) - f^{(n-1)}(0)

また、

fg=0xf(xu)g(u)du,x>0f * g = \int_0^x f(x-u)g(u) du, x > 0

について、

L[fg]=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=sinx,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)sn1f(0)sf(n2)(0)f(n1)(0),\mathcal L[f^{(n)}(x)] =s^nF(s)-s^{n-1}f(0)-\cdots-sf^{(n-2)}(0)-f^{(n-1)}(0),

以及卷积

(fg)(x)=0xf(xu)g(u)du(x>0)(f*g)(x)=\int_0^x f(x-u)g(u)\,du\qquad(x>0)

所满足的卷积定理

L[fg]=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
lny+1=x22+2x+C\ln|y+1| = \frac{x^2}{2} + 2x + C

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

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

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

  1. y+y=sinx,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=c1cosx+c2sinxy_h = c_1 \cos x + c_2 \sin x . 令 yp=Axcosx+Bxsinxy_p = Ax \cos x + Bx \sin x , 则 yp=AcosxAxsinx+Bsinx+Bxcosxy_p' = A \cos x - Ax \sin x + B \sin x + Bx \cos x , yp=AsinxAsinxAxcosx+Bcosx+BcosxBxsinx=2AsinxAxcosx+2BcosxBxsinxy_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 . 代入原方程,得到: (2AsinxAxcosx+2BcosxBxsinx)+(Axcosx+Bxsinx)=sinx(-2A \sin x - Ax \cos x + 2B \cos x - Bx \sin x) + (Ax \cos x + Bx \sin x) = \sin x 2Asinx+2Bcosx=sinx-2A \sin x + 2B \cos x = \sin x , 所以 A=12,B=0A = -\frac{1}{2}, B = 0 . 因此 yp=12xcosxy_p = -\frac{1}{2} x \cos x . 因此 y=c1cosx+c2sinx12xcosxy = c_1 \cos x + c_2 \sin x - \frac{1}{2} x \cos x . y(0)=c1=0y(0) = c_1 = 0 . y=c1sinx+c2cosx12cosx+12xsinxy' = -c_1 \sin x + c_2 \cos x - \frac{1}{2} \cos x + \frac{1}{2} x \sin x . y(0)=c212=0y'(0) = c_2 - \frac{1}{2} = 0 , 所以 c2=12c_2 = \frac{1}{2} . 因此 y=12sinx12xcosxy = \frac{1}{2} \sin x - \frac{1}{2} x \cos x .

  2. v=yv=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=C1ex+C2e3x+2x83v=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=ex+e3x+x283x+C3.y=e^{-x}+e^{-3x}+x^2-\frac83x+C_3.

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

y(x)=1+ex+e3x+x283x.\boxed{y(x)=1+e^{-x}+e^{-3x}+x^2-\frac83x}.