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

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

close all; clear all; clc; dbstop if error

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

t=;

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();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分别模拟卫星的静态和动态轨迹。

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

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

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

doc comet3查看帮助文档对comet3的介绍。

comet3输入擦除模式，例如，改变Figure绘图窗口的大小，comet3绘制的图形就会消失。