|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
求解这个微分方程。这里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
|
|