Optimal Approximation of Dynamical Systems with Rational Krylov Methods. Department of Mathematics, Virginia Tech, USA
|
|
- Brianna Chapman
- 5 years ago
- Views:
Transcription
1 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 University Workshop on Large-Scale Robust Optimization, August 31-September 2, 2005, Santa Fe, NM 1
2 Outline 1. Introduction and Problem Statement 2. Motivating Examples 3. Rational Krylov-Interpolation Framework 4. An Iterative Rational Krylov Algorithm 5. Inexact Solves in Krylov-based Model Reduction 3
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 = 4
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 5
5 Model reduction of Serkan from n=3 down to r=2 Cascades, Blacksburg, VA
6 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. 6
7 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) 7
8 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. 8
9 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) =
10 Motivating Example: Control Optimal Cooling of Steel Profiles in a Rolling Mill : Data: Peter Benner 10
11 Different steps in the production process require different temperatures of the raw material. To achieve high throughput, reduce the temperature as fast as possible to the required level before entering the next production phase. Cooling process by spraying cooling fluids on the surface Must be controlled so that material properties, such as durability or porosity, stay within given quality standards Modeled as boundary control of a two dimensional heat equation. A finite element discretization results in Eẋ(t) = Ax(t) + bu(t), y(t) = cx(t). n = 79, 841: nnz(a) : , nnz(e) :
12 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. Why to choose model reduction via rational interpolation? 12
13 Generically, any reduced model G r (s) can be obtained via interpolation. Interpolation points = Zeroes of G(s) G r (s). BUT: Prob-1: What is a good selection of interpolation points? Similar to polynomial approximation of complex functions. Recall: Trying to match the moments: c(σ k I A) (j+1) b Moments are extremely ill-conditioned Prob-2: Construct G r (s) without explicit moment computation Prob-2 easier to tackle using rational Krylov framework (Skelton et al. [1987], Grimme [1997]): 13
14 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, 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) s=σi Moment matching without explicit moment computation Still to answer: How to choose σ i? σ i = λ i (A) (Antoulas/G [2003]). Effective but not optimal. Does there exist an optimal selection? 14
15 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 Existence of a global minimal: Exists in the SISO case Not known for the MIMO case General approach: Find G r (s) that satisfies first-order necessary conditions: Wilson [1970], Meier and Luenburger [1967], Hyland and Bernstein [1985], Yan and Lam [1999],... 15
16 Framework of Wilson [1970] Given G r (s) = c r (si r A r ) 1 b r, define the error system G e (s) := G(s) G r (s) = c e (si A e ) 1 b e Let P e and Q e be the error gramians: P e = A e P e + P e A T e + b e b T e = 0, Q e A e + A T e Q e + c T e c e = 0 P 11 P 12 P T 12 P 22, Q e = Q 11 Q 12 Q T 12 Q 22 G e (s) 2 H 2 = c e P e c T e : = First-order necessary conditions: P T 12Q 12 + P 22 Q 22 = 0 Q T 12b + Q 22 b r = 0 c r P 22 cp 12 = 0. 16
17 Equivalently, V = P 12 P 1 22, Z = Q 12Q 1 22 and H 2 Iteration: A r = Z T AV, b r = Z T b, c r = cv. 1. Choose an initial G r (s) = c r (si r A r ) 1 b r. 2. Compute P e and Q e 3. Define V = P 12 P 1 22, Z = Q 12Q Let A r = Z T AV, b r = Z T b, c r = cv. 5. Return to Step 1. Two Lyapunov equations at each step. Similar framework by Hyland and Bernstein [1985] 17
18 Framework of Meier and Luenberger [1967] Let G r (s) = c r (si r A r ) 1 b r solves the optimal H 2 problem Let ˆλ i = λ i (A r ), i.e. the Ritz values. First-order conditions: G( ˆλ i ) = G r ( ˆλ i ), and d 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 18
19 Theorem: The two frameworks are equivalent. Proof: Starting point for Lyapunov Interpolation Framework: Lemma: (Gallivan et al. [2004], Antoulas/Sorensen [2002]) Let V solves AV + VA T r + bb T r = 0. Then, [ ] Ran(V) = Span ( ˆλ 1 I A) 1 b,, ( ˆλ r I A) 1 b. Starting point for Interpolation Lyapunov Framework: Model reduction via rational Krylov projection. 19
20 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 20
21 No methods guarantee convergence to global minimum. Question: Global minimum of a restricted H 2 minimization problem? Corollary: (Gaier 1980) Given stable G(s), and the stable reduced poles α 1,...,α r, define Ĝ(s) := β 0 + β 1 s + + β r s r (s α 1 )...(s α r ). Then G(s) Ĝ(s) H 2 is minimized if and only if G(s) = Ĝ(s) for s = α 1, α 2,..., α r. Upon convergence, IRKA minimizes the H 2 norm of the error system among all possible reduced models having the same reduced poles λ j. 21
22 Convergence? Understood better and better every day!!! A fixed point iteration: { } ({ }) (k+1) (k) σ i = f σ i Π (k+1) = h (Π (k)) Usual outcome is (numerical) convergence in 4 5 steps Convergence failure in rare circumstances. Newton Iteration Framework: Jacobian J: Sensitivity of λ i (A r ) wrt {σ i } Requires solving an r r generalized eigenvalue problem {σ i } (k+1) = {σ i } (k) (I + J) 1 ( {σ i } (k) + {λ i (A r )} (k)). 22
23 Stability? A r = Z T AV nonnormal reduced order model Reduced order stability not guaranteed in general. But, very hard to force convergence to unstable model (occasional unstable models can occur at intermediate stages) Fairly robust with respect to initial shift selection. Gugercin [CDC-2005]: Replace Z by QV(V T QV) 1 where implies stability. A T Q + QA + c T c = 0. 23
24 EXTREMELY small order benchmark examples Model r IRKA GFM OPM BTM FOM FOM FOM FOM Divergent FOM FOM Divergent FOM FOM FOM Divergent FOM FOM GFM: Gradient Flow Method of Yan and Lam [1999] OPM: Optimal Projection Method of Hyland and Bernstein [1985] BTM: Balanced Truncation Method of Moore [1981] FOM-1: n = 4, FOM-2: n = 7, FOM-3: n = 4, FOM-4: n = 2, 24
25 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 25
26 Part I: Conclusions and Future Work: Equivalence of first-order conditions for the H 2 problem Iterative Rational Krylov for optimal H 2 reduction First-order conditions while staying in Krylov framework No Lyapunov equations need to be solved Good H performance as well (Zolatorjov Problem (Beattie [2005])). Some open issues remain for convergence and stability. Newton s Iteration Formulation Application to controller reduction: Gugercin/Antoulas/Beattie [2005] Variations that guarantee stability (Gugercin [2005]) Find another name and acronym better than IRKA 26
27 Inexact Solves in Krylov-based Model Reduction Need for more detail and accuracy in the modeling stage System dimension n: O(10 6 ) or more (σi A)v = b cannot be solved directly Inexact solves need to be employed in constructing V and Z Questions: 1. What are the perturbation effects on interpolation? 2. Robustness with respect to the inexact solves? 3. What are the effective preconditioning, restarting strategies? 4. What is the effect on (the optimality of) the reduced model? 27
28 For simplicity, consider the one-sided projection, i.e. V = Z. Let ˆv j be an inexact solution for (σ j I A)v j = b (σ j I A)ˆv j b = δb j with δb j b ǫ Define δv j := ˆv j v j = (σ j I A) 1 δb j, and ˆK := [ (σ 1 I A) 1 b + δv 1, (σ r I A) 1 ] b + δv r. Inexact Krylov-based reduced model obtained by A r = V T A V, b r = V T b, c r = c V, where VT V = Ir. where V is an orthogonal basis for Range( ˆK) 28
29 Theorem: Given the above set-up, c r (σ j I r A r ) 1 b r = c(σ j I n A) 1 b + ε fwd = c(σ j I n A) 1 (b + b j ) where ε fwd = c [ (σ j I n A) 1 V(σ j I r A r ) 1 V T ] δb j. b j = [I n (σ j I n A)V(σ j I r A r ) 1 V T ]δb j. ε fwd : Forward error, b j : Backward error How well V(σ j I r A r ) 1 V T approximates (σ j I n A) 1 Expect optimal model to be robust with respect to inexact solves. Same analysis valid for the two-sided projection as well. 29
30 GMRES: 1. The same Krylov subspace for each (σ j I A)v j = b AW k = W k+1 Hk min σ j Ĩ H k b e 1 2. Span{v j } r j=1 is important, rather than each v j = min x ˆx 1,...,ˆx l (σ l+1 I A)x b 3. Two-sided case: BiCG,... Preconditioning: 1. If σ j is close to σ j+1, can re-use preconditioners for different linear systems 2. Cost of recomputing vs cost of using a close-by preconditioner 30
31 Inexact IRKA (I-IRKA) IRKA requires solving 2r linear systems at each step Expensive if n = O(10 6 ) Recall: {σ j } converge fast Use the solution from the previous step as an initial guess for the next step Expect faster convergence for a fixed tolerance value Optimal reduced model: Expect robustness 31
32 Example: Optimal Cooling of Steel Profiles ( P. Benner) G(s) = c(se A) 1 b, n = 20, 209 Bad shift selection: σ i = logspace( 8, 4, 6) r = 6 via Rational Krylov (RK) and Inexact-RK (I-RK). I-RK uses GMRES with tol = Bode Plots of H(s), H 1 (s) (RK), and H 2 (s) (I IRKA) H(s) H 1 (s) H 2 (s) Error Errors in the computed first moments Relative error Absolute error H(jw) 10 3 H(s) H 1 (s) = H(s) H 2 (s) = H 1 (s) H 2 (s) = Error Errors in the computed second moments 10 0 Absolute error Relative error frequency (rad/sec) σ 32
33 Optimal {σ i } obtained via IRKA Use these {σ i } in I-RK. I-RK uses GMRES with tol = Bode Plots of H(s), H 1 (s) (IRKA), and H 2 (s) (I IRKA) H(s) H 1 (s) H 2 (s) 10 5 Errors in the computed first moments Relative error 10 2 Error Absolute error H(jw) H(s) H 1 (s) = H(s) H 1 (s) = H 1 (s) H 2 (s) = Error Errors in the computed second moments 10 0 Relative error Absolute error frequency (rad/sec) σ 33
34 Same model with n = 79, 841 (Finer discretization) r = 6 via IRKA and I IRKA (tol = ) IRKA: Initial guess from the previous step H(jw) Amplitude Bode plots of H(s), H 1 (s) and H 2 (s) H(s) H 1 (s) H 2 (s) frequency (rad/sec) Number of GMRES steps Total Number of GMRES steps at the k th Iteration Iteration Index H(s) H 1 (s) = H(s) H 2 (s) = , H 1 (s) H 2 (s) =
35 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 35
36 Alexei Nikolaevich Krylov [1931] "On the numerical solution of the equation by which, in technical matters, frequencies of small oscillations of material systems are determined" to compute the characteristic polynomial coefficients.
37 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. 37
38 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). = 38
39 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. 39
40 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 40
41 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 41
42 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 42
43 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 43
44 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) 44
45 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. 45
46 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) 46
Rational Krylov Methods for Model Reduction of Large-scale Dynamical Systems
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
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
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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-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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationRealization-independent H 2 -approximation
Realization-independent H -approximation Christopher Beattie and Serkan Gugercin Abstract Iterative Rational Krylov Algorithm () of [] is an effective tool for tackling the H -optimal model reduction problem.
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 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 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 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 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 informationComputing Transfer Function Dominant Poles of Large Second-Order Systems
Computing Transfer Function Dominant Poles of Large Second-Order Systems Joost Rommes Mathematical Institute Utrecht University rommes@math.uu.nl http://www.math.uu.nl/people/rommes joint work with Nelson
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 informationarxiv: v1 [cs.sy] 17 Nov 2015
Optimal H model approximation based on multiple input/output delays systems I. Pontes Duff a, C. Poussot-Vassal b, C. Seren b a ISAE SUPAERO & ONERA b ONERA - The French Aerospace Lab, F-31055 Toulouse,
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 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 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 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 informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. MATRIX ANAL. APPL. Vol. 38, No. 2, pp. 401 424 c 2017 Society for Industrial and Applied Mathematics PRECONDITIONED MULTISHIFT BiCG FOR H 2 -OPTIMAL MODEL REDUCTION MIAN ILYAS AHMAD, DANIEL B.
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 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 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 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 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 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 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 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 informationTopic # Feedback Control Systems
Topic #20 16.31 Feedback Control Systems Closed-loop system analysis Bounded Gain Theorem Robust Stability Fall 2007 16.31 20 1 SISO Performance Objectives Basic setup: d i d o r u y G c (s) G(s) n control
More informationModel Reduction for Dynamical Systems
MAX PLANCK INSTITUTE Otto-von-Guericke Universitaet Magdeburg Faculty of Mathematics Summer term 2015 Model Reduction for Dynamical Systems Lecture 10 Peter Benner Lihong Feng Max Planck Institute for
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 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 informationMultivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006
Multivariable Control Lecture 3 Description of Linear Time Invariant Systems John T. Wen September 7, 26 Outline Mathematical description of LTI Systems Ref: 3.1-3.4 of text September 7, 26Copyrighted
More informationLMIs for Observability and Observer Design
LMIs for Observability and Observer Design Matthew M. Peet Arizona State University Lecture 06: LMIs for Observability and Observer Design Observability Consider a system with no input: ẋ(t) = Ax(t), x(0)
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 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 informationPole placement control: state space and polynomial approaches Lecture 2
: state space and polynomial approaches Lecture 2 : a state O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename -based November 21, 2017 Outline : a state
More informationHigh Performance Nonlinear Solvers
What is a nonlinear system? High Performance Nonlinear Solvers Michael McCourt Division Argonne National Laboratory IIT Meshfree Seminar September 19, 2011 Every nonlinear system of equations can be described
More informationEEE 184: Introduction to feedback systems
EEE 84: Introduction to feedback systems Summary 6 8 8 x 7 7 6 Level() 6 5 4 4 5 5 time(s) 4 6 8 Time (seconds) Fig.. Illustration of BIBO stability: stable system (the input is a unit step) Fig.. step)
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 informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
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 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 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 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 informationRobust Control 3 The Closed Loop
Robust Control 3 The Closed Loop Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /2/2002 Outline Closed Loop Transfer Functions Traditional Performance Measures Time
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
More informationDownloaded 05/27/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 35, No. 5, pp. B1010 B1033 c 2013 Society for Industrial and Applied Mathematics MODEL REDUCTION OF DESCRIPTOR SYSTEMS BY INTERPOLATORY PROJECTION METHODS SERKAN GUGERCIN, TATJANA
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 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 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 informationH 2 optimal model approximation by structured time-delay reduced order models
H 2 optimal model approximation by structured time-delay reduced order models I. Pontes Duff, Charles Poussot-Vassal, C. Seren 27th September, GT MOSAR and GT SAR meeting 200 150 100 50 0 50 100 150 200
More informationStability of CL System
Stability of CL System Consider an open loop stable system that becomes unstable with large gain: At the point of instability, K( j) G( j) = 1 0dB K( j) G( j) K( j) G( j) K( j) G( j) =± 180 o 180 o Closed
More informationPreconditioned inverse iteration and shift-invert Arnoldi method
Preconditioned inverse iteration and shift-invert Arnoldi method Melina Freitag Department of Mathematical Sciences University of Bath CSC Seminar Max-Planck-Institute for Dynamics of Complex Technical
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 informationPreconditioned multishift BiCG for H 2 -optimal model reduction. Mian Ilyas Ahmad, Daniel B. Szyld, and Martin B. van Gijzen
Preconditioned multishift BiCG for H 2 -optimal model reduction Mian Ilyas Ahmad, Daniel B. Szyld, and Martin B. van Gijzen Report 12-06-15 Department of Mathematics Temple University June 2012. Revised
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 informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationModel Order Reduction via Matlab Parallel Computing Toolbox. Istanbul Technical University
Model Order Reduction via Matlab Parallel Computing Toolbox E. Fatih Yetkin & Hasan Dağ Istanbul Technical University Computational Science & Engineering Department September 21, 2009 E. Fatih Yetkin (Istanbul
More informationNetwork Clustering for SISO Linear Dynamical Networks via Reaction-Diffusion Transformation
Milano (Italy) August 28 - September 2, 211 Network Clustering for SISO Linear Dynamical Networks via Reaction-Diffusion Transformation Takayuki Ishizaki Kenji Kashima Jun-ichi Imura Kazuyuki Aihara Graduate
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 informationTopic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback
Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall
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 informationRecent advances in approximation using Krylov subspaces. V. Simoncini. Dipartimento di Matematica, Università di Bologna.
Recent advances in approximation using Krylov subspaces V. Simoncini Dipartimento di Matematica, Università di Bologna and CIRSA, Ravenna, Italy valeria@dm.unibo.it 1 The framework It is given an operator
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 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 informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
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 informationAlternative correction equations in the Jacobi-Davidson method
Chapter 2 Alternative correction equations in the Jacobi-Davidson method Menno Genseberger and Gerard Sleijpen Abstract The correction equation in the Jacobi-Davidson method is effective in a subspace
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 22: H 2, LQG and LGR Conclusion To solve the H -optimal state-feedback problem, we solve min γ such that γ,x 1,Y 1,A n,b n,c n,d
More information