matlab Apollo卫星的运动轨迹已知Apollo卫星的运动轨迹 满足下面的方程u0=1/82.43;u1=1-u0,r1=sqrt((x+u0)^2+y^2),r2=sqrt((x-u1)^2+y^2),DDx=2*Dy+x-u1*(x+u0)/r1^3-u0*(x-u1)/r2^3,DDy=-2*Dx+y-u0*y/r1^3-u1*y/r2^3,x(0)=1.2,Dx(0)=0,y(0)=0,D
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/01 03:57:33
![matlab Apollo卫星的运动轨迹已知Apollo卫星的运动轨迹 满足下面的方程u0=1/82.43;u1=1-u0,r1=sqrt((x+u0)^2+y^2),r2=sqrt((x-u1)^2+y^2),DDx=2*Dy+x-u1*(x+u0)/r1^3-u0*(x-u1)/r2^3,DDy=-2*Dx+y-u0*y/r1^3-u1*y/r2^3,x(0)=1.2,Dx(0)=0,y(0)=0,D](/uploads/image/z/7107851-11-1.jpg?t=matlab+Apollo%E5%8D%AB%E6%98%9F%E7%9A%84%E8%BF%90%E5%8A%A8%E8%BD%A8%E8%BF%B9%E5%B7%B2%E7%9F%A5Apollo%E5%8D%AB%E6%98%9F%E7%9A%84%E8%BF%90%E5%8A%A8%E8%BD%A8%E8%BF%B9+%E6%BB%A1%E8%B6%B3%E4%B8%8B%E9%9D%A2%E7%9A%84%E6%96%B9%E7%A8%8Bu0%3D1%2F82.43%3Bu1%3D1-u0%2Cr1%3Dsqrt%28%28x%2Bu0%29%5E2%2By%5E2%29%2Cr2%3Dsqrt%28%28x-u1%29%5E2%2By%5E2%29%2CDDx%3D2%2ADy%2Bx-u1%2A%28x%2Bu0%29%2Fr1%5E3-u0%2A%28x-u1%29%2Fr2%5E3%2CDDy%3D-2%2ADx%2By-u0%2Ay%2Fr1%5E3-u1%2Ay%2Fr2%5E3%2Cx%280%29%3D1.2%2CDx%280%29%3D0%2Cy%280%29%3D0%2CD)
matlab Apollo卫星的运动轨迹已知Apollo卫星的运动轨迹 满足下面的方程u0=1/82.43;u1=1-u0,r1=sqrt((x+u0)^2+y^2),r2=sqrt((x-u1)^2+y^2),DDx=2*Dy+x-u1*(x+u0)/r1^3-u0*(x-u1)/r2^3,DDy=-2*Dx+y-u0*y/r1^3-u1*y/r2^3,x(0)=1.2,Dx(0)=0,y(0)=0,D
matlab Apollo卫星的运动轨迹
已知Apollo卫星的运动轨迹 满足下面的方程
u0=1/82.43;u1=1-u0,
r1=sqrt((x+u0)^2+y^2),
r2=sqrt((x-u1)^2+y^2),
DDx=2*Dy+x-u1*(x+u0)/r1^3-u0*(x-u1)/r2^3,
DDy=-2*Dx+y-u0*y/r1^3-u1*y/r2^3,
x(0)=1.2,Dx(0)=0,y(0)=0,Dy(0)=-1.04935751,
其中 ,试在初值 下分别用ode45命令和Simulink进行求解,并绘制出Apollo位置的 轨迹
matlab Apollo卫星的运动轨迹已知Apollo卫星的运动轨迹 满足下面的方程u0=1/82.43;u1=1-u0,r1=sqrt((x+u0)^2+y^2),r2=sqrt((x-u1)^2+y^2),DDx=2*Dy+x-u1*(x+u0)/r1^3-u0*(x-u1)/r2^3,DDy=-2*Dx+y-u0*y/r1^3-u1*y/r2^3,x(0)=1.2,Dx(0)=0,y(0)=0,D
function dx=appollo(t,x)
mu=1/82.45;
mustar=1-mu;
r1=sqrt((x(1)+mu)^2+x(3)^2);
r2=sqrt((x(1)-mustar)^2+x(3)^2);
dx=[x(2)
2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3
x(4)
-2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3];
--------------------------------------------------------------------------------------
x0=[1.2;0;0;-1.04935751];%x0(i)对应与xi的初值
options=odeset('reltol',1e-8);%该命令的另一种写法是options=odeset;options.reltol=1e-8;
tic
[t,y]=ode45(@appollo,[0,20],x0,options);%t是时间点,y的第i列对应xi的值,t和y的行数相同
toc
plot(y(:,1),y(:,3))%绘制x1和x3,也就是x和y的图形
title('Appollo卫星运动轨迹')
xlabel('X')
ylabel('Y')
Elapsed time is 0.597879 seconds.