On the condition number of the critically-scaled Laguerre Unitary Ensemble

We consider the Laguerre Unitary Ensemble (aka, Wishart Ensemble) of sample covariance matrices $A = XX^*$, where $X$ is an $N \times n$ matrix with iid standard complex normal entries. Under the scaling $n = N + \lfloor \sqrt{ 4 c N} \rfloor$, $c > 0$ and $N \rightarrow \infty$, we show that the rescaled fluctuations of the smallest eigenvalue, largest eigenvalue and condition number of the matrices $A$ are all given by the Tracy--Widom distribution ($\beta = 2$). This scaling is motivated by the study of the solution of the equation $Ax=b$ using the conjugate gradient algorithm, in the case that $A$ and $b$ are random: For such a scaling the fluctuations of the halting time for the algorithm are empirically seen to be universal.

[1]  J. Neumann,et al.  Numerical inverting of matrices of high order. II , 1951 .

[2]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[3]  G. Szegö Recent Advances and Open Questions on the Asymptotic Expansions of Orthogonal Polynomials , 1959 .

[4]  S. Kaniel Estimates for Some Computational Techniques - in Linear Algebra , 1966 .

[5]  T. Sugiyama On the Distribution of the Largest Latent Root and the Corresponding Latent Vector for Principal Component Analysis , 1966 .

[6]  V. Marčenko,et al.  DISTRIBUTION OF EIGENVALUES FOR SOME SETS OF RANDOM MATRICES , 1967 .

[7]  P. Krishnaiah,et al.  On the exact distribution of the smallest root of the wishart matrix using zonal polynomials , 1971 .

[8]  B. Simon Trace ideals and their applications , 1979 .

[9]  A. Edelman Eigenvalues and condition numbers of random matrices , 1988 .

[10]  A. Greenbaum Behavior of slightly perturbed Lanczos and conjugate-gradient recurrences , 1989 .

[11]  Athanassios S. Fokas,et al.  The isomonodromy approach to matric models in 2D quantum gravity , 1992 .

[12]  P. Forrester The spectrum edge of random matrix ensembles , 1993 .

[13]  C. Tracy,et al.  Level-spacing distributions and the Airy kernel , 1992, hep-th/9211141.

[14]  Stephanos Venakides,et al.  Asymptotics for polynomials orthogonal with respect to varying exponential weights , 1997 .

[15]  T. H. Baker,et al.  Random matrix ensembles with an effective extensive external charge , 1998 .

[16]  Stephanos Venakides,et al.  Strong asymptotics of orthogonal polynomials with respect to exponential weights , 1999 .

[17]  K. Johansson Shape Fluctuations and Random Matrices , 1999, math/9903134.

[18]  P. Deift Orthogonal Polynomials and Random Matrices: A Riemann-Hilbert Approach , 2000 .

[19]  I. Johnstone On the distribution of the largest eigenvalue in principal components analysis , 2001 .

[20]  W. Van Assche,et al.  The Riemann-Hilbert approach to strong asymptotics for orthogonal polynomials on [-1,1] , 2001 .

[21]  M. Vanlessen,et al.  Strong Asymptotics of Laguerre-Type Orthogonal Polynomials and Applications in Random Matrix Theory , 2005 .

[22]  Roderick Wong,et al.  Global asymptotic expansions of the Laguerre polynomials—a Riemann–Hilbert approach , 2008, Numerical Algorithms.

[23]  Ronald F. Boisvert,et al.  NIST Handbook of Mathematical Functions , 2010 .

[24]  Yang Chen,et al.  PDEs satisfied by extreme eigenvalues distributions of GUE and LUE , 2011, 1102.0402.

[25]  Danning Li,et al.  Approximation of Rectangular Beta-Laguerre Ensembles and Large Deviations , 2013, 1309.3882.

[26]  Dan Dai,et al.  Critical Edge Behavior and the Bessel to Airy Transition in the Singularly Perturbed Laguerre Unitary Ensemble , 2013, 1309.4354.

[27]  P. Deift,et al.  Universality in numerical computations with random data , 2014, Proceedings of the National Academy of Sciences.

[28]  T. Trogdon,et al.  Riemann-Hilbert Problems, Their Numerical Solution, and the Computation of Nonlinear Special Functions , 2015 .