求教状态方程

parapara2639

求六自由度，带阻尼系统在周期激振力作用下的响应：方程是：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)
可以运行会有很多错误，好像矩阵的维数不对应，不知道是状态方程写的不对还是初始值不对？请高手指点一下

yangzj

可以用lsim命令直接求解