[論文レビュー] On the condition number of the critically-scaled Laguerre Unitary Ensemble

Percy Deift, Govind Menon+1 more|arXiv (Cornell University)|Jul 2, 2015

Random Matrices and Applications参考文献 22被引用数 4

ひとこと要約

本稿では、$ n = N + \lfloor \sqrt{4cN} \rfloor $ という臨界スケーリング下で、ラゲル・ユニタリーエンsemble（LUE）— すなわちランダムなワイシャルト行列 $ A = XX^* $ — の条件数 $ \kappa $ が、$ N \to \infty $ の下でトレーシー・ワイドマン分布（$ \beta=2 $）に分布収束することを確立する。この結果は、ラゲル固有多項式の漸近解析をリーマン＝ヒルベルト問題を用いて行い、極端な固有値のフラクチュエーション（最小固有値、最大固有値、および条件数）が、このスケーリング下ですべて同じ普遍的トレーシー・ワイドマン則に従うことを示している。このスケーリングは、共役勾配法の停止時刻の普遍性に由来する。

ABSTRACT

We consider the Laguerre Unitary Ensemble (aka, Wishart Ensemble) of sample covariance matrices $A = XX^*$, where $X$ is an $N imes n$ matrix with iid standard complex normal entries. Under the scaling $n = N + \lfloor \sqrt{ 4 c N} floor$, $c > 0$ and $N ightarrow \infty$, we show that the rescaled fluctuations of the smallest eigenvalue, largest eigenvalue and condition number of the matrices $A$ are all given by the Tracy--Widom distribution ($β= 2$). This scaling is motivated by the study of the solution of the equation $Ax=b$ using the conjugate gradient algorithm, in the case that $A$ and $b$ are random: For such a scaling the fluctuations of the halting time for the algorithm are empirically seen to be universal.

研究の動機と目的

• ラゲル・ユニタリーエンsemble（LUE）の条件数 $ \kappa = \lambda_{\max}/\lambda_{\min} $ の極限分布を、$ n = N + \lfloor \sqrt{4cN} \rfloor $ という臨界スケーリング下で確立すること。
• 最小固有値、最大固有値、および条件数のフラクチュエーションが、すべて $ N \to \infty $ の極限でトレーシー・ワイドマン分布（$ \beta=2 $）に収束することを示すこと。
• このスケーリングと、ランダムな線形方程式系 $ Ax = b $ を解く際の共役勾配法の停止時刻の経験的普遍性との関連を明らかにすること。
• リーマン＝ヒルベルト解析を用いてラゲル多項式の漸近挙動を厳密に導出し、極端な固有値極限定理を支持すること。

提案手法

• LUEの固有値密度に不可欠なラゲル多項式の漸近挙動を解析するためのリーマン＝ヒルベルト問題技法の使用。
• 古典的なリーマン＝ヒルベルト解析を応用し、LUEの重み関数に関連する直交多項式の精密な漸近挙動を導出すること。
• LUEのグローバル固有値密度の導出と、ラゲル重みおよび直交多項式系との関連の確立。
• 導出された多項式の漸近挙動を用いた極端な固有値の漸近解析。特にスペクトルのソフトエッジおよびハードエッジ付近のエッジスケーリングに注目する。
• 正確なスケーリング $ \alpha = \lfloor \sqrt{4cN} \rfloor $ を用いて、$ \lambda_{\min} $、$ \lambda_{\max} $、$ \kappa $ の同時極限分布がトレーシー・ワイドマン則によって記述されることを確立する。
• エアリー核を含むフレドホルム行列式を用いて、トレーシー・ワイドマン分布 $ F_2(s) = \det(I - \mathcal{K}_{\text{Ai}}|_{L^2((s,\infty))}) $ を表現すること。

実験結果

リサーチクエスチョン

• RQ1 臨界スケーリング $ n = N + \lfloor \sqrt{4cN} \rfloor $ の下で、LUEの条件数はトレーシー・ワイドマン分布に収束するか？
• RQ2 このスケーリング領域において、極端な固有値フラクチュエーション（最小固有値および最大固有値）は漸近的にどのように振る舞うか？
• RQ3 このスケーリングと、ランダム線形系における共役勾配法の停止時刻の普遍性との関連は何か？
• RQ4 $ \lambda_{\min} $、$ \lambda_{\max} $、$ \kappa $ の同時分布は、$ N \to \infty $ の極限で、単一の普遍的法則によって記述可能か？
• RQ5 $ \alpha \sim \sqrt{N} $ スケーリングが、条件数の重尾的挙動を回避するために果たす役割は何か？

主な発見

• 適切にスケーリングされた最小固有値 $ \lambda_{\min} $ は、トレーシー・ワイドマン分布に分布収束する：$ \lim_{N\to\infty} \mathbb{P}\left( \frac{c - \lambda_{\min}}{c \alpha^{-2/3} 2^{2/3}} \leq t \right) = F_2(t) $。
• 同じスケーリング下で最大固有値 $ \lambda_{\max} $ もトレーシー・ワイドマン則に収束する：$ \lim_{N\to\infty} \mathbb{P}\left( \frac{\lambda_{\max} - \nu}{\nu^{1/3} 2^{2/3}} \leq t \right) = F_2(t) $、ここで $ \nu = 4N + 2\alpha + 2 $。
• 臨界スケーリング下で条件数 $ \kappa = \lambda_{\max}/\lambda_{\min} $ もトレーシー・ワイドマン分布に収束する：$ \lim_{N\to\infty} \mathbb{P}\left( \frac{\kappa - 4N/c}{4c^{-4/3} N^{2/3}} \leq t \right) = F_2(t) $。
• 臨界スケーリング $ \alpha = \lfloor \sqrt{4cN} \rfloor $ により、$ \lambda_{\min} $ がゼロから遠ざけられ、条件数の重尾的挙動が排除され、$ \kappa $ の純粋なトレーシー・ワイドマンフラクチュエーションが実現される。
• リーマン＝ヒルベルト問題によるラゲル多項式の漸近解析が、極端な固有値極限定理を導くために必要な基盤を提供する。
• 結果は、ランダムな線形系において共役勾配法の停止時刻の経験的普遍性を説明する：臨界スケーリングにより、収束時間の非退化で普遍的な極限が得られる。