[論文レビュー] $\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階変分方程式の解の検証付き包摂を計算するための高次陰解法 $Χ^1$-HO を提案する。高次テイラー予測子とヘルミート＝オブレシュコフに基づく補正子を組み合わせることで、$C^1$-Lohner法を改善し、より鋭い境界を得ることができ、ローラン系におけるカオス的不変集合のコンピュータ支援的証明（正のトポロジカルエントロピーを有する）を可能にする。

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法を改善し、区間包摂における過大評価を低減することを目的とする。
• 常微分方程式系におけるカオス的ダイナミクスの厳密なコンピュータ支援的証明を可能にすることを目的とする。
• 陰補正を用いた高次精度数値積分により、鋭く信頼性の高い境界を提供することを目的とする。

提案手法

• アルゴリズムは、時間ステップごとに解を進めるために高次テイラー法を予測子として用いる。
• ヘルミート＝オブレシュコフ補間に基づく陰補正子を用いて予測された解を精緻化し、局所切捨て誤差を低減する。
• 区間演算を用いてすべての数値誤差を厳密に境界付けることで、検証付き包摂を保証する。
• 解の包摂における $C^1$-滑らかさを維持することで、厳密な力学系解析に不可欠な性質を確保する。
• 予測子・補正子戦略と高次精度を組み合わせることで収束性を向上させ、過大評価を低減する。
• 本手法は、特にカオス的不変集合の検出を目的とした力学系におけるコンピュータ支援的証明に適用可能であるように設計されている。

実験結果

リサーチクエスチョン

• RQ1 高次陰解法は、既存の $C^1$-Lohner法と比較して、変分常微分方程式の解の検証付き包摂における過大評価を低減できるか？
• RQ2 $Χ^1$-HOアルゴリズムの向上した精度は、常微分方程式系におけるカオス的ダイナミクスの厳密な検出を可能にするか？
• RQ3 本アルゴリズムを用いて、ローラン系における一様双曲的不変集合の存在を検証できるか？
• RQ4 ローラン系の不変集合に制限されたカオス的ダイナミクスのトポロジカルエントロピーは何か？
• RQ5 境界の鋭さと計算コストの観点から、本アルゴリズムは先行手法と比較してどのように性能を発揮するか？

主な発見

• $Χ^1$-HOアルゴリズムは、$C^1$-Lohner法と比較して、変分解の検証付き包摂をより鋭く得ており、過大評価が低減している。
• 本アルゴリズムは、ローラン系の変分方程式について、厳密な境界を計算するのに成功した。
• コンピュータ支援的証明により、ローラン系に吸引子集合が存在することが確認された。
• その吸引子には、カオス的ダイナミクスを有する一様双曲的不変部分集合が含まれている。
• この部分集合上のダイナミクスは、有限型の部分シフトと位相的に共役であり、正のトポロジカルエントロピーを確認した。
• 本手法により、検証付き数値計算を用いたローラン系におけるカオス的挙動の初の厳密な検証が可能になった。