转子跌落仿真程序

huangziyuan212

求解这个微分方程。这里Fn，Ft都含有xy的函数
我需要画出横坐标x(-0.2,0.2),纵坐标y(-0.2,0.2)的图。可以我画出之后图是空白的。希望大侠能帮小弟看一下程序问题
M函数为 Guiji_fun.m
function f=Guiji_fun1(t,y)
         g=9.8;
         k=3.378e8;
         b=10/9;
         c=1000;
         m=1.1;
         Jr=1.9*10^(-4);
         s=1.25*10^(-4);
         r=5.875*10^(-3);
         u=1;
         F(1)=k*(sqrt(y(1)^2+y(2)^2)-s)^b+c*((y(2)*y(1)+y(3)*y(4))/(sqrt((y(1))^2)+(y(3))^2));
         F(2)=u*F(1);%F(t)
f=[y(2);-(F(1)*cos(atan(y(2)/y(1)))+F(2)*sin(atan(y(2)/y(1))))/m;y(4);-g+(F(2)*cos(atan(y(2)/y(1)))-F(1)*sin(atan(y(2)/y(1))))/m;y(6);-(F(2)*r)/Jr];


窗口执行
>> y0=[0,0,0,0,0,6000];
>> [t,y]=ode45(@Guiji_fun1,[-0.2,0.2],y0);
>> plot(y(:,1),y(:,3))
Warning: Imaginary parts of complex X and/or Y arguments ignored