转子跌落仿真程序
求解这个微分方程。这里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]