跳到主要内容

東京工業大学 環境・社会理工学院 融合理工学系 2018年8月実施 微分方程

Author

思齐塾, 祭音Myyura

Description

  1. Solve the following equations using real numbers.

(1) f(t)f(t)6f(t)=0,f(0)=1,f(0)=2f''(t) - f'(t) - 6f(t) = 0, \quad f(0) = 1, \quad f'(0) = 2

(2) f(t)4f(t)+4f(t)=e2t,f(0)=0,f(0)=3f''(t) - 4f'(t) + 4f(t) = e^{2t}, \quad f(0) = 0, \quad f'(0) = 3

(3) f(t)+9f(t)=6cos3t,f(0)=1,f(0)=3f''(t) + 9f(t) = 6 \cos 3t, \quad f(0) = 1, \quad f'(0) = 3

题目描述

用实数形式求下列三个二阶常微分方程初值问题的精确解:

  1. 求解初值问题:

    f(t)f(t)6f(t)=0,f(0)=1,f(0)=2.f''(t)-f'(t)-6f(t)=0,\qquad f(0)=1,\qquad f'(0)=2.
  2. 求解初值问题:

    f(t)4f(t)+4f(t)=e2t,f(0)=0,f(0)=3.f''(t)-4f'(t)+4f(t)=e^{2t},\qquad f(0)=0,\qquad f'(0)=3.
  3. 求解初值问题:

    f(t)+9f(t)=6cos3t,f(0)=1,f(0)=3.f''(t)+9f(t)=6\cos 3t,\qquad f(0)=1,\qquad f'(0)=3.

Kai

1.

(1)

特性方程式は

r2r6=(r3)(r+2)=0r^2-r-6=(r-3)(r+2)=0

なので、 f=C1e3t+C2e2tf=C_1e^{3t}+C_2e^{-2t} である。初期条件は

C1+C2=1,3C12C2=2C_1+C_2=1,\qquad 3C_1-2C_2=2

を与える。したがって C1=4/5, C2=1/5C_1=4/5,\ C_2=1/5 であり、

f(t)=45e3t+15e2t\boxed{f(t)=\frac45e^{3t}+\frac15e^{-2t}}

である。

(2)

左辺は (D2)2f(D-2)^2f である。 f=e2tg(t)f=e^{2t}g(t) とおけば

(D2)2f=e2tg=e2t,(D-2)^2f=e^{2t}g''=e^{2t},

したがって g=1g''=1 である。よって g=t2/2+C1t+C2g=t^2/2+C_1t+C_2 である。 f(0)=0f(0)=0 より C2=0C_2=0 、また f=e2t(2g+g)f'=e^{2t}(2g+g')f(0)=3f'(0)=3 より C1=3C_1=3 である。ゆえに

f(t)=e2t(3t+t22)\boxed{f(t)=e^{2t}\left(3t+\frac{t^2}{2}\right)}

となる。

(3)

同次解は C1cos3t+C2sin3tC_1\cos3t+C_2\sin3t である。右辺が固有振動数と共振するため、特解を fp=Atsin3tf_p=At\sin3t とおく。このとき

fp+9fp=6Acos3tf_p''+9f_p=6A\cos3t

なので A=1A=1 である。したがって

f=C1cos3t+C2sin3t+tsin3t.f=C_1\cos3t+C_2\sin3t+t\sin3t.

f(0)=1f(0)=1 より C1=1C_1=1f(0)=3f'(0)=3 より 3C2=33C_2=3 、すなわち C2=1C_2=1 である。よって

f(t)=cos3t+sin3t+tsin3t\boxed{f(t)=\cos3t+\sin3t+t\sin3t}

である。各結果は実数だけで表され、直接微分して代入すると方程式と初期条件の双方を満たす。