请高手帮忙谢谢:一个多自由的系统
我要处理的一个多自由的系统，所要求的三个速度和位移，我的m文件：
F=
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=;
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=;
tol=0.001;
trace=0;
=ode23('zhendongjuzhen2',,x0,tol,trace);
d=x(:,1);
v=x(:,2) 你的m文件明显不对头，先看看有关m文件的知识吧
[ 本帖最后由 咕噜噜 于 2007-5-10 13:13 编辑 ]
可以指导一下吗
我的m文件哪里有问题，可以指导一下吗
页:
[1]