📚 课程目录第2章时间系统与坐标系统

2. 时间系统与坐标系统

轨道确定与状态估计的第一步，不是写动力学方程，而是把“什么时候”和“在哪里”这两个基准统一起来。同一颗卫星的位置矢量，若时间尺度差 1 秒，沿轨方向就会偏差约 7.5 km；若坐标框架选错，误差可能达到数千公里。因此本节先把时间尺度与坐标系统这两套基准讲清楚，再给出它们之间的转换方法。

一、时间尺度：从原子时到质心动力学时

航天器轨道计算中常用的时间尺度可以按“均匀性”和“参考点”两条线索理解。均匀性最好的是原子时，与地球自转、太阳系运动无关；而世界时、恒星时则与地球自转直接绑定，均匀性较差但天文意义明确。

TAI（国际原子时）：由全球数百台铯、铷原子钟加权平均得到，是连续、均匀的时间尺度，无闰秒。它是其他时间尺度的“底层基准”。
UTC（协调世界时）：在 TAI 基础上通过闰秒调整，使其与地球自转（UT1）相差不超过 0.9 s。UTC 是日常与测控系统使用的时间，但不连续，闰秒时刻会出现 23:59:60。
TT（地球时）：原称地球动力学时 TDT，定义在 geocentric（地心）参考系下，与 TAI 保持固定差：TT = TAI + 32.184 s。它用于描述地球附近卫星的动力学方程。
TDB（质心动力学时）：定义在太阳系质心参考系下，与 TT 相差周期性项，量级约 1.6 ms。深空探测、行星际轨道必须使用 TDB。
儒略日 JD / 简化儒略日 MJD：以公元前 4713 年 1 月 1 日 12:00 UT 为起点的连续日计数。MJD = JD − 2400000.5，起点为 1858 年 11 月 17 日 0 时。儒略日把日期与时刻统一为一个实数，便于做时间差和插值。

它们之间的换算关系可以概括为：UTC 与 TAI 之间相差整数闰秒 ΔAT；TAI 与 TT 相差 32.184 s；TT 与 TDB 相差周期项。工程上常用近似：

TT  = TAI + 32.184 s
TAI = UTC + ΔAT   （ΔAT 为当前闰秒数，如 2017 年后为 37 s）
TDB ≈ TT + 0.001657·sin(g)  （g 为地球平近点角，单位 rad）

例如，某测站 2024-01-01 00:00:00 UTC 收到观测数据，此时 ΔAT = 37 s，则 TAI = 00:00:37，TT = 00:01:09.184。若直接拿 UTC 去积分动力学方程，一天内就会积累约 69 s 的系统偏差，对应低轨卫星沿轨约 500 km 的位置误差。

儒略日计算中，2000 年 1 月 1 日 12:00 TT 对应 JD 2451545.0，这是 J2000 历元的定义时刻，也是绝大多数星历和坐标转换的参考起点。

二、坐标系统：从惯性系到地固系

坐标系统解决“在哪里”的问题。按参考点可分为地心系与质心系；按是否随地球自转可分为惯性系与地固系。

J2000（J2000.0 地心惯性系）：原点在地心，X 轴指向 J2000 时刻的平春分点，Z 轴指向 J2000 平北极，Y 轴构成右手系。它是轨道动力学积分最常用的惯性参考系，星历、轨道根数通常都在此系下给出。
ECI（地心惯性系）：广义上指不随地球自转的参考系，J2000 是其中最常用的一种。严格说还有瞬时平赤道系、真赤道系等，区别在于岁差、章动是否改正。
ECEF（地心地固系）：随地球一起自转，Z 轴指向瞬时自转轴，X 轴指向本初子午线与赤道交点。测站坐标、地面轨迹、遥感影像定位都在 ECEF 下描述。
轨道根数（开普勒六根数）：半长轴 a、偏心率 e、倾角 i、升交点赤经 Ω、近地点幅角 ω、真近点角 ν（或平近点角 M）。它把位置和速度转化为几何直观的六个参数，便于判断轨道形状、朝向和卫星当前所处位置。

ECI 与 ECEF 之间的转换核心是地球自转角。设格林尼治视恒星时为 GAST，则位置矢量满足：

r_ECEF = R3(GAST) · r_ECI
R3(θ) = [ cosθ  sinθ  0
         -sinθ  cosθ  0
          0     0     1 ]

速度转换还需加上自转引起的附加项：v_ECEF = R3(GAST)·(v_ECI − ω×r_ECI)，其中 ω 为地球自转角速度矢量。若忽略这一项，低轨卫星速度会有约 0.46 km/s 的误差。

三、转换链路与实操要点

实际轨道确定中，时间与坐标的转换通常按以下链路进行：观测数据以 UTC 记录 → 转 TAI → 转 TT → 转 TDB（深空时）→ 在 J2000 下积分动力学 → 输出 ECI 状态 → 按 GAST 转 ECEF 用于地面比对。每一步都要保留足够位数，儒略日建议用双精度，避免大数相减损失精度。

下面这张图给出时间尺度与坐标系统之间的对应关系和转换方向，可作为工程实现的速查参考。

四、一个可复算的示例

设某卫星在 J2000 下的位置 r_ECI = [6878.137, 0, 0] km，速度 v_ECI = [0, 7.668, 0] km/s，对应约 500 km 高度圆轨道。取 GAST = 30°，地球自转角速度 ω = 7.292115×10⁻⁵ rad/s，则：

R3(30°) = [ 0.8660  0.5000  0
           -0.5000  0.8660  0
            0       0       1 ]

r_ECEF = [ 0.8660×6878.137, -0.5000×6878.137, 0 ]
       = [ 5956.5, -3439.1, 0 ] km

ω × r_ECI = [0, 0, 7.292115e-5] × [6878.137, 0, 0]
          = [0, 0.5016, 0] km/s

v_ECEF = R3(30°) · ([0, 7.668, 0] − [0, 0.5016, 0])
       = R3(30°) · [0, 7.1664, 0]
       = [ 3.5832, 6.2063, 0 ] km/s

可以看到，仅因地球自转，同一时刻的 ECEF 速度与 ECI 速度在数值和方向上都有明显差异。若在定轨中把 ECEF 观测直接当作 ECI 使用，残差会呈现随 GAST 变化的系统性趋势，而不是随机噪声。

常见错误：把 UTC 直接当作 TT 使用；在 ECI 与 ECEF 转换时只转位置不转速度；用单精度存储儒略日导致毫秒级抖动。这些都会在定轨残差中表现为系统性偏差，而不是白噪声。

掌握本节内容后，面对一组观测数据，应能明确回答：它记录在哪个时间尺度、对应哪个坐标框架、需要经过哪些旋转和时标改正才能进入统一的 J2000/TT 计算环境。这是后续建立观测方程和估计器状态量的前提。