chriswonder 发表于 2008-5-18 15:53

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

我用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 发表于 2008-5-18 16:03

原帖由 chriswonder 于 2008-5-18 15:53 发表 http://www.chinavib.com/forum/images/common/back.gif
我用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 发表于 2008-5-18 16:08

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

也没有zero这个命令

[ 本帖最后由 sigma665 于 2008-5-18 16:09 编辑 ]
页: [1]
查看完整版本: 请教我编的M函数老出问题请看下是什么原因