Abstract:We establish a new perturbation theory for orthogonal polynomials using a Riemann--Hilbert approach and consider applications in numerical linear algebra and random matrix theory. This new approach shows that the orthogonal polynomials with respect to two measures can be effectively compared using the difference of their Stieltjes transforms on a suitably chosen contour. Moreover, when two measures are close and satisfy some regularity conditions, we use the theta functions of a hyperelliptic Riemann surface to derive explicit and accurate expansion formulae for the perturbed orthogonal polynomials. In contrast to other approaches, a key strength of the methodology is that estimates can remain valid as the degree of the polynomial grows. The results are applied to analyze several numerical algorithms from linear algebra, including the Lanczos tridiagonalization procedure, the Cholesky factorization and the conjugate gradient algorithm. As a case study, we investigate these algorithms applied to a general spiked sample covariance matrix model by considering the eigenvector empirical spectral distribution and its limits. For the first time, we give precise estimates on the output of the algorithms, applied to this wide class of random matrices, as the number of iterations diverges. In this setting, beyond the first order expansion, we also derive a new mesoscopic central limit theorem for the associated orthogonal polynomials and other quantities relevant to numerical algorithms.

Blockwise Perturbation Theory for Block P-Cyclic Stochastic Matrices

Blockwise perturbation theory for nearly uncoupled Markov chains and its application

Blockwise Perturbation Theory for Markov Chains

The Local Stochastic Stability for Complex Networks under Pinning Control.

Perturbation and Inverse Problems of Stochastic Matrices

Probabilistic perturbation bounds of matrix decompositions

Invariant subspace perturbation of a matrix with Jordan blocks

A stochastic perturbation analysis of the QR decomposition and its applications

A Riemann--Hilbert approach to the perturbation theory for orthogonal polynomials: Applications to numerical linear algebra and random matrix theory

A note on entrywise perturbation theory for Markov chains

Entrywise perturbation theory for rate matrices of GI/M/1 type Markov chains

A Note on Perturbation Estimations for Spectral Projectors

Perturbative power series for block diagonalisation of Hermitian matrices

Perturbation analysis of the matrix equation X - \sum_{i=1}^m A_i^* X^{p_i} A_i = Q

Solutions and improved perturbation analysis for the matrix equation X-A^{*}X^{-p}A=Q (p>0)

Perturbation of the Drazin Inverse for Matrices with Equal Eigenprojections at Zero

Limit distribution of partial transposition of block random matrices

Entrywise relative perturbation bounds for exponentials of essentially non-negative matrices

On Kemeny's constant and stochastic complement

Application of the Schur-Cohn Theorem to the precise convergence domain for a p-cyclic SOR iteration matrix

Random $p$-adic matrices with fixed zero entries and the Cohen--Lenstra distribution