Generalized Rational Krylov Decompositions with an Application to Rational Approximation

نویسندگان

  • Mario Berljafa
  • Stefan Güttel
چکیده

Generalized rational Krylov decompositions are matrix relations which, under certain conditions, are associated with rational Krylov spaces. We study the algebraic properties of such decompositions and present an implicit Q theorem for rational Krylov spaces. Transformations on rational Krylov decompositions allow for changing the poles of a rational Krylov space without recomputation, and two algorithms are presented for this task. Using such transformations we develop a rational Krylov method for rational least squares fitting. Numerical experiments indicate that the proposed method converges fast and robustly. A MATLAB toolbox with implementations of the presented algorithms and experiments is provided.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Nleigs: a Class of Robust Fully Rational Krylov Methods for Nonlinear Eigenvalue Problems∗

A new rational Krylov method for the efficient solution of nonlinear eigenvalue problems, A(λ)x = 0, is proposed. This iterative method, called fully rational Krylov method for nonlinear eigenvalue problems (abbreviated as NLEIGS), is based on linear rational interpolation and generalizes the Newton rational Krylov method proposed in [R. Van Beeumen, K. Meerbergen, and W. Michiels, SIAM J. Sci....

متن کامل

Rational Krylov approximation of matrix functions: Numerical methods and optimal pole selection∗

Matrix functions are a central topic of linear algebra, and problems of their numerical approximation appear increasingly often in scientific computing. We review various rational Krylov methods for the computation of large-scale matrix functions. Emphasis is put on the rational Arnoldi method and variants thereof, namely, the extended Krylov subspace method and the shift-and-invert Arnoldi met...

متن کامل

Computing Approximate (Symmetric Block) Rational Krylov Subspaces without Explicit Inversion

It has been shown that approximate extended Krylov subspaces can be computed –under certain assumptions– without any explicit inversion or system solves. Instead the necessary products A−1v are obtained in an implicit way retrieved from an enlarged Krylov subspace. In this paper this approach is generalized to rational Krylov subspaces, which contain besides poles at infinite and zero also fini...

متن کامل

Application of Vector-Valued Rational Approximations to the Matrix Eigenvalue Problem and Connections with Krylov Subspace Methods

Let F(z) be a vector-valued function F C CN, which is analytic at z 0 and meromorphic in a neighborhood of z 0, and let its Maclaurin series be given. In a recent work [J. Approx. Theory, 76 (1994), pp. 89-111] by the author, vector-valued rational approximation procedures for F(z) that are based on its Maclaurin series, were developed, and some of their convergence properties were analyzed in ...

متن کامل

Superlinear convergence of the rational Arnoldi method for the approximation of matrix functions

A superlinear convergence bound for rational Arnoldi approximations to functions of matrices is derived. This bound generalizes the well-known superlinear convergence bound for the CG method to more general functions with finite singularities and to rational Krylov spaces. A constrained equilibrium problem from potential theory is used to characterize a max-min quotient of a nodal rational func...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • SIAM J. Matrix Analysis Applications

دوره 36  شماره 

صفحات  -

تاریخ انتشار 2015