Abstract:Rational approximation appears in many contexts throughout science and engineering, playing a central role in linear systems theory, special function approximation, and many others. There are many existing methods for solving the rational approximation problem, from fixed point methods like the Sanathanan-Koerner iteration and Vector Fitting, to partial interpolation methods like Adaptive Anderson Antoulas (AAA). While these methods can often find rational approximations with a small residual norm, they are unable to find optimizers with respect to a weighted l2 norm with a square dense weighting matrix. Here we develop a nonlinear least squares approach constructing rational approximations with respect to this norm. We explore this approach using two parameterizations of rational functions: a ratio of two polynomials and a partial fraction expansion. In both cases, we show how we can use Variable Projection (VARPRO) to reduce the dimension of the optimization problem. As many applications seek a real rational approximation that can be described as a ratio of two real polynomials, we show how this constraint can be enforced in both parameterizations. Although this nonlinear least squares approach often converge to suboptimal local minimizers, we find this can be largely mitigated by initializing the algorithm using the poles of the AAA algorithm applied to the same data. This combination of initialization and nonlinear least squares enables us to construct rational approximants using dense and potentially ill-conditioned weight matrices such as those that appear as a step in new H2 model reduction algorithm recently developed by the authors.

n-th Root Optimal Rational Approximants to Functions with Polar Singular Set

Multivariate rational approximation of functions with curves of singularities

Exponential node clustering at singularities for rational approximation, quadrature, and PDEs

A Geometric Approach to Polynomial and Rational Approximation

Resolution of singularities by rational functions

Riemann-Hilbert Characterisation of Rational Functions with a General Distribution of Poles on the Extended Real Line Orthogonal with Respect to Varying Exponential Weights: Multi-Point Padé Approximants and Asymptotics

Rational function approximation with normalized positive denominators

Unitary rational best approximations to the exponential function

Adelic perturbation of rational functions and applications

Rational minimax approximation via adaptive barycentric representations

Effective irrationality measures for real and $p$-adic roots of rational numbers close to $1$, with an application to parametric families of Thue-Mahler equations

A Rational Approximant for the Digamma Function

Meromorphic Functions with a Polar Asymptotic Value

Self-similar continued root approximants

optimal rational approximation on general domains

Representation of conformal maps by rational functions

Least Squares Rational Approximation

Quantitatively hyper-positive real rational functions II

Computation of polynomial and rational approximations in complex domains by the -method

Rational inner functions on the polydisk -- a survey

On sets of rational functions which locally represent all of $\mathbb{Q}$