1. 拆题:极地特快到底要解什么
这道题全名叫“极地特快:极轨道卫星的特殊解算”,编号是二体问题章节里的习题4.2。刚拿到时我以为就是把倾角设成90度然后套公式,实际做下来才发现,它真正要考验的是“从经典轨道根数到笛卡尔状态矢量,再通过数值积分做轨道预报”这一整条链路。二体运动方程看起来只有一行,但在极轨道场景下,坐标变换、开普勒方程、守恒量验证这些步骤里任何一个符号弄反,结果都会差得离谱。
先明确一下任务。题目给定一颗近地极轨卫星在初始时刻的轨道根数,要求在理想二体模型下计算它的初始位置速度,然后用数值方法积分一段时间,最后验证轨道能量和角动量是否守恒。我为了让过程可复现,把参数取成一套很常见的低轨遥感卫星量级:
| 轨道根数 | 数值 | 说明 |
|---|---|---|
| 半长轴 a | 6878.137 km | 相当于地面高度约500 km的近地轨道 |
| 偏心率 e | 0.001 | 近圆轨道,但不是完全圆轨道 |
| 轨道倾角 i | 90.0° | 纯极轨道,轨道面经过极区 |
| 升交点赤经 Ω | 0.0° | 初始节点方向取ECI的+X轴 |
| 近地点幅角 ω | 0.0° | 近地点正好落在升交点上 |
| 真近点角 ν | 0.0° | t0时刻卫星刚好在近地点 |
纯