Stability preserving post-processing methods applied in the Loewner framework
|
|
- Meryl Watson
- 5 years ago
- Views:
Transcription
1 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 Gosea and Athanasios C. Antoulas Jacobs University Bremen and Rice University Houston SPI 2016 Torino, Italy May 11, 2016
2 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 2 / 20 Outline 1 Introduction and Motivation 2 Loewner Framework 3 Post-Processing Methods 4 Numerical Examples
3 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 3 / 20 Introduction Model order reduction (MOR) is used to transform large, complex models of time dependent processes into smaller, simpler models that are still capable of representing accurately the behavior of the original process under a variety of conditions. Consider interpolatory model reduction methods: we seek reduced models whose transfer function matches that of the original system at selected interpolation points. Reduce dimension from n = dim(σ) to k = dim(ˆσ) { Eẋ(t) = Ax(t) + Bu(t) n k {Ê ˆx(t) = ˆx(t) + ˆBu(t) Σ : ˆΣ : y(t) = Cx(t) y(t) = Ĉ ˆx(t)
4 Main tool: Loewner framework data driven MOR method. Uses measured or computed data (e.g. measurements of the frequency response of a to-be approximated system) instead of the system matrices (E, A, B and C). Constructs reduced models based on a rank revealing factorization of appropriately constructed matrices. Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 4 / 20
5 Magnitude Main tool: Loewner framework data driven MOR method. Uses measured or computed data (e.g. measurements of the frequency response of a to-be approximated system) instead of the system matrices (E, A, B and C). Constructs reduced models based on a rank revealing factorization of appropriately constructed matrices. In some cases the resulting reduced systems may not be stable! Frequency response Original model Loewner model Frequency (ω) Poles of the Loewner reduced model Very good approximation but unstable model find a stable model which is an optimal or sub-optimal approximant of the original reduced Loewner model (in the H norm). Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 4 / 20
6 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 5 / 20 Loewner Framework Measured data: {(ω k, S k ) k = {1,, N}}, ω k C, S k C p m. { right data : (λ i, r i, w i ), i = {1,, ρ} Partition the data: ; left data : (µ j, l j, v j ), j = {1,, ν} Objective: Find H(s) C p m s.t. H(λ i )r i = w i, lj H(µ j) = vj
7 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 5 / 20 Loewner Framework Measured data: {(ω k, S k ) k = {1,, N}}, ω k C, S k C p m. { right data : (λ i, r i, w i ), i = {1,, ρ} Partition the data: ; left data : (µ j, l j, v j ), j = {1,, ν} Objective: Find H(s) C p m s.t. H(λ i )r i = w i, lj H(µ j) = vj The Loewner matrix L C ν ρ : v 1 r 1 l 1 w 1 µ 1 λ 1... L = v ν r 1 l ν w 1 µ ν λ 1... v 1 r ρ l 1 w ρ µ 1 λ ρ v ν r ρ l ν w ρ µ ν λ ρ V = ( v 1 v 2 v ν ) The shifted Loewner matrix L σ C ν ρ : µ 1 v 1 r 1 λ 1 l 1 w 1 µ 1 λ 1... L σ = µ ν v ν r 1 λ 1 l ν w 1 µ ν λ 1... W = ( w 1 w 2... w ρ ) µ 1 v 1 r ρ λ ρl 1 w ρ µ 1 λ ρ µ ν v ν r ρ λ ρl ν w ρ µ ν λ ρ
8 on Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 6 / 20 Theorem If (L σ, L) is regular, then {Ê = L, Â = L σ, ˆB = V, Ĉ = W } is a realization of the data. Hence, H(z) = W (L σ zl) 1 V is the required interpolant.
9 Theorem If (L σ, L) is regular, then {Ê = L, Â = L σ, ˆB = V, Ĉ = W } is a realization of the data. Hence, H(z) = W (L σ zl) 1 V is the required interpolant. Minimal amount of data: ( ) ( E A - L - Lσ B C V W ) Redundant data: ( ) ( E A -Y LX -Y L σ X B C Y V W X ) H(s) = W (L σ sl) 1 V Ĥ(s) = WX [Y (L σ sl)x ] 1 Y V Truncate at k = rank(l) SVD : [L L σ ] = Y Σ l X, [ L L σ ] = Ỹ Σ r X Theorem The quadruple {Ê = -Y LX, Â = -Y L σ X, ˆB = Y V, Ĉ = WX }, is the realization of an approximate data interpolant. Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 6 / 20
10 on Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 7 / 20 Reduction to order 2: L = L σ = [ H(µ1) H(λ ] 1) H(µ 1) H(λ 2 µ 1 λ 1 µ 1 λ 2 H(µ 2) H(λ 1) H(µ 2) H(λ 2) µ 2 λ 1 µ [ 2 λ 2 µ1h(µ 1) λ 1H(λ 1) µ 1H(µ 1) λ 2H(λ 2 µ 1 λ 1 µ 1 λ 2) µ 2H(µ 2) λ 1H(λ 1) µ 2H(µ 2) λ 2H(λ 2) µ 2 λ 1 µ 2 λ 2 [ ] H(µ1 ) V = H(µ 2 ) W = [ H(λ 1 ) H(λ 2 ) ] ] Interpolate 2k moments of the initial linear system (k=2): H(µ 1 ), H(µ 2 ), H(λ 1 ), H(λ 2 )
11 Post-Processing Methods Consider the following classes of linear systems Σ 0 n,p,m = {(E, A, B, C, D) Σ n,p,m ρ(e, A) ir = } Σ + n,p,m = {(E, A, B, C, D) Σ n,p,m ρ(e, A) C + } Σ n,p,m = {(E, A, B, C, D) Σ n,p,m ρ(e, A) C, } Σ + n,p,m and Σ n,p,m are the sets of the stable and antistable systems. Let Σ Σ 0 n,p,m. Find two systems so that Σ (Σ + Σ ) where: Σ + = (E +, A +, B +, C +, D) Σ + n +,p,m, Σ = (E, A, B, C, 0) Σ n,p,m. Define the space of all rational matrix transfer functions with bounded H norm (where G = sup ω R G(iω) 2 ):. RH = {G R(s) p m, ir D(G), G < } Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 8 / 20
12 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 9 / 20 Post-Processing Methods Definition (The (AP ) Problem) Given an unstable system Σ Σ 0 n,p,m with transfer function G, find a stable system ˆΣ whose transfer function is the best ñ N Σ ñ,p,m approximation of G in the space RH. G(Σ) G(ˆΣ) = inf G Σ G Σ Σ Σ ñ N ñ,p,m solving (AP ) for an unstable system Σ is equivalent to solving (AP ) for its antistable part Σ +. For standard systems Σ, the reversed is also known as an Nehari Problem.
13 Theorem (Main Result) The unstable reduced order Loewner model is decomposed: Σ L = Σ }{{} {E, A, B, C } Σ }{{} + {E +, A +, B +, C +} Σ 0 k,p,m Let the infinite gramians P + and Q + and σ 1 the largest Hankel singular value (i.e. σ 1 = max ( eig(p + Q + ) ). Take γ σ 1. R + = Q + E + P + E T + γ 2 I, Ê + = E T + R +, ˆB + = E + Q + B + Ĉ + = C + P + E T +, Â + = A T +R + C T + Ĉ+ It follows that: ˆΣ L = Σ }{{} {E, A, B, C } ˆΣ+ }{{} {Ê +, Â +, ˆB +, Ĉ +} Σ ˆk,p,m and the inequality holds: σ 1 Σ L ˆΣ L min(γ, Σ + ) Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 10 / 20
14 1 Σ = stable part of Σ L. 2 Σ opt = optimal stable approx. of Σ L in RH. 3 Σ γ sopt = sub-optimal stable approx. of Σ L in RH (for γ > σ 1 ). 4 Σ flp = flipping the unstable poles approx. of Σ L. Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 11 / 20
15 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 11 / 20 1 Σ = stable part of Σ L. 2 Σ opt = optimal stable approx. of Σ L in RH. 3 Σ γ sopt = sub-optimal stable approx. of Σ L in RH (for γ > σ 1 ). 4 Σ flp = flipping the unstable poles approx. of Σ L. The Loewner model transfer function is given (σ 1 = σ 2 = 1 6 ) H L (s) = } s {{ + 3 } s 2 1 }{{ s 1 } H (s) H +(s) H (s) = 1 s+3, Σ L Σ H = 1 3 H opt (s) = 1 s+3 1 6, Σ L Σ opt H = 1 6 H flp (s) = 1 s s+2 1 s+1, Σ L Σ flp H = 2 3 H sopt (s) = 1 s+3 s 6(8s 2 +25s+16), Σ L Σ sopt H = , γ = 7 6
16 on Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 12 / 20 Scattering parameters from EM simulation (Three coupled lines) Data acquired by courtesy of Prof. Dr. Grivet-Talocia, Politecnico di Torino (presentation held at MORML 16, Sttutgart Causality Verification for Data-Driven Macromodeling ) given sample pairs {(ω k, S k })
17 Scattering parameters from EM simulation (Three coupled lines) Data acquired by courtesy of Prof. Dr. Grivet-Talocia, Politecnico di Torino (presentation held at MORML 16, Sttutgart Causality Verification for Data-Driven Macromodeling ) given sample pairs {(ω k, S k }) 10 0 Singular values of the Loewner matrix Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 12 / 20
18 Scattering parameters from EM simulation (Three coupled lines) Data acquired by courtesy of Prof. Dr. Grivet-Talocia, Politecnico di Torino (presentation held at MORML 16, Sttutgart Causality Verification for Data-Driven Macromodeling ) given sample pairs {(ω k, S k }) 10 0 Singular values of the Loewner matrix Poles of the Loewner reduced model Original data Loewner model Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 12 / 20
19 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 13 / 20 Scattering parameters from EM simulation (Three coupled lines) Use truncation order k = 14. In this case, the Loewner model has 6 unstable poles (4 real and 2 complex conjugates) G(iω) 0.04 Magnitude Original data Loewner: Σ L Optimal H : Σ opt Suboptimal H : Σ γ sopt Flipped: Σ flp Frequency(ω)
20 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 14 / 20 Scattering parameters from EM simulation (Three coupled lines) Σ L Σ opt Σ L Σ γ sopt Σ L Σ Σ L Σ flp Frequency response of the error systems ΣL Σopt ΣL Σ γ sopt ΣL Σflp Magnitude Frequency(ω) Notice that as ω, the norm of the optimal error system stays constant while the norm of the sub-optimal error system 0.
21 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 15 / 20 Scattering parameters from EM simulation (Three coupled lines) Vary the offset γ σ 1 within (10 5, 10 1 ); 10 1 H norm of the error systems 10 2 Σ Σ opt H Σ Σ H Σ Σ γ sopt H Σ Σ flp H γ σ 1 lim γ σ1 Σ L Σ γ sopt = Σ L Σ opt = σ 1 lim γ Σ L Σ γ sopt = Σ L Σ = Σ +
22 7135 th order power system [FRM08] For the II nd experiment, consider a stable descriptor power system of dimension n = 7135 with two inputs and two outputs, i.e. m = p = Singular values of the Loewner matrix By analyzing the decay of the singular values, choose for instance reduction order k = 82. The reduced Loewner model has 9 unstable poles. Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston / 20
23 7135 th order power system [FRM08] 80 Number of unstable poles for different reduction orders Reduction order (k) σ1(g(iω)) σ2(g(iω)) Magnitude 10 2 Magnitude Frequency(ω) σ3(g(iω)) Original samples Loewner: ΣL Suboptimal H : Σ γ sopt Flipped: Σflp 10 1 Frequency(ω) σ4(g(iω)) Magnitude Frequency(ω) Magnitude Frequency(ω) on Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 17 / 20
24 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 18 / th order power system [FRM08] Magnitude 10 0 σ1 Magnitude100 σ Frequency(ω) σ3 ΣL Σ γ sopt H ΣL Σflp H 10 1 Frequency(ω) σ4 Magnitude 10 2 Magnitude 10 1 Frequency(ω) 10 1 Frequency(ω)
25 Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston 19 / 20 M. Koehler, On the closest stable descriptor system in the respective spaces RH 2 and RH. Linear Algebra and its Applications, 443(0):34-49, F. Freitas, J. Rommes, N. Martinsem, Gramian-Based Reduction Method Applied to Large Sparse Power System Descriptor Models IEEE Transactions on Power Systems, Vol. 23, Issue 3, August 2008, pp
26 Thank you for your attention! Any questions? Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016Houston / 20
Model 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 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 informationCase study: Approximations of the Bessel Function
Case study: Approximations of the Bessel Function DS Karachalios, IV Gosea, Q Zhang, AC Antoulas May 5, 218 arxiv:181339v1 [mathna] 31 Dec 217 Abstract The purpose of this note is to compare various approximation
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 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 informationOptimal RH 2 - and RH -Approximation of Unstable Descriptor Systems
Optimal RH - and RH -Approximation of Unstable Descriptor Systems Marcus Köhler August, 1 arxiv:13.595v [math.oc] 1 Aug 1 Stability perserving is an important topic in approximation of systems, e.g. model
More informationData-driven model order reduction of linear switched systems
Data-driven model order reduction of linear switched systems IV Gosea, M Petreczky, AC Antoulas arxiv:171205740v1 mathna 15 Dec 2017 Abstract The Loewner framework for model reduction is extended to the
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 informationModel reduction of coupled systems
Model reduction of coupled systems Tatjana Stykel Technische Universität Berlin ( joint work with Timo Reis, TU Kaiserslautern ) Model order reduction, coupled problems and optimization Lorentz Center,
More informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview and some new results Thanos Antoulas Rice University email: aca@rice.edu URL: www.ece.rice.edu/ aca LinSys2008, Sde Boker, 15-19 September 2008 Thanos
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 informationApproximation of the Linearized Boussinesq Equations
Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan
More informationKrylov Techniques for Model Reduction of Second-Order Systems
Krylov Techniques for Model Reduction of Second-Order Systems A Vandendorpe and P Van Dooren February 4, 2004 Abstract The purpose of this paper is to present a Krylov technique for model reduction of
More 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 informationLarge-scale and infinite dimensional dynamical systems approximation
Large-scale and infinite dimensional dynamical systems approximation Igor PONTES DUFF PEREIRA Doctorant 3ème année - ONERA/DCSD Directeur de thèse: Charles POUSSOT-VASSAL Co-encadrant: Cédric SEREN 1 What
More 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 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 informationComparison of Model Reduction Methods with Applications to Circuit Simulation
Comparison of Model Reduction Methods with Applications to Circuit Simulation Roxana Ionutiu, Sanda Lefteriu, and Athanasios C. Antoulas Department of Electrical and Computer Engineering, Rice University,
More 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 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 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 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 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 informationGramian Based Model Reduction of Large-scale Dynamical Systems
Gramian Based Model Reduction of Large-scale Dynamical Systems Paul Van Dooren Université catholique de Louvain Center for Systems Engineering and Applied Mechanics B-1348 Louvain-la-Neuve, Belgium vdooren@csam.ucl.ac.be
More 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 informationQualification of tabulated scattering parameters
Qualification of tabulated scattering parameters Stefano Grivet Talocia Politecnico di Torino, Italy IdemWorks s.r.l. stefano.grivet@polito.it 4 th IEEE Workshop on Signal Propagation on Interconnects
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 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 informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview and some new results Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/ aca Electrical and Computer Engineering
More informationParallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation
Parallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation Peter Benner 1, Enrique S. Quintana-Ortí 2, Gregorio Quintana-Ortí 2 1 Fakultät für Mathematik Technische Universität
More informationInfinite dimensional model analysis through data-driven interpolatory techniques
Infinite dimensional model analysis through data-driven interpolatory techniques... applications to time-delayed and irrational functions C. Poussot-Vassal 20th Conference of the International Linear Algebra
More informationH 2 optimal model approximation by structured time-delay reduced order models
H 2 optimal model approximation by structured time-delay reduced order models I. Pontes Duff, Charles Poussot-Vassal, C. Seren 27th September, GT MOSAR and GT SAR meeting 200 150 100 50 0 50 100 150 200
More 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 informationDESIGN OF OBSERVERS FOR SYSTEMS WITH SLOW AND FAST MODES
DESIGN OF OBSERVERS FOR SYSTEMS WITH SLOW AND FAST MODES by HEONJONG YOO A thesis submitted to the Graduate School-New Brunswick Rutgers, The State University of New Jersey In partial fulfillment of the
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 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 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 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 informationLQR and H 2 control problems
LQR and H 2 control problems Domenico Prattichizzo DII, University of Siena, Italy MTNS - July 5-9, 2010 D. Prattichizzo (Siena, Italy) MTNS - July 5-9, 2010 1 / 23 Geometric Control Theory for Linear
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 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 informationL 2 -optimal model reduction for unstable systems using iterative rational
L 2 -optimal model reduction for unstable systems using iterative rational Krylov algorithm Caleb agruder, Christopher Beattie, Serkan Gugercin Abstract Unstable dynamical systems can be viewed from a
More information6 Optimal Model Order Reduction in the Hankel Norm
6 Optimal Model Order Reduction in the Hankel Norm In this lecture, we discuss optimal model reduction in the Hankel norm. There is no known solution to the optimal model reduction problem in the H -norm.
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 informationThe Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals
The Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals Timo Hämäläinen Seppo Pohjolainen Tampere University of Technology Department of Mathematics
More informationVISA: Versatile Impulse Structure Approximation for Time-Domain Linear Macromodeling
VISA: Versatile Impulse Structure Approximation for Time-Domain Linear Macromodeling Chi-Un Lei and Ngai Wong Dept. of Electrical and Electronic Engineering The University of Hong Kong Presented by C.U.
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More informationarxiv: v1 [cs.sy] 17 Nov 2015
Optimal H model approximation based on multiple input/output delays systems I. Pontes Duff a, C. Poussot-Vassal b, C. Seren b a ISAE SUPAERO & ONERA b ONERA - The French Aerospace Lab, F-31055 Toulouse,
More 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 informationWavelet-Based Passivity Preserving Model Order Reduction for Wideband Interconnect Characterization
Wavelet-Based Passivity Preserving Model Order Reduction for Wideband Interconnect Characterization Mehboob Alam, Arthur Nieuwoudt, and Yehia Massoud Rice Automated Nanoscale Design Group Rice University,
More informationarxiv: v1 [math.na] 12 Jan 2012
Interpolatory Weighted-H 2 Model Reduction Branimir Anić, Christopher A. Beattie, Serkan Gugercin, and Athanasios C. Antoulas Preprint submitted to Automatica, 2 January 22 arxiv:2.2592v [math.na] 2 Jan
More informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
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 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 informationGramian Based Model Reduction of Large-scale Dynamical Systems
Gramian Based Model Reduction of Large-scale Dynamical Systems Paul Van Dooren Université catholique de Louvain Center for Systems Engineering and Applied Mechanics B-1348 Louvain-la-Neuve, Belgium vdooren@csamuclacbe
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 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 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 informationInterpolatory Model Reduction of Large-scale Dynamical Systems
Interpolatory Model Reduction of Large-scale Dynamical Systems Athanasios C. Antoulas, Christopher A. Beattie, and Serkan Gugercin Abstract Large-scale simulations play a crucial role in the study of a
More informationAN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL MOTIVATING EXAMPLE INVERTED PENDULUM
Controls Lab AN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL Eduardo Gildin (UT ICES and Rice Univ.) with Thanos Antoulas (Rice ECE) Danny Sorensen (Rice
More informationRobust Control of Time-delay Systems
Robust Control of Time-delay Systems Qing-Chang Zhong Distinguished Lecturer, IEEE Power Electronics Society Max McGraw Endowed Chair Professor in Energy and Power Engineering Dept. of Electrical and Computer
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 informationOn the continuity of the J-spectral factorization mapping
On the continuity of the J-spectral factorization mapping Orest V. Iftime University of Groningen Department of Mathematics and Computing Science PO Box 800, 9700 AV Groningen, The Netherlands Email: O.V.Iftime@math.rug.nl
More informationPOLITECNICO DI TORINO Repository ISTITUZIONALE
POLITECNICO DI TORINO Repository ISTITUZIONALE High-Performance Passive Macromodeling Algorithms for Parallel Computing Platforms Original High-Performance Passive Macromodeling Algorithms for Parallel
More informationarxiv: v1 [math.oc] 17 Oct 2014
SiMpLIfy: A Toolbox for Structured Model Reduction Martin Biel, Farhad Farokhi, and Henrik Sandberg arxiv:1414613v1 [mathoc] 17 Oct 214 Abstract In this paper, we present a toolbox for structured model
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 informationBalanced truncation for linear switched systems
Adv Comput Math (218 44:1845 1886 https://doi.org/1.17/s1444-18-961-z Balanced truncation for linear switched systems Ion Victor Gosea 1 Mihaly Petreczky 2 Athanasios C. Antoulas 1,3,4 Christophe Fiter
More informationRealization-independent H 2 -approximation
Realization-independent H -approximation Christopher Beattie and Serkan Gugercin Abstract Iterative Rational Krylov Algorithm () of [] is an effective tool for tackling the H -optimal model reduction problem.
More 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 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 informationNumerical Computation of Structured Complex Stability Radii of Large-Scale Matrices and Pencils
51st IEEE onference on Decision and ontrol December 10-13, 01. Maui, Hawaii, USA Numerical omputation of Structured omplex Stability Radii of Large-Scale Matrices and Pencils Peter Benner and Matthias
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII
CONTENTS VOLUME VII Control of Linear Multivariable Systems 1 Katsuhisa Furuta,Tokyo Denki University, School of Science and Engineering, Ishizaka, Hatoyama, Saitama, Japan 1. Linear Multivariable Systems
More informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview Thanos Antoulas joint work with Dan Sorensen aca@rice.edu URL: www-ece.rice.edu/ aca CITI Lecture, Rice, February 2005 Model reduction of large-scale
More 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 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 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 informationSigned Measures and Complex Measures
Chapter 8 Signed Measures Complex Measures As opposed to the measures we have considered so far, a signed measure is allowed to take on both positive negative values. To be precise, if M is a σ-algebra
More informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More informationModern Control Systems
Modern Control Systems Matthew M. Peet Illinois Institute of Technology Lecture 18: Linear Causal Time-Invariant Operators Operators L 2 and ˆL 2 space Because L 2 (, ) and ˆL 2 are isomorphic, so are
More informationTechnische Universität Berlin
; Technische Universität Berlin Institut für Mathematik arxiv:1610.03262v2 [cs.sy] 2 Jan 2017 Model Reduction for Large-scale Dynamical Systems with Inhomogeneous Initial Conditions Christopher A. Beattie
More 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 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 informationRegularization and Inverse Problems
Regularization and Inverse Problems Caroline Sieger Host Institution: Universität Bremen Home Institution: Clemson University August 5, 2009 Caroline Sieger (Bremen and Clemson) Regularization and Inverse
More informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
More informationEL2520 Control Theory and Practice
So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal
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 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 informationUpper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza
Aalborg Universitet Upper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza Published in: 7th IFAC Symposium on Robust Control Design DOI (link
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 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 informationPseudospectra and Nonnormal Dynamical Systems
Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 1 Overview of the Course These lectures
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 informationRational interpolation, minimal realization and model reduction
Rational interpolation, minimal realization and model reduction Falk Ebert Tatjana Stykel Technical Report 371-2007 DFG Research Center Matheon Mathematics for key technologies http://www.matheon.de Rational
More informationComputation of a Compact State Space Model for an Adaptive Spindle Head Configuration with Piezo Actuators using Balanced Truncation
Prod. Eng. Res. Devel. manuscript No. (will be inserted by the editor) Computation of a Compact State Space Model for an Adaptive Spindle Head Configuration with Piezo Actuators using Balanced Truncation
More informationModel Boundary Approximation Method as a Unifying Framework for Balanced Truncation and Singular Perturbation Approximation
1 Model Boundary Approximation Method as a Unifying Framework for Balanced Truncation and Singular Perturbation Approximation Philip E. Paré, David Grimsman, Alma T. Wilson, Mark K. Transtrum and Sean
More informationModel Reduction using a Frequency-Limited H 2 -Cost
Technical report from Automatic Control at Linköpings universitet Model Reduction using a Frequency-Limited H 2 -Cost Daniel Petersson, Johan Löfberg Division of Automatic Control E-mail: petersson@isy.liu.se,
More informationJ-SPECTRAL FACTORIZATION
J-SPECTRAL FACTORIZATION of Regular Para-Hermitian Transfer Matrices Qing-Chang Zhong zhongqc@ieee.org School of Electronics University of Glamorgan United Kingdom Outline Notations and definitions Regular
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 informationPOLITECNICO DI TORINO Repository ISTITUZIONALE
POLITECNICO DI TORINO Repository ISTITUZIONALE Sensitivity-based weighting for passivity enforcement of linear macromodels in power integrity applications Original Sensitivity-based weighting for passivity
More informationOptimal Input Elimination
SICE Journal of Control, Measurement, and System Integration, Vol. 11, No. 2, pp. 100 104, March 2018 Optimal Input Elimination Kazuhiro SATO Abstract : This paper studies an optimal inputs elimination
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 information