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.

The Covariance Extension Equation: A Riccati-Type Approach to Analytic Interpolation

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

Multivariable Analytic Interpolation with Complexity Constraints: A Modified Riccati Approach

The Covariance Extension Equation Revisited

A Homotopy Continuation Solution of the Covariance Extension Equation

A Convex Optimization Approach to the Rational Covariance Extension Problem

On the Convergence of the Accelerated Riccati Iteration Method.

On the well-posedness of the rational covariance extension problem

Using covariance extension equation to solve the Nevanlinna-Pick interpolation with degree constraint

A Convex Optimization Approach to the Covariance Extension Problem with Degree Constraint

The Geometry of Positive Real Functions with Applications to the Rational Covariance Extension Problem

A complete parameterization of all positive rational extensions of a covariance sequence

The Circulant Rational Covariance Extension Problem: the Complete Solution

On The Multivariate Circulant Rational Covariance Extension Problem

Multivariable simultaneous stabilization: A modified Riccati approach

From Finite Covariance Windows to Modeling Filters: A Convex Optimization Approach

A general solution to the simultaneous stabilization problem by analytic interpolation

Least Squares Rational Approximation

Rational Empirical Interpolation Methods with Applications

Homotopy for Rational Riccati Equations Arising in Stochastic Optimal Control.

Matrix-valued Nevanlinna-Pick Interpolation with Complexity Constraint: an Optimization Approach