Model Order Reduction for Parameter-Varying Systems
|
|
- Joy Griffin
- 5 years ago
- Views:
Transcription
1 Model Order Reduction for Parameter-Varying Systems Xingang Cao Promotors: Wil Schilders Siep Weiland Supervisor: Joseph Maubach Eindhoven University of Technology CASA Day November 1, 216 Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
2 Introduction Overview 1 Introduction 2 Model Order Reduction 3 Model Order Reduction for Linear Parameter-Varying Systems 4 Conclusions and Future Work Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
3 Introduction Self Introduction Russia Biography: 199 born in Baotou Mongolia BE in Electrical Engineering Tongji University, Shanghai MSc. In Systems and Control TU Eindhoven Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
4 Introduction Modeling & Control of Thermal Effects in Printing Systems Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
5 Introduction Model Order Reduction in a Broad Field Systems and control, Lyapunov, Truncated Balanced Realization Scientific computing, numerical MOR, Krylov methods, linear algebra, tensors Pade-via-Lanczos and PRIMA Passivity preserving Structure preserving Linear and nonlinear problems Parameterized methods Tensor analysis Large-scale Lyapunov systems Balanced realizations Observability and controllability Gramians Port-Hamiltonian systems Hankel singular values Projection methods Mathematical modeling, behavioral modeling, models via data MOR at operator level Reduced basis methods Karhunen-Loeve expansions Neural networks Vector fitting Behavioral models Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
6 Model Order Reduction Overview 1 Introduction 2 Model Order Reduction 3 Model Order Reduction for Linear Parameter-Varying Systems 4 Conclusions and Future Work Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
7 Model Order Reduction What is Model Order Reduction? Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
8 Model Order Reduction Abstract Problem Formulation Given a dynamical system Σ : {ẋ = f (x, u) y = g(x, u) with x(t) R n, u(t) R m and y(t) R p. Find another dynamical system with x(t) R r and Σ : { x = f ( x, u) ŷ = ĝ( x, u) r n. Constraints: Approximation error small, i.e., Σ Σ ɛ; Preserve characteristics, e.g., stability, dissipativity, etc; Computational efficient for simulation and control design. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
9 Model Order Reduction Why Model Order Reduction? Realistic models have high complexity (n ), but we want to perform simulations fast and with reliable outcomes; enable model based control design; cope with limited memory storage capacity; perform online tasks (optimization, prediction) efficiently; enable implementation on embedded systems. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
10 Model Order Reduction How to do Model Order Reduction? Model Reduction by Projection: Signal projection: V : R n R r x V x Residual: R = V x f (V x, u) Residual projection: W R = W V x W f (V x, u) z y -.1 x Trivial projection: 3D to 2D. Galerkin projection: W = V & V V = I r. Petrov Galerkin projection: W V & W V = I r. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
11 Model Order Reduction How to Find the Projections? Consider a linear time-invariant (LTI) system { Eẋ = Ax + Bu Σ : y = Cx + Du Take the Laplace transformation L : ẋ(t) sx (s), u(t) U(s) and y(t) Y (s) Σ defines an input-output mapping G(s) = C(sE A) 1 B + D. Expand G(s) = G(s + σ) around a point s, G(σ) = C (s E A) 1 E (s }{{} E A) 1 B }{{} i= M B σ i Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
12 Model Order Reduction How to Find Projections (Cont d)? Theorem If W and V span the left and right principal subspaces of the approximation error is small. M = (s E A) 1 E, Compute W and V by Krylov subspace methods W := span{w } = span{c, M C,..., (M ) q 1 C } V := span{v } = span{ B, M B,..., M q 2 B} dim(w) = dim(v) = r. Remark For B, C R n 1 (single-input single-output system), r = q 1 = q 2. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
13 Model Order Reduction for Linear Parameter-Varying Systems Overview 1 Introduction 2 Model Order Reduction 3 Model Order Reduction for Linear Parameter-Varying Systems 4 Conclusions and Future Work Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
14 Model Order Reduction for Linear Parameter-Varying Systems Linear Parameter-Varying Systems Linear parameter-varying (LPV) system is { E(p)ẋ = A(p)x + B(p)u Σ(p) : y = C(p)x + D(p)u where p(t) : R + P Rµ. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
15 Model Order Reduction for Linear Parameter-Varying Systems Track the Parameterized Projection Subspaces Problem Track/Estimate the principal subspaces of M(p) = (s E(p) A(p)) 1 E(p) with starting vectors C(p) and B(p). Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
16 Model Order Reduction for Linear Parameter-Varying Systems Geodesics on Grassmann Manifold Definition Let Y(t) : I S be a smooth, injective curve and Ẏ(t) for all t R. There exists a smooth vector field X : S T S such that Ẏ(t) = X (Y(t)) for all t I. If ẎẎ := ( X X )(Y(t)) =, then Y(t) is a geodesic. On Grassmann manifold, it satisfies Ÿ(t) + Y(t)(Ẏ(t) Ẏ(t)) =, t [, 1]. Given Y() = S and Ẏ() = H, Y(t) = ( ) ( ) cos(σ S Z U t) Z, sin(σ t) S 1 : = Y(1). H thin svd == U Σ V, t [, 1]; Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
17 Model Order Reduction for Linear Parameter-Varying Systems Data-Driven Subspace Tracking Step 1. Sample parameter p and get {p i } N i=1 ; Step 2. Compute the projections {W } N i=1 and {V }N i=1 ; Step 3. Define ψ i = vec(w i Wi ) and φ i = vec(v i Vi ); Step 4. Assume that ψ i and φ i satisfy Σ ψ : ψ i+1 = A ψ ψ i + B ψ p i+1, Σ φ : φ i+1 = A φ ψ i + B φ p i+1 ; Data Prediction φ 2 φ 1 φ A φ A φ φ 1 φ 2 B φ B φ p 1 p p 1 p 2 Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
18 Model Order Reduction for Linear Parameter-Varying Systems A Small-scale Example Example Consider the subspaces 1 W(p(t)) = span{, 1 V(p(t)) = span{, cos(p(t)) }, sin(p(t)) cos(p(t) + π 32 ) }. sin(p(t) + π 32 ) where e i, i = 1, 2, 3 R 3 and p(t) [ 1 1, π 4 π 32 ]. Take 2 random examples and apply the data-driven method to obtain a 3rd-order model. When the parameter is varying linearly from 1 1 to, compare the approximated and the real subspace representations. π 4 π 32 Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
19 Model Order Reduction for Linear Parameter-Varying Systems Approximation of W (p) W 11 = 1 Wreal Wapprox 1.5 W 12 = Wreal Wapprox 1.5 W 13 = Wreal Wapprox p W 21 = 1 Wreal.5 Wapprox p p W 22 = cos(p) Wreal Wapprox p p W 23 = sin(p) Wreal Wapprox p Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
20 Model Order Reduction for Linear Parameter-Varying Systems Approximation of V (p) V 11 = 1 Vreal Vapprox 1.5 V 12 = Vreal Vapprox 1.5 V 13 = Vreal Vapprox p V 21 = 1 Vreal.5 Vapprox p p V 22 = cos(p+π/32) 1 Vreal Vapprox p p V 23 = sin(p+π/32) Vreal Vapprox p Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
21 Conclusions and Future Work Overview 1 Introduction 2 Model Order Reduction 3 Model Order Reduction for Linear Parameter-Varying Systems 4 Conclusions and Future Work Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
22 Conclusions and Future Work Conclusions 1 Data-driven method: allows to track the parameter-varying subspaces; 2 Large-scale problems: off-line cost for data-driven method is high; 3 Data-driven method: faster at online phase than direct method. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
23 Conclusions and Future Work Future Work Formulate the LPV model order reduction problem as {W (p(t)), V (p(t))} = argmax Trace(W M(p(t))V ) W,V Gr(r, R n ) st. W V = I or W W = I, V V = I, W and V starts with specific vector B and C. where Gr(r, R n ) := {span(a) A R n r, rank(a) = r}. Solve it by optimization on matrix manifold methods. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
24 Conclusions and Future Work Questions? Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
25 Conclusions and Future Work References Schilders, Van der Vorst, Rommes. Model Order Reduction: Theory, Research Aspects and Applications. Springer, 28. Absil, Mahony, Sepulchre. Optimization Algorithms on Matrix Manifolds. Princeton University Press, 29. Benner, Gugercin, Willcox. A survey of projection-based model reduction methods for parametric dynamical systems. SIAM review, 215. Proctor, Brunton, Kutz. Dynamic Mode Decomposition with Control. SIAM Journal on Applied Dynamical Systems, 216. Xingang Cao (CASA TU/e) CASA Day 216 November November 1, / 25
Parametric 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 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 informationLinear and Nonlinear Matrix Equations Arising in Model Reduction
International Conference on Numerical Linear Algebra and its Applications Guwahati, January 15 18, 2013 Linear and Nonlinear Matrix Equations Arising in Model Reduction Peter Benner Tobias Breiten Max
More 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 informationFast Frequency Response Analysis using Model Order Reduction
Fast Frequency Response Analysis using Model Order Reduction Peter Benner EU MORNET Exploratory Workshop Applications of Model Order Reduction Methods in Industrial Research and Development Luxembourg,
More informationModel Reduction for Linear Dynamical Systems
Summer School on Numerical Linear Algebra for Dynamical and High-Dimensional Problems Trogir, October 10 15, 2011 Model Reduction for Linear Dynamical Systems Peter Benner Max Planck Institute for Dynamics
More informationModel Reduction 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 informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
More informationModel Reduction of 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 informationRecent developments for MOR in the electronics industry
Recent developments for MOR in the electronics industry Wil Schilders Reduced Order Models in Computational Science and Engineering Aachen, January 30-31, 2014 Full proposal oc-2013-1-15312 for a new COST
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 informationAre thermodynamical systems port-hamiltonian?
Are thermodynamical systems port-hamiltonian? Siep Weiland Department of Electrical Engineering Eindhoven University of Technology February 13, 2019 Siep Weiland (TUE) Are thermodynamical systems port-hamiltonian?
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 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 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 informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
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 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 informationA Black-Box Method for Parametric Model Order Reduction based on Matrix Interpolation with Application to Simulation and Control
A Black-Box Method for Parametric Model Order Reduction based on Matrix Interpolation with Application to Simulation and Control Matthias Martin Georg Geuß Vollständiger Abdruck der von der Fakultät für
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 informationPassive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form
Passive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form Boyuan Yan, Sheldon X.-D. Tan,PuLiu and Bruce McGaughy Department of Electrical Engineering, University of
More 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 informationLinear Matrix Inequality (LMI)
Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F 1 + + x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the
More informationWeighted balanced realization and model reduction for nonlinear systems
Weighted balanced realization and model reduction for nonlinear systems Daisuke Tsubakino and Kenji Fujimoto Abstract In this paper a weighted balanced realization and model reduction for nonlinear systems
More informationarxiv: v1 [math.ds] 11 Jul 2016
Model reduction of linear time-varying systems with applications for moving loads Maria Cruz Varona and Boris Lohmann Preprint: June 15, 2016 arxiv:1607.02846v1 [math.ds] 11 Jul 2016 Abstract In this paper
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 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 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 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 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 informationIFT 6760A - Lecture 1 Linear Algebra Refresher
IFT 6760A - Lecture 1 Linear Algebra Refresher Scribe(s): Tianyu Li Instructor: Guillaume Rabusseau 1 Summary In the previous lecture we have introduced some applications of linear algebra in machine learning,
More informationCME 345: MODEL REDUCTION - Projection-Based Model Order Reduction
CME 345: MODEL REDUCTION - Projection-Based Model Order Reduction Projection-Based Model Order Reduction Charbel Farhat and David Amsallem Stanford University cfarhat@stanford.edu 1 / 38 Outline 1 Solution
More informationFocus was on solving matrix inversion problems Now we look at other properties of matrices Useful when A represents a transformations.
Previously Focus was on solving matrix inversion problems Now we look at other properties of matrices Useful when A represents a transformations y = Ax Or A simply represents data Notion of eigenvectors,
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 informationLinear System Theory
Linear System Theory Wonhee Kim Lecture 3 Mar. 21, 2017 1 / 38 Overview Recap Nonlinear systems: existence and uniqueness of a solution of differential equations Preliminaries Fields and Vector Spaces
More informationEmpirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems
Empirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems Antoni Ras Departament de Matemàtica Aplicada 4 Universitat Politècnica de Catalunya Lecture goals To review the basic
More informationNonlinear Control Lecture 7: Passivity
Nonlinear Control Lecture 7: Passivity Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 7 1/26 Passivity
More informationTwo-sided Eigenvalue Algorithms for Modal Approximation
Two-sided Eigenvalue Algorithms for Modal Approximation Master s thesis submitted to Faculty of Mathematics at Chemnitz University of Technology presented by: Supervisor: Advisor: B.sc. Patrick Kürschner
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 informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Introduction to Nonlinear Controllability and Observability Hanz Richter Mechanical Engineering Department Cleveland State University Definition of
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationSpline Interpolation Based PMOR
Universitaet Bremen Zentrum fuer Technomathematik December 14, 2009 Outline 1 Statement of the problem: Parametric Model Order Reduction Interpolation with B-Splines 2 3 Statement of the problem Statement
More informationLINEAR ALGEBRA W W L CHEN
LINEAR ALGEBRA W W L CHEN c W W L Chen, 994, 28. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,
More informationA comparison of model reduction techniques from structural dynamics, numerical mathematics and systems and control
A comparison of model reduction techniques from structural dynamics, numerical mathematics and systems and control B. Besselink a, A. Lutowska b, U. Tabak c, N. van de Wouw a, H. Nijmeijer a, M.E. Hochstenbach
More informationLinear systems the finite dimensional case
Linear systems the finite dimensional case Johann Baumeister Goethe University, Frankfurt, Germany Rio de Janeiro / October 2017 Outline Controlled and observed systems Stability Observability Controllablity
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 informationEXAM. Exam #2. Math 2360 Summer II, 2000 Morning Class. Nov. 15, 2000 ANSWERS
EXAM Exam # Math 6 Summer II Morning Class Nov 5 ANSWERS i Problem Consider the matrix 6 pts A = 6 4 9 5 7 6 5 5 5 4 The RREF of A is the matrix R = A Find a basis for the nullspace of A Solve the homogeneous
More informationModel Order Reduction by Parameter-Varying Oblique Projection
Model Order Reduction by Parameter-Varying Oblique Projection Julian heis, Peter Seiler and Herbert Werner Abstract A method to reduce the dynamic order of linear parameter-varying (LPV) systems in grid
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 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 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 informationEquilibrium points: continuous-time systems
Capitolo 0 INTRODUCTION 81 Equilibrium points: continuous-time systems Let us consider the following continuous-time linear system ẋ(t) Ax(t)+Bu(t) y(t) Cx(t)+Du(t) The equilibrium points x 0 of the system
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 informationNECESSARY AND SUFFICIENT CONDITIONS FOR STATE-SPACE NETWORK REALIZATION. Philip E. Paré Masters Thesis Defense
NECESSARY AND SUFFICIENT CONDITIONS FOR STATE-SPACE NETWORK REALIZATION Philip E. Paré Masters Thesis Defense INTRODUCTION MATHEMATICAL MODELING DYNAMIC MODELS A dynamic model has memory DYNAMIC MODELS
More informationChapter 1 Vector Spaces
Chapter 1 Vector Spaces Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 110 Linear Algebra Vector Spaces Definition A vector space V over a field
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 informationIterative methods for symmetric eigenvalue problems
s Iterative s for symmetric eigenvalue problems, PhD McMaster University School of Computational Engineering and Science February 11, 2008 s 1 The power and its variants Inverse power Rayleigh quotient
More informationECEN 605 LINEAR SYSTEMS. Lecture 8 Invariant Subspaces 1/26
1/26 ECEN 605 LINEAR SYSTEMS Lecture 8 Invariant Subspaces Subspaces Let ẋ(t) = A x(t) + B u(t) y(t) = C x(t) (1a) (1b) denote a dynamic system where X, U and Y denote n, r and m dimensional vector spaces,
More informationMethods for eigenvalue problems with applications in model order reduction
Methods for eigenvalue problems with applications in model order reduction Methoden voor eigenwaardeproblemen met toepassingen in model orde reductie (met een samenvatting in het Nederlands) Proefschrift
More informationOrder Reduction of a Distributed Parameter PEM Fuel Cell Model
Order Reduction of a Distributed Parameter PEM Fuel Cell Model María Sarmiento Carnevali 1 Carles Batlle Arnau 2 Maria Serra Prat 2,3 Immaculada Massana Hugas 2 (1) Intelligent Energy Limited, (2) Universitat
More informationBalanced realization and model order reduction for nonlinear systems based on singular value analysis
Balanced realization and model order reduction for nonlinear systems based on singular value analysis Kenji Fujimoto a, and Jacquelien M. A. Scherpen b a Department of Mechanical Science and Engineering
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 informationAccepted to appear in SIAM Review, 2015
Accepted to appear in SIAM Review, 2015 A SURVEY OF PROJECTION-BASED MODEL REDUCTION METHODS FOR PARAMETRIC DYNAMICAL SYSTEMS PETER BENNER, SERKAN GUGERCIN, AND KAREN WILLCOX Abstract. Numerical simulation
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 informationDynamical systems: basic concepts
Dynamical systems: basic concepts Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 2014-2015
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationThe amount of work to construct each new guess from the previous one should be a small multiple of the number of nonzeros in A.
AMSC/CMSC 661 Scientific Computing II Spring 2005 Solution of Sparse Linear Systems Part 2: Iterative methods Dianne P. O Leary c 2005 Solving Sparse Linear Systems: Iterative methods The plan: Iterative
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Methods for Nonlinear Systems Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu 1 / 65 Outline 1 Nested Approximations 2 Trajectory PieceWise Linear (TPWL)
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
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 informationBlock Bidiagonal Decomposition and Least Squares Problems
Block Bidiagonal Decomposition and Least Squares Problems Åke Björck Department of Mathematics Linköping University Perspectives in Numerical Analysis, Helsinki, May 27 29, 2008 Outline Bidiagonal Decomposition
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationMath 314/ Exam 2 Blue Exam Solutions December 4, 2008 Instructor: Dr. S. Cooper. Name:
Math 34/84 - Exam Blue Exam Solutions December 4, 8 Instructor: Dr. S. Cooper Name: Read each question carefully. Be sure to show all of your work and not just your final conclusion. You may not use your
More informationMATH 323 Linear Algebra Lecture 12: Basis of a vector space (continued). Rank and nullity of a matrix.
MATH 323 Linear Algebra Lecture 12: Basis of a vector space (continued). Rank and nullity of a matrix. Basis Definition. Let V be a vector space. A linearly independent spanning set for V is called a basis.
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 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 informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report June 2009
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 09-21 June 2009 Model order reduction for nonlinear IC models with POD by A. Verhoeven, M. Striebel, E.J.W.
More informationModel Order Reduction of Continuous LTI Large Descriptor System Using LRCF-ADI and Square Root Balanced Truncation
, July 1-3, 15, London, U.K. Model Order Reduction of Continuous LI Large Descriptor System Using LRCF-ADI and Square Root Balanced runcation Mehrab Hossain Likhon, Shamsil Arifeen, and Mohammad Sahadet
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016 1. Please write your 1- or 2-digit exam number on
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 informationSJÄLVSTÄNDIGA ARBETEN I MATEMATIK
SJÄLVSTÄNDIGA ARBETEN I MATEMATIK MATEMATISKA INSTITUTIONEN, STOCKHOLMS UNIVERSITET Model Reduction for Piezo-Mechanical Systems Using Balanced Truncation av Mohammad Monir Uddin 2011 - No 5 Advisor: Prof.
More informationMathematics for Control Theory
Mathematics for Control Theory Outline of Dissipativity and Passivity Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials Only as a reference: Charles A. Desoer
More informationLyapunov functions on nonlinear spaces S 1 R R + V =(1 cos ) Constructing Lyapunov functions: a personal journey
Constructing Lyapunov functions: a personal journey Lyapunov functions on nonlinear spaces R. Sepulchre -- University of Liege, Belgium Reykjavik - July 2013 Lyap functions in linear spaces (1994-1997)
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 informationReduced Basis Method for Parametric
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
More informationSampling and Low-Rank Tensor Approximations
Sampling and Low-Rank Tensor Approximations Hermann G. Matthies Alexander Litvinenko, Tarek A. El-Moshely +, Brunswick, Germany + MIT, Cambridge, MA, USA wire@tu-bs.de http://www.wire.tu-bs.de $Id: 2_Sydney-MCQMC.tex,v.3
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationThe Newton-ADI Method for Large-Scale Algebraic Riccati Equations. Peter Benner.
The Newton-ADI Method for Large-Scale Algebraic Riccati Equations Mathematik in Industrie und Technik Fakultät für Mathematik Peter Benner benner@mathematik.tu-chemnitz.de Sonderforschungsbereich 393 S
More informationBALANCING-RELATED MODEL REDUCTION FOR LARGE-SCALE SYSTEMS
FOR LARGE-SCALE SYSTEMS Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Recent Advances in Model Order Reduction Workshop TU Eindhoven, November 23,
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 informationLecture 7: Positive Semidefinite Matrices
Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.
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 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 informationMODEL ORDER REDUCTION
MODEL ORDER REDUCTION Approximate yet accurate surrogates for large-scale simulation KAREN E. WILLCOX Science at Extreme Scales: Where Big Data Meets Large-Scale Computing Tutorials Institute for Pure
More informationStability of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
More informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 5. Input-Output Stability DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Input-Output Stability y = Hu H denotes
More information