轨道参数概述

概述

Python 是模拟简单二体轨道（如卫星绕地球运行）的理想选择。（它同样适用于复杂的三体轨道，但本课程不会涉及三体轨道。）你需要能够建模二维和三维中的基本圆形和椭圆轨道。你将利用 Matplotlib 的 FuncAnimation() 函数来动态展示二维轨道，显示卫星在靠近地球时速度增加、远离地球时速度减小的现象。此外，还可以使用 Matplotlib 的三维绘图功能绘制三维轨道，以更直观地展示不同轨道参数的影响。

求解二体问题需要某种常微分方程的数值方法来计算时间和位置。这些方法通常依赖于数值算法（如牛顿法、龙格-库塔法、开普勒法等），用于数值求解位置随时间的变化。你可以自行编写微分方程求解器，但也可以直接使用具备该功能的现有库。本课程中，你将使用 PyAstronomy 库中的 KeplerEllipse 类来生成卫星在一段时间内的 x、y 和 z 坐标，然后用 Matplotlib 绘制这些点。PyAstronomy 的文档可在本章末尾的脚注 [1] 中找到。

轨道力学基础

轨道力学的具体原理超出了本课程范围；不过，这里将为你提供一个快速入门。一颗卫星围绕某天体（如地球）的轨道可以通过以下参数进行描述和计算：

• 半长轴与半短轴 (a)：分别表示从椭圆中心到椭圆的最大和最小距离。通常仅使用半长轴作为参数。对于圆形轨道，半长轴与半短轴相等，此时称为轨道半径。

• 偏心率 (e)：描述轨道形状。偏心率为 0 表示轨道为完美圆形；介于 0 和 1 之间的值表示椭圆轨道。例如，偏心率为 0.50 的椭圆轨道意味着半长轴是中心到焦点距离的两倍。偏心率为 1 时表示抛物线轨道，即卫星仅绕地球一次后便不再返回。偏心率大于 1 则为双曲线轨道。

• 轨道倾角 (i)：轨道平面与地球赤道平面之间的夹角。倾角为 0° 的轨道位于赤道上方；倾角为 90° 的轨道经过南北极。通常，轨道倾角与发射地点的纬度一致：从佛罗里达州卡纳维拉尔角（北纬30度）发射的卫星，其轨道倾角约为30度。

• 近心点与远心点：近心点是卫星轨道上距离被绕行天体最近的点，远心点则是最远的点。

  • “-apsis” 是通用后缀，适用于任意天体，因此每条轨道都有近心点；
  • “-gee” 是地球专用后缀，所有绕地球运行的卫星都具有近地点（perigee）；
  • “-helion” 是太阳专用后缀，所有绕太阳运行的卫星都具有近日点（perihelion）。

• 近心点幅角 (w)：轨道上与赤道延长线交点和轨道上离地球最近点（近地点）之间的夹角。对于倾角为 0° 的圆形轨道，近心点幅角为 0°。该参数本质上是将轨道上离地球最近的点沿轨道移动到其他位置。

• 升交点赤经 (omega)：从某一参考方向（如太阳方向）到升交点的夹角，升交点是指轨道“上升”部分与赤道延长线的交点。该参数可用于“旋转”轨道相对于地球的方向。

不必担心完全理解所有参数。可视化正是帮助你掌握这些概念的关键。在了解了轨道力学的基本知识后，我们进入 Python 实践环节。你将首先绘制一个卫星绕地球运行的二维静态图像，接着实现该轨道的动态动画，然后创建同一轨道的三维图像，最后实现三维轨道的动态播放。