如何解决ode解微分方程出现虚解情况？

lhy

微分方程如下：不知是参数问题还是方程问题，用ode45解微分方程总会出现虚数解的情况，请各位高手指点~~~
dx=[x(2);...
    -2*ksai/lan*x(2)-k11*(x(1)-x(5))-Fx1/C1/m1/wr^2;...
    x(4);...
    -2*ksai/lan*x(4)+k21*x(1)-(k21+4*k22)*x(3)+4*k22*x(5)+e2*cos(t);...
    x(6);...
    -2*ksai/lan*x(6)+4*k32*x(2)-(4*k32+k33)*x(5)+k33*x(7)+e3*cos(t+fai);...
    x(8);...
    -2*ksai/lan*x(8)+k43*(x(5)-x(7))-Fx4/C2/m4/wr^2; ...
     x(10);...
    -2*ksai/lan*x(10)-k11*(x(9)-x(13))-Fy1/C1/m1/wr^2-g/C1/wr^2;...
    x(12);...
    -2*ksai/lan*x(12)+k21*x(9)-(k21+4*k22)*x(11)+4*k22*x(13)+e2*sign(t)-g/C1/wr^2;...
    x(14);...
    -2*ksai/lan*x(14)+4*k32*x(10)-(4*k32+k33)*x(13)+k33*x(15)+e3*sign(t+fai)-g/C2/wr^2;...
    x(16);...
    -2*ksai/lan*x(16)+k43*(x(13)-x(15))-Fy4/C2/m4/wr^2-g/C2/wr^2; ...
  ];

dongyuanzhixing

看不懂……  看着就想睡觉……