作者Honor1984 (希望愿望成真)
看板Grad-ProbAsk
标题Re: [理工] [工数]-O.D.E.
时间Tue Dec 15 22:10:22 2009
※ 引述《smallprawn (水中瑕)》之铭言:
: 2 2
: y''+y=1-(t / π ) 0≦t≦π
: I.C. y(0)=y'(0)=0
: 这个我知道前面 y = c1 cos t + c2 sin t
: h
[解一]
y' = (-c1 sint + c2 cost) + (c1' cost + c2' sint)
^^^^^^^^^^^^^^^^^^^^^ = 0
y'' = (-c1 cost - c2 sint) + (-c1'sint +c2' cost)
代回微分方程式
c1'sint + c2'cost = 1 - (t / π )^2
c1'cost + c2'sint = 0
c1(0) = 0
c2(0) = 0
分别解出c1' c2'
=> 解出 c1 c2
请记得这个做法
高阶也是一样
公式解就是这样来的
只是这题不一定好算
---------------------------------
[解二]
y_p为二阶多项式 a_2t^2 + a_1t + a_0
代回原式 = 2a_2 + a_2t^2 + a_1t + a_0 = 1 - (t / π )^2
比一下系数
y_p = -(t/π)^2 + 1 + 2(1/π^2)
----------------------------------
[解三]
y_p
= [D^2 + 1]^(-1) [1-(t / π )^2]
= [1 - D^2][1-(t / π )^2]
= 1-(t / π )^2 + (2/π^2)
-------------------------------------
解出y_p後
y(x) = y_p + y_h满足起始条件y(0) = y'(0) = 0
= -(1+(2 / π^2 ))cos t - (t /π)^2 + 1 + 2/π^2
: 但是後面要怎麽解?!
: 2 2 2 2
: 可得到答案 y(t)=-(1+(2 / π ))cos t - t /π +1 +2/π
: 再次麻烦高手指教了!!感谢!!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.124.105.128
1F:推 smallprawn:请问解一..y'跟y''後面为什麽多了c1'....c2'....?! 12/15 22:33
2F:推 honestonly:参数变异法 c1跟c2是假设的多项式函数 12/15 22:37
3F:推 smallprawn:嗯嗯~懂了!!感谢你@@!! 12/15 22:46