[論文レビュー] The polynomial eigenvalue problem is well conditioned for random inputs

Carlos Beltrán, Diego Armentano|arXiv (Cornell University)|Jun 19, 2017

Matrix Theory and Algorithms参考文献 22被引用数 4

ひとこと要約

本稿では、入力行列がi.i.d.標準複素ガウス分布に従う場合、多項式固有値問題（PEVP）は平均的に良好に条件付けられていることを確立する。ランダム行列理論と代数幾何学を用いて、正確な期待二乗条件数が(d+1)n²−1/dであることを導出しており、n や d が大きくても数値的解法が通常安定することを示している。

ABSTRACT

We compute the exact value of the squared condition number for the polynomial eigenvalue problem, when the input matrices have entries coming from the standard complex Gaussian distribution, showing that in general this problem is quite well conditioned.

研究の動機と目的

• 入力行列が確率的である場合の多項式固有値問題（PEVP）の平均的条件付けを特定すること。
• 関連問題の広範な研究にもかかわらず、PEVPの平均的条件数に関する理論的結果が不足しているという問題を解決すること。
• 期待二乗条件数を計算することにより、大規模問題に対してもPEVPが通常良好に条件付けられていることを確立すること。
• 代数幾何学とランダム行列理論を用いて、他の固有値問題への分析を一般化すること。
• PEVPの数値的解法が実際の応用で高精度に計算可能であり、安定しているという理論的根拠を提供すること。

提案手法

• 文献[16, Th. 4.2]で導出されたPEVPの条件数の公式を用いる。この公式は固有値、固有ベクトル、行列多項式の微分に依存する。
• 複素射影空間上でガウス分布を用いた代数的固有値問題の期待条件数に関する一般論である定理2.1を適用する。
• 代数多様体上の積分を計算するために、コアリティック・フォーミュラとノーマル・ジャコビアンを用いる。幾何学的構造と期待値計算を結びつける。
• 条件数がスケーリングに対して不変（μ(tA, (α,β)) = μ(A, (α,β)))であるという事実を活用し、分散を1に正規化できる。
• 解多様体の構造とその余次元を分析し、理論的枠組みの適用可能性を検証する。
• キーパターンの積分を計算するための補題A.1を用い、期待二乗条件数の正確な評価を可能にする。

実験結果

リサーチクエスチョン

• RQ1 入力行列がi.i.d.標準複素ガウス分布に従う場合、多項式固有値問題の期待二乗条件数は何か？
• RQ2 特定の例では悪条件になり得るが、n や d が大きくても、PEVPは平均的に良好に条件付けられていると言えるか？
• RQ3 文献[16]の理論的枠組みを拡張し、PEVPや関連問題の正確な期待条件数を計算できるか？
• RQ4 PEVPの平均的条件付けは、標準的または一般化固有値問題と比べてどう異なるか？
• RQ5 同様の手法を、固有値に代数的制約を課した他の構造的固有値問題に適用できるか？

主な発見

• i.i.d.標準複素ガウス分布に従う入力行列の下で、すべての固有値および確率的入力について平均をとったPEVPの期待二乗条件数は、正確に(d+1)n²−1/dである。
• 最大条件数μ_max(A)について、E[μ_max(A)²] ≤ ((d+1)n²−1)nが成り立つ。これは、最悪の固有値に対しても通常良好に条件付けられていることを示している。
• 対数期待条件数は(3/2)log n + (1/2)log(d+1)で有界であり、n や d が大きくても高精度な解法が可能であることを示している。
• 条件数のスケーリング不変性により、入力分布の分散σ²が任意の値であっても結果は成り立つ。
• この手法は、2n個の解を持つスパースPEVPや2n個の解を持つ三項問題など、他の固有値問題へも一般化可能であり、それぞれの期待二乗条件数は(3n²+3n)/4および3n²となる。
• 解集合がほとんど確実に有限であること、および関連する代数的多様体が定理2.1に必要な幾何的条件を満たしていることの検証により、理論的枠組みの妥当性が裏付けられている。