[论文解读] Gaussian fluctuations for products of random matrices

Vadim Gorin, Yi Sun|arXiv (Cornell University)|Dec 16, 2018

Random Matrices and Applications被引用 7

一句话总结

该论文在右酉不变性条件下，建立了 $N\times N$ 随机矩阵乘积中奇异值全局波动的中心极限定理。通过多变量贝塞尔生成函数及贝塞尔函数的渐近分析，证明了波动收敛到高斯场：当 $M$ 固定且 $N\to\infty$ 时为对数相关，当 $M,N\to\infty$ 联合时具有白噪声分量。

ABSTRACT

We study global fluctuations for singular values of $M$-fold products of several right-unitarily invariant $N imes N$ random matrix ensembles. As $N o \infty$, we show the fluctuations of their height functions converge to an explicit Gaussian field, which is log-correlated for $M$ fixed and has a white noise component for $M o \infty$ jointly with $N$. Our technique centers on the study of the multivariate Bessel generating functions of these spectral measures, for which we prove a central limit theorem for global fluctuations via certain conditions on the generating functions. We apply our approach to a number of ensembles, including square roots of Wishart, Jacobi, and unitarily invariant positive definite matrices with fixed spectrum, using a detailed asymptotic analysis of multivariate Bessel functions to verify the necessary conditions.

研究动机与目标

• 理解当 $N\to\infty$ 时，$N\times N$ 随机矩阵乘积中奇异值的全局波动。
• 为这类矩阵乘积的对数谱高度函数建立中心极限定理。
• 根据 $M$（矩阵因子个数）和 $N$ 的不同情形，表征极限高斯场结构——对数相关或白噪声。
• 通过多变量贝塞尔函数的渐近分析，验证中心极限定理所需的条件。
• 将该框架应用于具体系综，包括威沙特（Wishart）、雅可比（Jacobi）以及具有固定谱的酉不变正定矩阵。

提出的方法

• 使用多变量贝塞尔生成函数 $\phi_X(s)$ 作为关键工具，以编码随机矩阵的谱测度。
• 对解析函数应用对称化程序，以处理不变性并推导矩条件。
• 采用图示演算追踪生成函数展开中的组合结构。
• 依赖 Cholesky 分解 $X = R^*R$，将 $\phi_X(s)$ 表示为对角线元素 $R_{kk}$ 的期望，从而与 $S$-变换建立联系。
• 利用 Gelfand-Tsetlin 多面体和积分公式，对多变量贝塞尔函数进行详细的渐近分析。
• 通过验证协方差结构和使用累积量法证明高斯性，确立归一化波动收敛到高斯场。

实验结果

研究问题

• RQ1 当 $N\to\infty$ 时，$M$ 重随机矩阵乘积的奇异值全局波动行为如何？
• RQ2 在 $M$ 固定、$N\to\infty$ 的情形下，对数谱高度函数的极限高斯场结构是什么？
• RQ3 当 $M$ 和 $N$ 联合趋于无穷时，极限高斯场如何变化？
• RQ4 极限高斯场在何种条件下表现出对数相关性？
• RQ5 能否从 Cholesky 分解元素的渐近行为中恢复极限谱测度的 $S$-变换？

主要发现

• 当 $M$ 固定时，高度函数 $\mathcal{H}_N(t)$ 的波动在 $N\to\infty$ 时收敛到具有对数相关协方差核 $K(t,s) = -\frac{1}{2\pi^2}\log|t-s| + O(1)$ 的高斯场。
• 当 $M,N\to\infty$ 联合时，$M^{1/2}(\mathcal{H}_N(t) - \mathbb{E}[\mathcal{H}_N(t)])$ 的波动收敛到具有白噪声分量的高斯场，表明在小尺度上波动不相关。
• 通过基于多变量贝塞尔生成函数及其渐近性质的创新框架，建立了全局波动的中心极限定理。
• $S$-变换作为 $\mathbb{E}[2\log R_{\lfloor tN\rfloor,\lfloor tN\rfloor}]$ 的极限被恢复，将 Cholesky 分解与自由概率联系起来。
• 该方法适用于广泛系综，包括威沙特（Wishart）、雅可比（Jacobi）以及具有固定谱的酉不变正定矩阵。
• 对数相关性的关键技术条件是谱测度极限形状的光滑性，这确保了 $S$-变换行为良好。