## Dimension reduction for large-scale stochastic systems with non-zero initial states and controlled diffusion

Martin Redmann ∗

August 2, 2024

## Abstract

In this paper, we establish new strategies to reduce the dimension of large-scale controlled stochastic differential equations with non-zero initial states. The first approach transforms the original setting into a stochastic system with zero initial states. This transformation naturally leads to equations with controlled diffusion. A detailed analysis of dominant subspaces and bounds for the reduction error is provided in this controlled diffusion framework. Subsequently, we introduce a reduced system for the original framework and prove an a-priori error bound for the first ansatz. This bound involves so-called Hankel singular values that are linked to a new pair of Gramians. A second strategy is presented that is based on the idea of reducing control and initial state dynamics separately. Here, different Gramians are used in order to derive a reduced model and their relation to dominant subspaces are pointed out. We also show an a posteriori error bound for the second approach involving two types of Hankel singular values.

Keywords: model order reduction · stochastic systems · controlled diffusion · nonzero initial states · Gramians · matrix (in)equalities · error bounds MSC classification: 60H10 · 60H35 · 60J65 · 65C30 · 68Q25 · 93E03

## 1 Introduction

The analysis and control of large-scale stochastic systems is a significant area of research due to their wide applicability in engineering, finance, biology, and other fields. These systems are often described by controlled stochastic differential equations (SDEs) which can be computationally demanding to solve, particularly when the system size is large. Dimension reduction techniques offer a viable solution to mitigate this complexity by approximating high-dimensional systems with lower-dimensional counterparts while preserving essential system characteristics.

Traditional system-theoretical dimension reduction methods, such as balanced truncation [3, 6, 7, 15] and the iterative rational Krylov algorithm [10, 18], are well-established

∗ Martin Luther University Halle-Wittenberg, Institute of Mathematics, Theodor-Lieser-Str. 5, 06120 Halle (Saale), Germany, Email: martin.redmann@mathematik.uni-halle.de .

/negationslash for both deterministic and stochastic systems with zero initial states. However, many practical systems operate under non-zero initial conditions, making these classical approaches less effective. For linear systems of ordinary differential equation, there have been several approaches to overcome this issue [1, 2, 11, 21]. Having an initial state x 0 = 0 in a stochastic setting causes an enormous increase in complexity as millions of system evaluation might be required for each single x 0 (e.g. in a Monte-Carlo simulation). For that reason, first techniques for SDEs with non-zero initial data have been investigated [4, 19]. However, these approaches come with error bounds either depending on the terminal time T (exploding as T →∞ ) or on singular values of the so-called error system being practically less useful. This paper addresses a gap by developing new and very general strategies with efficient error bound for dimension reduction tailored to large-scale stochastic systems with non-zero initial states.

Our first approach involves the transformation of the original stochastic system with non-zero initial states into an equivalent system with zero initial states. This transformation aims to simplify the analysis and reduction process. However, it leads to SDEs with controlled diffusion terms, a setting, where no model reduction method is available so far. Within this framework, we conduct a detailed analysis of the dominant subspaces and provide bounds for the reduction error based on a new pair of Gramians. These bounds are crucial for ensuring that the reduced system accurately captures the behavior of the original system. One of the main results is an a-priori error bound that involves the truncated Hankel singular values, linked to the novel pair of Gramians, associated to the transformed system with controlled diffusion. Besides the result on the error bound, the advantage of our first approach is that it is very flexible and general. It, e.g., covers the two works [1, 21] at the same time when setting the diffusion equal to zero (deterministic case). The costs of the first approach are the computation of the newly proposed Gramians being solutions of linear matrix inequalities. Solving such inequalities in high dimensions is very challenging which can be seen as a drawback as long as no efficient solver is available.

Secondly, we propose to separate the control and initial state dynamics within the dimension reduction procedure. By treating these components individually, we employ two structurally different reachability Gramians to derive a reduced model. This method leverages the relationship between these Gramians and the dominant subspaces, allowing for a more nuanced reduction approach. Additionally, we establish an a posteriori error bound for this second method, providing further insight into the accuracy of the reduced model. This error bound involves the two different Hankel singular values of the control and initial state dynamics, demonstrating their critical role in both reduction procedures. The benefit of this approach is the flexibility of approximating the two underlying dynamics, with possibly different reduction potential, by reduced order systems not having the same dimensions. On the other hand, it might be computationally more involved as two pairs of Gramians have to be determined. Let us also mention that the second ansatz can be viewed as a generalization of the work on model order reduction for deterministic systems in [2].

Overall, our work contributes to the field of efficiently solving large-scale stochastic

systems with non-zero initial states. In detail, Section 2 introduces the stochastic setting with controlled drift and diffusion as well as a general initial condition. It sketches the two strategies mentioned above and delivers the required transformation of the original equations into a system with zero initial states. Section 3 provides a comprehensive theory on Gramian-based model reduction for such systems received after the transformation. The main contributions are the characterization of dominant subspace based a new pair of Gramians and an a-priori error bound. The same theory is established for uncontrolled systems with non-zero initial data in Section 4. The results of Sections 3 and 4 are the basis for the theory of dimension reduction for the original system. It is provided in Section 5 showing effective error bounds for our two approaches based on various truncated Hankel singular values.

## 2 Setting, notation and goal

Let ( Ω , F , ( F t ) t ∈ [0 ,T ] , P ) 1 be a filtered probability space on which every stochastic process appearing in this paper is defined. Given an R q -valued Wiener process w = [ w 1 . . . w q ] /latticetop , we assume that it is ( F t ) t ∈ [0 ,T ] -adapted and its increments w t ( + h ) -w t ( ) are independent of F t for t, h ≥ 0 and t + h ≤ T . We denote its covariance matrix by K = ( k ij ) i,j =1 ,...,d , so that E [ w t w t ( ) ( ) /latticetop ] = Kt holds. We study the following large-scale stochastic system with a non-zero initial state as well as a controlled drift and diffusion:

$$d x ( t ) & = [ A x ( t ) + B u ( t ) ] d t + \sum _ { i = 1 } ^ { q } [ N _ { i } x ( t ) + M _ { i } u ( t ) ] d w _ { i } ( t ), \quad x ( 0 ) = x _ { 0 } = X _ { 0 } v, \quad ( 1 a ) \\ y ( t ) & = C x ( t ) + D u ( t ), \quad t \in [ 0, T ],$$

$$y ( t ) = C x ( t ) + D u ( t ), \quad t \in [ 0, T ],$$

where A,N i ∈ R n × n , B,M i ∈ R n × m , C ∈ R p × n , D ∈ R p × m , X 0 ∈ R n × d and v ∈ R d is a generic vector. The image of the matrix X 0 represents the space of initial conditions that we consider. Moreover, the ( F t ) t ∈ [0 ,T ] -adapted stochastic control u takes values in R m and satisfies ‖ u ‖ 2 L 2 T := E ∫ T 0 ‖ u t ( ) ‖ 2 2 dt &lt; ∞ , where ‖·‖ 2 denotes the Euclidean norm. We assume that (1a) is mean square asymptotically stable for the rest of the paper. That is E ‖ x t ( ) ‖ 2 2 → 0 for u ≡ 0 and any x 0 ∈ R n , see, e.g., [12, 14]. We introduce an auxiliary variable ˜ associated to (1a) defined by the following system x

$$d \tilde { x } ( t ) & = \tilde { A } \tilde { x } ( t ) d t + \sum _ { i = 1 } ^ { q } \tilde { N } _ { i } \tilde { x } ( t ) d w _ { i } ( t ), \quad \tilde { x } ( 0 ) = x _ { 0 } = X _ { 0 } v, \\ \tilde { y } ( t ) & = C \tilde { x } ( t ).$$

$$\tilde { y } ( t ) = C \tilde { x } ( t ).$$

Here, ˜ A,N ˜ i ∈ R n × n are matrices that are arbitrary and that are chosen according to a desired ansatz. In fact, (2) is an uncontrolled version of (1), but mean square asymptotic stability is not necessarily assumed. However, both systems share the same initial state

1 ( F t ) t ∈ [0 ,T ] is right continuous and complete.

and output matrix C . We can now introduce a variable z t ( ) := x t ( ) -˜( ) having a zero x t initial condition. The corresponding system is

$$d z ( t ) = \left [ A z ( t ) + \left [ B \ A - \tilde { A } \right ] \left [ \begin{smallmatrix} u ( t ) \\ \vec { x } ( t ) \end{smallmatrix} \right ] \right ] d t + \sum _ { i = 1 } ^ { q } \left [ N _ { i } z ( t ) + \left [ M _ { i } \ \ N _ { i } - \tilde { N } _ { i } \right ] \left [ \begin{smallmatrix} u ( t ) \\ \vec { x } ( t ) \end{smallmatrix} \right ] \right ] d w _ { i } ( t ), \\ r \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$y ( t ) = C z ( t ) + \left [ D \ \ C \right ] \begin{bmatrix} u ( t ) \\ \tilde { x } ( t ) \end{bmatrix} = C z ( t ) + D u ( t ) + \tilde { y } ( t ), & \quad ( 3 b ) \\ \text{with } z ( 0 ) = 0 \text{ and } t \in [ 0, T ]. \ In \text{ fact, we aim to follow two approaches using the auxiliary}$$

with z (0) = 0 and t ∈ [0 , T ]. In fact, we aim to follow two approaches using the auxiliary variable ˜ in different ways. x These are discussed in the next two subsections.

## 2.1 First approach

Given X 0 , the matrices ˜ , A N ,... , N ˜ 1 ˜ q in (2) are chosen, so that the auxiliary variable ˜ x takes values in a (low-dimensional) subspace im[ V 0 ] for each x 0 ∈ im[ X 0 ], where V 0 ∈ R n × r 0 is a full rank matrix with r 0 ≤ n . In addition, we assume V 0 to be given explicitly without computational cost. We can now represent ˜( ) = x t V u 0 0 ( ) t v for some suitable R r 0 × d -valued process u 0 which we assume to be easily accessible. Hence, we interpret u v 0 to be another externally given control process. Inserting this into (3) yields

$$d z ( t ) & = \left [ A z ( t ) + \tilde { B } \tilde { u } ( t ) \right ] d t + \sum _ { i = 1 } ^ { q } \left [ N _ { i } z ( t ) + \tilde { M } _ { i } \tilde { u } ( t ) \right ] d w _ { i } ( t ), \quad z ( 0 ) = 0, \quad ( 4 \tt a ) \\ u ( t ) & = C z ( t ) + \tilde { D } \tilde { u } ( t ), \quad t \in [ 0, T ],$$

$$y ( t ) = C z ( t ) + \tilde { D } \tilde { u } ( t ), \quad t \in [ 0, T ],$$

where we define

$$( \tilde { B }, \tilde { M } _ { i }, \tilde { D }, \tilde { u } ( t ) ) \coloneqq \left ( \left [ B \ \ ( A - \tilde { A } ) V _ { 0 } \right ], \left [ M _ { i } \ \ ( N _ { i } - \tilde { N } _ { i } ) V _ { 0 } \right ], \left [ D \ \ C V _ { 0 } \right ], \left [ \ u _ { 0 } ( t ) v \right ] \right ). \\ \text{Consequently,} ( 3 ) \, \text{has become a system of the form} \ ( 1 ), \, \text{with a zero initial state, updated}$$

Consequently, (3) has become a system of the form (1), with a zero initial state, updated control matrices ˜ B,M ,D ˜ i ˜ and new control ˜. u The benefit of this first approach is that we save the step of reducing the dynamics corresponding to the initial state. This is due to the flexibility of choosing ˜ A and ˜ N i which is used to find a low-dimensional underlying structure without computational effort. However, the input space is enlarged in comparison to the original system, since we include the information regarding the initial state there. Notice that we have ˜( ) = Φ( ) x t ˜ t x 0 , where ˜ is the fundamental matrix. Φ In case all matrices ˜ and A ˜ N i commute it is Φ( ) = e ˜ t ( ˜ A -0 5 . ∑ q i,j =1 ˜ N N k i ˜ j ij ) + t ∑ q i =1 ˜ N w i i ( ) t . If we now set ˜ A = -αI and ˜ N i = γ I i for α &gt; 0 and γ i ∈ R , we obtain ˜( ) x t = e ( - -α 0 5 . ∑ q i,j =1 γ γ k i j ij ) + t ∑ q i =1 γ w i i ( ) t x 0 . Since x 0 = X v 0 , we obtain that V 0 = X 0 and u 0 ( ) = e t ( - -α 0 5 . ∑ q i,j =1 γ γ k i j ij ) + t ∑ q i =1 γ w i i ( ) t I d × d . Setting γ i = 0 for all i = 1 , . . . , q leads to an approach for ordinary differential equations that was studied in [21]. The most trivial choice is ˜ = A ˜ N i = 0 leading to ˜( ) = x t x 0 = X v 0 . Hence, V 0 = X 0 and u 0 ≡ I d × d . Such a scenario was studied in the deterministic case in [1]. However, this results in a

non-stable system (2) having a few disadvantages, e.g., an error estimate increasing in time. Certainly, we can also allow for using V 0 = I and u 0 ( ) = Φ( ) t ˜ t X 0 if ˜ Φ is explicitly available making u v 0 ≡ ˜ the additional control variable. x However, this leads to a large number of inputs, so that the system might be hard to reduce.

Remark 2.1. The structure of (1) is very natural to consider when following the ansatz of this subsection. Even though we might start with an uncontrolled system (1) , we have a controlled system in (4) which means that we potentially encounter a change in the structure. The main challenge is to develop a model reduction scheme for (4) or equivalently (1) with x (0) = 0 , since these systems have an input term in the diffusion. So far, there is no dimension reduction technique available for such a setting.

## 2.2 Second approach

We set ˜ = A A and ˜ N i = N i meaning that ˜ as well as its solution space are not explicitly x given in general. In this case, (3) is a system like in (1) having a zero initial state, but the output equation (3b) contains another variable ˜ that, in contrast to (4), is not known y a-priori:

$$d z ( t ) & = [ A z ( t ) + B u ( t ) ] d t + \sum _ { i = 1 } ^ { q } [ N _ { i } z ( t ) + M _ { i } u ( t ) ] d w _ { i } ( t ), \quad z ( 0 ) = 0, \quad & ( 5 a ) \\ y ( t ) - \tilde { y } ( t ) & = C z ( t ) + D u ( t ), \quad t \in [ 0, T ].$$

$$y ( t ) - \tilde { y } ( t ) = C z ( t ) + D u ( t ), \quad t \in [ 0, T ].$$

Now, computing y requires further effort and solving for ˜ is expensive. y Therefore, we approximate ˜ by applying model order reduction to (2). y

Remark 2.2. Following the ansatz of this subsection, we have to apply model reduction twice. The first step is to reduce the control dynamics. This is system (5) or equivalently system (1) with x (0) = 0 . Subsequently, dimension reduction is applied to (2) in order approximate ˜ y entering (5b) . This strategy can be viewed as a generalization of the work in [2] to stochastic differential equations.

## 3 Model reduction for (1) with x 0 = 0

This section is the foundation for the approaches described in Subsections 2.1 and 2.2. Throughout this section, we assume that x 0 = 0 in (1).

## 3.1 Gramians and dominant subspaces

We begin with identifying state directions that are of low relevance in (1a) while assuming that x 0 = 0 below. For that reason, we define a Gramian based on the following

operators:

$$L ( X ) \coloneqq A ^ { \top } X + X A + \sum _ { \dots } ^ { q } N _ { i } ^ { \top } X N _ { j } k _ { i j },$$

$$L ( X ) & \coloneqq A ^ { \top } X + X A + \sum _ { i, j = 1 } ^ { q } N _ { i } ^ { \top } X N _ { j } k _ { i j }, \\ S ( X ) & \coloneqq X B + \sum _ { i, j = 1 } ^ { q } N _ { i } ^ { \top } X M _ { j } k _ { i j }, \\ U ( X ) & \coloneqq \gamma I - \sum _ { i, j = 1 } ^ { q } M _ { i } ^ { \top } X M _ { j } k _ { i j }, \quad \gamma > 0, \\ X \in \mathbb { R } ^ { n \times n }. \text{ Although } U \text{ depends on the parameter } \gamma, t$$

for some generic X ∈ R n × n . Although U depends on the parameter γ , this is omitted in the rest of the paper. Below, we write M 1 ≥ M 2 ( M &gt; M 1 2 ) for two matrices M ,M 1 2 if M 1 -M 2 is symmetric positive semidefinite (definite). We are now able to state the following definition.

Definition 3.1. Given that x 0 = 0 , a reachability Gramian for (1) is a symmetric matrix P &gt; 0 satisfying U P ( -1 ) &gt; 0 and

$$L ( P ^ { - 1 } ) + S ( P ^ { - 1 } ) U ( P ^ { - 1 } ) ^ { - 1 } S ( P ^ { - 1 } ) ^ { \top } \leq 0.$$

Referring to [6], P can be viewed as a generalization of a so-called type II Gramian studied in that paper. Let us first ensure that a Gramian P according to Definition 3.1 exists.

Proposition 3.2. Let the solution of the uncontrolled equation (1a) be mean square asymptotically stable. Then, (6) has a solution P &gt; 0 with U P ( -1 ) &gt; 0 .

Proof. The underlying mean square asymptotic stability for system (1) is equivalent to the existence of a matrix X &gt; 0, so that L X ( ) = -Y &lt; 0, see [8, 17]. Given an /epsilon1 &gt; 0, we conclude that L /epsilon1X ( ) = -/epsilon1Y . We immediately see that U /epsilon1X ( ) = γI -/epsilon1 ∑ q i,j =1 M XM k /latticetop i j ij &gt; 0 for a sufficiently small /epsilon1 exploiting that ∑ q i,j =1 M XM k /latticetop i j ij ≥ 0. Let λ ≥ 0 denote an arbitrary eigenvalue of ∑ q i,j =1 M XM k /latticetop i j ij . Then, 1 γ -/epsilon1λ is an eigenvalue of U /epsilon1X ( ) -1 . If λ max is the largest out of all these eigenvalues, we obtain

$$\epsilon Y - S ( \epsilon X ) U ( \epsilon X ) ^ { - 1 } S ( \epsilon X ) ^ { \top } & \geq \epsilon Y - \epsilon ^ { 2 } S ( X ) S ( X ) ^ { \top } \frac { 1 } { \gamma - \epsilon \lambda _ { \max } } \\ & \geq I ( \epsilon \lambda _ { \min } ^ { Y } - \epsilon ^ { 2 } \lambda _ { \max } ^ { S S ^ { \top } } \frac { 1 } { \gamma - \epsilon \lambda _ { \max } } ) > 0$$

for /epsilon1 &gt; 0 small enough, where λ Y min &gt; 0 is the smallest and λ SS /latticetop max ≥ 0 the largest eigenvalue of Y and S X S X ( ) ( ) /latticetop , respectively. Using this estimate yields L /epsilon1X ( ) &lt; -S /epsilon1X ( ) U /epsilon1X ( ) -1 S /epsilon1X ( ) /latticetop for a sufficiently small /epsilon1 . For that reason, P = ( /epsilon1X ) -1 is a Gramian according to Definition 3.1.

Since P &gt; 0 is symmetric, we can find a basis of eigenvectors ( p k ) providing a representation for the state variable in (1a):

$$x ( t ) = \sum _ { k = 1 } ^ { n } \langle x ( t ), p _ { k } \rangle _ { 2 } \, p _ { k }.$$

The next theorem answers the question which directions p k can be neglected in (1a).

Theorem 3.3. Let P be as in Definition 3.1 and x be the solution of (1a) . Moreover, let ( p , λ k k ) be an associated basis of eigenvectors giving us (7) together with the eigenvalues λ k &gt; 0 . With x 0 = 0 , we have

$$\sup _ { t \in [ 0, T ] } \mathbb { E } \langle x ( t ), p _ { k } \rangle _ { 2 } ^ { 2 } \leq \lambda _ { k } \gamma \left \| u \right \| _ { L _ { T } ^ { 2 } } ^ { 2 },$$

where γ &gt; 0 enters the definition of the operator U .

Proof. We apply Ito's product rule and subsequently take the expectation, see Lemma A.1 for details. Hence,

$$\frac { d } { d t } \mathbb { E } \left [ x ( t ) ^ { \top } P ^ { - 1 } x ( t ) \right ] = & 2 \mathbb { E } \left [ x ( t ) ^ { \top } P ^ { - 1 } [ A x ( t ) + B u ( t ) ] \right ] \\ & + \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ \left ( N _ { i } x ( t ) + M _ { i } u ( t ) \right ) ^ { \top } P ^ { - 1 } ( N _ { j } x ( t ) + M _ { j } u ( t ) ) \right ] k _ { i j }.$$

We integrate this equation over [0 , t ] with t ≤ T yielding

$$\mathbb { E } \left [ x ( t ) ^ { \top } P ^ { - 1 } x ( t ) \right ] & = \int _ { 0 } ^ { t } \mathbb { E } \left [ x ( s ) ^ { \top } L ( P ^ { - 1 } ) x ( s ) \right ] + 2 \mathbb { E } \left [ x ( s ) ^ { \top } P ^ { - 1 } B u ( s ) \right ] \\ & \quad + \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ 2 x ( s ) ^ { \top } N _ { i } ^ { \top } P ^ { - 1 } M _ { j } u ( s ) + u ( s ) ^ { \top } M _ { i } ^ { \top } P ^ { - 1 } M _ { j } u ( s ) \right ] k _ { i j } d s \\ & = \int _ { 0 } ^ { t } \mathbb { E } \left [ x ( s ) ^ { \top } L ( P ^ { - 1 } ) x ( s ) \right ] + 2 \mathbb { E } \left [ x ( s ) ^ { \top } S ( P ^ { - 1 } ) u ( s ) \right ] \\ & \quad + \gamma \mathbb { E } \left [ u ( s ) ^ { \top } u ( s ) \right ] - \mathbb { E } \left [ u ( s ) ^ { \top } U ( P ^ { - 1 } ) u ( s ) \right ] d s.$$

We insert (6) into the above equation and obtain

$$\text{$\mathbb{E}$} \left [ x ( t ) ^ { \top } P ^ { - 1 } x ( t ) \right ] \leq \int _ { 0 } ^ { t } - \mathbb { E } \left [ x ( s ) ^ { \top } S ( P ^ { - 1 } ) U ( P ^ { - 1 } ) ^ { - 1 } S ( P ^ { - 1 } ) ^ { \top } x ( s ) \right ] \\ & + 2 \mathbb { E } \left [ x ( s ) ^ { \top } S ( P ^ { - 1 } ) u ( s ) \right ] + \gamma \mathbb { E } \| u ( s ) \| _ { 2 } ^ { 2 } - \mathbb { E } \left [ u ( s ) ^ { \top } U ( P ^ { - 1 } ) u ( s ) \right ] d s \\ = \mathbb { E } \int _ { 0 } ^ { t } - \| U ( P ^ { - 1 } ) ^ { - \frac { 1 } { 2 } } S ( P ^ { - 1 } ) ^ { \top } x ( s ) \| _ { 2 } ^ { 2 } - \| U ( P ^ { - 1 } ) ^ { \frac { 1 } { 2 } } u ( s ) \| _ { 2 } ^ { 2 } d s \\ & + \mathbb { E } \int _ { 0 } ^ { t } 2 \langle S ( P ^ { - 1 } ) ^ { \top } x ( s ), u ( s ) \rangle _ { 2 } d s + \gamma \| u \| _ { L _ { t } ^ { 2 } } ^ { 2 } \\ = \gamma \| u \| _ { L _ { t } ^ { 2 } } ^ { 2 } - \mathbb { E } \int _ { 0 } ^ { t } \| U ( P ^ { - 1 } ) ^ { - \frac { 1 } { 2 } } S ( P ^ { - 1 } ) ^ { \top } x ( s ) - U ( P ^ { - 1 } ) ^ { \frac { 1 } { 2 } } u ( s ) \| _ { 2 } ^ { 2 } d s \\ \leq \gamma \| u \| _ { L _ { T } ^ { 2 } } ^ { 2 } \\ \text{exploiting that $U(P^{-1})$ is positive definite.  The statement follows by the observation} \\ \text{+} \dots \text{*}$$

exploiting that U P ( -1 ) is positive definite. The statement follows by the observation that

$$\langle x ( t ), p _ { k } \rangle _ { 2 } ^ { 2 } & \leq \lambda _ { k } \, \sum _ { i = 1 } ^ { n } \lambda _ { i } ^ { - 1 } \langle x ( t ), p _ { i } \rangle _ { 2 } ^ { 2 } = \lambda _ { k } \left \| \sum _ { i = 1 } ^ { n } \lambda _ { i } ^ { - \frac { 1 } { 2 } } \langle x ( t ), p _ { i } \rangle _ { 2 } \, p _ { i } \right \| _ { 2 } ^ { 2 } = \lambda _ { k } \left \| P ^ { - \frac { 1 } { 2 } } x ( t ) \right \| _ { 2 } ^ { 2 } \\ & = \lambda _ { k } \, x ( t ) ^ { \top } P ^ { - 1 } x ( t ).$$

This concludes the proof.

Theorem 3.3 tells us that all eigenspaces corresponding to small eigenvalues of P can be removed only causing a small error. This will be exploited in the dimension reduction procedure later. However, we need to be able to solve for P . Therefore, let us briefly comment on the computation of this Gramian below.

Remark 3.4. Solving (6) seems to be very tricky and is a bottleneck of our approach as long as no efficient solver is available for such matrix inequalities. In fact, a straight forward approach is to solve for P -1 instead of computing P directly. This can be done by reformulating (6) as a linear matrix inequality based on Schur complement conditions for the definiteness of a matrix. This yields

$$\begin{bmatrix} L ( P ^ { - 1 } ) & S ( P ^ { - 1 } ) \\ S ( P ^ { - 1 } ) ^ { \top } & - U ( P ^ { - 1 } ) \end{bmatrix} \leq 0, \quad P ^ { - 1 }, U ( P ^ { - 1 } ) > 0, \\.$$

maximizing, e.g., the trace of P -1 . Here, the hope is to find a P -1 , such that P has a large number of small eigenvalues. This is beneficial within the dimension reduction approach due to Theorem 3.3. Now, a linear matrix inequality solver (e.g. YALMIP and MOSEK [13, 16]) can find a solution to the maximization of tr( P ) subject to (8) .

Below, we aim to investigate which directions in the state variable x (with zero initial condition) are of low relevance for the quantity of interest y in (1b). In this context, let us introduce a second Gramian.

Definition 3.5. Given x 0 = 0 , an observability Gramian for (1) is a symmetric matrix Q ≥ 0 satisfying

$$L ( Q ) + C ^ { \top } C \leq 0.$$

Choosing an arbitrary t 0 ∈ [0 , ∞ ), we denote the solution to (1a) by x t ( ; t 0 , x 0 , u ) for the moment with t 0 representing the initial time and t ≥ t 0 . Further, if t 0 = 0, we set x t ( ) := x t ( ; 0 , x 0 , u ). Exploiting the linearity of the solution in x 0 and u , we obtain

$$x ( t ) = x ( t ; 0, x _ { 0 }, u ) = x ( t ; t _ { 0 }, x ( t _ { 0 } ), u ) = x ( t ; t _ { 0 }, x ( t _ { 0 } ), 0 ) + x ( t ; t _ { 0 }, 0, u ). \quad ( 1 0 )$$

Based on (10), we aim to answer which directions in x t ( 0 ) are less relevant for the output y . Exploiting that an observability Gramian according to Definition 3.5 is symmetric, we choose a basis of eigenvectors ( q k ) of Q for R n yielding

$$x ( t _ { 0 } ) & = \sum _ { k = 1 } ^ { n } \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } \, q _ { k }. & ( 1 1 ) \\ \bigcap \text{and} \, \min \, \text{the in the initial state on one more more than}.$$

Inserting (11) into (10) and using the linearity in the initial state once more, we have

$$y ( t ) = \sum _ { k = 1 } ^ { n } \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } C x ( t ; t _ { 0 }, q _ { k }, 0 ) + C x ( t ; t _ { 0 }, 0, u ) & & ( 1 2 ) \\ > t _ { - } \, \text{ In order to answer which direction } a. \, \text{ can be notated in } u \, \text{ } w \, \text{ formula}.$$

for t ≥ t 0 . In order to answer which directions q k can be neglected in y , we formulate the following.

Proposition 3.6. Let ( q k ) be eigenvectors of Q that represent an orthonormal basis of R n with associated eigenvalues ( µ k ) . Let x t ( ; t 0 , q k , 0) be the solution of (1a) with x 0 = q k , u ≡ 0 and t 0 being the initial time. Then, we have

Proof. We obtain that

$$\ddots \stackrel { \lambda } { \mathbb { E } } \int _ { t _ { 0 } } ^ { \infty } \left \| \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } C x ( t ; t _ { 0 }, q _ { k }, 0 ) \right \| _ { 2 } ^ { 2 } d t \leq \mu _ { k } \mathbb { E } \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } ^ { 2 }. \\ \text{We obtain that } \right \| _ { \lambda } \lambda \right \|$$

$$& \mathcal { E } \left \| \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } C x ( t ; t _ { 0 }, q _ { k }, 0 ) \right \| _ { 2 } ^ { 2 } = \mathbb { E } \left [ \mathbb { E } \left \{ \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } ^ { 2 } \right \| C x ( t ; t _ { 0 }, q _ { k }, 0 ) \right \| _ { 2 } ^ { 2 } \left | \mathcal { F } _ { t _ { 0 } } \right \} \right ] \\ & \quad = \mathbb { E } \left [ \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } ^ { 2 } \mathbb { E } \left \{ \left \| C x ( t ; t _ { 0 }, q _ { k }, 0 ) \right \| _ { 2 } ^ { 2 } \left | \mathcal { F } _ { t _ { 0 } } \right \} \right ] = \mathbb { E } \left [ \langle x ( t _ { 0 } ), q _ { k } \rangle _ { 2 } ^ { 2 } \right ] \mathbb { E } \left \| C x ( t ; t _ { 0 }, q _ { k }, 0 ) \right \| _ { 2 } ^ { 2 } \\ & \text{exploiting that } x ( t _ { 0 } ) \text{ is } \mathcal { F } _ { t _ { 0 } } \text{-measurable and } x ( t ; t _ { 0 }, q _ { k }, 0 ) \text{ is independent of } \mathcal { F } _ { t _ { 0 } }. \text{ Using} \\ & \text{the arguments in the proof of Theorem 3.3 based on Ito's product rule (see also Lemma}$$

∥ ∥ ∣ ∥ ∥ exploiting that x t ( 0 ) is F t 0 -measurable and x t ( ; t 0 , q k , 0) is independent of F t 0 . Using the arguments in the proof of Theorem 3.3 based on Ito's product rule (see also Lemma A.1), we obtain

$$\frac { \, ^ { \prime \prime } \, \cdots \, } { d t } \mathbb { E } \left [ x ( t ; t _ { 0 }, q _ { k }, 0 ) ^ { \top } Q x ( t ; t _ { 0 }, q _ { k }, 0 ) \right ] = & 2 \mathbb { E } \left [ x ( t ; t _ { 0 }, q _ { k }, 0 ) ^ { \top } Q A x ( t ; t _ { 0 }, q _ { k }, 0 ) \right ] \\ & + \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ \left ( N _ { i } x ( t ; t _ { 0 }, q _ { k }, 0 ) \right ) ^ { \top } Q N _ { j } x ( t ; t _ { 0 }, q _ { k }, 0 ) \right ] k _ { i j }.$$

We integrate this equation over [ t 0 , t ] giving us

$$\mathbb { E } \left [ x ( t ; t _ { 0 }, q _ { k }, 0 ) ^ { \top } Q x ( t ; t _ { 0 }, q _ { k }, 0 ) \right ] & = q _ { k } ^ { \top } Q q _ { k } + \int _ { t _ { 0 } } ^ { t } \mathbb { E } \left [ x ( s ; t _ { 0 }, q _ { k }, 0 ) ^ { \top } L ( Q ) x ( s ; t _ { 0 }, q _ { k }, 0 ) \right ] d s. \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

We insert (9) into this equation resulting in

$$\int _ { t _ { 0 } } ^ { t } \mathbb { E } \left \| C x ( s ; t _ { 0 }, q _ { k }, 0 ) \right \| _ { 2 } ^ { 2 } d s \leq \mu _ { k } \\ \text{xy taking } t \to \infty.$$

and the result follows by taking t →∞ .

Proposition 3.6 together with (12) tell us that directions q k of x t ( 0 ) in (11) contribute very little to y t ( ), t ≥ t 0 for each t 0 ∈ [0 , ∞ ) in case µ k is small. Such directions will be neglected in the reduced dynamics.

## 3.2 Balancing transformation and reduced order model

Below, we further suppose that x 0 = 0. Using a regular matrix T ∈ R n × n , we introduce a balanced state variable by x b = Tx . This transformation is chosen in a particular way in order to simultaneously diagonalize P and Q introduced in Definitions 3.1 and 3.5. Such a T always exists under the condition that Q &gt; 0. As a consequence of the diagonalization, the associated bases of eigenvectors ( p k ) and ( q k ) will be the canonical basis of R n . Referring to Theorem 3.3 and Proposition 3.6 this means that components in the transformed state variable x b turn out to be of low relevance if the associated diagonal entries of the diagonalized Gramians are small. We now insert x = T -1 x b into (1) leading to a stochastic system with coefficients

$$( A _ { b }, B _ { b }, M _ { i, b }, N _ { i, b }, C _ { b } ) \coloneqq ( T A T ^ { - 1 }, T B, T M _ { i }, T N _ { i } T ^ { - 1 }, C T ^ { - 1 } )$$

instead of the original ones ( A,B,M ,N ,C i i ), i.e.,

$$d x _ { b } ( t ) & = [ A _ { b } x _ { b } ( t ) + B _ { b } u ( t ) ] d t + \sum _ { i = 1 } ^ { q } [ N _ { i, b } x _ { b } ( t ) + M _ { i, b } u ( t ) ] d w _ { i } ( t ), \quad & ( 1 3 a ) \\ y ( t ) & = C _ { h } x _ { b } ( t ) + D u ( t ), \quad x _ { b } ( 0 ) = 0, \quad t \in [ 0, T ]. \quad & ( 1 3 b )$$

$$y ( t ) = C _ { b } x _ { b } ( t ) + D u ( t ), \ \ x _ { b } ( 0 ) = 0, \ \ t \in [ 0, T ].$$

Systems (1) and (13) have the same input-output behaviour and are therefore considered to be equivalent. They share properties such as asymptotic stability. However, the Gramians are different. The particular choice of T as well as the form of the Gramians of (13) are stated in the following proposition.

Proposition 3.7. Suppose that T is an invertible matrix. Given the Gramians P and Q of the original system (1) with x 0 = 0 according to Definitions 3.1 and 3.5. Then, the Gramians of (13) are P b = TPT /latticetop and Q b = T -/latticetop QT -1 . Suppose that Q &gt; 0 , we

can find a simultaneous diagonalization, i.e., P b = Q b = Σ = diag( σ , . . . , σ 1 n ) using the balancing transformation

$$T = \Sigma ^ { \frac { 1 } { 2 } } U ^ { \top } L _ { P } ^ { - 1 },$$

where P = L L P /latticetop P and L QL /latticetop P P = U Σ 2 U /latticetop is a spectral factorization with an orthogonal U .

Proof. We multiply (6) with T -/latticetop from the left and with T -1 from the right. This results in

$$T ^ { - \top } L ( P ^ { - 1 } ) T ^ { - 1 } + T ^ { - \top } S ( P ^ { - 1 } ) U ( P ^ { - 1 } ) ^ { - 1 } S ( P ^ { - 1 } ) ^ { \top } T ^ { - 1 } \leq 0.$$

We look at each term individually and observe

$$T ^ { - \top } L ( P ^ { - 1 } ) T ^ { - 1 } & = T ^ { - \top } A ^ { \top } T ^ { \top } T ^ { - \top } P ^ { - 1 } T ^ { - 1 } + T ^ { - \top } P ^ { - 1 } T ^ { - 1 } T A T ^ { - 1 } \\ & \quad + \sum _ { i, j = 1 } ^ { q } T ^ { - \top } N _ { i } ^ { \top } T ^ { \top } T ^ { - \top } P ^ { - 1 } T ^ { - 1 } T N _ { j } T ^ { - 1 } k _ { i j } \\ & = A ^ { \top } P ^ { - 1 } + P ^ { - 1 } A _ { i j } + \sum _ { j = 1 } ^ { q } N ^ { \top } _ { j } P ^ { - 1 } N _ { j } \colon k _ { j j }.$$

$$= A _ { b } ^ { \top } P _ { b } ^ { - 1 } + P _ { b } ^ { - 1 } A _ { b } + \sum _ { i, j = 1 } ^ { q } N _ { i, b } ^ { \top } P _ { b } ^ { - 1 } N _ { j, b } k _ { i j }$$

as well as

$$& \text{as well as} \\ & T ^ { - \top } S ( P ^ { - 1 } ) = T ^ { - \top } P ^ { - 1 } T ^ { - 1 } T B + \sum _ { i, j = 1 } ^ { q } T ^ { - \top } N _ { i } ^ { \top } T ^ { \top } T ^ { - \top } P ^ { - 1 } T ^ { - 1 } T M _ { j } k _ { i j } \\ & \quad = P _ { b } ^ { - 1 } B _ { b } + \sum _ { i, j = 1 } ^ { q } N _ { i, b } ^ { \top } P _ { b } ^ { - 1 } M _ { j, b } k _ { i j } \\ & \text{and} \, U ( P ^ { - 1 } ) = \gamma I - \sum _ { i \colon } ^ { q } \,, \, M _ { i \colon } ^ { \top } T ^ { \top } T ^ { - \top } P ^ { - 1 } T ^ { - 1 } T M _ { i } k _ { i i } = \gamma I - \sum _ { i \colon } ^ { q } \,, \, M _ { i \colon } ^ { \top } P _ { i \colon } ^ { - 1 } M _ { i \colon } k _ { i i }.$$

and U P ( -1 ) = γI -∑ q i,j =1 M T /latticetop i /latticetop T -/latticetop P -1 T -1 TM k j ij = γI -∑ q i,j =1 M P /latticetop i,b -1 b M k j,b ij . This shows that P b is the reachability Gramian of (13). Let us now multiply (9) with T -/latticetop from the left and with T -1 from the right. Then,

$$T ^ { - \top } L ( Q ) T ^ { - 1 } + C _ { b } ^ { \top } C _ { b } \leq 0$$

shows that Q b is the observability Gramian of (13) taking into account that T -/latticetop L Q T ( ) -1 = A Q /latticetop b b + Q A b b + ∑ q i,j =1 N /latticetop i,b Q N b j,b k ij . Since P and Q are invertible, the same is true for L P and Σ. Therefore, T and its inverse T -1 = L U P Σ -1 2 are well-defined. Inserting this into the definitions of P b and Q b provides that both Gramians are equal to Σ. This completes the proof.

Since L QL /latticetop P P has the same eigenvalues as PQ it becomes clear that the Hankel singular values (HSVs) σ i are given by σ i = √ λ i ( PQ ), where λ i ( ) · denotes the i th eigenvalue of a matrix with real spectrum. We partition the balanced state x b = [ x 1 x 2 ] ,

x 1 ∈ R r , and the diagonal matrix of HSVs Σ = diag(Σ 1 , Σ ), where Σ 2 1 = diag( σ , . . . , σ 1 r ) contains the large and Σ 2 = diag( σ r +1 , . . . , σ n ) the small HSVs, r &lt; n . Accordingly, we partition the balanced realization as follows:

$$A _ { b } & = \begin{bmatrix} A _ { 1 1 } & A _ { 1 2 } \\ A _ { 2 1 } & A _ { 2 2 } \end{bmatrix}, \quad B _ { b } = \begin{bmatrix} B _ { 1 } \\ B _ { 2 } \end{bmatrix}, \quad M _ { i, b } \ = \begin{bmatrix} M _ { i, 1 } \\ M _ { i, 2 } \end{bmatrix}, \quad N _ { i, b } = \begin{bmatrix} N _ { i, 1 1 } & N _ { i, 1 2 } \\ N _ { i, 2 1 } & N _ { i, 2 2 } \end{bmatrix}, \quad C _ { b } = \begin{bmatrix} C _ { 1 } & C _ { 2 } \end{bmatrix}.$$

We set x 2 ≡ 0 in the equations of x 1 and y , since x 2 is of low relevance. Truncating the equation for x 2 leads to the following reduced order model

$$d x _ { r } ( t ) & = [ A _ { 1 1 } x _ { r } ( t ) + B _ { 1 } u ( t ) ] d t + \sum _ { i = 1 } ^ { q } [ N _ { i, 1 1 } x _ { r } ( t ) + M _ { i, 1 } u ( t ) ] d w _ { i } ( t ), \quad & ( 1 6 a ) \\ y _ { r } ( t ) & = C _ { 1 } x _ { r } ( t ) + D u ( t ), \quad x _ { r } ( 0 ) = 0, \quad t \in [ 0, T ].$$

$$y _ { r } ( t ) = C _ { 1 } x _ { r } ( t ) + D u ( t ), \ \ x _ { r } ( 0 ) = 0, \ \ t \in [ 0, T ].$$

For simplicity of the notation, let us for the remainder of this subsection assume that the original system is already balanced, i.e., x = x b as well as ( A,B,M ,N ,C i i ) = ( A , B , M b b i,b , N i,b , C b ). Consequently, we deal with matrix inequalities of the following form

$$L ( \Sigma ^ { - 1 } ) + S ( \Sigma ^ { - 1 } ) U ( \Sigma ^ { - 1 } ) ^ { - 1 } S ( \Sigma ^ { - 1 } ) ^ { \top } \leq 0, \ \text{ and } \ L ( \Sigma ) + C ^ { \top } C \leq 0. \quad \ ( 1 7 )$$

We now derive the same type of matrix inequalities for the reduced system (16).

Proposition 3.8. Given the balanced realization ( A,B,M ,N ,C i i ) and introducing the reduced operators

$$L _ { 1 1 } ( X ) \coloneqq A _ { 1 1 } ^ { \top } X + & X A _ { 1 1 } + \sum _ { i, j = 1 } ^ { q } N _ { i, 1 1 } ^ { \top } X N _ { j, 1 1 } k _ { i j }, \\ \underline { q } \quad -$$

$$U _ { 1 } ( X ) \coloneqq \gamma I - \sum _ { i, j = 1 } ^ { q } M _ { i, 1 1 } ^ { \top } X M _ { j, 1 1 } k _ { i j }, \quad \gamma > 0, \\.$$

$$S _ { 1 } ( X ) \coloneqq X B _ { 1 } + \sum _ { \begin{subarray} { c } i, j = 1 \\ \frac { q } { 2 } \end{subarray} } ^ { c, j - 1 } N _ { i, 1 1 } ^ { \top } X M _ { j, 1 1 } k _ { i j }, \\ \underline { q } \quad -$$

for a generic X ∈ R r × r , it holds that U 1 (Σ -1 1 ) &gt; 0 and

$$L _ { 1 1 } ( \Sigma _ { 1 } ^ { - 1 } ) + S _ { 1 } ( \Sigma _ { 1 } ^ { - 1 } ) U _ { 1 } ( \Sigma _ { 1 } ^ { - 1 } ) ^ { - 1 } S _ { 1 } ( \Sigma _ { 1 } ^ { - 1 } ) ^ { \top } \leq 0$$

as well as

$$L _ { 1 1 } ( \Sigma _ { 1 } ) \leq C _ { 1 } ^ { \top } C _ { 1 }.$$

This means that Σ 1 is a Gramian of the balanced reduced system (16) .

Proof. By the definiteness criterion of Schur, the first inequality in (17) is equivalent to

$$\begin{bmatrix} - L ( \Sigma ^ { - 1 } ) & S ( \Sigma ^ { - 1 } ) \\ S ( \Sigma ^ { - 1 } ) ^ { \top } & U ( \Sigma ^ { - 1 } ) \end{bmatrix} \geq 0.$$

We multiply this inequality with [ ] v z /latticetop from the left and with [ ] v z from the right, where v ∈ R n and z ∈ R m yielding

$$- v ^ { \top } L ( \Sigma ^ { - 1 } ) v + 2 v ^ { \top } S ( \Sigma ^ { - 1 } ) z + z ^ { \top } U ( \Sigma ^ { - 1 } ) z \geq 0.$$

We use the partition in (15) and the definition of the reduced operators. Now, setting v = [ v 1 0 ] , where v 1 ∈ R r , we obtain

$$& \quad \cdot \, \bar { \ } l \, \bar { \ } l } \\ 0 & \leq - v _ { 1 } ^ { \top } \left ( L _ { 1 1 } ( \Sigma _ { 1 } ^ { - 1 } ) + \sum _ { i, j = 1 } ^ { q } N _ { i, 2 1 } ^ { \top } \Sigma _ { 2 } ^ { - 1 } N _ { j, 2 1 } k _ { i j } \right ) v _ { 1 } \\ & \quad + 2 v _ { 1 } ^ { \top } \left ( S _ { 1 } ( \Sigma _ { 1 } ^ { - 1 } ) + \sum _ { i, j = 1 } ^ { q } N _ { i, 2 1 } ^ { \top } \Sigma _ { 2 } ^ { - 1 } M _ { j, 2 k _ { i j } } \right ) z + z ^ { \top } \left ( U _ { 1 } ( \Sigma _ { 1 } ^ { - 1 } ) - \sum _ { i, j = 1 } ^ { q } M _ { i, 2 } ^ { \top } \Sigma _ { 2 } ^ { - 1 } M _ { j, 2 k _ { i j } } \right ) z. \\ \quad \text{It remains to show that } - v _ { 1 } ^ { \top } \Sigma _ { 1 } ^ { \top } \cdot \cdot \cdot N _ { i, 2 1 } ^ { \top } \Sigma _ { 2 } ^ { - 1 } N _ { i, 2 k _ { i j } } \right ) z -$$

It remains to show that -v /latticetop 1 ∑ q i,j =1 N /latticetop i, 21 Σ -1 2 N j, 21 k ij v 1 +2 v /latticetop 1 ∑ q i,j =1 N /latticetop i, 21 Σ -1 2 M k z j, 2 ij -z /latticetop ∑ q i,j =1 M /latticetop i, 2 Σ -1 2 M k z j, 2 ij ≤ 0, since

$$- v _ { 1 } ^ { \top } L _ { 1 1 } ( \Sigma _ { 1 } ^ { - 1 } ) v _ { 1 } + 2 v _ { 1 } ^ { \top } S _ { 1 } ( \Sigma _ { 1 } ^ { - 1 } ) z + z ^ { \top } U _ { 1 } ( \Sigma _ { 1 } ^ { - 1 } ) z \geq 0$$

for all v 1 ∈ R r and z ∈ R m is equivalent to (18). In order to show the missing step for (18), we observe that

$$\ v e \, \text{that} \\ \sum _ { i, j = 1 } ^ { q } N _ { i, 2 1 } ^ { \top } \Sigma _ { 2 } ^ { - 1 } N _ { j, 2 1 } k _ { i j } & = N _ { 2 1 } ^ { \top } ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) N _ { 2 1 }, \\ \sum _ { i, j = 1 } ^ { q } N _ { i, 2 1 } ^ { \top } \Sigma _ { 2 } ^ { - 1 } M _ { j, 2 } k _ { i j } & = N _ { 2 1 } ^ { \top } ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) M _ { 2 }, \\ \sum _ { i, j = 1 } ^ { q } M _ { i, 2 } ^ { \top } \Sigma _ { 2 } ^ { - 1 } M _ { j, 2 } k _ { i j } & = M _ { 2 } ^ { \top } ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) M _ { 2 }, \\ [ N _ { 1, 2 1 } ^ { \top } \ \dots \ \, N _ { q, 2 1 } ^ { \top } ] \ \text{and} \, M _ { 2 } ^ { \top } \coloneqq [ M _ { 1, 2 1 } ^ { \top } \ \dots \ \, M _ { q, 2 } ^ { \top } ]. \ \text{Since} \, K \\ \quad. \quad. \quad. \quad.$$

where N /latticetop 21 := [ N /latticetop 1 21 , . . . N /latticetop q, 21 ] and M /latticetop 2 := [ M /latticetop 1 2 , . . . M /latticetop q, 2 ] . Since K = ( k ij ) and Σ 2 are positive (semi-)definite, we have that K ⊗ Σ -1 2 is positive (semi-)definite. Therefore, we obtain

$$& - v _ { 1 } ^ { \top } N _ { 2 1 } ^ { \top } ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) N _ { 2 1 } v _ { 1 } + 2 v _ { 1 } ^ { \top } N _ { 2 1 } ^ { \top } ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) M _ { 2 } z - z ^ { \top } M _ { 2 } ^ { \top } ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) M _ { 2 } z \\ & = - \left \| ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) ^ { \frac { 1 } { 2 } } N _ { 2 1 } v _ { 1 } - ( K \otimes \Sigma _ { 2 } ^ { - 1 } ) ^ { \frac { 1 } { 2 } } M _ { 2 } z \right \| _ { 2 } ^ { 2 } \leq 0$$

finally providing (18). We show (19) using partition (15) when evaluating the left upper block of the second inequality in (17). This yields

$$L _ { 1 1 } ( \Sigma _ { 1 } ) + \sum _ { i, j = 1 } ^ { q } N _ { i, 2 1 } ^ { \top } \Sigma _ { 2 } N _ { j, 2 1 } k _ { i j } \leq - C _ { 1 } ^ { \top } C _ { 1 }.$$

Since ∑ q i,j =1 N /latticetop i, 21 Σ 2 N j, 21 k ij = N /latticetop 21 ( K ⊗ Σ 2 ) N 21 ≥ 0, the proof is complete.

We are ready to prove a bound for the error between the outputs of (1) and (16).

Theorem 3.9. Let y be the output of (1) with x (0) = 0 , for which we assume to have a balanced realization, e.g., the assumptions of Proposition 3.7 hold true. Given the r -dimensional reduced system (16) with output y r and x r (0) = 0 . Then, for all u ∈ L 2 T , we have

$$\| y - y _ { r } \| _ { L _ { T } ^ { 2 } } \leq 2 \gamma ^ { \frac { 1 } { 2 } } \sum _ { k = r + 1 } ^ { n } \sigma _ { k } \| u \| _ { L _ { T } ^ { 2 } }.$$

Proof. Below, we continue making use of the partition in (15) while assuming that our original system is already balanced. We add a zero line to the reduced system (16) yielding

$$d \left [ \begin{smallmatrix} x _ { r } ( t ) \end{smallmatrix} \right ] = & \left [ A \left [ \begin{smallmatrix} x _ { r } ( t ) \end{smallmatrix} \right ] + B u ( t ) - \left [ \begin{smallmatrix} 0 \\ v _ { 0 } ( t ) \end{smallmatrix} \right ] d t + \sum _ { i = 1 } ^ { q } \left [ N _ { i } \left [ \begin{smallmatrix} x _ { r } ( t ) \end{smallmatrix} \right ] + M _ { i } u ( t ) - \left [ \begin{smallmatrix} 0 \\ v _ { i } ( t ) \end{smallmatrix} \right ] \end{smallmatrix} \right ] d w _ { i } ( t ),$$

where v 0 = A 21 x r + B u 2 and v i = N i, 21 x r + M u i, 2 . We introduce x -( ) := t x t ( ) -[ x r ( ) t 0 ] and x + ( ) := t x t ( ) + [ x r ( ) t 0 ] . Combining (1) with (21) gives us

$$d x _ { - } ( t ) = [ A x _ { - } ( t ) + \left [ \begin{smallmatrix} 0 \\ v _ { 0 } ( t ) \end{smallmatrix} \right ] ] d t + \sum _ { \cdot \cdot } ^ { q } \left [ N _ { i } x _ { - } ( t ) + \left [ \begin{smallmatrix} 0 \\ v _ { i } ( t ) \end{smallmatrix} \right ] d w _ { i } ( t ),$$

$$d x _ { - } ( t ) & = [ A x _ { - } ( t ) + \left [ \begin{smallmatrix} 0 \\ v _ { 0 } ^ { 0 } ( t ) \end{smallmatrix} \right ] d t + \sum _ { i = 1 } ^ { q } \left [ N _ { i } x _ { - } ( t ) + \left [ \begin{smallmatrix} 0 \\ v _ { i } ^ { 0 } ( t ) \end{smallmatrix} \right ] d w _ { i } ( t ), \\ d x _ { + } ( t ) & = [ A x _ { + } ( t ) + 2 B u ( t ) - \left [ \begin{smallmatrix} 0 \\ v _ { 0 } ( t ) \end{smallmatrix} \right ] d t + \sum _ { i = 1 } ^ { q } \left [ N _ { i } x _ { + } ( t ) + 2 M _ { i } u ( t ) - \left [ \begin{smallmatrix} 0 \\ v _ { i } ( t ) \end{smallmatrix} \right ] d w _ { i } ( t ). \\ \\ \intertext { v _ { i } \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{smallmatrix} }$$

We begin with the case of Σ 2 = σI , σ &gt; 0, which is the foundation for showing the error bound in full generality. In order to proof the bound, we need to derive equations for

quadratic forms of x -and x + . First, Ito's product rule, see Lemma A.1, yields

$$\text{quadratic forms of } x _ { - } \text{ and } x _ { + }. \text{ First, Ito} ^ { \prime } s \text{ product rule, see Lemma A.1, yields} \\ \mathbb { E } \left [ x _ { - } ( t ) ^ { \top } \Sigma _ { x _ { - } } ( t ) \right ] & = 2 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { - } ( s ) ^ { \top } \Sigma [ A x _ { - } ( s ) + \left [ v _ { 0 } ^ { 0 } ( s ) } \right ] \right ] d s \\ & \quad + \int _ { 0 } ^ { t } \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ \left ( N _ { i } x _ { - } ( s ) + \left [ v _ { i } ^ { 0 } ( s ) \right ] \right ) ^ { \top } \Sigma ( N _ { j } x _ { - } ( s ) + \left [ v _ { j } ^ { 0 } ( s ) } \right ] \right ) \right ] k _ { i j } d s \\ & = \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { - } ( s ) ^ { \top } L ( \Sigma ) x _ { - } ( s ) \right ] d s + 2 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { - } ( s ) ^ { \top } \Sigma \left [ v _ { 0 } ^ { 0 } ( s ) \right ] \right ] d s \\ & \quad + \int _ { 0 } ^ { t } \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ \left ( 2 N _ { i } x _ { - } ( s ) + \left [ v _ { i } ^ { 0 } ( s ) \right ] \right ) ^ { \top } \Sigma \left [ v _ { j } ^ { 0 } ( s ) \right ] \right ] k _ { i j } d s. \\ \text{We see that } x _ { - } ( s ) ^ { \top } \Sigma _ { \substack { \left [ v _ { 0 } ^ { 0 } ( s ) \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline{ \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline \sim \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim} \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \over { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \overline { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim } \\ \over { \sim }$$

We see that x -( s ) /latticetop Σ [ 0 v 0 ( s ) ] = x 2 ( s ) /latticetop Σ 2 v 0 ( s ) and ( 2 N x i -( s ) + [ 0 v i ( s ) ] ) /latticetop Σ [ 0 v j ( s ) ] = (2 [ N i, 21 N i, 22 ] ( x s ( ) -[ x r ( s ) 0 ] )+ v i ( s )) /latticetop Σ 2 v j ( s ) = (2 [ N i, 21 N i, 22 ] x s ( ) -v i ( s )+2 M u s i, 2 ( )) /latticetop Σ 2 v j ( s ). Consequently, inserting (9) into the above identify for E [ x -( ) t /latticetop Σ x -( ) t ] leads to

$$\mathbb { E } \left [ x _ { - } ( t ) ^ { \top } \Sigma x _ { - } ( t ) \right ] & \leq \int _ { 0 } ^ { t } - \mathbb { E } \left \| C x _ { - } ( s ) \right \| _ { 2 } ^ { 2 } d s + 2 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { 2 } ( s ) ^ { \top } \Sigma _ { 2 } v _ { 0 } ( s ) \right ] d s \\ & \quad + \int _ { 0 } ^ { t } \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ ( 2 \left [ \, N _ { i, 2 1 } \ N _ { i, 2 2 } \, \right ] x ( s ) - v _ { i } ( s ) + 2 M _ { i, 2 } u ( s ) ) ^ { \top } \Sigma _ { 2 } v _ { j } ( s ) \right ] k _ { i j } d s. \\ \text{Since $Cx$} \ = u - u _ { \cdot } \, \text{we obtain that}$$

Since Cx -= y -y r , we obtain that

$$\L = \dots \\ \| y - y _ { r } \| _ { L _ { t } ^ { 2 } } ^ { 2 } & \leq 2 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { 2 } ( s ) ^ { \top } \Sigma _ { 2 } v _ { 0 } ( s ) \right ] d s \\ & \quad + \int _ { 0 } ^ { t } \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ ( 2 \left [ N _ { i, 2 1 } \, N _ { i, 2 2 } \right ] x ( s ) - v _ { i } ( s ) + 2 M _ { i, 2 } u ( s ) ) ^ { \top } \Sigma _ { 2 } v _ { j } ( s ) \right ] k _ { i j } d s. \\ \intertext { W o n h r t h o a c e m m t i o n } \Sigma _ { \lambda } - \tau \, \text{and avnliit that } \Sigma ^ { q } \quad \text{and} \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \t }$$

We apply the assumption Σ 2 = σI and exploit that ∑ q i,j =1 v i ( s ) /latticetop v j ( s k ) ij ≥ 0. Hence,

$$\| y - y _ { r } \| _ { L _ { t } ^ { 2 } } ^ { 2 } & \leq 2 \sigma \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { 2 } ( s ) ^ { \top } v _ { 0 } ( s ) \right ] d s \\ & \quad + \sigma \int _ { 0 } ^ { t } \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ ( 2 \left [ N _ { i, 2 1 } \ N _ { i, 2 2 } \right ] x ( s ) + v _ { i } ( s ) + 2 M _ { i, 2 } u ( s ) ) ^ { \top } v _ { j } ( s ) \right ] k _ { i j } d s. \\ \text{We aim to combine inequality (22) with a second one that we obtain from the next step}$$

We aim to combine inequality (22) with a second one that we obtain from the next step. This is, finding an estimate for a quadratic form of x + . Ito's product rule (compare with Lemma A.1) leads to

$$\mathbb { E } \left [ x _ { + } ( t ) ^ { \top } \Sigma ^ { - 1 } x _ { + } ( t ) \right ] = 2 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } \Sigma ^ { - 1 } [ A x _ { + } ( s ) + 2 B u ( s ) - \left [ \begin{smallmatrix} 0 \\ v _ { 0 } ( s ) \end{smallmatrix} \right ] \right ] d s$$

$$& M. \, R e d m a n \\ & + \int _ { i j = 1 } ^ { t } \sum _ { \substack { \alpha \\ \alpha } } \mathbb { E } \left [ ( N _ { i x + } ( s ) + 2 M _ { i } u ( s ) - \left [ \sum _ { \substack { \alpha \\ \alpha } } ^ { 0 } \right ] \right ) ^ { \top } \Sigma ^ { - 1 } ( N _ { j x - } ( s ) + 2 M _ { i } u ( s ) - \left [ \sum _ { j } ^ { 0 } \right ] ) \right ] k _ { i j } d s \\ & = \int _ { i j = 1 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } L ( \Sigma ^ { - 1 } ) x _ { + } ( s ) \right ] d s \\ & \quad - 2 \int _ { i j = 1 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } \Sigma ^ { - 1 } \left [ \sum _ { \substack { \alpha \\ \alpha } } ^ { 0 } ( s ) \right ] \right ] d s + 4 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } \Sigma ^ { - 1 } B u ( s ) \right ] d s \\ & \quad + \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } L ( \Sigma ^ { - 1 } ) x _ { + } ( s ) \right ] d s + 4 \int _ { 0 } ^ { t } \mathbb { E } \left [ u ( s ) ^ { \top } \sum _ { i, j = 1 } ^ { \alpha } M _ { i } ^ { t } \Sigma ^ { - 1 } M _ { j } k _ { i j } u ( s ) \right ] d s \\ & \quad + 4 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } \sum _ { i, j = 1 } ^ { \alpha } N _ { i } ^ { t } \Sigma ^ { - 1 } M _ { j } k _ { i j } u ( s ) \right ] d s + 4 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } \Sigma ^ { - 1 } B u ( s ) \right ] d s \\ & \quad - \int _ { i, j = 1 } ^ { t } \mathbb { E } \left [ ( 2 N _ { i x + } ( s ) + 4 M _ { i } u ( s ) - \left [ \sum _ { i } ^ { 0 } \right ] \right ) ^ { \top } \Sigma ^ { - 1 } \left [ \sum _ { i } ^ { 0 } \left ] \right ] k _ { i j } d s \\ & \quad - 2 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { + } ( s ) ^ { \top } \Sigma ^ { - 1 } \left [ \sum _ { \substack { \alpha \\ \alpha } } ^ { 0 } \right ] \right ] d s. \\ & \quad \text{We observe that } x _ { + } ( s ) ^ { \top } \Sigma ^ { - 1 } \left [ \sum _ { \substack { \alpha \\ \alpha } } ^ { 0 } \right ] = \, x _ { 2 } ( s ) ^ { \top } \Sigma ^ { - 1 } v _ { 0 } ( s ) \text{ and } \, ( 2 N _ { i x + } ( s ) + 4 M _ { i } u ( s ) - \left [ \sum _ { i } ^ { 0 } \right ] \right ] \text{$\alpha$} - 1 $,\ell \colon = ( 2 N _ { i x + } ( s ) + \left [ x _ { + } ( s ) \right ] \, \vdash \, \mathcal { + } M _ { i } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { v } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { u } \colon \vec { w }$$

We observe that x + ( s ) /latticetop Σ -1 [ 0 v 0 ( s ) ] = x 2 ( s ) /latticetop Σ -1 2 v 0 ( s ) and (2 N x i + ( s ) + 4 Mu s i ( ) -[ 0 v i ( s ) ] ) /latticetop Σ -1 [ 0 v i ( s ) ] = (2[ N i, 21 N i, 22 ] ( x s ( ) + [ x r ( s ) 0 ] ) + 4 M u s i, 2 ( ) -v i ( s )) /latticetop Σ -1 2 v j ( s ) = (2 [ N i, 21 N i, 22 ] x s ( ) + v i ( s ) + 2 M u s i, 2 ( )) /latticetop Σ -1 2 v j ( s ). Using the definitions of the operators U and S , we therefore have

$$\mathbb { E } \left [ x _ { + } ( t ) ^ { \top } \Sigma ^ { - 1 } x _ { + } ( t ) \right ] & = \mathbb { E } \int _ { 0 } ^ { t } x _ { + } ( s ) ^ { \top } L ( \Sigma ^ { - 1 } ) x _ { + } ( s ) - 4 u ( s ) ^ { \top } U ( \Sigma ^ { - 1 } ) u ( s ) + 4 x _ { + } ( s ) ^ { \top } S ( \Sigma ^ { - 1 } ) u ( s ) d s \\ & \quad - \int _ { 0 } ^ { t } \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ ( 2 \left [ N _ { i, 2 1 } \ N _ { i, 2 2 } \right ] x ( s ) + v _ { i } ( s ) + 2 M _ { i, 2 u } ( s ) ) ^ { \top } \Sigma ^ { - 1 } _ { 2 } v _ { j } ( s ) \right ] k _ { i j } d s \\ & \quad - 2 \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { 2 } ( s ) ^ { \top } \Sigma ^ { - 1 } _ { 2 } v _ { 0 } ( s ) \right ] d s + 4 \gamma \int _ { 0 } ^ { t } \mathbb { E } \| u ( s ) \| _ { 2 } ^ { 2 } d s.$$

We multiply (20) with [ x + ( s ) -2 u s ( ) ] /latticetop from the left and with [ x + ( s ) -2 u s ( ) ] from the right in order to obtain -x + ( s ) /latticetop L (Σ -1 ) x + ( s ) -4 x + ( s ) /latticetop S (Σ -1 ) u s ( ) + 4 u s ( ) /latticetop U (Σ -1 ) u s ( ) ≥ 0.

This yields

$$- \dots \cdots \cdots \\ 4 \gamma \| u \| _ { L ^ { 2 } _ { t } } ^ { 2 } & \geq \sigma ^ { - 1 } \int _ { 0 } ^ { t } \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ \left ( 2 \left [ \, N _ { i, 2 1 } \ N _ { i, 2 2 } \, \right ] x ( s ) + v _ { i } ( s ) + 2 M _ { i, 2 } u ( s ) \right ) ^ { \top } v _ { j } ( s ) \right ] k _ { i j } d s \\ & \quad + 2 \sigma ^ { - 1 } \int _ { 0 } ^ { t } \mathbb { E } \left [ x _ { 2 } ( s ) ^ { \top } v _ { 0 } ( s ) \right ] d s \\ \text{exploiting that $\Sigma_{2}=\sigma I$. Combining this estimate with (22), we see that}$$

exploiting that Σ 2 = σI . Combining this estimate with (22), we see that

$$\| y - y _ { r } \| _ { L ^ { 2 } _ { t } } \leq 2 \sigma \gamma ^ { \frac { 1 } { 2 } } \| u \| _ { L ^ { 2 } _ { t } }.$$

Let us turn our attention to the general case. We can bound the L 2 T -error of the dimension reduction procedure as follows

$$\| y - y _ { r } \| _ { L ^ { 2 } _ { T } } & \leq \| y - y _ { n - 1 } \| _ { L ^ { 2 } _ { T } } + \sum _ { i = r + 1 } ^ { n - 1 } \| y _ { k } - y _ { k - 1 } \| _ { L ^ { 2 } _ { T } } \,, \\ \cdot \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \end{cases}$$

where y k is the output of the reduced model with dimension k ∈ { r, . . . , n - } 1 . By (23) with t = T , we obtain ‖ y -y n -1 ‖ L 2 T ≤ 2 σ γ n 1 2 ‖ u ‖ L 2 T . In the next reduction step from y n -1 to y n -2 we exploit Proposition 3.8. It tells us that y n -1 is the output of a balanced system with Gramian diag( σ , . . . , σ 1 n -1 ). Therefore, we can apply (23) once more. We continue this procedure until we arrive at the reduced system of order r and find in every step that ‖ y k -y k -1 ‖ L 2 T ≤ 2 σ γ k 1 2 ‖ u ‖ L 2 T concluding this proof.

Theorem 3.9 is one of our main results and will be the basis for the error estimates that we provide for the approaches presented in Subsections 2.1 and 2.2. It contains an a-priori error bound that links the truncated HSVs to the error of the dimension reduction. Let us now turn our attention to a model reduction technique that is needed for the ansatz of Subsection 2.2.

## 4 Model reduction for (1) with u ≡ 0

We proceed with the uncontrolled case as the second ansatz of Section 2.2 requires a dimension reduction approach for (2) with ˜ = A A and ˜ N i = N i . We begin with the definition of a reachability Gramian having a very different structure than the one in Definition 3.1.

Definition 4.1. Given that u ≡ 0 , a reachability Gramian for (1) is a symmetric matrix P ≥ 0 satisfying

$$L ^ { * } ( \mathbf P ) + X _ { 0 } X _ { 0 } ^ { \top } = 0,$$

where L ∗ is the adjoint operator of L given by L ∗ ( X ) := AX + XA /latticetop + ∑ q i,j =1 N XN i /latticetop j k ij for some generic matrix X of suitable dimension.

/negationslash

Using the terminology created in [6], P has the form of a type I Gramian. This is, e.g., used in [5, 7] in a control system context, whereas here we have u ≡ 0 and x 0 = 0. Let us show why it is relevant in out setting.

Proposition 4.2. Let P be as in Definition 4.1 and x be the mean square asymptotically stable solution of (1a) with u ≡ 0 . Moreover, let ( p k , λ k ) be an associated basis of eigenvectors for R n with the eigenvalues λ k ≥ 0 . Then,

$$\int _ { 0 } ^ { \infty } \mathbb { E } \left | \langle x ( t ), \mathbf p _ { k } \rangle _ { 2 } \right | ^ { 2 } d t \leq \lambda _ { k } \left \| v \right \| _ { 2 } ^ { 2 }. \\ \dots \left. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right. \right.$$

Proof. Let Φ be the fundamental solution of (1a). We obtain that

$$\int _ { 0 } ^ { \infty } \mathbb { E } \left | \langle x ( t ), \mathbf p _ { k } \rangle _ { 2 } \right | ^ { 2 } d t & = \int _ { 0 } ^ { \infty } \mathbb { E } \left | \langle \Phi ( t ) X _ { 0 } v, \mathbf p _ { k } \rangle _ { 2 } \right | ^ { 2 } d t = \int _ { 0 } ^ { \infty } \mathbb { E } \left | \langle v, X _ { 0 } ^ { \top } \Phi ( t ) ^ { \top } \mathbf p _ { k } \rangle _ { 2 } \right | ^ { 2 } d t \\ & \leq \mathbf p _ { k } ^ { \top } \int _ { 0 } ^ { \infty } \mathbb { E } \left [ \Phi ( t ) X _ { 0 } X _ { 0 } ^ { \top } \Phi ( t ) ^ { \top } \right ] d t \, \mathbf p _ { k } \left \| v \right \| _ { 2 } ^ { 2 }. \\ \text{It is a well-known fact that } \mathbf P \coloneqq \int _ { 0 } ^ { \infty } \mathbb { E } \left [ \Phi ( t ) X _ { 0 } X _ { 0 } ^ { \top } \Phi ( t ) ^ { \top } \right ] d t \, \text{solves } ( 2 4 ), \, \sec, \, e. g., \, [ 8 ].$$

It is a well-known fact that P := ∫ ∞ 0 E [ Φ( ) t X X 0 /latticetop 0 Φ( ) t /latticetop ] dt solves (24), see, e.g., [8]. Therefore, we find that p /latticetop k ∫ ∞ 0 E [ Φ( ) t X X 0 /latticetop 0 Φ( ) t /latticetop ] dt p k = λ k which concludes the proof.

Remark 4.3. Let ¯ P be a solution of the inequality in (24) . Then, we have L ∗ ( ¯ P -P ) = -Y ≤ 0 for some Y ≥ 0 . Using the stochastic representation of this matrix equation, we see that ¯ P -P ≥ 0 allowing to replace P by ¯ P in Proposition 4.2.

When u ≡ 0, there is no change in the arguments used below Definition 3.5. Therefore, the observability Gramians are the same regardless of having an uncontrolled setting or not. However, for simplicity we assume that we do not use a more general inequality but rather the solution to the equation. For that reason, we consider the following definition.

Definition 4.4. Given that u ≡ 0 , an observability Gramian of system (1) is a symmetric matrix Q ≥ 0 satisfying

$$L ( \mathbf Q ) + C ^ { \top } C = 0.$$

We aim to diagonalize P and Q simultaneously for the same reason explained in Section 3.2. The associated state space transformation is here denoted by T defining the balanced state variable x b = T x satisfying

$$d x _ { b } ( t ) & = \mathbf A _ { b } x _ { b } ( t ) d t + \sum _ { i = 1 } ^ { q } \mathbf N _ { i, b } x _ { b } ( t ) d w _ { i } ( t ), & ( 2 6 \mathbf a ) \\ y ( t ) & = \mathbf C _ { b } x _ { b } ( t ), \quad \mathbf x _ { b } ( 0 ) = \mathbf X _ { 0. b } v, \quad t \in [ 0, T ], & ( 2 6 b )$$

$$y ( t ) = \mathbf C _ { b } \mathbf x _ { b } ( t ), \ \mathbf x _ { b } ( 0 ) = \mathbf X _ { 0, b } v, \ \ t \in [ 0, T ],$$

with the same output as (1) ( u ≡ 0) and coefficients

$$( \mathbf A _ { b }, \mathbf N _ { i, b }, \mathbf C _ { b }, \mathbf X _ { 0, b } ) \coloneqq ( \mathbf T A \mathbf T ^ { - 1 }, \mathbf T N _ { i } \mathbf T ^ { - 1 }, \mathbf C \mathbf T ^ { - 1 }, \mathbf T X _ { 0 } ).$$

Given P Q , &gt; 0, the balancing transformation T exists and has exactly the same structure like T in Proposition 3.7. Just the underlying Gramians are different. The structure is identical because the Gramians of (26) are of the form P b = TPT /latticetop and Q b = T -/latticetop QT -1 again. The only difference in the proof is that the result is here obtained by multiplying (24) by T from the left and with T /latticetop from the right. The same type of partition

$$& \mathbf A _ { b } = \begin{bmatrix} \mathbf A _ { 1 1 } & \mathbf A _ { 1 2 } \\ \mathbf A _ { 2 1 } & \mathbf A _ { 2 2 } \end{bmatrix}, \quad \mathbf X _ { 0, b } = \begin{bmatrix} \mathbf X _ { 0, 1 } \\ \mathbf X _ { 0, 2 } \end{bmatrix}, \quad \mathbf N _ { i, b } = \begin{bmatrix} \mathbf N _ { i, 1 1 } & \mathbf N _ { i, 1 2 } \\ \mathbf N _ { i, 2 1 } & \mathbf N _ { i, 2 2 } \end{bmatrix}, \quad \mathbf C _ { b } = \begin{bmatrix} \mathbf C _ { 1 } & \mathbf C _ { 2 } \end{bmatrix} \\ & \text{like in } / 1 5 \, \text{load to the download execution}$$

like in (15) leads to the reduced system

$$d x _ { \mathbf r } ( t ) & = \mathbf A _ { 1 1 } x _ { \mathbf r } ( t ) d t + \sum _ { i = 1 } ^ { q } \mathbf N _ { i, 1 1 } x _ { \mathbf r } ( t ) d w _ { i } ( t ), & ( 2 7 \mathbf a ) \\ \mathbf y _ { \mathbf r } ( t ) & = \mathbf C _ { 1 } x _ { \mathbf r } ( t ), \quad x _ { \mathbf r } ( 0 ) = \mathbf X _ { 0. 1 } v, \quad t \in [ 0, T ]. & ( 2 7 \mathbf b )$$

$$y _ { \mathbf r } ( t ) = \mathbf C _ { 1 } \mathbf x _ { \mathbf r } ( t ), \ \ x _ { \mathbf r } ( 0 ) = \mathbf X _ { 0, 1 } v, \ \ t \in [ 0, T ].$$

Here, x r ( ) t ∈ R r , where r is the number of small HSVs θ i = √ λ i ( PQ ). The dimension of the reduced matrices in (27) is accordingly.

/negationslash

Remark 4.5. It is interesting to notice that the model reduction procedure to lower the dimension of the initial state dynamics is less complex than the one for the dynamics of the control. This is due to the choice of Gramian in each scenario. If x 0 = 0 but the control is present, the algorithm of Section 3 requires the computation of a solution P &gt; 0 satisfying a linear matrix inequality (6) . The scheme in this section, where u ≡ 0 and x 0 = 0 , uses solutions P and Q of matrix equations (24) and (25) . Therefore, it is much easier to determine P and Q . We refer to [9, 20] for potential strategies to solve such matrix equations.

It is not hard to find a first error estimate between the output y of (1) with u ≡ 0 and y r . Using the fundamental matrix Φ r of (27), this is

$$\| y - \mathbf y _ { \mathbf r } \| _ { L ^ { 2 } _ { T } } ^ { 2 } & = \int _ { 0 } ^ { T } \mathbb { E } \, \| y ( t ) - \mathbf y _ { \mathbf r } ( t ) \| _ { 2 } ^ { 2 } \, d t = \int _ { 0 } ^ { T } \mathbb { E } \, \| C \Phi ( t ) X _ { 0 } v - C _ { 1 } \Phi _ { \mathbf r } ( t ) X _ { 0, 1 } v \| _ { 2 } ^ { 2 } \, d t \\ & \leq \int _ { 0 } ^ { \infty } \mathbb { E } \, \| C \Phi ( t ) X _ { 0 } v - C _ { 1 } \Phi _ { \mathbf r } ( t ) X _ { 0, 1 } \| _ { F } ^ { 2 } \, d t \, \| v \| _ { 2 } ^ { 2 } \,. \\ \nu _ { 1 } \, \dots \dots \dots \nu _ { 2 } \, \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{array}.$$

It remains to analyze the integral ∫ ∞ 0 E ‖ C Φ( ) t X v 0 -C Φ 1 r ( ) t X 0 1 , ‖ 2 F dt which is finite due to the mean square asymptotic stability of (1). This was done in [7, Section 4.2], where the integral has been expressed as a term depending on the matrix Θ 2 = diag( θ r +1 , . . . , θ n ) of truncated HSVs. We exploit this result in the following theorem.

Theorem 4.6. Let y be the output of (1) with u ≡ 0 and given that the Gramians P and Q associated to this uncontrolled case are positive definite. Suppose that y r is the

output of the r -dimensional reduced system (27) with x r (0) = 0 . Then, for all u ∈ L 2 T , we have

$$\| y - \mathbf y _ { \mathbf r } \| _ { L _ { T } ^ { 2 } } ^ { 2 } \leq \mathbf r ( \Theta _ { 2 } \mathcal { W } ) \, \| v \| _ { 2 } ^ { 2 } \,,$$

where W = B B 2 /latticetop 2 +2 ˆ P A 2 /latticetop 21 + ∑ q i,j =1 (2 [ N i, 21 N i, 22 ] ˆ P -N i, 21 P r ) N /latticetop j, 21 k ij . The matrices ˆ = P [ ˆ P 1 ˆ P 2 ] and P r are the solutions of

$$A _ { 1 1 } P _ { r } + P _ { r } A _ { 1 1 } ^ { \top } + \sum _ { i, j = 1 } ^ { \prime } N _ { i, 1 1 } P _ { r } N _ { j, 1 1 } ^ { \top } k _ { i j } = - X _ { 0, 1 } X _ { 0, 1 } ^ { \top },$$

$$A _ { b } \hat { P } + \hat { P } A _ { 1 1 } ^ { \top } + \sum _ { i, j = 1 } ^ { \nu } N _ { i, b } \hat { P } N _ { j, 1 1 } ^ { \top } k _ { i j } = - X _ { 0, b } X _ { 0, 1 } ^ { \top }.$$

$$\mathfrak { U } _ { 1 1 } ^ { \top } + \sum _ { i, j = 1 } ^ { q } N _ { i, 1 1 } \mathbf P _ { r } N _ { j, 1 1 } ^ { \top } k _ { i j } & = - X _ { 0, 1 } X _ { 0, 1 } ^ { \top }, \\ \mathbb { P } \mathbf A _ { 1 1 } ^ { \top } + \sum _ { i, j = 1 } ^ { q } N _ { i, b } \hat { \mathbf P } \mathbf N _ { j, 1 1 } ^ { \top } k _ { i j } & = - X _ { 0, b } X _ { 0, 1 } ^ { \top }. \\ \end{cases} \, \text{and in natural to Promotion } 4 6 \text{ in } [ 7 ] \text{ for}$$

Proof. We refer to Section 4.2 and in particular to Proposition 4.6 in [7] for a proof of this statement, where uncorrelated L´vy noise was considered. The proof is identical for e correlated Wiener noise.

Theorem 4.6 contains a bound that involves truncated HSVs and hence links these values to the approximation error. However, in contrast to Theorem 3.9, it is not computable a-priori. We are now ready to discuss our two techniques used for reducing the dimension of (1).

## 5 Model reduction for general systems (1)

We conclude this paper with formulating error bounds for the model reduction strategies proposed in Subsections 2.1 and 2.2. Here, we make use of the results established in Sections 3 and 4.

## 5.1 The approach according to Subsection 2.1

The dimension reduction procedure for (1) can be conducted using the approach of Section 3, which first of all covers the special case of x 0 = 0. However, we can transform (1) into a system with a zero initial state, that is given in (4). The price that needs to be paid is the enlargement of the control space. In fact, the information concerning the initial state is moved to the new control ˜( ) = u t [ u t ( ) u 0 ( ) t v ] . Moreover, all control matrices are updated, meaning that we replace B,M i and D by ˜ B,M ˜ i and ˜ . D We can now apply the model reduction technique of Section 3 to (4). It is based on Gramians ˜ P,Q &gt; 0 satisfying ˜ ( ˜ U P -1 ) &gt; 0 and

$$L ( \tilde { P } ^ { - 1 } ) + \tilde { S } ( \tilde { P } ^ { - 1 } ) \tilde { U } ( \tilde { P } ^ { - 1 } ) ^ { - 1 } \tilde { S } ( \tilde { P } ^ { - 1 } ) ^ { \top } \leq 0 \quad \text{and} \quad L ( Q ) + C ^ { \top } C \leq 0, \quad \text{(28)}$$

where the updated operators are defined by

$$\tilde { S } ( X ) & = X \tilde { B } + \sum _ { i, j = 1 } ^ { q } N _ { i } ^ { \top } X \tilde { M } _ { j } k _ { i j } \quad \text{and} \quad \tilde { U } ( X ) & = \tilde { \gamma } I - \sum _ { i, j = 1 } ^ { q } \tilde { M } _ { i } ^ { \top } X \tilde { M } _ { j } k _ { i j } \quad ( \tilde { \gamma } > 0 ).$$

We determine the balancing transformation ˜ T according to (14) leading to balanced coefficients that we partition as follows

$$\begin{array} { c c c }. \\ \tilde { T } A \tilde { T } ^ { - 1 } = \begin{bmatrix} \tilde { A } _ { 1 1 } & * \\ * & * \end{bmatrix}, & \tilde { T } \tilde { B } = \begin{bmatrix} \tilde { B } _ { 1 } \\ * \end{bmatrix}, & \tilde { T } \tilde { M } _ { i } = \begin{bmatrix} \tilde { M } _ { i, 1 } \\ * \end{bmatrix}, \\ \tilde { T } N _ { i } \tilde { T } ^ { - 1 } = \begin{bmatrix} \tilde { N } _ { i, 1 1 } & * \\ * & * \end{bmatrix}, & C \tilde { T } ^ { - 1 } = \begin{bmatrix} \tilde { C } _ { 1 } & * \end{bmatrix}. \\ \text{quently, the reduced system is of the form}$$

Consequently, the reduced system is of the form

$$d \tilde { x } _ { r } ( t ) & = \left [ \tilde { A } _ { 1 1 } \tilde { x } _ { r } ( t ) + \tilde { B } _ { 1 } \tilde { u } ( t ) \right ] d t + \sum _ { i = 1 } ^ { q } \left [ \tilde { N } _ { i, 1 1 } \tilde { x } _ { r } ( t ) + \tilde { M } _ { i, 1 } \tilde { u } ( t ) \right ] d w _ { i } ( t ), \\ \tilde { y } _ { r } ( t ) & = \tilde { C } _ { 1 } \tilde { x } _ { r } ( t ) + \tilde { D } \tilde { u } ( t ), \quad \tilde { x } _ { r } ( 0 ) = 0.$$

By Theorem 3.9, we know that the error between the original output y in (1b) and ˜ y r can be bound as follows

$$\| y - \tilde { y } _ { r } \| _ { L _ { T } ^ { 2 } } & \leq 2 \tilde { \gamma } ^ { \frac { 1 } { 2 } } \sum _ { k = r + 1 } ^ { n } \tilde { \sigma } _ { k } \| \tilde { u } \| _ { L _ { T } ^ { 2 } }, \\ -$$

where ˜ σ k = √ λ k ( ˜ PQ ) are the HSVs of the transformed system (4). Exploiting that ‖ ˜ u ‖ 2 L 2 T = ‖ u ‖ 2 L 2 T + ‖ u v 0 ‖ 2 L 2 T , we obtain the final error estimate that we formulate in a theorem.

Theorem 5.1. Let y be the the quantity of interest in (1) and given that the Gramian Q is positive definite. Let ˜ y r be the output of the reduced system (29) . Then, we have

$$\| y - \tilde { y } _ { r } \| _ { L ^ { 2 } _ { T } } & \leq 2 \tilde { \gamma } ^ { \frac { 1 } { 2 } } \sum _ { k = r + 1 } ^ { n } \tilde { \sigma } _ { k } \sqrt { \mathbb { E } \int _ { 0 } ^ { T } \| u ( t ) \| _ { 2 } ^ { 2 } d t + \mathbb { E } \int _ { 0 } ^ { T } \| u _ { 0 } ( t ) \| _ { 2 } ^ { 2 } d t \| v \| _ { 2 } ^ { 2 } }, \\ \cdot \quad \tilde { \ } \cdot \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \quad \tilde { \ }. \end{$$

where ˜ σ k , k ∈ { 1 , . . . , n } , are the HSVs of (4) with underlying Gramians solving (28) as well as u ∈ L 2 T and v ∈ R d . Moreover, ˜ γ &gt; 0 is a parameter charaterizing the operator ˜ U and hence the Gramian ˜ P . The process u 0 taking values in R r 0 × d is fixed and describes the underlying low-dimensional subspace of the auxiliary system (2) .

We found an a-priori error estimate for our dimension reduction procedure in Theorem 5.1. In particular, the truncated HSVs ˜ σ r +1 , . . . , σ ˜ n deliver a nice criterion for the choice of a suitable r . Let us briefly recall two choices for ˜ and A ˜ N i in (2) mentioned in

Subsection 2.1. First, we can fix ˜ = A ˜ N i = 0 yielding r 0 = d , V 0 = X 0 and consequently a constant u 0 being the identity of size d × d . Therefore, we obtain E ∫ T 0 ‖ u 0 ( ) t ‖ 2 2 dt = T meaning that the bound is increasing in the terminal time. This illustrates the drawback of exploiting a non-stable system (2), since the error bound explodes as T → ∞ . This is an indicator for a bad approximation quality on long time scales. Secondly, we can use multiples of the identify, i.e., ˜ A = -αI and ˜ N i = γ I i , where γ i ∈ R . The parameter α can actually also be arbitrary, but for simplicity we set α &gt; 0. This yields r 0 = d , V 0 = X 0 and u 0 ( ) = e t ( - -α 0 5 . ∑ q i,j =1 γ γ k i j ij ) + t ∑ q i =1 γ w i i ( ) t I d × d . We calculate that E ‖ u 0 ( ) t ‖ 2 2 = e βt , where β = -2 α + ∑ q i,j =1 γ γ k i j ij . If β = 0, we have the same effect as in the first case. Else, we see that E ∫ T 0 ‖ u 0 ( ) t ‖ 2 2 dt = 1 β (e βT -1). Choosing γ i , so that β &gt; 0 leads to an exponential growth in the error bound and therefore seems to be very non-optimal (unstable system (2)). With β &lt; 0 (mean square asymptotically stable (2)), a small error expression can be ensured and passing to the limit as T → ∞ is possible in Theorem 5.1 if u is square integrable on Ω × [0 , ∞ ), too.

We discuss an alternative to the reduced system (29) in the next section.

## 5.2 The approach according to Subsection 2.2

In contrast to Section 2.1, model reduction needs to applied twice when applying the alternative discussed below. First of all, the control and secondly the initial state dynamics are tackled with a dimension reduction technique. We can simply write the output y of the original system (1) as y = ( y -˜) + ˜, where ˜ is the output of (2) with y y y ˜ = A A and ˜ N i = N i . Now, we can approximate the output ( y -˜) of y (5) by y r associated to the reduced system (16). This approach comes with an error bounded by an expression stated in Theorem 3.9. Moreover, we approximate ˜ by the output y y r of the reduced dynamics in (27). Here, the error is charaterized by Theorem 4.6. In total, we obtain the following error estimate that we formulate in a theorem.

Theorem 5.2. Let y be the quantity of interest in (1) and y r , y r the reduced order outputs of systems (16) and (27) , respectively. We assume that the Gramians P and Q associated to (1) with u ≡ 0 are positive definite. Then, we have

$$\| y - y _ { r } - \mathbf y _ { \mathbf r } \| _ { L _ { T } ^ { 2 } } & \leq 2 \gamma ^ { \frac { 1 } { 2 } } \sum _ { k = r + 1 } ^ { n } \sigma _ { k } \| u \| _ { L _ { T } ^ { 2 } } + ( \mathbf r ( \Theta _ { 2 } W ) ) ^ { \frac { 1 } { 2 } } \left \| v \right \| _ { 2 }, \\ \overbrace { \sim \, \Xi \sim } \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

where σ i = √ λ i ( PQ ) are the HSVs of (1) with x 0 = 0 and Θ = diag( 2 θ r +1 , . . . , θ n ) is the matrix of truncated HSVs θ i = √ λ i ( PQ ) of (1) with u ≡ 0 . The matrix W is a weight specified in Theorem 4.6, γ &gt; 0 charaterizes the operator U , u ∈ L 2 T and v ∈ R d .

Proof. By the triangle inequality, we obtain that

$$\| y - y _ { r } - \mathbf y _ { \mathbf r } \| _ { L ^ { 2 } _ { T } } \leq \| ( y - \tilde { y } ) - y _ { r } \| _ { L ^ { 2 } _ { T } } + \| \tilde { y } - \mathbf y _ { \mathbf r } \| _ { L ^ { 2 } _ { T } } \,,$$

where ˜ is the output of (2) with y ˜ = A A and ˜ N i = N i . The statement follows by applying Theorems 3.9 and 4.6.

We can also pass to the limit as T → ∞ in Theorem 5.2 if lim T →∞ ‖ u ‖ L 2 T exists. Moreover, notice that the reduced dimensions r and r may differ. This makes sense as the potential in reducing the control and the initial state dynamics can be very different. This is an advantage of this approach as, in this regard, it is more flexible than the one of Subsection 2.1. However, a second pair of Gramians P and Q being solutions to (24) and (25) needs to be computed. At least this second pair is easier to compute than P and Q , see Remark 4.5. Another drawback in comparison to the ansatz of Subsection 2.1 is that the bound cannot be computed a-priori. On the other hand, Theorem 5.2 shows nicely how truncated HSVs of control and initial state dynamics are related to the error of the reduction.

## A Supporting lemma

Below, we introduce a lemma that is based on Ito's product rule.

Lemma A.1. Suppose that A,B ,... , B 1 q are R n -valued ( F t ) t ∈ [0 ,T ] -adapted processes with A being almost surely Lebesgue integrable and B i satisfying E ∫ T 0 ‖ B t i ( ) ‖ 2 2 dt &lt; ∞ . Let w = [ w 1 . . . w q ] /latticetop be a Wiener process with covariance matrix K = ( k ij ) . If x is an Ito process given by

$$d x ( t ) & = A ( t ) d t + B ( t ) d w ( t ) = A ( t ) d t + \sum _ { i = 1 } ^ { q } B _ { i } ( t ) d w _ { i } ( t ), \\ \intertext { n \quad \uparrow n } \quad \intertext { n \quad \in 1 \quad n + 1 \cdots \ \dots \ \iota \dots }$$

where B = [ B 1 . . . B q ] . Then, we have

$$\frac { d } { d t } \mathbb { E } \left [ x ( t ) ^ { \top } x ( t ) \right ] & = 2 \mathbb { E } \left [ x ( t ) ^ { \top } A ( t ) \right ] + \sum _ { i, j = 1 } ^ { q } \mathbb { E } \left [ B _ { i } ( t ) ^ { \top } B _ { j } ( t ) \right ] k _ { i j }. \\ \intertext {..... } \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \. & \$$

Proof. Ito's product rule yields d x t ( ( ) /latticetop x t ( ) ) = d x t ( ( ) /latticetop ) x t ( )+ x t ( ) /latticetop d x t ( ( ) ) + d x t ( ( ) /latticetop ) d x t ( ( ) ) . We insert the representation of x into this equation providing

$$d ( x ( t ) ^ { \top } x ( t ) ) & = A ( t ) ^ { \top } x ( t ) d t + ( B ( t ) d w ( t ) ) ^ { \top } x ( t ) + x ( t ) ^ { \top } A ( t ) d t + x ( t ) ^ { \top } B ( t ) d w ( t ) \\ & \quad + ( A ( t ) ^ { \top } d t + \sum _ { i = 1 } ^ { q } B _ { i } ( t ) ^ { \top } d w _ { i } ( t ) ) \left ( A ( t ) d t + \sum _ { j = 1 } ^ { q } B _ { j } ( t ) d w _ { j } ( t ) \right ). \\ \text{The result now follows by the rule } d w _ { i } ( t ) d w _ { i } ( t ) = k _ { i }. d t \text{ and by applying the expected}$$

The result now follows by the rule dw i ( ) t dw j ( ) = t k ij dt and by applying the expected value to the above identify exploiting the fact that the Ito integral has mean zero.

## Acknowledgments

MR is supported by the DFG via the individual grant 'Low-order approximations for large-scale problems arising in the context of high-dimensional PDEs and spatially discretized SPDEs'- project number 499366908.

## References

- [1] U. Baur, P. Benner, and L. Feng. Model Order Reduction for Linear and Nonlinear Systems: A System-Theoretic Perspective. Arch. Comput. Methods Eng. , 21(4):331358, 2014.
- [2] C. Beattie, S. Gugercin, and V. Mehrmann. Model reduction for systems with inhomogeneous initial conditions. Syst. Control Lett. , 99:99-106, 2017.
- [3] S. Becker and C. Hartmann. Infinite-dimensional bilinear and stochastic balanced truncation with error bounds. Math. Control Signals, Syst. , 31:1-37, 2019.
- [4] S. Becker, C. Hartmann, M. Redmann, and L. Richter. Error bounds for model reduction of feedback-controlled linear stochastic dynamics on Hilbert spaces. Stoch. Process. Their Appl. , 149:107-141, 2022.
- [5] P. Benner and T. Damm. Lyapunov equations, energy functionals, and model order reduction of bilinear and stochastic systems. SIAM J. Control Optim. , 49(2):686711, 2011.
- [6] P. Benner, T. Damm, and Y. R. Rodriguez Cruz. Dual pairs of generalized Lyapunov inequalities and balanced truncation of stochastic linear systems. IEEE Trans. Autom. Control , 62(2):782-791, 2017.
- [7] P. Benner and M. Redmann. Model Reduction for Stochastic Systems. Stoch. PDE: Anal. Comp. , 3(3):291-338, 2015.
- [8] T. Damm. Rational Matrix Equations in Stochastic Control. Lecture Notes in Control and Information Sciences 297. Berlin: Springer, 2004.
- [9] T. Damm. Direct methods and ADI-preconditioned Krylov subspace methods for generalized Lyapunov equations. Numer. Linear Algebra Appl. , 15(9):853-871, 2008.
- [10] S. Gugercin, A. C. Antoulas, and C. Beattie. H 2 Model Reduction for Large-Scale Linear Dynamical Systems. SIAM J. Matrix Anal. Appl. , 30(2):609-638, 2008.
- [11] M. Heinkenschloss, T. Reis, and A. C. Antoulas. Balanced truncation model reduction for systems with inhomogeneous initial conditions. Automatica , 47(3):559-564, 2011.
- [12] R. Z. Khasminskii. Stochastic stability of differential equations , volume 66 of Stochastic Modelling and Applied Probability . Springer, Heidelberg, second edition, 2012.
- [13] J. L¨fberg. o YALMIP: A Toolbox for Modeling and Optimization in MATLAB. In In Proceedings of the CACSD Conference , Taipei, Taiwan, 2004.
- [14] X. Mao. Stochastic Differential Equations and Applications (Second Edition) . Woodhead Publishing, 2007.

- [15] B. Moore. Principal component analysis in linear systems: Controllability, observability, and model reduction. IEEE Trans. Autom. Control , 26(1):17-32, 1981.
- [16] MOSEK ApS. The MOSEK optimization toolbox for MATLAB manual. Version 10.0. , 2022.
- [17] M. Redmann. Type II singular perturbation approximation for linear systems with L´ evy noise. SIAM J. Control Optim. , 56(3):2120-2158, 2018.
- [18] M. Redmann and M. A. Freitag. Optimization based model order reduction for stochastic systems. Appl. Math. Comput. , Volume 398, 2021.
- [19] M. Redmann and I. Pontes Duff. Full state approximation by Galerkin projection reduced order models for stochastic and bilinear systems. Appl. Math. Comput. , 490, 2022.
- [20] V. Simoncini. Computational methods for linear matrix equations. SIAM Rev. , 38(3):377-441, 2016.
- [21] M. Voigt and C. Schr¨der. Balanced truncation model reduction with a priori error o bounds for LTI systems with nonzero initial value. J. Comput. Appl. Math. , 420, 2023.
Martin Redmann
Dimension reduction for large-scale stochastic systems with non-zero initial states and controlled diffusion
In this paper, we establish new strategies to reduce the dimension of
large-scale controlled stochastic differential equations with non-zero initial
states. The first approach transforms the original setting into a stochastic
system with zero initial states. This transformation naturally leads to
equations with controlled diffusion. A detailed analysis of dominant subspaces
and bounds for the reduction error is provided in this controlled diffusion
framework. Subsequently, we introduce a reduced system for the original
framework and prove an a-priori error bound for the first ansatz. This bound
involves so-called Hankel singular values that are linked to a new pair of
Gramians. A second strategy is presented that is based on the idea of reducing
control and initial state dynamics separately. Here, different Gramians are
used in order to derive a reduced model and their relation to dominant
subspaces are pointed out. We also show an a posteriori error bound for the
second approach involving two types of Hankel singular values.