2. 二体问题与开普勒轨道

二体问题是轨道动力学中最基础、也最重要的可解析模型。它假设空间中只有两个质点，彼此仅通过万有引力相互作用，忽略其他天体摄动、大气阻力、太阳光压与推力。虽然真实航天任务几乎不可能完全满足这一假设，但二体模型给出的开普勒轨道是理解轨道形状、周期、能量与后续摄动分析的基准。本节从惯性系下的运动方程出发，建立开普勒六根数与轨道几何的对应关系，并给出用数值积分验证解析结果的最小实现路径。

惯性系下的运动方程

取地心惯性系（ECI）作为参考系，原点位于地球质心，坐标轴在惯性空间中保持指向稳定。设航天器位置矢量为 r，速度矢量为 v，地球引力常数为 μ = GM。二体问题的运动方程为：

dr/dt = v
dv/dt = -μ · r / |r|³

这里 r 的方向由航天器指向远离地心，因此引力加速度方向与 r 相反，负号不可省略。该方程是二阶非线性常微分方程，解析解对应圆锥曲线轨道。数值积分时通常把它写成一阶状态方程，状态量为六维向量 [x, y, z, vx, vy, vz]，便于用 RK4 或 ode45 求解。

需要特别说明参照约定：本文所有位置、速度均在 ECI 下描述，角度以地心为顶点、以赤道面为基准面。若切换到轨道平面内的近焦点坐标系，则需明确说明，避免方向与正负号混淆。

开普勒六根数

二体轨道由六个独立常数完全确定，称为开普勒六根数。它们分别描述轨道的大小、形状、空间朝向和航天器在轨道上的位置：

• 半长轴 a：决定轨道能量与周期，椭圆轨道 a > 0。
• 偏心率 e：描述轨道形状，0 ≤ e < 1 为椭圆，e = 0 为圆。
• 轨道倾角 i：轨道平面与赤道面的夹角，0° ≤ i ≤ 180°。
• 升交点赤经 Ω：从春分点方向到升交点方向的夹角。
• 近地点幅角 ω：从升交点方向到近地点方向的夹角，在轨道平面内度量。
• 真近点角 ν：从近地点方向到航天器当前位置的夹角，随时间变化。

前五个根数确定轨道几何与空间姿态，第六个根数确定航天器在轨道上的瞬时位置。若已知某时刻的 r 与 v，可通过矢量运算反解出六根数；反之，由六根数也可生成任意时刻的位置速度。

轨道几何与近地点/远地点

椭圆轨道的几何关系由半长轴 a 与偏心率 e 决定。近地点半径与远地点半径分别为：

r_p = a(1 - e)
r_a = a(1 + e)

近地点是航天器离地心最近的位置，速度最大；远地点是离地心最远的位置，速度最小。由活力公式可得：

v_p = sqrt(μ(1 + e) / (a(1 - e)))
v_a = sqrt(μ(1 - e) / (a(1 + e)))

轨道半通径 p = a(1 - e²)，它决定轨道在焦点处的“宽度”。这些量共同刻画了轨道的几何形态，也是判断轨道是否满足任务约束（如覆盖、回归、寿命）的直接依据。

周期与能量

开普勒第三定律给出椭圆轨道周期：

T = 2π · sqrt(a³ / μ)

周期只与半长轴有关，与偏心率无关。这意味着只要半长轴相同，无论轨道是圆还是扁椭圆，周期都相同。比机械能由半长轴唯一决定：

ε = -μ / (2a)

椭圆轨道 ε < 0，抛物线 ε = 0，双曲线 ε > 0。能量为负是轨道被地球引力束缚的标志。数值积分时，可以用能量守恒作为积分精度的一个检查量：若积分结果中比机械能随时间明显漂移，说明步长过大或积分器精度不足。

最简轨道模型与数值积分验证

下面给出一个可直接运行的最小验证流程。以近地圆轨道为例，取 μ = 398600.4418 km³/s²，轨道高度 500 km，则 r = 6878.137 km，圆轨道速度 v = sqrt(μ/r) ≈ 7.613 km/s。初始状态设为：

r0 = [6878.137, 0, 0] km
v0 = [0, 7.613, 0] km/s

用 RK4 以 10 s 步长积分一个周期 T ≈ 5676.9 s，得到末端位置速度。验证要点有三：其一，末端位置与初始位置的距离应远小于轨道半径，理想情况下回到起点；其二，比机械能 ε = v²/2 - μ/r 应保持在 -μ/(2a) 附近，相对漂移量可控制在 1e-8 量级；其三，由末端状态反解六根数，半长轴与初始设定一致，偏心率接近 0。

通过上述流程，二体模型从解析公式到数值实现形成闭环：运动方程给出状态演化规律，六根数给出轨道几何与瞬时位置，周期与能量给出全局约束，数值积分则验证解析关系是否在离散计算中保持。掌握这一最小模型后，后续加入 J2 摄动、大气阻力或三体引力时，只需在加速度项上叠加相应摄动，而轨道根数的变化与能量漂移仍可作为判断摄动影响的重要观察量。