如何解决ode解微分方程出现虚解情况？
微分方程如下：不知是参数问题还是方程问题，用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; ...
]; 看不懂……看着就想睡觉……
页:
[1]