7-2.m

来自「一本很不错的教程,适合初学者,希望对各位爱好者有一定的帮助.」· M 代码 · 共 24 行

a=12;b=9;T0=2*pi;						            
%T0是轨道周期
T=5*T0;dt=pi/100;t=[0:dt:T]';
f=sqrt(a^2-b^2);									
%地球与另一焦点的距离
th=12.5*pi/180;									
%卫星轨道与x-y平面的倾角
E=exp(-t/20);										
%轨道收缩率
x=E.*(a*cos(t)-f);
y=E.*(b*cos(th)*sin(t));
z=E.*(b*sin(th)*sin(t));
plot3(x,y,z,'g')									    
%画全程轨线
hold on,sphere(20),                                 
%画地球
axis off
title '卫星返回地球示例'
x1=-18*R0;x2=6*R0;y1=-12*R0;y2=12*R0;z1=-6*R0;z2=6*R0;
axis([x1 x2 y1 y2 z1 z2])						   
%确定坐标范围
comet3(x,y,z,0.02),hold off		                   
%画运动轨线