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.

Perturbation Analysis for Best Approximation and the Polar Factor by Subunitary Matrices

The Equivalence Relations Between the Ditzian–Totik Moduli of Smoothness and Best Polynomial Approximation in the Besov Spaces

Fe b 20 10 PARTIAL DIFFERENCE SETS FROM pARY WEAKLY REGULAR BENT FUNCTIONS AND QUADRATIC FORMS

Perturbation Analysis for Subspace Decomposition with Applications in Subspace-Based Algorithms

Componentwise perturbation analysis for the generalized Schurdecomposition

Polar Coding for Arbitrary Channels Based on Gaussian Decomposition

Generalized Gaussian Approximation and Its Application in Polar Decoding.

Matrix perturbation analysis of methods for extracting singular values from approximate singular subspaces

LIL and the Approximation of Rectangular Sums of B-valued Random Variables when Extreme Terms Are Excluded

Probabilistic perturbation bounds of matrix decompositions

General Algorithms for the Polar Decomposition and Strains

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

PERTURBATION THEORY FOR THE LU AND QR FACTORIZATIONS

The polarization hierarchy for polynomial optimization over convex bodies, with applications to nonnegative matrix rank

A stochastic perturbation analysis of the QR decomposition and its applications

Making matrices better: Geometry and topology of polar and singular value decomposition

Perturbation Analysis for Qr Factorization of Unitary-Symmetric Matrix

n-th Root Optimal Rational Approximants to Functions with Polar Singular Set

Approximation and quadrature by weighted least squares polynomials on the sphere

Bounds on Polarization Problems on Compact Sets via Mixed Integer Programming

New rigorous perturbation bounds for the LU and QR factorizations