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

linesegment3:=odeplot(Ns1,[t,KepAbs],0..40,numpoints=1000, labels=[t,Lo]): #绘动量矩模变化图一。 linesegment3; #显示动量矩模变化图一。 linesegment4:=odeplot(Ns2,[x(t),y(t),z(t)],0..40,numpoints=1000): #绘轨迹曲线二。 linesegment4; #显示轨迹曲线二。 linesegment5:=odeplot(Ns2,[KepVec[2],KepVec[3]],0..40, numpoints=1000,labels=[Ly,Lz]): #绘动量矩投影图二。 linesegment5; #显示动量矩投影图二。 linesegment6:=odeplot(Ns1,[t,KepAbs],0..40,numpoints=1000, labels=[t,Lo]): #绘动量矩模变化图二。 linesegment6; #显示动量矩模变化图二。 （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):