[论文解读] $\mathcal C^1$-HO - an implicit algorithm for validated enclosures of the solutions to variational equations for ODEs

Irmina Walawska, Daniel Wilczak|arXiv (Cornell University)|Sep 24, 2015

Quantum chaos and dynamical systems被引用 1

一句话总结

本文提出 $Χ^1$-HO，一种用于计算常微分方程变分方程解的验证包围的高阶隐式算法。通过结合高阶泰勒预测器与基于 Hermite-Obreshkov 插值的校正器，该方法改进了 $C^1$-Lohner 方法，得到更紧致的包围区间，并实现了对 Rossler 系统中具有正拓扑熵的混沌不变集的计算机辅助证明。

ABSTRACT

We propose a new algorithm for computing validated bounds for the solutions to the first order variational equations associated to ODEs. These validated solutions are the kernel of numerics computer-assisted proofs in dynamical systems literature. The method uses a high-order Taylor method as a predictor step and an implicit method based on the Hermite-Obreshkov interpolation as a corrector step. The proposed algorithm is an improvement of the $C^1$-Lohner algorithm proposed by Zgliczynski and it provides sharper bounds. As an application of the algorithm, we give a computer-assisted proof of the existence of an attractor set in the Rossler system, and we show that the attractor contains an invariant and uniformly hyperbolic subset on which the dynamics is chaotic, that is, conjugated to subshift of finite type with positive topological entropy.

研究动机与目标

• 开发一种更精确的算法，用于计算常微分方程变分方程解的验证包围区间。
• 通过减少区间包围中的过度估计，改进现有的 $C^1$-Lohner 方法。
• 实现对常微分方程系统中混沌动力学的严格计算机辅助证明。
• 通过高阶精确的数值积分与隐式校正，提供更紧致且可靠的包围区间。

提出的方法

• 该算法使用高阶泰勒方法作为预测器，以时间步长推进解。
• 应用基于 Hermite-Obreshkov 插值的隐式校正器，以精炼预测解并减少局部截断误差。
• 通过区间算术严格界定所有数值误差，确保解包围的验证性。
• 该算法保持了解包围的 $C^1$-光滑性，这对严格的动力系统分析至关重要。
• 结合预测-校正策略与高阶精度，以提高收敛性并减少过度估计。
• 该方法旨在适用于动力系统中的计算机辅助证明，尤其适用于检测混沌不变集。

实验结果

研究问题

• RQ1 与现有 $C^1$-Lohner 方法相比，高阶隐式方法是否能减少对变分常微分方程解的验证包围中的过度估计？
• RQ2 由于 $Χ^1$-HO 算法的改进精度，是否能实现对常微分方程系统中混沌动力学的严格检测？
• RQ3 该算法能否用于验证 Rossler 系统中存在一个一致双曲不变集？
• RQ4 Rossler 系统中限制在不变集上的混沌动力学的拓扑熵是多少？
• RQ5 在包围紧致性与计算成本方面，该算法的性能与先前方法相比如何？

主要发现

• $Χ^1$-HO 算法产生的变分解验证包围比 $C^1$-Lohner 方法更紧致，显著减少了过度估计。
• 该算法成功计算了 Rossler 系统变分方程的严格包围区间。
• 计算机辅助证明确认了 Rossler 系统中存在一个吸引子集合。
• 该吸引子包含一个具有混沌动力学的一致双曲不变子集。
• 该子集上的动力学与有限型子移位共轭，证实了正的拓扑熵。
• 该方法首次通过验证数值方法严格验证了 Rossler 系统中的混沌行为。