BALANCING-RELATED MODEL REDUCTION FOR LARGE-SCALE SYSTEMS
|
|
- Florence Norman
- 5 years ago
- Views:
Transcription
1 FOR LARGE-SCALE SYSTEMS Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Recent Advances in Model Order Reduction Workshop TU Eindhoven, November 23, 2007
2 Overview 1 Model Reduction for Linear Systems Application Areas Goals 2 The Basic Ideas Balancing-Related Model Reduction Numerical Algorithms for 3 Solving Large-Scale Matrix ADI Method for Lyapunov Newton s Method for AREs 4 of Complexity < O(n 3 ) 5 Systems with a Large Number of Terminals 6 7 8
3 Linear Systems Application Areas Goals Dynamical Systems with Σ : states x(t) R n, inputs u(t) R m, outputs y(t) R p. { ẋ(t) = f (t, x(t), u(t)), x(t0 ) = x 0, y(t) = g(t, x(t), u(t)),
4 Model Reduction for Dynamical Systems Linear Systems Application Areas Goals Original System Σ : j ẋ(t) = f (t, x(t), u(t)), y(t) = g(t, x(t), u(t)). states x(t) R n, inputs u(t) R m, outputs y(t) R p. Reduced-Order System bσ : j ˆx(t) = b f (t, ˆx(t), u(t)), ŷ(t) = bg(t, ˆx(t), u(t)). states ˆx(t) R r, r n inputs u(t) R m, outputs ŷ(t) R p. Goal: y ŷ < tolerance u for all admissible input signals.
5 Model Reduction for Dynamical Systems Linear Systems Application Areas Goals Original System Σ : j ẋ(t) = f (t, x(t), u(t)), y(t) = g(t, x(t), u(t)). states x(t) R n, inputs u(t) R m, outputs y(t) R p. Reduced-Order System bσ : j ˆx(t) = b f (t, ˆx(t), u(t)), ŷ(t) = bg(t, ˆx(t), u(t)). states ˆx(t) R r, r n inputs u(t) R m, outputs ŷ(t) R p. Goal: y ŷ < tolerance u for all admissible input signals.
6 Model Reduction for Dynamical Systems Linear Systems Application Areas Goals Original System Σ : j ẋ(t) = f (t, x(t), u(t)), y(t) = g(t, x(t), u(t)). states x(t) R n, inputs u(t) R m, outputs y(t) R p. Reduced-Order System bσ : j ˆx(t) = b f (t, ˆx(t), u(t)), ŷ(t) = bg(t, ˆx(t), u(t)). states ˆx(t) R r, r n inputs u(t) R m, outputs ŷ(t) R p. Goal: y ŷ < tolerance u for all admissible input signals.
7 Model Reduction for Linear Systems Linear Systems Application Areas Goals Linear, Time-Invariant (LTI) Systems f (t, x, u) = Ax + Bu, A R n n, B R n m, g(t, x, u) = Cx + Du, C R p n, D R p m. Linear Systems Frequency Domain Application of Laplace transformation (x(t) x(s), ẋ(t) sx(s)) to linear system with x(0) = 0: sx(s) = Ax(s) + Bu(s), yields I/O-relation in frequency domain: ( y(s) = G is the transfer function of Σ. C(sI n A) 1 B + D }{{} =:G(s) y(s) = Bx(s) + Du(s), ) u(s)
8 Model Reduction for Linear Systems Linear Systems Application Areas Goals Linear, Time-Invariant (LTI) Systems f (t, x, u) = Ax + Bu, A R n n, B R n m, g(t, x, u) = Cx + Du, C R p n, D R p m. Linear Systems Frequency Domain Application of Laplace transformation (x(t) x(s), ẋ(t) sx(s)) to linear system with x(0) = 0: sx(s) = Ax(s) + Bu(s), yields I/O-relation in frequency domain: ( y(s) = G is the transfer function of Σ. C(sI n A) 1 B + D }{{} =:G(s) y(s) = Bx(s) + Du(s), ) u(s)
9 Model Reduction for Linear Systems Linear Systems Application Areas Goals Linear, Time-Invariant (LTI) Systems f (t, x, u) = Ax + Bu, A R n n, B R n m, g(t, x, u) = Cx + Du, C R p n, D R p m. Linear Systems Frequency Domain Application of Laplace transformation (x(t) x(s), ẋ(t) sx(s)) to linear system with x(0) = 0: sx(s) = Ax(s) + Bu(s), yields I/O-relation in frequency domain: ( y(s) = G is the transfer function of Σ. C(sI n A) 1 B + D }{{} =:G(s) y(s) = Bx(s) + Du(s), ) u(s)
10 Model Reduction for Linear Systems Linear Systems Application Areas Goals Problem Approximate the dynamical system ẋ = Ax + Bu, A R n n, B R n m, y = Cx + Du, C R p n, D R p m. by reduced-order system ˆx = ˆx + ˆBu,  R r r, ˆB R r m, ŷ = Ĉ ˆx + ˆDu, Ĉ R p r, ˆD R p m, of order r n, such that y ŷ = Gu Ĝu G Ĝ u < tolerance u. = Approximation problem: min order ( Ĝ) r G Ĝ.
11 Model Reduction for Linear Systems Linear Systems Application Areas Goals Problem Approximate the dynamical system ẋ = Ax + Bu, A R n n, B R n m, y = Cx + Du, C R p n, D R p m. by reduced-order system ˆx = ˆx + ˆBu,  R r r, ˆB R r m, ŷ = Ĉ ˆx + ˆDu, Ĉ R p r, ˆD R p m, of order r n, such that y ŷ = Gu Ĝu G Ĝ u < tolerance u. = Approximation problem: min order ( Ĝ) r G Ĝ.
12 Application Areas Linear Systems Application Areas Goals Feedback Control controllers designed by LQR/LQG, H 2, H methods are LTI systems of order n, but technological implementation needs order 10. Optimization/open-loop control time-discretization of already large-scale systems leads to huge number of equality constraints in mathematical program. Microelectronics verification of VLSI/ULSI chip design requires high number of simulations for different input signals, various effects due to progressive miniaturization lead to large-scale systems of differential(-algebraic) equations (order 10 8 ). MEMS/Microsystem design smart system integration needs compact models for efficient coupled simulation....
13 Goals Linear Systems Application Areas Goals Automatic generation of compact models. Satisfy desired error tolerance for all admissible input signals, i.e., want y ŷ < tolerance u u L 2 (R, R m ). = Need computable error bound/estimate! Preserve physical properties: stability (poles of G in C, i.e., Λ (A) C ), minimum phase (zeroes of G in C ), passivity: t u(τ)t y(τ) dτ 0 t R, u L 2 (R, R m ). ( system does not generate energy ).
14 Goals Linear Systems Application Areas Goals Automatic generation of compact models. Satisfy desired error tolerance for all admissible input signals, i.e., want y ŷ < tolerance u u L 2 (R, R m ). = Need computable error bound/estimate! Preserve physical properties: stability (poles of G in C, i.e., Λ (A) C ), minimum phase (zeroes of G in C ), passivity: t u(τ)t y(τ) dτ 0 t R, u L 2 (R, R m ). ( system does not generate energy ).
15 Goals Linear Systems Application Areas Goals Automatic generation of compact models. Satisfy desired error tolerance for all admissible input signals, i.e., want y ŷ < tolerance u u L 2 (R, R m ). = Need computable error bound/estimate! Preserve physical properties: stability (poles of G in C, i.e., Λ (A) C ), minimum phase (zeroes of G in C ), passivity: t u(τ)t y(τ) dτ 0 t R, u L 2 (R, R m ). ( system does not generate energy ).
16 Goals Linear Systems Application Areas Goals Automatic generation of compact models. Satisfy desired error tolerance for all admissible input signals, i.e., want y ŷ < tolerance u u L 2 (R, R m ). = Need computable error bound/estimate! Preserve physical properties: stability (poles of G in C, i.e., Λ (A) C ), minimum phase (zeroes of G in C ), passivity: t u(τ)t y(τ) dτ 0 t R, u L 2 (R, R m ). ( system does not generate energy ).
17 Goals Linear Systems Application Areas Goals Automatic generation of compact models. Satisfy desired error tolerance for all admissible input signals, i.e., want y ŷ < tolerance u u L 2 (R, R m ). = Need computable error bound/estimate! Preserve physical properties: stability (poles of G in C, i.e., Λ (A) C ), minimum phase (zeroes of G in C ), passivity: t u(τ)t y(τ) dτ 0 t R, u L 2 (R, R m ). ( system does not generate energy ).
18 Goals Linear Systems Application Areas Goals Automatic generation of compact models. Satisfy desired error tolerance for all admissible input signals, i.e., want y ŷ < tolerance u u L 2 (R, R m ). = Need computable error bound/estimate! Preserve physical properties: stability (poles of G in C, i.e., Λ (A) C ), minimum phase (zeroes of G in C ), passivity: t u(τ)t y(τ) dτ 0 t R, u L 2 (R, R m ). ( system does not generate energy ).
19 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Idea: A system Σ, realized by (A, B, C, D), is called balanced, if solutions P, Q of the Lyapunov equations AP + PA T + BB T = 0, A T Q + QA + C T C = 0, satisfy: P = Q = diag(σ 1,..., σ n ) with σ 1 σ 2... σ n > 0. {σ 1,..., σ n } are the Hankel singular values (HSVs) of Σ. Compute balanced realization of the system via state-space transformation T : (A, B, C, D) (TAT 1, TB, CT 1, D)»» A11 A 12 B1 =, A 21 A 22 B 2 (Â, ˆB, Ĉ, ˆD) = (A 11, B 1, C 1, D)., ˆ C 1 C 2, D «
20 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Idea: A system Σ, realized by (A, B, C, D), is called balanced, if solutions P, Q of the Lyapunov equations AP + PA T + BB T = 0, A T Q + QA + C T C = 0, satisfy: P = Q = diag(σ 1,..., σ n ) with σ 1 σ 2... σ n > 0. {σ 1,..., σ n } are the Hankel singular values (HSVs) of Σ. Compute balanced realization of the system via state-space transformation T : (A, B, C, D) (TAT 1, TB, CT 1, D)»» A11 A 12 B1 =, A 21 A 22 B 2 (Â, ˆB, Ĉ, ˆD) = (A 11, B 1, C 1, D)., ˆ C 1 C 2, D «
21 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Idea: A system Σ, realized by (A, B, C, D), is called balanced, if solutions P, Q of the Lyapunov equations AP + PA T + BB T = 0, A T Q + QA + C T C = 0, satisfy: P = Q = diag(σ 1,..., σ n ) with σ 1 σ 2... σ n > 0. {σ 1,..., σ n } are the Hankel singular values (HSVs) of Σ. Compute balanced realization of the system via state-space transformation T : (A, B, C, D) (TAT 1, TB, CT 1, D)»» A11 A 12 B1 =, A 21 A 22 B 2 (Â, ˆB, Ĉ, ˆD) = (A 11, B 1, C 1, D)., ˆ C 1 C 2, D «
22 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Idea: A system Σ, realized by (A, B, C, D), is called balanced, if solutions P, Q of the Lyapunov equations AP + PA T + BB T = 0, A T Q + QA + C T C = 0, satisfy: P = Q = diag(σ 1,..., σ n ) with σ 1 σ 2... σ n > 0. {σ 1,..., σ n } are the Hankel singular values (HSVs) of Σ. Compute balanced realization of the system via state-space transformation T : (A, B, C, D) (TAT 1, TB, CT 1, D)»» A11 A 12 B1 =, A 21 A 22 B 2 (Â, ˆB, Ĉ, ˆD) = (A 11, B 1, C 1, D)., ˆ C 1 C 2, D «
23 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Motivation: HSV are system invariants: they are preserved under T and determine the energy transfer given by the Hankel map H : L 2 (, 0) L 2 (0, ) : u y +.
24 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Motivation: HSV are system invariants: they are preserved under T and determine the energy transfer given by the Hankel map H : L 2 (, 0) L 2 (0, ) : u y +. In balanced coordinates... energy transfer from u to y + : E := sup u L 2 (,0] x(0)=x 0 y(t) T y(t) dt 0 0 u(t) T u(t) dt = 1 x 0 2 n σj 2 x0,j 2 j=1
25 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Motivation: HSV are system invariants: they are preserved under T and determine the energy transfer given by the Hankel map H : L 2 (, 0) L 2 (0, ) : u y +. In balanced coordinates... energy transfer from u to y + : E := sup u L 2 (,0] x(0)=x 0 y(t) T y(t) dt 0 0 u(t) T u(t) dt = 1 x 0 2 n σj 2 x0,j 2 j=1 = Truncate states corresponding to small HSVs = complete analogy to best approximation via SVD!
26 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Implementation: SR Method 1 Compute Cholesky factors of the solutions of the Lyapunov equations, P = S T S, Q = R T R. 2 Compute SVD (thanks, Gene!) 3 Set SR T = [ U 1, U 2 ] [ Σ1 Σ 2 ] [ V T 1 V T 2 ]. W = R T V 1 Σ 1/2 1, V = S T U 1 Σ 1/ Reduced model is (W T AV, W T B, CV, D).
27 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Implementation: SR Method 1 Compute Cholesky factors of the solutions of the Lyapunov equations, P = S T S, Q = R T R. 2 Compute SVD (thanks, Gene!) 3 Set SR T = [ U 1, U 2 ] [ Σ1 Σ 2 ] [ V T 1 V T 2 ]. W = R T V 1 Σ 1/2 1, V = S T U 1 Σ 1/ Reduced model is (W T AV, W T B, CV, D).
28 The Basic Ideas The Basic Ideas Balancing- Related MR Numerical Algorithms for Implementation: SR Method 1 Compute Cholesky factors of the solutions of the Lyapunov equations, P = S T S, Q = R T R. 2 Compute SVD (thanks, Gene!) 3 Set SR T = [ U 1, U 2 ] [ Σ1 Σ 2 ] [ V T 1 V T 2 ]. W = R T V 1 Σ 1/2 1, V = S T U 1 Σ 1/ Reduced model is (W T AV, W T B, CV, D).
29 Balancing-Related Model Reduction The Basic Ideas Balancing- Related MR Numerical Algorithms for Basic Principle of Given positive semidefinite matrices P = S T S, Q = R T R, compute balancing state-space transformation so that P = Q = diag(σ 1,..., σ n ) = Σ, σ 1... σ n 0, and truncate corresponding realization at size r with σ r > σ r+1.
30 Balancing-Related Model Reduction The Basic Ideas Balancing- Related MR Numerical Algorithms for Basic Principle of Given positive semidefinite matrices P = S T S, Q = R T R, compute balancing state-space transformation so that P = Q = diag(σ 1,..., σ n ) = Σ, σ 1... σ n 0, and truncate corresponding realization at size r with σ r > σ r+1. Classical (BT) Mullis/Roberts 76, Moore 81 P = controllability Gramian of system given by (A, B, C, D). Q = observability Gramian of system given by (A, B, C, D). P, Q solve dual Lyapunov equations AP + PA T + BB T = 0, A T Q + QA + C T C = 0.
31 Balancing-Related Model Reduction The Basic Ideas Balancing- Related MR Numerical Algorithms for Basic Principle of Given positive semidefinite matrices P = S T S, Q = R T R, compute balancing state-space transformation so that P = Q = diag(σ 1,..., σ n ) = Σ, σ 1... σ n 0, and truncate corresponding realization at size r with σ r > σ r+1. LQG (LQGBT) Jonckheere/Silverman 83 P/Q = controllability/observability Gramian of closed-loop system based on LQG compensator. P, Q solve dual algebraic Riccati equations (AREs) 0 = AP + PA T PC T CP + B T B, 0 = A T Q + QA QBB T Q + C T C.
32 Balancing-Related Model Reduction The Basic Ideas Balancing- Related MR Numerical Algorithms for Basic Principle of Given positive semidefinite matrices P = S T S, Q = R T R, compute balancing state-space transformation so that P = Q = diag(σ 1,..., σ n ) = Σ, σ 1... σ n 0, and truncate corresponding realization at size r with σ r > σ r+1. Stochastic (BST) Desai/Pal 84, Green 88 P = controllability Gramian of system given by (A, B, C, D), i.e., solution of Lyapunov equation AP + PA T + BB T = 0. Q = observability Gramian of right spectral factor of power spectrum of system given by (A, B, C, D), i.e., solution of ARE  T Q + Q + QB W (DD T ) 1 B T W Q + C T (DD T ) 1 C = 0, where  := A B W (DD T ) 1 C, B W := BD T + PC T.
33 Balancing-Related Model Reduction The Basic Ideas Balancing- Related MR Numerical Algorithms for Basic Principle of Given positive semidefinite matrices P = S T S, Q = R T R, compute balancing state-space transformation so that P = Q = diag(σ 1,..., σ n ) = Σ, σ 1... σ n 0, and truncate corresponding realization at size r with σ r > σ r+1. Positive-Real (PRBT) Green 88 Based on positive-real equations, related to positive real (Kalman-Yakubovich-Popov-Anderson) lemma. P, Q solve dual AREs 0 = ĀP + PĀ T + PC T R 1 CP + B R 1 B T, 0 = ĀT Q + QĀ + QB R 1 B T Q + C T R 1 C, where R = D + D T, Ā = A B R 1 C.
34 Balancing-Related Model Reduction The Basic Ideas Balancing- Related MR Numerical Algorithms for Basic Principle of Given positive semidefinite matrices P = S T S, Q = R T R, compute balancing state-space transformation so that P = Q = diag(σ 1,..., σ n ) = Σ, σ 1... σ n 0, and truncate corresponding realization at size r with σ r > σ r+1. Other Balancing-Based Methods Bounded-real balanced truncation (BRBT) based on bounded real lemma [Opdenacker/Jonckheere 88]; H balanced truncation (HinfBT) closed-loop balancing based on H compensator [Mustafa/Glover 91]. Both approaches require solution of dual AREs. Frequency-weighted versions of the above approaches.
35 Balancing-Related Model Reduction Properties The Basic Ideas Balancing- Related MR Numerical Algorithms for Guaranteed preservation of physical properties like stability (all), passivity (PRBT), minimum phase (BST). Computable error bounds, e.g., BT: G G r 2 LQGBT: G G r 2 BST: G G r n j=r+1 n j=r+1 n j=r+1 σ BT j, σ LQG j q 1+(σ LQG j ) 2 1+σ BST j 1 σ BST j 1 G, Can be combined with singular perturbation approximation for steady-state performance. Computations can be modularized.
36 Numerical Algorithms for General misconception: complexity O(n 3 ) true for several implementations! (e.g., Matlab, SLICOT). The Basic Ideas Balancing- Related MR Numerical Algorithms for
37 Numerical Algorithms for General misconception: complexity O(n 3 ) true for several implementations! (e.g., Matlab, SLICOT). Algorithmic ideas from numerical linear algebra (since 1997): The Basic Ideas Balancing- Related MR Numerical Algorithms for
38 Numerical Algorithms for The Basic Ideas Balancing- Related MR Numerical Algorithms for General misconception: complexity O(n 3 ) true for several implementations! (e.g., Matlab, SLICOT). Algorithmic ideas from numerical linear algebra (since 1997): Instead of Gramians P, Q compute S, R R n k, k n, such that P SS T, Q RR T. Compute S, R with problem-specific Lyapunov/Riccati solvers of low complexity directly.
39 Numerical Algorithms for The Basic Ideas Balancing- Related MR Numerical Algorithms for General misconception: complexity O(n 3 ) true for several implementations! (e.g., Matlab, SLICOT). Algorithmic ideas from numerical linear algebra (since 1997): Instead of Gramians P, Q compute S, R R n k, k n, such that P SS T, Q RR T. Compute S, R with problem-specific Lyapunov/Riccati solvers of low complexity directly. need solver for large-scale matrix equations which computes S, R directly!
40 Solving Large-Scale Matrix Large-Scale Algebraic Lyapunov and Riccati ADI for Lyapunov Newton s Method for AREs General form for A, G = G T, W = W T R n n given and X R n n unknown: 0 = L(X ) := A T X + XA + W, 0 = R(X ) := A T X + XA XGX + W. In large scale applications, typically n = (= unknowns!), A has sparse representation (A = M 1 K for FEM), G, W low-rank with G, W {BB T, C T C}, where B R n m, m n, C R p n, p n. Standard (eigenproblem-based) O(n 3 ) methods are not applicable!
41 Solving Large-Scale Matrix Large-Scale Algebraic Lyapunov and Riccati ADI for Lyapunov Newton s Method for AREs General form for A, G = G T, W = W T R n n given and X R n n unknown: 0 = L(X ) := A T X + XA + W, 0 = R(X ) := A T X + XA XGX + W. In large scale applications, typically n = (= unknowns!), A has sparse representation (A = M 1 K for FEM), G, W low-rank with G, W {BB T, C T C}, where B R n m, m n, C R p n, p n. Standard (eigenproblem-based) O(n 3 ) methods are not applicable!
42 ADI Method for Lyapunov ADI for Lyapunov Newton s Method for AREs For A R n n stable, B R n m (w n), consider Lyapunov equation AX + XA T = BB T. ADI Iteration: [Wachspress 88] (A + p k I )X (j 1)/2 = BB T X k 1 (A T p k I ) (A + p k I )X k T = BB T X (j 1)/2 (A T p k I ) with parameters p k C and p k+1 = p k if p k R. For X 0 = 0 and proper choice of p k : lim k X k = X superlinear. Re-formulation using X k = Y k Y T k yields iteration for Y k...
43 ADI Method for Lyapunov ADI for Lyapunov Newton s Method for AREs For A R n n stable, B R n m (w n), consider Lyapunov equation AX + XA T = BB T. ADI Iteration: [Wachspress 88] (A + p k I )X (j 1)/2 = BB T X k 1 (A T p k I ) (A + p k I )X k T = BB T X (j 1)/2 (A T p k I ) with parameters p k C and p k+1 = p k if p k R. For X 0 = 0 and proper choice of p k : lim k X k = X superlinear. Re-formulation using X k = Y k Y T k yields iteration for Y k...
44 ADI Method for Lyapunov ADI for Lyapunov Newton s Method for AREs For A R n n stable, B R n m (w n), consider Lyapunov equation AX + XA T = BB T. ADI Iteration: [Wachspress 88] (A + p k I )X (j 1)/2 = BB T X k 1 (A T p k I ) (A + p k I )X k T = BB T X (j 1)/2 (A T p k I ) with parameters p k C and p k+1 = p k if p k R. For X 0 = 0 and proper choice of p k : lim k X k = X superlinear. Re-formulation using X k = Y k Y T k yields iteration for Y k...
45 ADI Method for Lyapunov ADI for Lyapunov Newton s Method for AREs For A R n n stable, B R n m (w n), consider Lyapunov equation AX + XA T = BB T. ADI Iteration: [Wachspress 88] (A + p k I )X (j 1)/2 = BB T X k 1 (A T p k I ) (A + p k I )X k T = BB T X (j 1)/2 (A T p k I ) with parameters p k C and p k+1 = p k if p k R. For X 0 = 0 and proper choice of p k : lim k X k = X superlinear. Re-formulation using X k = Y k Y T k yields iteration for Y k...
46 Factored ADI Iteration Lyapunov equation AX + XA T = BB T. Setting X k = Y k Y T k, some algebraic manipulations = ADI for Lyapunov Newton s Method for AREs Algorithm [Penzl 97, Li/White 02, B./Li/Penzl 99/ 07] V 1 p 2Re (p 1)(A + p 1I ) 1 B, Y 1 V 1 FOR j = 2, 3,... q V k Re (pk ) `Vk 1 Re (p k 1 (p ) k + p k 1 )(A + p k I ) 1 V k 1, Y k rrqr `ˆ Y k 1 V k % column compression At convergence, Y kmax Yk T max X, where range (Y kmax ) = range ([ V 1... V kmax ]), Vk = C n m.
47 Factored ADI Iteration Lyapunov equation AX + XA T = BB T. Setting X k = Y k Y T k, some algebraic manipulations = ADI for Lyapunov Newton s Method for AREs Algorithm [Penzl 97, Li/White 02, B./Li/Penzl 99/ 07] V 1 p 2Re (p 1)(A + p 1I ) 1 B, Y 1 V 1 FOR j = 2, 3,... q V k Re (pk ) `Vk 1 Re (p k 1 (p ) k + p k 1 )(A + p k I ) 1 V k 1, Y k rrqr `ˆ Y k 1 V k % column compression At convergence, Y kmax Yk T max X, where range (Y kmax ) = range ([ V 1... V kmax ]), Vk = C n m.
48 Newton s Method for AREs [Kleinman 68, Mehrmann 91, Lancaster/Rodman 95, B./Byers 94/ 98, B. 97, Guo/Laub 99] ADI for Lyapunov Newton s Method for AREs Consider 0 = R(X ) = C T C + A T X + XA XBB T X. Frechét derivative of R(X ) at X : R X : Z (A BBT X ) T Z + Z(A BB T X ). Newton-Kantorovich method: ( ) 1 X j+1 = X j R X R(Xj j ), j = 0, 1, 2,... Newton s method (with line search) for AREs FOR j = 0, 1,... 1 A j A BB T X j =: A BK j. 2 Solve the Lyapunov equation A T j N j + N j A j = R(X j ). 3 X j+1 X j + t j N j. END FOR j
49 Newton s Method for AREs [Kleinman 68, Mehrmann 91, Lancaster/Rodman 95, B./Byers 94/ 98, B. 97, Guo/Laub 99] ADI for Lyapunov Newton s Method for AREs Consider 0 = R(X ) = C T C + A T X + XA XBB T X. Frechét derivative of R(X ) at X : R X : Z (A BBT X ) T Z + Z(A BB T X ). Newton-Kantorovich method: ( ) 1 X j+1 = X j R X R(Xj j ), j = 0, 1, 2,... Newton s method (with line search) for AREs FOR j = 0, 1,... 1 A j A BB T X j =: A BK j. 2 Solve the Lyapunov equation A T j N j + N j A j = R(X j ). 3 X j+1 X j + t j N j. END FOR j
50 Newton s Method for AREs [Kleinman 68, Mehrmann 91, Lancaster/Rodman 95, B./Byers 94/ 98, B. 97, Guo/Laub 99] ADI for Lyapunov Newton s Method for AREs Consider 0 = R(X ) = C T C + A T X + XA XBB T X. Frechét derivative of R(X ) at X : R X : Z (A BBT X ) T Z + Z(A BB T X ). Newton-Kantorovich method: ( ) 1 X j+1 = X j R X R(Xj j ), j = 0, 1, 2,... Newton s method (with line search) for AREs FOR j = 0, 1,... 1 A j A BB T X j =: A BK j. 2 Solve the Lyapunov equation A T j N j + N j A j = R(X j ). 3 X j+1 X j + t j N j. END FOR j
51 Newton s Method for AREs [Kleinman 68, Mehrmann 91, Lancaster/Rodman 95, B./Byers 94/ 98, B. 97, Guo/Laub 99] ADI for Lyapunov Newton s Method for AREs Consider 0 = R(X ) = C T C + A T X + XA XBB T X. Frechét derivative of R(X ) at X : R X : Z (A BBT X ) T Z + Z(A BB T X ). Newton-Kantorovich method: ( ) 1 X j+1 = X j R X R(Xj j ), j = 0, 1, 2,... Newton s method (with line search) for AREs FOR j = 0, 1,... 1 A j A BB T X j =: A BK j. 2 Solve the Lyapunov equation A T j N j + N j A j = R(X j ). 3 X j+1 X j + t j N j. END FOR j
52 Newton s Method for AREs Properties and Implementation ADI for Lyapunov Newton s Method for AREs Convergence for K 0 stabilizing: A j = A BK j = A BB T X j is stable j 0. lim j R(X j ) F = 0 (monotonically). lim j X j = X 0 (locally quadratic). Need large-scale Lyapunov solver; here, ADI iteration: linear systems with dense, but sparse+low rank coefficient matrix A j : A j = A B K j = sparse m m n = efficient inversion using Sherman-Morrison-Woodbury formula: (A BK j ) 1 = (I n + A 1 B(I m K j A 1 B) 1 K j )A 1. BUT: X = X T R n n = n(n + 1)/2 unknowns!
53 Newton s Method for AREs Properties and Implementation ADI for Lyapunov Newton s Method for AREs Convergence for K 0 stabilizing: A j = A BK j = A BB T X j is stable j 0. lim j R(X j ) F = 0 (monotonically). lim j X j = X 0 (locally quadratic). Need large-scale Lyapunov solver; here, ADI iteration: linear systems with dense, but sparse+low rank coefficient matrix A j : A j = A B K j = sparse m m n = efficient inversion using Sherman-Morrison-Woodbury formula: (A BK j ) 1 = (I n + A 1 B(I m K j A 1 B) 1 K j )A 1. BUT: X = X T R n n = n(n + 1)/2 unknowns!
54 Newton s Method for AREs Properties and Implementation ADI for Lyapunov Newton s Method for AREs Convergence for K 0 stabilizing: A j = A BK j = A BB T X j is stable j 0. lim j R(X j ) F = 0 (monotonically). lim j X j = X 0 (locally quadratic). Need large-scale Lyapunov solver; here, ADI iteration: linear systems with dense, but sparse+low rank coefficient matrix A j : A j = A B K j = sparse m m n = efficient inversion using Sherman-Morrison-Woodbury formula: (A BK j ) 1 = (I n + A 1 B(I m K j A 1 B) 1 K j )A 1. BUT: X = X T R n n = n(n + 1)/2 unknowns!
55 Newton s Method for AREs Properties and Implementation ADI for Lyapunov Newton s Method for AREs Convergence for K 0 stabilizing: A j = A BK j = A BB T X j is stable j 0. lim j R(X j ) F = 0 (monotonically). lim j X j = X 0 (locally quadratic). Need large-scale Lyapunov solver; here, ADI iteration: linear systems with dense, but sparse+low rank coefficient matrix A j : A j = A B K j = sparse m m n = efficient inversion using Sherman-Morrison-Woodbury formula: (A BK j ) 1 = (I n + A 1 B(I m K j A 1 B) 1 K j )A 1. BUT: X = X T R n n = n(n + 1)/2 unknowns!
56 Low-Rank Newton-ADI for AREs ADI for Lyapunov Newton s Method for AREs Re-write Newton s method for AREs A T j (X j + N j ) }{{} A T j A T j N j + N j A j = R(X j ) + (X j + N j ) }{{} =X j+1 A j = C T C X j BB T X j }{{} =X j+1 =: W j W T j Set X j = Z j Zj T for rank (Z j ) n = ( Zj+1 Zj+1 T ) ( + Zj+1 Zj+1 T ) Aj = W j Wj T Factored Newton Iteration [B./Li/Penzl 99/ 07] Solve Lyapunov equations for Z j+1 directly by factored ADI iteration and use sparse + low-rank structure of A j.
57 Low-Rank Newton-ADI for AREs ADI for Lyapunov Newton s Method for AREs Re-write Newton s method for AREs A T j (X j + N j ) }{{} A T j A T j N j + N j A j = R(X j ) + (X j + N j ) }{{} =X j+1 A j = C T C X j BB T X j }{{} =X j+1 =: W j W T j Set X j = Z j Zj T for rank (Z j ) n = ( Zj+1 Zj+1 T ) ( + Zj+1 Zj+1 T ) Aj = W j Wj T Factored Newton Iteration [B./Li/Penzl 99/ 07] Solve Lyapunov equations for Z j+1 directly by factored ADI iteration and use sparse + low-rank structure of A j.
58 of Complexity < O(n 3 ) Parallelization: Efficient parallel algorithms based on matrix sign function. Complexity O(n 3 /q) on q-processor machine. library PLiCMR. (B./Quintana-Ortí/Quintana-Ortí since 1999) Formatted Arithmetic: For special problems from PDE control use implementation based on hierarchical matrices and matrix sign function method (Baur/B.), complexity O(n log 2 (n)r 2 ).
59 of Complexity < O(n 3 ) Parallelization: Efficient parallel algorithms based on matrix sign function. Complexity O(n 3 /q) on q-processor machine. library PLiCMR. (B./Quintana-Ortí/Quintana-Ortí since 1999) Formatted Arithmetic: For special problems from PDE control use implementation based on hierarchical matrices and matrix sign function method (Baur/B.), complexity O(n log 2 (n)r 2 ).
60 of Complexity < O(n 3 ) Sparse : Sparse implementation using sparse Lyapunov solver (ADI+MUMPS/SuperLU). Complexity O(n(k 2 + r 2 )). : + Matlab toolbox LyaPack (Penzl 99), + library SpaRed with WebComputing interface. (Badía/B./Quintana-Ortí/Quintana-Ortí since 03)
61 Systems with a Large Number of Terminals Efficient BT implementations are based on assumption n m, p. For on-chip clock distribution networks, power grids, wide buses, this assumption is not justified; here,, e.g., m = p = n 2, n 4. Cure: BT can easily be combined with SVDMOR [Feldmann/Liu 04]: for G(s) = C(sE A) 1 B, let G(s 0) = C(s 0E A) 1 B = ˆ»» Σ 1 V T U 1 U 1 2 Σ 2 V2 T U 1Σ 1V T 1 (rank-k approximation), so that G(s 0) U 1Σ 1V T 1 2 = σ k+1. Now define B := BV 1, C := U T 1 C, then G(s) U 1 B(sE A) 1 C V1 T. {z } =: G(s) Now apply BT to G(s) Ĝ(s), then G(s) U1Ĝ(s)V T 1.
62 Systems with a Large Number of Terminals Efficient BT implementations are based on assumption n m, p. For on-chip clock distribution networks, power grids, wide buses, this assumption is not justified; here,, e.g., m = p = n 2, n 4. Cure: BT can easily be combined with SVDMOR [Feldmann/Liu 04]: for G(s) = C(sE A) 1 B, let G(s 0) = C(s 0E A) 1 B = ˆ»» Σ 1 V T U 1 U 1 2 Σ 2 V2 T U 1Σ 1V T 1 (rank-k approximation), so that G(s 0) U 1Σ 1V T 1 2 = σ k+1. Now define B := BV 1, C := U T 1 C, then G(s) U 1 B(sE A) 1 C V1 T. {z } =: G(s) Now apply BT to G(s) Ĝ(s), then G(s) U1Ĝ(s)V T 1.
63 Systems with a Large Number of Terminals Efficient BT implementations are based on assumption n m, p. For on-chip clock distribution networks, power grids, wide buses, this assumption is not justified; here,, e.g., m = p = n 2, n 4. Cure: BT can easily be combined with SVDMOR [Feldmann/Liu 04]: for G(s) = C(sE A) 1 B, let G(s 0) = C(s 0E A) 1 B = ˆ»» Σ 1 V T U 1 U 1 2 Σ 2 V2 T U 1Σ 1V T 1 (rank-k approximation), so that G(s 0) U 1Σ 1V T 1 2 = σ k+1. Now define B := BV 1, C := U T 1 C, then G(s) U 1 B(sE A) 1 C V1 T. {z } =: G(s) Now apply BT to G(s) Ĝ(s), then G(s) U1Ĝ(s)V T 1.
64 Systems with a Large Number of Terminals Efficient BT implementations are based on assumption n m, p. For on-chip clock distribution networks, power grids, wide buses, this assumption is not justified; here,, e.g., m = p = n 2, n 4. Cure: BT can easily be combined with SVDMOR [Feldmann/Liu 04]: for G(s) = C(sE A) 1 B, let G(s 0) = C(s 0E A) 1 B = ˆ»» Σ 1 V T U 1 U 1 2 Σ 2 V2 T U 1Σ 1V T 1 (rank-k approximation), so that G(s 0) U 1Σ 1V T 1 2 = σ k+1. Now define B := BV1, C := U T 1 C, then G(s) U 1 B(sE A) 1 C V1 T. {z } =: G(s) Now apply BT to G(s) Ĝ(s), then G(s) U1Ĝ(s)V T 1.
65 Systems with a Large Number of Terminals Efficient BT implementations are based on assumption n m, p. For on-chip clock distribution networks, power grids, wide buses, this assumption is not justified; here,, e.g., m = p = n 2, n 4. Cure: BT can easily be combined with SVDMOR [Feldmann/Liu 04]: for G(s) = C(sE A) 1 B, let G(s 0) = C(s 0E A) 1 B = ˆ»» Σ U 1 U 1 V T 1 2 Σ 2 V2 T U 1Σ 1V T 1 (rank-k approximation), so that G(s 0) U 1Σ 1V T 1 2 = σ k+1. Now define B := BV1, C := U T 1 C, then G(s) U 1 B(sE A) 1 C V1 T. {z } =: G(s) Now apply BT to G(s) Ĝ(s), then G(s) U1Ĝ(s)V T 1. multi-point truncated SVD, error bounds Master s thesis A. Schneider, 2007.
66 Systems with a Large Number of Terminals Efficient BT implementations are based on assumption n m, p. For on-chip clock distribution networks, power grids, wide buses, this assumption is not justified; here,, e.g., m = p = n 2, n 4. Cure: BT can easily be combined with SVDMOR [Feldmann/Liu 04]: for G(s) = C(sE A) 1 B, let G(s 0) = C(s 0E A) 1 B = ˆ»» Σ U 1 U 1 V T 1 2 Σ 2 V2 T U 1Σ 1V T 1 (rank-k approximation), so that G(s 0) U 1Σ 1V T 1 2 = σ k+1. Now define B := BV 1, C := U T 1 C, then G(s) U 1 B(sE A) 1 C V1 T. {z } =: G(s) Now apply BT to G(s) Ĝ(s), then G(s) U 1 Ĝ(s)V T 1. Alternative for medium-size m: superposition of reduced-order SIMO models using Padé-type approximation [Feng/B./Rudnyi 07].
67 Optimal Control: Cooling of Steel Profiles Mathematical model: boundary control for linearized 2D heat equation c ρ t x = λ x, ξ Ω λ n x = κ(u k x), ξ Γ k, 1 k 7, x = 0, ξ Γ7. n = m = 7, p = 6. FEM Discretization, different models for initial mesh (n = 371), 1, 2, 3, 4 steps of mesh refinement n = 1357, 5177, 20209, Source: Physical model: courtesy of Mannesmann/Demag. Math. model: Tröltzsch/Unger 1999/2001, Penzl 1999, Saak 2003.
68 Optimal Control: Cooling of Steel Profiles n = 1357, Absolute Error BT model computed with sign function method, MT w/o static condensation, same order as BT model.
69 Optimal Control: Cooling of Steel Profiles n = 1357, Absolute Error n = 79841, Absolute error BT model computed with sign function method, MT w/o static condensation, same order as BT model. BT model computed using SpaRed, computation time: 8 min.
70 MEMS: Microthruster Co-integration of solid fuel with silicon micromachined system. Goal: Ignition of solid fuel cells by electric impulse. Application: nano satellites. Thermo-dynamical model, ignition via heating an electric resistance by applying voltage source. Design problem: reach ignition temperature of fuel cell w/o firing neighboring cells. Spatial FEM discretization of thermo-dynamical model linear system, m = 1, p = 7. Source: The Oberwolfach Benchmark Collection Courtesy of C. Rossi, LAAS-CNRS/EU project Micropyros.
71 MEMS: Microthruster axial-symmetric 2D model FEM discretization using linear (quadratic) elements n = 4, 257 (11, 445) m = 1, p = 7. Reduced model computed using SpaRed, modal truncation using ARPACK, and Z. Bai s PVL implementation.
72 MEMS: Microthruster axial-symmetric 2D model FEM discretization using linear (quadratic) elements n = 4, 257 (11, 445) m = 1, p = 7. Reduced model computed using SpaRed, modal truncation using ARPACK, and Z. Bai s PVL implementation. Relative error n = 4, 257
73 MEMS: Microthruster axial-symmetric 2D model FEM discretization using linear (quadratic) elements n = 4, 257 (11, 445) m = 1, p = 7. Reduced model computed using SpaRed, modal truncation using ARPACK, and Z. Bai s PVL implementation. Relative error n = 4, 257 Relative error n = 11, 445
74 MEMS: Microthruster axial-symmetric 2D model FEM discretization using linear (quadratic) elements n = 4, 257 (11, 445) m = 1, p = 7. Reduced model computed using SpaRed, modal truncation using ARPACK, and Z. Bai s PVL implementation. Frequency Response BT/PVL
75 MEMS: Microthruster axial-symmetric 2D model FEM discretization using linear (quadratic) elements n = 4, 257 (11, 445) m = 1, p = 7. Reduced model computed using SpaRed, modal truncation using ARPACK, and Z. Bai s PVL implementation. Frequency Response BT/PVL Frequency Response BT/MT
76 MEMS: Microgyroscope (Butterfly Gyro) Vibrating micro-mechanical gyroscope for inertial navigation. Rotational position sensor. By applying AC voltage to electrodes, wings are forced to vibrate in anti-phase in wafer plane. Coriolis forces induce motion of wings out of wafer plane yielding sensor data. Source: The Oberwolfach Benchmark Collection Courtesy of D. Billger (Imego Institute, Göteborg), Saab Bofors Dynamics AB.
77 MEMS: Butterfly Gyro FEM discretization of structure dynamical model using quadratic tetrahedral elements (ANSYS-SOLID187) n = 34, 722, m = 1, p = 12. Reduced model computed using SpaRed, r = 30.
78 MEMS: Butterfly Gyro FEM discretization of structure dynamical model using quadratic tetrahedral elements (ANSYS-SOLID187) n = 34, 722, m = 1, p = 12. Reduced model computed using SpaRed, r = 30. Frequency Response Analysis
79 MEMS: Butterfly Gyro FEM discretization of structure dynamical model using quadratic tetrahedral elements (ANSYS-SOLID187) n = 34, 722, m = 1, p = 12. Reduced model computed using SpaRed, r = 30. Frequency Response Analysis Hankel Singular Values
80 Robust Control: Heat Flow in Copper Boundary control problem for 2D heat flow in copper on rectangular domain; control acts on two sides via Robins BC. FDM n = 4496, m = 2; 4 sensor locations p = 4. Numerical ranks of BT Gramians are 68 and 124, respectively, for LQG BT both have rank 210. Computed reduced-order model: r = 10. Source: COMPl eib v1.1,
81 Robust Control: Heat Flow in Copper Boundary control problem for 2D heat flow in copper on rectangular domain; control acts on two sides via Robins BC. FDM n = 4496, m = 2; 4 sensor locations p = 4. Numerical ranks of BT Gramians are 68 and 124, respectively, for LQG BT both have rank 210. Computed reduced-order model: r = 10. Source: COMPl eib v1.1,
82 Robust Control: Heat Flow in Copper Boundary control problem for 2D heat flow in copper on rectangular domain; control acts on two sides via Robins BC. FDM n = 4496, m = 2; 4 sensor locations p = 4. Numerical ranks of BT Gramians are 68 and 124, respectively, for LQG BT both have rank 210. Computed reduced-order model: r = 10. Source: COMPl eib v1.1,
83 Model Reduction SLICOT Model and Controller Reduction Toolbox Stand-alone Matlab Toolbox based on the Subroutine Library in Systems and Control Theory SLICOT which provides > 400 Fortran 77 routines for systems and control related computations; other Matlab Toolboxes: SLICOT Basic Systems and Control Toolbox, SLICOT System Identification Toolbox. Much enhanced functionality compared to Matlab s Control Toolboxes, in particular coprime factorization methods, frequency-weighted BT methods, controller reduction. Maintained by NICONET e.v. Distributed by Synoptio GmbH, Berlin. For more information, visit and
84 Model Reduction LyaPack Matlab toolbox for computations involving large-scale Lyapunov equations with sparse coefficient matrix A and low-rank constant term; contains BT implementation and dominant subspace approximation method. Available as additional software (no registration necessary) from New version (algorithmic improvements, easier to use) in progress. PLiCMR, SpaRed Parallel implementations of BT and some balancing-related methods for dense and sparse linear systems. MorLAB Model-reduction LABoratory, contains Matlab functions for balancing-related methods, available from
85 1 G. Obinata and B.D.O. Anderson. Model Reduction for Control System Design. Springer-Verlag, London, UK, Z. Bai. Krylov subspace techniques for reduced-order modeling of large-scale dynamical systems. Appl. Numer. Math, 43(1 2):9 44, R. Freund. Model reduction methods based on Krylov subspaces. Acta Numerica, 12: , P. Benner, E.S. Quintana-Ortí, and G. Quintana-Ortí. State-space truncation methods for parallel model reduction of large-scale systems. Parallel Comput., 29: , P. Benner. Numerical linear algebra for model reduction in control and simulation. GAMM Mitt., 29(2): , P. Benner, V. Mehrmann, and D. Sorensen (editors). Dimension Reduction of Large-Scale Systems. Lecture Notes in Computational Science and Engineering, Vol. 45, Springer-Verlag, Berlin/Heidelberg, Germany, A.C. Antoulas. Lectures on the Approximation of Large-Scale Dynamical Systems. SIAM Publications, Philadelphia, PA, P. Benner, R. Freund, D. Sorensen, and A. Varga (editors). Special issue on Order Reduction of Large-Scale Systems. Linear Algebra Appl., Vol. 415, June 2006.
BALANCING-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 informationSYSTEM-THEORETIC METHODS FOR MODEL REDUCTION OF LARGE-SCALE SYSTEMS: SIMULATION, CONTROL, AND INVERSE PROBLEMS. Peter Benner
SYSTEM-THEORETIC METHODS FOR MODEL REDUCTION OF LARGE-SCALE SYSTEMS: SIMULATION, CONTROL, AND INVERSE PROBLEMS Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität
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 informationModel Reduction for Dynamical Systems
Otto-von-Guericke Universität Magdeburg Faculty of Mathematics Summer term 2015 Model Reduction for Dynamical Systems Lecture 8 Peter Benner Lihong Feng Max Planck Institute for Dynamics of Complex Technical
More informationPassivity Preserving Model Reduction for Large-Scale Systems. Peter Benner.
Passivity Preserving Model Reduction for Large-Scale Systems Peter Benner Mathematik in Industrie und Technik Fakultät für Mathematik Sonderforschungsbereich 393 S N benner@mathematik.tu-chemnitz.de SIMULATION
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 informationModel Reduction for Linear Dynamical Systems
Summer School on Numerical Linear Algebra for Dynamical and High-Dimensional Problems Trogir, October 10 15, 2011 Model Reduction for Linear Dynamical Systems Peter Benner Max Planck Institute for Dynamics
More informationMODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS
MODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS Ulrike Baur joint work with Peter Benner Mathematics in Industry and Technology Faculty of Mathematics Chemnitz University of Technology
More informationEfficient Implementation of Large Scale Lyapunov and Riccati Equation Solvers
Efficient Implementation of Large Scale Lyapunov and Riccati Equation Solvers Jens Saak joint work with Peter Benner (MiIT) Professur Mathematik in Industrie und Technik (MiIT) Fakultät für Mathematik
More informationBalancing-Related Model Reduction for Large-Scale Systems
Balancing-Related Model Reduction for Large-Scale Systems Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz D-09107 Chemnitz benner@mathematik.tu-chemnitz.de
More informationS N. hochdimensionaler Lyapunov- und Sylvestergleichungen. Peter Benner. Mathematik in Industrie und Technik Fakultät für Mathematik TU Chemnitz
Ansätze zur numerischen Lösung hochdimensionaler Lyapunov- und Sylvestergleichungen Peter Benner Mathematik in Industrie und Technik Fakultät für Mathematik TU Chemnitz S N SIMULATION www.tu-chemnitz.de/~benner
More informationModel 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 informationParallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation
Parallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation Peter Benner 1, Enrique S. Quintana-Ortí 2, Gregorio Quintana-Ortí 2 1 Fakultät für Mathematik Technische Universität
More informationBalanced Truncation Model Reduction of Large and Sparse Generalized Linear Systems
Balanced Truncation Model Reduction of Large and Sparse Generalized Linear Systems Jos M. Badía 1, Peter Benner 2, Rafael Mayo 1, Enrique S. Quintana-Ortí 1, Gregorio Quintana-Ortí 1, A. Remón 1 1 Depto.
More informationMoving Frontiers in Model Reduction Using Numerical Linear Algebra
Using Numerical Linear Algebra Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory Group Magdeburg, Germany Technische Universität Chemnitz
More informationModel reduction of coupled systems
Model reduction of coupled systems Tatjana Stykel Technische Universität Berlin ( joint work with Timo Reis, TU Kaiserslautern ) Model order reduction, coupled problems and optimization Lorentz Center,
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 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 informationParametrische Modellreduktion mit dünnen Gittern
Parametrische Modellreduktion mit dünnen Gittern (Parametric model reduction with sparse grids) Ulrike Baur Peter Benner Mathematik in Industrie und Technik, Fakultät für Mathematik Technische Universität
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 informationBalancing-Related Model Reduction for Parabolic Control Systems
Balancing-Related Model Reduction for Parabolic Control Systems Peter Benner Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany (e-mail: benner@mpi-magdeburg.mpg.de). Abstract:
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 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 informationFast Frequency Response Analysis using Model Order Reduction
Fast Frequency Response Analysis using Model Order Reduction Peter Benner EU MORNET Exploratory Workshop Applications of Model Order Reduction Methods in Industrial Research and Development Luxembourg,
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationModel order reduction of electrical circuits with nonlinear elements
Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a
More informationAccelerating Model Reduction of Large Linear Systems with Graphics Processors
Accelerating Model Reduction of Large Linear Systems with Graphics Processors P. Benner 1, P. Ezzatti 2, D. Kressner 3, E.S. Quintana-Ortí 4, Alfredo Remón 4 1 Max-Plank-Institute for Dynamics of Complex
More informationFactorized Solution of Sylvester Equations with Applications in Control
Factorized Solution of Sylvester Equations with Applications in Control Peter Benner Abstract Sylvester equations play a central role in many areas of applied mathematics and in particular in systems and
More informationNumerical Solution of Descriptor Riccati Equations for Linearized Navier-Stokes Control Problems
MAX PLANCK INSTITUTE Control and Optimization with Differential-Algebraic Constraints Banff, October 24 29, 2010 Numerical Solution of Descriptor Riccati Equations for Linearized Navier-Stokes Control
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
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 informationModel reduction of nonlinear circuit equations
Model reduction of nonlinear circuit equations Tatjana Stykel Technische Universität Berlin Joint work with T. Reis and A. Steinbrecher BIRS Workshop, Banff, Canada, October 25-29, 2010 T. Stykel. Model
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 informationParallel Solution of Large-Scale and Sparse Generalized algebraic Riccati Equations
Parallel Solution of Large-Scale and Sparse Generalized algebraic Riccati Equations José M. Badía 1, Peter Benner 2, Rafael Mayo 1, and Enrique S. Quintana-Ortí 1 1 Depto. de Ingeniería y Ciencia de Computadores,
More informationGramian-Based Model Reduction for Data-Sparse Systems CSC/ Chemnitz Scientific Computing Preprints
Ulrike Baur Peter Benner Gramian-Based Model Reduction for Data-Sparse Systems CSC/07-01 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific Computing Preprints ISSN 1864-0087 (1995
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 informationModel Order Reduction of Electrical Circuits with Nonlinear Elements
Model Order Reduction of Electrical Circuits with Nonlinear Elements Tatjana Stykel and Technische Universität Berlin July, 21 Model Order Reduction of Electrical Circuits with Nonlinear Elements Contents:
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
More informationFluid flow dynamical model approximation and control
Fluid flow dynamical model approximation and control... a case-study on an open cavity flow C. Poussot-Vassal & D. Sipp Journée conjointe GT Contrôle de Décollement & GT MOSAR Frequency response of an
More information6 OUTPUT FEEDBACK DESIGN
6 OUTPUT FEEDBACK DESIGN When the whole sate vector is not available for feedback, i.e, we can measure only y = Cx. 6.1 Review of observer design Recall from the first class in linear systems that a simple
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 informationIn this article we discuss sparse matrix algorithms and parallel algorithms, Solving Large-Scale Control Problems
F E A T U R E Solving Large-Scale Control Problems EYEWIRE Sparsity and parallel algorithms: two approaches to beat the curse of dimensionality. By Peter Benner In this article we discuss sparse matrix
More informationModel order reduction of large-scale dynamical systems with Jacobi-Davidson style eigensolvers
MAX PLANCK INSTITUTE International Conference on Communications, Computing and Control Applications March 3-5, 2011, Hammamet, Tunisia. Model order reduction of large-scale dynamical systems with Jacobi-Davidson
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 informationComparison of methods for parametric model order reduction of instationary problems
Max Planck Institute Magdeburg Preprints Ulrike Baur Peter Benner Bernard Haasdonk Christian Himpe Immanuel Maier Mario Ohlberger Comparison of methods for parametric model order reduction of instationary
More informationProblem set 5 solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION Problem set 5 solutions Problem 5. For each of the stetements below, state
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
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 informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
More informationKrylov Techniques for Model Reduction of Second-Order Systems
Krylov Techniques for Model Reduction of Second-Order Systems A Vandendorpe and P Van Dooren February 4, 2004 Abstract The purpose of this paper is to present a Krylov technique for model reduction of
More informationIntro. Computer Control Systems: F9
Intro. Computer Control Systems: F9 State-feedback control and observers Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 21 dave.zachariah@it.uu.se F8: Quiz! 2 / 21 dave.zachariah@it.uu.se
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationPrinciple Component Analysis and Model Reduction for Dynamical Systems
Principle Component Analysis and Model Reduction for Dynamical Systems D.C. Sorensen Virginia Tech 12 Nov 2004 Protein Substate Modeling and Identification PCA Dimension Reduction Using the SVD Tod Romo
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 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 informationSolving Large Scale Lyapunov Equations
Solving Large Scale Lyapunov Equations D.C. Sorensen Thanks: Mark Embree John Sabino Kai Sun NLA Seminar (CAAM 651) IWASEP VI, Penn State U., 22 May 2006 The Lyapunov Equation AP + PA T + BB T = 0 where
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 informationBDF Methods for Large-Scale Differential Riccati Equations
BDF Methods for Large-Scale Differential Riccati Equations Peter Benner Hermann Mena Abstract We consider the numerical solution of differential Riccati equations. We review the existing methods and investigate
More informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 WeC10.4 An iterative SVD-Krylov based method for model reduction
More informationRiccati Equations and Inequalities in Robust Control
Riccati Equations and Inequalities in Robust Control Lianhao Yin Gabriel Ingesson Martin Karlsson Optimal Control LP4 2014 June 10, 2014 Lianhao Yin Gabriel Ingesson Martin Karlsson (LTH) H control problem
More informationStructure preserving Krylov-subspace methods for Lyapunov equations
Structure preserving Krylov-subspace methods for Lyapunov equations Matthias Bollhöfer, André Eppler Institute Computational Mathematics TU Braunschweig MoRePas Workshop, Münster September 17, 2009 System
More informationOn Linear-Quadratic Control Theory of Implicit Difference Equations
On Linear-Quadratic Control Theory of Implicit Difference Equations Daniel Bankmann Technische Universität Berlin 10. Elgersburg Workshop 2016 February 9, 2016 D. Bankmann (TU Berlin) Control of IDEs February
More informationModel Order Reduction for Parameter-Varying Systems
Model Order Reduction for Parameter-Varying Systems Xingang Cao Promotors: Wil Schilders Siep Weiland Supervisor: Joseph Maubach Eindhoven University of Technology CASA Day November 1, 216 Xingang Cao
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationIterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations
Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University
More informationIntroduction to Model Order Reduction
Introduction to Model Order Reduction Lecture 1: Introduction and overview Henrik Sandberg Kin Cheong Sou Automatic Control Lab, KTH ACCESS Specialized Course Graduate level Ht 2010, period 1 1 Overview
More informationSystem Identification by Nuclear Norm Minimization
Dept. of Information Engineering University of Pisa (Italy) System Identification by Nuclear Norm Minimization eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain Systems
More informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 1964 1986 www.elsevier.com/locate/laa An iterative SVD-Krylov based method for model reduction of large-scale dynamical
More informationReduced-order Model Based on H -Balancing for Infinite-Dimensional Systems
Applied Mathematical Sciences, Vol. 7, 2013, no. 9, 405-418 Reduced-order Model Based on H -Balancing for Infinite-Dimensional Systems Fatmawati Department of Mathematics Faculty of Science and Technology,
More information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
More informationModel reduction of interconnected systems
Model reduction of interconnected systems A Vandendorpe and P Van Dooren 1 Introduction Large scale linear systems are often composed of subsystems that interconnect to each other Instead of reducing the
More informationOrder Reduction of a Distributed Parameter PEM Fuel Cell Model
Order Reduction of a Distributed Parameter PEM Fuel Cell Model María Sarmiento Carnevali 1 Carles Batlle Arnau 2 Maria Serra Prat 2,3 Immaculada Massana Hugas 2 (1) Intelligent Energy Limited, (2) Universitat
More informationModel reduction of large-scale dynamical systems
Model reduction of large-scale dynamical systems Lecture III: Krylov approximation and rational interpolation Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/
More informationModel Reduction of Inhomogeneous Initial Conditions
Model Reduction of Inhomogeneous Initial Conditions Caleb Magruder Abstract Our goal is to develop model reduction processes for linear dynamical systems with non-zero initial conditions. Standard model
More informationDECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION
Journal of ELECTRICAL ENGINEERING, VOL. 58, NO. 6, 2007, 307 312 DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Szabolcs Dorák Danica Rosinová Decentralized control design approach based on partial
More informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
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 informationParametric Model Order Reduction for Linear Control Systems
Parametric Model Order Reduction for Linear Control Systems Peter Benner HRZZ Project Control of Dynamical Systems (ConDys) second project meeting Zagreb, 2 3 November 2017 Outline 1. Introduction 2. PMOR
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 informationModel reduction of large-scale systems by least squares
Model reduction of large-scale systems by least squares Serkan Gugercin Department of Mathematics, Virginia Tech, Blacksburg, VA, USA gugercin@mathvtedu Athanasios C Antoulas Department of Electrical and
More informationTowards high performance IRKA on hybrid CPU-GPU systems
Towards high performance IRKA on hybrid CPU-GPU systems Jens Saak in collaboration with Georg Pauer (OVGU/MPI Magdeburg) Kapil Ahuja, Ruchir Garg (IIT Indore) Hartwig Anzt, Jack Dongarra (ICL Uni Tennessee
More informationStructure preserving model reduction of network systems
Structure preserving model reduction of network systems Jacquelien Scherpen Jan C. Willems Center for Systems and Control & Engineering and Technology institute Groningen Reducing dimensions in Big Data:
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationAN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL MOTIVATING EXAMPLE INVERTED PENDULUM
Controls Lab AN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL Eduardo Gildin (UT ICES and Rice Univ.) with Thanos Antoulas (Rice ECE) Danny Sorensen (Rice
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 informationOn the numerical solution of large-scale linear matrix equations. V. Simoncini
On the numerical solution of large-scale linear matrix equations V. Simoncini Dipartimento di Matematica, Università di Bologna (Italy) valeria.simoncini@unibo.it 1 Some matrix equations Sylvester matrix
More informationSyreNe System Reduction for Nanoscale IC Design
System Reduction for Nanoscale Max Planck Institute for Dynamics of Complex Technical Systeme Computational Methods in Systems and Control Theory Group Magdeburg Technische Universität Chemnitz Fakultät
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 informationMPC/LQG for infinite-dimensional systems using time-invariant linearizations
MPC/LQG for infinite-dimensional systems using time-invariant linearizations Peter Benner 1 and Sabine Hein 2 1 Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstr. 1, 39106 Magdeburg,
More informationNUMERICAL ANALYSIS AND SYSTEMS THEORY
Int. J. Appl. Math. Comput. Sci., 2001, Vol.11, No.5, 1025 1033 NUMERICAL ANALYSIS AND SYSTEMS THEORY Stephen L. CAMPBELL The area of numerical analysis interacts with the area of control and systems theory
More informationKrylov Subspace Methods for Nonlinear Model Reduction
MAX PLANCK INSTITUT Conference in honour of Nancy Nichols 70th birthday Reading, 2 3 July 2012 Krylov Subspace Methods for Nonlinear Model Reduction Peter Benner and Tobias Breiten Max Planck Institute
More informationEmpirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems
Empirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems Antoni Ras Departament de Matemàtica Aplicada 4 Universitat Politècnica de Catalunya Lecture goals To review the basic
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationDownloaded 05/27/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. MATRIX ANAL. APPL. Vol. 31, No. 2, pp. 272 288 c 2009 Society for Industrial and Applied Mathematics INEXACT KLEINMAN NEWTON METHOD FOR RICCATI EQUATIONS F. FEITZINGER, T. HYLLA, AND E. W. SACHS
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
More information