轨道预报基础：时间系统与坐标系统
做轨道预报这些年，我最大的体会就是——时间系统和坐标系统是地基，地基不稳，上面算得再漂亮也是白搭。说白了，你连卫星现在几点、在哪儿都不知道，怎么预报它明天在哪儿？
这一章我们就把这两个基础打牢。内容不复杂，但细节很多。我建议你耐心看完，尤其是那些容易踩坑的地方。
时间系统：为什么会有这么多时间？
你可能觉得奇怪，时间不就是时间吗？为什么航天领域要搞出好几种？
嗯，这里有个核心矛盾：地球自转不均匀。你想想看，我们日常用的时间（比如北京时间）是基于地球自转定义的。但地球自转在变慢，而且不是匀速变慢。这就导致了一个问题——如果你用基于地球自转的时间来做轨道预报，误差会越来越大。
所以航天领域引入了几个关键的时间系统：
- 世界时（UT1）：直接基于地球自转。我刚开始做项目时，以为UT1就是格林尼治时间，其实不完全对。UT1是天文观测得到的真实地球自转时间。
- 国际原子时（TAI）：基于原子钟，极其稳定。说白了，原子时不管地球转得快还是慢，它只管自己嘀嗒嘀嗒走。
- 协调世界时（UTC）：这是TAI和UT1的折中方案。UTC用原子时的秒长，但通过闰秒来保持和UT1的差异不超过0.9秒。
- 力学时（TDB/TDT）：用于太阳系动力学计算。做深空探测时，这个时间系统就很重要了。
核心要点：轨道预报中，我们通常使用UTC作为输入输出时间，但内部计算必须转换到统一的连续时间系统（如TAI或TDB）。
时间转换：一个容易出错的环节
我在项目中遇到过好几次，预报结果和实测数据对不上，最后发现是时间转换出了问题。尤其是闰秒，很多人会忽略。
注意：UTC到TAI的转换要考虑闰秒。目前（2024年）UTC比TAI慢37秒。这个值不是固定的，每次闰秒都会变。
举个例子，如果你要计算某颗卫星在UTC时间2024-01-01 00:00:00的位置，实际计算时应该先转成TAI：
TAI = UTC + 37秒（当前闰秒偏移）
TDT = TAI + 32.184秒为什么TDT还要加32.184秒？这是历史原因，记住就行。我个人习惯是把这些偏移量写成一个配置文件，每次更新闰秒时改一个地方就好。
坐标系统：卫星在哪儿？
时间搞定了，接下来是坐标。航天器轨道预报常用的坐标系统有几种，每种都有自己的用途。
地心惯性系（J2000.0）
这是最常用的惯性坐标系。原点在地心，X轴指向春分点，Z轴指向北极。说白了，这个坐标系在空间中是固定的，不随地球自转。
为什么需要惯性系？因为牛顿定律只在惯性系中成立。你算轨道运动，必须用惯性系。
地心地固系（ITRS）
这个坐标系和地球固连，随地球一起转。我们平时说的经纬度、海拔，都是在这个坐标系下定义的。
我记得有一次，同事直接把GPS接收机输出的WGS84坐标当成惯性系坐标来算轨道，结果预报误差大得离谱。嗯，这就是坐标系没搞清楚的典型错误。
轨道坐标系
这个坐标系跟着卫星走。比如LVLH坐标系（Local Vertical Local Horizontal），Z轴指向地心，X轴指向速度方向。做相对运动分析时特别好用。
小技巧：做交会对接或者编队飞行时，用轨道坐标系比用惯性系直观得多。我建议你熟练掌握LVLH和惯性系之间的转换。
坐标转换：绕不开的矩阵运算
不同坐标系之间的转换，本质上就是旋转矩阵的乘法。这里我给出一个常用的转换流程：
J2000.0 → 瞬时真天球 → 瞬时真地球 → ITRS每一步都涉及一些参数：
- 岁差和章动（J2000.0到瞬时真天球）
- 地球自转（瞬时真天球到瞬时真地球）
- 极移（瞬时真地球到ITRS）
你可能会问，这些参数怎么获取？我一般用IERS（国际地球自转服务）发布的公报。网上有现成的数据文件，直接下载就行。
知识体系总览
下面这张图把时间系统和坐标系统的核心逻辑串起来了。你可以看到，时间系统和坐标系统在轨道预报中是交织在一起的。
实际应用中的注意事项
讲完了理论，说说实际中容易出问题的地方。
- 时间精度：轨道预报对时间精度要求很高。1毫秒的时间误差，对于低轨卫星可能造成几米的位置误差。我建议你至少使用微秒级的时间精度。
- 坐标系选择：近地轨道用J2000.0就够了。但如果是高轨或者深空任务，可能需要考虑更精确的坐标系，比如ICRS。
- 转换参数更新：岁差、章动、极移这些参数不是固定不变的。我每年都会更新一次IERS公报数据。
避坑指南：我曾经因为忘记更新闰秒偏移量，导致某次轨道预报误差达到几十公里。排查了整整两天才发现问题。从那以后，我每次做预报前都会先检查时间系统的配置。
好了，时间系统和坐标系统就讲到这里。内容不多，但都是基础中的基础。你把这些搞清楚了，后面学轨道预报算法就会顺畅很多。