[論文レビュー] Computing Enclosures of Overdetermined Interval Linear Systems

Jaroslav Horáček, Milan Hladík|arXiv (Cornell University)|Apr 17, 2013

Numerical Methods and Algorithms参考文献 8被引用数 10

ひとこと要約

本稿では、変数より方程式が多い過決定区間線形系（OILS）における解集合の包含域を計算するための複数の手法を提示・比較している。区間解析を用い、Rohn法、最小二乗法、線形計画法、ガウス・エリミネーションなどの手法を評価した結果、区間半径が狭い場合にはRohn法が最も高速であることが判明した一方、Verifylssは計算時間の増加を伴うが、よりきめの細かい包含域を提供した。

ABSTRACT

This work considers special types of interval linear systems - overdetermined systems. Simply said these systems have more equations than variables. The solution set of an interval linear system is a collection of all solutions of all instances of an interval system. By the instance we mean a point real system that emerges when we independently choose a real number from each interval coefficient of the interval system. Enclosing the solution set of these systems is in some ways more difficult than for square systems. The main goal of this work is to present various methods for solving overdetermined interval linear systems. We would like to present them in an understandable way even for nonspecialists in a field of linear systems. The second goal is a numerical comparison of all the methods on random interval linear systems regarding widths of enclosures, computation times and other special properties of methods.

研究の動機と目的

• 研究分野を問わず幅広い研究者に適した、過決定区間線形系（OILS）における解集合の包含域を計算する既存手法の包括的で理解しやすい概要を提供すること。
• ランダムに生成されたOILSインスタンスに対して、Rohn法、最小二乗法、線形計画法、反復法、ガウス・エリミネーションといったさまざまな手法の性能を評価・比較すること。
• 各手法の長所・短所および実用的適用可能性を特定すること。特に、包含域のきびつさ、計算時間、不適合性の検出能力に注目する。
• 単純または「無知な」手法が、区間半径やシステムサイズの変化に伴って、より洗練された手法を上回るケースを浮き彫りにすること。
• システムのサイズ、係数の範囲、望ましい包含域の鋭さといった特性に応じて、研究者が最適な手法を選択できるようにガイドすること。

提案手法

• 本稿では、ガウス・エリミネーション（GE）、事前条件付きガウス・エリミネーション（GEpre）、Rohnの反復法（Rohn）、最小二乗法（Verifylss）、線形計画法（LP）という5つの主要な手法を評価している。これらは区間ヘッドの計算に用いられている。
• 各手法は、次元（5×3から1000×663まで）と区間半径（≤10⁻³および≤10⁻⁵）が異なるランダムに生成された過決定区間線形系に適用されている。
• 包含域のきびつさは、得られた区間解集合の幅によって測定され、幅が小さいほど性能が優れているとされる。
• 計算時間は各手法間で記録・比較されており、特にVerifylssとRohnの実行時間比に注目している。
• Rohn法は、中点行列と半径行列を用いた反復的精錬を実行し、区間系の構造を活かして包含域を効率的に計算する。
• Verifylssは、検証付き最小二乗ソルバーを適用して鋭い包含域を計算する。一方、LPは計算コストが許容される場合には正確な区間ヘッドを計算するために用いられる。

実験結果

リサーチクエスチョン

• RQ1 さまざまなシステムサイズと区間半径において、過決定区間線形系に対してどの手法が最もきめの細かい包含域を提供するか？
• RQ2 計算時間は各手法でどのように変化するのか？また、速度と包含域の鋭さのトレードオフはどのようになるか？
• RQ3 単純またはヒューリスティックな手法（例：Rohn法）が、Verifylss や LP といったより複雑な検証付き手法を上回る条件は何か？
• RQ4 区間半径や中点範囲の変化に対して、各手法の感度はどの程度か？これは信頼性や性能に影響を与えるか？
• RQ5 Rohnのような反復法は、他のソルバーの事前処理として効果的に使用できるか？この戦略は、時間的・精度的観点からいつ有効になるか？

主な発見

• Rohn法は、テストされたすべてのシステムサイズで最も高速であり、特にmが大きくnが小さいシステム（例：894×8）ではVerifylssに比べて著しく計算時間が短い（比：423.7）。
• VerifylssはRohn法よりもきめの細かい包含域を生成するが、計算時間も著しく長くなるため、大規模システムにはあまり適さない。
• 最大半径≤10⁻³のシステムでは、200×170システムにおいてRohn法の平均時間は0.6632秒、Verifylssは9.8570秒で、14倍以上も遅い。
• VerifylssとRohnの計算時間比は、システムのスパarsity（例：894×8：423.7）が高くなると急激に増加し、Rohn法の過決定系における優位性が顕著に現れる。
• ガウス・エリミネーション（GEpre）は、速度面で著しく劣っており、大規模問題には推奨されないが、不適合性の検出には役立つ可能性がある。
• 各手法の性能は、区間半径や中点範囲の変化に対して相対的に感度が低く、さまざまな係数分布に対しても頑健であると示唆されている。