2. 轨道机动数学基础
轨道机动的规划与执行，本质上是在一个可预测的动力学框架内，回答三个问题：航天器现在在哪里、要到哪里去、中间需要施加多大的速度增量。要回答这三个问题，必须先建立一套统一的数学描述。本节围绕二体问题、轨道六根数、状态矢量与轨道根数的互转，以及时间系统与坐标系统的选用，梳理机动计算所依赖的基础工具。
一、二体问题：机动计算的动力学底座
在轨道机动规划中，最常用的简化模型是二体问题：只考虑中心天体与航天器之间的引力，忽略其他摄动。其运动方程为
r̈ = -μ · r / |r|³其中 r 为航天器相对中心天体的位置矢量，μ 为中心天体引力常数。以地球为例，μ ≈ 3.986004418×10¹⁴ m³/s²。二体问题给出两个关键结论：轨道是圆锥曲线；航天器在轨道上的机械能守恒。机动计算中，速度增量 Δv 通常通过比较机动前后两条轨道的速度矢量得到，而二体模型正是这一比较的基准。
需要强调的是，二体模型是“基准模型”而非“真实模型”。地球非球形引力、大气阻力、日月三体引力、太阳光压等摄动会使实际轨道偏离二体预测。机动规划时，通常先用二体模型快速求解，再通过摄动修正或数值积分进行精化。因此，二体问题不是可选项，而是所有机动计算的起点。
二、轨道六根数：轨道的几何与时间描述
二体轨道可以用六个独立参数完整描述，称为轨道六根数：
- 半长轴 a：决定轨道能量与周期，是轨道“大小”的度量。
- 偏心率 e：决定轨道形状，e=0 为圆轨道，0<e<1 为椭圆轨道。
- 轨道倾角 i：轨道平面与参考平面（通常为赤道面）的夹角。
- 升交点赤经 Ω：从参考方向（春分点）到升交点的角度，决定轨道平面在空间中的方位。
- 近地点幅角 ω：从升交点到近地点的角度，决定轨道在平面内的朝向。
- 真近点角 f：航天器当前相对近地点的角度，决定航天器在轨道上的瞬时位置。
前五个根数确定轨道几何，第六个根数确定航天器在轨道上的位置。机动规划中，改变轨道几何通常对应改变 a、e、i、Ω、ω 中的若干项，而改变航天器在轨道上的位置则对应改变 f。例如，共面机动主要改变 a 和 e，而改变轨道平面则需要改变 i 和 Ω。
三、状态矢量与轨道根数互转
在数值计算中，航天器状态通常用位置矢量 r 和速度矢量 v 表示，称为状态矢量。状态矢量与轨道根数之间的互转是机动计算的核心操作。
由状态矢量求轨道根数，基本步骤如下：
- 计算轨道能量：
ε = |v|²/2 - μ/|r|，由a = -μ/(2ε)得到半长轴。 - 计算角动量矢量：
h = r × v，其方向确定轨道平面法向。 - 计算偏心率矢量：
e = ((|v|² - μ/|r|) · r - (r·v) · v) / μ，其模为偏心率，方向指向近地点。 - 由
h与参考平面法向的夹角得到倾角 i，由h在参考平面上的投影得到升交点赤经 Ω。 - 由
e与升交点方向的夹角得到近地点幅角 ω，由r与e的夹角得到真近点角 f。
由轨道根数求状态矢量，则先在轨道平面内计算位置和速度，再通过三次旋转变换（Ω、i、ω）转到惯性系。这一过程在机动规划中反复出现：给定初始轨道根数，计算当前状态矢量；施加 Δv 后得到新速度矢量；再由新状态矢量反算新轨道根数，判断是否满足目标。
四、时间系统与坐标系统在机动计算中的选用
机动计算离不开时间与坐标的统一定义。时间系统方面，常用的有：
- UTC：协调世界时，用于地面测控与任务时间标注。
- TAI：国际原子时，与 UTC 相差整数秒的闰秒。
- TDB：质心力学时，用于高精度轨道动力学计算。
- JD/MJD：儒略日与简化儒略日，便于连续时间计算。
在机动规划中，若只做量级估算，UTC 与 TDB 的差异可以忽略；但若涉及长时间预报或高精度定轨，必须统一到 TDB 或指定时间尺度，并明确历元。
坐标系统方面，常用的有：
- J2000 地心惯性系：以 J2000.0 历元的平赤道和平春分点为基准，是轨道动力学计算的标准惯性系。
- 地心固连系：随地球自转，用于地面站位置、测控可见性计算。
- 轨道坐标系：以航天器为原点，用于描述相对运动与姿态。
- 月心/日心惯性系：用于深空机动与借力飞行。
选用原则是：动力学积分与机动 Δv 计算在惯性系中进行；涉及地面站、测控弧段、目标点经纬度时，转换到固连系；涉及交会对接、编队飞行时，使用轨道坐标系描述相对状态。坐标转换必须明确旋转顺序与符号约定，避免因参照系不一致导致 Δv 方向错误。
五、机动计算中的典型工作流
把上述工具串起来，一次典型的轨道机动计算可以这样组织：
- 给定初始时刻
t₀和初始轨道根数，转换为 J2000 惯性系下的状态矢量(r₀, v₀)。 - 根据目标轨道，确定需要改变的根数，例如将半长轴从 7000 km 提升到 8000 km。
- 在二体模型下求解所需速度增量 Δv，得到机动后速度
v₁ = v₀ + Δv。 - 由
(r₀, v₁)反算新轨道根数，检查是否满足目标。 - 若需要高精度，将
(r₀, v₁)送入含摄动的数值积分器，在 TDB 时间尺度下预报，并迭代修正 Δv。
这一流程中，状态矢量与轨道根数的互转是“翻译器”，时间系统是“标尺”，坐标系统是“舞台”。三者缺一不可。实际工程中，Δv 的求解往往不是一次完成，而是通过打靶法或优化算法迭代，但每一次迭代都离不开本节所述的基础转换与统一约定。
掌握这些数学基础后，后续的共面机动、平面改变、交会对接等具体机动方式，都可以归结为在统一的时间与坐标框架下，对状态矢量施加 Δv 并重新解释轨道根数的过程。