跳到主要内容

九州大学 システム情報科学府 情報学専攻・情報知能工学専攻・電気電子工学専攻 共通 2012年度入学 数学 問2(微分方程式)

Author​

思齐塾, 祭音Myyura

Description​

2つの関数 x(t)x(t) , y(t)y(t) について、次の連立微分方程式を解け.

{dxdt=−2x−y+cos⁡tdydt=−dxdt−6x\begin{cases} \frac{dx}{dt} = -2x - y + \cos t \\ \frac{dy}{dt} = - \frac{dx}{dt} - 6x \end{cases}

题目描述​

设 x(t)x(t)、y(t)y(t) 为关于 tt 的函数,求下列非齐次联立微分方程的通解:

{dxdt=−2x−y+cos⁡t,dydt=−dxdt−6x.\begin{cases} \dfrac{dx}{dt}=-2x-y+\cos t,\\[2mm] \dfrac{dy}{dt}=-\dfrac{dx}{dt}-6x. \end{cases}

Kai​

Let's rewrite the system of equations as:

{dxdt=−2x−y+cos⁡tdydt=−dxdt−6x\begin{cases} \frac{dx}{dt} = -2x - y + \cos t \\ \frac{dy}{dt} = - \frac{dx}{dt} - 6x \end{cases}

Substitute the first equation into the second one:

dydt=−(−2x−y+cos⁡t)−6x=2x+y−cos⁡t−6x=−4x+y−cos⁡t\frac{dy}{dt} = - (-2x - y + \cos t) - 6x = 2x + y - \cos t - 6x = -4x + y - \cos t

Now we have:

{dxdt=−2x−y+cos⁡tdydt=−4x+y−cos⁡t\begin{cases} \frac{dx}{dt} = -2x - y + \cos t \\ \frac{dy}{dt} = -4x + y - \cos t \end{cases}

Differentiate the first equation with respect to t:

d2xdt2=−2dxdt−dydt−sin⁡t\frac{d^2x}{dt^2} = -2\frac{dx}{dt} - \frac{dy}{dt} - \sin t

Substitute the first and second equations into this one:

d2xdt2=−2(−2x−y+cos⁡t)−(−4x+y−cos⁡t)−sin⁡t=4x+2y−2cos⁡t+4x−y+cos⁡t−sin⁡t=8x+y−cos⁡t−sin⁡t\frac{d^2x}{dt^2} = -2(-2x - y + \cos t) - (-4x + y - \cos t) - \sin t = 4x + 2y - 2\cos t + 4x - y + \cos t - \sin t = 8x + y - \cos t - \sin t

From the first equation, we have y=−dxdt−2x+cos⁡ty = -\frac{dx}{dt} - 2x + \cos t . Substitute this into the above equation:

d2xdt2=8x+(−dxdt−2x+cos⁡t)−cos⁡t−sin⁡t=6x−dxdt−sin⁡t\frac{d^2x}{dt^2} = 8x + (- \frac{dx}{dt} - 2x + \cos t) - \cos t - \sin t = 6x - \frac{dx}{dt} - \sin t
d2xdt2+dxdt−6x=−sin⁡t\frac{d^2x}{dt^2} + \frac{dx}{dt} - 6x = - \sin t

This is a second-order linear non-homogeneous differential equation. Let's find the homogeneous solution first:

The characteristic equation is r2+r−6=0r^2 + r - 6 = 0 , which factors to (r+3)(r−2)=0(r+3)(r-2) = 0 . Thus, r1=−3r_1 = -3 and r2=2r_2 = 2 .

The homogeneous solution is xh(t)=c1e−3t+c2e2tx_h(t) = c_1e^{-3t} + c_2e^{2t} .

Now, let's find a particular solution. Since the right-hand side is −sin⁡t-\sin t , we can assume a particular solution of the form xp(t)=Acos⁡t+Bsin⁡tx_p(t) = A\cos t + B\sin t .

dxpdt=−Asin⁡t+Bcos⁡t\frac{dx_p}{dt} = -A\sin t + B\cos t
d2xpdt2=−Acos⁡t−Bsin⁡t\frac{d^2x_p}{dt^2} = -A\cos t - B\sin t

Substitute these into the differential equation:

(−Acos⁡t−Bsin⁡t)+(−Asin⁡t+Bcos⁡t)−6(Acos⁡t+Bsin⁡t)=−sin⁡t(-A\cos t - B\sin t) + (-A\sin t + B\cos t) - 6(A\cos t + B\sin t) = -\sin t
(−A+B−6A)cos⁡t+(−B−A−6B)sin⁡t=−sin⁡t(-A + B - 6A)\cos t + (-B - A - 6B)\sin t = -\sin t
(−7A+B)cos⁡t+(−A−7B)sin⁡t=−sin⁡t(-7A + B)\cos t + (-A - 7B)\sin t = -\sin t

So we have the system:

{−7A+B=0−A−7B=−1\begin{cases} -7A + B = 0 \\ -A - 7B = -1 \end{cases}

From the first equation, B=7AB = 7A . Substitute this into the second equation:

−A−7(7A)=−1  ⟹  −A−49A=−1  ⟹  −50A=−1  ⟹  A=150-A - 7(7A) = -1 \implies -A - 49A = -1 \implies -50A = -1 \implies A = \frac{1}{50}

Thus, B=7A=750B = 7A = \frac{7}{50} .

The particular solution is xp(t)=150cos⁡t+750sin⁡tx_p(t) = \frac{1}{50}\cos t + \frac{7}{50}\sin t .

The general solution for x(t) is x(t)=c1e−3t+c2e2t+150cos⁡t+750sin⁡tx(t) = c_1e^{-3t} + c_2e^{2t} + \frac{1}{50}\cos t + \frac{7}{50}\sin t .

To find y(t), we use y=−dxdt−2x+cos⁡ty = -\frac{dx}{dt} - 2x + \cos t :

dxdt=−3c1e−3t+2c2e2t−150sin⁡t+750cos⁡t\frac{dx}{dt} = -3c_1e^{-3t} + 2c_2e^{2t} - \frac{1}{50}\sin t + \frac{7}{50}\cos t
y(t)=−(−3c1e−3t+2c2e2t−150sin⁡t+750cos⁡t)−2(c1e−3t+c2e2t+150cos⁡t+750sin⁡t)+cos⁡ty(t) = -(-3c_1e^{-3t} + 2c_2e^{2t} - \frac{1}{50}\sin t + \frac{7}{50}\cos t) - 2(c_1e^{-3t} + c_2e^{2t} + \frac{1}{50}\cos t + \frac{7}{50}\sin t) + \cos t
y(t)=3c1e−3t−2c2e2t+150sin⁡t−750cos⁡t−2c1e−3t−2c2e2t−125cos⁡t−725sin⁡t+cos⁡ty(t) = 3c_1e^{-3t} - 2c_2e^{2t} + \frac{1}{50}\sin t - \frac{7}{50}\cos t - 2c_1e^{-3t} - 2c_2e^{2t} - \frac{1}{25}\cos t - \frac{7}{25}\sin t + \cos t
y(t)=c1e−3t−4c2e2t−1350sin⁡t+4150cos⁡ty(t) = c_1e^{-3t} - 4c_2e^{2t} - \frac{13}{50}\sin t + \frac{41}{50}\cos t

Therefore:

{x(t)=c1e−3t+c2e2t+150cos⁡t+750sin⁡ty(t)=c1e−3t−4c2e2t−1350sin⁡t+4150cos⁡t\begin{cases} x(t) = c_1e^{-3t} + c_2e^{2t} + \frac{1}{50}\cos t + \frac{7}{50}\sin t \\ y(t) = c_1e^{-3t} - 4c_2e^{2t} - \frac{13}{50}\sin t + \frac{41}{50}\cos t \end{cases}