$\mathcal{H}_2$-optimal model reduction of linear quadratic-output systems by multivariate rational interpolation
By: Sean Reiter , Ion Victor Gosea , Igor Pontes Duff and more
Potential Business Impact:
Makes complex systems simpler for computers.
This paper addresses the $\mathcal{H}_2$-optimal approximation of linear dynamical systems with quadratic-output functions, also known as linear quadratic-output systems. Our major contributions are threefold. First, we derive interpolation-based first-order optimality conditions for the linear quadratic-output $\mathcal{H}_2$ minimization problem. These conditions correspond to the mixed-multipoint tangential interpolation of the full-order linear- and quadratic-output transfer functions, and generalize the Meier-Luenberger optimality framework for the $\mathcal{H}_2$-optimal model reduction of linear time-invariant systems. Second, given the interpolation data, we show how to enforce these mixed-multipoint tangential interpolation conditions explicitly by Petrov-Galerkin projection of the full-order model matrices. Third, to find the optimal interpolation data, we build on this projection framework and propose a generalization of the iterative rational Krylov algorithm for the $\mathcal{H}_2$-optimal model reduction of linear quadratic-output systems, called LQO-IRKA. Upon convergence, LQO-IRKA produces a reduced linear quadratic-output system that satisfies the interpolatory optimality conditions. The method only requires solving shifted linear systems and matrix-vector products, thus making it suitable for large-scale problems. Numerical examples are included to illustrate the effectiveness of the proposed method.
Similar Papers
First and Second Order Optimal $\mathcal{H}_2$ Model Reduction for Linear Continuous-Time Systems
Optimization and Control
Makes computer models simpler and more accurate.
An iterative tangential interpolation framework for model reduction of MIMO systems
Systems and Control
Makes big computer models smaller and faster.
Interpretable and flexible non-intrusive reduced-order models using reproducing kernel Hilbert spaces
Computational Engineering, Finance, and Science
Makes computer models faster and easier to understand.