H2-optimal model order reduction using hyperbolic geometry
Abstract
We revisit the problem of approximating a linear time-invariant LTI model by a reduced-order model (ROM) in an H2-optimal sense. A widely used method for computing such ROMs is the iterative rational Krylov algorithm (IRKA). However, in practice, the stability of the resulting ROM is not generally guaranteed. In this work, we develop a novel framework that can be used to obtain H2-optimal models of single-input single-output LTI systems while enforcing stability of the ROM. We use Malmquist-Takenaka orthogonal expansions in the corresponding H2-spaces and their connection to hyperbolic geometry to develop an equilibrium function. We show that the poles of an H2-optimal ROM transfer function are defined by the zeros of this equilibrium function. The resulting nonlinear root-finding problem is solved using a Newton-type method formulated on the products of Poincaré disk manifolds, where updates are performed along hyperbolic geodesics ensuring that the stability constraint is satisfied at every iteration. We present numerical experiments on benchmark systems for both discrete- and continuous-time LTI MOR problems, including examples where IRKA produces unstable reduced models while the proposed method preserves stability throughout the iteration.
Related
- Exponential Turnpike in Optimal Boundary Control of Multi-Dimensional First-Order Hyperbolic Systems
- Truncated ℋ2 Optimal Model Order Reduction for Bilinear Quadratic Output Systems
- Porteous Meets Walsh: Counting Critical Points for ℋ2-Norm Model Order Reduction Using Algebraic Geometry Behind Multiparameter Eigenvalue Problems
- Norm-optimal model order reduction in reproducing kernel Hilbert spaces
- Hyperbolic periodic orbits and optimal 1/t convergence of the Lax--Oleinik semigroup
BibTeX
@misc{angino2026h2optimal,
title = {{$H_2$-optimal model order reduction using hyperbolic geometry}},
author = {Andrea Angino and Tamás Dózsa and Matthias Voigt},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.16249},
eprint = {2610.16249},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.16249},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - H_2-optimal model order reduction using hyperbolic geometry
AU - Andrea Angino
AU - Tamás Dózsa
AU - Matthias Voigt
PY - 2026
DA - 2026/10/12
PB - arXiv
DO - 10.48550/arxiv.2610.16249
UR - https://arxiv.org/abs/2610.16249
AB - We revisit the problem of approximating a linear time-invariant LTI model by a reduced-order model (ROM) in an H_2-optimal sense. A widely used method for computing such ROMs is the iterative rational Krylov algorithm (IRKA). However, in practice, the stability of the resulting ROM is not generally guaranteed. In this work, we develop a novel framework that can be used to obtain H_2-optimal models of single-input single-output LTI systems while enforcing stability of the ROM. We use Malmquist-Takenaka orthogonal expansions in the corresponding H_2-spaces and their connection to hyperbolic geometry to develop an equilibrium function. We show that the poles of an H_2-optimal ROM transfer function are defined by the zeros of this equilibrium function. The resulting nonlinear root-finding problem is solved using a Newton-type method formulated on the products of Poincaré disk manifolds, where updates are performed along hyperbolic geodesics ensuring that the stability constraint is satisfied at every iteration. We present numerical experiments on benchmark systems for both discrete- and continuous-time LTI MOR problems, including examples where IRKA produces unstable reduced models while the proposed method preserves stability throughout the iteration.
ER -