3. 卫星轨道基础与位置速度计算

观测几何分析的第一步，是把“卫星在哪里、朝哪里、动得多快”变成可计算的量。没有统一的时空基准，后续的可见性判断、地面覆盖、侧摆角与成像窗口都无从谈起。本节围绕开普勒轨道根数、二体问题、摄动影响、星历表使用与轨道递推展开，目标是让读者能够从一组轨道参数出发，算出任意时刻卫星在惯性系中的位置与速度，并理解这些结果在观测几何中的使用边界。

一、开普勒轨道根数：轨道的六个自由度

在二体假设下，卫星相对地球的运动由六个独立参数完全确定，即开普勒轨道根数。它们不是随意选取的坐标，而是分别描述轨道形状、空间朝向和卫星在轨道上的瞬时位置：

• 半长轴 a：决定轨道能量与周期，是轨道尺度的核心参数，单位通常为 km。
• 偏心率 e：描述轨道形状，e=0 为圆轨道，0<e<1 为椭圆轨道。
• 轨道倾角 i：轨道面与赤道面的夹角，决定卫星能覆盖的纬度范围。
• 升交点赤经 Ω：从春分点方向量到升交点的角度，确定轨道面在惯性空间中的旋转位置。
• 近地点幅角 ω：在轨道面内，从升交点量到近地点的角度，确定椭圆长轴的方向。
• 真近点角 ν：卫星当前相对近地点的角度，是随时间变化的“位置参数”。

前五个根数在短时间内可视为常量，第六个随轨道运动持续变化。因此轨道递推的本质，就是求解 ν 随时间的变化，再把它映射为惯性系中的位置和速度。

二、二体问题：从根数到位置速度

二体问题忽略除地球中心引力外的所有作用力，其运动方程可解析求解。工程上常用平近点角 M 作为时间代理量：

M = n · (t - t0) + M0
n = sqrt(μ / a³)          # 平均角速度，μ 为地球引力常数
E - e·sinE = M            # 开普勒方程，解出偏近点角 E
ν = 2·atan2( sqrt(1+e)·sin(E/2), sqrt(1-e)·cos(E/2) )
r = a·(1 - e·cosE)        # 地心距

其中 μ≈398600.4418 km³/s²。开普勒方程是超越方程，实际计算用牛顿迭代，初值取 E≈M，通常 3～5 次即可收敛到 1e-10 rad 量级。得到 r 与 ν 后，先在轨道面内构造位置和速度分量，再经三次旋转（近地点幅角、轨道倾角、升交点赤经）转到地心惯性系（ECI）。这一步是观测几何的“时空基准”来源：所有地面站、目标点和视线方向，最终都要统一到同一惯性系或地固系中比较。

三、摄动影响：二体之外的修正

真实卫星并不只受中心引力。地球非球形（主要是 J2 项）、大气阻力、太阳光压和日月三体引力都会使轨道偏离二体解。对观测几何影响最直接的是 J2 项，它引起升交点赤经和近地点幅角的长期漂移，量级可达每天零点几度到几度。若忽略这一项，太阳同步轨道的降交点地方时会在数天内明显偏移，成像光照条件随之改变。

工程上通常不直接对二体解做复杂积分，而是采用“平均根数 + 长期项修正”的方式：先按二体公式递推，再叠加 J2 引起的 Ω、ω、M 线性变化率。对于低轨卫星，大气阻力使半长轴缓慢衰减，轨道周期缩短，长期递推必须定期用新星历重置，否则位置误差会快速累积。

四、星历表使用与轨道递推

星历表是轨道计算的“事实来源”。常见格式包括 TLE（两行根数）和 SP3 精密星历。TLE 给出的是平均根数，配合 SGP4 模型使用；SP3 直接给出等间隔的 ECI 或 ECEF 位置，精度更高但数据量更大。使用 TLE 时要注意：

• 历元时刻必须与根数严格对应，不能把不同历元的根数混用。
• TLE 的参考系是 TEME，与 J2000 惯性系存在差异，做精密几何分析时需转换。
• TLE 的有效期通常为几天，超过一周后位置误差可能达到数十公里。

轨道递推的典型流程是：读取星历 → 解析根数 → 计算目标时刻的 M → 解 E → 求 ν、r → 旋转到 ECI → 按需转换到 ECEF 或站心系。下面给出一个可参考的递推参数表，用于教学示例（数值为合理假设，非实测）：

按上述参数，平均角速度 n≈0.001047 rad/s，轨道周期约 100 分钟。递推 30 分钟后，M≈1.885 rad，解出 E 后可得 r≈7000 km 量级，位置在 ECI 中约转过 180°。若叠加 J2 长期项，Ω 在一天内约漂移 1° 左右，这正是太阳同步轨道保持光照条件稳定的关键。

总结来看，开普勒根数提供了轨道的完整描述，二体问题给出了可解析的位置速度计算路径，摄动修正弥补了真实力模型的偏差，星历表与递推则把这一切落到具体时刻。只有把这条链路走通，观测几何中的可见性、覆盖、侧摆和成像窗口才具备可靠的时空基准。