Maple理论力学 李银山 第六部分 第24章新.pptVIP

（a）轨迹 (b) 动量矩向量 (c) 动量矩的模随 时间的变化 图24-5 （a）轨迹 (b) 动量矩向量 (c) 动量矩的模随 时间的变化 图24-6 例 24-5(人造地球卫星) 将一个人造地球卫星放在绕地球的轨道上。人造卫星在近地点（最靠近地球的点）的速度为 。近地点的高度为 。求人造卫星的轨道并验证Kepler第二定律的有效性。 已知： ， ， 。 求：轨迹。 解：●建模 此时卫星的轨道是平面的，所以我们在平面上用 坐标解这个问题。 ①利用Newton引力定律，建立人造地球卫星动力学方程； ②绘人造地球卫星轨迹图； ③计算动量矩的误差，验证Kepler第二定律的有效性； ④绘出时程曲线。 ⑤计算周期。我们可以观察到，轨道的周期变化在第326和327时间段之间结束。我们用二等分法确定轨道的周期。 ● Maple程序 restart; #清零。 eq1:=diff(x(t),t)=vx(t): # 方向速度与位移的关系。 eq2:=diff(y(t),t)=vy(t): # 方向速度与位移的关系。 eq3:=diff(vx(t),t)=-G*Mz*x(t)/(x(t)^2+y(t)^2)^(3/2): #人造地球卫星动力学方程在 方向的投影。 eq4:=diff(vy(t),t)=-G*Mz*y(t)/(x(t)^2+y(t)^2)^(3/2): #人造地球卫星动力学方程在 方向的投影。 G:=6.67*10^(-11): Mz:=6*10^24: #引力常数和地球质量。 Inic:=x(0)=7*10^6, vx(0)=0, y(0)=0, vy(0)=9*10^3: #初始条件。 Digits:=15: #计算精度。 var:=x(t),y(t),vx(t),vy(t): #相空间变量。 Ns:=dsolve({eq1,eq2,eq3,eq4,Inic},{var},numeric): #微分方程求解。 MofI(t):=x(t)*vy(t)-y(t)*vx(t): #计算动量矩。 MofI0:=0.63*10^11: #计算初动量矩。 with(plots): #加载绘图库。 p1:=polarplot(6378*10^3,phi=0..2*Pi): #绘地球形状图。 p2:=odeplot(Ns,[x(t),y(t)],t=0..16000,numpoints=400,thickness=2): #绘人造地球卫星轨迹图。 display({p1,p2},labels=[x,y],scaling=constrained); #合并图形。 p3:=odeplot(Ns,[t,(MofI(t)- MofI0)/ MofI0],t=0..16000, numpoints=400): #绘动量矩的相对偏差图。 p3; #显示动量矩的相对偏差。 p4:=odeplot(Ns,[t,x(t)],t=0..16000,numpoints=400): #绘横坐标位移图。 p4; #显示绘横坐标位移图。 p5:=odeplot(Ns,[t,y(t)],t=0..16000,numpoints=400): #绘纵坐标位移图。 p5;