exm051422_1.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 20 行

%exm051422_1.m
G=6.672e-11;ME=5.97e24;vy0=4000;	x0=-4.2e7;t0=0;tf=60*60*24*9;
tspan=[t0,tf];Y0=[x0;0;0;vy0];
[t,YY,Te,Ye,Ie]=ode45('DYDt3',[],[],[],G,ME,tspan,Y0);	%		<3>
X=YY(:,1);Y=YY(:,2);
shg,plot(X,Y,'b','Linewidth',2);hold on
text(0,6e7,'轨道','Color','b')		%产生蓝色文字注释
axis('image');							%保证x、y轴等长刻度，且坐标框恰包容图形
%在三个事件发生点上画标记
plot(Ye(1,1),0.4e7+Ye(1,2),'r^','MarkerSize',10)
plot(Ye(2,1),0.4e7+Ye(2,2),'bv','MarkerSize',10)
plot(Ye(3,1),-0.4e7+Ye(3,2),'b^','MarkerSize',10)
%把轨道的半周期和全周期标在图上
text(0.8*Ye(3,1),-2e7+Ye(3,2),['t3=' sprintf('%6.0f',Te(3))])
text(0.8*Ye(2,1),1.5e7+Ye(2,2),['t2=' sprintf('%6.0f',Te(2))])
%在x-y坐标上画地球
[XE,YE,ZE] = sphere(10);RE=0.64e7;XE=RE*XE;YE=RE*YE;ZE=0*ZE;
mesh(XE,YE,ZE)
text(1e7,1e7,'地球','Color','r'),	hold off		%产生红色文字注释