跳到主要内容

九州大学 システム情報科学府 共通 2011年8月実施 微分方程

Author

思齐塾, 祭音Myyura

Description

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

{dxdt=2xy+costdydt=dxdt6x\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=2xy+cost,dydt=dxdt6x.\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=2xy+costdydt=dxdt6x\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=(2xy+cost)6x=2x+ycost6x=4x+ycost\frac{dy}{dt} = - (-2x - y + \cos t) - 6x = 2x + y - \cos t - 6x = -4x + y - \cos t

Now we have:

{dxdt=2xy+costdydt=4x+ycost\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=2dxdtdydtsint\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(2xy+cost)(4x+ycost)sint=4x+2y2cost+4xy+costsint=8x+ycostsint\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=dxdt2x+costy = -\frac{dx}{dt} - 2x + \cos t . Substitute this into the above equation:

d2xdt2=8x+(dxdt2x+cost)costsint=6xdxdtsint\frac{d^2x}{dt^2} = 8x + (- \frac{dx}{dt} - 2x + \cos t) - \cos t - \sin t = 6x - \frac{dx}{dt} - \sin t
d2xdt2+dxdt6x=sint\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+r6=0r^2 + r - 6 = 0 , which factors to (r+3)(r2)=0(r+3)(r-2) = 0 . Thus, r1=3r_1 = -3 and r2=2r_2 = 2 .

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

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

dxpdt=Asint+Bcost\frac{dx_p}{dt} = -A\sin t + B\cos t
d2xpdt2=AcostBsint\frac{d^2x_p}{dt^2} = -A\cos t - B\sin t

Substitute these into the differential equation:

(AcostBsint)+(Asint+Bcost)6(Acost+Bsint)=sint(-A\cos t - B\sin t) + (-A\sin t + B\cos t) - 6(A\cos t + B\sin t) = -\sin t
(A+B6A)cost+(BA6B)sint=sint(-A + B - 6A)\cos t + (-B - A - 6B)\sin t = -\sin t
(7A+B)cost+(A7B)sint=sint(-7A + B)\cos t + (-A - 7B)\sin t = -\sin t

So we have the system:

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

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

A7(7A)=1    A49A=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)=150cost+750sintx_p(t) = \frac{1}{50}\cos t + \frac{7}{50}\sin t .

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

To find y(t), we use y=dxdt2x+costy = -\frac{dx}{dt} - 2x + \cos t :

dxdt=3c1e3t+2c2e2t150sint+750cost\frac{dx}{dt} = -3c_1e^{-3t} + 2c_2e^{2t} - \frac{1}{50}\sin t + \frac{7}{50}\cos t
y(t)=(3c1e3t+2c2e2t150sint+750cost)2(c1e3t+c2e2t+150cost+750sint)+costy(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)=3c1e3t2c2e2t+150sint750cost2c1e3t2c2e2t125cost725sint+costy(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)=c1e3t4c2e2t1350sint+4150costy(t) = c_1e^{-3t} - 4c_2e^{2t} - \frac{13}{50}\sin t + \frac{41}{50}\cos t

Therefore:

{x(t)=c1e3t+c2e2t+150cost+750sinty(t)=c1e3t4c2e2t1350sint+4150cost\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}