Abstract:The Quadrature-based Balanced Truncation (QuadBT) framework of <a class="link-https" data-arxiv-id="2104.01006" href="https://arxiv.org/abs/2104.01006">arXiv:2104.01006</a> is a "non-intrusive" reformulation of balanced truncation; a classical projection-based model-order reduction technique for linear systems. QuadBT is non-intrusive in the sense that it builds approximate balanced reduced-order models entirely from system response data (e.g., transfer function measurements) without the need to reference an explicit state-space realization of the underlying full-order model. In this work, we generalize and extend QuadBT to other types of balanced truncation model reduction. Namely, we develop non-intrusive implementations for balanced stochastic truncation, positive-real balanced truncation, and bounded-real balanced truncation. We show that the data-driven construction of these balanced reduced-order models requires sampling certain spectral factors associated with the system of interest. Numerical examples are included in each case to validate our approach.
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is: how to extend the model order - reduction method based on balanced truncation (BT) to other types of balanced order - reduction models, and achieve non - invasive (data - driven) implementation in these extensions. Specifically, the authors hope to develop a general framework that enables the construction of various types of balanced order - reduction models from the input - output data of the system without referring to the explicit state - space implementation of the system.
### Main contributions of the paper
1. **Extension of the existing data - driven balanced truncation method**:
- Starting from Quadrature - based Balanced Truncation (QuadBT), the authors extended it to other types of balanced truncation methods, such as Balanced Stochastic Truncation (BST), Positive - Real Balanced Truncation (PRBT), and Bounded - Real Balanced Truncation (BRBT).
2. **Proposed non - invasive implementation method**:
- These extended balanced truncation methods can be implemented through system response data (for example, transfer function measurements), without the need to access the internal dynamics of the system (i.e., the state - space form). This method is called "non - invasive" or "data - driven".
3. **Solved the problem of sampling required for different balanced order - reduction models**:
- The authors detailed what types of data need to be sampled for different balanced order - reduction models. Specifically, they showed how to implement these data - driven order - reduction models by sampling certain spectral factors related to the system.
### Key technical details
- **Calculation of Gramians**: The core of the balanced truncation method is Gramians (observability Gramian and controllability Gramian), and these matrices are obtained by solving Lyapunov equations or Riccati equations.
- **Square - root algorithm**: For numerical stability, the square - root algorithm is usually used to avoid direct calculation of Gramians. The square - root algorithm achieves balancing and order - reduction through singular value decomposition (SVD).
- **Numerical integration and quadratic form**: In the data - driven method, the approximation of Gramians is achieved through numerical integration, in which a quadratic - form (quadrature - based) method is used. The specific quadratic - form factors are extracted from the frequency - domain response of the transfer function.
### Formula summary
The key formulas involved in the paper include:
- **Lyapunov equations**:
\[
A^\top Q+QA + C^\top C = 0
\]
\[
AP + PA^\top+BB^\top = 0
\]
- **Riccati equations** (used for variants such as BST, PRBT, BRBT):
\[
A^\top QW+QWA+(C - B^\top WQW)^\top(DD^\top)^{-1}(C - B^\top WQW) = 0
\]
\[
A^\top QM+QMA+(C - B^\top QM)^\top(D + D^\top)^{-1}(C - B^\top QM) = 0
\]
\[
A^\top QJ+QJA+C^\top C+(B^\top QJ + D^\top C)^\top(I_m - D^\top D)^{-1}(B^\top QJ + D^\top C) = 0
\]
- **Transfer function**:
\[
G(s)=C(sI_n - A)^{-1}B + D
\]
- **Quadratic - form factor**:
\[
\tilde{U}_X=\rh