2. 时空基准与坐标系统

轨道测量数据处理的第一步，不是解算轨道，而是回答一个更基础的问题：“这个测量值是在哪个坐标系、哪个时刻、以什么时间尺度表达的？”如果时空基准不统一，后续的几何距离、速度、光行时改正都会出现系统性偏差。本节围绕惯性系、地固系、轨道坐标系与时间系统四条主线，说明它们的定义、用途以及在实际归算中的统一与转换方法。

一、四类基准的分工

惯性系（CIS）是牛顿力学成立、无旋转的参考框架。航天器动力学积分、轨道根数、摄动计算都在惯性系中完成。工程上常用 J2000.0 地心惯性系：原点在地心，X 轴指向 J2000.0 时刻的平春分点，Z 轴指向该时刻的平北极，Y 轴按右手系确定。它的关键性质是“不随地球自转”，因此卫星位置矢量在惯性系中变化平滑，适合做动力学建模。

地固系（CTS/ECEF）随地球一起旋转，原点同样在地心，Z 轴指向瞬时自转轴（协议地极 CTP），X 轴指向格林尼治子午面与赤道面的交线。测站坐标、地面天线指向、地球重力场模型、地图投影都建立在地固系上。测站位置在地固系中是常值，这正是它相对惯性系的优势。

轨道坐标系是随航天器运动的局部坐标系，常见形式有：以径向、横向、法向构成的 RSW 系；以切向、法向、径向构成的 TNW 系；以及用于姿态描述的星体坐标系。它把“绝对位置”转化为“相对轨道的方向关系”，便于描述测站可见性、天线指向和相对几何。

时间系统则解决“何时”的问题。TAI 是原子时连续均匀；UTC 通过闰秒与地球自转对齐，是地面测控时钟的常用尺度；TT 与 TAI 相差 32.184 s，用于动力学方程；UT1 反映地球真实自转相位，与 UTC 之差即 DUT1；GPS 时与 TAI 有固定偏差。测量数据归算中，必须把各类时标统一到同一时间尺度，并区分“时刻”与“时间间隔”。

二、坐标转换的实际操作

从地固系到惯性系的转换，标准链路是：CTS → 极移 → 自转 → 岁差 → 章动 → CIS。工程实现中通常写成矩阵连乘：

r_CIS = P · N · R · W · r_CTS
其中：
  W：极移矩阵，由极移量 (xp, yp) 构成
  R：地球自转矩阵，由格林尼治真恒星时 GAST 构成
  N：章动矩阵，由章动角 Δψ、Δε 构成
  P：岁差矩阵，由岁差角 ζ、θ、z 构成

以某测站为例，假设其地固坐标为 (X, Y, Z) = (-2148744.5, 4426641.4, 4044655.9) m（示例数据），在某一时刻计算得到 GAST = 1.234 rad，极移量 xp = 0.05″、yp = 0.30″。先构造 W 矩阵修正极移，再构造 R 矩阵绕 Z 轴旋转 GAST，得到瞬时真赤道地心系下的坐标，最后经岁差、章动矩阵转到 J2000.0 惯性系。每一步的旋转矩阵都是正交矩阵，转换不改变矢量长度，只改变方向表达。

反过来，从惯性系到轨道坐标系，需要先由卫星位置和速度计算轨道面法向：W = (r × v) / |r × v|，径向 R = r / |r|，横向 S = W × R。这样得到的 RSW 系中，测站相对卫星的矢量可以分解为径向、横向、法向三个分量，直接用于判断可见性和计算几何精度因子。

三、时间系统的统一与转换

测量数据中常见的时间标签有 UTC、GPS 时、本地时。归算时统一到 TT 或 TAI 是稳妥做法。转换关系为：

TAI = UTC + 闰秒数（当前为 37 s）
TT  = TAI + 32.184 s
GPS 时 = TAI - 19 s
UT1 = UTC + DUT1（DUT1 由 IERS 公报给出，量级 |DUT1| < 0.9 s）

需要注意两个易错点：其一，闰秒是阶跃量，跨闰秒时刻的数据必须按事件表逐段处理，不能简单用线性公式；其二，UT1 与 UTC 的差异影响地球自转矩阵 R 的精度，对于高精度测距数据，DUT1 的 0.1 s 误差会带来约 46 m 的赤道方向位置偏差，必须使用 IERS 发布的最终或快速公报值。

四、统一基准的检查清单

在实际归算流程中，建议按以下顺序核对基准一致性：

• 时间标签：确认是 UTC、GPS 时还是本地时，是否已扣除闰秒，是否与测站时钟同步。
• 测站坐标：确认是地固系还是惯性系，是否包含极移和板块运动改正。
• 卫星位置：确认是惯性系还是地固系，是否与测站坐标处于同一框架。
• 光行时：信号传播时间约 0.07～0.13 s，需用发射时刻的卫星位置和接收时刻的测站位置迭代求解。
• 方向约定：RSW 与 TNW 的轴定义不同，全文和图示必须保持一致，切换时明确说明。

只有当时、空、坐标系三者都统一到同一基准，测量数据才能进入轨道确定和参数估计环节。这一步做扎实，后续的残差分析才有意义。