Gramian Based Model Reduction of Large-scale Dynamical Systems
|
|
- Gladys Muriel Cunningham
- 6 years ago
- Views:
Transcription
1 Gramian Based Model Reduction of Large-scale Dynamical Systems Paul Van Dooren Université catholique de Louvain Center for Systems Engineering and Applied Mechanics B-1348 Louvain-la-Neuve, Belgium Origin and motivation Ex: Storm surge forecasting Typical techniques (Gramians) Linear time-invariant Linear time-varying Non-linear Numerical/Algorithmic issues (Krylov) 1
2 Storm surge forecasting in the North Sea see Verlaan-Heemink 97 depth Wick n North Shields Lowestoft Sheerness Dover Huibertgat Harlingen Delfzijl Den Helder IJmuiden-II HoekvanHolland Vlissingen m Problem Using measurements predict the state of the North Sea variables in order to operate the sluices in due time (6h.) Solution x(t). = [h(t), v x (t), v y (t)] satisfies the shallow water equations x(t)/ t = F (x(t), w(t)) y(t) = G(x(t), v(t)) with measurements y(t) and noise processes v(.), w(.) = estimate and predict ˆx(t) using Kalman filtering 2
3 Some data : very few measurements (x s and + s) water level current HF radar station Scheveningen Noordwijk Hook of Holland but discretized state is very large-scale (6. variables) y [km] x [km] 3
4 Nevertheless it works... True water levels Estimated water levels ZUNOWAK True water levels (t=1) ZUNOWAK Estimated water levels (t=1) GERMANY 4 GERMANY 2 ENGLAND THE NETHERLANDS 2 ENGLAND THE NETHERLANDS BELGIUM FRANCE RIJKSWATERSTAAT RIKZ BELGIUM FRANCE RIJKSWATERSTAAT RIKZ ZUNOWAK True water levels (t=3) ZUNOWAK Estimated water levels (t=3) GERMANY GERMANY 2 ENGLAND THE NETHERLANDS 2 ENGLAND THE NETHERLANDS 1.4 BELGIUM FRANCE RIJKSWATERSTAAT RIKZ BELGIUM FRANCE RIJKSWATERSTAAT RIKZ ZUNOWAK True water levels (t=8) ZUNOWAK Estimated water levels (t=8) GERMANY GERMANY 2.2 ENGLAND.4.6 THE NETHERLANDS 2 ENGLAND THE NETHERLANDS BELGIUM FRANCE RIJKSWATERSTAAT RIKZ BELGIUM FRANCE RIJKSWATERSTAAT RIKZ Reconstruction works well around estuarium 4
5 Visualisation of computed variance of the error std water level [m] n m Standard deviation of filter using 8 measurement locations std water level [m] n m Standard deviation of filter using measurement locations (1 used for validation) 5
6 Suppose a good discretization x(.) R N is given Dynamical systems modeled via explicit equations y 1 (.) y 2 (.). y p (.) x(t) R N, N >> m, p u 1 (.) u 2 (.). u m (.) discretize = continuous-time { ẋ(t) = G(x(t), u(t)) y(t) = H(x(t), u(t)) discrete-time { x(k + 1) = G(x(k), u(k)) y(k) = H(x(k), u(k)) linearize { ẋ(t)=a(t)x(t)+b(t)u(t) y(t) = C(t)x(t)+D(t)u(t) { x(k+1) = A(k)x(k)+B(k)u(k) y(k) = C(k)x(k)+D(k)u(k) { ẋ(t) = Ax(t) + Bu(t) y(t) = Cx(t) + Du(t) freeze time { x(k + 1) = Ax(k) + Bu(k) y(k) = Cx(k) + Du(k) Many control problems require N 3 /( t) operations. Because of cubic complexity in N model reduction 6
7 Linear Time Invariant Systems Given large model {A NN, B Nm, C pn } { ẋ(t) = Ax(t) + Bu(t) y(t) = Cx(t) + Du(t), u(t) R m, y(t) R p, x(t) R N, N << m, p find small model {Ânn, ˆB nm, Ĉpn } with n << N { ˆx(t) = ˆx(t) + ˆBu(t) ŷ(t) = Ĉ ˆx(t) + ˆDu(t), driven by the same input u(t) with small error y(t) ŷ(t) Model reduction = find a smaller model, i.e. n << N : approximation problem stability is important measure is important 7
8 How to capture the essence of the system? Transfer functions and norms H(s) = C(sI N A) 1 B + D, Ĥ(s) = Ĉ(sI n Â) 1 ˆB + ˆD, are p m rational matrices try to match frequency responses Frequency response Frequency response Full order model Reduced order model Theory : by minimizing their difference using H(.) Ĥ(.) balanced truncation (Moore 81). = sup σ max {H(jω) Ĥ(jω)} ω optimal Hankel norm approximation (Adamian-Arov-Krein 71, Glover 9) interpolation (Gragg-Lindquist 83) Other references : Gallivan-Grimme-VanDooren, Jaimoukha-Kasanelly, Villemagne-Skelton, Boley, Craig, Freund-Feldman, Sorensen-Antoulas,... 8
9 Why. norm? Fourier transforms of signals : u f (ω) = Fu(t), y f (ω) = Fy(t), ŷ f (ω) = Fŷ(t) yields y f (ω) = H(jω)u f (ω), ŷ f (ω) = Ĥ(jω)u f(ω). and hence a bound for e(t). = [y(t) ŷ(t)] : Fe(t) = e f (ω) = [H(jω) Ĥ(jω)]u f(ω). Minimize worst case error e f (ω) 2 for u f (ω) 2 = 1 by minimizing H(.) Ĥ(.). = sup H(jω) Ĥ(jω) 2, ω but this is a difficult norm to handle! 9
10 Use Hankel norm instead of. approximations Consider the mapping past inputs = future outputs Continuous-time y(t) = Ce A(t τ) Bu(τ)dτ = Ce At y(t) = Ce At x(), x() = e Aτ Bu( τ)dτ, e Aτ Bu( τ)dτ Continuous input u(t) Continuous output y(t) Discrete-time y(k) = CA (k j) Bu(j) = CA k A j Bu( j), y(k) = CA k x(), x() = A j Bu( j) Discrete input u(k) Discrete output y(k) 1
11 N N Gramians derived from the Hankel map From y([, )) = Ox(), x() = Cu((, )), define the dual maps O : y([, )) x(), C : x() u((, )) and the (observability and controllability) Gramians G o. = O O, G c. = CC Continuous-time G o = + (Ce At ) T (Ce At )dt, G c = + (e At B)(e At B) T dt, Discrete-time G o = + (CA k ) T (CA k ), G c = + (A k B)(A k B) T, Gramians can be viewed as energy functions G c for past inputs x() G o for x() future outputs Perform ordered eigendecomposition Λ = T 1 (G c G o )T (with λ n >> λ n+1 ) and project on the first n coordinates : T 1 AT =. [ ] A11 A 12, T 1 B. [ ] B1 =, CT =. [ ] C A 21 A 22 B 1 C 2. 2 {Â, ˆB, Ĉ}. = {A 11, B 1, C 1 } 11
12 Interpolating the frequency response H(ω) seems a good idea since Parseval s theorem implies and G o = 1 2π G c = 1 2π + G o = 1 2π G c = 1 2π ( jωi A T ) 1 C T C(jωI A) 1 dω, (jωi A) 1 BB T ( jωi A T ) 1 dω. (e jω I A T ) 1 C T C(e jω I A) 1 (e jω I A) 1 BB T (e jω I A T ) 1. What technique to use? The discrete-time case : y u 1 C y 1 y 2 = OC u 2 u 3 = CA [ CA 2 B AB A 2 B... ] suggests to use Krylov spaces! K j (M, R) = Im { R, MR, M 2 R,..., M j 1 R } u 1 u 2 u 3. 12
13 Rational interpolation and moment matching Let X and Y define a projector (Y T X = I n ) and {Â, ˆB, Ĉ, D} = { Y T AX, Y T B, CX, D } Taylor series of H(s). = C(sI A) 1 B + D around H(s) = H + H 1 s 1 + H 2 s 2 +, where the moments H i are equal to : H = D, H i = CA i 1 B, i = 1, 2,... The reduced order model Ĥ(s). = Ĉ(sI Â) 1 ˆB + ˆD has a similar expansion with moments Ĥi : Ĥ(s) = Ĥ + Ĥ1s 1 + Ĥ2s 2 +, Ĥ = ˆD, Ĥ i = ĈÂi 1 ˆB, i = 1, 2,... Moment matching of both models (Padé approximation) is obtained by using Krylov spaces (Gragg-Lindquist 83) 13
14 Theorem : Let m = p, Y T X = I n and assume ImX = Im [ B, AB, A 2 B,... A k 1 B ], ] ImY = Im [C T, A T C T, A 2T C T,... A (k 1)T C T then the first 2k moments match : H j = Ĥj, j = 1,..., 2k. Rational Krylov methods extend this to several points σ i : multipoint (Padé) approximations are obtained by just using modified Krylov spaces : (Gallivan-Grimme-VanDooren 97) ImX = i K ji {(A σ i I) 1, (A σ i I) 1 B} ImY = i K ji {(A T σ i I) 1, (A T σ i I) 1 C T } in the above theorem 14
15 Comparison Optimal and rational approximations 15th order approximation of 12 th order CD player Frequency responses in db 5 1 Frequency responses in db th order Hankel norm approximation th order balanced truncation 1 5 Frequency responses in db th order RK with points 1,2,2,2,2,2,2,2,2,3,4,4,5,5,Inf Legend : Hankel norm - - Balanced truncation Rational Krylov Errors T (.) ˆT (.) ln( T (.) ˆT (.) ) Hankel Balanc Rat. Kr Rational approximation looks better on a logarithmic scale 15
16 Time-varying systems Approximate the (discrete) time-varying systems { x(k + 1) = A(k)x(k) + B(k)u(k) y(k) = C(k)x(k) + D(k)u(k) by a lower order models of same type. We notice that y(k) C(k) y(k+1) y(k+2) = C(k+1)A(k) C(k+2)A(k+1)A(k) x(k),.. x(k) = [ B(k 1) A (k 1) B (k 2) A (k 1) A (k 2) B (k 3)... ] u(k 1) u(k 2) u(k 3). which again suggests Krylov. Use low rank approximations G o (k) S o (k)s T o (k), G c (k) S c (k)s T c (k), where S o (k) and S c (k) are N n matrices Such approximations are obtained e.g. by keeping only the n dominant singular vectors at step k of the Krylov recurrence : S c (k) = SV D n [B(k), A(k)S c (k 1)] 16
17 Let s go back to the Storm Surge example The important matrix in this Kalman filtering problem is P (k) = E{(x(k) ˆx(k))(x(k) ˆx(k)) T } which represents the error covariance of the estimation. If we approximate P (k) S(k)S(k) T, S(k) R N n, then we obtain the recurrence S(k + 1) SV D n [ R(k) C(k)S(k) A(k)S(k) B(k)Q(k) where the new factor S(k + 1) stays of rank n by a projection (using the SVD). The only big matrix involved here is the sparse matrix A(k) which is multiplied with only n columns Verlaan-Heemink 97 ] 17
18 Time-varying linearized problems Consider the discrete-time system { x(k + 1) = G(x(k), u(k)) y(k) = H(x(k), u(k)). One could linearize along a nominal trajectory (x(k), u(k)) and get A(.), B(.), C(.), D(.) from Taylor expansion of G(.,.), H(.,.) Simpler idea (POD) : (Holmes-Lumley-Berkooz 96) Use the energy function G. = k f k=k i x(k)x(k) T. From x(k + 1) = A(k)x(k) with initial conditions x(k i ), we have x(k) = Φ(k, k i )x(k i ). Therefore G looks like a Gramian : k f (Φ(k, k i )x(k i ))(Φ(k, k i )x(k i )) T. k=k i Now project on its dominant subspace (POD) 18
19 Example : Use POD in CVD reactor (Ly-Tran 99) Schematic representation of a horizontal quartz reactor in a steel confinement shell Compute state trajectories for one typical input : Snap shots of typical states Ten dominant states 19
20 Concluding remarks large-scale is typically sparse time-stepping (simulation) is cheaper than control (optimization) find an energy function that is cheap and project on its dominant features Future work find error bounds on the fly incorporate projections in closed loop Further reading P. Van Dooren, Gramian based model reduction of large-scale dynamical systems, in Numerical Analysis 1999, Chapman and Hall, pp , CRC Press, London, 2. M. Verlaan and A. Heemink, Tidal flow forecasting using reduced rank square root filters, Stochastic Hydrology and Hydraulics, 11, pp , See also SIAM short course notes on vdooren/ 2
Gramian Based Model Reduction of Large-scale Dynamical Systems
Gramian Based Model Reduction of Large-scale Dynamical Systems Paul Van Dooren Université catholique de Louvain Center for Systems Engineering and Applied Mechanics B-1348 Louvain-la-Neuve, Belgium vdooren@csamuclacbe
More informationProjection Methods for Model Reduction of Large-scale Systems
Projection Methods for Model Reduction of Large-scale Systems Kyle A. Gallivan School of Computational Science and Information Technology Florida State University Support: NSF, DOE, IBM Corporation, Belgian
More informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview Thanos Antoulas joint work with Dan Sorensen aca@rice.edu URL: www-ece.rice.edu/ aca CITI Lecture, Rice, February 2005 Model reduction of large-scale
More informationFourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics
Fourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics K. Willcox and A. Megretski A new method, Fourier model reduction (FMR), for obtaining stable, accurate, low-order models
More 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 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 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 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 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 informationModule 08 Observability and State Estimator Design of Dynamical LTI Systems
Module 08 Observability and State Estimator Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha November
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 information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 7: State-space Models Readings: DDV, Chapters 7,8 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology February 25, 2011 E. Frazzoli
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 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 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 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 informationProjection of state space realizations
Chapter 1 Projection of state space realizations Antoine Vandendorpe and Paul Van Dooren Department of Mathematical Engineering Université catholique de Louvain B-1348 Louvain-la-Neuve Belgium 1.0.1 Description
More 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 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 informationStability preserving post-processing methods applied in the Loewner framework
Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 1 / 20 Stability preserving post-processing methods applied in the Loewner framework Ion Victor
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 informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
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 information16.30 Estimation and Control of Aerospace Systems
16.30 Estimation and Control of Aerospace Systems Topic 5 addendum: Signals and Systems Aeronautics and Astronautics Massachusetts Institute of Technology Fall 2010 (MIT) Topic 5 addendum: Signals, Systems
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More informationSystem Identification by Nuclear Norm Minimization
Dept. of Information Engineering University of Pisa (Italy) System Identification by Nuclear Norm Minimization eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain Systems
More informationModule 07 Controllability and Controller Design of Dynamical LTI Systems
Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October
More informationModel Reduction of Linear Systems, an Interpolation point of View
Université Catholique De Louvain Faculté des sciences appliquées Département d ingénierie mathématique Center for Systems Engineering and Applied Mechanics Model Reduction of Linear Systems, an Interpolation
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 informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationObservability and state estimation
EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time observability
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More informationLecture 19 Observability and state estimation
EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time
More informationComparison of four state observer design algorithms for MIMO system
Archives of Control Sciences Volume 23(LIX), 2013 No. 2, pages 131 144 Comparison of four state observer design algorithms for MIMO system VINODH KUMAR. E, JOVITHA JEROME and S. AYYAPPAN A state observer
More informationBALANCING AS A MOMENT MATCHING PROBLEM
BALANCING AS A MOMENT MATCHING PROBLEM T.C. IONESCU, J.M.A. SCHERPEN, O.V. IFTIME, AND A. ASTOLFI Abstract. In this paper, we treat a time-domain moment matching problem associated to balanced truncation.
More informationProblem set 5 solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION Problem set 5 solutions Problem 5. For each of the stetements below, state
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationPrinciple Component Analysis and Model Reduction for Dynamical Systems
Principle Component Analysis and Model Reduction for Dynamical Systems D.C. Sorensen Virginia Tech 12 Nov 2004 Protein Substate Modeling and Identification PCA Dimension Reduction Using the SVD Tod Romo
More 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 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 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 informationSubspace-based Identification
of Infinite-dimensional Multivariable Systems from Frequency-response Data Department of Electrical and Electronics Engineering Anadolu University, Eskişehir, Turkey October 12, 2008 Outline 1 2 3 4 Noise-free
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 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 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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
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 informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
More informationModel Reduction (Approximation) of Large-Scale Systems. Introduction, motivating examples and problem formulation Lecture 1
Lecture outines Motivating examples Norm and linear algebra reminder Introduction Summary Model Reduction (Approximation) of Large-Scale Systems Introduction, motivating examples and problem formulation
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
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 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 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 informationPost flood analysis. Demands of the Storm Surge Warning Service (SVSD) of the Netherlands. Annette Zijderveld/ Martin Verlaan
Post flood analysis Demands of the Storm Surge Warning Service (SVSD) of the Netherlands Annette Zijderveld/ Martin Verlaan Objective Bringing in coastal issues to the FEWS user days (second attempt) Demonstrating
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 informationBalancing of Lossless and Passive Systems
Balancing of Lossless and Passive Systems Arjan van der Schaft Abstract Different balancing techniques are applied to lossless nonlinear systems, with open-loop balancing applied to their scattering representation.
More 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 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 informationComputation and Application of Balanced Model Order Reduction
Computation and Application of Balanced Model Order Reduction D.C. Sorensen Collaborators: M. Heinkenschloss, A.C. Antoulas and K. Sun Support: NSF and AFOSR MIT Nov 2007 Projection Methods : Large Scale
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 informationAPPROXIMATE SOLUTION OF A SYSTEM OF LINEAR EQUATIONS WITH RANDOM PERTURBATIONS
APPROXIMATE SOLUTION OF A SYSTEM OF LINEAR EQUATIONS WITH RANDOM PERTURBATIONS P. Date paresh.date@brunel.ac.uk Center for Analysis of Risk and Optimisation Modelling Applications, Department of Mathematical
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 informationLinear System Theory. Wonhee Kim Lecture 1. March 7, 2018
Linear System Theory Wonhee Kim Lecture 1 March 7, 2018 1 / 22 Overview Course Information Prerequisites Course Outline What is Control Engineering? Examples of Control Systems Structure of Control Systems
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 informationRational 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 informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More 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 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 informationDATA ASSIMILATION STUDY USING THE DUTCH CONTINENTAL SHELF MODEL WITH FULL MEASUREMENT
DATA ASSIMILATION STUDY USING THE DUTCH CONTINENTAL SHELF MODEL WITH FULL MEASUREMENT JULIUS H. SUMIHAR AND MARTIN VERLAAN,2 Faculty of Electrical Engineering, Mathematics and Computer Science, Department
More informationModel reduction of large-scale systems
Model reduction of large-scale systems Lecture IV: Model reduction from measurements Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/ aca International School,
More informationStabilisation of Linear Time-Invariant Systems subject to Output Saturation
Stabilisation of Linear Time-Invariant Systems subject to Output Saturation by Gijs Hilhorst, BSc A Thesis in Partial Fulfillment of the Requirement for the Degree of Master of Science at the Department
More informationDiscrete-time linear systems
Automatic Control Discrete-time linear systems Prof. Alberto Bemporad University of Trento Academic year 2-2 Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year 2-2 / 34 Introduction
More informationBalancing-Related Model Reduction for Large-Scale Systems
Balancing-Related Model Reduction for Large-Scale Systems Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz D-09107 Chemnitz benner@mathematik.tu-chemnitz.de
More 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 informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture State space models, 1st part: Model: Sec. 10.1 The
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 informationDESIGNING A KALMAN FILTER WHEN NO NOISE COVARIANCE INFORMATION IS AVAILABLE. Robert Bos,1 Xavier Bombois Paul M. J. Van den Hof
DESIGNING A KALMAN FILTER WHEN NO NOISE COVARIANCE INFORMATION IS AVAILABLE Robert Bos,1 Xavier Bombois Paul M. J. Van den Hof Delft Center for Systems and Control, Delft University of Technology, Mekelweg
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 information14 - Gaussian Stochastic Processes
14-1 Gaussian Stochastic Processes S. Lall, Stanford 211.2.24.1 14 - Gaussian Stochastic Processes Linear systems driven by IID noise Evolution of mean and covariance Example: mass-spring system Steady-state
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More informationSystems and Control Theory Lecture Notes. Laura Giarré
Systems and Control Theory Lecture Notes Laura Giarré L. Giarré 2017-2018 Lesson 5: State Space Systems State Dimension Infinite-Dimensional systems State-space model (nonlinear) LTI State Space model
More informationLinear dynamical systems with inputs & outputs
EE263 Autumn 215 S. Boyd and S. Lall Linear dynamical systems with inputs & outputs inputs & outputs: interpretations transfer function impulse and step responses examples 1 Inputs & outputs recall continuous-time
More informationSparsity in system identification and data-driven control
1 / 40 Sparsity in system identification and data-driven control Ivan Markovsky This signal is not sparse in the "time domain" 2 / 40 But it is sparse in the "frequency domain" (it is weighted sum of six
More informationA Novel Scheme for Positive Real Balanced Truncation
A Novel Scheme for Positive Real Balanced runcation Kari Unneland, Paul Van Dooren and Olav Egeland Abstract Here, model reduction, based on balanced truncation, of stable and passive systems will be considered.
More 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 informationROBUST PASSIVE OBSERVER-BASED CONTROL FOR A CLASS OF SINGULAR SYSTEMS
INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 5 Number 3-4 Pages 480 487 c 2009 Institute for Scientific Computing and Information ROBUST PASSIVE OBSERVER-BASED CONTROL FOR A CLASS OF
More informationSolving Dynamic Equations: The State Transition Matrix
Overview Solving Dynamic Equations: The State Transition Matrix EGR 326 February 24, 2017 Solutions to coupled dynamic equations Solutions to dynamic circuits from EGR 220 The state transition matrix Discrete
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 informationAdvanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification
Advanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification 1. Consider the time series x(k) = β 1 + β 2 k + w(k) where β 1 and β 2 are known constants
More informationKalman Filter and Parameter Identification. Florian Herzog
Kalman Filter and Parameter Identification Florian Herzog 2013 Continuous-time Kalman Filter In this chapter, we shall use stochastic processes with independent increments w 1 (.) and w 2 (.) at the input
More informationData Driven Discrete Time Modeling of Continuous Time Nonlinear Systems. Problems, Challenges, Success Stories. Johan Schoukens
1/51 Data Driven Discrete Time Modeling of Continuous Time Nonlinear Systems Problems, Challenges, Success Stories Johan Schoukens fyuq (,, ) 4/51 System Identification Data Distance Model 5/51 System
More informationWhy do we need a model?
Why do we need a model? For simulation (study the system output for a given input) Ex. Thermal study of a space shuttle when it enters the atmosphere (étude thermique concernant la rentrée dans l atmosphère
More informationRural/Urban Migration: The Dynamics of Eigenvectors
* Analysis of the Dynamic Structure of a System * Rural/Urban Migration: The Dynamics of Eigenvectors EGR 326 April 11, 2019 1. Develop the system model and create the Matlab/Simulink model 2. Plot and
More informationLecture 2. Linear Systems
Lecture 2. Linear Systems Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering April 6, 2015 1 / 31 Logistics hw1 due this Wed, Apr 8 hw due every Wed in class, or my mailbox on
More informationS. Gugercin and A.C. Antoulas Department of Electrical and Computer Engineering Rice University
Proceedings of the 39" IEEE Conference on Decision and Control Sydney, Australia December, 2000 A Comparative Study of 7 Algorithms for Model Reduct ion' S. Gugercin and A.C. Antoulas Department of Electrical
More informationIntro. Computer Control Systems: F8
Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2
More informationSquare-Root Algorithms of Recursive Least-Squares Wiener Estimators in Linear Discrete-Time Stochastic Systems
Proceedings of the 17th World Congress The International Federation of Automatic Control Square-Root Algorithms of Recursive Least-Squares Wiener Estimators in Linear Discrete-Time Stochastic Systems Seiichi
More informationOrder Reduction for Large Scale Finite Element Models: a Systems Perspective
Order Reduction for Large Scale Finite Element Models: a Systems Perspective William Gressick, John T. Wen, Jacob Fish ABSTRACT Large scale finite element models are routinely used in design and optimization
More information