[論文レビュー] A Lohner-type algorithm for control systems and ordinary differential inclusions

Tomasz Kapela, Piotr Zgliczyński|ArXiv.org|Dec 6, 2007

Numerical methods for differential equations参考文献 9被引用数 4

ひとこと要約

本稿では、微分方程式の摂動を扱えるように $C^0$-Lohner 法を拡張することで、制御系および常微分包含の到達集合を厳密に計算する Lohner 形アルゴリズムを提示する。この手法は、成分ごとの推定と対数ノルムを用いて不確実性下での解の境界を定め、Poincaré 写像の厳密な計算を可能にし、$ u \in [0.032, 0.127]$ の範囲で Kuramoto-Sivashinsky PDE に複数の周期的軌道の存在を証明する。このアプローチは CAPD ライブラリに実装されており、有界な摂動に対して従来の手法を上回る精度を達成している。

ABSTRACT

We describe a Lohner-type algorithm for the computation of rigorous upper bounds for reachable set for control systems, solutions of ordinary differential inclusions and perturbations of ODEs.

研究の動機と目的

• 制御系および有界な不確実性を伴う微分包含における到達集合を厳密に計算する数値的手法の開発を目的とする。
• $C^0$-Lohner アルゴリズムを摂動付き常微分方程式および誤差境界を保証する非-autonomous 系に拡張することを目的とする。
• 非-autonomous 常微分方程式および小さな遅れを伴う系に対して、断面における Poincaré 写像の正確な計算を可能とすることを目的とする。
• Kuramoto-Sivashinsky 方程式のような散逸型 PDE における周期的軌道の存在を証明するための計算フレームワークを提供することを目的とする。
• ベンチマーク問題に対して、対数ノルム (LN) 法と成分ごとの (CW) 法の2つの推定戦略を比較・検証することを目的とする。

提案手法

• 区間演算とテイラー型モデルを用いて、$|y(t)| < \epsilon$ を満たす摂動付き常微分方程式 $x'(t) = f(x(t), y(t))$ の解の厳密な包摂を計算する。
• 片側リプシッツ条件に基づく成分ごとの推定と対数ノルムを組み合わせ、不確実性下での解の成長に対する鋭い境界を導出する。
• 主な革新点として、対数ノルムと片側リプシッツ定数の間の新しい不等式を提案し、境界の鋭さを向上させる。
• 時間ステップごとに区間ニュートン型補正と区間ハロー伝播を用いて、解の包摂を維持する。
• 断面 $\alpha(x) = 0$ に到達するまで軌道を追跡することで、Poincaré 写像の計算を可能にする。
• 実装は CAPD ライブラリに統合されており、$C^0$-Lohner 法と拡張された成分ごとの推定戦略の両方をサポートする。

実験結果

リサーチクエスチョン

• RQ1 制御系および有界な摂動を伴う微分包含に対して、Lohner 形アルゴリズムを厳密に拡張できるか？
• RQ2 摂動付き常微分方程式において、成分ごとの (CW) 法と対数ノルム (LN) 法の推定戦略は、精度と計算コストの観点でどのように比較できるか？
• RQ3 本アルゴリズムは、非-autonomous 常微分方程式および小さな遅れを伴う系に対して、厳密な Poincaré 写像を計算できるか？
• RQ4 本アルゴリズムは、Kuramoto-Sivashinsky PDE において複数の周期的軌道の存在を証明する上で、どの程度の性能を示すか？
• RQ5 摂動の大きさがベクトル場の成分の大きさを上回る場合でも、本手法は信頼性のある包摂を提供できるか？

主な発見

• 摂動付き調和振動子（$\epsilon = 0.1$）において、成分ごとの (CW) 法が対数ノルム (LN) 法よりも鋭い包摂を生成し、時間ステップを 100 から 100,000 に改善した際、直径が 1.620 から 0.820 に減少した。
• Rössler 系に $[-10^{-4}, 10^{-4}]^3$ の摂動を加えた場合、断面 $x=0, x'>0$ における厳密な包摂の直径は (0.619, 0.322, 0.003) に達した。
• Kuramoto-Sivashinsky PDE に対して、$\nu \in [0.032, 0.127)$ の範囲で、安定な軌道と不安定な軌道（1つの不安定方向を有する）を含む複数の周期的軌道の存在を証明した。
• 摂動が $\epsilon = 10$ に達しても、包摂は有効であり、直径は比例的に増大したが、計算可能でかつ厳密なままであった。
• すべてのテスト問題において、CW 法が LN 法を一貫して上回り、特に大きな時間ステップや摂動サイズの下で顕著に優れた精度を示した。
• 調和振動子の $T = 2\pi$ 後の解の直径推定において、細かい時間離散化により LN 法で約 1.620、CW 法で約 0.820 に収束した。