跳到主要内容

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

Author

思齐塾, 祭音Myyura

Description

以下の連立微分方程式の解 x(t)x(t) , y(t)y(t) を求めよ。 ただし、初期値は x(0)=0x(0)=0 , y(0)=1y(0) = 1 とする。

{dxdt=4y6e3tsin2tdydt=x3e3tcos2t\begin{cases} \frac{dx}{dt} = 4y - 6e^{-3t}\sin{2t} \\ \frac{dy}{dt} = -x - 3e^{-3t}\cos{2t} \end{cases}

题目描述

求下列线性微分方程组的解 x(t)x(t)y(t)y(t)

{dxdt=4y6e3tsin2t,dydt=x3e3tcos2t,\begin{cases} \dfrac{dx}{dt}=4y-6e^{-3t}\sin 2t,\\ \dfrac{dy}{dt}=-x-3e^{-3t}\cos 2t, \end{cases}

并满足初值

x(0)=0,y(0)=1.x(0)=0,\qquad y(0)=1.

Kai

Let X(s)=L[x(t)]X(s) = L[x(t)] and Y(s)=L[y(t)]Y(s) = L[y(t)] . Taking the Laplace transform of the given system of differential equations, we get:

sX(s)x(0)=4Y(s)6L[e3tsin2t]sX(s) - x(0) = 4Y(s) - 6L[e^{-3t}\sin{2t}] sY(s)y(0)=X(s)3L[e3tcos2t]sY(s) - y(0) = -X(s) - 3L[e^{-3t}\cos{2t}]

Using the initial conditions x(0)=0x(0) = 0 and y(0)=1y(0) = 1 , and the Laplace transforms L[e3tsin2t]=2(s+3)2+4L[e^{-3t}\sin{2t}] = \frac{2}{(s+3)^2 + 4} and L[e3tcos2t]=s+3(s+3)2+4L[e^{-3t}\cos{2t}] = \frac{s+3}{(s+3)^2 + 4} , we have:

sX(s)=4Y(s)12(s+3)2+4sX(s) = 4Y(s) - \frac{12}{(s+3)^2 + 4} sY(s)1=X(s)3(s+3)(s+3)2+4sY(s) - 1 = -X(s) - \frac{3(s+3)}{(s+3)^2 + 4}

Multiplying the second equation by 4 gives: 4sY(s)4=4X(s)12(s+3)(s+3)2+44sY(s) - 4 = -4X(s) - \frac{12(s+3)}{(s+3)^2 + 4}

From the first equation, Y(s)=sX(s)4+3(s+3)2+4Y(s) = \frac{sX(s)}{4} + \frac{3}{(s+3)^2 + 4} . Substituting this into the above equation:

4s(sX(s)4+3(s+3)2+4)4=4X(s)12(s+3)(s+3)2+44s(\frac{sX(s)}{4} + \frac{3}{(s+3)^2 + 4}) - 4 = -4X(s) - \frac{12(s+3)}{(s+3)^2 + 4} s2X(s)+12s(s+3)2+44=4X(s)12(s+3)(s+3)2+4s^2X(s) + \frac{12s}{(s+3)^2 + 4} - 4 = -4X(s) - \frac{12(s+3)}{(s+3)^2 + 4} (s2+4)X(s)=412s(s+3)2+412(s+3)(s+3)2+4(s^2 + 4)X(s) = 4 - \frac{12s}{(s+3)^2 + 4} - \frac{12(s+3)}{(s+3)^2 + 4} (s2+4)X(s)=412(2s+3)(s2+6s+13)(s^2 + 4)X(s) = 4 - \frac{12(2s+3)}{(s^2 + 6s + 13)} X(s)=4(s2+6s+13)12(2s+3)(s2+4)(s2+6s+13)=4s2+24s+5224s36(s2+4)(s2+6s+13)=4s2+16(s2+4)(s2+6s+13)=4(s2+4)(s2+4)(s2+6s+13)=4s2+6s+13=4(s+3)2+4X(s) = \frac{4(s^2 + 6s + 13) - 12(2s+3)}{(s^2 + 4)(s^2 + 6s + 13)} = \frac{4s^2 + 24s + 52 - 24s - 36}{(s^2 + 4)(s^2 + 6s + 13)} = \frac{4s^2 + 16}{(s^2 + 4)(s^2 + 6s + 13)} = \frac{4(s^2+4)}{(s^2 + 4)(s^2 + 6s + 13)} = \frac{4}{s^2 + 6s + 13} = \frac{4}{(s+3)^2 + 4}

Therefore, X(s)=4(s+3)2+4=22(s+3)2+4X(s) = \frac{4}{(s+3)^2 + 4} = 2\frac{2}{(s+3)^2 + 4} , so x(t)=2e3tsin2tx(t) = 2e^{-3t}\sin{2t} .

Y(s)=sX(s)4+3(s+3)2+4=s(4(s+3)2+4)4+3(s+3)2+4=s+3(s+3)2+4Y(s) = \frac{sX(s)}{4} + \frac{3}{(s+3)^2 + 4} = \frac{s(\frac{4}{(s+3)^2 + 4})}{4} + \frac{3}{(s+3)^2 + 4} = \frac{s + 3}{(s+3)^2 + 4} So y(t)=e3tcos2ty(t) = e^{-3t}\cos{2t} .

Therefore, x(t)=2e3tsin2tx(t) = 2e^{-3t}\sin{2t} and y(t)=e3tcos2ty(t) = e^{-3t}\cos{2t} .

Verification: dxdt=2(3e3tsin2t+2e3tcos2t)=6e3tsin2t+4e3tcos2t\frac{dx}{dt} = 2(-3e^{-3t}\sin{2t} + 2e^{-3t}\cos{2t}) = -6e^{-3t}\sin{2t} + 4e^{-3t}\cos{2t} 4y6e3tsin2t=4e3tcos2t6e3tsin2t4y - 6e^{-3t}\sin{2t} = 4e^{-3t}\cos{2t} - 6e^{-3t}\sin{2t}

dydt=3e3tcos2t2e3tsin2t\frac{dy}{dt} = -3e^{-3t}\cos{2t} - 2e^{-3t}\sin{2t}

x3e3tcos2t=2e3tsin2t3e3tcos2t-x - 3e^{-3t}\cos{2t} = -2e^{-3t}\sin{2t} - 3e^{-3t}\cos{2t} x(0)=2e0sin0=0x(0) = 2e^0\sin{0} = 0 y(0)=e0cos0=1y(0) = e^0\cos{0} = 1