Bernhard Beckermann
The Condition Number of real Vandermonde, Krylov and positive definite Hankel matrices
Key words : Condition number, Vandermonde matrix, Krylov matrix, Hankel matrix, Hilbert matrix.
Classifications: AMS(MOS): 15A12, 65F35; CR: G1.3, G1.6.

Abstract

We show that the Euclidean condition number of any positive definite Hankel matrix of order $ngeq 3$ may be bounded from below by $gamma^{n-1}/(16n)$ with $gamma=exp(4 cdot Catalan/pi) approx 3.210$, and that this bound may be improved at most by a factor $8 gamma n$. Similar estimates are given for the class of real Vandermonde matrices, the class of row scaled real Vandermonde matrices, and the class of Krylov matrices with Hermitian argument. Improved bounds are derived for the case where the abscissa or eigenvalues are included in a given real interval. Our findings confirm that all such matrices --- including for instance the famous Hilbert matrix --- are ill--conditioned already for ``moderate'' order.

As application, we describe implications of our results for the numerical condition of various tasks in Numerical Analysis such as polynomial and rational interpolation at real nodes, determination of real roots of polynomials, computation of coefficients of orthogonal polynomials, or the iterative solution of linear Hermitian systems of equations.