2. 轨道力学基础与相对运动描述

交会接近轨迹设计的所有公式，最终都要落到一个具体的问题上：追踪航天器相对于目标航天器，此刻在哪里、以什么速度运动、下一刻会到哪里。要回答这个问题，必须先建立一套自洽的力学模型和坐标约定。本节从二体问题出发，依次给出轨道六根数、三类坐标系以及相对位置速度的坐标转换方法，为后续的接近走廊、V-bar/R-bar 机动和 CW 方程建模提供统一的数学底座。

一、二体问题：一切轨道运动的起点

二体问题假设空间中只有两个质点，彼此仅受万有引力作用，忽略其他天体摄动、大气阻力、太阳光压和推力。设地球质量 M、航天器质量 m，且 m≪M，则航天器在地心惯性系中的运动方程为：

r̈ = -μ · r / |r|³
其中 μ = GM ≈ 3.986004418×10¹⁴ m³/s²（地球引力常数）
r 为地心指向航天器的位置矢量，|r| 为地心距

该方程存在解析解，轨道是一条圆锥曲线。对交会任务而言，追踪航天器与目标航天器通常都处于近圆或小偏心率轨道，因此我们最关心的是圆轨道与近圆轨道情形。二体问题的两个守恒量——比机械能 ε 和比角动量 h——决定了轨道的大小与形状：

ε = v²/2 - μ/r      （v 为惯性速度大小）
h = |r × v|

由 ε 可判断轨道类型：ε<0 为椭圆，ε=0 为抛物线，ε>0 为双曲线。交会任务中追踪器与目标器一般都在椭圆轨道上，且 ε 接近圆轨道值 -μ/(2a)。

二、轨道六根数：把状态矢量翻译成可读参数

任意时刻的位置矢量 r 和速度矢量 v 共 6 个分量，完全确定轨道。但工程上更习惯用 6 个几何与时间参数描述，即轨道六根数：

• 半长轴 a：决定轨道能量与周期，T = 2π√(a³/μ)。
• 偏心率 e：决定轨道形状，e=0 为圆，0<e<1 为椭圆。
• 轨道倾角 i：轨道面与赤道面的夹角，0°≤i≤180°。
• 升交点赤经 Ω：从春分点方向量到升交点的角度，确定轨道面在惯性空间中的方位。
• 近地点幅角 ω：从升交点量到近地点的角度，确定椭圆在轨道面内的朝向。
• 真近点角 f（或平近点角 M、偏近点角 E）：确定航天器在轨道上的瞬时位置。

六根数与状态矢量的互换是交会轨迹设计的基本功。例如由 r、v 求 a 和 e：

a = 1 / (2/r - v²/μ)
e = √(1 + 2εh²/μ²)   （ε、h 同上）

在交会任务中，目标器的六根数通常由地面测控给出，追踪器则需根据任务规划反算所需六根数，再通过变轨实现。两器的 a、e、i、Ω 差异决定了相对运动的长期漂移趋势，这是后续接近策略设计的重要输入。

三、三类坐标系：从惯性到相对

交会轨迹设计涉及三个层次的坐标系，必须明确各自定义与转换关系。

地心惯性系（ECI，记为 O-XYZ）：原点在地心，X 轴指向春分点，Z 轴指向北极，Y 轴与 X、Z 构成右手系。该系不随地球自转，是牛顿力学定律成立的参考系，二体问题方程即在此系中建立。

轨道系（记为 O-xyz）：原点在航天器质心，x 轴沿地心指向航天器的方向（径向），z 轴沿轨道面法向（即 h 方向），y 轴与 x、z 构成右手系，指向速度方向附近。轨道系随航天器沿轨道运动，是描述轨道面内几何关系最自然的坐标系。

相对运动系（记为 T-xyz，又称 LVLH 系）：原点在目标器质心，x 轴沿目标器径向，z 轴沿目标器轨道面法向，y 轴沿目标器速度方向（即 V-bar 方向）。追踪器相对目标器的位置速度均在此系中描述。CW 方程、V-bar 接近、R-bar 接近等概念都建立在该系上。

三者的关系是：ECI 是绝对参考，轨道系是单星局部参考，相对运动系是双星相对参考。设计接近轨迹时，先在相对运动系中规划期望的相对位置速度，再转换到 ECI 中求解追踪器的绝对状态，最后反算其六根数并设计变轨。

四、相对位置速度的坐标转换方法

设目标器在 ECI 中的位置为 r_T、速度为 v_T，追踪器为 r_C、v_C。相对运动系 T-xyz 的三轴单位矢量可由目标器状态构造：

x̂ = r_T / |r_T|                    （径向）
ẑ = (r_T × v_T) / |r_T × v_T|      （轨道法向）
ŷ = ẑ × x̂                          （沿迹方向，右手系）

由这三个单位矢量组成从 ECI 到 LVLH 的旋转矩阵 R：

R = [ x̂ᵀ ; ŷᵀ ; ẑᵀ ]   （3×3，行向量为三轴单位矢量）

则相对位置与相对速度在 LVLH 系中表示为：

ρ_LVLH = R · (r_C - r_T)
ρ̇_LVLH = R · (v_C - v_T) - ω × ρ_LVLH

其中 ω 为目标器轨道角速度矢量，ω = (r_T × v_T)/|r_T|²。第二式中的 -ω×ρ 项不可省略，它体现了 LVLH 系自身在惯性空间中旋转带来的牵连速度。若忽略该项，在近地轨道上会引入约 1.1×10⁻³ rad/s 量级的角速度误差，对应百米级相对距离下约 0.1 m/s 的速度偏差，足以影响接近轨迹的精度评估。

反向转换同样常用：已知期望的 ρ_LVLH 和 ρ̇_LVLH，求追踪器在 ECI 中的绝对状态：

r_C = r_T + Rᵀ · ρ_LVLH
v_C = v_T + Rᵀ · (ρ̇_LVLH + ω × ρ_LVLH)

得到 r_C、v_C 后，即可按第二节的公式反算追踪器的六根数，作为变轨设计的初始条件。

五、一个可复算的数值示例

设目标器在 500 km 圆轨道上，取某时刻 ECI 状态（示例假设值）：

r_T = [6878.137, 0, 0] km
v_T = [0, 7.613, 0] km/s
μ = 3.986004418×10⁵ km³/s²

由定义可得 x̂=[1,0,0]，ẑ=[0,0,1]，ŷ=[0,1,0]，即 R 为单位阵。若追踪器位于目标器后方 200 m、下方 50 m，即 ρ_LVLH=[-0.05, -0.2, 0] km，则 r_C = r_T + ρ = [6878.087, -0.2, 0] km。若追踪器相对速度 ρ̇_LVLH=[0, 0.01, 0] km/s，目标器角速度 ω=√(μ/|r_T|³)≈1.098×10⁻³ rad/s，则 v_C = v_T + ρ̇ + ω×ρ ≈ [0, 7.623, 0] km/s。由此可反算追踪器半长轴 a≈6878.1 km，与目标器接近，说明该相对状态对应一次小幅接近机动。

本节建立的二体模型、六根数、三类坐标系与转换公式，构成了交会轨迹设计的完整数学入口。后续的相对运动方程、接近走廊约束和机动规划，都将在这套约定下展开。