huangziyuan212 发表于 2013-9-13 22:30

转子跌落仿真程序

求解这个微分方程。这里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=;


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



页: [1]
查看完整版本: 转子跌落仿真程序