1、第一,模拟下图所示的卫星返回地球轨迹。
2、第二,启动MATLAB,新建脚本(Ctrl+N),在脚本编辑区输入如下代码:close a盟敢势袂ll; clear all; c造婷用痃lc; dbstop if errora=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 onsphere(20);axis equal; axis([-15,3,-7,7,-1.5,1.5]);axis offtitle('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个轨道周期。