作者ytyty (该换个版潜水了™ )
看板Grad-ProbAsk
标题Re: [理工] [工数]-拉式常系数ode
时间Fri Jan 29 18:46:21 2010
※ 引述《smallprawn (水中瑕)》之铭言：
: t-1, 1<t<2
: 3.y''+4y'+3y=r(t) y(0)=2, y'(0)=-2, r(t)={
: 0, otherwise
: 请问这题要怎麽算?!
: 麻烦哩!!
r(t)=(t-1)[H(t-1)-H(t-2)]
y''+4y'+3y=r(t)
令￡[y(t)]=Y(s)
原式取Laplace transform得
2 -s -2s
[s Y-sy(0)-y'(0)]+4[sY-y(0)]+3Y=e ￡[(t-1)+1] - e ￡[(t-1)+2]
2 -s 1 -2s 1 1
(s+4s+3)Y-2s+2-8=e ── - e (── + ──)
s^2 s^2 s
1 -s -2s 1 -2s
(s+1)(s+3)Y= 2s+6 + ──(e - e ) - ── e
s^2 s
2 1 -s -2s 1 -2s
Y= ── + ─────── (e - e ) - ─────── e
s+1 s^2(s+1)(s+3) s(s+1)(s+3)
1 A B C D
令 ─────── = ── + ─── + ─── + ───
s^2(s+1)(s+3) s s^2 s+1 s+3
去分母得
1 = As(s+1)(s+3) + B(s+1)(s+3) + Cs^2(s+3) + Ds^2(s+1)
s=-1代入得 1=C[(-1)^2]*2 , C=1/2
s=-3代入得 1=D[(-3)^2]*(-2), D=-1/18
s=0 代入得 1=B*1*3, B=1/3
等号左式之s^3之系数为0，等号右式之s^3之系数为A+C+D
所以 A+C+D=0 即A+(1/2)+(-1/18)=0, A=-4/9
1 E F G
令 ─────── = ── + ─── + ───
s(s+1)(s+3) s s+1 s+3
去分母得
1 = E(s+1)(s+3) + Fs(s+3) + Gs(s+1)
s=-1代入得 1=F(-1)*2 , F=-1/2
s=-3代入得 1=G(-3)*(-2), G=1/6
s=0 代入得 1=E*1*3, E=1/3
所以
2 -4/9 1/3 1/2 -1/18 -s -2s
Y= ── +( ── + ─── + ─── + ───) (e - e )
s+1 s s^2 s+1 s+3
1/3 -1/2 1/6 -2s
-(── + ─── + ─── )e
s s+1 s+3
2 -4/9 1/3 1/2 -1/18 -s
= ── +( ── + ─── + ─── + ───) e
s+1 s s^2 s+1 s+3
1/9 -1/3 -1/9 -2s
+(── + ─── + ───)e
s s^2 s+3
-1 -t -(t-1) -3(t-1)
y(t)=￡ [Y(s)] = 2e +[-4/9+(1/3)(t-1)+(1/2)e -(1/18)e ]H(t-1)
-3(t-2)
+[1/9-(1/3)(t-2)-(1/9)e ]H(t-2)
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.114.244.49
1F：推 smallprawn:m(_ _)m 辛苦了..感谢.. 01/29 19:58