$\mathcal{H}_2$-Optimal Model Reduction Using Projected Nonlinear Least Squares

Abstract

In many applications throughout science and engineering, model reduction plays an important role replacing expensive large-scale linear dynamical systems by inexpensive reduced order models that capture key features of the original, full order model. One approach to model reduction is to find reduced order models that are locally optimal approximations in the H 2 norm, an approach taken by the Iterative Rational Krylov Algorithm (IRKA) and several others. Here we introduce a new approach for H 2 -optimal model reduction using the projected nonlinear least squares framework previously introduced in [J. M. Hokanson, SIAM J. Sci. Comput. 39 (2017), pp. A3107-A3128]. At each iteration, we project the H 2 optimization problem onto a finite-dimensional subspace yielding a weighted least rational approximation problem. Subsequent iterations append this subspace such that the least squares rational approximant asymptotically satisfies the first order necessary conditions of the original, H 2 optimization problem. This enables us to build reduced order models with similar error in the H 2 norm as competing methods but using far fewer evaluations of the expensive, full order model. Moreover, our new algorithm only requires access to the transfer function of the full order model, unlike IRKA which requires a state-space representation or TF-IRKA which requires both the transfer function and its derivative. This application of projected nonlinear least squares to the H 2 -optimal model reduction problem suggests extensions of this approach related model reduction problems.