Chris Beattie and Serkan Gugercin
|
|
- Andrew Dalton
- 5 years ago
- Views:
Transcription
1 Krylov-based model reduction of second-order systems with proportional damping Chris Beattie and Serkan Gugercin Department of Mathematics, Virginia Tech, USA 44th IEEE Conference on Decision and Control and European Control Conference CDC-ECC 2005 Seville, Spain, December
2 Introduction Consider the second-order system of the form Mẍ(t) + Gẋ(t) + Kx(t) = Bu(t), y(t) = Cx(t) H(s) = C(Ms 2 + Gs + K) 1 B M,G,K R n n, and B,C T R n. In many cases, n is too large for efficient simulation and control Generate, for some r n, an r th order reduced system M r ẍ r (t) + G r ẋ r (t) + K r x r (t) = B r u(t), y r (t) = C r x r (t) M r,g r,k r R r r, and B r,c T r R r y r (t) approximates y(t) for a wide range of inputs u(t). 2
3 H(s) := Krylov-based model reduction E q(t) = Aq(t) + Bu(t) H(s) = C(sE A) 1 B y(t) = Cq(t) Find H r (s) = C r (se r A r ) 1 B r so that d j H(s) ds j = dj H r (s) s=σk ds j, s=σk K : Number of interpolation points J : Number of interpolation conditions at each σ k. k = 1,...,K j = 0,...,J 1. E r := Z T EV, A r := Z T AV, B r := Z T B, C r := CV 3
4 Given F C N N and g C N, define the Krylov subspace: K J (F,g) = span{g, Fg, F 2 g,, F J 1 g } Construct V and Z such that Ran(V) = span {K J1 (F 1,g 1 ),, K JK (F K,g K )}, and Ran(Z) = span { K JK+1 (F K+1,g K+1 ),, K J2K (F 2K,g 2K ) }, where F i = (σ i E A) 1 E, g i = (σ i E A) 1 B, for i = 1,...,K F i = (σ i E A) T E T, g i = (σ i E A) T C T, for i = K + 1,...,2K H r (s) satisfies the interpolation conditions 4
5 Krylov reduction of second-order systems Convert Mẍ(t) + Gẋ(t) + Kx(t) = Bu(t), y(t) = Cx(t) into E q = Aq(t) + Bu(t), y(t) = Cq(t) where E = I 0 0 M, A = K 0 I αm βk, B = 0 B, C T = CT 0. Apply first-order Krylov reduction and transfer back. Not always possible and structure is lost = Apply reduction directly in the second-order system framework Find W R n r such that W T W = I r and M r = W T MW, G r = W T GW, K r = W T KW, B r = W T B, and C r = CW 5
6 Second-order Krylov Subspaces (Bai[2003]) Given F 1, F 2 C n n and g C n, second-order Krylov subspace: K (2) J (F 1, F 2,g) = span{r 0, r 1,, r J 1 } where r 0 = g, r 1 = F 1 r 0, r j = F 1 r j 1 + F 2 r j 2, for j 2. For second-order reduction with interpolation point σ: W r = K (2) r ( K G, KM,r 0 ) where K = σ 2 M + σg + K, D = 2σM + G, r0 = K 1 B Reduction directly in the second-order system framework Reduced model matches the first r moments at σ. 6
7 Second-order systems with proportional damping Mẍ(t) + (αm + βk)ẋ(t) + Kx(t) = Bu(t), where α, β > 0, and αβ < 1. y(t) = Cx(t) E = First-order equivalent system: H(s) = C (se A) 1 B with I 0 0 M, A = K 0 I αm βk, B = 0 B, C T = CT 0. To match the first r moments at σ, construct V 2n r = span { g, Fg,, F r 1 g } where F = (σe A) 1 E and g = (σe A) 1 B. 7
8 Theorem: Given the above set-up, V 2n r W n r W n r where W n r = K r (K 1 σ M, K 1 σ B) with K σ = σ 2 M + σ(αm + βk) + K. Model reduction directly in the second-order framework using Ran(W r ) = W with W T r W r = I r, W r R n r H r (s) matches r moments at σ. Second-order reduction using regular first-order Krylov subspace W r If σ = 0, W n r = span{k 1 B,, (K 1 M) r 1 K 1 B} If σ =, W n r = span{m 1 B,, (M 1 K) r 1 M 1 B} 8
9 Proof: g = (σe A) 1 B = K 1 σ B σk 1 σ B and K 1 σ B W. Let v = v 1 and v 1,v 2 K p (K 1 σ M, K 1 σ B). v 2 The next vector in W is v = (σe A) 1 Ev = v 1 v 2 = κ 1v 1 + κ 2 K 1 σ Mv 1 + κ 3 K 1 σ Mv 2 κ 4 v 1 + κ 5 K 1 σ Mv 1 + κ 6 K 1 σ Mv 2 v 1, v 2 K p+1 (K 1 σ M, K 1 σ B) Moment matching: Second-order reduction using Ran(W) = W is equivalent to first-order reduction with V = Z = W 0. 0 W 9
10 Approximation by Interpolation (Beattie[2004]) Best uniform approximation H r (s) that minimizes max ω R H(jω) H r (jω) should make H(jω) H r (jω) constant as r ( near-circularity of best uniform rational approximation error (Trefethen, 1981) related to Chebyshev equioscillation theorem. Best uniform approximations are hard to calculate. Interpolants are easy to calculate. How to choose the shifts? Pick interpolation points carefully to recover a good uniform approximant. log H(z) H r (z) has positive singularities at system eigenvalues. negative singularities at interpolation points. 10
11 Pick interpolation points to balance the contours of log H(z) H r (z) (makes log H(z) H r (z) nearly constant along the imaginary axis) Interpolation at Ritz values mirrored across the imaginary axis is good. Mirror Ritz values: optimal choice for H 2 minimization as well (G./Antoulas/Beattie [2004]) If system spectra is structured (e.g., circular), interpolation at equivalent lumped charge locations will have much the same effect. 11
12 Shift Selection for Proportional Damping Mẍ + (αm + βk)ẋ + Kx = p. Proposition: All damped eigenvalues (the system poles) are on circle with center: 1 β, radius: 1 αβ β and on ray (, 1 β ]. Distribution depends on undamped natural frequencies, but usual elastic vibration models lead to distributions that are close to equilibrium condenser distributions Interpret log H(z) H r (z) as potential function associated with charge distribution (poles have net charge of +1, interpolation points have net charge of 1). 12
13 Only ONE shift is necessary - replace aggregate of interpolation points (negative charge distribution) with single α shift (an equivalent lumped charge) at σ = β. Optimal choice for condenser distribution of system poles; pretty good choice for most K and M). α Regular first-order Krylov subspace using a single shift: β. 13
14 Pick α, β (0, 1) Exact Condenser Distribution K = α β 2 1 αβ 1 αβ αβ αβ αβ 1 αβ M = 2+ 1 αβ 1 αβ αβ αβ αβ 1 αβ. G = αm + β K, α = β = 0.05 B = C T = [ ] T. 14
15 Reduction from n = 2000 to r = 30 using a single shift 20 Pole locations for exact condenser distribution 10 1 H error vs interpolation point Imag 0 5 H error Real σ σ = α β = 1 is the optimal shift. 15
16 A 1-D Beam Model: n = α = 1/10, β = 1/500, B = e 1, C = e T 200. H error vs interpolation point 500 Distributiion of Observed and Condenser Poles Observed Poles Condenser Poles 10 1 Imag G(s) G r (s) r=5 r=10 r= Real σ σ = α β = : Very close to being optimal. 16
17 Another 1-D Beam Model: n = 200 and α = β = 1/300, B = C T = e 1. Compare with balanced truncation and other shift selections Reduction done in the first-order framework to r (1) = 40 (r = 20) Amplitude bode plots for the beam model 10 4 Amplitude bode plots of the error systems for the beam model H H BT 10 4 H σ* 10 5 H( jw ) H BT H(s) H error ( jw ) H H σ* freq (rad/sec) freq (rad/sec) 17
18 Amplitude bode plots for the beam model 3.6 Beam with n=200, α = β = 1/ σ = σ = 5 H( jw ) σ = 1 σ i, i=1:40 log 10 H er Hinf Observed Convergence rate: Expected Convergence rate: σ = freq (rad/sec) r 18
19 Conclusions Considered Mẍ + (αm + βk)ẋ + Kx = Bu(t), y(t) = Cx(t) Second-order Krylov reduction using first-order Krylov subspaces Pole locations for proportional damping Equivalent single shift (lumped charge) for aggregate of interpolation points (negative charge distribution) Optimal single shift σ = for condenser distribution Close to optimal in general Future work: Extensions to other types of damping. β α 19
Krylov-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 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 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 informationOptimal Approximation of Dynamical Systems with Rational Krylov Methods. Department of Mathematics, Virginia Tech, USA
Optimal Approximation of Dynamical Systems with Rational Krylov Methods Serkan Güǧercin Department of Mathematics, Virginia Tech, USA Chris Beattie and Virginia Tech. jointly with Thanos Antoulas Rice
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 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 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 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 informationarxiv: v1 [math.na] 12 Jan 2012
Interpolatory Weighted-H 2 Model Reduction Branimir Anić, Christopher A. Beattie, Serkan Gugercin, and Athanasios C. Antoulas Preprint submitted to Automatica, 2 January 22 arxiv:2.2592v [math.na] 2 Jan
More informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 WeC10.4 An iterative SVD-Krylov based method for model reduction
More 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 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 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 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 informationCollocated versus non-collocated control [H04Q7]
Collocated versus non-collocated control [H04Q7] Jan Swevers September 2008 0-0 Contents Some concepts of structural dynamics Collocated versus non-collocated control Summary This lecture is based on parts
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 informationOn Dominant Poles and Model Reduction of Second Order Time-Delay Systems
On Dominant Poles and Model Reduction of Second Order Time-Delay Systems Maryam Saadvandi Joint work with: Prof. Karl Meerbergen and Dr. Elias Jarlebring Department of Computer Science, KULeuven ModRed
More informationK(s +2) s +20 K (s + 10)(s +1) 2. (c) KG(s) = K(s + 10)(s +1) (s + 100)(s +5) 3. Solution : (a) KG(s) = s +20 = K s s
321 16. Determine the range of K for which each of the following systems is stable by making a Bode plot for K = 1 and imagining the magnitude plot sliding up or down until instability results. Verify
More informationEfficient computation of transfer function dominant poles of large second-order dynamical systems
Chapter 6 Efficient computation of transfer function dominant poles of large second-order dynamical systems Abstract This chapter presents a new algorithm for the computation of dominant poles of transfer
More informationParametric Order Reduction of Proportionally Damped Second Order Systems
RT Lehrstuhl für Regelungstechnik Technical Reports on Automatic Control Vol. TRAC-1, Oct. 2006 Available online at: www.rt.mw.tum.de TECHNISCHE UNIVERSITÄT MÜNCHEN Parametric Order Reduction of Proportionally
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 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 informationSchool of Mechanical Engineering Purdue University
Case Study ME375 Frequency Response - 1 Case Study SUPPORT POWER WIRE DROPPERS Electric train derives power through a pantograph, which contacts the power wire, which is suspended from a catenary. During
More informationLecture 9 Infinite Impulse Response Filters
Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9
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 informationDynamics of structures
Dynamics of structures 1.2 Viscous damping Luc St-Pierre October 30, 2017 1 / 22 Summary so far We analysed the spring-mass system and found that its motion is governed by: mẍ(t) + kx(t) = 0 k y m x x
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Serkan Gugercin, Tatjana Stykel, Sarah Wyatt Model Reduction of Descriptor Systems by Interpolatory Projection Methods Preprint Nr. 3/213 29.
More informationAutomatic Order Reduction for Finite Element Models
Automatic Order Reduction for Finite Element Models Andreas Greiner, Jan Lienemann, Evgenii Rudnyi, and Jan G. Korvink IMTEK, Albert-Ludwigs-Universität Freiburg Geroges-Köhler-Allee 103, D-79110 Freiburg,
More informationME 375 EXAM #1 Friday, March 13, 2015 SOLUTION
ME 375 EXAM #1 Friday, March 13, 2015 SOLUTION PROBLEM 1 A system is made up of a homogeneous disk (of mass m and outer radius R), particle A (of mass m) and particle B (of mass m). The disk is pinned
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 informationSTABILITY ANALYSIS TECHNIQUES
ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response
More informationL 2 -optimal model reduction for unstable systems using iterative rational
L 2 -optimal model reduction for unstable systems using iterative rational Krylov algorithm Caleb agruder, Christopher Beattie, Serkan Gugercin Abstract Unstable dynamical systems can be viewed from a
More informationModel 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 informationExercises Lecture 15
AM1 Mathematical Analysis 1 Oct. 011 Feb. 01 Date: January 7 Exercises Lecture 15 Harmonic Oscillators In classical mechanics, a harmonic oscillator is a system that, when displaced from its equilibrium
More informationInexact Solves in Interpolatory Model Reduction
Inexact Solves in Interpolatory Model Reduction Sarah Wyatt Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements for the
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using
More informationAn Interpolation-Based Approach to the Optimal H Model Reduction
An Interpolation-Based Approach to the Optimal H Model Reduction Garret M. Flagg Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment
More informationPositioning Servo Design Example
Positioning Servo Design Example 1 Goal. The goal in this design example is to design a control system that will be used in a pick-and-place robot to move the link of a robot between two positions. Usually
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 20) Final Examination December 9, 20 Name: Kerberos Username: Please circle your section number: Section Time 2 am pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
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 informationInexact Solves in Krylov-based Model Reduction
Inexact Solves in Krylov-based Model Reduction Christopher A. Beattie and Serkan Gugercin Abstract We investigate the use of inexact solves in a Krylov-based model reduction setting and present the resulting
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using the
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 informationOverview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros)
Overview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros) J. McNames Portland State University ECE 222 Bode Plots Ver.
More informationEE Experiment 11 The Laplace Transform and Control System Characteristics
EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal
More informationMAS107 Control Theory Exam Solutions 2008
MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve
More informationDownloaded 05/27/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 35, No. 5, pp. B1010 B1033 c 2013 Society for Industrial and Applied Mathematics MODEL REDUCTION OF DESCRIPTOR SYSTEMS BY INTERPOLATORY PROJECTION METHODS SERKAN GUGERCIN, TATJANA
More information6.003: Signals and Systems. CT Fourier Transform
6.003: Signals and Systems CT Fourier Transform April 8, 200 CT Fourier Transform Representing signals by their frequency content. X(jω)= x(t)e jωt dt ( analysis equation) x(t)= 2π X(jω)e jωt dω ( synthesis
More information6.003: Signals and Systems. Applications of Fourier Transforms
6.003: Signals and Systems Applications of Fourier Transforms November 7, 20 Filtering Notion of a filter. LTI systems cannot create new frequencies. can only scale magnitudes and shift phases of existing
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 informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationLike bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.
Inversion of the z-transform Focus on rational z-transform of z 1. Apply partial fraction expansion. Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform. Let X(z)
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 information6.003: Signals and Systems. CT Fourier Transform
6.003: Signals and Systems CT Fourier Transform April 8, 200 CT Fourier Transform Representing signals by their frequency content. X(jω)= x(t)e jωt dt ( analysis equation) x(t)= X(jω)e jωt dω ( synthesis
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
More informationModeling and Experimentation: Mass-Spring-Damper System Dynamics
Modeling and Experimentation: Mass-Spring-Damper System Dynamics Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin July 20, 2014 Overview 1 This lab is meant to
More informationDynamics of Structures
Dynamics of Structures Elements of structural dynamics Roberto Tomasi 11.05.2017 Roberto Tomasi Dynamics of Structures 11.05.2017 1 / 22 Overview 1 SDOF system SDOF system Equation of motion Response spectrum
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 22: The Nyquist Criterion Overview In this Lecture, you will learn: Complex Analysis The Argument Principle The Contour
More informationProblem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013
EE 56: Digital Control Systems Problem Set 3: Solution Due on Mon 7 th Oct in class Fall 23 Problem For the causal LTI system described by the difference equation y k + 2 y k = x k, () (a) By first finding
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 information6.003 Homework #6 Solutions
6.3 Homework #6 Solutions Problems. Maximum gain For each of the following systems, find the frequency ω m for which the magnitude of the gain is greatest. a. + s + s ω m = w This system has poles at s
More informationkx m x B N 1 C L, M Mg θ
.004 MODELING DYNAMICS AND CONTROL II Spring 00 Solutions to Problem Set No. 7 Problem 1. Pendulum mounted on elastic support. This problem is an execise in the application of momentum principles. Two
More informationTopic # Feedback Control Systems
Topic #20 16.31 Feedback Control Systems Closed-loop system analysis Bounded Gain Theorem Robust Stability Fall 2007 16.31 20 1 SISO Performance Objectives Basic setup: d i d o r u y G c (s) G(s) n control
More informationTopic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback
Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall
More informationSolutions for homework 5
1 Section 4.3 Solutions for homework 5 17. The following equation has repeated, real, characteristic roots. Find the general solution. y 4y + 4y = 0. The characteristic equation is λ 4λ + 4 = 0 which has
More informationSection 4.9; Section 5.6. June 30, Free Mechanical Vibrations/Couple Mass-Spring System
Section 4.9; Section 5.6 Free Mechanical Vibrations/Couple Mass-Spring System June 30, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This Session: (1) Free
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 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 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 informationTransform Solutions to LTI Systems Part 3
Transform Solutions to LTI Systems Part 3 Example of second order system solution: Same example with increased damping: k=5 N/m, b=6 Ns/m, F=2 N, m=1 Kg Given x(0) = 0, x (0) = 0, find x(t). The revised
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More informationFrequency Response of Linear Time Invariant Systems
ME 328, Spring 203, Prof. Rajamani, University of Minnesota Frequency Response of Linear Time Invariant Systems Complex Numbers: Recall that every complex number has a magnitude and a phase. Example: z
More informationLectures Chapter 10 (Cutnell & Johnson, Physics 7 th edition)
PH 201-4A spring 2007 Simple Harmonic Motion Lectures 24-25 Chapter 10 (Cutnell & Johnson, Physics 7 th edition) 1 The Ideal Spring Springs are objects that exhibit elastic behavior. It will return back
More informationChapter 2 SDOF Vibration Control 2.1 Transfer Function
Chapter SDOF Vibration Control.1 Transfer Function mx ɺɺ( t) + cxɺ ( t) + kx( t) = F( t) Defines the transfer function as output over input X ( s) 1 = G( s) = (1.39) F( s) ms + cs + k s is a complex number:
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 2011) Final Examination December 19, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 1 pm 4 2 pm Grades will be determined by the correctness of your answers
More informationME scope Application Note 28
App Note 8 www.vibetech.com 3/7/17 ME scope Application Note 8 Mathematics of a Mass-Spring-Damper System INTRODUCTION In this note, the capabilities of ME scope will be used to build a model of the mass-spring-damper
More informationClustering-Based Model Order Reduction for Multi-Agent Systems with General Linear Time-Invariant Agents
Max Planck Institute Magdeburg Preprints Petar Mlinarić Sara Grundel Peter Benner Clustering-Based Model Order Reduction for Multi-Agent Systems with General Linear Time-Invariant Agents MAX PLANCK INSTITUT
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More informationChapter 14 Oscillations. Copyright 2009 Pearson Education, Inc.
Chapter 14 Oscillations Oscillations of a Spring Simple Harmonic Motion Energy in the Simple Harmonic Oscillator Simple Harmonic Motion Related to Uniform Circular Motion The Simple Pendulum The Physical
More informationNormalization of Complex Mode Shapes by Truncation of the. Alpha-Polynomial. A thesis submitted to the. Graduate School
Normalization of Complex Mode Shapes by Truncation of the Alpha-Polynomial A thesis submitted to the Graduate School of the University of Cincinnati in partial fulfillment of the requirements for the degree
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 informationIssues in Interpolatory Model Reduction: Inexact Solves, Second-order Systems and DAEs
Issues in Interpolatory Model Reduction: Inexact Solves, Second-order Systems and DAEs Sarah Wyatt Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial
More 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 informationAutomatic Control 2. Model reduction. Prof. Alberto Bemporad. University of Trento. Academic year
Lecture: Automatic Control 2 Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 17 Lecture:
More informationChapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx
Chapter 1 Lecture Notes Chapter 1 Oscillatory Motion Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx When the mass is released, the spring will pull
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method
.. AERO 422: Active Controls for Aerospace Vehicles Frequency Response- Method Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. ... Response to
More informationStructure-preserving tangential interpolation for model reduction of port-hamiltonian Systems
Structure-preserving tangential interpolation for model reduction of port-hamiltonian Systems Serkan Gugercin a, Rostyslav V. Polyuga b, Christopher Beattie a, and Arjan van der Schaft c arxiv:.3485v2
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationTopic # Feedback Control
Topic #5 6.3 Feedback Control State-Space Systems Full-state Feedback Control How do we change the poles of the state-space system? Or,evenifwecanchangethepolelocations. Where do we put the poles? Linear
More informationIntroduction to Controls
EE 474 Review Exam 1 Name Answer each of the questions. Show your work. Note were essay-type answers are requested. Answer with complete sentences. Incomplete sentences will count heavily against the grade.
More informationDynamic circuits: Frequency domain analysis
Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.011: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 2010 ANSWER BOOKLET
More informationEE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley
University of California Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences Professors : N.Morgan / B.Gold EE225D Digital Filters Spring,1999 Lecture 7 N.MORGAN
More informationNetwork Clustering for SISO Linear Dynamical Networks via Reaction-Diffusion Transformation
Milano (Italy) August 28 - September 2, 211 Network Clustering for SISO Linear Dynamical Networks via Reaction-Diffusion Transformation Takayuki Ishizaki Kenji Kashima Jun-ichi Imura Kazuyuki Aihara Graduate
More information6.003 (Fall 2011) Quiz #3 November 16, 2011
6.003 (Fall 2011) Quiz #3 November 16, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 3 1 pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
More informationSTABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential
More information