exm07431_2.m

来自「书名应为《掌握和精通MATLAB》附有完整的源程序」· M 代码 · 共 25 行

%exm07431_2.m

shg;
R0=1;													%以地球半径为一个单位
a=12*R0;b=9*R0;
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')									%画全程轨线
[X,Y,Z]=sphere(30);X=R0*X;Y=R0*Y;Z=R0*Z;	%获得单位球坐标
grid on,hold on,
surf(X,Y,Z),shading  interp					%画地球
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])							%确定坐标范围
view([117 37]),										%设视角、画运动轨线
comet3(x,y,z,0.02),
hold off