[論文レビュー] The polynomial method for random matrices

N. Raj Rao, Alan Edelman|ArXiv.org|Jan 16, 2006

Random Matrices and Applications被引用数 11

ひとこと要約

本稿では、ランダム行列の多項式法を導入し、Stieltjes変換を符号化する2変数多項式方程式を用いて、極限固有値分布の正確な計算を可能にする。主な貢献は、Wigner や Wishart および関連アンサンブルに応用可能な、代数的ランダム行列のモーメント、密度、自由畳み込みを効率的に計算するための体系的フレームワークを提供することであり、RMTool ソフトウェアで実装されている。

ABSTRACT

We define a class of "algebraic" random matrices. These are random matrices for which the Stieltjes transform of the limiting eigenvalue distribution function is algebraic, i.e., it satisfies a (bivariate) polynomial equation. The Wigner and Wishart matrices whose limiting eigenvalue distributions are given by the semi-circle law and the Marcenko-Pastur law are special cases. Algebraicity of a random matrix sequence is shown to act as a certificate of the computability of the limiting eigenvalue density function. The limiting moments of algebraic random matrix sequences, when they exist, are shown to satisfy a finite depth linear recursion so that they may often be efficiently enumerated in closed form. In this article, we develop the mathematics of the polynomial method which allows us to describe the class of algebraic matrices by its generators and map the constructive approach we employ when proving algebraicity into a software implementation that is available for download in the form of the RMTool random matrix "calculator" package. Our characterization of the closure of algebraic probability distributions under free additive and multiplicative convolution operations allows us to simultaneously establish a framework for computational (non-commutative) "free probability" theory. We hope that the tools developed allow researchers to finally harness the power of the infinite random matrix theory.

研究の動機と目的

• Wigner や Marčēnko-Pastur のような古典的分布を超える広範なランダム行列クラスの極限固有値分布を計算する体系的手法の開発。
• Stieltjes変換が2変数多項式方程式を満たす代数的ランダム行列—これにより、モーメントの効率的列挙と分布計算が可能であることを確立。
• 多項式変換を用いて、自由加法的および乗法的畳み込み演算を代数的枠組みで構成的に定式化。
• RMTool と呼ばれるソフトウェアパッケージにこの手法を実装し、複雑な行列アンサンブルの固有値密度およびモーメントの実用的計算を可能に。
• 代数的分布が自由畳み込みに関して閉じていることの特徴づけを通じ、自由確率論を統合的かつ拡張的に展開。

提案手法

• Stieltjes変換 $ m(z) $ が2変数多項式方程式 $ L_{\textrm{mz}}(m,z) = 0 $ を満たす行列を代数的ランダム行列と定義し、$ z \in \mathbb{C}^+ $ を満たす。
• 結果式と代数的消去法を用いて、独立な代数的ランダム行列の和および積の生成多項式から、それらのための多項式方程式を導出。
• 表7に示す操作則（例：'Add Atomic Wishart' および 'Multiply Wishart'）を用い、自由加法的および乗法的畳み込みの下で多項式方程式を変換。
• 代数的行列のモーメントが有限深さの線形再帰に従うという事実を活用し、モーメントの閉形式列挙を効率的に行う。
• 理論的枠組みを RMTool パッケージを通じてソフトウェア実装に変換し、入力された多項式から固有値密度およびモーメントを計算可能に。
• D-有限級数と組合せ的対応（例：Catalan数、Riordan数、Schröder数）を用い、行列のモーメントと既知の整数列を結びつける。

実験結果

リサーチクエスチョン

• RQ1 ランダム行列系列の極限固有値分布は、そのStieltjes変換が満たす2変数多項式方程式によって代数的に特徴づけられるか？
• RQ2 自由加法的および乗法的畳み込み演算が、それらの定義多項式の2変数多項式にどのように体系的に写像されるか？
• RQ3 代数的ランダム行列系列の極限モーメントが有限深さの線形再帰によって計算可能であるための条件は何か？
• RQ4 多項式法を用いて、Wigner や Wishart のような標準的タイプを超える複雑な行列アンサンブルの正確な固有値密度およびモーメントを計算できる範囲はどの程度か？
• RQ5 有限の行列実現から最小次数の多項式を同定できるか？また、有効な確率分布を定義するための多項式係数の許容可能な制約は何か？

主な発見

• Wigner 行列および Wishart 行列の極限固有値分布のStieltjes変換は、それぞれ $ m^2 + z m + 1 = 0 $ および $ m^2 z - (-2z + 1)m + 2 = 0 $ という2変数多項式方程式を満たす。
• 独立な Wigner 行列と Wishart 行列の和は、そのStieltjes変換が $ m^3 + (z+2)m^2 - (-2z+1)m + 2 = 0 $ を満たす極限分布をもたらし、'Add Atomic Wishart' の規則により導出される。
• 独立な Wigner 行列と Wishart 行列の積は、そのStieltjes変換が $ m^4 z^2 - 2m^3 z + m^2 + 4m z + 4 = 0 $ を満たす極限分布をもたらし、'Multiply Wishart' の規則により導出される。
• $ \mathbf{W}_N(1) - \mathbf{I}_N $ の極限モーメントは Riordan 数であり、$ 2\mathbf{W}_N(0.5) $ のものは large Schr"{o}der 数、$ 4\mathbf{W}_N(0.125) $ のものは small Schr"{o}der 数である。
• 代数的ランダム行列のクラスは、自由加法的および乗法的畳み込みに関して閉じており、既知のものから新しい代数的アンサンブルを体系的に構成可能である。
• RMTool ソフトウェアパッケージは多項式法を実装し、定義多項式から直接固有値密度およびモーメントを計算可能にしている。