MATLAB模拟卫星返回地球&三维彗星轨道图comet3

1、第一，模拟下图所示的卫星返回地球轨迹。

2、第二，启动MATLAB，新建脚本（Ctrl+N），在脚本编辑区输入如下代码：

close all; clear all; clc; dbstop if error

a=12; b=9; T0=2*pi; dt=pi/100;

t=[0:dt:5*T0];

f=sqrt(a^2-b^2);

theta=12.5*pi/180;

E=exp(-t/20);

x=E.*(a*cos(t)-f); y=E.*(b*cos(theta)*sin(t)); z=E.*(b*sin(theta)*sin(t));

plot3(x,y,z,'--k','LineWidth',1); hold on

sphere(20); 

axis equal; axis([-15,3,-7,7,-1.5,1.5]);axis off

title('A satellite returns to the Earth')

xlabel('x'); ylabel('y'); zlabel('z')

comet3(x,y,z,0.001); hold off

其中， T0是卫星轨道周期；

f是地球到另一个焦点的距离；

theta是卫星轨道与x,y平面倾角；

E是卫星轨道的收缩率；

plot3和comet3分别模拟卫星的静态和动态轨迹。

3、第三，保存和运行上述脚本，得到如下卫星返回地球的模拟图。

4、第四，关于comet3，可以在命令行输入doc comet3，然后回车，查看帮助文档对其介绍。

5、第五，通过改变t=[0:dt:5*T0]，增加或减少卫星返回地球时的轨道数量，其中T0=2*pi是轨道周期。例如改为t=[0:dt:10*T0]，卫星返回地球时就要经历10个轨道周期。