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 计算第三体项后，可在不损失精度的情况下放大步长[631]。
  • 地月空间第三体摄动评估：地球邻域的太阳、月球摄动。月球邻域的地球、太阳摄动，两类问题都满足 r≪ρjr\ll\rho_j，Battin-Giorgi 是默认选择。
  • 星载导航滤波：实时计算需求下，避免相消意味着可用更短字长（如 32 位浮点）获得同等精度，对星载计算机友好。
  • 教学参考：Battin《航天动力学的数学方法》1999 年修订版与 2018 年中译本 §8.4.1 给出完整推导[187]。Sanna 等[572] 在 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 法内部直接调用。

相关概念

参考文献

  • Giorgi, 1964, Una Formulazione Caratteristica del Metodo di Encke in Vista dell'Applicazione Numerica, Scuola di Ingegneria Aerospaziale, Roma
完善页面
最近更新:
Contributors: cislunarspace, ouyangjiahong