[論文レビュー] METHODS FOR VERIFIED SOLUTIONS TO CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATIONS

Tayyebe Haqiri, Federico Poloni|arXiv (Cornell University)|Sep 7, 2015

Numerical Methods and Algorithms参考文献 26被引用数 3

ひとこと要約

本稿では、連続時間代数リカッチ方程式（CAREs）の安定化解を計算する検証済み数値的手法を提示する。Krawczyk法の強化版と、AD Iにインspiredされた固定点反復に基づく、O(n³)の直接検証技術が用いられている。この手法は数学的に厳密な誤差境界を保証し、標準ベンチマークにおいて競争力のある性能を示す。

ABSTRACT

We describe a procedure based on the Krawczyk method to compute a verified enclosure for the stabilizing solution of a continuous- time algebraic Riccati equation AX + XA + Q = XGX, building on the work of (B. Hashemi, SCAN 2012) and adding several modifications to the Krawczyk procedure. Moreover, we describe a new O(n 3 ) direct method for verification, based on a fixed-point formulation of the equation inspired by the ADI procedure. The resulting methods are tested on a number of standard benchmark examples, and are competitive with the state-of-the-art methods for the same problem.

研究の動機と目的

• 連続時間代数リカッチ方程式（CAREs）の安定化解に対する数学的に検証された包含域を提供し、浮動小数点演算下での計算の正しさを保証すること。
• 検証数値計算分野における先行研究を基盤とし、CAREsに適した信頼性と収束性を向上させる改良版Krawczyk法を拡張すること。
• AD I法にインspiredされた固定点構造を活用し、O(n³)の計算複雑度を持つ新しい直接検証手法を開発することで、大規模問題への検証解のスケーラビリティを実現すること。
• 提案手法を標準ベンチマーク例に対してテスト・検証し、最先端の検証済みおよび非検証済みソルバーと比較して、競争力のある性能を示すこと。

提案手法

• 本稿では、区間演算とニュートン型補正を組み合わせた、変更されたKrawczyk法を用いて、CARE AX + XA + Q = XGX の安定化解の検証済み包含域を計算する。
• 区間境界における過大評価を低減し、収束性を向上させるために、改良されたプリコンディショニングと反復的精錬をKrawczyk手順に導入する。
• AD I（交互方向陰解法）反復の構造に着想を得た固定点定式化に基づく、新しい直接検証手法を提案する。これにより、効率的かつ検証可能な解の計算が可能になる。
• 固定点定式化により、O(n³)のアルゴリズム的複雑度が達成され、反復的検証手法と比較して大幅にスケーラビリティが向上する。
• 計算誤差の数学的保証を得るために、区間演算を用いて解の包含域を生成し、真の解を数学的に確実に含むことを保証する。
• 後退誤差解析と区間ニュートン技術を統合した検証フレームワークにより、計算された包含域内に安定化解が存在し、一意であることが確認される。

実験結果

リサーチクエスチョン

• RQ1 Krawczyk法は、連続時間代数リカッチ方程式の安定化解に対する検証済み包含域を効果的に計算するために、適切に適合および強化可能か？
• RQ2 AD I法にインスパイアされた固定点構造を活用することで、CAREsに対してO(n³)の複雑度を持つ直接検証手法を構築可能か？
• RQ3 提案手法の精度、信頼性、計算コストは、既存の最先端ソルバーと比較してどの程度優れているか？
• RQ4 Krawczyk手順に対する改良が、検証数値計算におけるCAREsの過大評価をどの程度低減し、収束性を向上させるか？

主な発見

• 強化されたKrawczyk法は、標準実装と比較して収束性が向上し、過大評価が低減した検証済み包含域を効果的に計算できた。
• 提案されたO(n³)の直接検証手法は、反復的検証手法と比較してスケーラブルで効率的な代替手段を提供し、大規模問題への検証解の実現可能性を高めた。
• 標準ベンチマーク例において、提案手法は検証済みの結果を達成し、解の精度と計算コストの観点で競争力のある性能を示した。
• AD Iにインスパイアされた反復に基づく固定点定式化により、数学的厳密性を維持しながら計算オーバーヘッドを最小限に抑えた、強固で効率的な検証フレームワークが実現された。
• 区間演算とニュートン型補正の組み合わせにより、計算された解の包含域が真の安定化解を数学的に保証して含むことが確認された。
• 結果として、検証済みCAREs解は、高い信頼性と競争力のある効率性を併せ持つことが確認され、検証数値線形代数分野における最先端技術が前進した。