[問題] ODE45 解法器已回收
主程式
[t,x]=ode45('sspace',0:(L0/V0)/200:L0/V0,[0 0 0 0 0 0 0 0])
這個ODE我已經先解出 t與x的矩陣
假設t算出來是201X1的矩陣，x是201X8的矩陣
之後我又要解第2組
[t1,x1]=ode45('sspace1',0:(L0/V0)/200:L0/V0,[0 0 0 0 0 0 0 0]);
因為我第2組的ode需要用到第一組算出來的x，也就是第2組跑第一個時間時
我要用到x的第一列，跑第2個時間時要用到第2列
-----------------------------------------------------------------
以下為 sspace1.m 的一部分
function dx=sspace1(t1,X1)
global x
.....
.....
F_matrix1=[0,-M0*g*W_n1(1,V0*t1),-M0*g*W_n1(2,V0*t1),-M0*g*W_n1(3,V0*t1)]'
-C_matrix1*x(第t1列,5:8)'-K_matrix1*x(第t1列,1:4)'
......
當時間跑到t1時，我的列數想要是第一組t的位置 t=t1
我上面的第t1列要怎麼寫?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.57.150.174