2. 参考坐标系与时间系统

在卫星姿态指向与观测几何分析中,任何“指向”都必须回答两个问题:相对于谁、在什么时刻。坐标系定义了空间基准,时间系统定义了时间基准,二者共同构成观测几何计算的“地基”。本节先建立五个常用坐标系的定义与转换关系,再梳理四种常用时间系统及其在观测几何中的角色。

一、五个参考坐标系的定义

地心惯性系(ECI):以地球质心为原点,Z轴指向J2000.0历元平春分点方向的地球自转轴(近似指向北极),X轴指向J2000.0平春分点,Y轴与X、Z构成右手系。ECI不随地球自转,是描述卫星轨道运动、太阳/月球方向、惯性敏感器测量的自然坐标系。轨道六根数、星历(如TLE、精密星历)通常都在ECI下给出。

地固系(ECEF):同样以地球质心为原点,Z轴指向协议地球极(CTP),X轴指向本初子午线与赤道交点,Y轴构成右手系。ECEF随地球自转,地面站坐标、地形、目标点经纬高都在此系下定义。ECI与ECEF之间的转换由地球自转角(GAST)和极移矩阵决定,是观测几何中“卫星—地面”链路计算的关键一步。

轨道系(Orbital Frame):原点在卫星质心,Z轴指向地心(或天底方向),X轴沿轨道速度方向,Y轴沿负轨道法向(即轨道面法向的反方向),构成右手系。轨道系是描述对地观测、对日定向、轨道机动等任务的中间坐标系,卫星对地指向通常先在轨道系中表达,再转到星体系。

星体系(Satellite Body Frame):固连于卫星本体,原点在卫星质心,三轴与卫星结构基准轴一致(如+X为对地面、+Z为对日面等,具体由卫星设计定义)。星体系是姿态确定与控制的直接对象,星敏感器、陀螺、反作用飞轮等都在此系下安装和测量。

本体系(Local Frame):也称当地轨道坐标系或局部水平坐标系,原点在卫星质心,Z轴指向天顶(或地心),X轴指向当地北向,Y轴指向东向。本体系常用于描述卫星相对于当地水平面的姿态,如对地成像时的侧摆角、俯仰角,以及地面站可见性判断中的方位角、俯仰角。

二、坐标系之间的转换链

观测几何分析中,一条典型转换链为:目标点(ECEF)→ 卫星位置(ECI)→ 轨道系 → 星体系 → 本体系。每一步转换都对应一个旋转矩阵或四元数,且必须明确旋转顺序与正负号约定。例如,ECI到ECEF的转换可写为:

r_ECEF = Rz(-GAST) · r_ECI

其中GAST为格林尼治真恒星时。ECEF到轨道系的转换需要卫星位置和速度矢量,先构造轨道系三轴单位矢量,再组装旋转矩阵。轨道系到星体系的转换由卫星姿态确定系统给出,通常以四元数或欧拉角形式输出。星体系到本体系的转换则依赖卫星位置对应的地理经纬度。

转换链中每一步的旋转矩阵必须正交,且行列式为+1。实际实现时建议用四元数或方向余弦矩阵,避免欧拉角奇异点。

三、时间系统的定义与角色

UTC(协调世界时):以原子时秒长为基础,通过闰秒与UT1(地球自转时)保持接近。UTC是地面站、任务规划、遥测遥控的常用时间标签,但因其含闰秒,不适合直接用于高精度轨道计算。

TAI(国际原子时):由全球原子钟加权平均得到,均匀连续,无闰秒。TAI与UTC的关系为:TAI = UTC + 闰秒数(当前为37秒)。TAI是许多航天器时间系统的内部基准。

TT(地球时):原称地球动力学时,是描述太阳系天体运动的均匀时间尺度。TT = TAI + 32.184秒。在精密定轨、星历计算中,TT用于消除引力场时间延迟影响。

GPS时:由GPS系统维持,起点为1980年1月6日0时UTC,无闰秒。GPS时 = TAI - 19秒。GPS接收机输出的时间标签通常为GPS时,需转换为UTC或TAI才能与地面任务时间对齐。

四、时间系统在观测几何中的角色

观测几何计算中,时间系统的作用体现在三个层面:第一,卫星位置和速度是时间的函数,ECI下的星历插值必须使用均匀时间尺度(TAI或TT),否则闰秒会导致位置误差;第二,地球自转角GAST是UT1的函数,而UT1与UTC的差异(DUT1)需通过IERS公报获取,用于ECI到ECEF的精确转换;第三,地面站与卫星的可见性判断、成像时刻的太阳高度角计算,都需要统一到同一时间基准上。

例如,某卫星在UTC 2025-01-01 00:00:00时刻过境,若直接将该UTC标签用于ECI星历插值,忽略闰秒,则卫星位置在沿轨方向可能产生约37秒×7.5 km/s ≈ 277 km的误差,完全不可接受。正确做法是先将UTC转换为TAI或TT,再进行星历插值。

坐标系转换链与时间系统角色 地心惯性系 ECI (J2000) 地固系 ECEF 轨道系 Orbital Frame 星体系 Body Frame GAST r,v 姿态 本体系 Local Frame 经纬度 时间系统及其在观测几何中的角色 UTC 含闰秒,任务标签 TAI 均匀连续,内部基准 TT 星历计算,精密定轨 GPS时 接收机输出,无闰秒 +37s +32.184s -19s 地面站/任务规划 星历插值基准 引力场时间延迟 接收机时间标签 ⚠ UTC直接用于ECI插值:37s × 7.5km/s ≈ 277km沿轨误差

五、实际计算中的注意事项

在编写观测几何计算程序时,建议遵循以下顺序:首先将任务时间标签统一转换为TAI或TT;然后利用IERS提供的DUT1和极移参数,将ECI星历转换到ECEF;接着根据卫星位置和速度构造轨道系;最后结合姿态数据完成到星体系和本体系的转换。每一步转换都应保留中间结果,便于排查误差来源。

对于教学示例,可假设某卫星在UTC 2025-01-01 00:00:00时刻的ECI位置为(7000, 0, 0) km,速度为(0, 7.5, 0) km/s。若忽略闰秒,直接将该UTC标签用于插值,则卫星位置在沿轨方向产生约277 km误差。正确做法是先将UTC转换为TAI(+37秒),再进行插值,此时位置误差可控制在米级以内。

坐标系转换中,旋转矩阵的乘法顺序不可交换。ECI→ECEF→轨道系与ECI→轨道系→ECEF的结果不同,必须严格按照物理意义确定顺序。