数值延拓（Numerical Continuation）

本文作者：天疆说

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

定义

数值延拓（numerical continuation）是跟踪含参数非线性方程组[171][423]

F(x,λ)=0,x∈Rn, λ∈R\mathbf{F}(\mathbf{x},\lambda)=\mathbf{0},\qquad \mathbf{x}\in\mathbb{R}^n,\ \lambda\in\mathbb{R}

解曲线的一类数值方法。其基本思想是：在已知一个解 (x0,λ0)(\mathbf{x}_0,\lambda_0) 的前提下，以它为起点，沿解曲线逐步推进到参数 λ\lambda 的相邻取值，每一步以上一步的解作为下一步的初值猜测，从而扫出整条解分支，见 Seydel 2010 和 Allgower、Georg 1990。

在 圆型限制性三体问题（CR3BP） 中，F\mathbf{F} 通常是周期性条件（状态向量返回初值的打靶方程）加上一个约束，x\mathbf{x} 是初始状态与周期等自由变量，λ\lambda 是沿族推进的参数（雅可比常数 CC、振幅 AzA_z、近月点高度等）。延拓是 周期轨道族（Halo、Lyapunov、DRO、NRHO、butterfly 等）系统化计算与分岔分析的标准工具，由庞加莱首次提出并在 20 世纪 60 至 70 年代随数字计算普及而成熟[341][689]。

一般形式

设系统有 nn 个方程、n+1n+1 个未知量（x\mathbf{x} 含 nn 个分量，外加参数 λ\lambda），则解集在 (x,λ)∈Rn+1(\mathbf{x},\lambda)\in\mathbb{R}^{n+1} 空间中通常是一维曲线（分支）。延拓就是沿这条曲线作参数化为弧长 ss 的追踪：

F(x(s),λ(s))=0.\mathbf{F}(\mathbf{x}(s),\lambda(s))=\mathbf{0}.

对上式关于 ss 求导，得切向 (x˙,λ˙)(\dot{\mathbf{x}},\dot{\lambda}) 满足

Fx x˙+Fλ λ˙=0,\mathbf{F}_\mathbf{x}\,\dot{\mathbf{x}}+\mathbf{F}_\lambda\,\dot{\lambda}=\mathbf{0},

即雅可比 Fx\mathbf{F}_\mathbf{x} 在非退化点（满秩）处唯一确定切线方向（差一个归一化）。各类延拓算法的差异主要在如何选步进方向、如何校正、如何处理转向点。

主要变体

1. 自然参数延拓（Natural Parameter Continuation）

最简单的形式：直接以参数 λ\lambda 为推进变量，取序列 λ0<λ1<⋯\lambda_0<\lambda_1<\cdots，对每个 λk\lambda_k 用 微分修正 / 打靶法求 xk\mathbf{x}_k，初值取 xk−1\mathbf{x}_{k-1}。

F(xk,λk)=0,初值 xk(0)=xk−1.\mathbf{F}(\mathbf{x}_k,\lambda_k)=\mathbf{0},\quad \text{初值}\ \mathbf{x}_k^{(0)}=\mathbf{x}_{k-1}.

实现简单、对光滑且单调的分支有效。致命局限：当解曲线对 λ\lambda 出现转向点（fold / turning point，即 dλ/ds=0d\lambda/ds=0）时，λ\lambda 不再单调，方法失效。CR3BP 中周期-能量关系普遍存在转折，因此自然参数延拓只用于快速预扫，正式计算需用伪弧长。

2. 牛顿延拓（Newton Continuation）

自然参数延拓的初值迭代版本：每步用牛顿法解 F(xk,λk)=0\mathbf{F}(\mathbf{x}_k,\lambda_k)=\mathbf{0}，雅可比 Fx\mathbf{F}_\mathbf{x} 在 xk−1\mathbf{x}_{k-1} 处因子化后迭代

xk(j+1)=xk(j)−Fx−1F(xk(j),λk).\mathbf{x}_k^{(j+1)}=\mathbf{x}_k^{(j)}-\mathbf{F}_\mathbf{x}^{-1}\mathbf{F}(\mathbf{x}_k^{(j)},\lambda_k).

所谓牛顿指的是校正环节采用牛顿迭代，本质仍是自然参数延拓。同样无法穿越转向点。文献中有时把它与下文的牛顿同伦（构造同伦函数的特定方式，见 同伦方法）相混，两者并非一回事。

3. 伪弧长延拓（Pseudo-arclength Continuation）

为绕过转向点而设计，由 Keller 在 1977 年系统化，是当前主流算法（AUTO、MATCONT 等持续软件的核心）。思路是把 λ\lambda 也当成未知量，引入弧长约束补回欠定的一维自由度：

{F(x,λ)=0,x˙k−1 T(x−xk−1)+λ˙k−1(λ−λk−1)=Δs.\begin{cases}\mathbf{F}(\mathbf{x},\lambda)=\mathbf{0},\\[2pt] \dot{\mathbf{x}}_{k-1}^{\,T}(\mathbf{x}-\mathbf{x}_{k-1})+\dot{\lambda}_{k-1}(\lambda-\lambda_{k-1})=\Delta s.\end{cases}

第二式是上一步切向 (x˙k−1,λ˙k−1)(\dot{\mathbf{x}}_{k-1},\dot{\lambda}_{k-1}) 与位移向量的内积等于给定弧长步长 Δs\Delta s。这给雅可比增加一行，新雅可比

[FxFλx˙k−1 Tλ˙k−1]\begin{bmatrix}\mathbf{F}_\mathbf{x} & \mathbf{F}_\lambda \\ \dot{\mathbf{x}}_{k-1}^{\,T} & \dot{\lambda}_{k-1}\end{bmatrix}

在标准雅可比 Fx\mathbf{F}_\mathbf{x} 奇异的转向点处仍非奇异，使方法可以平滑地绕过折叠，见 Allgower、Georg 1990 和 Seydel 2010。

预测-校正（predictor-corrector）实现：每步先沿切向作预测 x~=xk−1+Δs x˙k−1\tilde{\mathbf{x}}=\mathbf{x}_{k-1}+\Delta s\,\dot{\mathbf{x}}_{k-1}，λ~=λk−1+Δs λ˙k−1\tilde{\lambda}=\lambda_{k-1}+\Delta s\,\dot{\lambda}_{k-1}，再用牛顿法联立上述两式校正回解曲线。在 CR3BP 周期轨道族的计算中，这是 Halo、Lyapunov、DRO 等族全分支扫描与分岔跟踪的标准做法，见 Doedel 等 2007、Galan-Vioque 等[321] 和 Zhang 2019。

4. 分段线性（Piecewise-Linear / Simplicial）延拓

不依赖导数、对 F\mathbf{F} 仅要求连续的鲁棒实现：在 (x,λ)(\mathbf{x},\lambda) 空间布单纯形剖分，跟踪完全标号单纯形来近似解曲线。优点是无需雅可比、对不可微问题也能工作，缺点是步进精度低、计算量大，远慢于预测-校正法，见 Allgower、Georg 1990 和 Haberkorn 等[370]。在轨道力学中很少使用，只在雅可比不可计算或同伦曲线高度不规则时作为后备。

转向点与分岔

延拓过程中真正有价值的发现通常发生在解曲线的奇点上：

  • 转向点（fold / saddle-node）：λ\lambda 沿曲线方向局部极值，意味着同一 λ\lambda 值对应多个解或解消失。CR3BP 周期-能量关系图上频繁出现，是族存在性边界的标志。

  • 分岔点（branching / bifurcation）：两条或多条解分支相交[582]。Halo 轨道族即从平面 Lyapunov 族经过 pitchfork 分岔产生（后者关于 z→−zz\to -z 对称性破缺）。DRO、NRHO、butterfly 等族之间的连接点也都通过延拓过程中检测雅可比零空间或 Floquet 乘子穿越 +1+1 来识别，见 Galan-Vioque 等 2014 和 Zhang 2019。

延拓扫描出整张分岔图后，再在各分支点处切换到新分支继续延拓，即可得完整的族系图谱。

CR3BP 中的延拓参数

地月空间周期轨道族延拓常用的参数包括：

参数典型族说明
雅可比常数 CC所有族CR3BP 守恒量。沿 CC 延拓直接对应能量层面的扫描
zz 向振幅 AzA_zHalo 族经典 Halo 族参数（Richardson 三阶近似亦以 AzA_z 为参）
xx 向振幅 AxA_xLyapunov 族平面振幅，自然推广至 Lissajous
周期 TTDRO 族DRO 族在很大周期范围内存在
近月点高度 hph_p月球轨道族工程上直观，但族存在性区间窄
初始 y˙0\dot{y}_0单值族通过固定一个分量降低自由度

固定步长 Δλ\Delta\lambda（或 Δs\Delta s）选得过大会导致预测点偏离真实解曲线、微分修正 不收敛。过小则计算量大。常用自适应策略：根据上一步牛顿迭代次数（迭代次数少则放大步长，多则缩小）或解曲线曲率（曲率大处减小步长）调整。

应用要点

  1. 族的系统化扫描。给定一条种子周期轨道（通常由 微分修正 求得），延拓可一次生成整族数千条轨道，避免逐个孤立的初值猜测。这是 Halo/Lyapunov/DRO/NRHO 族图谱的工业级生成方式，见 Zhang 2019。
  2. 模型间的过渡（model continuation）。把参数取为模型保真度（如从 CR3BP 到双圆四体再到星历 N 体），每步对当前模型求解平动点或周期轨道，得到保真度递增的解序列。这一思路由 Ren 等 2012 和 Dei Tos、Topputo[279] 系统化，是把 CR3BP 设计的轨道移植到真实星历环境的标配流程。
  3. 准周期不变环面族的延拓。在参考周期轨道的中心流形附近建立状态网格，施加频闪映射固定、沿环面无漂移、周期匹配等约束修正后，沿族的切向扰动延拓。这是 CR3BP 中准周期 Lissajous / quasi-Halo 族生成的标准方法，见 Capannolo 等[237] 和 Gómez 等 2001。
  4. 发射窗口与鲁棒性分析。把出发时刻、月面停留时长等工程参数作为延拓变量，从最优解出发扫描可行区间，评估窗口宽度与备份能力[36]。

相关概念

参考文献

  • Doedel et al., 2007, AUTO-07P: Continuation and Bifurcation Software for Ordinary Differential Equations, User's Guide, Concordia University, Montreal，主流延拓软件的实现与文档
完善页面
最近更新:
Contributors: cislunarspace, ouyangjiahong