中心流形（Center Manifold）
本文作者：天疆说
定义
中心流形是动力系统理论中描述平衡点附近中心方向运动的不变流形。在圆形限制性三体问题（CR3BP）中，共线平动点（拉格朗日点）的线性化方程特征值包含一对实根（对应鞍点型不稳定/稳定方向）和一对纯虚根（对应中心型振荡方向）。中心流形即由中心方向张成的二维不变流形，包含了平动点附近的全部拟周期轨道和 Lyapunov 轨道。
对于三角平动点（、），其线性化特征值全为纯虚根，整个局部运动空间即为中心流形，因此三角平动点本身即为动力系统理论中的中心点。
与鞍点和中心点的关系
从动力系统理论的角度，两类平动点对应不同的平衡点类型：
| 平动点类型 | 平衡点类型 | 特征值结构 | 存在的流形 |
|---|---|---|---|
| 共线平动点 | 鞍点 | 一对实根 + 一对纯虚根 | 中心流形 + 稳定/不稳定流形 |
| 三角平动点 | 中心点 | 两对共轭纯虚根 | 中心流形 |
共线平动点附近的动力学特征极其丰富：既存在中心流形结构，还存在与之关联的稳定流形和不稳定流形结构。
化简到中心流形的方法
将完整的六维相空间动力学化简到二维中心流形，是研究共线平动点附近运动的关键步骤。该过程涉及以下环节：
1. 线性化与特征值分解
在平动点附近对动力学方程进行线性化，得到 Jacobian 矩阵并计算特征值和特征向量。
2. 非线性坐标变换
通过非线性变换将系统从原始坐标转换到正规形坐标，将动力学方程分离为：
- 中心方向（振荡）部分
- 双曲方向（鞍点）部分
3. 中心流形哈密尔顿量计算
由于限制性三体问题是哈密顿系统，化简到中心流形后得到的约化哈密尔顿量保留了系统的保守性质。中心流形哈密尔顿量的计算涉及高阶项的递推求解。
改进方法
传统方法在高阶项计算中存在数值精度问题。改进方法包括：
- 对中心流形哈密尔顿量的递推计算过程进行优化
- 改进非线性坐标变换的求取算法，提高高阶变换系数的数值精度
在轨道设计中的应用
拟周期轨道计算
中心流形包含了共线平动点附近的全部拟周期轨道。通过在中心流形上选取适当的初始条件，可以计算得到满足任务要求的拟周期轨道。中心流形的降维特性将问题从六维相空间降至二维，大幅简化了轨道计算的复杂度。
基于投影的稳定保持
将动力学中心流形结构引入轨道控制方法设计中，可以得到基于投影到中心流形的稳定保持方法：
- 将航天器的实时状态投影到中心流形上
- 在二维中心流形上计算最优控制量
- 将控制量映射回全空间执行
这种方法充分利用了中心流形的降维特性，可在显著降低燃料消耗的基础上达到良好的稳定保持效果。
相关概念
参考文献
- 李言俊, 张科, 吕梅柏, 张汉清. 利用拉格朗日点的深空探测技术[M]. 西北工业大学出版社.
- Gómez G, Llibre J, Martínez R, et al. Dynamics and Mission Design Near Libration Points, Vol. I: Fundamentals[M]. World Scientific, 2001.
- Jorba A. A methodology for the numerical computation of normal forms, centre manifolds and first integrals of Hamiltonian systems[J]. Experimental Mathematics, 1999, 8(2): 155-195.
- 郑越. 地月三体系统低能转移轨道设计与控制[D]. 西北工业大学, 2019.