声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 760|回复: 2

[综合讨论] 请高手帮忙谢谢:一个多自由的系统

[复制链接]
发表于 2007-5-10 10:02 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我要处理的一个多自由的系统,所要求的三个速度和位移,
我的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)
回复
分享到:

使用道具 举报

发表于 2007-5-10 13:10 | 显示全部楼层
你的m文件明显不对头,先看看有关m文件的知识吧

[ 本帖最后由 咕噜噜 于 2007-5-10 13:13 编辑 ]
 楼主| 发表于 2007-5-10 15:03 | 显示全部楼层

可以指导一下吗

我的m文件哪里有问题,可以指导一下吗
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-12 21:36 , Processed in 0.054364 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表