[論文レビュー] Limit Distribution of Eigenvalues for Random Hankel and Toeplitz Band Matrices

Dang-Zheng Liu, Zhengdong Wang|ArXiv.org|Apr 20, 2009

Random Matrices and Applications参考文献 13被引用数 3

ひとこと要約

本稿は、行列サイズ $N \to \infty$ のとき、バンド幅 $b_N \to \infty$ で $b_N/N \to b \in [0,1]$ となるランダムなヘルミート・トーペリッツ行列および実対称ハンケル帯行列について、固有値分布の弱収束およびほとんど確実収束を確立する。極限分布 $\gamma_T(b)$ および $\gamma_H(b)$ は普遍的で対称的であり、それらのモーメントはペア分割上の積分として表現される。特に $\gamma_T(0)$ は正規分布であり、$\gamma_H(0)$ は $|x|e^{-x^2}$ である。$b > 0$ における異なる峰のモーメントが非正規的行動を確認する。この手法はトレースモーメント解析と、シフト作用素を用いたトーペリッツ行列の新規な線形分解に依存する。

ABSTRACT

Consider real symmetric, complex Hermitian Toeplitz and real symmetric Hankel band matrix models, where the bandwidth $b_{N} a \iy$ but $b_{N}/N o b$, $b\in [0,1]$ as $N o \infty$. We prove that the distributions of eigenvalues converge weakly to universal, symmetric distributions $γ_{_{T}}(b)$ and $γ_{_{H}}(b)$. In the case $b>0$ or $b=0$ but with the addition of $b_{N}\geq C N^{{1/2}+ε_{0}}$ for some positive constants $ε_{0}$ and $C$, we prove almost sure convergence. The even moments of these distributions are the sum of some integrals related to certain pair partitions. In particular, when the bandwidth grows slowly, i.e. $b=0$, $γ_{_{T}}(0)$ is the standard Gaussian distribution and $γ_{_{H}}(0)$ is the distribution $|x| \exp(-x^{2})$. In addition, from the fourth moments we know that the $γ_{_{T}}(b)$'s are different for different $b$'s, the $γ_{_{H}}(b)$'s different for different $b\in [0,{1/2}]$ and the $γ_{_{H}}(b)$'s different for different $b\in [{1/2},1]$.

研究の動機と目的

• ランダムなトーペリッツおよびハンケル帯行列の固有値の極限分布の存在および普遍性を確立すること。
• 極限スペクトル分布がバンド幅比 $b = \lim b_N/N$ にどのように依存するかを分析すること。
• ペア分割上の積分を用いて、極限分布のモーメントの明示的公式を導出すること。
• バンド幅 $b_N$ に対するやや緩い増加条件のもとで、ほとんど確実収束を証明すること。
• $b$ の異なる値に対して $\gamma_T(b)$ と $\gamma_H(b)$ の極限分布を、高次モーメントを用いて区別すること。

提案手法

• トーペリッツおよびハンケル帯行列を、確率的係数を伴う決定的シフト行列 $B$ と $F$ の線形結合としてモデル化すること。
• トレースモーメント法を用いて $\mathbb{E}[\mathrm{tr}(X_N^k)]$ を計算し、第4モーメントのフラクチュエーションを分析して収束を確立すること。
• 極限分布の $k$-次モーメントを、$\pi \in \mathcal{P}_2(2k)$ のペア分割上の積分の和として表現すること。
• $[0,1] \times [-1,1]^k$ 上の積分 $p_\pi(b)$ を定義し、変数 $x_l$ に関する線形形式の指示関数を含め、バンド幅依存性を捉えること。
• モーメント収束技術を適用し、高次キュムラントを制御することで、$b_N \geq C N^{1/2 + \epsilon_0}$ のもとでほとんど確実収束を証明すること。
• シフト作用素表現を用いてトーペリッツ行列の代数的構造を活用し、モーメント公式を導出すること。

実験結果

リサーチクエスチョン

• RQ1 ランダムなトーペリッツ帯行列の固有値分布は、$N \to \infty$ かつ $b_N/N \to b$ のとき、弱収束するか？
• RQ2 $\gamma_T(b)$ および $\gamma_H(b)$ の極限スペクトル分布は、すべての $b \in [0,1]$ に対して普遍的かつ対称的か？
• RQ3 $\gamma_T(b)$ および $\gamma_H(b)$ の第4モーメントは $b$ にどのように依存するか。これはそれらの相違性に何を示唆するか？
• RQ4 $b_N$ にどのような条件下で、固有値分布の収束が弱収束ではなく、ほとんど確実収束となるか？
• RQ5 インデックスのペア分割上の積分を用いて、極限分布を明示的に特徴づけられるか？

主な発見

• 極限固有値分布 $\gamma_T(b)$ は普遍的で対称的であり、$\gamma_T(0)$ は標準正規分布に一致する。
• ハンケル行列の極限分布 $\gamma_H(b)$ は $\gamma_H(0) = |x| e^{-x^2}$ を満たし、非正規的で重たい尾を持つ形状を示す。
• $b \in (0,1]$ のとき、$\gamma_T(b)$ は正規分布でない。これは $[0,1]$ 上で第4モーメント $m_4(\gamma_T(b))$ が厳密に減少することによって確認される。
• $\gamma_H(b)$ の第4モーメントは $[0,1/2]$ で厳密に減少し、$[1/2,1]$ で厳密に増加するため、各区間内で異なる $b$ に対して $\gamma_H(b)$ が区別可能であることが証明される。
• $\gamma_T(b)$ および $\gamma_H(b)$ の偶数次モーメントは、ペア分割上の積分の和として与えられ、$k=2$ について明示的に計算されている。
• $b_N \geq C N^{1/2 + \epsilon_0}$（$C>0$、$\epsilon_0>0$）の条件下で、行列要素のモーメント条件のもと、固有値分布のほとんど確実収束が確立されている。