Constructing Sobolev orthonormal rational functions via an updating procedure

Abstract

In this paper, we generate the recursion coefficients for rational functions with prescribed poles that are orthonormal with respect to a continuous Sobolev inner product. Using a rational Gauss quadrature rule, the inner product can be discretized, thus allowing a linear algebraic approach. The presented approach involves reformulating the problem as an inverse eigenvalue problem involving a Hessenberg pencil, where the pencil will contain the recursion coefficients that generate the sequence of Sobolev orthogonal rational functions. This reformulation is based on the connection between Sobolev orthonormal rational functions and the orthonormal bases for rational Krylov subspaces generated by a Jordan-like matrix. An updating procedure, introducing the nodes of the inner product one after the other, is proposed and the performance is examined through some numerical examples.