[問題] 解多次微分方程式已回收
我想問一下 我有一個方程式為 m*D2u+c*Du^a+k*u=0
我在matlab寫
eq='m*D2u+c*(Du)^a+k*u=0'
int='u(0)=0,Du(0)=0'
dsolve(eq,int,'t')
請問我還缺少啥，還是說要改成怎樣寫才求的出u??
感謝
另外 我若想求得對u微分兩次的數值
我的方程式假如為 20*D2u+2*Du+100*u=60
在M檔案裡我寫成
function du=mck(t,u)
du(1)=u(2);
du(2)=(-100*u(1)-2*u(2)+60)/20;
在主程式裡輸入
>> dt=1/40;n=250;tspan=linspace(0,(n-1)*dt,n);
>> [t,z]=ode45('mck',tspan,[0 0])
還要再輸入啥 才能得到微分兩次 也就是du(2)的值
感謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.13.113
※ 編輯: tarepanda275 來自: 112.104.53.214 (08/10 23:35)
※ 編輯: tarepanda275 來自: 112.104.53.214 (08/10 23:37)
※ 編輯: tarepanda275 來自: 112.104.53.214 (08/10 23:41)
※ 編輯: tarepanda275 來自: 112.104.53.214 (08/10 23:43)
※ 編輯: tarepanda275 來自: 112.104.53.214 (08/10 23:44)
推
08/11 00:17, , 1F
08/11 00:17, 1F
討論串 (同標題文章)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
