请教我编的M函数老出问题请看下是什么原因

chriswonder

我用ODE45做1个24维的2阶微分方程
function dy=weifen(t,x)
dy=zeros(48,1);
if t<=6
    A=c;
else
    A=zero(24,0);
dy(1:24)=x(25:48);
dy(25:48)=-inv(M)*M*A-inv(M)*k*x(1:24);
end
里面C是已知的矩阵,但是我的这个M函数在MATLAB里面运行的时候出现Error in ==> comput at 5
if (t<=6)不知道是怎么回事

ch_j1985

原帖由 chriswonder 于 2008-5-18 15:53 发表

                               
登录/注册后可看大图

我用ODE45做1个24维的2阶微分方程
function dy=weifen(t,x)
dy=zeros(48,1);
if t comput at 5
if (t

首先，传递的参数t是数组还是一个数字？
其次，A=zero(24,0);是不是有错误？
可以根据以上的两个问题，自己调试一下

sigma665

A=zero(24,0);
恩，这个有问题
24行，0列？

也没有zero这个命令

[ 本帖最后由 sigma665 于 2008-5-18 16:09 编辑 ]