声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1173|回复: 1

[综合] 求教状态方程

[复制链接]
发表于 2009-5-12 20:30 | 显示全部楼层 |阅读模式

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

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

x
求六自由度,带阻尼系统在周期激振力作用下的响应:方程是:MX''+CX'+KX=F(t)
     我的matlab程序如下:
     先定义函数:
       function xdot=odefun(t,x)
                C;
                K;
                M;  %阻尼,质量,刚度矩阵已知,元素就不写了
        F=[1*sin(10*t)  2*sin(10*t)  3*sin(10*t)  4*sin(10*t)  5*sin(10*t)  6*sin(10*t)]'  %激振力列阵
                xdot=[x(2);inv(M)*(F-C-K*x(1));
%***********************************************************************
              下面引用函数:
      t0=0;
           tfinal=3;
           x0=[zeros(6,1);zeros(6,1)];
           [t,x]=ode45(@odefun,[t0,tfinal],x0)
可以运行会有很多错误,好像矩阵的维数不对应,不知道是状态方程写的不对还是初始值不对?请高手指点一下
回复
分享到:

使用道具 举报

发表于 2009-5-13 13:28 | 显示全部楼层
可以用lsim命令直接求解
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-26 08:37 , Processed in 0.057848 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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