前阵子帮人排查一套地面站的轨道预报模块现象挺怪同一份两行根数注进去算出来的过境时刻能差二十多秒跟踪天线偶尔偏出半个波束宽度。查了两天坐标转换没问题校时也没问题最后发现是代码里把两行根数中那个平均运动字段直接当成开普勒根数里的平均角速度套二体公式反算了半长轴再用这条理想轨道做递推。误差就是这么一点点攒出来的。这件事说明一个很基础但极容易被跳过的点卫星六根数和两行根数是两套语言前者是几何描述后者是带摄动修正的拟合结果长得像脾气完全不同。这篇就把这两套东西拆开讲清楚包括每个量的物理含义、相互换算、代码实现以及我在实际项目里踩过的那些坑。不管你是刚接触轨道计算的学生还是要写定轨、预报、覆盖分析模块的工程师下面这些内容都能直接拿去用。1. 六个数字如何把一条轨道钉死1.1 半长轴定大小偏心率定形状开普勒根数的核心思路是把卫星看成在一个固定椭圆上跑这个椭圆由六个独立参数完全确定。其中两个描述形状和大小四个描述姿态和相位。半长轴 a 就是椭圆长轴的一半单位通常用公里。它直接决定轨道周期公式是 T 2π√(a³/μ)μ 是地球引力常数取 398600.4418 km³/s²。拿一条典型低轨遥感轨道举例a 6778.137 km代进去算6778.137³ ≈ 3.1139×10¹¹除以 μ 得 7.812×10⁵开方约 883.9 秒乘 2π 得到周期约 5554 秒也就是 92.6 分钟。这个数字和实际观测的低轨周期完全对得上说明公式没记错。偏心率 e 描述椭圆的扁的程度取值 0 到 1 之间是椭圆等于 0 是正圆大于等于 1 就变成抛物线或双曲线那已经是逃逸轨道了。工程上真正关心的往往是近地点和远地点高度换算很直接近地点半径 rp a(1-e)远地点半径 ra a(1e)各自减去地球平均半径 6378.137 km 就是高度。还是上面那组数e 0.001 时近地点高度约 392.7 km远地点高度约 406.3 km轨道高度起伏只有十几公里。如果换成 e 0.7 的椭圆轨道远地点半径是近地点的 5.67 倍那就是另一回事了这种轨道在转移段和某些科学探测任务里才见得到。这里有个新手常犯的错拿到 a 和 e 之后直接报轨道高度 6778 公里把半长轴当地心距用了。地心距是随位置变化的只有圆轨道上才处处等于 a。判断一个量到底是半径还是高度看它有没有减地球半径这一步错了后面全错。1.2 倾角与升交点赤经轨道平面在惯性系里的朝向确定了椭圆形状还要确定这个椭圆所在平面在空间中的朝向。两个角度就够倾角 i 和升交点赤经 Ω。倾角是轨道平面和赤道面的夹角范围 0 到 180 度。i 小于 90 度是顺行轨道卫星运动方向和地球自转同向大于 90 度是逆行轨道接近 90 度是极轨能覆盖全球。太阳同步轨道是个特殊例子典型倾角在 97 度到 98 度之间属于轻微逆行靠这个倾角配合轨道高度让升交点每天向东进动大约 0.9856 度正好跟上太阳在黄道上的视运动从而保证卫星每次过同一纬度时的地方时大致相同。这个性质对光学遥感极其关键因为成像需要稳定的光照条件。升交点赤经 Ω 是从春分点方向惯性系 X 轴沿赤道面量到升交点的角度范围 0 到 360 度。升交点就是卫星从南半球穿过赤道面进入北半球的交点。直观理解i 决定轨道面倾斜多少Ω 决定这个倾斜面绕赤道转了多少两个角合起来就把轨道面的法向量方向唯一确定了。用球坐标类比i 相当于极角Ω 相当于方位角。要特别注意Ω 不是常数。地球扁率 J2 项会让轨道面缓慢进动低轨卫星的升交点每天能漂好几度。很多人第一次写长期预报直接用初始 Ω 加上时间推下去几天后地面轨迹就完全错位问题就出在这里。这个进动的量级估算方法我放到第 5 章细讲。1.3 近地点幅角与真近点角椭圆上的两个角度标尺平面确定了椭圆在这个平面里还可以转卫星在椭圆上的位置也需要一个角度。这就是 ω 和 ν 的活儿。近地点幅角 ω 从升交点开始沿轨道运动方向量到近地点范围 0 到 360 度。它决定椭圆长轴在轨道面内的朝向换句话说近地点到底在轨道的哪个方位上。真近点角 ν 从近地点量到卫星当前位置也是 0 到 360 度它给出卫星此刻在椭圆上的确切位置。这里有个很好用的组合量纬度幅角 u ω ν表示从升交点量到卫星当前位置的角度。在近圆轨道上u 比单独的 ω 和 ν 稳定得多因为 e 接近 0 时近地点位置几乎是随机的ω 和 ν 会剧烈抖动但它们的和很稳。真近点角随时间不是匀速变化的这是开普勒第二定律的直接结果近地点附近跑得快远地点附近跑得慢。为了让时间推进变成线性运算引入了平近点角 M定义是 M n(t - tp)其中 n √(μ/a³) 是平均角速度tp 是过近地点时刻。M 和 ν 之间要通过偏近点角 E 搭桥关系是开普勒方程 M E - e·sinE这是个超越方程得迭代求解。E 和 ν 的换算用 tan(ν/2) √((1e)/(1-e))·tan(E/2)。这套换算我见过太多人写错尤其是忘了半角公式里的平方根或者把 E 和 ν 搞混。建议写完一定做闭环验证给定六个根数正向算出位置速度再反解回来六个量应该原样复现。1.4 为什么工程上要换一套轨道要素标准的开普勒六根数在数学上很漂亮但在工程实现里有明显短板。第一个短板是近圆轨道的病态问题。e 小于 10⁻⁴ 时近地点方向基本没有物理意义ω 和 ν 的数值会随机跳变虽然它们的和
