[论文解读] Numerical Computation of the Stable and Unstable Manifolds of Invariant Tori

Derin B. Wysham, James D. Meiss|arXiv (Cornell University)|Apr 26, 2005

Quantum chaos and dynamical systems参考文献 26被引用 4

一句话总结

该论文提出了一种基于幂法和广义特征值求解器的迭代数值方法，用于计算具有不可公度旋转矢量的微分同胚中不变环面的稳定与不稳定流形。通过将线性化动力系统重新表述为广义特征值问题，该方法高效地近似了首阶不变子空间，验证了基于Melnikov理论对同宿与异宿分岔的预测，适用于体积守恒映射。

ABSTRACT

We develop an iterative technique for computing the unstable and stable eigenfunctions of the invariant tori of diffeomorphisms. Using the approach of Jorba, the linearized equations are rewritten as a generalized eigenvalue problem. Casting the system in this light allows us to take advantage of the speed of eigenvalue solvers and create an efficient method for finding the first order approximations to the invariant manifolds of the torus. We present a numerical scheme based on the power method that can be used to determine the behavior normal to such tori, and give some examples of the application of the method. We confirm the qualitative conclusions of the Melnikov calculations of Lomelí and Meiss (2003) for a volume-preserving mapping.

研究动机与目标

• 开发一种高效数值算法，用于计算微分同胚中不变环面的稳定与不稳定流形。
• 将不变流形的迭代方法从不动点扩展到更高维、准周期不变环面。
• 为具有不可公度旋转矢量和法向双曲环面的系统提供一种计算上可行的方法。
• 通过Melnikov理论对异宿交叉的分析预测，验证该方法的有效性。

提出的方法

• 该方法采用Jorba提出的方法，将不变环面上的线性化动力系统重新表述为广义特征值问题。
• 应用幂法迭代计算由环面上的上积系统定义的主导特征函数。
• 当系统可约化时，通过坐标变换将准周期矩阵系统简化为常数矩阵形式。
• 通过线性化系统的特征空间计算不变流形的一阶近似。
• 该方法能有效处理一维特征空间，对中间模态采用迭代降维法，但存在数值不稳定性问题。
• 该方法应用于体积守恒映射，结果与Melnikov函数等值线对比以验证。

实验结果

研究问题

• RQ1 如何以高效率计算微分同胚中具有不可公度旋转矢量的不变环面的稳定与不稳定流形？
• RQ2 广义特征值问题在近似准周期上积系统的不变子空间中起什么作用？
• RQ3 如幂法等迭代方法能否被适配以计算可约化准周期系统的特征函数？
• RQ4 数值结果与Melnikov理论对异宿交叉的分析预测相比如何？
• RQ5 当旋转矢量变化时，稳定与不稳定流形交线的同调分岔如何演变？

主要发现

• 基于幂法的迭代方案成功计算了具有不可公度旋转矢量的不变环面的稳定与不稳定流形的一阶近似。
• 该方法验证了Lomelí与Meiss（2003）关于体积守恒映射中异宿交叉的存在性及其同调类型的定性预测。
• 当旋转矢量趋近于共振1/3时，观察到从(1,0)到(3,1)的同调分岔，与Melnikov函数分析一致。
• Melnikov函数的等值线图显示零水平集与观测到的交线曲线吻合，验证了数值结果。
• 该方法在处理一维特征空间时有效，但对多重或复特征值存在数值稳定性挑战，未来工作将解决此问题。
• 使用广义特征值求解器可借助现有高性能线性代数库实现高效计算。