推 iHakka:[y how]=simple(y)化簡後做diff(y,2)+9*y我得到sin4t 03/24 00:54
→ kevio:謝謝你的回答 我回去試試~~ 03/24 13:05
→ kevio:原來....是我算錯了 orz.... 03/25 12:01
我想用Matlab symbolic function 來解一個簡單ODE
y'' + 9y = sin (4t) , initial condition y=y'=0
用手算應該是 y= (-1/7) * sin (4t)
但Matlab 給我相當神奇的答案
Code
y=dsolve('D2y+9*y=sin(4*t)','y(0)=0', 'Dy(0)=0','t')
-->得到
y = (4*sin(3*t))/21 - sin(3*t)*(cos(7*t)/42 + cos(t)/6)
+ cos(3*t)*(sin(7*t)/42 - sin(t)/6)
然後我用diff 函式 做一階 二階微分 帶回去根本就無法讓等號成立
請問是哪邊出錯?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 111.255.170.128