[论文解读] The rank 1 real Wishart spiked model

M. Y. Mo|arXiv (Cornell University)|Jan 26, 2011

Random Matrices and Applications参考文献 68被引用 15

一句话总结

本文研究了秩1实Wishart尖刺模型，通过一种新颖的针对秩1矩阵的哈里什-钱德拉-伊茨伊克森-祖伯积分的围道积分公式，推导出在相变点处最大特征值的渐近分布。关键成果是基于Painlevé超越函数的新弗雷德霍姆行列式公式，给出了在大N极限下最大特征值分布的积分表示，该表示与以往基于随机算子的方法不同。

ABSTRACT

In this paper, we consider N-dimensional real Wishart matrices Y in the class $W_{\mathbb{R}}(Σ,M)$ in which all but one eigenvalues of $Σ$ is 1. Let the non-trivial eigenvalue of $Σ$ be $1+τ$, then as N, $M ightarrow\infty$, with $M/N=γ^2$ finite and non-zero, the eigenvalue distribution of $Y$ will converge into the Marchenko-Pastur distribution inside a bulk region. When $τ$ increases from zero, one starts to see a stray eigenvalue of Y outside of the support of the Marchenko-Pastur density. As the this stray eigenvalue leaves the bulk region, a phase transition will occur in the largest eigenvalue distribution of the Wishart matrix. In this paper we will compute the asymptotics of the largest eigenvalue distribution when the phase transition occur. We will first establish the results that are valid for all N and M and will use them to carry out the asymptotic analysis. In particular, we have derived a contour integral formula for the Harish-Chandra Itzykson-Zuber integral $\int_{O(N)}e^{- r(XgYg^T)}g^T\D g$ when X, Y are real symmetric and Y is a rank 1 matrix. This allows us to write down a Fredholm determinant formula for the largest eigenvalue distribution and analyze it using orthogonal polynomial techniques. As a result, we obtain an integral formula for the largest eigenvalue distribution in the large N limit characterized by Painleve transcendents. The approach used in this paper is very different from a recent paper (Bloemendal and Virag ArXiv:1011.1877), in which the largest eigenvalue distribution was obtained using stochastic operator method. In particular, the Painleve formula for the largest eigenvalue distribution obtained in this paper is new.

研究动机与目标

• 分析当 N, M → ∞ 且 M/N = γ² 为有限值时，秩1实Wishart尖刺模型中最大特征值的渐近行为。
• 推导当一个矩阵为秩1且实对称时，哈里什-钱德拉-伊茨伊克森-祖伯积分的围道积分公式。
• 利用正交多项式技术，建立最大特征值分布的弗雷德霍姆行列式公式。
• 利用Painlevé超越函数，表征在 τ = γ⁻¹ 处最大特征值分布的相变行为。
• 在大N极限下，提供一个与先前随机算子方法不同的最大特征值分布的新解析公式。

提出的方法

• 推导当一个矩阵为秩1且实对称时，实对称哈里什-钱德拉-伊茨伊克森-祖伯积分的围道积分表示，从而实现特征值统计的精确计算。
• 利用正交多项式技术分析最大特征值分布的弗雷德霍姆行列式表示。
• 为再生核构造黎曼-希尔伯特问题，通过单值性分析将其与Painlevé II方程联系起来。
• 建立再生函数与Painlevé II方程的Hastings-McLeod解之间的联系。
• 对黎曼-希尔伯特问题应用渐近分析，识别出次主导项以确保正确的归一化与可解性。
• 推导出在大N极限下最大特征值分布的积分公式，其参数由Painlevé超越函数表示。

实验结果

研究问题

• RQ1 在 τ = γ⁻¹ 的相变点，秩1实Wishart尖刺模型中最大特征值分布的行为如何？
• RQ2 当一个矩阵为秩1且实对称时，能否为哈里什-钱德拉-伊茨伊克森-祖伯积分推导出围道积分公式？
• RQ3 该模型中最大特征值分布的弗雷德霍姆行列式表示是什么？
• RQ4 在相变点，最大特征值的渐近行为与Tracy-Widom分布有何不同？
• RQ5 在大N极限下，最大特征值分布能否用Painlevé超越函数表示？

主要发现

• 推导出当一个矩阵为秩1且实对称时，哈里什-钱德拉-伊茨伊克森-祖伯积分的新围道积分公式，从而实现特征值统计的精确计算。
• 通过正交多项式技术和黎曼-希尔伯特问题分析，建立了最大特征值分布的弗雷德霍姆行列式公式。
• 最大特征值分布的相变发生在 τ = γ⁻¹，与先前的复数和四元数模型一致。
• 最大特征值的渐近分布由一个涉及Painlevé超越函数的新积分公式表征，该公式与Tracy-Widom分布不同。
• 该方法与先前的随机算子方法有本质不同，为最大特征值分布提供了基于Painlevé的新公式。
• 黎曼-希尔伯特问题的解由Airy函数的渐近行为以及Painlevé II方程的单值性结构唯一确定。