请高手帮忙谢谢:一个多自由的系统

猴王

我要处理的一个多自由的系统，所要求的三个速度和位移，
我的m文件：
F=[25;24;20]
  m1=31.2;
  m2=31.2;
  m3=31.2;
  k1=67.51;
  k2=67.51;
  k3=67.51;
  c1=0.01;
  c2=0.01;
  c3=0.01;
M=[ m1+m2,-m2,0;-m2,m2+m3,-m3;0,-m3,m3];
B=[ c1+c2,-c2,0;-c2,c2+c3,-c3;0,-c3,c3];
k=[k1+k2,-k2,0;-k2,k2+k3,-k3;0,-k3,k3];
A=zeros(2*3);
A(1:3,1:3)=zeros(3);
A(1:3,3+1:end)=eye(3);
A(3+1:end,1:3)=-inv(M)*K;
A(3+1:end,3+1:end)=-inv(M)*C;
B=zeros(2*3,1);
B(1:3)=zeros(3);
B(3+1:2*3)=F;
%f-作用力的向量，%n为系统的自由度
xdot=Ax+B
我的执行文件：
t0=0;
   tfinal=3;
   x0=[0,0];
   tol=0.001;
   trace=0;
   [t,x]=ode23('zhendongjuzhen2',[t0,tfinal],x0,tol,trace);
   d=x(:,1);
   v=x(:,2)

咕噜噜

你的m文件明显不对头，先看看有关m文件的知识吧

[ 本帖最后由 咕噜噜 于 2007-5-10 13:13 编辑 ]

猴王

我的m文件哪里有问题，可以指导一下吗