[論文レビュー] Sensitivity Analysis Using a Fixed Point Interval Iteration

Alexandre Goldsztejn|ArXiv.org|Nov 18, 2008

Matrix Theory and Algorithms参考文献 5被引用数 7

ひとこと要約

本稿では、パラメトリック方程式系の厳密な感度解析のための固定点区間反復法を提案する。この手法は、Hansen-Sengupta作用素に基づくもので、パラメータ領域全体にわたり解の存在と一意性を証明し、Krawczyk作用素よりも鋭い包摂を提供する。近似ノーマル解からのロバストな解集合の境界特定が可能である。

ABSTRACT

Proving the existence of a solution to a system of real equations is a central issue in numerical analysis. In many situations, the system of equations depend on parameters which are not exactly known. It is then natural to aim proving the existence of a solution for all values of these parameters in some given domains. This is the aim of the parametrization of existence tests. A new parametric existence test based on the Hansen-Sengupta operator is presented and compared to a similar one based on the Krawczyk operator. It is used as a basis of a fixed point iteration dedicated to rigorous sensibility analysis of parametric systems of equations.

研究の動機と目的

• パラメータが与えられた区間内に不確実に存在するパラメトリック方程式系の感度解析のための厳密な手法を開発すること。
• 初期解の包摂が不要な状況でも、パラメータ領域全体にわたり解の存在と一意性を証明する課題に対処すること。
• パラメトリック拡張と固定点反復を用いたHansen-Sengupta作用素を用いることで、既存の区間ベースの存在検証を改善すること。
• 事前の解の変動範囲の知識が不要な状況でも、近似ノーマル解のみを用いて解集合の実用的計算を可能とすること。

提案手法

• 本稿では、Hansen-Sengupta存在検証のパラメトリック拡張を提案し、与えられた区間領域内のすべてのパラメータ値に対して解の存在と一意性を検証する。
• 各反復で反復誤差を要因 δ > 1 で拡張する固定点区間反復を採用し、有限ステップで収縮を保証し、存在証明を可能にする。
• 反復には、変数ベクトルの定義域を ℝⁿ に拡張することで、過剰な定義域の交差を回避するHansen-Sengupta作用素を適用する。
• 区間演算を用いて微分の包摂と関数範囲の拡張を計算し、微分可能性を仮定しない厳密な境界を実現する。
• 収束性と包摂の鋭さを向上させるために、前処理行列 C を適用する。
• 計算された包摂が直前の反復内に厳密に含まれる場合に、アルゴリズムは終了し、パラメータ領域内すべてのパラメータに対して解の存在と一意性が確認される。

実験結果

リサーチクエスチョン

• RQ1 Hansen-Sengupta作用素に基づく固定点区間反復法は、初期解の包摂が不要な状況でも、パラメトリック方程式系の解集合を厳密に包摂できるか？
• RQ2 包摂の鋭さと収束速度の観点から、パラメトリックHansen-Sengupta作用素はKrawczyk作用素と比較してどのように異なるか？
• RQ3 δ > 1 における反復的拡張が、収縮の確保と感度解析における存在証明の実現に与える影響は何か？
• RQ4 近似ノーマル解からのみ出発する場合でも、本手法はパラメータ領域全体にわたり解の存在と一意性を信頼性高く証明できるか？
• RQ5 Lipschitzまたはスロープベースの微分包摂を用いることで、非微分可能または非滑らかな系はどのように処理されるか？

主な発見

• パラメトリックHansen-Sengupta作用素は、Krawczyk作用素よりも鋭い解包摂を生成し、反復5回目における相対的幅差は約20％であり、k → ∞ に伴い0に収束する。
• 例2では、初期反復がその直前の反復に厳密に含まれない場合（過剰濃縮のため）でも、4反復目で解の存在と一意性が正当に証明された。
• δ = 1.01 を用いた反復的拡張により、[x_{k+1}] ⊆ int[x_k] となる収縮条件が有限ステップで満たされ、厳密な存在証明が可能となった。
• アルゴリズムは固定点 [x_∞] に収束し、[x_∞] = [H]([x_∞], ℝⁿ) を満たす。この固定点は、[a] 内のすべてのパラメータに対して完全な解集合を含む。
• 本手法はロバストかつ実用的であり、近似ノーマル解のみを必要とし、事前の解の変動範囲の知識が不要である。
• 包摂の幅比の収束速度は約 12exp(−0.46k) であり、0への指数的収束を示している。