Battin-Giorgi 方法（Battin-Giorgi Method）

本文作者：天疆说

本站地址：https://cislunarspace.cn

定义

Battin-Giorgi 方法是第三体引力摄动加速度的一种数值稳定计算格式。给定中心天体 m1m_1（地球）、摄动天体 mjm_j（月球或太阳）与航天器 r⃗\vec r，最直接的第三体摄动加速度写为

a⃗3B=Gmj ⁣[ρ⃗j−r⃗∣ρ⃗j−r⃗∣3−ρ⃗j∣ρ⃗j∣3]=−Gmj ⁣[d⃗jdj3+ρ⃗jρj3],\vec a_{3B} = G m_j\!\left[\frac{\vec\rho_j-\vec r}{|\vec\rho_j-\vec r|^3}-\frac{\vec\rho_j}{|\vec\rho_j|^3}\right] = -G m_j\!\left[\frac{\vec d_j}{d_j^3}+\frac{\vec\rho_j}{\rho_j^3}\right],

其中 ρ⃗j\vec\rho_j 是摄动天体相对中心天体的位置、d⃗j=r⃗−ρ⃗j\vec d_j=\vec r-\vec\rho_j 是航天器相对摄动天体的位置。当 r≪ρjr\ll\rho_j（航天器紧贴中心天体）时，d⃗j/dj3\vec d_j/d_j^3 与 −ρ⃗j/ρj3-\vec\rho_j/\rho_j^3 是两个量级相近、方向相反的大矢量，直接相减会产生灾难性相消（catastrophic cancellation），有效位数急剧损失。Giorgi (1964) 给出、Battin (1999, §8.4.1) 整理为标准形式的算法，引入辅助变量 qq 把这两个大矢量的差写成单矢量上的小修正：

a⃗3B=−Gmjdj3[r⃗+f(qj) ρ⃗j],\vec a_{3B} = -\frac{G m_j}{d_j^3}\big[\vec r + f(q_j)\,\vec\rho_j\big],

qj=r⃗⋅(r⃗−2ρ⃗j)ρ⃗j⋅ρ⃗j=rρj ⁣(rρj−2cos⁡αj),f(q)=(1+q)3/2−1.q_j = \frac{\vec r\cdot(\vec r-2\vec\rho_j)}{\vec\rho_j\cdot\vec\rho_j} = \frac{r}{\rho_j}\!\left(\frac{r}{\rho_j}-2\cos\alpha_j\right),\qquad f(q) = (1+q)^{3/2}-1.

几何上 q≈−2(r/ρj)cos⁡αj+(r/ρj)2q\approx -2(r/\rho_j)\cos\alpha_j+(r/\rho_j)^2，当 r≪ρjr\ll\rho_j 时 ∣q∣≪1|q|\ll 1，相消问题被转移到对 qq 的良态多项式求值上。

f(q) 的稳定求值

直接按 (1+q)3/2−1(1+q)^{3/2}-1 计算 f(q)f(q) 在 qq 接近零时同样相消。Battin (1999, Eq. 8-60) 给出封闭稳定形式：

f(q)=q 3+3q+q21+(1+q)3/2.f(q) = q\,\frac{3+3q+q^2}{1+(1+q)^{3/2}}.

分子是 qq 的多项式，分母在 q=0q=0 时为 2，整个表达式对任何 qq 都不敏感。也可使用幂级数展开

f(q)=32q ⁣(1+14q−124q2+196q3−⋯ ),f(q) = \tfrac{3}{2}q\!\left(1+\tfrac{1}{4}q-\tfrac{1}{24}q^2+\tfrac{1}{96}q^3-\cdots\right),

在 ∣q∣≪1|q|\ll 1 时收敛迅速，但封闭形式在所有 qq 都给出最佳精度。

与摄动函数展开的关系

第三体摄动函数

Rj=Gmj ⁣(1dj−r⃗⋅ρ⃗jρj3)R_j = G m_j\!\left(\frac{1}{d_j}-\frac{\vec r\cdot\vec\rho_j}{\rho_j^3}\right)

可展开为 x=r/ρjx=r/\rho_j 的幂级数，其中 1/dj=(1+q)−1/2=∑kPk(cos⁡α) xk1/d_j=(1+q)^{-1/2}=\sum_k P_k(\cos\alpha)\,x^k——Legendre 多项式 PkP_k 自然出现。因此 qq 既是相消补救变量，也是经典 Legendre 展开的自变量。Battin 把这一点与 Cowell 摄动法、Encke 法、Hansen 法等摄动数值积分的统一表述联系起来。

应用场景

  • 高精度星历数值积分：地月转移、平动点轨道、月球低轨等场景都需在地球椭球、第三体摄动下长时间积分。Cowell 法直接对总加速度积分时步长受迫变小；改用 Battin-Giorgi 计算第三体项后，可在不损失精度的情况下放大步长。
  • 地月空间第三体摄动评估：地球邻域的太阳、月球摄动；月球邻域的地球、太阳摄动——两类问题都满足 r≪ρjr\ll\rho_j，Battin-Giorgi 是默认选择。
  • 星载导航滤波：实时计算需求下，避免相消意味着可用更短字长（如 32 位浮点）获得同等精度，对星载计算机友好。
  • 教学参考：Battin《航天动力学的数学方法》（1999 / 中译本 2018）§8.4.1 给出完整推导；Sanna et al. (2024) 在 Gateway 到低月轨的最优脉冲转移中援引该方法并把"第三体摄动加速度的稳定计算"列为基础假设之一。

易混点

  • Battin-Vaughan 算法：基于通用变量的 Lambert 问题求解器，与 Battin-Giorgi 完全不同。两者都源自 Battin 的工作，但前者是 Lambert 问题 的根求解，后者是摄动加速度的数值稳定化。
  • q 在 Lambert 问题中：通用变量 Lambert 算法里也定义了一个 q=sin⁡2(ΔE/2)q=\sin^2(\Delta E/2) 类型的辅助变量，与 Battin-Giorgi 的 qq 同名但意义不同。
  • 与 Encke 法的关系：Encke 法在密切轨道偏差方程中直接出现 1/d3−1/ρ31/d^3-1/\rho^3 类型的差，本质上同样面临相消；Battin-Giorgi 提供的是单一稳定的摄动加速度求值，可在 Encke 法内部直接调用。

相关概念

参考文献

  • Battin, R. H. (1999). An Introduction to the Mathematics and Methods of Astrodynamics, Revised Edition, §8.4.1. AIAA.（中译本：《航天动力学的数学方法（修订版）》，2018，第 8 章。）
  • Giorgi, L. (1964). Sur la détermination des positions apparentes des planètes. Astronomische Nachrichten, 339, 250–258.
  • Sanna, D. et al. (2024). Optimal impulsive orbit transfers from Gateway to low lunar orbit. Aerospace, 11(10), 460.
  • Vallado, D. A. (2022). Fundamentals of Astrodynamics and Applications, 5th ed., §9 (numerical treatment of perturbations).
完善页面
最近更新:
Contributors: ouyangjiahong, Ou Yang Jiahong