|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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)
可以运行会有很多错误,好像矩阵的维数不对应,不知道是状态方程写的不对还是初始值不对?请高手指点一下 |
|