[论文解读] On the stability of periodic orbits for differential systems in $\mathbb{R}^n$

Armengol Gasull, Héctor Giacomini+1 more|ArXiv.org|Oct 4, 2006

Advanced Differential Equations and Dynamical Systems参考文献 2被引用 3

一句话总结

本文提出了一种新颖的方法，通过利用 $n-1$ 个不变超曲面（其横截交集定义了周期轨道）来计算 $ ^n$ 中周期轨道的特征乘子，避免了使用标准的一阶变分方程。该方法将稳定性分析简化为求解一个二阶线性微分方程，在斯捷科夫刚体周期轨道上的应用表明，其计算效率高于传统方法。

ABSTRACT

We consider an autonomous differential system in $\mathbb{R}^n$ with a periodic orbit and we give a new method for computing the characteristic multipliers associated to it. Our method works when the periodic orbit is given by the transversal intersection of $n-1$ codimension one hypersurfaces and is an alternative to the use of the first order variational equations. We apply it to study the stability of the periodic orbits in several examples, including a periodic solution found by Steklov studying the rigid body dynamics.

研究动机与目标

• 开发一种不依赖一阶变分方程的周期轨道特征乘子计算替代方法。
• 解决当标准变分方法在计算上变得繁琐时，周期轨道稳定性分析的挑战。
• 将平面结果（将乘子与沿轨道的余因子积分关联）推广至高维系统，利用 $n-1$ 个不变超曲面。
• 将该方法应用于具体实例，包括刚体动力学中的斯捷科夫周期轨道，以评估其实际效用。

提出的方法

• 该方法假设周期轨道 $\Gamma$ 是由 $n-1$ 个余维数为一的超曲面 $f_i(\mathbf{x}) = 0$，$i=1,\dots,n-1$ 定义的横截交集。
• 引入一个余因子矩阵 $\mathbf{k}(\mathbf{x})$，满足 $D\mathbf{f}(\mathbf{x})\,\mathbf{X}(\mathbf{x}) = \mathbf{k}(\mathbf{x})\,\mathbf{f}(\mathbf{x})$，推广了平面余因子概念。
• 稳定性分析被简化为沿轨道求解一个 $n-1$ 阶线性系统，利用变分方程的基本解矩阵 $\mathbf{v}(t)$。
• 通过利用首次积分和对称性，系统进一步简化为关于 $v_4(t)$ 的二阶线性微分方程（16），显著简化了计算。
• 然后从该简化系统的单变量矩阵中推导出特征乘子。
• 该方法在斯捷科夫周期轨道上得到验证，其中将六阶变分系统简化为二阶方程。

实验结果

研究问题

• RQ1 能否在不求解完整一阶变分方程的情况下计算 $ ^n$ 中周期轨道的特征乘子？
• RQ2 如何利用不变超曲面的结构来简化周期轨道的稳定性分析？
• RQ3 余因子矩阵 $\mathbf{k}(\mathbf{x})$ 在将平面余因子方法推广至高维时起到什么作用？
• RQ4 该方法能否有效应用于复杂系统，如斯捷科夫刚体解？
• RQ5 首次积分与超曲面不变性如何降低变分系统的阶数？

主要发现

• 该方法通过将问题简化为二阶线性微分方程（16），成功计算了 $ ^n$ 中周期轨道的特征乘子，避免了在斯捷科夫情形下求解六阶系统。
• 对于斯捷科夫周期轨道，该方法将变分系统从六阶降低至二阶，显著简化了稳定性分析。
• 特征乘子由简化后二阶方程（16）的单变量矩阵导出，该方程通过函数 $\omega(t)$、$A_0(t)$ 和 $A_{nh}(t)$ 包含了动力学信息。
• 通过使用 $n-1$ 个不变超曲面和矩阵余因子 $\mathbf{k}(\mathbf{x})$，该方法将平面余因子公式推广至 $n$ 维，同时保持了指数积分结构。
• 通过简化系统对斯捷科夫轨道的稳定性进行了分析，可基于方程（16）解的行为判断轨道是否不稳定。
• 首次积分与超曲面不变性的使用使系统阶数降低了三阶，展示了该方法在对称系统中的高效性。