机器翻译，未经母语人士审校。以英文版为准。

轨道在GitHub上编辑来源：docs/AGENCY-ORBIT-VALIDATION.md

通过星历拟合进行力模型验证：方法与经独立核对的残差

Kshana的全力模型引擎（src/precise_od.rs）对真实的机构精密轨道产品进行拟合，所得残差如实、可引用，并带有提交哈希标记。这是路线图里程碑P4（“精密天体动力学：高阶重力、SRP（太阳辐射压）、对照机构数据集的验证”）的验证记录。

此处经独立核对的残差意味着什么#

仅当以下各项全部成立时才报告残差（设计文档docs/design/2026-06-09-precise-astrodynamics-design.md）：

  1. 力模型包含在目标精度下重要的每一项摄动：EGM2008（地球重力场模型2008）高阶地球位（截断于所述的阶次，d/o）、固体潮 + 海潮 + 大气潮汐（src/tides.rs，IERS（国际地球自转与参考系统服务）规范2010第6章）、日月第三体、带圆锥阴影的球形SRP + 估计的C_R、阻力（仅限LEO（低地球轨道））、Schwarzschild + Lense–Thirring广义相对论（GR）项（src/forces.rs）。
  2. 估计器是真正的Gauss–Newton批处理最小二乘，带有变分状态转移矩阵（与全弧段有限差分交叉核对，< 1e-6）、1/σ²观测加权，以及n倍标准差的粗差剔除。
  3. 坐标系/时间使用真实的IERS finals2000A地球定向参数（EOP）（UT1（世界时1，地球自转时间）−UTC（协调世界时）、极移；src/eop.rs），经由已验证的IAU（国际天文学联合会）2006/2000A CIO（天球中间零点）链（src/cio.rs）。SP3（标准产品3，精密轨道格式）GPS（全球定位系统）时 → TT（地球时），采用固定的51.184 s偏移（timescales::gps_to_tt）。
  4. 残差以RTN（径向/沿迹/法向）和3-D给出，含与不含经验加速度，并同时给出原始（未拟合）重叠。
  5. 每个数字都可复现（开放数据CI（持续集成）门禁）且可引用（提交哈希 + 数据集参考 + 夹具SHA-256（SHA：安全哈希算法））。

方法，按数据集#

  1. 解析SP3精密轨道；对选定的卫星，取每一个ITRF（国际地球参考框架）位置解。
  2. 将SP3的GPS历元转换为TT；在该历元从finals2000A求出(UT1, xₚ, yₚ)。
  3. 利用这些EOP经CIO链把每个ITRF位置解旋转到GCRS（地心天球参考系）——即惯性位置观测。动力学对地球位的地固旋转使用同一组 EOP，因此观测与力共用同一个坐标系。
  4. 设定历元状态（位置 = 第一个位置解；速度 = 2阶有限差分），并批量拟合[r, v, C_R]（Tier 1），然后再加入9个RTN每转一周（CPR）经验加速度（Tier 2，带先验约束）。
  5. 报告两种层级及原始重叠的拟合后RTN + 3-D RMS（均方根）。

结果#

Galileo中地球轨道——GREEN（< 5 m判据）#

  • 数据集：ESA/ESOC（ESA：欧洲空间局；ESOC：欧洲空间运行中心）最终多GNSS（GNSS：全球导航卫星系统）轨道ESA0MGNFIN，ITRF，5分钟采样，卫星E11（GSAT0101，Galileo在轨验证（IOV），标称MEO（中地球轨道）），2022-01-01。
  • 开放来源（无需登录）：ESA导航办公室镜像，navigation-office.esa.int/products/gnss-products/2190/。EOP：IERS datacenter.iers.org/data/9/finals2000A.all。
  • 验证提交： 66da3ff (tests/agency_galileo.rs).
  • 夹具（SHA-256）：SP3 e7297f4c…d3a24a3；EOP 6b781d36…cb2ed00f（tests/fixtures/agency/NOTICE.md）。
运行 弧段 d/o n_obs 3-D RMS RTN（R、T、N） C_R 备注
CI夹具，Tier 1（力 + C_R） 8 h 12 97 0.133 m 0.105、0.067、0.047 m 1.174 原始重叠78.7 km
CI夹具，Tier 2（+ 经验CPR） 8 h 12 97 0.070 m 0.048、0.044、0.027 m — 残差减半
全弧段派发，Tier 1（修复前†） 24 h 12 289 0.611 m 0.276、0.381、0.390 m 1.244 workflow_dispatch

上述CI夹具的数字在S2大气潮修复（xval/s2-air-tide-fix/）之后于0.36中重新记录：Galileo的Tier 1由0.132 → 0.133 m（C_R 1.1744 → 1.1743），Swarm-A的N分量变化1–2 mm；全弧段派发的数字没有重新运行。† 修复前：全弧段派发的数字记录于0.36的S2大气潮修复之前，尚未在修复后的代码上重新运行（整日产品由派发工作流获取，做出该修复的会话中无法获得）。在CI夹具上，该修复使Galileo的3-D RMS变化不到1 mm，使Swarm-A的分量变化1–2 mm，因此预计会有同一量级的变化；下一次workflow_dispatch运行会将其重新记录。

MEO场在8阶时重力已收敛（8 h的Tier-1结果在d/o 8、10和12时在毫米量级上完全一致），因此CI夹具截断在d/o-12可以忽略；派发作业运行完整的d/o-70。8 h拟合的纯力结果达到13 cm，加入经验层后为7 cm；完整的24 h弧段（SRP/日食的影响更大，动力学跨度更长）为61 cm——均远在5 m判据之内。

Swarm-A LEO——GREEN（< 5 m判据）#

  • 数据集：ESA Swarm Level-2简化动力学精密科学轨道SW_OPER_SP3ACOM_2_（RDOD_AR，GPS导出，ITRF / IGb14，约2 cm，代尔夫特理工大学处理），卫星Swarm-A（SP3标识L47，约430 km LEO），2022-01-01。
  • 开放来源（无需登录）：ESA Swarm分发服务器https://swarm-diss.eo.esa.int/ → Level2daily/Latest_baselines/POD/RD/Sat_A/（依据ESA数据政策开放）。EOP：同一IERS finals2000A 2022-001序列。
  • 验证提交： ceea70a (tests/agency_swarm.rs).
  • 夹具（SHA-256）：SP3 6cd84b78…acb733e（tests/fixtures/agency/NOTICE.md）。

LEO弧段在力模型中增加了大气阻力。由于密度模型是静态的分段指数模型，轨道分两层拟合：动力学层（仅估计历元状态；C_R固定为1，因为在LEO阻力占主导，而SRP在短弧段上难以分离），以及简化动力学层，它加入承载未建模阻力的每转一周经验加速度——这才是业务上有意义的LEO轨道。

运行 弧段 d/o n_obs 3-D RMS RTN（R、T、N） 备注
CI夹具，动力学（C_R=1） 3 h 70 181 2.687 m 0.925、2.522、0.045 m 残差 ≈ 纯沿迹（阻力）
CI夹具，简化动力学（+经验） 3 h 70 181 0.098 m 0.026、0.092、0.025 m 经验项吸收了阻力

动力学拟合达到5 m判据以内，残差几乎完全在沿迹方向——这是约430 km处典型的阻力特征。经验层吸收了该沿迹误差（2.52 → 0.09 m），相对ESA自己约2 cm的轨道，得到约10 cm的简化动力学拟合。整日、全阶的运行是被忽略的swarm_full_arc_dispatch（分发服务器通过其文件浏览器会话提供该产品，因此创始人下载当日的SP3，并让KSHANA_SWARM_SP3指向它）。带空间天气驱动的NRLMSISE-00（NRLMSISE：美国海军研究实验室质谱仪与非相干散射雷达扩展大气模型）是已注明的升级方案，可进一步收紧动力学层。

LRO月球——经核对，高于 5 m判据（如实）（P4 W4b）#

  • 数据集：真实的NASA/JPL（NASA：美国国家航空航天局；JPL：喷气推进实验室）月球勘测轨道器（LRO；NAIF（导航与辅助信息设施）−85）重建轨迹，来自JPL Horizons，ICRF（国际天球参考框架）中的几何月心状态矢量，2022-01-01，高度约98 km，1分钟采样（241个历元，4 h，约2圈）。使用Horizons文本矢量无需SPK/SPICE（SPK：行星星历核；SPICE：航天器、行星、仪器、C矩阵、事件）读取器。
  • 重力：GRAIL（重力恢复与内部结构实验室）GRGM660PRIM场（NASA戈达德太空飞行中心，660阶），截断到d/o 150，并在d/o 100处拟合，在月球固连的主轴坐标系中计算。
  • 开放来源（无需登录）：Horizons API（应用程序编程接口）ssd.jpl.nasa.gov/api/horizons.api；重力场经ICGEM（国际全球地球模型中心）icgem.gfz-potsdam.de。
  • 验证提交： 4fb82bb (tests/agency_lro.rs).
  • 夹具（SHA-256）：LRO 574e3518…d100f0；GRGM 0ff04184…f029977ae（tests/fixtures/agency/NOTICE.md）。

这是月心动力学——一个不同的力模型（src/lunar_od.rs）：月球固连坐标系中的GRGM场（IAU 2015平均地球取向src/lunar_frame.rs与固定的DE421（Development Ephemeris 421）平均地球 → 主轴（ME→PA）偏移复合而成），加上地球（月球轨道最主要的摄动）和太阳作为第三体，经由与地球数据集所用同一个通用精密Gauss–Newton估计器拟合（即precise_od::ForceModel trait）。

运行 弧段 d/o n_obs 3-D RMS RTN（R、T、N） 备注
CI夹具，动力学（仅状态） 4 h 100 241 12.05 m 2.91、10.11、5.88 m 原始重叠39.9 m
CI夹具，简化动力学（+经验1+2/转） 4 h 100 241 6.67 m 3.41、4.58、3.44 m 每转1次和2次的经验项吸收了沿迹误差

0.36更正：在0.36之前，DE421的ME→PA偏移是按错误的方向（且欧拉角顺序也错）施加的，这使解析的PA坐标系与PA相差约208″，而不是与之重合。现在它已按SPICE的定义转录（经与pxform('MOON_ME','MOON_PA')核对，一致到2.6e-16 rad，tests/lunar_me_pa_offset_spice_oracle.rs；docs/design/evidence/fut05b-me-pa-offset-sign/）。上面的数字是修复之后的；在此之前，原始重叠为53.8 m，动力学为12.6 m（3.07、10.19、6.81 m），简化动力学为6.65 m（3.49、4.56、3.35 m）。

这些残差高于 Galileo（0.13 m）和Swarm-A（0.10 m）所达到的5 m判据，并如实报告。

究竟是什么决定了下限——一次DE级别的交叉验证（及一处更正）

我们曾假设，限制因素是解析的月球取向和星历的保真度（IAU天平动级数，与JPL DE（Development Ephemeris）数值积分的MOON_PA相比精度为数十角秒；Montenbruck–Gill地球/太阳星历，与DE/SPICE核相比约0.3°）。为直接检验这一点，被排除在工作区之外的交叉验证crate xval/anise-lunar-od 仅将这两项输入换成DE级的输入——DE440（Development Ephemeris 440）月球主轴取向（moon_pa_de440_200625.bpc）和DE440星历（de440s.bsp），通过ANISE（姿态、导航、仪器、航天器、星历——纯Rust的SPICE核读取器）读取——并重新运行同一个估计器（提交ae96a29，核SHA-256 c1c7fee… / 60cd55a…）：

层级 解析式 DE-grade
原始重叠 53.8 m 41.5 m
动力学 12.6 m 12.0 m
简化动力学（1+2/转） 6.65 m 6.67 m

0.36说明：下文动力学层的大部分改善来自偏移误差，而不是DE级输入。在ME→PA偏移更正之后，解析方案给出原始39.9 m、动力学12.05 m、简化动力学6.67 m，而DE级为41.5、12.0和6.67 m：解析与DE级拟合现在在每一层上都一致。下表的0.30读数由此而来，保持原文所写。

该假设只对了一半，而且头条结论需要更正。DE级取向和星历确实改善了原始重叠（53.8 → 41.5 m）和动力学拟合（12.6 → 12.0 m）——所以解析取向/星历的误差是真实存在的，并确实限制了这些层级。但它们使简化动力学残差基本不变（6.65 → 6.67 m）：1+2每转经验层早已吸收了取向/星历误差，因此业务上约6.6 m的下限由别的因素决定——一种每转基函数吸收不了、DE级坐标系也消除不了的残差。最一致的解释是在4 h的短弧段上存在未建模的LRO非引力动力学（热再辐射、放气）；场的阶数已被排除（d/o 100和150时完全相同），积分器容差也已被排除（atol 1e-6与1e-9时完全相同）。建设性的推论是：Kshana精简、无需核文件的解析月球方案，对简化动力学（业务）月球轨道已经与DE级保真度相当——纯系数的默认方案在这方面没有留下任何余地。要越过5 m，需要对LRO非引力进行建模（针对该航天器的盒翼/热模型以及更长的多弧段拟合），而不是更好的坐标系。

诚实性约定#

  • < 5 m的“green”判据对Galileo MEO（动力学0.13 m）和Swarm-A LEO（动力学2.69 m / 简化动力学0.10 m）已达到。LRO月球按原样公布：动力学12.05 m / 简化动力学6.67 m（1+2每转经验）——高于该判据（0.36，在ME→PA偏移修复之后；此前为12.6 / 6.65 m）。上文的DE级交叉验证检验并更正了关于限制因素的说法：简化动力学的下限不是解析取向/星历（DE440核使其保持在6.67 m不变），因此它不是靠更好的坐标系来消除，而是靠对LRO非引力动力学建模——这是通往米级的、如实更新后的路径。
  • 对于LEO，动力学（仅状态、静态密度）和简化动力学（带经验加速度）两个层级始终分开报告，让读者看到经验项吸收了什么；简化动力学层才是业务轨道。MEO的经验/纯力层级同样如此分开。
  • 每个残差都在上文附有其提交哈希、数据集参考和夹具校验和。
  • 与真实机构真值核对过的数据集：3个中有3个（Galileo MEO、Swarm-A LEO、LRO月球）。达到< 5 m判据：3个中有2个（Galileo ✓，Swarm-A ✓；LRO为6.67 m，如实高于该判据）。