[論文レビュー] Higher Order Methods for Differential Inclusions

Sanja Gonzalez Zivanovic, Pieter Collins|arXiv (Cornell University)|Jun 28, 2012

Numerical methods for differential equations被引用数 5

ひとこと要約

本稿では、Fliess型展開と対数ノルムに基づく誤差推定を用いて、微分包含における到達集合の厳密な上界近似を計算する高次数値法を提示する。この手法は一次手法よりも良好な局所誤差境界を達成し、不確実性やノイズを含む入力を持つ系においても、成分ごとかつ一様な誤差制御を維持しながら、長時間にわたる統合の精度を向上させる。

ABSTRACT

We present a numerical method for rigorous over-approximation of a reachable set of differential inclusions. The method gives high-order error bounds for single step approximations and a uniform bound on the error over the finite time interval. We provide formulas for the local error based on Lipschitz constants and bounds on higher-order derivatives. The method is based on a Fliess-like expansion, and extends previous results by providing error estimates which are valid for all possible inputs.

研究の動機と目的

• 微分包含における到達集合の厳密な上界近似を計算する高次数値法の開発を目的とする。
• 制御系やノイズを含む系において、すべての可能な入力に対して有効な一様的で高次の誤差境界を提供することで、一次手法を拡張することを目的とする。
• コンパクトで凸値をとる右辺をもつ微分包含の長時間統合における精度と効率の向上を目的とする。
• 不確実性、不連続なダイナミクス、またはモデル還元誤差を含む系の信頼性のある検証を可能にすることを目的とする。
• 到達集合計算における局所誤差の精度と計算コストのトレードオフを調査することを目的とする。

提案手法

• 時間ステップごとの解の時間的変化をモデル化するために、Fliessに類似したべき級数展開を用いる。
• リプシッツ定数よりもタイトな成分ごとの誤差境界を導出するために、対数ノルムを適用する。これにより、安定なダイナミクスと不安定なダイナミクスを区別できる。
• 高階微分の上限を評価し、一様な入力バウンディングを用いることで、2次および3次誤差推定を計算する。
• 各時間区間における近似解と一様な誤差境界を組み合わせることで、到達集合の上界近似を構築する。
• グリッドベースの手法を避けるために、適応的時間刻み幅を用いた1段階統合スキームを採用し、効率を向上させる。
• 関数計算と自動誤差伝播を用いて、状態空間における到達集合の閉包を計算する。

実験結果

リサーチクエスチョン

• RQ1 多値ダイナミクスをもつ微分包含に、高次誤差推定を効果的に適用できるか？
• RQ2 リプシッツに基づく境界と比較して、対数ノルムに基づく境界はタイトさと適用可能性においてどのように異なるか？
• RQ3 高次到達集合計算において、局所誤差の精度と計算コストのトレードオフはいかなるものか？
• RQ4 ノイズを含む入力やモデル還元誤差を含む系に対し、この手法は信頼性のある上界近似を提供できるか？
• RQ5 不確実パラメータをもつロシュル方程式のような複雑な系において、この手法はどのように性能を発揮するか？

主な発見

• 摂動を受ける調和振動子に対して、3次誤差推定では h=0.01 のとき到達集合半径が 0.00100759 となり、[17] が2ノルム対数ノルムを用いて得た 0.00100502 と同等の結果を得た。
• h=0.25 のとき、2次誤差推定では半径が 0.0420586 となったが、[17] の 0.0284025 よりわずかに悪い結果であったが、h を小さくすると著しく改善された。
• a=5.7 で、v_i ∈ [-10⁻⁴, 10⁻⁴] のノイズを含むロシュル方程式に対して、h=0.005、T=11.1 の範囲で R(X₀,T) ≈ {0}×[-0.15572, 0.15391]×[-3.75926, -3.41772]×[0.03139, 0.03398] を計算した。
• 解析的誤差境界は e = 8.586×10⁻⁸ であり、閉包における高い精度を示した。
• 手計算による境界ではなく自動誤差推定を用いているにもかかわらず、[17] と同等の品質の結果を得た。
• この手法は、不確実な入力をもつ複雑でカオス的系に対しても、高次誤差推定が実用的かつ有効であることを示しており、厳密な上界近似を維持している。