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