Rational Krylov Methods for Model Reduction of Large-scale Dynamical Systems
|
|
- Toby Skinner
- 5 years ago
- Views:
Transcription
1 Rational Krylov Methods for Model Reduction of Large-scale Dynamical Systems Serkan Güǧercin Department of Mathematics, Virginia Tech, USA Aerospace Computational Design Laboratory, Dept. of Aeronautics & Astronautics, MIT, October 21, 2005, Cambridge, MA 1
2 Outline 1. Introduction and Problem Statement 2. Motivating Examples 3. Rational Krylov-Interpolation Framework 4. Current Trends in Rational Krylov Projection 5. Conclusions 2
3 Introduction Consider an n th order single-input/single-output system G(s): ẋ(t) = Ax(t) + bu(t) G(s) : G(s) = c(si n A) 1 b y(t) = cx(t) = n(s) d(s) u(t) R: input, x(t) R n : state, y(t) R: output A R n n, b,c T R n. Will assume R(λ i (A)) < 0 Need for improved accuracy = Include more details in the modeling stage In many applications, n is quite large, n O(10 6, 10 7 ), Untenable demands on computational resources = 3
4 Model Reduction Problem: Find ẋ r (t) = A r x r (t) + b r u(t) y r (t) = c r x r (t) G r (s) = c r (si r A r ) 1 b r where A r R r r, b r,c T r R r, with r n such that 1. y y r is small. 2. The procedure is computationally efficient. A B A r B r C D C r D r G r (s): used for simulation or designing a reduced-order controller 4
5 Model reduction through projection: a unifying framework. Construct Π = VZ T, where V,Z R n r with Z T V = I r : ẋ r = Z T AV }{{} :=A r x r (t) + Z T b }{{} u(t), y r (t) = }{{} cv x r (t) :=b r :=c r What is the approximation error e(t) := y(t) y r (t)? G(s): Associate a convolution operator S: S : u(t) y(t) = (Su)(t) = (g u)(t) = t g(t τ)u(τ)dτ. g(t) = ce At b for t 0: Impulse response. Transfer function: G(s) = (Lg)(s) = c(si A) 1 b. 5
6 The H Norm : 2-2 induced norm of S: G(s) H = sup u 0 y 2 u 2 = sup u 0 Su 2 = sup G(jw) u 2 2 w R G G r = Worst output error y(t) y r (t) 2 u(t) 2 = 1. The H 2 Norm : L 2 norm of g(t) in time domain: G(s) 2 H 2 = 0 trace[g T (t)g(t)]dt = 1 2π + trace[g (jw)g(jw)]dw G(s) G(jw) (db) G(s) freq (rad/sec) 6
7 Motivating Example: Simulation A Tunable Optical Filter: (Data: D. Hohlfeld, T. Bechtold, and H. Zappe) An optical filter, tunable by thermal means. Silicon-based fabrication. The thin-film filter: membrane to improve thermal isolation Wavelength tuning by thermal modulation of resonator optical thickness The device features low power consumption, high tuning speed and excellent optical performance. 7
8 Modeling: A simplified thermal model to analyze/simulate important thermal issues: 2D model and 3D model Meshed and discretized in ANSYS 6.1 by the finite element methods The Dirichlet boundary conditions at the bottom of the chip. A constant load vector corresponding to the constant input power of of 1 mw for 2D model and 10 mw for 3D model The output nodes located in the membrane Eẋ(t) = Ax(t) + bu(t), y(t) = cx(t) 2D: n = 1668, nnz(a) = 6209, nnz(e) = D: n = , nnz(a) = , nnz(e) =
9 Motivating Example: Control International Space Station (ISS) Modules: ISS: a complex structure of many modules. More than 40 Shuttle flights to complete the assembly. Each module modeled by 1000 state variables For ISS, controllers will be needed for many reasons: 1. to reduce the oscillatory motions; 2. to fix the orientation of the space station with respect to some desired direction Controllers to be implemented on-board Controllers of low complexity needed due to: hardware, radiation, throughput and testing issues 9
10 Flex Structure Variation During Assembly Stage 1R - Dec 99 Stage 4A - Mar 00 Stage 13A - May 02 Figure 1: ISS: Evolution of the frequency response as more modules are added (Data: Draper Labs) 10
11 SVD based approximation methods Singular Value Decomposition (SVD) provides low-rank approximation of matrices. What are the singular values of G(s)? Two Lyapunov Equations: AP + PA T + bb T = 0, A T Q + QA + c T c = 0, 0 < P : Reachability Gramian 0 < Q : Observability Gramian Hankel Singular Values: σ i (Σ) := λ i (PQ) for i = 1,...,n. States difficult to reach correspond to small eigenvalues of P. States difficult to observe correspond to small eigenvalues of Q. 11
12 Balanced Truncation Method States difficult to reach States difficult to observe? Make P = Q by a basis change: x = Tx, dett 0 Ā = TAT 1, b = Tb, c = ct 1. P = UU T, U T QU = RΣ 2 R T = T b := Σ 1/2 R T U 1 In the balanced basis: P = Q = Σ = diag(σ1,...,σ n ) Projection Framework: P = UU T, Q = LL T, U T L = WΣY T, Σ = diag(σ 1,...,σ n ). Z := LY 1 Σ 1/2 1, V := UW 1 Σ 1/2 1. VZ T : oblique projector G r (s) = A r b r c r 0 = ZT AV Z T b cv 0. 12
13 Properties of the reduced model : 1. G r (s) is asymptotically stable. 2. G(s) G r (s) 2 (σ r σ n ). Drawbacks in large-scale settings : Requires U and L 1. O(n 3 ) arithmetic operations 2. O(n 2 ) memory requirement Remedy: σ i (Σ), λ i (P) and λ i (Q) decay rapidly. = Find Ũ, L R n k s.t. P := ŨŨ P and Q := L L T Q. (Approximate) balancing using Ũ and L. O(n 2 k) arithmetic operations, O(nk) memory requirement. Penzl [2000], Li/White [2000]: LR-Smith, G/Sorensen/Antoulas [2002]: Modified LR-Smith, Willcox/Peraire [2002]: Balanced POD 13
14 Model Reduction via Interpolation Rational Interpolation: Given G(s), find G r (s) so that G r (s) interpolates G(s) and certain number of its derivatives at selected frequencies σ k in the complex plane ( 1) j j! d j G(s) ds j = ( 1)j s=σk j! d j G r (s) ds j, s=σk for k = 1,...,K, and j = 1,...,J ( 1)j j! d j G(s) ds j = c(σ k I A) (j+1) b: s=σk = j th moment of G(s) at σ k =: η σk (j). 14
15 For σ 1 = σ 2 =... = σ K = σ = η (j) σ = ca j 1 b : j th Markov Parameter Moment matching problem: Partial realization Solution: Lanczos and Arnoldi procedures For arbitrary σ η (j) σ = c(σi A) (j+1) b Moment matching problem: Rational Interpolation Solution: Rational Lanczos/Arnoldi procedures For multiple points = rational Krylov method 15
16 Moments are very ill-conditioned to compute Markov Parameters, σ = 10 2 Moments at σ = Ball et al. [1990], Antoulas et al. [1990], Pillage et al. [1993]: Moment matching problem with explicit moment computation. = Total loss of accuracy in large-scale settings Feldman and Freund [1995]: Padé via Lanczos ( Krylov setting) Moment matching without moment computation Iterative implementation 16
17 The Arnoldi Procedure for σ = 1. Given A R n n, b R n, and c T R n, 2. v 1 := b b, w := Av 1; α 1 := v T 1 w f 1 := w v 1 α 1 ; V 1 := (v 1 ); H 1 := (α 1 ) 3. For j = 1, 2,, r 1 β j := f j, v j+1 := f j β j V j+1 := (V j v j+1 ), Ĥ j = H j β j e T j w := Av j+1, h := Vj+1w, T f j+1 = w V j+1 h ) H j+1 := (Ĥj h AV = VH + fe T r, V T V = I r, V T f = 0, V R n r A r = H, b r = b e 1, c r = cv 17
18 G r (s) matches the first r Markov parameters of G(s) without computing them Only matrix-vector multiplications. No SVD, or Schur decomposition required Range(V) = Span ( [ b, Ab,, A r 1 b ] ) Lanczos Procedure: Two-sided Arnoldi. Range(V) = Span ( [ b, Ab,, A r 1 b ] ) Range(Z) = Span ( [ c T, A T c T,, (A T ) r 1 c T ] ) G r (s) matches the first 2r Markov parameters of G(s) without computing them For arbitrary σ, A (σi A) 1, b (σi A) 1 b, c c(σi A T ) 1 18
19 Approximation of a Building Model Los Angeles University Hospital: n = 48 to r = Bal. Red. 40 Σ H(jw) (db) σ = Σ σ= Σ Freq (rad/sec) 19
20 Rational Krylov Method Lanczos/Arnoldi: Moment matching at a single point. Large error around other frequencies Match the moments around various frequencies Better approximation over a broad range Solution by projection: First by Skelton et. al. [1985] Grimme [1997]: Efficient computation by Krylov projection 20
21 Given 2r interpolation points: {σ i } 2r i=1 Set V = Span [ (σ 1 I A) 1 b,, (σ r I A) 1 b ], and Z = Span [ (σ r+1 I A T ) 1 c T,, (σ 2r I A T ) 1 c T ], Make Z T V = I r. A r = Z T AV, b r = Z T b, c r = cv = G(σ i ) = G r (σ i ) for i = 1,...,2r Moment matching without explicit moment computation 21
22 Given r interpolation points: {σ i } r i=1 Set V = Span [ (σ 1 I A) 1 b,, (σ r I A) 1 b ], and Z = Span [ (σ 1 I A T ) 1 c T,, (σ r I A T ) 1 c T ], Make Z T V = I r. A r = Z T AV, b r = Z T b, c r = cv = G(σ i ) = G r (σ i ), and d ds G(s) = d s=σi ds G r(s) Moment matching without explicit moment computation s=σi Efficient computations of V and Z: Dual Rational Arnoldi, Rational Lanczos, Rational Power Krylov.(Grimme [1997]) 22
23 Why to choose model reduction via rational interpolation? van Dooren et al. [2004], Halevi [2005]: Generically, any reduced model G r (s) can be obtained via interpolation. Interpolation points = Zeroes of G(s) G r (s). Example: G(s) = G(s) G r (s) = 1 s 2 + 3s + 2 V = (σ 1 A) 1 b = A 1 b, BUT: and G r (s) = 0.2 s s(s 2) s s s σ 1 = 0, σ 2 = 2 Z = (σ 2 A) T c T = (2I A) T c T What is a good selection of interpolation points? Similar to polynomial approximation of complex functions. 23
24 Advantages of Krylov based methods: Iterative implementation Only matrix-vector multiplications and sparse LU decompositions O(r 2 n) arithmetic operation and O(nr) storage requirements Drawbacks: Stability is not guaranteed No global error bound Selection of σ i is an ad-hoc process. 24
25 Stability: Implicit Restart (Sorensen et al.[1995]) Let A r,b r,c r be obtained via r steps of Arnoldi algorithm. Let A r have one unstable eigenvalue µ, i.e., R(µ) > 0. Define b := (µi n A)b Run r 1 steps of Arnoldi process on the pair (A, b) The new reduced system will be asymptotically stable Matching the moments of Ĝ(s) = c(si A) 1 b Trade-off between guaranteed stability and exact moment matching. 25
26 Least-squares model reduction (G./Antoulas [2003]) Bilinear transformation and reduction in the discrete-time Combines the SVD and Krylov reduction V : Krylov-side, Z : SVD-side Stability and moment matching Fourier model reduction (Willcox/Megretski [2005]) Bilinear transformation and reduction in the discrete-time Stability and moment matching 26
27 An H 2 Error Expression (G./Antoulas [2002]) G(s) 2 H 2 = cpc T = b T Qb. φ i := G(s)(s λ i ) s=λi and φ j := G r (s)(s λ j ) s= λj. G(s) G r (s) 2 H 2 = n φ i (G( λ i ) G r ( λ i )) + i=1 r j=1 ) φ j (G r ( λ j ) G( λ j ). Error due to mismatch at λ i and λ j. Try to kill the first term σ i = λ i (A). Select λ i (A) with large φ i. 27
28 CD Player and Selection of σ i Dynamics of a CD Player with 120 states, m = p = 1. Reduce the order to r = 14 by Rational Krylov G r: Choosing σ i = λ i (A). G 1,G 2,G 3 : Arbitrary σ i. Best 3 among 2000 models. G r G 1 G 2 G 3 H H Table 1: H 2 and H errors for the CD Player Example Does there exist an optimal shift selection? 28
29 Optimal H 2 approximation Problem: Given a stable dynamical system G(s), find a reduced model G r (s) that satisfies G r (s) = arg G(s) Ĝ(s) H2. min deg(ĝ) = r Ĝ : stable First-order conditions: (Meier and Luenberger [1967]) d G( ˆλ i ) = G r ( ˆλ i ), and ds G(s) = d s= ˆλi ds G r(s) s= ˆλi Match the first two moments at the mirror images of the Ritz values. First-order conditions as interpolation. Rational Krylov Framework 29
30 For the H 2 problem, simply set σ i = ˆλ i ˆλ i NOT known a priori = Needs iterative rational steps An Iterative Rational Krylov Algorithm (IRKA): (G..Beattie/Antoulas [2004]) 1. Choose σ i for i = 1,..., r. 2. V = Span [ (σ 1 I A) 1 b,, (σ r I A) 1 b ], 3. Z = Span [ (σ 1 I A T ) 1 c T,, (σ r I A T ) 1 c T ], Z T V = I r. 4. while [relative change in σ j ] > ǫ (a) A r = Z T AV, (b) σ i λ i (A r ) for i = 1,..., r (c) V = Span [ (σ 1 I A) 1 b,, (σ r I A) 1 b ]. (d) Z = Span [ (σ 1 I A T ) 1 c T,, (σ r I A T ) 1 c T ], Z T V = I r. 5. A r = Z T AV, b r = Z T b, c r = cv Upon convergence, first-order conditions satisfied via Krylov projection framework, no Lyapunov solvers 30
31 ISS 12a Module n = Reduce to r = 2 : 2 : 60 Compare with balanced truncation Comparision between IRKA and BT 1 Evolution of IRKA for r=2 Real(σ 1 ) Relative H 2 error 10 1 BT Imag(σ 1 ) IRKA r: reduced order Relative H 2 error Number of iterations 31
32 Example: Optimal Cooling of Steel Profiles ( P. Benner) G(s) = c(se A) 1 b, n = 79, 841 r = 6 via IRKA 10 1 Bode Plots of G(s) and G r (s) (IRKA) G(s) G r (s) 10 2 H(jw) frequency (rad/sec) G(s) G 1 (s) =
33 van Dooren et al.[2005]: Model reduction of interconnected systems via rational Krylov Reduce each subsystem via interpolation Reduced-interconnection interpolates the full-model Bai [2003]: Krylov-reduction of second-order systems: Mẍ(t) + Gẋ(t) + Kx(t) = bu(t), y(t) = cx(t) A = G(s) : 0 I M 1 K M 1 G,B = q(t) = Aq(t) + B u(t) y(t) = Cq(t) 0 M 1 b,c T = ct Apply model reduction in the first-order state-space Transfer back to second-order form Structure is lost 0 T 33
34 find a matrix W R n r such that W T W = I r M r = W T MW, G r = W T GW, K r = W T KW, b r = W T b, and c r = cw M r ẍ r (t) + G r ẋ r (t) + K r x r (t) = b r u(t), y r (t) = c r x r (t) Second-order Krylov subspaces (Second-order Arnoldi Iteration): W spans the required Krylov subspace but constructed directly in the second-order framework Interpolation conditions hold. Antoulas/Sorensen [2003]: Passivity preserving Krylov-based model reduction G(jw) + (G(jw)) > 0. G(s) + G (s) = H(s)H (s) Interpolation points: Zeros of H(s) Passivity is preserved. 34
35 Conclusions and Future Work: Rational Krylov framework for model reduction In the SISO case, any G r (s) can be obtained via interpolation Match moments without moment computation Cost is reduced from O(n 3 ) to O(r 2 n) To guarantee stability Implicit restart Least-squares reduction (G./Antoulas) Fourier reduction (Willcox/Megretski) Optimal selection of interpolation points Interpolation condition for optimal H 2 approximation Iterative Rational Krylov Algorithm On-going and future work - Interconnected systems - Second-order systems - Controller reduction - Optimal approximation - Guaranteed error estimates - Balancing via Krylov 35
36 Controller reduction for large-scale systems Consider an n th order plant G(s) = c(si A) 1 b n th κ order stabilizing controller: K(s) = c K (si A K ) 1 b K + d K LQG, H control designs n κ = n (i) Complex hardware (ii) Degraded accuracy (iii) Degraded computational speed Obtain K r (s) of order r n κ to replace K(s) in the closed loop. 36
37 Controller reduction via frequency weighting Small open loop error K(s) K r (s) not enough. Minimize the weighted error: W o (s)(k(s) K r (s))w i (s). How to obtain the weights W o (s) and W i (s)? If K(s) and K r (s) have the same number of unstable poles and if [K(s) Kr (s)]]g(s)[i + G(s)K(s)] 1 < 1, or [I + G(s)K(s)] 1 G(s)[K(s) K r (s)] < 1, = K r (s) stabilizes G(s). = 37
38 For stability considerations: W i (s) = I and W o (s) = [I + G(s)K(s)] 1 G(s) or W o (s) = I and W i (s) = G(s)[I + G(s)K(s)] 1. To preserve closed-loop performance: W i (s) = [I + G(s)K(s)] 1 and W o (s) = [I + G(s)K(s)] 1 G(s). Solved by frequency-weighted balancing (Anderson and Liu [1989], Schelfhout and De Moor [1996], Varga and Anderson [2002] ). Requires solving two Lyapunov equations of order n + n κ. A i P + PA T i + b i b T i = 0, A T o Q + QA o + c T o c o = 0, A i,b i : K(s)W i (s), A o,c o : W o (s)k(s) Balance P and Q. 38
39 Controller-reduction via Krylov Projection How to modify IRKA for the controller reduction problem? Let W i (s) = I and W o (s) = [I + G(s)K(s)] 1 G(s) A K P + PA T K + b K b T K = 0 }{{} unweighted Lyapunov eq. A T wq + QA w + c T wc w = 0. }{{} weighted Lyapunov eq. Z = K(A T,C T, σ i ), and V = K(A,B, µ j ) Z and σ i : Reflect W o (s): the closed-loop information. σ i = jw i over the region where W o (jw) is dominant V and µ j : Obtained in an (optimal) open loop sense. µ j : From an iterative rational Krylov iteration 39
40 An Iterative Rational Krylov Iteration for Controller Reduction: 1. Choose σ i = jw i, for i = 1,..., r where w i is chosen to reflect W o (jw). 2. Z = Span [ (σ 1 I A T K ) 1 c T K (σ ri A T K ) 1 c T K ] with Z T Z = I r. 3. V = Z 4. while [relative change in µ j ] > ǫ (a) A r = Z T A K V, (b) µ j λ i (A r ) for j = 1,..., r (c) V = Span [ (µ 1 I A K ) 1 b K (µ r I A K ) 1 ] b K with Z T V = I r. 5. A r = Z T A K V, b r = Z T b K, c r = c K V Z K r (s) includes the closed loop information V K r (s) is optimal in a restricted H 2 sense Π = ZVT 40
41 International Space Station Module 1R: n = 270. G(s) is lightly damped Long-lasting oscillations. K(s) is designed to remove these oscillations. n κ = x 10 3 Impulse response of the open loop system 6 x 10 3 Impulse response of the closed loop system Amplitude 2 1 Amplitude time (sec) time (sec) Reduce the order to r = 19 using iterative Rational Krylov and to r = 23 using one-sided frequency weighted balancing 41
42 FWBR: Frequency-weighted balancing with W i (s) = I and W o (s) = [I + G(s)K(s)] 1 G(s). IRK-CL: Iterative Rational Krylov - Closed Loop version: σ i reflect the weight W o (s). 30 Bode Plot of W o (s) = T(s) = [ I + G(s)K(s)] 1 G(s) Singular Values (db) Frequency (rad/sec) σ i = j logspace( 1, 2, 10) rad/sec 42
43 Singular Values (db) ISS Example: Bode Plots of reduced closed loop systems T T IRK CL T FW20 T FW23 Relative Errors H error T T FW T T FW23 T T IRK CL Frequency (rad/sec) ISS Model: Error in closed loop systems T T FW20 Relative Errors H 2 error T T FW T T FW23 T T IRK CL Singular Values (db) T T FW23 Weighted Errors H 2 error W i (K K FW20 ) < T T IRK CL W i (K K FW23 ) < 1 W i (K K IRK CL ) < Frequency (rad/sec) 43
44 An Unstable Model: n=2000. K(s) of order n κ = 2000 stabilizes the model. 70 Impulse response of G(s) Amplitude Impulse Response of T(s) Amplitude time (sec) K(s) has four unstable poles. 44
45 Reduce the order to r = 14: Stabilizing controller K r (s) has 4 unstable poles as desired Bode Plots of T(s) and T r (s) T(s) T r (s) 0.5 Impulse Response T(s) T r (s) Singular Values Amplitude Singular Values 10 0 Bode Plots of K(s) and K r (s) 10 1 K(s) K r (s) Amplitude Response to u(t)= sin (4t) T(s) T r (s) Frequency (rad/sec) time (sec) 45
46 Part II: Conclusions and Future Work n >> 10 6 : Forces usage of Inexact Solves in Krylov-based reduction Perturbation effects: Backward and forward error analysis framework Good/Optimal shift selection robust with respect to inexact solves I-IRKA (Locally) optimal reduced models for n > 10 6 without user intervention Acceleration strategies Open issues: Global H 2 and/or H perturbation effects Modifications to GMRES, effective preconditioning strategies Scalable parallel versions A large-scale easy-to-use model reduction toolbox Modify the algorithms to fit into the framework of, e.g., Trilinos Implementation on Virginia Tech.-System X 46
47 Effect of Initialization: CD Player Model n = 120 to r = 8 via IRKA (Logarithmic grid) Logarithmic Grid Sweep 8 th Order Reduced Model Relative H 2 Error Imaginary Axis Real Axis
48 Starting point: H 2 error expression (Gugercin and Antoulas [2003]) G(s) with poles λ i and G r (s) with poles λ j φ i := G(s)(s λ i ) s=λi and φ j := G r (s)(s λ j ) s= λj. G(s) G r (s) 2 H 2 = n φ i (G( λ i ) G r ( λ i )) + i=1 r j=1 ) φ j (G r ( λ j ) G( λ j ). Open loop error due to mismatch at λ i and λ j 48
49 Given Σ = The Lanczos Procedure: SISO Case [ A C B D ] = v 1 := 1 B and z 1 := sgn(cb) C T CB CB AV = V T + fet r A T Z = ZT T + ge T r, V T Z = I r, V T g = Z T f = 0 Σ r = [ Z T AV CV Z T B D ] = [ T sgn(cb) CB e T 1 Σ r is unique and matches 2r Markov parameters CB e1 D ] Im(V ) = Im ( [ B AB A r 1 B ] ) =: Im (R r (A, B)) Im(Z) = Im ( [ C T A T C T (A T ) r 1 C T ] ) =: Im ( O r (A, C) T) 49
50 Define H r := η 1 η 2 η r η 2 η 3 η r η r η r+1 η 2r 1, where η i = CA i 1 B. Lanczos requires deth i 0 for i = 1,, r. Violated for many systems Look-Ahead Lanczos deth i 0 can be avoided if tridiagonal structure not required. Only requirement: deth r 0. Arnoldi Procedure: The r th order reduced model Σ r matches only r Markov parameters and is not unique. For arbitrary σ, A (σi A) 1, B (σi A) 1, C C(σI A T ) 1 50
Optimal Approximation of Dynamical Systems with Rational Krylov Methods. Department of Mathematics, Virginia Tech, USA
Optimal Approximation of Dynamical Systems with Rational Krylov Methods Serkan Güǧercin Department of Mathematics, Virginia Tech, USA Chris Beattie and Virginia Tech. jointly with Thanos Antoulas Rice
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 informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
An iterative SVD-Krylov based method for model reduction of large-scale dynamical systems Serkan Gugercin Department of Mathematics, Virginia Tech., Blacksburg, VA, USA, 24061-0123 gugercin@math.vt.edu
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 informationChris Beattie and Serkan Gugercin
Krylov-based model reduction of second-order systems with proportional damping Chris Beattie and Serkan Gugercin Department of Mathematics, Virginia Tech, USA 44th IEEE Conference on Decision and Control
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 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 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 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 informationInexact Solves in Krylov-based Model Reduction
Inexact Solves in Krylov-based Model Reduction Christopher A. Beattie and Serkan Gugercin Abstract We investigate the use of inexact solves in a Krylov-based model reduction setting and present the resulting
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 informationarxiv: v1 [math.na] 12 Jan 2012
Interpolatory Weighted-H 2 Model Reduction Branimir Anić, Christopher A. Beattie, Serkan Gugercin, and Athanasios C. Antoulas Preprint submitted to Automatica, 2 January 22 arxiv:2.2592v [math.na] 2 Jan
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 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 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 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 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 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 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 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 informationS. Gugercin and A.C. Antoulas Department of Electrical and Computer Engineering Rice University
Proceedings of the 39" IEEE Conference on Decision and Control Sydney, Australia December, 2000 A Comparative Study of 7 Algorithms for Model Reduct ion' S. Gugercin and A.C. Antoulas Department of Electrical
More informationApproximation of the Linearized Boussinesq Equations
Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan
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 informationInexact Solves in Interpolatory Model Reduction
Inexact Solves in Interpolatory Model Reduction Sarah Wyatt Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements for the
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 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 informationComparison of Model Reduction Methods with Applications to Circuit Simulation
Comparison of Model Reduction Methods with Applications to Circuit Simulation Roxana Ionutiu, Sanda Lefteriu, and Athanasios C. Antoulas Department of Electrical and Computer Engineering, Rice University,
More informationThe rational Krylov subspace for parameter dependent systems. V. Simoncini
The rational Krylov subspace for parameter dependent systems V. Simoncini Dipartimento di Matematica, Università di Bologna valeria.simoncini@unibo.it 1 Given the continuous-time system Motivation. Model
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 informationKrylov-based model reduction of second-order systems with proportional damping
Krylov-based model reduction of second-order systems with proportional damping Christopher A Beattie and Serkan Gugercin Abstract In this note, we examine Krylov-based model reduction of second order systems
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 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 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 informationAdaptive rational Krylov subspaces for large-scale dynamical systems. V. Simoncini
Adaptive rational Krylov subspaces for large-scale dynamical systems V. Simoncini Dipartimento di Matematica, Università di Bologna valeria@dm.unibo.it joint work with Vladimir Druskin, Schlumberger Doll
More informationStructure-preserving tangential interpolation for model reduction of port-hamiltonian Systems
Structure-preserving tangential interpolation for model reduction of port-hamiltonian Systems Serkan Gugercin a, Rostyslav V. Polyuga b, Christopher Beattie a, and Arjan van der Schaft c arxiv:.3485v2
More informationL 2 -optimal model reduction for unstable systems using iterative rational
L 2 -optimal model reduction for unstable systems using iterative rational Krylov algorithm Caleb agruder, Christopher Beattie, Serkan Gugercin Abstract Unstable dynamical systems can be viewed from a
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 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 informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More informationGramians of structured systems and an error bound for structure-preserving model reduction
Gramians of structured systems and an error bound for structure-preserving model reduction D.C. Sorensen and A.C. Antoulas e-mail:{sorensen@rice.edu, aca@rice.edu} September 3, 24 Abstract In this paper
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Serkan Gugercin, Tatjana Stykel, Sarah Wyatt Model Reduction of Descriptor Systems by Interpolatory Projection Methods Preprint Nr. 3/213 29.
More informationAn Interpolation-Based Approach to the Optimal H Model Reduction
An Interpolation-Based Approach to the Optimal H Model Reduction Garret M. Flagg Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment
More informationGramian Based Model Reduction of Large-scale Dynamical Systems
Gramian Based Model Reduction of Large-scale Dynamical Systems Paul Van Dooren Université catholique de Louvain Center for Systems Engineering and Applied Mechanics B-1348 Louvain-la-Neuve, Belgium vdooren@csam.ucl.ac.be
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 informationIssues in Interpolatory Model Reduction: Inexact Solves, Second-order Systems and DAEs
Issues in Interpolatory Model Reduction: Inexact Solves, Second-order Systems and DAEs Sarah Wyatt Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
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 informationOrder reduction of large scale second-order systems using Krylov subspace methods
Linear Algebra and its Applications 415 (2006) 385 405 www.elsevier.com/locate/laa Order reduction of large scale second-order systems using Krylov subspace methods Behnam Salimbahrami, Boris Lohmann Lehrstuhl
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 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 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 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 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 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 informationFourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics
Fourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics K. Willcox and A. Megretski A new method, Fourier model reduction (FMR), for obtaining stable, accurate, low-order models
More informationSolution of eigenvalue problems. Subspace iteration, The symmetric Lanczos algorithm. Harmonic Ritz values, Jacobi-Davidson s method
Solution of eigenvalue problems Introduction motivation Projection methods for eigenvalue problems Subspace iteration, The symmetric Lanczos algorithm Nonsymmetric Lanczos procedure; Implicit restarts
More informationOrder Reduction for Large Scale Finite Element Models: a Systems Perspective
Order Reduction for Large Scale Finite Element Models: a Systems Perspective William Gressick, John T. Wen, Jacob Fish ABSTRACT Large scale finite element models are routinely used in design and optimization
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 informationKrylov-based model reduction of second-order systems with proportional damping
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 TuA05.6 Krylov-based model reduction of second-order systems
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 informationPassive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form
Passive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form Boyuan Yan, Sheldon X.-D. Tan,PuLiu and Bruce McGaughy Department of Electrical Engineering, University of
More informationAn Internal Stability Example
An Internal Stability Example Roy Smith 26 April 2015 To illustrate the concept of internal stability we will look at an example where there are several pole-zero cancellations between the controller and
More informationSiMpLIfy: A Toolbox for Structured Model Reduction
SiMpLIfy: A Toolbox for Structured Model Reduction Martin Biel, Farhad Farokhi, and Henrik Sandberg ACCESS Linnaeus Center, KTH Royal Institute of Technology Department of Electrical and Electronic Engineering,
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More informationModel reduction of large-scale systems
Model reduction of large-scale systems Lecture IV: Model reduction from measurements Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/ aca International School,
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 informationADI-preconditioned FGMRES for solving large generalized Lyapunov equations - A case study
-preconditioned for large - A case study Matthias Bollhöfer, André Eppler TU Braunschweig Institute Computational Mathematics Syrene-MOR Workshop, TU Hamburg October 30, 2008 2 / 20 Outline 1 2 Overview
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 informationInterpolatory Model Reduction of Large-scale Dynamical Systems
Interpolatory Model Reduction of Large-scale Dynamical Systems Athanasios C. Antoulas, Christopher A. Beattie, and Serkan Gugercin Abstract Large-scale simulations play a crucial role in the study of a
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 informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview Thanos Antoulas joint work with Dan Sorensen aca@rice.edu URL: www-ece.rice.edu/ aca CITI Lecture, Rice, February 2005 Model reduction of large-scale
More informationA Novel Scheme for Positive Real Balanced Truncation
A Novel Scheme for Positive Real Balanced runcation Kari Unneland, Paul Van Dooren and Olav Egeland Abstract Here, model reduction, based on balanced truncation, of stable and passive systems will be considered.
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 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 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 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 informationarxiv: v3 [math.na] 6 Jul 2018
A Connection Between Time Domain Model Order Reduction and Moment Matching for LTI Systems arxiv:181.785v3 [math.na] 6 Jul 218 Manuela Hund a and Jens Saak a a Max Planck Institute for Dynamics of Complex
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 informationLMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System
LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz
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 informationEECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total
More informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
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 informationAn Optimum Fitting Algorithm for Generation of Reduced-Order Models
An Optimum Fitting Algorithm for Generation of Reduced-Order Models M.M. Gourary 1, S.G. Rusakov 1, S.L. Ulyanov 1, M.M. Zharov 1, B.J. Mulvaney 2 1 IPPM, Russian Academy of Sciences, Moscow 1523, e-mail:
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 informationTechnische Universität Berlin
; Technische Universität Berlin Institut für Mathematik arxiv:1610.03262v2 [cs.sy] 2 Jan 2017 Model Reduction for Large-scale Dynamical Systems with Inhomogeneous Initial Conditions Christopher A. Beattie
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 informationLarge-scale and infinite dimensional dynamical systems approximation
Large-scale and infinite dimensional dynamical systems approximation Igor PONTES DUFF PEREIRA Doctorant 3ème année - ONERA/DCSD Directeur de thèse: Charles POUSSOT-VASSAL Co-encadrant: Cédric SEREN 1 What
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationProjection of state space realizations
Chapter 1 Projection of state space realizations Antoine Vandendorpe and Paul Van Dooren Department of Mathematical Engineering Université catholique de Louvain B-1348 Louvain-la-Neuve Belgium 1.0.1 Description
More informationModule 07 Controllability and Controller Design of Dynamical LTI Systems
Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October
More informationControl Systems. Frequency domain analysis. L. Lanari
Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic
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 informationControl Systems Design
ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:
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 information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 22: Balanced Realization Readings: DDV, Chapter 26 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology April 27, 2011 E. Frazzoli
More informationSolution of eigenvalue problems. Subspace iteration, The symmetric Lanczos algorithm. Harmonic Ritz values, Jacobi-Davidson s method
Solution of eigenvalue problems Introduction motivation Projection methods for eigenvalue problems Subspace iteration, The symmetric Lanczos algorithm Nonsymmetric Lanczos procedure; Implicit restarts
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 information