2. 轨道动力学核心量复习

轨道异常识别本质上是一个“比较”问题：把当前观测或定轨结果，与按动力学模型推演出来的“应有状态”作比较，超出容许范围就判为异常。因此必须先厘清四个核心量之间的链条关系——瞬时位置速度、轨道六根数、摄动加速度、预报误差传播，它们共同决定了异常判据的门限从哪里来、为什么这样设。

一、位置速度与轨道六根数：同一状态的两种表达

在惯性系（本课程统一采用 J2000 地心赤道惯性系，位置记 r、速度记 v）中，航天器瞬时状态由 6 个分量完全确定。位置速度矢量是最“原始”的表达，直接对应测量量（测距、测速、测角）和数值积分状态；轨道六根数则是同一状态的几何化表达：半长轴 a、偏心率 e、轨道倾角 i、升交点赤经 Ω、近地点幅角 ω、真近点角 f（或平近点角 M）。

两者可互相转换，且是“一对一”的（除 e=0 或 i=0 的奇异情形）。这一点的工程意义在于：异常既可以用位置速度偏差描述，也可以用根数偏差描述，但两者的敏感度完全不同。例如 a 的微小误差会随时间累积成沿迹方向的位置漂移，而 i、Ω 的误差主要体现为法向偏差。因此判据设计时，往往需要把位置速度残差投影到径向、沿迹、法向（R/T/N）三轴，再分别与根数偏差对应起来看。

下表给出一个近地轨道示例（假设值，仅用于说明量级关系），展示同一状态在两种表达下的对应：

二、摄动加速度：预报误差的物理来源

二体模型下轨道是闭合的，但真实航天器始终受到摄动加速度作用。主要项包括：地球非球形引力（J2 项量级约 1e-3 g，是最大摄动）、大气阻力（低轨约 1e-6~1e-4 g，随高度急剧变化）、三体引力（日月，约 1e-6 g）、太阳光压（约 1e-7 g）。这些加速度虽然量级小，但长时间积分后会造成显著的位置偏差。

关键认识是：摄动加速度不是“误差”，而是模型的一部分。真正导致预报误差的是摄动模型的不完整与参数不准，例如大气密度模型在太阳活动高年误差可达 20%~50%，J2 高阶项截断、光压面质比估计偏差等。因此摄动加速度的建模精度，直接决定了预报误差的增长速率。

三、误差传播：从初始偏差到位置偏差

轨道预报误差传播可用状态转移矩阵 Φ(t, t0) 描述：初始状态偏差 δx0 经 δx(t) = Φ(t, t0) δx0 放大。对近地轨道，沿迹方向误差近似随时间线性增长，位置误差大致满足：

δr_along ≈ 3 · (n · δa) · Δt
其中 n 为平均角速度，δa 为半长轴初始误差，Δt 为预报时长

以 a=7000 km 为例，n≈0.0011 rad/s。若初始 δa=100 m，预报 1 天后沿迹误差约 3×0.0011×100×86400 ≈ 28.5 km。这说明半长轴误差是沿迹预报误差的主导因素，也解释了为什么定轨中 a 的精度要求最高。

误差传播还表现为协方差矩阵的演化：P(t) = Φ P0 Φᵀ + Q(t)，其中 Q 为过程噪声，代表摄动模型误差的累积。工程上常用“误差椭球”或“沿迹误差包络”来可视化预报不确定度。

四、异常判据：把误差传播变成可判断的门限

异常判据的核心逻辑是：观测残差若超出预报误差的统计范围，则判为异常。常用形式为马氏距离判据：

d² = (y - h(x))ᵀ · S⁻¹ · (y - h(x))
S = H P Hᵀ + R
若 d² > χ²(α, m)，则判为异常

其中 y 为观测，h(x) 为预报观测，H 为观测矩阵，P 为预报协方差，R 为观测噪声协方差，m 为观测维数，α 为置信水平。这个判据把“误差传播”与“异常判断”直接连起来：预报误差越大，门限越宽，越不容易误报；但门限过宽又会漏报真实异常。

实际工程中还会设置分级判据：例如沿迹偏差超过 3σ 为黄色预警，超过 6σ 为红色告警；同时结合根数偏差（如 δa 超过 500 m、δi 超过 0.05°）做交叉确认，避免单一指标误判。

五、四者关系总览

下面用一张关系图把本节四个核心量串起来：位置速度与轨道六根数是同一状态的双向表达；摄动加速度作用于状态，驱动其演化；误差传播把初始偏差和模型误差放大为预报不确定度；异常判据则基于预报不确定度设定门限，对观测残差做统计检验，最终输出正常/异常结论。