We show that the well-known Levinson algorithm for computing the inverse Cholesky factorization of positivedefinite Toeplitz matrices can be viewed as a special case of a more general process. The latter processprovides a very efficient implementation of the Arnoldi process when the underlying operator is isometric.This is analogous with the case of Hermitian operators where the Hessenberg matrix becomes tridiagonal andresults in the Hermitian Lanczos process. We investigate the structure of the Hessenberg matrices in theisometric case and show that simple modifications of them move all their eigenvalues to the unit circle. Theseeigenvalues are then interpreted as abscissas for analogs of Gaussian quadrature, now on the unit circleinstead of the real line. The trapezoidal rule appears as the analog of the Gauss-Legendre formula.
展开▼