作者kevio (成仔)
看板MATLAB
标题[讨论] 请问用dsolve解ODE
时间Sat Mar 23 00:59:50 2013
我想用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
1F：推 iHakka:[y how]=simple(y)化简後做diff(y,2)+9*y我得到sin4t 03/24 00:54
2F：→ kevio:谢谢你的回答 我回去试试~~ 03/24 13:05
3F：→ kevio:原来....是我算错了 orz.... 03/25 12:01