[問題] 解ODE的問題已回收
各位好:
我的原始方程式:[M][y'']+[k][y]=[f11]*V
整理成[y'']=[M]^-1[f(t)]-[M]^-1[k][y]
[M]: 6*6
[k]: 6*6
[f11] :6*1
所以是六條聯立解ODE ,我參考了書本上的範例(一階2*2)
以下都是我自己想的,實在是程式苦手~.~"
=========================================
我自己降階
x1=y ,x2=y'
x2'=[M]^-1[f(t)]-[M]^-1[k]x1
x1'=x2 這邊的x1 x2 x1' x2'他會自己判斷成矩陣吧?
那請問我副程式可以寫成這樣嗎 檔名pp.m
function xdot=pp(t,x)
globle w
mainprogram; %這一條是另一個程式,跑出M,K的參數用的
if t<.1
V=sin(w*t)
else
V=0 %希望在t<0.1實有外力作用
end
xdot(1)=x(2);
xdot(2)=M\f11*V-M\K*x(1)
xdot=[xdot(1);xdot(2)];
然後主程式解
globle w
w=5;
[t,x]=ode23('pp',[0,5],[0,0,0,0,0,0]) ==>這個初始條件在這邊我不太會放
[0,0]還是[0,0,0,0,0,0]?
plot(t,x)
我當然知道我跑不出來,我只是想請問各位高手我那邊觀念有問題?
單純一條ODE我會做,n*n矩陣我就不太會寫,我程式想全部用矩陣去寫
因為我的矩陣維度不算小,麻煩高手替小弟解惑一下,我想的好辛苦@@
還是說有高手寫過的程式可以借我參考!!
謝謝!!!!<(_ _)>
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.66.117
→
03/19 16:48, , 1F
03/19 16:48, 1F
→
03/20 22:03, , 2F
03/20 22:03, 2F
→
03/29 22:20, , 3F
03/29 22:20, 3F