Adaptive Rational Interpolation: Arnoldi and Lanczos-like Equations

Abstract

The Arnoldi and Lanczos algorithms, which belong to the class of Krylov subspace methods, are increasingly used for model reduction of large-scale systems. The standard versions of the algorithms tend to create reduced order models that poorly approximate low frequency dynamics. Rational Arnoldi and Lanczos algorithms produce reduced models that approximate dynamics at various frequencies. This paper tackles the issue of developing simple Arnoldi and Lanczos equations for the rational case. This allows a simple error analysis to be carried out for both algorithms and permits the development of computationally efficient model reduction algorithms, where the frequencies at which the dynamics are to be matched can be updated adaptively.