2. 二体问题回顾与轨道根数

二体问题是轨道力学的基石。在忽略其他天体摄动、把两个天体视为质点的前提下，牛顿万有引力给出了它们之间唯一的相互作用力。本节从万有引力定律出发，推导二体相对运动方程，并引出描述轨道几何与朝向的开普勒轨道根数，最后说明近焦点坐标系与各轨道要素的几何意义。

从万有引力到相对运动方程

设两质点质量分别为 m₁、m₂，位置矢量为 r₁、r₂，相对位置矢量 r = r₂ − r₁，距离 r = |r|。牛顿万有引力定律给出 m₂ 受到的引力：

F₂₁ = −G·m₁·m₂·r / r³

对两体分别写牛顿第二定律，得到两个耦合方程。将两式相减，并引入引力常数 μ = G(m₁ + m₂)，耦合项被消去，得到只含相对位置 r 的二阶微分方程：

r̈ = −μ·r / r³

这就是二体相对运动方程。它的意义在于：两体的相对运动等价于一个质量可忽略的质点在中心引力场 μ 中的运动。对地球卫星，μ ≈ 3.986×10⁵ km³/s²；对太阳系行星，μ 取太阳的引力常数。求解该方程可得轨道为圆锥曲线，其类型由机械能决定：能量为负对应椭圆（闭合轨道），为零对应抛物线，为正对应双曲线。

开普勒轨道根数

相对运动方程的解是一条平面圆锥曲线，但要在三维空间中唯一确定一条轨道，还需要六个独立参数。工程上最常用的是经典开普勒轨道根数（也称轨道要素），它们把轨道的形状、大小、空间朝向和航天器在轨道上的位置分开描述：

• 半长轴 a：决定轨道大小，椭圆轨道中为长轴的一半。
• 偏心率 e：决定轨道形状，e = 0 为圆，0 < e < 1 为椭圆，e ≥ 1 为抛物/双曲。
• 轨道倾角 i：轨道平面与赤道平面的夹角，范围 0°～180°。
• 升交点赤经 Ω：从春分点方向沿赤道量到升交点的角度，确定轨道平面绕地轴的旋转。
• 近地点幅角 ω：在轨道平面内，从升交点量到近地点的角度，确定椭圆在轨道平面内的朝向。
• 真近点角 ν（或平近点角 M）：从近地点量到航天器当前位置的角度，确定航天器在轨道上的瞬时位置。

前两个根数 (a, e) 描述轨道形状与大小，中间两个 (i, Ω) 描述轨道平面的空间姿态，ω 描述椭圆在平面内的指向，ν 描述航天器当前所在位置。这种“形状—姿态—位置”的分层结构，正是轨道根数便于工程使用的原因。

近焦点坐标系与几何意义

为把上述角度关系表达清楚，引入近焦点坐标系（PQW 系）：原点位于焦点（中心天体），P 轴指向近地点，Q 轴在轨道平面内与 P 轴垂直并指向运动方向，W 轴沿轨道面法向。在该系中，轨道方程可写成极坐标形式：

r = a(1 − e²) / (1 + e·cos ν)

当 ν = 0 时 r 最小，即近地点距离 r_p = a(1 − e)；当 ν = 180° 时 r 最大，即远地点距离 r_a = a(1 + e)。由此可直观看出偏心率对轨道“扁”的程度的影响。近焦点坐标系与地心惯性系之间的转换，由 ω、i、Ω 三个欧拉角依次旋转实现，这也是轨道根数与直角坐标互相转换的标准做法。

一个可动手的换算示例

假设某地球卫星的轨道根数为 a = 7000 km、e = 0.05、i = 45°、Ω = 60°、ω = 30°，当前真近点角 ν = 90°。先算近地点与远地点距离：r_p = 7000×(1−0.05) = 6650 km，r_a = 7000×(1+0.05) = 7350 km。再算当前地心距：

r = 7000×(1−0.05²)/(1+0.05×cos90°)
  = 7000×0.9975/1
  ≈ 6982.5 km

若要把该位置转到地心惯性系，可按 R₃(−Ω)·R₁(−i)·R₃(−ω) 的顺序把近焦点系中的位置矢量旋转到惯性系。实际操作中，先由 ν 求出近焦点系坐标 (r·cosν, r·sinν, 0)，再依次施加三次旋转即可。这类换算是轨道预报、可见性分析和地面站跟踪的起点，建议用一个小脚本把六个根数作为输入，输出惯性系位置速度，便于反复验证。

掌握二体相对运动方程与六个轨道根数后，就具备了把“引力”翻译成“轨道几何”的基本语言。后续的摄动分析、轨道机动与数值积分，都是在这个基准解之上做修正与扩展。