九州大学 システム情報科学府 共通 2011年8月実施 微分方程
Author
思齐塾, 祭音Myyura
Description
2つの関数 x(t) , y(t) について、次の連立微分方程式を解け.
{dtdx=−2x−y+costdtdy=−dtdx−6x
题目描述
设 x(t)、y(t) 为关于 t 的函数,求下列非齐次联立微分方程的通解:
⎩⎨⎧dtdx=−2x−y+cost,dtdy=−dtdx−6x.
Kai
Let's rewrite the system of equations as:
{dtdx=−2x−y+costdtdy=−dtdx−6x
Substitute the first equation into the second one:
dtdy=−(−2x−y+cost)−6x=2x+y−cost−6x=−4x+y−cost
Now we have:
{dtdx=−2x−y+costdtdy=−4x+y−cost
Differentiate the first equation with respect to t:
dt2d2x=−2dtdx−dtdy−sint
Substitute the first and second equations into this one:
dt2d2x=−2(−2x−y+cost)−(−4x+y−cost)−sint=4x+2y−2cost+4x−y+cost−sint=8x+y−cost−sint
From the first equation, we have y=−dtdx−2x+cost . Substitute this into the above equation:
dt2d2x=8x+(−dtdx−2x+cost)−cost−sint=6x−dtdx−sint
dt2d2x+dtdx−6x=−sint
This is a second-order linear non-homogeneous differential equation. Let's find the homogeneous solution first:
The characteristic equation is r2+r−6=0 , which factors to (r+3)(r−2)=0 . Thus, r1=−3 and r2=2 .
The homogeneous solution is xh(t)=c1e−3t+c2e2t .
Now, let's find a particular solution. Since the right-hand side is −sint , we can assume a particular solution of the form xp(t)=Acost+Bsint .
dtdxp=−Asint+Bcost
dt2d2xp=−Acost−Bsint
Substitute these into the differential equation:
(−Acost−Bsint)+(−Asint+Bcost)−6(Acost+Bsint)=−sint
(−A+B−6A)cost+(−B−A−6B)sint=−sint
(−7A+B)cost+(−A−7B)sint=−sint
So we have the system:
{−7A+B=0−A−7B=−1
From the first equation, B=7A . Substitute this into the second equation:
−A−7(7A)=−1⟹−A−49A=−1⟹−50A=−1⟹A=501
Thus, B=7A=507 .
The particular solution is xp(t)=501cost+507sint .
The general solution for x(t) is x(t)=c1e−3t+c2e2t+501cost+507sint .
To find y(t), we use y=−dtdx−2x+cost :
dtdx=−3c1e−3t+2c2e2t−501sint+507cost
y(t)=−(−3c1e−3t+2c2e2t−501sint+507cost)−2(c1e−3t+c2e2t+501cost+507sint)+cost
y(t)=3c1e−3t−2c2e2t+501sint−507cost−2c1e−3t−2c2e2t−251cost−257sint+cost
y(t)=c1e−3t−4c2e2t−5013sint+5041cost
Therefore:
{x(t)=c1e−3t+c2e2t+501cost+507sinty(t)=c1e−3t−4c2e2t−5013sint+5041cost