Reduced Basis Method for Parametric
|
|
- Ferdinand Clyde Payne
- 5 years ago
- Views:
Transcription
1 1/24 Reduced Basis Method for Parametric H 2 -Optimal Control Problems NUMOC2017 Andreas Schmidt, Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation, University of Stuttgart andreas.schmidt@mathematik.uni-stuttgart.de 22 nd June, 2017
2 Outline 1 Introduction 2 Reduced Basis Approximation 3 Numerical Results 4 Conclusion and Outlook 2/24
3 Introduction 3/24
4 Basic Overview A general control problem Goal: Set up a general framework for realistic control problems. Noise w 1 (t) Control u(t) System x(t) State z(t) Performance y(t) Measurement Noise w 2 (t) In a more mathematical way: LTI System x(t) = Ax(t) +B u u(t) +B w w 1 (t) y(t) = C y x(t) +D yw w 2 (t) z(t) = C z x(t) +D zu u(t) x(0) = x 0. State equation Measurements Performance output 4/24
5 Basic Overview A general control problem Goal: Set up a general framework for realistic control problems. Noise w 1 (t) Control u(t) System x(t) State z(t) Performance y(t) Measurement Noise w 2 (t) In a more mathematical way: LTI System x(t) = Ax(t) +B u u(t) +B w w 1 (t) y(t) = C y x(t) +D yw w 2 (t) z(t) = C z x(t) +D zu u(t) x(0) = x 0. State equation Measurements Performance output 4/24
6 Basic Overview A general control problem Goal: Set up a general framework for realistic control problems. Noise w 1 (t) Control u(t) System x(t) State z(t) Performance y(t) Measurement Noise w 2 (t) In a more mathematical way: LTI System x(t) = Ax(t) +B u u(t) +B w w 1 (t) y(t) = C y x(t) +D yw w 2 (t) z(t) = C z x(t) +D zu u(t) x(0) = x 0. State equation Measurements Performance output 4/24
7 Basic Overview A general control problem Goal: Set up a general framework for realistic control problems. Noise w 1 (t) Control u(t) System x(t) State z(t) Performance y(t) Measurement Noise w 2 (t) In a more mathematical way: LTI System x(t) = A(µ)x(t) +B u (µ)u(t) +B w (µ)w 1 (t) y(t) = C y (µ)x(t) +D yw (µ)w 2 (t) z(t) = C z (µ)x(t) +D zu (µ)u(t) x(0) = x 0 (µ). State equation Measurements Performance output 4/24
8 The Control Problem What should the controller do? w(t) System z(t) u(t) y(t) Controller Control objective Try to find a controller that maps the noisy measurements y(t) to control signals u(t) such that certain goals are achieved: System is stable Effect of noise w(t) on output z(t) is minimized in H 2 norm 1 G 2 = 2π G(iω) 2 F dω, where G(s) = C (si A ) 1 B 5/24
9 The Control Problem What should the controller do? w(t) u(t) Control objective System Controller y(t) z(t) Controller: x K = A K x K +B K y u = C K x K Closed-loop LTI: h = A h + Bw z = C h Try to find a controller that maps the noisy measurements y(t) to control signals u(t) such that certain goals are achieved: System is stable Effect of noise w(t) on output z(t) is minimized in H 2 norm 1 G 2 = 2π G(iω) 2 F dω, where G(s) = C (si A ) 1 B 5/24
10 The Control Problem What should the controller do? w(t) u(t) Control objective System Controller y(t) z(t) Controller: x K = A K x K +B K y u = C K x K Closed-loop LTI: h = A h + Bw z = C h Try to find a controller that maps the noisy measurements y(t) to control signals u(t) such that certain goals are achieved: System is stable Effect of noise w(t) on output z(t) is minimized in H 2 norm 1 G 2 = 2π G(iω) 2 F dω, where G(s) = C (si A ) 1 B 5/24
11 The Control Problem What should the controller do? w(t) u(t) Control objective System Controller y(t) z(t) Controller: x K = A K x K +B K y u = C K x K Closed-loop LTI: h = A h + Bw z = C h Try to find a controller that maps the noisy measurements y(t) to control signals u(t) such that certain goals are achieved: System is stable Effect of noise w(t) on output z(t) is minimized in H 2 norm 1 G 2 = 2π G(iω) 2 F dω, where G(s) = C (si A ) 1 B 5/24
12 The Control Problem What should the controller do? w(t) u(t) Control objective System Controller y(t) z(t) Controller: x K = A K x K +B K y u = C K x K Closed-loop LTI: h = A h + Bw z = C h Try to find a controller that maps the noisy measurements y(t) to control signals u(t) such that certain goals are achieved: System is stable Effect of noise w(t) on output z(t) is minimized in H 2 norm 1 G 2 = 2π G(iω) 2 F dω, where G(s) = C (si A ) 1 B 5/24
13 The Control Problem The solution Theorem 1 ( [Zhou, Doyle, Glover 1996]) Let P,Q R n n be the unique stabilizing solutions to the algebraic Riccati equations (AREs) A T P +PA PB u (D T zud zu ) 1 B T u P +C T z C z = 0 (1) AQ +QA T QCy T (D yw Dyw) T 1 C y Q +B w B T w = 0. (2) Define the control gain K := (D T zud zu ) 1 B T u P and the observer gain L := QC T y (D yw D T yw) 1. Then the H 2 optimal controller is given by u(t) := Kx K (t), (3) where x K (t) R n is the solution to the observer-equation: x K (t) = (A BK)x K (t) +L[y(t) C y x K (t)]. (4) K. Zhou, J. Doyle and K. Glover. Robust and Optimal Control /24
14 The Control Problem The solution Theorem 1 ( [Zhou, Doyle, Glover 1996]) Let P,Q R n n be the unique stabilizing solutions to the algebraic Riccati equations (AREs) A T P +PA PB u (D T zud zu ) 1 B T u P +C T z C z = 0 (1) AQ +QA T QCy T (D yw Dyw) T 1 C y Q +B w B T w = 0. (2) Define the control gain K := (D T zud zu ) 1 B T u P and the observer gain L := QC T y (D yw D T yw) 1. Then the H 2 optimal controller is given by u(t) := Kx K (t), (3) where x K (t) R n is the solution to the observer-equation: x K (t) = (A BK)x K (t) +L[y(t) C y x K (t)]. (4) K. Zhou, J. Doyle and K. Glover. Robust and Optimal Control /24
15 The Control Problem The solution Theorem 1 ( [Zhou, Doyle, Glover 1996]) Let P,Q R n n be the unique stabilizing solutions to the algebraic Riccati equations (AREs) A T P +PA PB u (D T zud zu ) 1 B T u P +C T z C z = 0 (1) AQ +QA T QCy T (D yw Dyw) T 1 C y Q +B w B T w = 0. (2) Define the control gain K := (D T zud zu ) 1 B T u P and the observer gain L := QC T y (D yw D T yw) 1. Then the H 2 optimal controller is given by u(t) := Kx K (t), (3) where x K (t) R n is the solution to the observer-equation: x K (t) = (A BK)x K (t) +L[y(t) C y x K (t)]. (4) K. Zhou, J. Doyle and K. Glover. Robust and Optimal Control /24
16 Reduced Basis Approximation 7/24
17 Reduced Basis Method in a Nutshell Calculation is split in offline and online phase Online efficiency through parameter separability Basis generation (usually) by Greedy procedure Many applications: Elliptic/parabolic PDEs, variational inequalities, dyn. systems, Kalman filter,... A. Patera and G. Rozza. Reduced Basis Approximation and a Posteriori Error Estimation for Parametrized Partial Differential Equations. To appear in (tentative) MIT Pappalardo Graduate Monographs in Mechanical Engineering B. Haasdonk. Reduced Basis Methods for Parametrized PDEs A Tutorial Introduction for Stationary and Instationary Problems. SimTech Preprint. Chapter to appear in P. Benner, A. Cohen, M. Ohlberger and K. Willcox (eds.): "Model Reduction and Approximation: Theory and Algorithms", SIAM, Philadelphia, IANS, University of Stuttgart, Germany, /24
18 Reduced Basis Method in a Nutshell Calculation is split in offline and online phase Online efficiency through parameter separability Basis generation (usually) by Greedy procedure Many applications: Elliptic/parabolic PDEs, variational inequalities, dyn. systems, Kalman filter,... A. Patera and G. Rozza. Reduced Basis Approximation and a Posteriori Error Estimation for Parametrized Partial Differential Equations. To appear in (tentative) MIT Pappalardo Graduate Monographs in Mechanical Engineering B. Haasdonk. Reduced Basis Methods for Parametrized PDEs A Tutorial Introduction for Stationary and Instationary Problems. SimTech Preprint. Chapter to appear in P. Benner, A. Cohen, M. Ohlberger and K. Willcox (eds.): "Model Reduction and Approximation: Theory and Algorithms", SIAM, Philadelphia, IANS, University of Stuttgart, Germany, /24
19 Reduced Basis Method in a Nutshell Calculation is split in offline and online phase Online efficiency through parameter separability Basis generation (usually) by Greedy procedure Many applications: Elliptic/parabolic PDEs, variational inequalities, dyn. systems, Kalman filter,... A. Patera and G. Rozza. Reduced Basis Approximation and a Posteriori Error Estimation for Parametrized Partial Differential Equations. To appear in (tentative) MIT Pappalardo Graduate Monographs in Mechanical Engineering B. Haasdonk. Reduced Basis Methods for Parametrized PDEs A Tutorial Introduction for Stationary and Instationary Problems. SimTech Preprint. Chapter to appear in P. Benner, A. Cohen, M. Ohlberger and K. Willcox (eds.): "Model Reduction and Approximation: Theory and Algorithms", SIAM, Philadelphia, IANS, University of Stuttgart, Germany, /24
20 Reduced Basis Method in a Nutshell Calculation is split in offline and online phase Online efficiency through parameter separability Basis generation (usually) by Greedy procedure Many applications: Elliptic/parabolic PDEs, variational inequalities, dyn. systems, Kalman filter,... A. Patera and G. Rozza. Reduced Basis Approximation and a Posteriori Error Estimation for Parametrized Partial Differential Equations. To appear in (tentative) MIT Pappalardo Graduate Monographs in Mechanical Engineering B. Haasdonk. Reduced Basis Methods for Parametrized PDEs A Tutorial Introduction for Stationary and Instationary Problems. SimTech Preprint. Chapter to appear in P. Benner, A. Cohen, M. Ohlberger and K. Willcox (eds.): "Model Reduction and Approximation: Theory and Algorithms", SIAM, Philadelphia, IANS, University of Stuttgart, Germany, /24
21 Reduction of the ARE How to do that? Recall: The ARE is a nonlinear matrix-valued equation for X(µ) R n n : A(µ) T X(µ) +X(µ)A(µ) X(µ)B(µ)R(µ) 1 B(µ) T X(µ) +C(µ) T C(µ) = 0. SVD/Eigenvalue decomposition X(µ) = V ΣV T σi Index i Projection leads to small N n dimensional ARE Idea: [S., Haasdonk 2017] Approximate X(µ) X(µ) := V X X N (µ)v T X A N (µ) T X N (µ)+x N (µ)a N (µ) X N (µ)b N (µ)r(µ) 1 B(µ) T X N (µ)+c N (µ) T C N (µ) = 0. A. Schmidt and B. Haasdonk. Reduced basis approximation of large scale parametric algebraic Riccati equations. ESAIM: Control, Optimisation and Calculus of Variations (Feb. 2017). V. Simoncini. Analysis of the Rational Krylov Subspace Projection Method for Large-Scale Algebraic Riccati Equations. SIAM Journal on Matrix Analysis and Applications (Jan. 2016). 9/24
22 Reduction of the ARE How to do that? Recall: The ARE is a nonlinear matrix-valued equation for X(µ) R n n : A(µ) T X(µ) +X(µ)A(µ) X(µ)B(µ)R(µ) 1 B(µ) T X(µ) +C(µ) T C(µ) = 0. SVD/Eigenvalue decomposition X(µ) = V ΣV T σi Index i Projection leads to small N n dimensional ARE Idea: [S., Haasdonk 2017] Approximate X(µ) X(µ) := V X X N (µ)v T X A N (µ) T X N (µ)+x N (µ)a N (µ) X N (µ)b N (µ)r(µ) 1 B(µ) T X N (µ)+c N (µ) T C N (µ) = 0. A. Schmidt and B. Haasdonk. Reduced basis approximation of large scale parametric algebraic Riccati equations. ESAIM: Control, Optimisation and Calculus of Variations (Feb. 2017). V. Simoncini. Analysis of the Rational Krylov Subspace Projection Method for Large-Scale Algebraic Riccati Equations. SIAM Journal on Matrix Analysis and Applications (Jan. 2016). 9/24
23 Reduction of the ARE How to do that? Recall: The ARE is a nonlinear matrix-valued equation for X(µ) R n n : A(µ) T X(µ) +X(µ)A(µ) X(µ)B(µ)R(µ) 1 B(µ) T X(µ) +C(µ) T C(µ) = 0. SVD/Eigenvalue decomposition X(µ) = V ΣV T Z(µ)Z(µ) T σi Index i Projection leads to small N n dimensional ARE Idea: [S., Haasdonk 2017] Approximate X(µ) X(µ) := V X X N (µ)v T X A N (µ) T X N (µ)+x N (µ)a N (µ) X N (µ)b N (µ)r(µ) 1 B(µ) T X N (µ)+c N (µ) T C N (µ) = 0. A. Schmidt and B. Haasdonk. Reduced basis approximation of large scale parametric algebraic Riccati equations. ESAIM: Control, Optimisation and Calculus of Variations (Feb. 2017). V. Simoncini. Analysis of the Rational Krylov Subspace Projection Method for Large-Scale Algebraic Riccati Equations. SIAM Journal on Matrix Analysis and Applications (Jan. 2016). 9/24
24 Reduction of the ARE How to do that? Recall: The ARE is a nonlinear matrix-valued equation for X(µ) R n n : A(µ) T X(µ) +X(µ)A(µ) X(µ)B(µ)R(µ) 1 B(µ) T X(µ) +C(µ) T C(µ) = 0. SVD/Eigenvalue decomposition X(µ) = V ΣV T Z(µ)Z(µ) T σi Index i Projection leads to small N n dimensional ARE Idea: [S., Haasdonk 2017] Approximate X(µ) X(µ) := V X X N (µ)v T X A N (µ) T X N (µ)+x N (µ)a N (µ) X N (µ)b N (µ)r(µ) 1 B(µ) T X N (µ)+c N (µ) T C N (µ) = 0. A. Schmidt and B. Haasdonk. Reduced basis approximation of large scale parametric algebraic Riccati equations. ESAIM: Control, Optimisation and Calculus of Variations (Feb. 2017). V. Simoncini. Analysis of the Rational Krylov Subspace Projection Method for Large-Scale Algebraic Riccati Equations. SIAM Journal on Matrix Analysis and Applications (Jan. 2016). 9/24
25 Reduction of the ARE How to do that? Recall: The ARE is a nonlinear matrix-valued equation for X(µ) R n n : A(µ) T X(µ) +X(µ)A(µ) X(µ)B(µ)R(µ) 1 B(µ) T X(µ) +C(µ) T C(µ) = 0. SVD/Eigenvalue decomposition X(µ) = V ΣV T Z(µ)Z(µ) T σi Index i Projection leads to small N n dimensional ARE Idea: [S., Haasdonk 2017] Approximate X(µ) X(µ) := V X X N (µ)v T X A N (µ) T X N (µ)+x N (µ)a N (µ) X N (µ)b N (µ)r(µ) 1 B(µ) T X N (µ)+c N (µ) T C N (µ) = 0. A. Schmidt and B. Haasdonk. Reduced basis approximation of large scale parametric algebraic Riccati equations. ESAIM: Control, Optimisation and Calculus of Variations (Feb. 2017). V. Simoncini. Analysis of the Rational Krylov Subspace Projection Method for Large-Scale Algebraic Riccati Equations. SIAM Journal on Matrix Analysis and Applications (Jan. 2016). 9/24
26 Low-Rank Factor Greedy Algorithm Algorithm 1: Low-Rank Factor Greedy Algorithm (LRFG) for the basis generation. Data: Initial basis matrix V 0, training set P train P, tolerance ε, inner tolerance tol i [0,1], error indicator (V,µ) Set V := V 0. while max µ Ptrain (V,µ) > ε do µ := argmax µ Ptrain (V,µ) Solve the full dimensional ARE for the low rank factor Z(µ ) Set Z := (I n VV T )Z(µ ) Ẑ = POD(Z,tol i ) Extend the current basis matrix V = (V,Ẑ) Return V. 10/24
27 Low-Rank Factor Greedy Algorithm Algorithm 2: Low-Rank Factor Greedy Algorithm (LRFG) for the basis generation. Data: Initial basis matrix V 0, training set P train P, tolerance ε, inner tolerance tol i [0,1], error indicator (V,µ) Set V := V 0. while max µ Ptrain (V,µ) > ε do µ := argmax µ Ptrain (V,µ) Solve the full dimensional ARE for the low rank factor Z(µ ) Set Z := (I n VV T )Z(µ ) Ẑ = POD(Z,tol i ) Extend the current basis matrix V = (V,Ẑ) Return V. 10/24
28 Low-Rank Factor Greedy Algorithm Algorithm 3: Low-Rank Factor Greedy Algorithm (LRFG) for the basis generation. Data: Initial basis matrix V 0, training set P train P, tolerance ε, inner tolerance tol i [0,1], error indicator (V,µ) Set V := V 0. while max µ Ptrain (V,µ) > ε do µ := argmax µ Ptrain (V,µ) Solve the full dimensional ARE for the low rank factor Z(µ ) Set Z := (I n VV T )Z(µ ) Ẑ = POD(Z,tol i ) Extend the current basis matrix V = (V,Ẑ) Return V. (W,µ) := R( P(µ)) F C(µ) T Q(µ)C(µ) F 10/24
29 Low-Rank Factor Greedy Algorithm Algorithm 4: Low-Rank Factor Greedy Algorithm (LRFG) for the basis generation. Data: Initial basis matrix V 0, training set P train P, tolerance ε, inner tolerance tol i [0,1], error indicator (V,µ) Set V := V 0. while max µ Ptrain (V,µ) > ε do µ := argmax µ Ptrain (V,µ) Solve the full dimensional ARE for the low rank factor Z(µ ) Set Z := (I n VV T )Z(µ ) Ẑ = POD(Z,tol i ) Extend the current basis matrix V = (V,Ẑ) Return V. (W,µ) := R( P(µ)) F C(µ) T Q(µ)C(µ) F 10/24
30 Low-Rank Factor Greedy Algorithm Algorithm 5: Low-Rank Factor Greedy Algorithm (LRFG) for the basis generation. Data: Initial basis matrix V 0, training set P train P, tolerance ε, inner tolerance tol i [0,1], error indicator (V,µ) Set V := V 0. while max µ Ptrain (V,µ) > ε do µ := argmax µ Ptrain (V,µ) Solve the full dimensional ARE for the low rank factor Z(µ ) Set Z := (I n VV T )Z(µ ) Ẑ = POD(Z,tol i ) Extend the current basis matrix V = (V,Ẑ) Return V. (W,µ) := R( P(µ)) F C(µ) T Q(µ)C(µ) F POD(X,tol i ) returns first l SVs, such that: li=1 σ 2 i tol i (> 0.99) rank(x) i=1 σ 2 i 10/24
31 What is the overall idea? The observer equation The state estimation involves an n-dimensional ODE: x K (t;µ) = (A(µ) B(µ)K(µ))x K (t;µ) +L(µ)[y(t) C y (µ)x K (t;µ)]. Apply POD by calculating snapshots for some µ P and inputs y(t): X := [x K (t 1 ;µ 1 ),x K (t 2 ;µ 1 ),...,x K (t l ;µ P )] Figure: Some bases for an advection diffusion equation. S. Volkwein. Lecture Notes: Proper Orthogonal Decomposition: Theory and Reduced-Order Modelling /24
32 What is the overall idea? The observer equation The state estimation involves an n-dimensional ODE: x K (t;µ) = (A(µ) B(µ)K(µ))x K (t;µ) +L(µ)[y(t) C y (µ)x K (t;µ)]. Apply POD by calculating snapshots for some µ P and inputs y(t): X := [x K (t 1 ;µ 1 ),x K (t 2 ;µ 1 ),...,x K (t l ;µ P )] Figure: Some bases for an advection diffusion equation. S. Volkwein. Lecture Notes: Proper Orthogonal Decomposition: Theory and Reduced-Order Modelling /24
33 The Reduced Controller Or: Why is now everything faster? Build three different bases V P,V Q and V xk (independently, offline). See[S., Haasdonk 2016]. Definition 2 (Reduced controller) Let P NP and Q NQ be the solutions of the reduced AREs and define the reduced gains K := (D T zud zu ) 1 B T u P, L := QC T y (D yw D T yw) 1, where and the control signal P = V P P NP V T P, Q = V Q Q NQ V T Q, û(t) := K x(t), where x(t) R N x satisfies the reduced observer equation x = V T x K (A B K LC y )V xk x +V T x K Ly A. Schmidt and B. Haasdonk. Reduced basis method for H2 optimal feedback control problems. IFAC-PapersOnLine (2016). 2nd IFAC Workshop on Control of Systems Governed by Partial Differential Equations CPDE /24
34 The Reduced Controller Or: Why is now everything faster? Build three different bases V P,V Q and V xk (independently, offline). See[S., Haasdonk 2016]. Definition 2 (Reduced controller) Let P NP and Q NQ be the solutions of the reduced AREs and define the reduced gains K := (D T zud zu ) 1 B T u P, L := QC T y (D yw D T yw) 1, where and the control signal P = V P P NP V T P, Q = V Q Q NQ V T Q, û(t) := K x(t), where x(t) R N x satisfies the reduced observer equation x = V T x K (A B K LC y )V xk x +V T x K Ly A. Schmidt and B. Haasdonk. Reduced basis method for H2 optimal feedback control problems. IFAC-PapersOnLine (2016). 2nd IFAC Workshop on Control of Systems Governed by Partial Differential Equations CPDE /24
35 Numerical Results 13/24
36 Distributed Control of Damped 2D Wave Consider for Ω := [0,1] 2 the following damped wave-equation with parameters µ 1,µ 2 [0.1,2] [1,100]. f tt 0.1 f + µ 1 f t = 10 1 Ωu u +3 1 Ωw w 1, t > 0, f (0,ξ;µ) = sin(ξ 1 π)sin(ξ 2 π), ξ Ω f t (0,ξ;µ) = 0, ξ Ω, s(t;µ) = 1 f (t,ξ;µ)dξ +0.05w 2 (t). Ω y Ω y Consider two performance outputs: o 1 (t) = 1 µ 2 f (t,ξ;µ)dξ, Ω y Ω y Ω w Ω u Ω y Output o1(t) o 2 (t) = 0.1u(t) Uncontrolled Controlled /24
37 Distributed Control of Damped 2D Wave Basis generation results The equations are discretized in space by using FD: n = 800. We construct the three bases: V P,V Q and V x. The observer system was simulated with y(t) = sin(t 2 ) with t [0,2π]. t full [s] t red [s] Basis size Feedback ARE P Observer ARE Q State Estimation ODE x Table: Basis generation results. Overall RB approximation Overall controller is ODE of dimension 25 instead of 800 and the speed up for the AREs is in the magnitude of (using care). By using parameter separability, e.g. A(µ) = Q A i=1 θ i(µ)a i the online simulation can be implemented independent of n. 15/24
38 Distributed Control of Damped 2D Wave Approximation of control signal Setting: Full basis for observer gain basis and state estimation basis x Control u(t) Full N P =5 N P =10 N P = Time t 16/24
39 Distributed Control of Damped 2D Wave Approximation of control signal Setting: Full basis for feedback and observer gain basis Control u(t) 2 Full 4 N x =5 N x =10 N x = Time t 17/24
40 Distributed Control of Damped 2D Wave Approximation of the H 2 -norm Relative maximum error over test set P test with P test = 50. G( ;µ) Ĝ( ;µ) 2 max. (5) µ P test G( ;µ) 2 Table: Relative error in the transfer functions from w to z in the closed loop systems. N Q N P e e e e e e e e e e e e e e e e-04 Surprise: No instabilities occured! See famous one-page article[doyle 78]. J. Doyle. Guaranteed margins for LQG regulators. Automatic Control, IEEE Transactions on (Aug. 1978). 18/24
41 Conclusion and Outlook 19/24
42 Conclusion and Outlook We have seen: Model order reduction for H 2 optimal control problems. Full parametric and realistic control setup Expensive due to two AREs and state estimation Model reduction for AREs and the state estimation Large speed-up What next? Stability considerations and error estimation Robustification? Same approach for H -control problems 20/24
43 Conclusion and Outlook We have seen: Model order reduction for H 2 optimal control problems. Full parametric and realistic control setup Expensive due to two AREs and state estimation Model reduction for AREs and the state estimation Large speed-up What next? Stability considerations and error estimation Robustification? Same approach for H -control problems 20/24
44 Thank you for your attention! Questions? Comments? 21/24
45 References I J. Doyle. Guaranteed margins for LQG regulators. Automatic Control, IEEE Transactions on 23.4 (Aug. 1978). B. Haasdonk. Reduced Basis Methods for Parametrized PDEs A Tutorial Introduction for Stationary and Instationary Problems. SimTech Preprint. Chapter to appear in P. Benner, A. Cohen, M. Ohlberger and K. Willcox (eds.): "Model Reduction and Approximation: Theory and Algorithms", SIAM, Philadelphia, IANS, University of Stuttgart, Germany, A. Patera and G. Rozza. Reduced Basis Approximation and a Posteriori Error Estimation for Parametrized Partial Differential Equations A. Schmidt and B. Haasdonk. Reduced basis method for H2 optimal feedback control problems. IFAC-PapersOnLine 49.8 (2016). 2nd IFAC Workshop on Control of Systems Governed by Partial Differential Equations CPDE A. Schmidt and B. Haasdonk. Reduced basis approximation of large scale parametric algebraic Riccati equations. ESAIM: Control, Optimisation and Calculus of Variations (Feb. 2017). V. Simoncini. Analysis of the Rational Krylov Subspace Projection Method for Large-Scale Algebraic Riccati Equations. SIAM Journal on Matrix Analysis and Applications 37.4 (Jan. 2016). 22/24
46 References II S. Volkwein. Lecture Notes: Proper Orthogonal Decomposition: Theory and Reduced-Order Modelling K. Zhou, J. Doyle and K. Glover. Robust and Optimal Control /24
47 The Control Problem Mathematical model of the controller Consider the open-loop system x = Ax +B u u +B w w z = C z x +D zu u y = C y x +D yw w Ansatz for controller w(t) u(t) System Controller y(t) z(t) x K = A K x K +B K y u = C K x K The interconnection (i.e. the closed-loop system) for h T = (x,x K ) can be written as h = A h + Bw, z = C h. 24/24
48 The Control Problem Mathematical model of the controller Consider the open-loop system x = Ax +B u u +B w w z = C z x +D zu u y = C y x +D yw w Ansatz for controller w(t) u(t) System Controller y(t) z(t) x K = A K x K +B K y u = C K x K The interconnection (i.e. the closed-loop system) for h T = (x,x K ) can be written as h = A h + Bw, z = C h. 24/24
49 The Control Problem Mathematical model of the controller Consider the open-loop system x = Ax +B u u +B w w z = C z x +D zu u y = C y x +D yw w Ansatz for controller w(t) u(t) System Controller y(t) z(t) x K = A K x K +B K y u = C K x K The interconnection (i.e. the closed-loop system) for h T = (x,x K ) can be written as h = A h + Bw, z = C h. 24/24
Reduced basis method for H 2 optimal feedback control problems
A. Schmidt a B. Haasdonk b Reduced basis method for H 2 optimal feedback control problems Stuttgart, November 215 a Institute for Applied Analysis and Numerical Simulation, University of Stuttgart, Pfaffenwaldring
More informationParameterized Partial Differential Equations and the Proper Orthogonal D
Parameterized Partial Differential Equations and the Proper Orthogonal Decomposition Stanford University February 04, 2014 Outline Parameterized PDEs The steady case Dimensionality reduction Proper orthogonal
More informationSchmidt, Andreas 1 and Haasdonk, Bernard 1
ESAIM: Control, Optimisation and Calculus of Variations URL: http://www.emath.fr/cocv/ Will be set by the publisher REDUCED BASIS APPROXIMATION OF LARGE SCALE PARAMETRIC ALGEBRAIC RICCATI EQUATIONS Schmidt,
More informationEmpirical Interpolation of Nonlinear Parametrized Evolution Operators
MÜNSTER Empirical Interpolation of Nonlinear Parametrized Evolution Operators Martin Drohmann, Bernard Haasdonk, Mario Ohlberger 03/12/2010 MÜNSTER 2/20 > Motivation: Reduced Basis Method RB Scenario:
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 informationReduced-Order Greedy Controllability of Finite Dimensional Linear Systems. Giulia Fabrini Laura Iapichino Stefan Volkwein
Universität Konstanz Reduced-Order Greedy Controllability of Finite Dimensional Linear Systems Giulia Fabrini Laura Iapichino Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 364, Oktober 2017 ISSN
More informationA-Posteriori Error Estimation for Parameterized Kernel-Based Systems
A-Posteriori Error Estimation for Parameterized Kernel-Based Systems Daniel Wirtz Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation University of Stuttgart 70569 Stuttgart, Pfaffenwaldring
More informationReduced Basis Methods for inverse problems
Reduced Basis Methods for inverse problems Dominik Garmatter dominik.garmatter@mathematik.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany Joint work with Bastian
More informationModel reduction for multiscale problems
Model reduction for multiscale problems Mario Ohlberger Dec. 12-16, 2011 RICAM, Linz , > Outline Motivation: Multi-Scale and Multi-Physics Problems Model Reduction: The Reduced Basis Approach A new Reduced
More informationLecture 9. Introduction to Kalman Filtering. Linear Quadratic Gaussian Control (LQG) G. Hovland 2004
MER42 Advanced Control Lecture 9 Introduction to Kalman Filtering Linear Quadratic Gaussian Control (LQG) G. Hovland 24 Announcement No tutorials on hursday mornings 8-9am I will be present in all practical
More informationA reduced basis method and ROM-based optimization for batch chromatography
MAX PLANCK INSTITUTE MedRed, 2013 Magdeburg, Dec 11-13 A reduced basis method and ROM-based optimization for batch chromatography Yongjin Zhang joint work with Peter Benner, Lihong Feng and Suzhou Li Max
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Proper Orthogonal Decomposition (POD) Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu 1 / 43 Outline 1 Time-continuous Formulation 2 Method of Snapshots
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
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 informationReduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems. Laura Iapichino Stefan Trenz Stefan Volkwein
Universität Konstanz Reduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems Laura Iapichino Stefan Trenz Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 347, Dezember 2015
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationI. D. Landau, A. Karimi: A Course on Adaptive Control Adaptive Control. Part 9: Adaptive Control with Multiple Models and Switching
I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 1 Adaptive Control Part 9: Adaptive Control with Multiple Models and Switching I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 2 Outline
More informationLecture 8. Applications
Lecture 8. Applications Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering May 8, 205 / 3 Logistics hw7 due this Wed, May 20 do an easy problem or CYOA hw8 (design problem) will
More informationOrder reduction numerical methods for the algebraic Riccati equation. V. Simoncini
Order reduction numerical methods for the algebraic Riccati equation V. Simoncini Dipartimento di Matematica Alma Mater Studiorum - Università di Bologna valeria.simoncini@unibo.it 1 The problem Find X
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More 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 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 information6.241 Dynamic Systems and Control
6.41 Dynamic Systems and Control Lecture 5: H Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 11, 011 E. Frazzoli (MIT) Lecture 5: H Synthesis May 11, 011
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationGreedy control. Martin Lazar University of Dubrovnik. Opatija, th Najman Conference. Joint work with E: Zuazua, UAM, Madrid
Greedy control Martin Lazar University of Dubrovnik Opatija, 2015 4th Najman Conference Joint work with E: Zuazua, UAM, Madrid Outline Parametric dependent systems Reduced basis methods Greedy control
More informationComparison of Feedback Controller for Link Stabilizing Units of the Laser Based Synchronization System used at the European XFEL
Comparison of Feedback Controller for Link Stabilizing Units of the Laser Based Synchronization System used at the European XFEL M. Heuer 1 G. Lichtenberg 2 S. Pfeiffer 1 H. Schlarb 1 1 Deutsches Elektronen
More informationSome a posteriori error bounds for reduced order modelling of (non-)parametrized linear systems
Max Planck Institute Magdeburg Preprints Lihong Feng Athanasios C. Antoulas Peter Benner Some a posteriori error bounds for reduced order modelling of (non-)parametrized linear systems MAX PLANCK INSTITUT
More informationKrylov Subspace Methods for Nonlinear Model Reduction
MAX PLANCK INSTITUT Conference in honour of Nancy Nichols 70th birthday Reading, 2 3 July 2012 Krylov Subspace Methods for Nonlinear Model Reduction Peter Benner and Tobias Breiten Max Planck Institute
More informationA POD Projection Method for Large-Scale Algebraic Riccati Equations
A POD Projection Method for Large-Scale Algebraic Riccati Equations Boris Kramer Department of Mathematics Virginia Tech Blacksburg, VA 24061-0123 Email: bokr@vt.edu John R. Singler Department of Mathematics
More informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationLinear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters
Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 204 Emo Todorov (UW) AMATH/CSE 579, Winter
More information5. Observer-based Controller Design
EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1
More informationReduced Basis Methods for Parametrized PDEs A Tutorial Introduction for Stationary and Instationary Problems
B. Haasdonk Reduced Basis Methods for Parametrized PDEs A Tutorial Introduction for Stationary and Instationary Problems Stuttgart, March 2014, revised version March 2016 Institute of Applied Analysis
More informationLecture 7 : Generalized Plant and LFT form Dr.-Ing. Sudchai Boonto Assistant Professor
Dr.-Ing. Sudchai Boonto Assistant Professor Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Quadratic Gaussian The state space
More informationLinear Quadratic Gausssian Control Design with Loop Transfer Recovery
Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Leonid Freidovich Department of Mathematics Michigan State University MI 48824, USA e-mail:leonid@math.msu.edu http://www.math.msu.edu/
More informationStructured State Space Realizations for SLS Distributed Controllers
Structured State Space Realizations for SLS Distributed Controllers James Anderson and Nikolai Matni Abstract In recent work the system level synthesis (SLS) paradigm has been shown to provide a truly
More informationMultiobjective PDE-constrained optimization using the reduced basis method
Multiobjective PDE-constrained optimization using the reduced basis method Laura Iapichino1, Stefan Ulbrich2, Stefan Volkwein1 1 Universita t 2 Konstanz - Fachbereich Mathematik und Statistik Technischen
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 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 informationOn Linear-Quadratic Control Theory of Implicit Difference Equations
On Linear-Quadratic Control Theory of Implicit Difference Equations Daniel Bankmann Technische Universität Berlin 10. Elgersburg Workshop 2016 February 9, 2016 D. Bankmann (TU Berlin) Control of IDEs February
More informationControl System Design
ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More 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 informationGreedy Control. Enrique Zuazua 1 2
Greedy Control Enrique Zuazua 1 2 DeustoTech - Bilbao, Basque Country, Spain Universidad Autónoma de Madrid, Spain Visiting Fellow of LJLL-UPMC, Paris enrique.zuazua@deusto.es http://enzuazua.net X ENAMA,
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationReduced Basis Method for Explicit Finite Volume Approximations of Nonlinear Conservation Laws
Proceedings of Symposia in Applied Mathematics Reduced Basis Method for Explicit Finite Volume Approximations of Nonlinear Conservation Laws Bernard Haasdonk and Mario Ohlberger Abstract. The numerical
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 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 informationNumerics and Control of PDEs Lecture 1. IFCAM IISc Bangalore
1/1 Numerics and Control of PDEs Lecture 1 IFCAM IISc Bangalore July 22 August 2, 2013 Introduction to feedback stabilization Stabilizability of F.D.S. Mythily R., Praveen C., Jean-Pierre R. 2/1 Q1. Controllability.
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 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 informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
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 informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationTowards parametric model order reduction for nonlinear PDE systems in networks
Towards parametric model order reduction for nonlinear PDE systems in networks MoRePas II 2012 Michael Hinze Martin Kunkel Ulrich Matthes Morten Vierling Andreas Steinbrecher Tatjana Stykel Fachbereich
More informationSYSTEMTEORI - KALMAN FILTER VS LQ CONTROL
SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL 1. Optimal regulator with noisy measurement Consider the following system: ẋ = Ax + Bu + w, x(0) = x 0 where w(t) is white noise with Ew(t) = 0, and x 0 is a stochastic
More informationStationary trajectories, singular Hamiltonian systems and ill-posed Interconnection
Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute
More informationL2 gains and system approximation quality 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION L2 gains and system approximation quality 1 This lecture discusses the utility
More informationTHE LOCALIZED REDUCED BASIS MULTISCALE METHOD
Proceedings of ALGORITMY 212 pp. 393 43 THE LOCALIZED REDUCED BASIS MULTISCALE METHOD FELIX ALBRECHT, BERNARD HAASDONK, SVEN KAULMANN OHLBERGER, AND MARIO Abstract. In this paper we introduce the Localized
More informationUsing Model Reduction techniques for simulating the heat transfer in electronic systems.
Using Model Reduction techniques for simulating the heat transfer in electronic systems. -Pramod Mathai, Dr. Benjamin Shapiro, University of Maryland, College Park Abstract: There is an increasing need
More informationIntroduction to Model Order Reduction
Introduction to Model Order Reduction Lecture 1: Introduction and overview Henrik Sandberg Kin Cheong Sou Automatic Control Lab, KTH ACCESS Specialized Course Graduate level Ht 2010, period 1 1 Overview
More informationH 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University Motivation Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process
More informationRobust Control 9 Design Summary & Examples
Robust Control 9 Design Summary & Examples Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/6/003 Outline he H Problem Solution via γ-iteration Robust stability via
More informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More informationChapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Observers,
More informationMTNS 06, Kyoto (July, 2006) Shinji Hara The University of Tokyo, Japan
MTNS 06, Kyoto (July, 2006) Shinji Hara The University of Tokyo, Japan Outline Motivation & Background: H2 Tracking Performance Limits: new paradigm Explicit analytical solutions with examples H2 Regulation
More informationNPTEL Online Course: Control Engineering
NPTEL Online Course: Control Engineering Ramkrishna Pasumarthy Assignment-11 : s 1. Consider a system described by state space model [ ] [ 0 1 1 x + u 5 1 2] y = [ 1 2 ] x What is the transfer function
More informationMarcus Pantoja da Silva 1 and Celso Pascoli Bottura 2. Abstract: Nonlinear systems with time-varying uncertainties
A NEW PROPOSAL FOR H NORM CHARACTERIZATION AND THE OPTIMAL H CONTROL OF NONLINEAR SSTEMS WITH TIME-VARING UNCERTAINTIES WITH KNOWN NORM BOUND AND EXOGENOUS DISTURBANCES Marcus Pantoja da Silva 1 and Celso
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 informationAn LQ R weight selection approach to the discrete generalized H 2 control problem
INT. J. CONTROL, 1998, VOL. 71, NO. 1, 93± 11 An LQ R weight selection approach to the discrete generalized H 2 control problem D. A. WILSON², M. A. NEKOUI² and G. D. HALIKIAS² It is known that a generalized
More informationComputer Problem 1: SIE Guidance, Navigation, and Control
Computer Problem 1: SIE 39 - Guidance, Navigation, and Control Roger Skjetne March 12, 23 1 Problem 1 (DSRV) We have the model: m Zẇ Z q ẇ Mẇ I y M q q + ẋ U cos θ + w sin θ ż U sin θ + w cos θ θ q Zw
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
More informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationLinear Matrix Inequalities in Robust Control. Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002
Linear Matrix Inequalities in Robust Control Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002 Objective A brief introduction to LMI techniques for Robust Control Emphasis on
More informationReduced Modeling in Data Assimilation
Reduced Modeling in Data Assimilation Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Challenges in high dimensional analysis and computation
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 5: Uncertainty and Robustness in SISO Systems [Ch.7-(8)] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 5:Uncertainty and Robustness () FEL3210 MIMO
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 informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationHankel Optimal Model Reduction 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and
More informationThe Behavioral Approach to Systems Theory
The Behavioral Approach to Systems Theory... and DAEs Patrick Kürschner July 21, 2010 1/31 Patrick Kürschner The Behavioral Approach to Systems Theory Outline 1 Introduction 2 Mathematical Models 3 Linear
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 Multiscale Lithium-Ion Battery Simulation
Model Reduction for Multiscale Lithium-Ion Battery Simulation Mario Ohlberger 1, Stephan Rave 1, and Felix Schindler 1 Applied Mathematics Münster, CMTC & Center for Nonlinear Science, University of Münster,
More informationLecture 2: Discrete-time Linear Quadratic Optimal Control
ME 33, U Berkeley, Spring 04 Xu hen Lecture : Discrete-time Linear Quadratic Optimal ontrol Big picture Example onvergence of finite-time LQ solutions Big picture previously: dynamic programming and finite-horizon
More informationSome recent developments on ROMs in computational fluid dynamics
Some recent developments on ROMs in computational fluid dynamics F. Ballarin 1, S. Ali 1, E. Delgado 2, D. Torlo 1,3, T. Chacón 2, M. Gómez 2, G. Rozza 1 1 mathlab, Mathematics Area, SISSA, Trieste, Italy
More informationEE C128 / ME C134 Feedback Control Systems
EE C128 / ME C134 Feedback Control Systems Lecture Additional Material Introduction to Model Predictive Control Maximilian Balandat Department of Electrical Engineering & Computer Science University of
More informationIntroduction to Model Order Reduction
KTH ROYAL INSTITUTE OF TECHNOLOGY Introduction to Model Order Reduction Lecture 1: Introduction and overview Henrik Sandberg, Bart Besselink, Madhu N. Belur Overview of Today s Lecture What is model (order)
More informationDiscussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough
Discussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough H.L. Trentelman 1 The geometric approach In the last
More informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
More informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
More informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More information6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski. Solutions to Problem Set 1 1. Massachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Solutions to Problem Set 1 1 Problem 1.1T Consider the
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 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 informationMathematics for Control Theory
Mathematics for Control Theory H 2 and H system norms Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials We will use: Michael Green and David Limebeer, Linear Robust
More informationBounding Stability Constants for Affinely Parameter-Dependent Operators
Bounding Stability Constants for Affinely Parameter-Dependent Operators Robert O Connor a a RWTH Aachen University, Aachen, Germany Received *****; accepted after revision +++++ Presented by Abstract In
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 14 January 2007 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More information