3. 二体问题与轨道根数

在航天器轨道确定中，二体问题是所有轨道计算的起点。它把地球（或中心天体）视为质点，把航天器视为另一个质点，忽略大气阻力、非球形引力、日月摄动等次要因素，只保留万有引力这一主项。虽然真实轨道并不严格满足二体假设，但二体解提供了轨道根数的定义、位置速度与根数之间的解析转换关系，也是后续数值定轨、摄动修正和状态估计的参考基准。理解二体问题，核心是掌握三件事：开普勒定律描述的几何规律、二体运动方程给出的动力学约束，以及轨道六根数如何唯一确定一条轨道。

3.1 开普勒定律与二体运动方程

开普勒三定律从几何角度概括了二体运动：第一定律指出航天器轨道是以中心天体为焦点的椭圆（圆、抛物线和双曲线是极限情形）；第二定律指出相等时间内矢径扫过的面积相等，等价于角动量守恒；第三定律指出轨道周期的平方与半长轴的立方成正比，即 \(T=2\pi\sqrt{a^3/\mu}\)，其中 \(\mu\) 为中心天体引力常数。对地球，\(\mu\approx3.986\times10^{14}\,\mathrm{m^3/s^2}\)。

从动力学看，二体运动方程为

r̈ = -μ r / |r|³

其中 r 为航天器相对中心天体的位置矢量。该方程是二阶非线性微分方程，但存在解析解。由角动量守恒可得轨道平面固定，由能量守恒可得轨道形状固定。定义机械能 \(\varepsilon=v^2/2-\mu/r\)，当 \(\varepsilon<0\) 时为椭圆轨道，\(\varepsilon=0\) 时为抛物线，\(\varepsilon>0\) 时为双曲线。对椭圆轨道，半长轴由 \(a=-\mu/(2\varepsilon)\) 确定。偏心率 \(e\) 则与角动量 \(h=|\mathbf{r}\times\mathbf{v}|\) 和能量共同决定：\(e=\sqrt{1+2\varepsilon h^2/\mu^2}\)。

3.2 轨道六根数及其几何含义

一条二体轨道在惯性空间中由六个独立参数唯一确定，称为轨道六根数：半长轴 \(a\)、偏心率 \(e\)、轨道倾角 \(i\)、升交点赤经 \(\Omega\)、近地点幅角 \(\omega\)、真近点角 \(f\)（或平近点角 \(M\)、偏近点角 \(E\)）。前两个决定轨道形状和大小，中间两个决定轨道平面在惯性空间中的姿态，第五个决定椭圆在轨道平面内的指向，第六个决定航天器在轨道上的瞬时位置。

以地心惯性系（ECI）为例，参考平面为赤道面，参考方向为春分点。倾角 \(i\) 是轨道面与赤道面的夹角，\(0^\circ\le i\le180^\circ\)；升交点赤经 \(\Omega\) 是从春分点沿赤道到升交点的角度；近地点幅角 \(\omega\) 是从升交点沿轨道到近地点的角度。这三个角度共同完成从惯性系到轨道面的旋转。真近点角 \(f\) 是近地点到航天器当前位置的夹角，随时间变化，是描述航天器在轨道上位置的量。

3.3 位置速度与轨道根数的相互转换

轨道确定中经常需要在位置速度矢量 \(\mathbf{r},\mathbf{v}\) 与轨道六根数之间来回转换。从根数到位置速度的流程较为直接：先由 \(a,e\) 计算半通径 \(p=a(1-e^2)\)，再由真近点角 \(f\) 计算轨道面内的位置和速度分量，最后通过三次旋转（绕 \(z\) 轴转 \(\Omega\)、绕 \(x\) 轴转 \(i\)、绕 \(z\) 轴转 \(\omega\)）得到惯性系下的 \(\mathbf{r},\mathbf{v}\)。轨道面内的位置分量为

r_pf = [p cos f / (1+e cos f), p sin f / (1+e cos f), 0]
v_pf = [-sqrt(μ/p) sin f, sqrt(μ/p)(e+cos f), 0]

从位置速度反求根数则按以下顺序进行：先由 \(h=|\mathbf{r}\times\mathbf{v}|\) 和能量 \(\varepsilon=v^2/2-\mu/r\) 求 \(a\) 和 \(e\)；再由 \(i=\arccos(h_z/h)\) 求倾角；由节点矢量 \(\mathbf{n}=\hat{z}\times\mathbf{h}\) 求 \(\Omega\)；由 \(\mathbf{n}\) 与偏心率矢量 \(\mathbf{e}\) 的夹角求 \(\omega\)；最后由 \(\mathbf{e}\) 与 \(\mathbf{r}\) 的夹角求 \(f\)。这一过程在定轨软件中通常封装为 rv2coe 和 coe2rv 两个函数。

3.4 近圆与近直线轨道的奇异性处理

经典六根数在两类特殊轨道上会出现奇异：当 \(e\approx0\)（近圆轨道）时，近地点方向不确定，\(\omega\) 和 \(f\) 无法独立定义；当 \(i\approx0\) 或 \(i\approx180^\circ\)（近赤道轨道）时，升交点方向不确定，\(\Omega\) 和 \(\omega\) 退化。实际任务中，太阳同步轨道、GEO 轨道、冻结轨道等都可能接近这些奇异条件，直接使用经典根数会导致数值不稳定。

工程上常用的处理方式有两种。第一种是改用非奇异根数，例如用 \(e\cos\omega\)、\(e\sin\omega\) 代替 \(e\) 和 \(\omega\)，用 \(\sin i\cos\Omega\)、\(\sin i\sin\Omega\) 代替 \(i\) 和 \(\Omega\)，这样在近圆和近赤道情况下仍然连续。第二种是采用春分点根数（Equinoctial Elements），包括半长轴 \(a\)、\(h=e\sin(\omega+\Omega)\)、\(k=e\cos(\omega+\Omega)\)、\(p=\tan(i/2)\sin\Omega\)、\(q=\tan(i/2)\cos\Omega\)、平经度 \(\lambda=\Omega+\omega+M\)。这组根数在 \(e=0\) 和 \(i=0\) 处均无奇异，适合长时间轨道递推和状态估计。

下面给出一个近圆轨道示例。设某航天器在 ECI 系中的位置和速度分别为 \(\mathbf{r}=[7000,0,0]\,\mathrm{km}\)、\(\mathbf{v}=[0,7.5,0.5]\,\mathrm{km/s}\)。由 \(h=|\mathbf{r}\times\mathbf{v}|\) 可得 \(h\approx52535\,\mathrm{km^2/s}\)，能量 \(\varepsilon\approx-21.4\,\mathrm{km^2/s^2}\)，进而 \(a\approx9310\,\mathrm{km}\)，\(e\approx0.02\)。由于 \(e\) 很小，若直接计算 \(\omega\) 和 \(f\)，两者会因近地点方向接近不确定而出现较大数值误差。此时可改用 \(e\cos\omega\) 和 \(e\sin\omega\) 作为状态量，或直接采用春分点根数进行后续滤波。该示例中的数值为教学假设，实际使用时需根据真实测量数据重新计算。