[問題] 請問用dsolve解ODE

看板MATLAB作者 (成仔)時間13年前 (2013/03/23 00:59), 編輯推噓1(102)
留言3則, 2人參與, 最新討論串1/1
我想用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

03/24 00:54, , 1F
[y how]=simple(y)化簡後做diff(y,2)+9*y我得到sin4t
03/24 00:54, 1F

03/24 13:05, , 2F
謝謝你的回答 我回去試試~~
03/24 13:05, 2F

03/25 12:01, , 3F
原來....是我算錯了 orz....
03/25 12:01, 3F
文章代碼(AID): #1HJ8u9IE (MATLAB)
文章代碼(AID): #1HJ8u9IE (MATLAB)