English

Maximum of the resolvent over matrices with given spectrum

Numerical Analysis 2018-02-27 v2 Functional Analysis Spectral Theory

Abstract

In numerical analysis it is often necessary to estimate the condition number CN(T)=TT1CN(T)=||T||_{} \cdot||T^{-1}||_{} and the norm of the resolvent (ζT)1||(\zeta-T)^{-1}||_{} of a given n×nn\times n matrix TT. We derive new spectral estimates for these quantities and compute explicit matrices that achieve our bounds. We recover the well-known fact that the supremum of CN(T)CN(T) over all matrices with T1||T||_{} \leq1 and minimal absolute eigenvalue r=mini=1,...,nλi>0r=\min_{i=1,...,n}|\lambda_{i}|>0 is the Kronecker bound 1rn\frac{1}{r^{n}}. This result is subsequently generalized by computing the corresponding supremum of (ζT)1||(\zeta-T)^{-1}||_{} for any ζ1|\zeta| \leq1. We find that the supremum is attained by a triangular Toeplitz matrix. This provides a simple class of structured matrices on which condition numbers and resolvent norm bounds can be studied numerically. The occuring Toeplitz matrices are so-called model matrices, i.e. matrix representations of the compressed backward shift operator on the Hardy space H2H_2 to a finite-dimensional invariant subspace.

Keywords

Cite

@article{arxiv.1501.07007,
  title  = {Maximum of the resolvent over matrices with given spectrum},
  author = {Oleg Szehr and Rachid Zarouf},
  journal= {arXiv preprint arXiv:1501.07007},
  year   = {2018}
}