ESA PB05 式 Brouwer–Lyddane 平均根数轨道传播器（J2–J6 + 日月）

传播器用 RK4 积分平均根数 (a, e, i, Ω, ω, M)，右端项是对所选深度带谐场的 Gauss 变分率轨道平均——即 Brouwer–Lyddane/PB05 Delaunay 正则化的数值等价，并用闭式一阶 J2 率交叉校验（Ω̇ = −3/2·J2·n·(Re/p)²·cos i，ω̇ = 3/4·J2·n·(Re/p)²·(5cos²i−1)）。日月选项叠加四极阶双重平均第三体项（同时对卫星与摄动轨道平均），能再现 GEO 约 0.75–0.95°/年的经典倾角增长，并携带 Kozai–Lidov 临界互倾角 arccos(√(3/5)) = 39.23°。短周期 J2 偏差由瞬时率积分得到，可加回采样密切根数（平均→密切反变换）。报告还运行同力模型的独立 RK4 笛卡尔数值传播并给出两者偏差、地面轨迹共振分类、冻结轨道（J3 长周期）参数与给定几何的太阳同步倾角解。