Pith. sign in

REVIEW 1 major objections 6 minor 2 cited by

$\mathcal{H}_2$-optimal model reduction of linear quadratic-output systems by multivariate rational interpolation

T0 review · 1 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash

Pith's one-line read An $\mathcal{H}_2$-optimal reduced linear quadratic-output model must tangentially interpolate the full system's linear and quadratic transfer functions at the mirror images of its own poles.

desk verdict Genuinely extends H2-optimal interpolation theory to LQO systems, but the main proof has a repairable yet central class-membership gap. read the letter →

arxiv 2505.03057 v2 pith:W3C4CKLD submitted 2025-05-05 math.NA cs.NAcs.SYeess.SYmath.DSmath.OC

classification math.NAcs.NAcs.SYeess.SYmath.DSmath.OC MSC 34C2041A0549K1565J0565F9993A1593C1093C80
keywords modelreductionH2-optimalitylinearquadratic-outputsystemstangentialinterpolationmultivariaterationaliterativeKrylovalgorithmPetrov-GalerkinprojectionVolterrakernels
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper derives first-order necessary conditions for a reduced model of a linear quadratic-output (LQO) system, meaning linear state dynamics with an output made of a linear term plus a quadratic term $M(x\otimes x)$, to minimize the $\mathcal{H}_2$ approximation error. The conditions say that at the mirror images of the reduced model's poles, the reduced linear and quadratic transfer functions must tangentially match the full-order ones, and certain weighted sums of the two must be matched in both value and derivative. These mixed-multipoint tangential interpolation conditions generalize the classical Meier-Luenberger conditions for linear systems. The paper also shows how to enforce all conditions by Petrov-Galerkin projection and gives an iterative rational Krylov algorithm, LQO-IRKA, that converges to models satisfying them. If correct, this turns $\mathcal{H}_2$-optimal quadratic-output model reduction into a multivariate rational interpolation problem that can be solved with shifted linear solves.

What carries the argument

The load-bearing object is the pole-residue expansion of the reduced transfer functions, $\tilde G_1(s)=\sum_j c_j b_j^T/(s-\lambda_j)$ and $\tilde G_2(s_1,s_2)=\sum_{j,k} m_{j,k}(b_j\otimes b_k)^T/((s_1-\lambda_j)(s_2-\lambda_k))$. Theorem 2.1 uses these expansions to express the $\mathcal{H}_2$ inner product and norm of an LQO system as finite evaluations of $G_1$ and $G_2$ at the mirrored poles $-\lambda_j$. Lemma 2.1's symmetry identities for $G_2$, inherited from the commutation matrix, let the authors group the left-tangential and derivative conditions into the compact form (29c)-(29d). Theorem 3.2 then constructs projection bases $V,W$ whose columns are the shifted solves of (31) and (32), so that Petrov-Galerkin projection enforces all $3r+r^2$ interpolation conditions simultaneously.

What would settle it

One could settle Theorem 3.1 by computing, for a small LQO example, a global $\mathcal{H}_2$-optimal order-$r$ reduced model with simple poles and evaluating whether (29a)-(29d) fail at its mirrored poles; any violation would refute the necessity claim. Alternatively, a stable full-order LQO system whose LQO-IRKA iterates converge to a model that still violates one of the four conditions would show the algorithm does not deliver the claimed optimality certificate.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is Theorem 3.1: if an asymptotically stable order-$r$ LQO system $\tilde G$ with simple poles minimizes the squared $\mathcal{H}_2$ error against a full-order LQO system $G$, then the transfer functions $\tilde G_1$ and $\tilde G_2$ satisfy the four tangential interpolation conditions (29a)-(29d). In particular, $\tilde G_1$ and $\tilde G_2$ individually right-interpolate $G_1$ and $G_2$ at $-\lambda_k$ along residue directions, and linear combinations of the two transfer functions are interpolated in the left-tangential Lagrange and bi-tangential Hermite senses at all pairs of mirrored poles. This gives the LQO analogue of the Meier-Luenberger characterization and, because the proof only uses Hardy-space membership of the full-order functions, the paper notes the conditions extend beyond the specific LQO realization.

Load-bearing premise

The argument assumes the $\mathcal{H}_2$-optimal reduced model can be chosen with simple poles and that small pole/residue perturbations stay within the class of order-$r$ asymptotically stable LQO systems; the paper does not prove that a minimizer with simple poles always exists, and LQO-IRKA's convergence is reported empirically rather than proved.

Editorial extensions

If this is right

  • Every $\mathcal{H}_2$-optimal reduced LQO model with simple poles is a tangential interpolant, so the quest for optimal reduced models can be reposed as choosing poles and residue directions that satisfy (29a)-(29d).
  • Petrov-Galerkin projection with the bases (31) and (32) enforces all the optimality conditions at once, guaranteeing that any converged LQO-IRKA iterate meets the first-order necessary conditions.
  • LQO-IRKA requires only shifted linear system solves and matrix-vector products, making the framework applicable to large-scale systems where balancing-based LQO methods need costly Lyapunov solves.
  • When the quadratic output term is zero, conditions (29) reduce to the classical linear Meier-Luenberger interpolation conditions, so the paper's results contain the linear theory as a special case.
  • Because the $\mathcal{H}_2$ error bounds the time-domain $L^\infty$ output error via (20), satisfying these optimality conditions also controls worst-case output accuracy over time for finite-energy inputs.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The proof of Theorem 3.1 does not require the full-order model to have LQO structure, so the same interpolation conditions should characterize $\mathcal{H}_2$-optimal rational approximants for systems with second-order or delay dynamics; whether projection can construct them is left open by the paper.
  • No convergence proof for LQO-IRKA is given; a natural test is whether iterates always reach a stationary point from arbitrary initial data, and whether the simple-pole restriction ever excludes the true minimizer.
  • The mixed conditions couple every pair of poles through the $m_{k,\ell}$ residues, so a large-scale implementation may need to exploit sparsity or low-rank structure in the sums over $\ell$; the paper does not address this cost.
  • A data-driven analogue could identify the interpolation data directly from frequency samples of $G_1$ and $G_2$, bypassing the state-space realization entirely.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 6 minor

Summary. This manuscript studies H2-optimal model reduction for linear systems with quadratic outputs (LQO systems). The main theoretical result, Theorem 3.1, states that any asymptotically stable order-r LQO reduced model with simple poles that minimizes the H2 error (22) must satisfy the four groups of tangential interpolation conditions (29a)-(29d), generalizing the Meier-Luenberger conditions from LTI systems. Theorem 3.2 gives explicit Petrov-Galerkin bases V and W that enforce these interpolation conditions at prescribed data, and Algorithm 4.1 (LQO-IRKA) iteratively updates the interpolation data using the current reduced model's poles and residue directions. The paper includes a detailed proof of the pole-residue H2 inner product formula (Theorem 2.1), a real-valued basis construction (Lemma 4.1), and numerical experiments on a 1D advection-diffusion problem with a quadratic cost, with code and data released on Zenodo.

Significance. The interpolatory optimality conditions are a natural and non-obvious generalization of the classical H2 interpolation theory, and the projection-based enforcement in Theorem 3.2 gives a constructive route to reduced models satisfying them. The numerical results suggest LQO-IRKA is competitive with balanced truncation for LQO systems and robust to initialization. Strengths include the detailed proofs of Theorem 2.1 and Theorem 3.2, the careful treatment of real-valued bases, the explicit complexity statement (shifted linear solves only), and the availability of reproducible code. The main caveat is the gap in the proof of Theorem 3.1 discussed below; if repaired, the paper would be a solid contribution to the model reduction literature.

major comments (1)
  1. [Appendix B / Theorem 3.1] The proof of Theorem 3.1 does not establish that the perturbed systems used to derive the necessary conditions lie in the admissible class of real order-r LQO systems. For (29b), equation (53) perturbs only the (j,k)-th quadratic residue direction m_{j,k}; when j≠k this violates the symmetry constraint m_{j,k}=m_{k,j} of equation (25), so ˇG is not of the form (2) and the sub-optimality inequality (30) cannot be invoked. For (29c), equation (54) perturbs only b_k, and for (29d), equation (58) shifts only λ_k; in both cases, for complex-conjugate pole pairs the necessary conjugate update of b_{k+1} or λ_{k+1} is omitted, so the perturbed system is complex-valued rather than a real admissible LQO system. These are not merely formalities: the contradiction argument proves stationarity over a larger unconstrained pole-residue class, which does not imply stationarity over the constrained LQO set. The theorem is likely repairable by using symmetric two-index perturbations for (29b) and real-structure-preserving perturbations for (29c)-(29d), but the proof of the central claim must be revised.
minor comments (6)
  1. [Proof of Theorem 3.2] In the sentence following equation (37), "to prove (33b)" should read "to prove (33c)", since the argument establishes the left-tangential Lagrange condition.
  2. [Lemma 4.1] The block matrix Q stated as (1/√2)[1 1; i −i] does not correctly relate W and W_p; the transformation should be (1/2)[1 −i; 1 i] up to column scaling, and "orthogonal" should be "unitary".
  3. [Section 5.2] The numerical computation of the relative H2 error in (47) is not described; the authors should state whether a full eigendecomposition of the order-3000 system or a Lyapunov solve was used.
  4. [Section 5.3] Only a single benchmark is tested, and no comparison is made with the Wilson/gramian-based H2-optimal method of [33] or the Riemannian BFGS method of [44]; such a comparison would strengthen the practical claims.
  5. [Section 2.3 / Eq. (17)] The notation \(\overline{G_1(s)}\) for conjugation of coefficients only is confusing; suggest defining it explicitly or using a different symbol.
  6. [Section 4.2.2] The absence of a convergence proof for LQO-IRKA is acknowledged, but the authors should state explicitly that convergence and attainment of (29) are empirical observations rather than guaranteed properties.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: Theorem 3.1's interpolatory optimality conditions are obtained by an independent variational perturbation argument; self-citations are contextual and non-load-bearing.

full rationale

The paper's central claim, Theorem 3.1, is not circular. The mixed-multipoint tangential interpolation conditions (29a)-(29d) are derived in Appendix B from the sub-optimality inequality (30) by explicit first-order perturbations of the reduced-model poles and residue directions, using the H2 inner-product and norm formulae (26)-(27) that are proved independently in Theorem 2.1 and Appendix A. The optimality conditions are not assumed as inputs; they are consequences of stationarity of the squared H2 error. The paper explicitly recognizes the apparent circularity of needing the optimal reduced model to set the interpolation data ('Of course, this assumes having access to the optimal reduced model. We resolve this circular causality issue in the section.') and resolves it in Section 4 with the LQO-IRKA fixed-point iteration, which updates interpolation data from the previous iterate. This is the standard IRKA self-consistency device, not a circular derivation: convergence to a fixed point satisfies independently derived necessary conditions. The projection construction in Theorem 3.2 is explicit and self-contained; the citation to the authors' prior work ([33, Theorem 3.2], 'any H2-optimal approximation ... is necessarily obtained via a Petrov-Galerkin projection') is contextual and not load-bearing, since Theorem 3.2 proves directly how to enforce the interpolation conditions by constructing V and W. Self-citations to [9,22,23,24] provide background and comparison algorithms, not the theorem's proof. Numerical benchmarking against the full-order time-domain output and against LQO-BT is external evidence, not a fit. The manuscript does flag genuine limitations: convergence of LQO-IRKA is reported empirically rather than proved (Section 4.2.2), and asymptotic stability is not guaranteed by the iteration. In addition, the Appendix B perturbation proof has a class-membership gap: perturbing a single residue m_{j,k} in (53), or applying a single complex pole shift in (58), may leave the class of real-valued LQO systems satisfying the symmetry constraint (25), so inequality (30) is invoked for systems not shown to be admissible. This is a correctness and rigor issue, not a circularity issue: the claimed conditions are not equivalent to the proof's inputs by construction. No step in the derivation reduces to a fitted parameter renamed as a prediction, and no uniqueness theorem is imported from the authors' prior work to force the choice.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No free parameters or invented entities appear in the derivation. The theory rests on standard complex analysis and on domain assumptions about stability, real-valued realizations, and simple poles of the reduced model. The algorithmic inputs (initial interpolation data, tolerance, order r) are user choices rather than fitted parameters.

assumptions (5)
  • domain assumption The full-order system (1) is asymptotically stable and E is nonsingular
    Assumed at the start of Section 1; required for the transfer functions G1 and G2 to lie in the Hardy spaces used to define the H2 norm.
  • domain assumption The reduced model has simple poles
    Assumed in Section 2.4 and Theorem 3.1; enables the pole-residue expansion (23) and the residue-based H2 inner product formula in Theorem 2.1. The paper notes repeated poles rarely appear in practice but does not prove the optimal minimizer has simple poles.
  • domain assumption Real-valued state-space realization, so G1(s)=G1(s) and G2(s1,s2)=G2(s1,s2)
    Used in the perturbation proofs in Appendix B to evaluate inner products via Theorem 2.1; this is standard for real dynamical systems.
  • standard math Residue theorem and contour integration results apply to the bivariate transfer functions
    The proof of Theorem 2.1 in Appendix A uses ML-estimates and the residue theorem to evaluate the double integrals; these are standard results in complex analysis.
  • domain assumption Perturbations of the reduced model parameters produce valid order-r LQO systems
    The proof of Theorem 3.1 perturbs poles and residue directions and implicitly assumes the resulting transfer functions correspond to some order-r LQO system; this holds because any function of the form (23) with the symmetry (25) is realizable, but the paper does not explicitly state this.

how reviews work

0 comments
Cite this review

Pith. "Pith review of $\mathcal{H}_2$-optimal model reduction of linear quadratic-output systems by multivariate rational interpolation." pith.science (2026). https://pith.science/paper/W3C4CKLD

@misc{pith2026250503057,
  author       = {Pith},
  title        = {Pith review of: $\mathcalH_2$-optimal model reduction of linear quadratic-output systems by multivariate rational interpolation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/W3C4CKLD}},
  note         = {Machine review of arXiv:2505.03057}
}
abstract

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 interpolatory 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 optimal interpolation data, we show how to enforce the interpolatory optimality conditions explicitly by Petrov-Galerkin projection of the full-order model. 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 reduced linear quadratic-output systems that satisfy 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.

Figures

Figures reproduced from arXiv: 2505.03057 by the authors.

Figure 1
Figure 1. Output magnitudes and pointwise relative errors [PITH_FULL_IMAGE:figures/full_fig_p024_1.png] view at source ↗
Figure 2
Figure 2. Relative H2 errors of the intermediate reduced models computed by LQO-IRKAeigs and LQO-IRKAimag for the first 50 iterations. LQO-IRKAeigs LQO-IRKAimag LQO-BT interponeStep,eigs interponeStep,imag Run time (s) 58.23 s 56.31 s 72.61 s 0.41 s 0.44 s Iteration count 124 110 N/A N/A N/A [PITH_FULL_IMAGE:figures/full_fig_p025_2.png] view at source ↗
Figure 3
Figure 3. Relative H2 errors (47) due to the hierarchy of reduced models for orders r = 2, 4, . . . , 30. LQO-IRKAimag converge to the same local minimum for each order of reduction. 6. Conclusion We have presented a novel H2-optimality framework for the approximation of linear quadratic￾output systems (1) based on multivariate rational interpolation. In Theorem 3.1, we derive first￾order optimality conditions; these amount t… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. $\mathcal{H}_\infty$ model order reduction for quadratic output systems

    math.OC 2025-05 conditional novelty 7.0 of 10

    Introduces an H-infinity norm for linear systems with quadratic output and an optimization-based reduced-order modeling algorithm that minimizes it, with a structure-preserving variant for port-Hamiltonian systems.

  2. L2-L2-gain bounds for quadratic output systems

    math.OC 2026-07 unverdicted novelty 5.0 of 10

    Derives explicit L2-L2-gain bound for quadratic-output LTI systems; equals L2-norm of bivariate transfer function on anti-diagonal when output is purely state-quadratic, and is obtained by solving linear matrix equations.

Reference graph

Works this paper leans on

54 extracted references · 26 canonical work pages · cited by 2 Pith papers

  1. [33]

    H2 optimal model reduction of linear systems with multiple quadratic output s

    Sean Reiter, Igor Pontes Duff, Ion Victor Gosea, and Serka n Gugercin. H2 optimal model reduction of linear systems with multiple quadratic output s. e-prints 2405.05951, arXiv, 2024. URL: https://arxiv.org/abs/2405.05951. Preprint. May 7, 2025 S. Reiter, I. V. Gosea, I. Pontes Duff, S. Gugercin: H2-optimal MOR of LQO systems 42

  2. [44]

    H2 optimal model reduction of linear dy- namical systems with quadratic output by the Riemannian BFG S method

    Ping Yang, Zhao-Hong Wang, and Yao-Lin Jiang. H2 optimal model reduction of linear dy- namical systems with quadratic output by the Riemannian BFG S method. Mathematics and Computers in Simulation , 2025. doi:10.1016/j.matcom.2025.03.021

  3. [1]

    Antoulas

    Athanasios C. Antoulas. Approximation of Large-Scale Dynamical Systems . SIAM, Philadel- phia, PA, 2005. doi:10.1137/1.9780898718713

  4. [2]

    Antoulas, Christopher A

    Athanasios C. Antoulas, Christopher A. Beattie, and Ser kan G¨ u˘ gercin.Interpolatory Methods for Model Reduction . Computational Science & Engineering. SIAM, Philadelphia , PA, 2020. doi:10.1137/1.9781611976083

  5. [3]

    Quirin Aumann and Steffen W. R. Werner. Structured model or der reduction for vibro- acoustic problems using interpolation and balancing metho ds. J. Sound Vib. , 543:117363,

  6. [4]

    Energy-based approx imation of linear systems with polyno- mial outputs

    Linus Balicki and Serkan Gugercin. Energy-based approx imation of linear systems with polyno- mial outputs. e-prints 2409.19730, arXiv, 2024. URL: https://arxiv.org/abs/2409.19730

  7. [5]

    Model order r eduction for linear and nonlinear systems: a system-theoretic perspective

    Ulrike Baur, Peter Benner, and Lihong Feng. Model order r eduction for linear and nonlinear systems: a system-theoretic perspective. Archives of Computational Methods in Engineering , 21(4):331–358, 2014. doi:10.1007/s11831-014-9111-2

  8. [6]

    Inexact solves in interpo- latory model reduction

    Christopher Beattie, Serkan Gugercin, and Sarah Wyatt. Inexact solves in interpo- latory model reduction. Linear Algebra and its Applications , 436(8):2916–2943, 2012. doi:10.1016/j.laa.2011.07.015

Show all 54 references
  1. [7]

    Interpolation-based H2-model reduction of bilinear con- trol systems

    Peter Benner and Tobias Breiten. Interpolation-based H2-model reduction of bilinear con- trol systems. SIAM Journal on Matrix Analysis and Applications , 33(3):859–885, 2012. doi:10.1137/110836742. Preprint. May 7, 2025 S. Reiter, I. V. Gosea, I. Pontes Duff, S. Gugercin: H2-op...

  2. [8]

    H2-quasi-optimal model order reduction for quadratic-bilinear control systems

    Peter Benner, Pawan Goyal, and Serkan Gugercin. H2-quasi-optimal model order reduction for quadratic-bilinear control systems. SIAM Journal on Matrix Analysis and Applications , 39(2):983–1032, 2018. doi:10.1137/16M1098280

  3. [9]

    Gramians, energy functionals, and bal- anced truncation for linear dynamical systems with quadrat ic outputs

    Peter Benner, Pawan Goyal, and Igor Pontes Duff. Gramians, energy functionals, and bal- anced truncation for linear dynamical systems with quadrat ic outputs. IEEE Transactions on Automatic Control, 67(2):886–893, 2021. doi:10.1109/TAC.2021.3086319

  4. [10]

    Sorensen

    Peter Benner, Volker Mehrmann, and Danny C. Sorensen. Dimension Reduction of Large- Scale Systems, volume 45 of Lectures Notes in Computional Science and Engineering . Springer, Berlin, Heidelberg, 2005. doi:10.1007/3-540-27909-1

  5. [11]

    Model Reduc- tion and Approximation: Theory and Algorithms

    Peter Benner, Mario Ohlberger, Albert Cohen, and Karen Willcox. Model Reduc- tion and Approximation: Theory and Algorithms . SIAM, Philadelphia, PA, 2017. doi:10.1137/1.9781611974829

  6. [12]

    Fourier Transforms

    Salomon Bochner and Komaravolu Chandrasekharan. Fourier Transforms. Number 19. Prince- ton University Press, 1949. doi:10.1515/9781400882243

  7. [13]

    Kronecker products and matrix calculus in system theory

    John Brewer. Kronecker products and matrix calculus in system theory. IEEE Transactions on circuits and systems , 25(9):772–781, 1978. doi:10.1109/TCS.1978.1084534

  8. [14]

    Krylov subspace model order reduction of l inear dynamical systems with quadratic output

    Yan-Ping Bu. Krylov subspace model order reduction of l inear dynamical systems with quadratic output. Transactions of the Institute of Measurement and Control , 47(5):827–838,

  9. [15]

    h2-norm optimal model reduction for large scale discrete dynamical MIMO systems

    Angelika Bunse-Gerstner, Dorota Kubalinska, Georg Vo ssen, and Daniel Wilczek. h2-norm optimal model reduction for large scale discrete dynamical MIMO systems. Journal of Compu- tational and Applied Mathematics , 233(5):1202–1216, 2010. Special Issue Dedicated to Willi am B. G...

  10. [16]

    Interpolation-based model order reduction for quadratic-bilinear systems and H2 optimal approximation

    Xingang Cao, Joseph Maubach, Wil Schilders, and Siep We iland. Interpolation-based model order reduction for quadratic-bilinear systems and H2 optimal approximation. In Realization and Model Reduction of Dynamical Systems: A Festschrift in H onor of the 70th Birthday of Thanos...

  11. [17]

    Diaz, Matthias Heinkenschloss, Ion Victo r Gosea, and Athanasios C

    Alejandro N. Diaz, Matthias Heinkenschloss, Ion Victo r Gosea, and Athanasios C. An- toulas. Interpolatory model reduction of quadratic-bilin ear dynamical systems with quadratic-bilinear outputs. Advances in Computational Mathematics , 49(6):1–28, 2023. doi:10.1007/s10444-023-10096-2

  12. [18]

    Realization independent sin- gle time-delay dynamical model interpolation and H2-optimal approximation

    Igor Pontes Duff, Charles Poussot-Vassal, and C´ edric Se ren. Realization independent sin- gle time-delay dynamical model interpolation and H2-optimal approximation. In 2015 54th IEEE Conference on Decision and Control (CDC) , pages 4662–4667. IEEE, 2015. doi:10.1109/CDC.2015.7402946

  13. [19]

    Multipoint Volterra series interpolation and H2 optimal model reduction of bilinear systems

    Garret Flagg and Serkan Gugercin. Multipoint Volterra series interpolation and H2 optimal model reduction of bilinear systems. SIAM Journal on Matrix Analysis and Applications , 36(2):549–579, 2015. doi:10.1137/130947830. Preprint. May 7, 2025 S. Reiter, I. V. Gosea, I. Pontes...

  14. [20]

    Interpolation Methods for the Model Reduction of Bilinear Sy stems

    Garret Michael Flagg. Interpolation Methods for the Model Reduction of Bilinear Sy stems. Dissertation, Virginia Tech, 2012. doi:10919/27521

  15. [21]

    Complex Analysis

    Theodore Gamelin. Complex Analysis . Springer Science & Business Media, New York, NY,

  16. [22]

    Antoulas

    Ion Victor Gosea and Athanasios C. Antoulas. A two-side d iterative framework for model reduction of linear systems with quadratic output. In 2019 IEEE 58th Conference on Decision and Control (CDC) , pages 7812–7817. IEEE, 2019. doi:10.1109/CDC40024.2019.9030025

  17. [23]

    Data-driven mode ling of linear dynamical systems with quadratic output in the AAA framework

    Ion Victor Gosea and Serkan Gugercin. Data-driven mode ling of linear dynamical systems with quadratic output in the AAA framework. Journal of Scientific Computing , 91(1):16,

  18. [24]

    Antoulas, and Christop her Beattie

    Serkan Gugercin, Athanasios C. Antoulas, and Christop her Beattie. H2 model reduction for large-scale linear dynamical systems. SIAM Journal on Matrix Analysis and Applications , 30(2):609–638, 2008. doi:10.1137/06066612

  19. [25]

    Energy matching in reduced passive and port-Hamiltonian systems

    Tobias Holicki, Jonas Nicodemus, Paul Schwerdtner, an d Benjamin Unger. Energy matching in reduced passive and port-Hamiltonian systems. e-prints 2309.05778, arXiv, 2023. URL: https://arxiv.org/abs/2309.05778

  20. [26]

    The commutation matri x: some properties and applica- tions

    Jan R Magnus and Heinz Neudecker. The commutation matri x: some properties and applica- tions. The Annals of Statistics , 7(2):381–394, 1979. doi:10.1214/aos/1176344621

  21. [27]

    Control of port-Ha miltonian differential-algebraic sys- tems and applications

    Volker Mehrmann and Benjamin Unger. Control of port-Ha miltonian differential-algebraic sys- tems and applications. Acta Numerica, 32:395–515, 2023. doi:10.1017/S0962492922000083

  22. [28]

    Approximation of linear c onstant systems

    Lewis Meier and D Luenberger. Approximation of linear c onstant systems. IEEE Transactions on Automatic Control , 12(5):585–588, 1967. doi:10.1109/TAC.1967.1098680

  23. [29]

    Balanced truncation of descriptor systems with a quadratic output

    Jennifer Przybilla, Igor Pontes Duff, Pawan Goyal, and Pe ter Benner. Balanced truncation of descriptor systems with a quadratic output. e-prints 240 2.14716, arXiv, 2024. URL: https://arxiv.org/abs/2402.14716

  24. [30]

    Energy-based model order reduction for l inear stochastic Galerkin systems of second order

    Roland Pulch. Energy-based model order reduction for l inear stochastic Galerkin systems of second order. PAMM, 23(3):e202300038, 2023. doi:10.1002/pamm.20230003833

  25. [31]

    Balanced truncation for model order reduction of linear dy- namical systems with quadratic outputs

    Roland Pulch and Akil Narayan. Balanced truncation for model order reduction of linear dy- namical systems with quadratic outputs. SIAM Journal on Scientific Computing , 41(4):A2270– A2295, 2019. doi:10.1137/17M1148797

  26. [32]

    H2-optimal model reduction of linear quadratic-output systems by multivariate ration al interpolation

    Sean Reiter. Code, data, and results for numerical expe riments in “H2-optimal model reduction of linear quadratic-output systems by multivariate ration al interpolation” (version 1.0), May

  27. [34]

    Sean Reiter and Steffen W. R. Werner. Interpolatory model order reduction of large-scale dynamical systems with root mean squared error measures. e- prints 2403.08894, arXiv, 2024. URL: https://arxiv.org/abs/2403.08894

  28. [35]

    Nonlinear System Theory

    Wilson John Rugh. Nonlinear System Theory . Johns Hopkins University Press, Baltimore, MD, 1981. ISBN: O-8018-2549-0, Web version prepared in 2002

  29. [36]

    Balanced truncation of linear systems with quadratic outpu ts in limited time and frequency intervals

    Qiu-Yan Song, Umair Zulfiqar, Zhi-Hua Xiao, Mohammad Mo nir Uddin, and Victor Sreeram. Balanced truncation of linear systems with quadratic outpu ts in limited time and frequency intervals. e-prints 2402.11445, arXiv, 2024. URL: https://arxiv.org/abs/2402.11445

  30. [37]

    doi:10.5281/zenodo.15319961

  31. [38]

    Model reduction for dynamical systems with quadratic output

    Roel Van Beeumen, Katrien Van Nimmen, Geert Lombaert, a nd Karl Meerbergen. Model reduction for dynamical systems with quadratic output. International Journal for Numerical Methods in Engineering , 91(3):229–248, 2012. doi:10.1002/nme.4255

  32. [39]

    Port-Hamiltonian systems: an int roductory survey

    Arjan van der Schaft. Port-Hamiltonian systems: an int roductory survey. In International congress of mathematicians , pages 1339–1365. European Mathematical Society Publishi ng House (EMS Ph), 2006. doi:10.4171/022-3/65

  33. [40]

    H2-optimal model re- duction of MIMO systems

    Paul van Dooren, Kyle A Gallivan, and P-A Absil. H2-optimal model re- duction of MIMO systems. Applied Mathematics Letters , 21(12):1267–1273, 2008. doi:10.1016/j.aml.2007.09.015

  34. [41]

    Gallivan, and P.-A

    Paul van Dooren, Kyle A. Gallivan, and P.-A. Absil. H2-optimal model reduction with higher- order poles. SIAM Journal on Matrix Analysis and Applications , 31(5):2738–2753, 2010. doi:10.1137/080731591

  35. [42]

    Model reduction b y balanced truncation of linear systems with a quadratic output

    Roel Van Beeumen and Karl Meerbergen. Model reduction b y balanced truncation of linear systems with a quadratic output. In AIP Conference Proceedings, volume 1281, pages 2033–

  36. [43]

    David A. Wilson. Optimum solution of model-reduction p roblem. In Proceedings of the Institution of Electrical Engineers , volume 117, pages 1161–1165. IET, 1970. doi:10.1049/piee.1970.0227

  37. [45]

    Using Krylov-Pad´ e model o rder reduction for accelerating design optimization of structures and vibrations in the fre quency domain

    Yao Yue and Karl Meerbergen. Using Krylov-Pad´ e model o rder reduction for accelerating design optimization of structures and vibrations in the fre quency domain. International Journal for Numerical Methods in Engineering , 90(10):1207–1232, 2012. doi:10.1002/nme.3357

  38. [46]

    Accelerating optimizatio n of parametric linear sys- tems by model order reduction

    Yao Yue and Karl Meerbergen. Accelerating optimizatio n of parametric linear sys- tems by model order reduction. SIAM Journal on Optimization , 23(2):1344–1370, 2013. doi:10.1137/120869171. Preprint. May 7, 2025 S. Reiter, I. V. Gosea, I. Pontes Duff, S. Gugercin: H2-optimal MO...

  39. [47]

    H2-optimal model reduction of linear quadratic output system s in finite frequency range

    Umair Zulfiqar, Zhi-Hua Xiao, Qiu-Yan Song, Mohammad Mo nir Uddin, and Victor Sreeram. H2-optimal model reduction of linear quadratic output system s in finite frequency range. e- prints 2408.07939, arXiv, 2024. URL: https://arxiv.org/abs/2408.07939

  40. [48]

    Steffen W. R. Werner. Structure-Preserving Model Reduction for Mechanical Syst ems. Disser- tation, Otto-von-Guericke-Universit¨ at, Magdeburg, Germany, 2021. doi:10.25673/38617

  41. [54]

    Time-limitedH2-optimal model order reduction of linear systems with quadr atic outputs

    Umair Zulfiqar, Zhi-Hua Xiao, Qiu-Yan Song, Mohammad Mo nir Uddin, and Victor Sreeram. Time-limitedH2-optimal model order reduction of linear systems with quadr atic outputs. e- prints 2408.05965, arXiv, 2024. URL: https://arxiv.org/abs/2408.05965. Preprint. May 7, 2025

  42. [2003]

    doi:10.1007/978-0-387-21607-2

  43. [2022]

    doi:10.1007/s10915-022-01771-5

  44. [2023]

    doi:10.1016/j.jsv.2022.117363

  45. [2025]

    doi:10.1177/01423312241257298

  46. [2036]

    doi:10.1063/1.3498345

    American Institute of Physics, 2010. doi:10.1063/1.3498345

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.