Dual pairs of generalized Lyapunov inequalities and balanced truncation of stochastic linear systems
|
|
- Jayson Logan
- 5 years ago
- Views:
Transcription
1 Dual pairs of generalized Lyapunov inequalities and balanced truncation of stochastic linear systems Peter Benner, obias Damm, and Yolanda Rocio Rodriguez Cruz arxiv:54.55v [math.ds 8 Apr 5 Abstract We consider two approaches to balanced truncation of stochastic linear systems, which follow from different generalizations of the reachability Gramian of deterministic systems. Both preserve mean-square asymptotic stability, but only the second leads to a stochastic H -type bound for the approximation error of the truncated system. Index erms generalized Lyapunov equation, model order reduction, balanced truncation, stochastic linear system, asymptotic mean square stability 5A4, 93A5, 93B36, 93B4, 93D5, 93E5, INRODUCION Optimization and (feedback) control of dynamical systems is often computationally infeasible for high dimensional plant models. herefore, one tries to reduce the order of the system, so that the input-output mapping is still computable with sufficient accuracy, but at considerably smaller cost than for the original system, [, [, [3, [4, [5. o guarantee the desired accuracy, computable error bounds are required. Moreover, system properties which are relevant in the context of control system design like asymptotic stability need to be preserved. It has long been known that for linear time-invariant (LI) systems the method of balanced truncation preserves asymptotic stability and provides an error bound for the L - induced input-output norm, that is the H -norm of the associated transfer function, see [6, [7. When considering model order reduction of more general system classes, it is natural to try to extend this approach. his has been worked out for descriptor systems in [8, for time-varying systems in [9, [, [, for bilinear systems in [, [3, [4 and general nonlinear systems e.g. in [5. Yet another generaliztion of LI systems is obtained considering dynamics driven by noise processes. his leads to the class of stochastic systems, which have been considered in a system theoretic context e.g. in [6, [7, [8. Quite recently, balanced truncation has also been described for linear stochastic systems of Itô type in [4, [9, [. Already the formulation of the method leads to two different variants that are equivalent in the deterministic case, but not so for stochastic systems. It is natural to ask which of the above mentioned properties of balanced truncation also hold for these variants. he aim of this paper is to answer this question. P. Benner is with the Max Planck Institute for Dynamics of Complex echnical Systems, Sandtorstr., 396 Magdeburg, Germany benner@mpi-magdeburg.mpg.de.. Damm and Y.R. Rodriguez Cruz are with University of Kaiserslautern, Department of Mathematics, Kaiserslautern, Germany, damm@mathematik.uni-kl.de, rodrigue@mathematik.uni-kl.de Let us first recapitulate balanced truncation for linear deterministic control systems of the form ẋ = Ax + Bu, y = Cx, σ(a) C. () Here A R n n, B R n m, C R p n, and x(t) R n, y(t) R p and u(t) R m are the state, output, and input of the system, respectively. Moreover σ(a) denotes the spectrum of A and C the open left half complex plane. Let L A : X A X + XA denote the Lyapunov operator and L A : X AX + XA its adjoint with respect to the Frobenius inner product. hen σ(a) C if and only if there exists a positive definite solution X of the Lyapunov inequality L A (X) <, by Lyapunov s classical stability theorem, see e.g. [. Balanced truncation means truncating a balanced realization. his realization is obtained by a state space transformation computed from the Gramians P and Q, which solve the dual pair of Lyapunov equations L A (Q) = A Q + QA = C C, L A(P ) = AP + P A = BB, or more generally the inequalities (a) (b) L A (Q) C C, L A(P ) BB. (3) hese (in)equalities are essential in the characterization of stability, controllability and observability of system (). If det P, the inequalities (3) can be written as L A (Q) C C, L A (P ) = P A + A P P BB P. (4a) (4b) In the present paper we discuss extensions of (3) and (4) for stochastic linear systems. As indicated above, the equivalent formulations (3) and (4) lead to different generalizations, if we consider Itô-type stochastic systems of the form dx = Ax dt + Nx dw + Bu dt, y = Cx, (5) where A, B, C are as in () and N R n n. System (5) is asymptotically mean-square stable (e.g. [, [3, [8), if and only if there exists a positive definite solution X of the generalized Lyapunov inequality (L A + Π N )(X) = A X + XA + N XN <.
2 Here Π N : X N XN and Π N : X NXN. his stability criterion indicates that in the stochastic context the generalized Lyapunov operator L A + Π N takes over the role of L A. Substituting L A by L A +Π N in (3) and (4), we obtain two different dual pairs of generalized Lyapunov inequalities. We call them type I: (L A + Π N )(Q) = A Q + QA + N QN C C, (6a) (L A + Π N ) (P ) = AP + P A + NP N BB, (6b) and type II: (L A + Π N )(Q) = A Q + QA + N QN C C, (L A + Π N )(P ) = A P + P A + N P N P BB P. (7a) (7b) Note that (6) corresponds to (3) in the sense that L A (P ) has been replaced by (L A + Π N ) (P ), while (7) corresponds to (4), where L A (P ) has been replaced by (L A + Π N )(P ). In general (if N and P do not commute), the inequalities (6b) and (7b) are not equivalent. At first glance it is not clear which generalization is more appropriate. If the system is asymptotically mean-square stable and certain observability and reachability conditions are fulfilled, then for both types there are solutions Q, P >. By a suitable state space-transformation, it is possible to balance the system such that Q = P = Σ > is diagonal. Consequently, the usual procedure of balanced truncation can be applied to reduce the order of (5). For simplicity, let us refer to this as type I or type II balanced truncation. Under natural assumptions, this reduction preserves meansquare asymptotic stability. For type I, this nontrivial fact has been proven in [4. Moreover, in [, an H -error bound has been provided. However, different from the deterministic case, there is no H -type error bound in terms of the truncated entries in Σ. his will be shown in Example I.3. In contrast, for type II, an H -type error bound has been obtained in [9. In the present paper, as one of our main contributions, we show in heorem II. that type II balanced truncation also preserves mean-square asymptotic stability. he proof differs significantly from the one given for type I. Using this result, we are able to give a more compact proof of the error bound, heorem II.4, which exploits the stochastic bounded real lemma [7. We illustrate our results by analytical and numerical examples in Section IV. I. YPE I BALANCED RUNCAION Consider a stochastic linear control system of Itô-type dx = Ax dt + k N j x dw j + Bu dt, y = Cx, (8) j= where w j = (w j (t)) t R+ are uncorrelated zero mean real Wiener processes on a probability space (Ω, F, µ) with respect to an increasing family (F t ) t R+ of σ-algebras F t F (e.g. [5, [6). o simplify the notation, we only consider the case k = and set w = w, N = N. But all results can immediately be generalized for k >. Let L w(r +, R q ) denote the corresponding space of nonanticipating stochastic processes v with values in R q and norm ( ) v( ) L := E v(t) dt <, w where E denotes expectation. Let the homogeneous equation dx = Ax dt + Nx dw be asymptotically mean-square-stable, i.e. E( x(t) ) t, for all solutions x. hen, by heorem A. the equations A Q + QA + N QN = C C, AP + P A + NP N = BB, have unique solutions Q and P. Under suitable observability and controllability conditions, Q and P are nonsingular. A similarity transformation (A, N, B, C) (S AS, S NS, S B, CS) of the system implies the contragredient transformation as (Q, P ) (S QS, S P S ). Choosing e.g. S = LV Σ /, with Cholesky factorizations LL = P, R R = Q and a singular value decomposition RL = UΣV, we obtain S = Σ / U R and S QS = S P S = Σ = diag(σ,..., σ n ). After suitable partitioning [ Σ Σ =, S = [ S Σ S, S = a truncated system is given in the form [ (A, N, B, C ) = ( AS, NS, B, CS ). he following result has been proven in [4. heorem I. Let A, N R n n satisfy σ(i A + A I + N N) C. For a block-diagonal matrix Σ = diag(σ, Σ ) > with σ(σ ) σ(σ ) =, assume that A Σ + ΣA + N ΣN and AΣ + ΣA + NΣN. hen, with the usual partitioning of A and N, we have σ(i A + A I + N N ) C. Its implication for mean-square stability of the truncated system is immediate. Corollary I. Consider an asymptotically mean square stable stochastic linear system dx = Ax dt + Nx dw.
3 3 Assume that a matrix Σ = diag(σ, Σ ) is given as in heorem I. and A and N are partitioned accordingly. hen the truncated system dx r = A x r dt + N x r dw is also asymptotically mean square stable. If the diagonal entries of Σ are small, it is expected that the truncation error is small. In fact this is supported by an H -error bound obtained in [. Additionally, however, from the deterministic situation (see [6, [), one would also hope for an H -type error bound of the form y y r L w (R +,R p )? α(trace Σ ) u L w (R +,R m ) (9) with some number α >. he following example shows that no such general α exists. [ Example I.3 Let A = a with a >, N = [, B = [, C = [. Solving (6) with equality, we get P = [ 4a a [ 4a, Q = with σ(p Q) = { 8a, 8a 4 } so that Σ = diag(σ, σ ), where σ = 8a and σ = is balanced by the transformation S = 8a [ a /. he system /4. [ hen CS = so that Cr = for the truncated /4 system of order. hus the output of the reduced system is y r, and the truncation error L L r is equal to the stochastic H -norm (see [7) of the original system, L = sup y L w. x()=, u L = w We show now that this norm is equal to a = aσ. hus, depending on a, the ratio of the truncation error and trace Σ = σ can be arbitrarily large. According to the stochastic bounded real lemma, heorem A.5, L is the infimum over all γ so that the Riccati inequality < A X + XA + N XN C C γ XBB X () [ = x + x 3 γ x (a + )x γ x x (a + )x γ x x a x 3 γ x [ x x <. possesses a solution X = x x 3 If a given matrix X satisfies this condition, then so does the same matrix with x replaced by. Hence we can assume that x =, and end up with the two conditions x 3 < a and (after multiplying the upper left entry with γ ) > x + γ x γ x 3 = (x + γ ) γ (γ + x 3 ) > (x + γ ) γ (γ a ). hus necessarily γ > a, i.e. γ > a. his already proves that L a = aσ, which suffices to disprove the existence of a general bound α in (9). aking infima, it is easy to show that indeed L = a. II. YPE II BALANCED RUNCAION We now consider the inequalities (7). Lemma II. Assume that dx = Ax dt + Nx dw is asymptotically mean-square-stable. hen inequality (7b) is solvable with P >. Proof: By heorem A., for a given Y <, there exists a P >, so that A P + P A + N P N = Y. hen P = ε P, for sufficiently small ε >, satisfies A P + P A + N P N = εy < ε P BB P so that (7b) holds even in the strict form. It is easy to see that like in the previous section a state space transformation (A, N, B, C) (S AS, S NS, S B, CS) leads to a contragredient transformation Q S QS, P S P S of the solutions. hat is, Q and P satisfy (7a) and (7b), if and only if S QS and S P S do so for the transformed data. As before, we can assume the system to be balanced with Q = P = Σ = diag(σ I,..., σ ν I) =, () [ Σ Σ where σ >... > σ ν > and σ(σ ) = {σ,..., σ r }, σ(σ ) = {σ r+,..., σ ν }. Hence, we will now assume (after balancing) that a diagonal matrix Σ as in () is given which satisfies A Σ + ΣA + N ΣN C C, A Σ + Σ A + N Σ N Σ BB Σ. Partitioning A, N, B, C like Σ, we write the system as (a) (b) dx = (A x + A x ) dt + (N x + N x ) dw + B u dt dx = (A x + A x ) dt + (N x + N x ) dw + B u dt y = C x + C x. he reduced system obtained by truncation is dx r = A x r + N x r dw + B u dt, y r = C x r. he index r is the number of different singular values σ j that have been kept in the reduced system. In the following subsections, we consider matrices [ [ A A A = N N, N =, A A N N Σ = diag(σ, Σ ) as in (), and equations of the form A Σ + ΣA + N ΣN = C C A Σ + Σ A + N Σ N = B B (3b) (3a) with arbitrary right-hand sides C C and B B.
4 4 For convenience, we write out the blocks of these equations explicitly: A Σ + Σ A + N Σ N A Σ + Σ A + N Σ N A Σ + Σ A + N Σ N = N Σ N C C (4) = N Σ N C C (5) = N Σ N C C (6) A Σ + Σ A + N Σ N = N Σ N B B (7) A Σ + Σ A + N Σ N = N Σ N B B (8) A Σ + Σ A + N Σ N = N Σ N B B (9) A. Preservation of asymptotic stability he following theorem is the main new result of this paper. heorem II. Let A and N be given such that σ(i A + A I + N N) C. () Assume further that for a block-diagonal matrix Σ = diag(σ, Σ ) > with σ(σ ) σ(σ ) =, we have A Σ + ΣA + N ΣN and (a) A Σ + Σ A + N Σ N. hen, with the usual partitioning of A and N, we have (b) σ(i A + A I + N N ) C. () Again we have an immediate interpretation in terms of meansquare stability of the truncated system. Corollary II.3 Consider an asymptotically mean square stable stochastic linear system dx = Ax dt + Nx dw. Assume that a matrix Σ = diag(σ, Σ ) is given as in heorem II. and A and N are partitioned accordingly. hen the truncated system dx r = A x r dt + N x r dw is also asymptotically mean square stable. Proof of heorem II.: Note that the inequalities () are equivalent to the equations (4) (9) with appropriate righthand sides C C and B B. By way of contradiction, we assume that () does not hold. hen by heorem A.3, there exist V, V, α such that A V + V A + N V N = αv. (3) aking the scalar product of the equation (4) with V, we obtain α trace(σ V ) whence α = and C V =, N V = by Corollary A.4. Hence ( A Σ + Σ A + N Σ N ) V =. (4) Analogously, we have B V = by (5). In particular, from N V =, we get ([ ) (L A + Π V N) [ A V + V A = + N V N V A + N V N A V + N V N N V N [ V A =. A V We will show that A V =, which implies σ(i A + A I + N N) (5) in contradiction to (), and thus finishes the proof. We first show that Im V is invariant under A and N. o this end let V z =. hen by (3), = z ( A V + V A + N V N ) z = z N V N z, whence also V N z =, i.e. N z Ker V. From this, we have = ( A V + V A + N V N ) z = V A z, implying A z Ker V. hus A Ker V Ker V and N Ker V Ker V. Since Ker V = (Im V ), it follows further that Im V is invariant under A and N. Let V = V V, where V has full column rank, i.e. det V V. hen by the invariance, there exist square matrices X and Y, such that It follows that A V = V X and N V = V Y. = A V V + V V A + N V V N = V (X + X + Y Y )V, whence X + X + Y Y =. Moreover, from (4), we get A Σ V = Σ A V N Σ N V = Σ V X N Σ V Y. (6) Using this substitution in the following computation, we obtain V Σ ( A Σ + Σ A + NΣ N ) Σ V = V Σ (A Σ V ) + (A Σ V ) Σ V + V Σ NΣ N Σ V = V Σ 3 V X X V Σ 3 V (7) V Σ NΣ V Y Y V Σ N Σ V + V Σ NΣ N Σ V.
5 5 We will show that the right hand side has nonnegative trace. his then implies that the whole term vanishes. Note that trace(y V Σ 3 V Y ) = trace(v Σ 3 V Y Y ) = trace( V Σ 3 V (X + X )) = trace( V Σ 3 V X X V Σ 3 V ). aking the trace in (7), we have ( trace Y V Σ 3 V Y V Σ NΣ V Y ) Y V Σ N Σ V + V Σ NΣ N Σ V [ [ V Y V Y = trace M. V where [ M = V Σ 3 Σ N Σ Σ NΣ Σ NΣ N Σ. he matrix M is positive semidefinite, because the upper left block is positive definite, and the corresponding Schur complement Σ N Σ N Σ Σ N Σ Σ 3 Σ N Σ = vanishes. Hence [ Σ 3 Σ N Σ Σ NΣ Σ NΣ N Σ [ V Y V = implying via the first block row that N Σ V = Σ V Y. From (7), using also (6) again, we thus have = ( A Σ + Σ A + N Σ N ) Σ V = Σ V X N Σ V Y + Σ A Σ V + N Σ V Y = Σ V X + Σ A Σ V, i.e. A Σ V = Σ V X. It follows that for arbitrary k N, the eigenvector V in (3) can be replaced by because Σ k V Σ k = Σ k V V Σ k = Σ ( V X + X + Y Y ) V Σ ( = A Σ V V Σ ) ( + Σ V V Σ ) A ( + N Σ V V Σ ) N. Induction leads to ( = A Σ k V V Σ k ) ( + Σ k V V ( + N Σ k V V Σ k ) N. Σ k ) A As above, we conclude that N Σ k V =, C Σ k V =, and B Σ k V =. Multiplying (5) with Σ (k ) V and (8) with Σ k V, we get A Σ k V + Σ A Σ (k ) V + NΣ k V Y =, A Σ k V + Σ A Σ k V + NΣ k V Y =. Hence (after multiplication with Σ ), for all k, we have Σ A Σ (k ) ( V = Σ A Σ k V + NΣ k V Y ) = A Σ k V. Applying this identity repeatedly, we get A Σ k V = Σ k A V for all k N. If µ is the minimal polynomial of Σ, then σ(σ ) σ(σ ) = implies det µ(σ ) and = A µ(σ )V = µ(σ )A V, whence A V = and also A V =. Hence we obtain the contradiction (5). B. Error estimate he following theorem has been proven in [9 using LMI-techniques. Exploiting the stability result in the previous subsection, we can give a slightly more compact proof based on the stochastic bounded real lemma, heorem A.6. heorem II.4 Let A and N satisfy σ(i A + A I + N N) C. Assume furthermore that for Σ = diag(σ, Σ ) > with Σ = diag(σ r+ I,..., σ ν I) and σ(σ ) σ(σ ) =, the following Lyapunov inequalities hold, A Σ + ΣA + N ΣN C C, A Σ + Σ A + N Σ N Σ BB Σ. If x() = and x r () =, then for all >, it holds that y y r L w ([, ) (σ r σ ν ) u L w ([, ). Proof: We adapt a proof for deterministic systems e.g. [, heorem 7.9. In the central argument we treat the case where Σ = σ ν I and show that y y ν L w [, σ ν u L w [,. (8) From (4) and (7), we can see that also A Σ + Σ A + N Σ N C C, A Σ + Σ A + N Σ N Σ B B Σ. Hence we can repeat the above argument to remove σ ν,..., σ r+ successively. By the triangle inequality we find that ν y y r L w [, y j+ y j L w [, j=r (σ ν σ r+ ) u L w [,. which then concludes the proof. o prove (8), we make use of the stochastic bounded real lemma. In the following let r = ν and consider the error system defined by dx e = A e x e dt + N e x e dw + B e u dt, y e = C e x e = y y r,
6 6 where x A A x e = x, A e = A A, x r A N N B N e = N N, B e = B, N B C e = [ C C C. Applying the state space transformation x x x r = x x r x x + x r = we obtain the transformed system I r I r I n r I r I r } {{ } =S x x x r A A Ã e = S A e S = A A A, A A N N Ñ e = S N e S = N N N, N N B e = S B B, Ce = C e S = [ C C. B By heorem A.6, we have L e σ ν, if the Riccati inequality, satisfies (9). Partitioning R σν (X) = we have R R R 3 R R R 3 R 3 R 3 R 33 R = A Σ + Σ A + N Σ N + σ ν N N + C C = A Σ + Σ A + N Σ N + N Σ N + C C σ ν N N R = A Σ + σ ν A + N Σ N + σ ν N N + C C R 3 = σ ν N N R = σ ν (A + A + NN ) + NΣ N + σνn Σ N + C C + B B = A Σ + Σ A + NΣ N + NΣ N + C C + σν(a Σ + Σ A + NΣ N + NΣ N + Σ B B Σ ) R 3 = σν(σ A + NΣ N ) + σ ν (A + NN ) + σ ν Σ B B = σν(σ A + NΣ N + A Σ + NΣ N + Σ B B Σ ) R 33 = σν(a Σ + Σ A + NΣ N ) + σ ν N N + σνσ B B Σ = σν(a Σ + Σ A + NΣ N + Σ B B Σ + N Σ N ) σ ν With the permutation matrix J = [ I I N N we define M = J(A Σ + Σ A + N Σ N + Σ BB Σ )J, where M by (3b). Using (4) (9), we have [ A R σν (X) = Σ + ΣA + N ΣN + C C σ N N [ ν + σ N N ν, M which is inequality (9)., R γ (X) = Ã e X + XÃe + Ñ e XÑe + C e C e + 4σν X B e B e X (9) possesses a solution X. We will show now that the blockdiagonal matrix X = diag(σ, Σ, σ νσ ) = diag(σ, σ ν I, σ νσ ) > Example II.5 Let the system (A, N, B, C) and Q be as in Example I.3. he matrix [ + p P = >, where < p, p satisfies inequality (7b). As in Example I.3, we have L r = for the corresponding reduced system of order, so that the truncation error again is a, independently of p,. On the other hand we have σ = min σ(p Q) = 4a ( + p) 8a,
7 7 with equality for p. heorem II.4 thus gives the sharp error bound σ = a. Note, that there is no P > satisfying the equation (7b). he previous example illustrates the problem of optimizing over all solutions of inequality (7b). III. NUMERICAL EXAMPLES o compare the reduction methods we need to compute Q, P from (6) or (7). Instead of the inequalities (6a), (6b), (7a) we can consider the corresponding equations, for which quite efficient algorithms have been developed recently, e.g. [7, [8, [9, [3. hese also allow for a low-rank approximation of the solutions. In contrast we cannot replace (7b) by the corresponding equation, because this may not be solvable (see Example II.5). Even worse, we do not have any solvability or uniqueness criteria nor reliable algorithms. herefore, in general, we have to work with the inequality (7b), which is solvable according to Lemma II., but of course not uniquely solvable. In view of our application, we aim at a solution P of (7b), so that (some of) the eigenvalues of P Q are particularly small, since they provide the error bound. Choosing a matrix Y < and a very small ε along the lines of the proof of Lemma II. can be contrary to this aim. Hence some optimization over all solutions of (7b) is required. Note also that a matrix P > satisfies (7b), if and only if it satisfies the linear matrix inequality (LMI) [ P A + AP + BB P N. (3) NP P hus, LMI optimal solution techniques are applicable. However, their complexity will be prohibitive for large-scale problems. herefore further research for alternative methods to solve (7b) adequately is required. By L and L r, we always denote the original and the r- th order approximated system. he stochastic H -type norm L L r is computed by a binary search of the infimum of all γ such that the Riccati inequality () is solvable. he latter is solved via a Newton iteration as in [8. Finally, the Lyapunov equations () are solved by preconditioned Krylov subspace methods described in [7. Unfortunately, for small γ, i.e. for small approximation errors, this method of computing the error runs into numerical problems, because () contains the term γ. his apparently leads to cancellation phenomena in the Newton iteration, if e.g. γ < 7. herefore we mainly concentrate on cases where the error is larger, that is we make r sufficiently small. A. ype II can be better than type I In many examples we observe that type II reduction gives a valid error bound, but the approximation error still is better with type I. his, however, is not always true, as the example (A, N, B, C ) = ([, [, [ 3, [ 3 ) shows. It can easily be verified that the type I Lyapunov equations (6) are solved by [ [ Q = and P = he type II inequalities (7) are e.g. solved by [ [ Q = and P = 3 3. Reduction to order r = gives the following error bounds and approximation errors for both types: σ L L I II As we see, the type I approximation error is larger than both the truncated singular value and the type II approximation error. B. An electrical ladder network with perturbed inductance As our first example with a physical background, we take up the electrical ladder network described in [3, consisting of n/ sections with a capacitor C, inductor L and two resistors R and R as depicted here. V R =. L =. I C =. R = But following e.g. [3, we assume that the inductance L is subject to stochastic perturbations. For simplicity, we replace the inverse L formally by L + ẇ in all sections. Here L =. and ẇ is white noise of a certain intensity σ, where we set σ =. E.g. for n = 6, we have the system matrices A = CR L C R R L(R+ R) R C(R+ R) R L(R+ R) C(R+ R) L C R R L(R+ R) R C(R+ R) R R R R+ R R+ R N = R R R R+ R R+ R R B = [ CR C = [. R R L(R+ R) C(R+ R) L For larger n, the band structure of A and N is extended periodically. o see the behaviour of our two methods, we C R L
8 8 reduce from order n = to the orders r =, 3, 5,..., 9, and compute both the theoretical bounds and the actual approximation errors in the H -norm. he results are shown in the following figure. bound I error I bound II error II As we can see, the upper error bound fails for type I, but is correct for type II. Nevertheless, judging from the H error, neither of the types seems to be preferable over the other. D. Summary Clearly, higher dimensional examples are required to get more insight. o this end a more sophisticated method for the solution of (3) is needed. With general purpose LMI-software on a standard Laptop, we hardly got higher than n = In this example, for both types the bounds hold, and for all reduced orders, type I gives a better approximation than type II. C. A heat transfer problem As another example we consider a stochastic modification of the heat transfer problem described in [4. On the unit square Ω = [, the heat equation x t = x is given with Dirichlet condition x = u j, j =,, 3 on three of the boundary edges and a stochastic Robin condition n x = (/ + ẇ)x on the fourth edge (where ẇ stands for white noise). A standard 5- point finite difference discretization on a grid leads to a modified Poisson matrix A R and corresponding matrices [ N R and B R 3. We use the input u and choose the average temperature as the output, i.e. C = [,...,. We apply balanced truncation of type I and type II. For type II, an LMI-solver (MALAB R function mincx) is used to compute P as a solution of the LMI (3) which minimizes trace P or trace P Q. In the following two figures, we compare the reduced systems of order r = for both types. he left figure shows the decay of the singular values. Since the LMI-solver was called with tolerance level 9, only the first about 5 singular values for type II have the correct order of magnitude. he right figure shows the approximation error y(t) y r (t) over a given time interval. For both types it has the same order of magnitude. In fact, for many examples we have observed both methods to yield very similar results. Decay of singular values 5 5 type type x 5 Output error, n=, r=.5.5 type I type II We have computed the estimated error norm and the actual approximation error for both types: j= σ j L L j= σ j L L I 4.66e 6 9.3e 6.e e 9 II.75e e 6.7e 8 9.7e 9 IV. CONCLUSIONS We have compared two types of balanced truncation for stochastic linear systems, which are related to different Gramian type matrices P and Q. he following table collects properties of these reduction methods. ype I II Def. of P, Q (6) (7) Stability? Yes, [4 Yes, hm. II. H -bound? Yes, [ no result H -bound? No, Ex. I.3 Yes, hm. II.4 or [9 he main contributions of this paper are the preservation of asymptotic stability for type II balanced truncation proved in heorem II. and the new proof of the H error bound in heorem II.4. he efficient solution of (7b) is an open issue and requires further research. he same is true for the computation of the stochastic H -norm. APPENDIX ASYMPOIC MEAN SQUARE SABILIY Consider the stochastic linear system of Itô-type dx = Ax dt + Nx dw, (3) where w = (w(t)) t R+ is a zero mean real Wiener process on a probability space (Ω, F, µ) with respect to an increasing family (F t ) t R+ of σ-algebras F t F (e.g. [5, [6). Let L w(r +, R q ) denote the corresponding space of non-anticipating stochastic processes v with values in R q and norm ( ) v( ) L := E v(t) dt <, w where E denotes expectation. By definition, system (3) is asymptotically mean-square-stable, if E( x(t) ) t, for all initial conditions x() = x. We have the following version of Lyapunov s matrix theorem, see [3. Here denotes the Kronecker product. heorem A. he following are equivalent. (i) System (3) is asymptotically mean-square stable. (ii) max{rλ λ σ(a I + I A + N N)} < (iii) Y > : X > : A X + XA + N XN = Y (iv) Y > : X > : A X + XA + N XN = Y (v) Y : X : A X + XA + N XN = Y
9 9 Remark A. he theorem (like all other results in this paper) carries over to systems dx = Ax dt + k N j x dw j j= with more than one noise term, and many more equivalent criteria can be provided, see e.g. [33 or [8, heorem he following theorem does not require any stability assumptions (see [8, heorem 3..3). It is central in the analysis of mean-square stability. heorem A.3 Let α = max{rλ λ σ(a I + I A + N N)}. hen there exists a nonnegative definite matrix V, such that (L A + Π N )(V ) = AV + V A + NV N = αv. We also note a simple consequence of this theorem [4, Corollary 3.. Here Y, V = trace(y V ) is the Frobenius inner product for symmetric matrices. Corollary A.4 Let α, V as in the theorem. For given Y assume that X > : L A (X) + Π N (X) Y. (3) hen α. Moreover, if α = then Y V = V Y =. HE SOCHASIC BOUNDED REAL LEMMA Now let us consider system (5) with input u and output y. If system (3) is asymptotically mean-square stable, then (5) defines an input output operator L : u y from L w(r, R m ) to L w(r, R p ), see [7. By L we denote the induced operator norm, which is an analogue of the deterministic H -norm. It can be characterized by the stochastic bounded real lemma. heorem A.5 [7 For γ >, the following are equivalent. (i) System (3) is asymptotically mean-square stable and L < γ. (ii) here exists a negative definite solution X < to the Riccati inequality A X + XA + N XN C C γ XBB X >. (iii) here exists a positive definite solution X > to the Riccati inequality A X + XA + N XN + C C + γ XBB X <. We have stated the obviously equivalent formulations (ii) and (iii) to avoid confusion arising from different formulations in the literature. Under additional assumptions also non-strict versions can be formulated. he following sufficient criterion is given in [8, Corollary..3 (where also the signs are changed). Unlike in the previous theorem, here asymptotic mean-square stability is assumed at the outset. heorem A.6 Assume that (3) is asymptotically stable in mean-square. If there exists a nonnegative definite matrix X, satisfying A X + XA + N XN + C C + γ XBB X, then L γ. REFERENCES [ G. Obinata, B. Anderson, Model Reduction for Control System Design, Springer,. [ A. C. Antoulas, Approximation of large-scale dynamical systems, Vol. 6 of Advances in Design and Control, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 5. [3 P. Benner, V. Mehrmann, D. C. Sorensen (Eds.), Dimension Reduction of Large-Scale Systems, Vol. 45 of Lecture Notes in Computational Science and Engineering, Springer-Verlag, 5. [4 W. H. Schilders, H. A. van der Vorst, J. Rommes (Eds.), Model Order Reduction: heory, Research Aspects and Applications, Vol. 3 of Mathematics in Industry, Springer-Verlag, 8. [5 U. Baur, P. Benner, L. Feng, Model order reduction for linear and nonlinear systems: A system-theoretic perspective, Arch. Comput. Method. E. (4) (4) [6 B. C. Moore, Principal component analysis in linear systems: controllability, observability, and model reduction, IEEE rans. Autom. Control AC-6 (98) 7 3. [7 L. Pernebo, L. M. Silverman, Model reduction via balanced state space representations, IEEE rans. Autom. Control AC-7 () (98) [8. Stykel, Analysis and numerical solution of generalized Lyapunov equations, Ph.D. thesis, echnische Universität Berlin (). [9 A. Shokoohi, L. M. Silverman, P. M. Van Dooren, Linear time-variable systems: balancing and model reduction, IEEE rans. Automat. Contr. AC-8 (8) (983) 8 8. [ E. Verriest,. Kailath, On generalized balanced realizations, IEEE rans. Automat. Contr. AC-8 (8) (983) [ H. Sandberg, A. Rantzer, Balanced truncation of linear time-varying systems, IEEE rans. Automat. Contr. 49 () (4) 7 9. [ S. A. Al-Baiyat, M. Bettayeb, U. M. Al-Saggaf, New model reduction scheme for bilinear systems, Int. J. Systems Sci. 5 (994) [3 W. S. Gray, J. Mesko, Energy functions and algebraic Gramians for bilinear systems, in: Preprints of the 4th IFAC Nonlinear Control Systems Design Symposium, Enschede, he Netherlands, 998, pp [4 P. Benner,. Damm, Lyapunov equations, energy functionals, and model order reduction of bilinear and stochastic systems, SIAM J. Control Optim. 49 () () [5 J. M. A. Scherpen, Balancing for nonlinear systems, Syst. Control Lett. () (993) [6 W. M. Wonham, Random differential equations in control theory, in: A.. Bharucha-Reid (Ed.), Probab. Methods Appl. Math., Vol., Academic Press, New York - London, 97, pp. 3. [7 D. Hinrichsen, A. J. Pritchard, Stochastic H, SIAM J. Control Optim. 36 (5) (998) [8. Damm, Rational Matrix Equations in Stochastic Control, no. 97 in Lecture Notes in Control and Information Sciences, Springer, 4. [9. Damm, P. Benner, Balanced truncation for stochastic linear systems with guaranteed error bound, in: Proceedings of MNS-4, Groningen, he Netherlands, 4, pp [ M. Redmann, P. Benner, Model reduction for stochastic systems, Preprint MPIMD/4-3, Max Planck Institute Magdeburg (4). [ F. R. Gantmacher, he heory of Matrices (Vol. II), Chelsea, New York, 959. [ D. L. Kleinman, On the stability of linear stochastic systems, IEEE rans. Autom. Control AC-4 (969) [3 R. Z. Khasminskij, Stochastic Stability of Differential Equations, Sijthoff & Noordhoff, Alphen aan den Rijn, NL, 98. [4 P. Benner,. Damm, M. Redmann, Y. Rocio Rodriguez Cruz, Positive operators and stable truncation, Linear Algebra Appl. doi:.6/j.laa.4..5, in press. published electronically, Dec. 3, 4. [5 L. Arnold, Stochastic Differential Equations: heory and Applications. ranslation., John Wiley and Sons Inc., New York etc., 974. [6 B. Oeksendal, Stochastic Differential Equations, 5th Edition, Springer- Verlag, 998.
10 [7. Damm, Direct methods and ADI-preconditioned Krylov subspace methods for generalized Lyapunov equations, Numer. Lin. Alg. Appl. 5 (9) (8) [8 P. Benner,. Breiten, Low rank methods for a class of generalized Lyapunov equations and related issues., Numer. Math. 4 (3) (3) [9 D. Kressner, P. Sirković, Greedy low-rank methods for solving general linear matrix equations, echnical report, ANCHP, MAHICSE, EPF Lausanne, Switzerland (4). [3 S. Shank, V. Simoncini, D. Szyld, Efficient low-rank solutions of generalized Lyapunov equations, Report 4--, Department of Mathematics emple University, Philadelphia, PA 9 (4). [3 S. Gugercin, A. Antoulas, A survey of model reduction by balanced truncation and some new results, Int. J. Control 77 (8) (4) [3 V. A. Ugrinovskii, I. R. Petersen, Absolute stabilization and minimax optimal control of uncertain systems with stochastic uncertainty, SIAM J. Control Optim. 37 (4) (999) 89. [33 H. Schneider, Positive operators and an inertia theorem, Numer. Math. 7 (965) 7.
Model reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
More informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationBalancing of Lossless and Passive Systems
Balancing of Lossless and Passive Systems Arjan van der Schaft Abstract Different balancing techniques are applied to lossless nonlinear systems, with open-loop balancing applied to their scattering representation.
More informationKrylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations
Krylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations Peter Benner and Tobias Breiten Abstract We discuss Krylov-subspace based model reduction
More informationIterative Solution of a Matrix Riccati Equation Arising in Stochastic Control
Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control Chun-Hua Guo Dedicated to Peter Lancaster on the occasion of his 70th birthday We consider iterative methods for finding the
More informationImplicit Volterra Series Interpolation for Model Reduction of Bilinear Systems
Max Planck Institute Magdeburg Preprints Mian Ilyas Ahmad Ulrike Baur Peter Benner Implicit Volterra Series Interpolation for Model Reduction of Bilinear Systems MAX PLANCK INSTITUT FÜR DYNAMIK KOMPLEXER
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
More informationRow and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems
ROMANIAN JOURNAL OF INFORMATION SCIENCE AND TECHNOLOGY Volume 19, Numbers 1-2, 2016, 127 136 Row and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems Octavian PASTRAVANU,
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationLinear and Nonlinear Matrix Equations Arising in Model Reduction
International Conference on Numerical Linear Algebra and its Applications Guwahati, January 15 18, 2013 Linear and Nonlinear Matrix Equations Arising in Model Reduction Peter Benner Tobias Breiten Max
More informationTechnical Notes and Correspondence
1108 IEEE RANSACIONS ON AUOMAIC CONROL, VOL. 47, NO. 7, JULY 2002 echnical Notes and Correspondence Stability Analysis of Piecewise Discrete-ime Linear Systems Gang Feng Abstract his note presents a stability
More informationISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION
ISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION Harald K. Wimmer 1 The set of all negative-semidefinite solutions of the CARE A X + XA + XBB X C C = 0 is homeomorphic
More informationWeighted balanced realization and model reduction for nonlinear systems
Weighted balanced realization and model reduction for nonlinear systems Daisuke Tsubakino and Kenji Fujimoto Abstract In this paper a weighted balanced realization and model reduction for nonlinear systems
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of
More informationCANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM
CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationw T 1 w T 2. w T n 0 if i j 1 if i = j
Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -
More informationON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION
Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th
More informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationBlock companion matrices, discrete-time block diagonal stability and polynomial matrices
Block companion matrices, discrete-time block diagonal stability and polynomial matrices Harald K. Wimmer Mathematisches Institut Universität Würzburg D-97074 Würzburg Germany October 25, 2008 Abstract
More informationBALANCING AS A MOMENT MATCHING PROBLEM
BALANCING AS A MOMENT MATCHING PROBLEM T.C. IONESCU, J.M.A. SCHERPEN, O.V. IFTIME, AND A. ASTOLFI Abstract. In this paper, we treat a time-domain moment matching problem associated to balanced truncation.
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More informationonly nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr
The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands
More informationThroughout these notes we assume V, W are finite dimensional inner product spaces over C.
Math 342 - Linear Algebra II Notes Throughout these notes we assume V, W are finite dimensional inner product spaces over C 1 Upper Triangular Representation Proposition: Let T L(V ) There exists an orthonormal
More informationarxiv: v1 [math.oc] 29 Apr 2017
BALANCED TRUNCATION MODEL ORDER REDUCTION FOR QUADRATIC-BILINEAR CONTROL SYSTEMS PETER BENNER AND PAWAN GOYAL arxiv:175.16v1 [math.oc] 29 Apr 217 Abstract. We discuss balanced truncation model order reduction
More informationClustering-Based Model Order Reduction for Multi-Agent Systems with General Linear Time-Invariant Agents
Max Planck Institute Magdeburg Preprints Petar Mlinarić Sara Grundel Peter Benner Clustering-Based Model Order Reduction for Multi-Agent Systems with General Linear Time-Invariant Agents MAX PLANCK INSTITUT
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH
ON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH V. FABER, J. LIESEN, AND P. TICHÝ Abstract. Numerous algorithms in numerical linear algebra are based on the reduction of a given matrix
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More informationQUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY
QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics
More informationPeter C. Müller. Introduction. - and the x 2 - subsystems are called the slow and the fast subsystem, correspondingly, cf. (Dai, 1989).
Peter C. Müller mueller@srm.uni-wuppertal.de Safety Control Engineering University of Wuppertal D-497 Wupperta. Germany Modified Lyapunov Equations for LI Descriptor Systems or linear time-invariant (LI)
More informationOn a quadratic matrix equation associated with an M-matrix
Article Submitted to IMA Journal of Numerical Analysis On a quadratic matrix equation associated with an M-matrix Chun-Hua Guo Department of Mathematics and Statistics, University of Regina, Regina, SK
More informationOrder reduction numerical methods for the algebraic Riccati equation. V. Simoncini
Order reduction numerical methods for the algebraic Riccati equation V. Simoncini Dipartimento di Matematica Alma Mater Studiorum - Università di Bologna valeria.simoncini@unibo.it 1 The problem Find X
More information3118 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 57, NO. 12, DECEMBER 2012
3118 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 57, NO 12, DECEMBER 2012 ASimultaneousBalancedTruncationApproachto Model Reduction of Switched Linear Systems Nima Monshizadeh, Harry L Trentelman, SeniorMember,IEEE,
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationSUCCESSIVE POLE SHIFTING USING SAMPLED-DATA LQ REGULATORS. Sigeru Omatu
SUCCESSIVE POLE SHIFING USING SAMPLED-DAA LQ REGULAORS oru Fujinaka Sigeru Omatu Graduate School of Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Sakai, 599-8531 Japan Abstract: Design of sampled-data
More informationCramér-Rao Bounds for Estimation of Linear System Noise Covariances
Journal of Mechanical Engineering and Automation (): 6- DOI: 593/jjmea Cramér-Rao Bounds for Estimation of Linear System oise Covariances Peter Matiso * Vladimír Havlena Czech echnical University in Prague
More informationMemoryless output feedback nullification and canonical forms, for time varying systems
Memoryless output feedback nullification and canonical forms, for time varying systems Gera Weiss May 19, 2005 Abstract We study the possibility of nullifying time-varying systems with memoryless output
More informationA Perron-type theorem on the principal eigenvalue of nonsymmetric elliptic operators
A Perron-type theorem on the principal eigenvalue of nonsymmetric elliptic operators Lei Ni And I cherish more than anything else the Analogies, my most trustworthy masters. They know all the secrets of
More informationContour integral solutions of Sylvester-type matrix equations
Contour integral solutions of Sylvester-type matrix equations Harald K. Wimmer Mathematisches Institut, Universität Würzburg, 97074 Würzburg, Germany Abstract The linear matrix equations AXB CXD = E, AX
More informationStability and Inertia Theorems for Generalized Lyapunov Equations
Published in Linear Algebra and its Applications, 355(1-3, 2002, pp. 297-314. Stability and Inertia Theorems for Generalized Lyapunov Equations Tatjana Stykel Abstract We study generalized Lyapunov equations
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 14: LMIs for Robust Control in the LF Framework ypes of Uncertainty In this Lecture, we will cover Unstructured,
More informationThe Bock iteration for the ODE estimation problem
he Bock iteration for the ODE estimation problem M.R.Osborne Contents 1 Introduction 2 2 Introducing the Bock iteration 5 3 he ODE estimation problem 7 4 he Bock iteration for the smoothing problem 12
More informationc 1998 Society for Industrial and Applied Mathematics
SIAM J MARIX ANAL APPL Vol 20, No 1, pp 182 195 c 1998 Society for Industrial and Applied Mathematics POSIIVIY OF BLOCK RIDIAGONAL MARICES MARIN BOHNER AND ONDŘEJ DOŠLÝ Abstract his paper relates disconjugacy
More informationH 2 -optimal model reduction of MIMO systems
H 2 -optimal model reduction of MIMO systems P. Van Dooren K. A. Gallivan P.-A. Absil Abstract We consider the problem of approximating a p m rational transfer function Hs of high degree by another p m
More informationZero controllability in discrete-time structured systems
1 Zero controllability in discrete-time structured systems Jacob van der Woude arxiv:173.8394v1 [math.oc] 24 Mar 217 Abstract In this paper we consider complex dynamical networks modeled by means of state
More informationLINEAR-QUADRATIC OPTIMAL CONTROL OF DIFFERENTIAL-ALGEBRAIC SYSTEMS: THE INFINITE TIME HORIZON PROBLEM WITH ZERO TERMINAL STATE
LINEAR-QUADRATIC OPTIMAL CONTROL OF DIFFERENTIAL-ALGEBRAIC SYSTEMS: THE INFINITE TIME HORIZON PROBLEM WITH ZERO TERMINAL STATE TIMO REIS AND MATTHIAS VOIGT, Abstract. In this work we revisit the linear-quadratic
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2017 LECTURE 5
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 17 LECTURE 5 1 existence of svd Theorem 1 (Existence of SVD) Every matrix has a singular value decomposition (condensed version) Proof Let A C m n and for simplicity
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationA Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations
A Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations Peter Benner Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory
More informationOHSx XM511 Linear Algebra: Solutions to Online True/False Exercises
This document gives the solutions to all of the online exercises for OHSx XM511. The section ( ) numbers refer to the textbook. TYPE I are True/False. Answers are in square brackets [. Lecture 02 ( 1.1)
More informationThe Newton-ADI Method for Large-Scale Algebraic Riccati Equations. Peter Benner.
The Newton-ADI Method for Large-Scale Algebraic Riccati Equations Mathematik in Industrie und Technik Fakultät für Mathematik Peter Benner benner@mathematik.tu-chemnitz.de Sonderforschungsbereich 393 S
More informationA Symmetry Approach for Balanced Truncation of Positive Linear Systems
A Symmetry Approach for Balanced Truncation of Positive Linear Systems Grussler, Christian; Damm, Tobias Published in: IEEE Control Systems Society Conference 2012, Maui 2012 Link to publication Citation
More informationOptimization Based Output Feedback Control Design in Descriptor Systems
Trabalho apresentado no XXXVII CNMAC, S.J. dos Campos - SP, 017. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics Optimization Based Output Feedback Control Design in
More informationj=1 x j p, if 1 p <, x i ξ : x i < ξ} 0 as p.
LINEAR ALGEBRA Fall 203 The final exam Almost all of the problems solved Exercise Let (V, ) be a normed vector space. Prove x y x y for all x, y V. Everybody knows how to do this! Exercise 2 If V is a
More informationModel Order Reduction of Continuous LTI Large Descriptor System Using LRCF-ADI and Square Root Balanced Truncation
, July 1-3, 15, London, U.K. Model Order Reduction of Continuous LI Large Descriptor System Using LRCF-ADI and Square Root Balanced runcation Mehrab Hossain Likhon, Shamsil Arifeen, and Mohammad Sahadet
More informationLecture notes: Applied linear algebra Part 1. Version 2
Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and
More informationUpper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza
Aalborg Universitet Upper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza Published in: 7th IFAC Symposium on Robust Control Design DOI (link
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
More informationLinear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,
More informationROBUST STABILITY TEST FOR UNCERTAIN DISCRETE-TIME SYSTEMS: A DESCRIPTOR SYSTEM APPROACH
Latin American Applied Research 41: 359-364(211) ROBUS SABILIY ES FOR UNCERAIN DISCREE-IME SYSEMS: A DESCRIPOR SYSEM APPROACH W. ZHANG,, H. SU, Y. LIANG, and Z. HAN Engineering raining Center, Shanghai
More informationPrediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate
www.scichina.com info.scichina.com www.springerlin.com Prediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate WEI Chen & CHEN ZongJi School of Automation
More information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationAn Observation on the Positive Real Lemma
Journal of Mathematical Analysis and Applications 255, 48 49 (21) doi:1.16/jmaa.2.7241, available online at http://www.idealibrary.com on An Observation on the Positive Real Lemma Luciano Pandolfi Dipartimento
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationKTH. Access to the published version may require subscription.
KTH This is an accepted version of a paper published in IEEE Transactions on Automatic Control. This paper has been peer-reviewed but does not include the final publisher proof-corrections or journal pagination.
More informationLMI based output-feedback controllers: γ-optimal versus linear quadratic.
Proceedings of the 17th World Congress he International Federation of Automatic Control Seoul Korea July 6-11 28 LMI based output-feedback controllers: γ-optimal versus linear quadratic. Dmitry V. Balandin
More informationSolvability of Linear Matrix Equations in a Symmetric Matrix Variable
Solvability of Linear Matrix Equations in a Symmetric Matrix Variable Maurcio C. de Oliveira J. William Helton Abstract We study the solvability of generalized linear matrix equations of the Lyapunov type
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
More informationALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA
ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA Kent State University Department of Mathematical Sciences Compiled and Maintained by Donald L. White Version: August 29, 2017 CONTENTS LINEAR ALGEBRA AND
More informationWhere is matrix multiplication locally open?
Linear Algebra and its Applications 517 (2017) 167 176 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Where is matrix multiplication locally open?
More informationThe Important State Coordinates of a Nonlinear System
The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative
More informationStat 159/259: Linear Algebra Notes
Stat 159/259: Linear Algebra Notes Jarrod Millman November 16, 2015 Abstract These notes assume you ve taken a semester of undergraduate linear algebra. In particular, I assume you are familiar with the
More informationOn Positive Real Lemma for Non-minimal Realization Systems
Proceedings of the 17th World Congress The International Federation of Automatic Control On Positive Real Lemma for Non-minimal Realization Systems Sadaaki Kunimatsu Kim Sang-Hoon Takao Fujii Mitsuaki
More informationWe first repeat some well known facts about condition numbers for normwise and componentwise perturbations. Consider the matrix
BIT 39(1), pp. 143 151, 1999 ILL-CONDITIONEDNESS NEEDS NOT BE COMPONENTWISE NEAR TO ILL-POSEDNESS FOR LEAST SQUARES PROBLEMS SIEGFRIED M. RUMP Abstract. The condition number of a problem measures the sensitivity
More informationMatrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein
Matrix Mathematics Theory, Facts, and Formulas with Application to Linear Systems Theory Dennis S. Bernstein PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Special Symbols xv Conventions, Notation,
More informationMatrix Inequalities by Means of Block Matrices 1
Mathematical Inequalities & Applications, Vol. 4, No. 4, 200, pp. 48-490. Matrix Inequalities by Means of Block Matrices Fuzhen Zhang 2 Department of Math, Science and Technology Nova Southeastern University,
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationThe model reduction algorithm proposed is based on an iterative two-step LMI scheme. The convergence of the algorithm is not analyzed but examples sho
Model Reduction from an H 1 /LMI perspective A. Helmersson Department of Electrical Engineering Linkoping University S-581 8 Linkoping, Sweden tel: +6 1 816 fax: +6 1 86 email: andersh@isy.liu.se September
More informationNetwork Reconstruction from Intrinsic Noise: Non-Minimum-Phase Systems
Preprints of the 19th World Congress he International Federation of Automatic Control Network Reconstruction from Intrinsic Noise: Non-Minimum-Phase Systems David Hayden, Ye Yuan Jorge Goncalves Department
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
More informationVector and Matrix Norms. Vector and Matrix Norms
Vector and Matrix Norms Vector Space Algebra Matrix Algebra: We let x x and A A, where, if x is an element of an abstract vector space n, and A = A: n m, then x is a complex column vector of length n whose
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
More informationSTA141C: Big Data & High Performance Statistical Computing
STA141C: Big Data & High Performance Statistical Computing Numerical Linear Algebra Background Cho-Jui Hsieh UC Davis May 15, 2018 Linear Algebra Background Vectors A vector has a direction and a magnitude
More informationRank-one LMIs and Lyapunov's Inequality. Gjerrit Meinsma 4. Abstract. We describe a new proof of the well-known Lyapunov's matrix inequality about
Rank-one LMIs and Lyapunov's Inequality Didier Henrion 1;; Gjerrit Meinsma Abstract We describe a new proof of the well-known Lyapunov's matrix inequality about the location of the eigenvalues of a matrix
More informationSome New Results on Lyapunov-Type Diagonal Stability
Some New Results on Lyapunov-Type Diagonal Stability Mehmet Gumus (Joint work with Dr. Jianhong Xu) Department of Mathematics Southern Illinois University Carbondale 12/01/2016 mgumus@siu.edu (SIUC) Lyapunov-Type
More informationSimultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer
Preprints of the 19th World Congress The International Federation of Automatic Control Simultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer Fengming Shi*, Ron J.
More informationCanonical lossless state-space systems: staircase forms and the Schur algorithm
Canonical lossless state-space systems: staircase forms and the Schur algorithm Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics School of Mathematical Sciences Projet APICS Universiteit
More informationHankel Optimal Model Reduction 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and
More information