2. 二体问题与轨道根数回顾
摄动分析的所有结论，最终都要落到一个基准上：如果航天器只受中心天体引力、且天体为质点，它的运动就是二体问题。摄动建模的本质，是把真实受力写成"二体引力 + 摄动加速度"，因此必须先明确二体解的形式、轨道根数的定义，以及状态矢量与根数之间的双向转换关系。本节把这些内容整理成后续误差传播可直接引用的初始条件框架。
一、二体问题的解与六个积分常数
在惯性系中，设中心天体引力常数为 μ，航天器位置矢量为 r、速度矢量为 v，二体运动方程为：
d²r/dt² = -μ · r / |r|³该方程存在六个独立积分常数，正好对应轨道六根数。其中：
- 比角动量矢量 h = r × v，方向垂直于轨道面，模长 h = |r||v|cos(飞行路径角)，它给出轨道倾角 i 和升交点赤经 Ω；
- 偏心率矢量 e = (v × h)/μ − r/|r|，模长即偏心率 e，方向指向近地点，给出近地点幅角 ω；
- 轨道能量 ε = v²/2 − μ/r，由它得到半长轴 a = −μ/(2ε)；
- 真近点角 ν 由位置在轨道面内的角度确定，是唯一的随时间变化的量。
需要强调：a、e、i、Ω、ω 在二体问题下是常数，ν 随时间单调变化。摄动的作用正是让前五个量缓慢漂移，而 ν 的变化率也偏离开普勒值。误差传播分析中，初始时刻的这六个量就是状态初值。
二、轨道六根数的几何含义与约定
轨道六根数分为三组：
| 分组 | 根数 | 几何含义 | 取值范围 |
|---|---|---|---|
| 形状 | a | 半长轴，决定轨道尺度与周期 | a > 0（椭圆） |
| 形状 | e | 偏心率，决定轨道扁圆程度 | 0 ≤ e < 1（椭圆） |
| 定向 | i | 轨道面与赤道面夹角 | 0° ≤ i ≤ 180° |
| 定向 | Ω | 春分点方向到升交点的夹角 | 0° ≤ Ω < 360° |
| 定向 | ω | 升交点到近地点的夹角 | 0° ≤ ω < 360° |
| 位置 | ν（或 M、E） | 近地点到当前航天器的夹角 | 0° ≤ ν < 360° |
这里采用地心赤道惯性系（ECI）作为参考系，X 轴指向春分点，Z 轴指向北极，Y 轴按右手法则确定。所有角度均从参考方向逆时针量取。后续摄动分析中若切换到轨道面坐标系或地固系，必须显式说明，避免符号混乱。
三、状态矢量与轨道根数的双向转换
实际工作中，初始条件常以位置速度矢量给出，而摄动分析更习惯用根数表达。两者转换是本节的核心操作。
由根数求状态矢量的流程为：先由 a、e 求半通径 p = a(1 − e²)；再由 ν 求轨道面内的位置和速度分量：
r = p / (1 + e·cosν)
x_p = r·cosν, y_p = r·sinν
vx_p = -√(μ/p)·sinν
vy_p = √(μ/p)·(e + cosν)然后依次绕 Z 轴转 ω、绕 X 轴转 i、绕 Z 轴转 Ω，得到 ECI 系下的 r 和 v。旋转矩阵顺序不能颠倒，这是最容易出错的地方。
由状态矢量求根数则按以下顺序计算：
- h = r × v，i = arccos(h_z / h)；
- 节点矢量 n = k × h，Ω = arccos(n_x / n)，若 n_y < 0 则 Ω = 360° − Ω；
- e = |e|，a = 1 / (2/r − v²/μ)；
- ω = arccos(n·e / (n·e))，若 e_z < 0 则 ω = 360° − ω；
- ν = arccos(e·r / (e·r))，若 r·v < 0 则 ν = 360° − ν。
当 e ≈ 0（近圆轨道）时，ω 和 ν 的几何意义退化，此时应改用近地点幅角与真近点角之和（纬度幅角 u = ω + ν）来描述；当 i ≈ 0（赤道轨道）时，Ω 退化，应改用真经度。这些奇异情形在摄动分析中会反复出现，必须提前建立处理习惯。
四、从二体基准到摄动初始条件
把二体解作为参考轨道，真实运动可写成：
r(t) = r_2body(t) + δr(t)
v(t) = v_2body(t) + δv(t)其中 δr、δv 是摄动引起的偏差。误差传播方程（如状态转移矩阵方法）需要以初始时刻的 δr₀、δv₀ 为输入，而这些偏差往往来自初始根数误差：
δr₀ = (∂r/∂a)δa + (∂r/∂e)δe + (∂r/∂i)δi
+ (∂r/∂Ω)δΩ + (∂r/∂ω)δω + (∂r/∂ν)δν因此，根数到状态矢量的雅可比矩阵是误差传播的起点。实际工程中，初始根数误差通常来自定轨残差或注入参数精度，例如 a 的误差约几十米、角度量误差约 0.001°～0.01° 量级。把这些误差代入上述线性化关系，即可得到初始位置速度偏差，作为后续摄动积分的初值。
把上述关系固化下来，后续每一节的摄动模型都可以按同一套流程处理：先写出摄动加速度，再投影到轨道根数变化率方程（如拉格朗日行星方程或高斯方程），最后用数值积分或状态转移矩阵把初始误差传播到目标时刻。本节建立的根数—矢量转换与初始条件约定，就是这条链路的第一环。