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.

Stability-Preserving Rational Approximation Subject to Interpolation Constraints

Stable Rational Approximation in the Context of Interpolation and Convex Optimization

On Degree-Constrained Analytic Interpolation with Interpolation Points Close to the Boundary.

The Inverse Problem of Analytic Interpolation with Degree Constraint

The Inverse Problem of Analytic Interpolation with Degree Constraint and Weight Selection for Control Synthesis

Analytic interpolation with degree constraint: a constructive theory with applications to control and signal processing

Least Squares Rational Approximation

Passivity-preserving Model Reduction by Analytic Interpolation

Multivariable Analytic Interpolation with Complexity Constraints: A Modified Riccati Approach

Rational Approximation of Some Limited Smooth Functions

On Analytic Interpolation with Non-Classical Constraints for Solving Problems in Robust Control.

Rational Empirical Interpolation Methods with Applications

CONSTRUCTIVE THEORY OF MULTIVARIATE RATIONAL INTERPOLATION OF DEGREE ZERO WITH ALGORITHM

The Linear Pencil Approach to Rational Interpolation

An algorithm for two‐variable rational interpolation suitable for matrix manipulations with the evaluation–interpolation method

Interpolation via weighted $l_1$ minimization

Stabilized interpolation using radial basis functions augmented with selected radial polynomials

On Exact and Approximate Interpolation of Sparse Rational Functions

A Modified Riccati Approach to Analytic Interpolation with Applications to System Identification and Robust Control

On the convergence of polynomial approximation of rational functions

On Interpolation Approximation: Convergence Rates for Polynomial Interpolation for Functions of Limited Regularity.