4. 时间系统与坐标系统
在卫星任务仿真中,时间与坐标是两条贯穿始终的“基准线”。任何一条轨道根数、一次姿态机动、一段测控可见性判断,最终都要落到某个具体时刻和某个具体坐标系里。如果时间基准不统一,或者坐标系定义混用,仿真结果会出现米级甚至公里级的偏差,而且这种偏差往往不会报错,只会悄悄地把结论带偏。本节把常用的四种时间系统和四类坐标系梳理清楚,并给出它们之间的转换关系与实现要点。
一、四种时间系统:从原子时到航天工程时
时间系统的差异,本质上是“秒长定义”和“起点定义”的差异。理解这一点,就能理解为什么需要这么多时间系统。
UTC(协调世界时)是我们日常和地面测控最常用的时间。它的秒长来自原子钟,但通过闰秒(leap second)使其时刻与地球自转(UT1)保持接近,差值控制在 0.9 秒以内。UTC 的麻烦之处在于闰秒不可预测,历史上大约每 1~2 年插入一次,因此 UTC 不是均匀时间尺度,不能直接用于高精度轨道递推。
TAI(国际原子时)是连续、均匀的原子时,没有闰秒。TAI 与 UTC 的固定关系为:
TAI = UTC + 闰秒总数(截至2024年为37 s)
例如:UTC 2024-01-01 00:00:00 → TAI 2024-01-01 00:00:37
TT(地球时)是用于描述地球附近天体运动的均匀时间尺度,与 TAI 相差一个常数:
TT = TAI + 32.184 s
TT = UTC + 37 s + 32.184 s = UTC + 69.184 s(2024年)
在 J2000 历元下,TT 时刻为 2000-01-01 12:00:00 TT,对应 UTC 为 2000-01-01 11:58:55.816。这个 69.184 秒的差值在低轨卫星约 7.8 km/s 的速度下,对应约 540 km 的位置偏差,绝不能忽略。
GPS 时(GPST)由 GPS 系统维持,起点为 1980-01-06 00:00:00 UTC,秒长与 TAI 一致,但不插入闰秒。因此:
GPST = TAI - 19 s
GPST = UTC + 闰秒总数 - 19 s(2024年为 UTC + 18 s)
GPS 接收机输出的时间标签通常是 GPST,而星历中的参考时刻也以 GPST 表示。仿真中若把 GPST 当作 UTC 使用,会引入 18 秒误差,对低轨卫星相当于约 140 km 的沿轨偏差。
二、四类坐标系:从惯性系到星体固连系
坐标系的选择取决于要描述的对象和运动特性。惯性系适合描述轨道运动,地固系适合描述地面目标与星下点,LVLH 适合描述相对姿态与相对运动。
J2000(J2000.0 平赤道平春分点惯性系)是最常用的准惯性系。原点在地球质心,X 轴指向 J2000 历元(2000-01-01 12:00:00 TT)的平春分点,Z 轴指向该历元的平北极,Y 轴按右手定则确定。它近似惯性,适合作为轨道积分和星历表达的基准。
ECI(地心惯性系)是一个统称,J2000 是其中一种具体实现。工程中常见的还有 TOD(真赤道真春分点)和 MOD(平赤道平春分点)。ECI 与 J2000 的差别主要来自岁差、章动和极移的处理方式。在一般仿真中,若精度要求不高,可直接用 J2000 代替 ECI。
ECEF(地心地固系)随地球一起旋转,原点在地球质心,Z 轴指向协议地球极(CTP),X 轴指向本初子午线与赤道交点,Y 轴按右手定则。ECEF 用于描述地面站、星下点、地面轨迹和地球引力场模型。ECI 到 ECEF 的转换需要用到地球自转角(GAST):
r_ECEF = Rz(GAST) · r_ECI
其中 Rz(GAST) 为绕 Z 轴旋转 GAST 角度的旋转矩阵
LVLH(局部水平局部垂直系)又称轨道坐标系,原点在卫星质心,Z 轴指向地心,X 轴沿轨道速度方向(或指向飞行方向),Y 轴与轨道面法向一致,构成右手系。LVLH 用于描述对地观测姿态、交会对接相对运动和天线指向。ECI 到 LVLH 的转换由卫星位置和速度矢量决定:
Z_LVLH = -r / |r|
Y_LVLH = -(r × v) / |r × v|
X_LVLH = Y_LVLH × Z_LVLH
三、转换链与实现要点
实际仿真中,时间与坐标的转换通常按以下链路组织:
- 输入 UTC 时刻,先转换为 TAI,再转换为 TT 或 GPST;
- 用 TT 时刻计算 J2000 下的轨道位置和速度;
- 用 GAST 将 J2000 位置旋转到 ECEF,得到星下点和地面可见性;
- 用卫星位置和速度构造 LVLH 基向量,将目标方向投影到 LVLH 下判断姿态。
下面这张图给出了四种时间系统与四类坐标系之间的转换关系,以及它们在仿真链路中的典型位置。
四、工程实现中的常见陷阱
第一,闰秒处理。UTC 到 TAI 的转换需要查闰秒表,不能简单用固定值。仿真中若跨越闰秒时刻,必须使用 IERS 发布的闰秒文件。第二,GAST 计算。GAST 由 GMST 加上赤经章动得到,GMST 又依赖 UT1,而 UT1 与 UTC 的差值(DUT1)需要从 IERS 公报获取。若忽略 DUT1,误差可达 0.9 秒,对应约 400 米的地面位置偏差。第三,LVLH 定义歧义。不同文献中 LVLH 的 X 轴可能指向速度方向,也可能指向地心反方向,使用前必须确认定义。第四,坐标系旋转方向。ECI 到 ECEF 是绕 Z 轴旋转 GAST,方向为从 ECI 到 ECEF 的被动旋转,矩阵形式为 Rz(GAST),不要写成 Rz(-GAST)。
总结来说,时间系统的核心是“秒长是否均匀、起点是否连续”,坐标系统的核心是“原点、轴向和随动关系”。在仿真中,先固定时间基准,再固定坐标基准,最后按“UTC→TAI→TT/GPST→J2000→ECEF→LVLH”的链路组织转换,就能保证整条仿真链路的一致性和可复现性。