exm041421_1.m

来自「这本书是matlab经典书籍」· M 代码 · 共 14 行

global G ME						%在主程序中定义全局变量传递参数		<1>
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];					%按式(4.14.2.1-6)给定初值向量
[t,YY]=ode45('DYDt',tspan,Y0);										%<8>
X=YY(:,1);						%输出Y的第一列是位移数据x(t)
Y=YY(:,2);						%输出Y的第二列是位移数据y(t)
plot(X,Y,'b','Linewidth',2); hold on
axis('image')						%保证x、y轴等长刻度，且坐标框恰包容图形
[XE,YE,ZE] = sphere(10);		%产生单位球面数据
RE=0.64e7;						%地球半径
XE=RE*XE;YE=RE*YE;ZE=0*ZE;		%坐标纸上的地球平面数据
mesh(XE,YE,ZE),hold off			%绘地球示意图