parapara2639 发表于 2009-5-12 20:30

求教状态方程

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

yangzj 发表于 2009-5-13 13:28

可以用lsim命令直接求解
页: [1]
查看完整版本: 求教状态方程