[論文レビュー] Numerical computation of formal solutions to interval linear systems of equations

Sergey P. Shary|arXiv (Cornell University)|Mar 15, 2019

Numerical Methods and Algorithms被引用数 4

ひとこと要約

本稿では、行列分割に基づく定常反復解法を用いて、クラスカー区間算術における区間線形システムの形式的解を数値的に計算するための数値的手法を提示する。提案されたARMSplit法は、大規模なシステムに対し高い精度（最大13桁の有効数字）を達成し、収束特性を分析・他の部分微分ニュートン法と比較することで、特異的または悪条件な行列に対しても頑健性を示す。

ABSTRACT

The work is devoted to the development of numerical methods for computing "formal solutions" of interval systems of linear algebraic equations. These solutions are found in Kaucher interval arithmetic, which extends and completes the classical interval arithmetic algebraically. The need to solve these problems naturally arises, for example, in inner and outer estimation of various solution sets to interval linear systems of equations. The work develops two approaches to the construction of stationary iterative methods for computing the formal solutions that are based on splitting the matrix of the system. We consider their convergence and implementation issues, compare with the other approaches to computing formal solutions.

研究の動機と目的

• クラスカー区間算術における区間線形システムの形式的解を効率的に計算するための数値的手法の開発。
• 行列分割に基づく定常反復法の収束性および実装の分析。
• ARMSplit法と他の手法（部分微分ニュートン法など）の性能および適用範囲の比較。
• 特異的または悪条件な区間線形システム、特に大規模システムに対して、これらの手法の頑健性の評価。

提案手法

• 本稿では、クラスカー算術における区間線形システムを解くために、行列分割技術を用いて定常反復法を構築する。
• ARMSplitアルゴリズムは、区間行列の点分割に基づいて開発され、形式的解の反復的近似を可能にする。
• 反復スキームの収束解析が行われ、クラスカー算術における収束の理論的条件が導出された。
• 比較のため、収束性を向上させるためにダミング係数τ = 1を用いた部分微分ニュートン法が適用された。
• 数値実験では、三重対角行列やネイマイヤー行列を含む標準的なテストセットからの区間行列が使用され、右辺は正確な解の摂動から導出された。
• 実装は、クラスカー算術演算においてシステムを正確に満たす区間ベクトルの計算に焦点を当てた。

実験結果

リサーチクエスチョン

• RQ1 どのようにしてクラスカー区間算術における区間線形システムの形式的解を効率的に計算できるか？
• RQ2 行列分割に基づく定常反復法の収束特性は、クラスカー算術においてどのように特徴づけられるか？
• RQ3 形式的解の計算において、ARMSplit法は部分微分ニュートン法と比較して、精度および収束速度においてどのように異なるか？
• RQ4 ARMSplitのような反復法は、大規模または特異的区間線形システムに対しても、形式的解を信頼性高く計算できるか？
• RQ5 行列の構造および条件数は、形式的解計算における反復法の収束行動にどのように影響を与えるか？

主な発見

• 2×2システムでは、ARMSplit法は10反復後に3桁の有効数字、20反復後に6桁の有効数字を達成し、従来の反復法と同等の効率性を示した。
• 40×40三重対角区間システムでは、ARMSplit法は元の行列および双対化された行列の両方に対して16反復で12〜13桁の正確な有効数字を達成した。
• ネイマイヤー40×40システムでは、ARMSplit法は40反復後に約10⁻⁸の精度に達し、行列が特異であっても形式的解を正常に計算できた。
• 収束しない7×7システムでは、部分微分ニュートン法はτ = 1を用いて9反復で収束し、収束性の頑健性が優れていることが示された。
• 7×7行列の(7,7)成分を[8,82]に狭めた場合、ARMSplit法は収束が遅く、5桁の有効数字を得るのには約100反復を要した。
• 部分微分ニュートン法は、収束速度および適用範囲の両面でARMSplit法を上回り、特に反復法が発散または収束が遅い場合に顕著に優れた性能を示した。