[論文レビュー] Three Problems Related to the Eigenvalues of Complex Non-central Wishart Matrices with a Rank-1 Mean

Prathapasinghe Dharmawansa|arXiv (Cornell University)|Jun 27, 2013

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

ひとこと要約

本稿では、ランク1の平均をもつ $ n \times n $ 複素非中心 Wishart 行列の最小固有値の累積分布関数（c.d.f.）について、新しい行列式に基づく表現を導出する。この表現は、$ m-n+1 $ 次の行列に依存することを示している。Wang [44] が提示した直交多項式技術およびコーシー積分表示を活用することで、$ \text{tr}(\mathbf{W})/\lambda_{\min} $ の分布および $ \det[z\mathbf{I}_n + \mathbf{W}]^{-1} $ の平均についても特徴づけられ、スムーズ化解析や確率行列理論への応用が可能な正確な解析的表現が得られる。

ABSTRACT

Recently, D. Wang has devised a new contour integral based method to simplify certain matrix integrals. Capitalizing on that approach, we derive a new expression for the probability density function (p.d.f.) of the joint eigenvalues of a complex non-central Wishart matrix with a rank-1 mean. The resulting functional form in turn enables us to use powerful classical orthogonal polynomial techniques in solving three problems related to the non-central Wishart matrix. To be specific, for an $n imes n$ complex non-central Wishart matrix $\mathbf{W}$ with $m$ degrees of freedom ($m\geq n$) and a rank-1 mean, we derive a new expression for the cumulative distribution function (c.d.f.) of the minimum eigenvalue ($λ_{\min}$). The c.d.f. is expressed as the determinant of a square matrix, the size of which depends only on the difference $m-n$. This further facilitates the analysis of the microscopic limit for the minimum eigenvalue which takes the form of the determinant of a square matrix of size $m-n$ with the Bessel kernel. We also develop a moment generating function based approach to derive the p.d.f. of the random variable $\frac{ ext{tr}(\mathbf{W})}{λ_{\min}}$, where $ ext{tr}(\cdot)$ denotes the trace of a square matrix. This random quantity is of great importance in the so-called smoothed analysis of Demmel condition number. Finally, we find the average of the reciprocal of the characteristic polynomial $\det[z\mathbf{I}_n+\mathbf{W}],\; |\arg z|

研究の動機と目的

• ランク1の平均をもつ複素非中心 Wishart 行列の最小固有値の c.d.f. に対する新しい解析的表現を導出すること。
• スムーズ化解析における Demmel 条件数に関連する比 $ \text{tr}(\mathbf{W})/\lambda_{\min} $ の分布を特徴づけること。
• 複素 Wishart 行列に対する特性多項式の逆行列の統計的平均を計算すること。
• $ m-n $ のみに依存する行列式表現として $ \lambda_{\min} $ の c.d.f. を提示し、微視的極限解析を可能にする。

提案手法

• Wang [44] が開発した、ランク1行列を含む2変数ハイパージェオメトリック関数のコーシー積分表示を用いる。
• 古典的直交多項式技術（特にラゲール多項式 $ L_j^{(\alpha)}(x) $）を用いて、固有値分布の表現を簡略化する。
• Khatri [31] の研究に基づく、連関固有値確率密度関数の行列式表現を、コーシー積分法により再定式化する。
• $ \lambda_{\min} $ の c.d.f. を $ (m-n+1) \times (m-n+1) $ 行列の行列式として導出し、漸近的解析を可能にする。
• モーメント母関数技術を用いて $ \text{tr}(\mathbf{W})/\lambda_{\min} $ の確率密度関数を導出。積分変換および直交多項式の恒等式を活用する。
• 極限処理および行列式の恒等式を適用して、$ \det[z\mathbf{I}_n + \mathbf{W}]^{-1} $ の平均を評価し、ラゲール多項式の積分を含む閉形式表現が得られる。

実験結果

リサーチクエスチョン

• RQ1 ランク1の平均をもつ複素非中心 Wishart 行列の最小固有値の c.d.f. は、$ m-n $ のみに依存するサイズの行列式として表現可能か？
• RQ2 比 $ \text{tr}(\mathbf{W})/\lambda_{\min} $ の正確な分布は何か？ また、これは Demmel 条件数のスムーズ化解析とどのように関係するか？
• RQ3 特性多項式 $ \det[z\mathbf{I}_n + \mathbf{W}] $ の逆行列の統計的平均は何か？ そして、閉形式で表現可能か？
• RQ4 最小固有値分布の微視的極限はどのように振る舞い、極限でどのようなカーネル構造が出現するか？

主な発見

• $ \lambda_{\min} $ の c.d.f. は、$ m-n $ が小さいとき、標準的な $ n \times n $ 行列式よりもよりコンactな形で、$ (m-n+1) \times (m-n+1) $ 行列の行列式として表現される。
• 最小固有値分布の微視的極限は、サイズ $ m-n $ の行列式として、ベッセルカーネルをとる形を取り、普遍的な端効果を示す。
• モーメント母関数のアプローチにより、$ \text{tr}(\mathbf{W})/\lambda_{\min} $ の正確な確率密度関数が得られ、これは Demmel 条件数のスムーズ化解析にとって重要である。
• $ \det[z\mathbf{I}_n + \mathbf{W}]^{-1} $ の平均は、ラゲール多項式と有理関数を含む積分の組み合わせとして導出され、直交多項式の恒等式を用いて閉形式表現が得られる。
• 導出された $ \lambda_{\min} $ の c.d.f. は、既知の結果と同等であることが示され、漸近的および数値的解析を容易にする形で表現されている。
• Wang [44] のコーシー積分法が、確率行列理論における複雑な行列積分表現の簡略化に有効であることが分析によって確認された。