辛结构与哈密顿正规形

本文作者：天疆说

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

辛矩阵与辛变换

在标准相空间 R2n\mathbb{R}^{2n}，坐标 z=(q,p)\mathbf{z}=(\mathbf{q},\mathbf{p})，定义块矩阵 J=(0I−I0)J=\begin{pmatrix}0&I\\-I&0\end{pmatrix}。实方阵 MM 满足 M ⁣⊤JM=JM^{\!\top}JM=J 则为辛矩阵，所有这样的矩阵构成辛群 Sp(2n,R)Sp(2n,\mathbb{R})。辛矩阵行列式为 1，特征值在 λ↦λˉ\lambda\mapsto\bar\lambda、λ↦1/λ\lambda\mapsto 1/\lambda 下对称[490]，见该书第 2.4 节。光滑变量替换 z=Z(ζ)\mathbf{z}=\mathbf{Z}(\boldsymbol{\zeta}) 若其 Jacobi 处处辛，则为辛变换（正则变换）。它保持 Hamilton 方程的形式不变。

任何 Hamilton 流的状态转移矩阵都是辛矩阵，这是平动点轨道单值矩阵具有 saddle×center×center 特征值结构的根源，也是 庞加莱截面 保持相体积的原因。辛几何是其微分几何抽象：辛流形携带闭的、非退化的 2-形式 ω\omega。Darboux 定理保证所有辛流形局部都形如 (R2n,ω0)(\mathbb{R}^{2n},\omega_0)。

泊松括号

相空间上光滑函数 F,GF,G 的泊松括号定义为

{F,G}=∇F ⁣⊤J ∇G=∂F∂q ⁣⊤∂G∂p−∂F∂p ⁣⊤∂G∂q.\{F,G\}=\nabla F^{\!\top}J\,\nabla G=\frac{\partial F}{\partial\mathbf{q}^{\!\top}}\frac{\partial G}{\partial\mathbf{p}}-\frac{\partial F}{\partial\mathbf{p}^{\!\top}}\frac{\partial G}{\partial\mathbf{q}}.

它是双线性、反对称且满足 Jacobi 恒等式的，使光滑函数构成 Lie 代数。Hamilton 方程写作 F˙={F,H}\dot F=\{F,H\}。FF 守恒当且仅当 {F,H}=0\{F,H\}=0。nn 个相互对合的独立积分使系统可积。泊松括号在辛变量替换下不变，这是 Meyer & Offin（2017）定理 2.6.3 的代数内容，也是正则变换成为哈密顿力学自然变量的原因。

生成函数与李变换

正则变换可由生成函数编码。在 Hori 1966、Deprit 1969 与 Meyer 的李级数方法中[281][391]，自治 Hamilton 量 GG 生成一个单参数正则变换族 ζ↦Zt(ζ)\boldsymbol{\zeta}\mapsto\mathbf{Z}_t(\boldsymbol{\zeta})，对任意函数 FF 的作用为 F↦F+{G,F}+12{G,{G,F}}+⋯F\mapsto F+\{G,F\}+\tfrac12\{G,\{G,F\}\}+\cdots，即李变换。逐阶选择 GG 消去 HH 中不需要的项即得正规形[343]。Hori/Deprit 方法给出了一个递归的、便于计算机代数实现的形式，在天体力学中广泛使用。

正规形流程与同调方程

考虑平衡点（如共线平动点）附近的 Hamilton 量 H=H2+∑k≥3HkH=H_2+\sum_{k\geq 3}H_k，HkH_k 为 kk 次齐次多项式。正规形寻找一个由生成 Hamilton 量 G=∑k≥3GkG=\sum_{k\geq 3}G_k 的时间一映射给出的正则变换 ζ↦Z(ζ)\boldsymbol{\zeta}\mapsto\mathbf{Z}(\boldsymbol{\zeta})，使新变量下不需要的项逐阶消失[344]。第 kk 阶的决定方程为同调方程

H~k=Hk+{Gk,H2}+(低阶已知项),\widetilde H_k = H_k + \{G_k,H_2\} + (\text{低阶已知项}),

未知量为 GkG_k 与新的系数多项式 H~k\widetilde H_k。HkH_k 的每个单项式求解 GkG_k 时产生形如 k⋅ω\mathbf{k}\cdot\boldsymbol{\omega}（频率点积）的分母。若 k⋅ω\mathbf{k}\cdot\boldsymbol{\omega} 很小但不为零，GkG_k 的系数极大（即小除数问题），形式级数可能发散。共线平动点处的双曲频率 λ\lambda 为非共振项分母提供下界，故有限阶化简是良定义的[244]。

Birkhoff 与 Birkhoff-Gustavson 正规形

Birkhoff 正规形（BNF）消去所有非共振单项式，仅保留与 H2H_2 对合的项。非共振情形下化简后的 Hamilton 量只依赖作用量，使截断系统可积[203]，相关论述见 Meyer & Offin 2017 年第 10 章。Birkhoff-Gustavson 正规形处理共振情形，保留共振单项式，适用于空间 CR3BP（两中心频率接近 1:1 共振，产生 halo 轨道，Hill 情形下更高阶共振在 57 阶附近出现）。

偏正规形、约化与 Moser 定理

偏正规形（PNF）只消去不稳定（双曲）项，保留中心流形动力学。配合中心流形 WcW^c 约化，得到 4 维（或平面 2 维）保守 Hamilton 系统，捕捉 L1,2L_{1,2} 邻域的 Lissajous、halo、拟 halo 轨道族而不必处理双曲方向[415]。正规形约化对庞加莱映射做同样的事：化简辛映射以读出稳定性与分岔结构。

Moser 1958 年的定理给出 Birkhoff 正规形在椭圆平衡点附近收敛的充分条件，形式级数此时描述真实的不变曲线，为平动点邻域的稳定性论断与不变流形参数化提供依据[497]。

应用要点

  • L1/L2L_1/L_2 处的正规形计算到 15–35 阶，可给出 halo 与 Lissajous 轨道的准确半解析近似，作为微分修正的初值。
  • 小除数决定实际可达阶数。低阶共振附近需改用共振正规形。
  • Hori/Deprit 方案的计算机代数实现（符号处理器）是主力工具。只需中心动力学时使用偏正规形。

相关概念

完善页面
最近更新:
Contributors: cislunarspace, ouyangjiahong