Lecture Note #6 (Chap.10)
|
|
- Sydney Charles
- 6 years ago
- Views:
Transcription
1 System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryoo Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of models Model reduction Model order reduction in a linear system Balanced realization Pade approximation Moment matching Continued fraction approximation Model approximation of a nonlinear system by a linear one Linearization Describing function analysis Approximation of the nonlinear system by ignoring higher-order harmonics
2 Linearization Nonlinear differential equation: x f( xu, ) Linearization: x f + f ( x, u )( x x ) + f ( x, u )( u u ) Discretization Sampling time: h xt () Fx() t + Gu() t x + Φx Γu + yt () Cx() t y Cx where Fh h Fs Φ e and Γ e Gds Some heuristic model reduction methods (For linear system only) H( z). z /(.7z.8 z ). z /(.8 z )(+. z ) Polynomial truncation H ( z ).3 z /(.7 z ) Method of dominating poles H () z. z /(.8 z ) Pole-zero cancellations 3 (). /(. ) x u H z z + z Cancel (-.8z - ) with (-.78z - ) Balanced Realization and Model Reduction he Reachability Gramian x x u n S: + Φ Γ + x R y Cx What states can be reached with a given input energy assuming that x? N Given finite input energy: Ju ( ) euu uu E States: N N [ N N x ][ u N-] N Φ Γu Φ Φ Γ u ψ NU N Input Sequence for desired x N : U ( ψ ψ ) ψ x N N N N N N N N N N Ju ( ) uu U U x P x E N N N N N where P ψ ψ ΦΓΓ ( Φ ) Bound on the reachable states at time N P N : Reachability Gramian P N ΦPNΦ + ΓΓ + ΦPΦ P + ΓΓ Lyapunov equation
3 he Observability Gramian What states energy is necessary for all u in order to obtain a specified output energy? Specified output energy: J ( y ) e y y E yy Φ Φ (, ) ( ) Q Φ CCΦ Q : Observability Gramian J y y y xq x E y Cx C x C x u where ( ) Q Φ QΦ + CC Φ QΦ Q+ CC Lyapunov equation ( Φ QΦ + CC ( Φ ) + CC Φ + + CC Φ ) Reachability Gramian decides the minimum input energy to drive a system from x to x N. Observability Gramian decides the maximum output energy that the initial state x can generate. Balanced Realization Reachability and observability Gramians P and Q define matrices that describe the sensitivity of the input-output map in different directions of state space. Consider a state-space transformation z x, then the system becomes z+ Φ z + Γu P S : z P y C z Qz Q Different state-space realization may result in different Gramians Balanced realization: P z Q z Pz Qz Σ diag( σ, σ) (where σi λi( PQ)) Algorithm using Cholesyfactorization Hanel singular value Q QQ QPQ UΣ U ( UU I) Σ U Q Σ Σ Σ If some σ i s are relatively small, the Σ can be reduced and the resulting transformation produces abalanced model resuction.
4 Example. (Balanced model reduction) S : H() z. z /(.7z.8 z ) Controllable canonical realization.7.8 x+ x u + y Balanced realization.57.7 P Q Q Σ z z + u (. ) y Σ x From Lyapunov equation ( ) z Reduced model Since first singular value is dominant, second variable can be reduced as if it has no dynamics (z + z )..55 Σ z + Φ Φ z Γ.69 u z Φ Φ + + z Γ z Φ z + Φ z + Γ u z ( I Φ ) ( Φ z + Γ u ) z Φ z + Φ ( I Φ ) ( Φ z + Γ u ) + Γ u + ( Φ + Φ ( I Φ ) Φ ) z + ( Γ + Φ ( I Φ ) Γ ) u + ( Φ) ( Φ + Γ ) y Cz C I z u ( C + C ( I Φ ) Φ ) z + C ( I Φ ) Γ u z.7976z u y.4478z +.95u
5 Balanced model reduction for continuous systems x A A x B x u x A A + x B y ( C C) + Du x x is assumed to have no dynamics x and x A A x A B u his method preserves the essential low-frequency properties (gain). he reduced-order state-space model x ( A A A A ) x + ( B A A B) u y ( C CA A) x + ( D CA B) u Other choices x, x Ax + Bu y Cx + Du x α x x ( A A ( αi A ) A ) x + ( B A ( αi A ) B) u, Comparison Direct transmission term (non-zero even if D) y ( C C ( αi A ) A ) x + ( D C ( αi A ) B ) u x Gred () G() x G ( ) ( ) x red G αx Gred ( α) G( α) Example. (Rohr s system) x x+.338 u y x ( ) Σ /(s+) Reduced second-order model x x+ u y ( ) x+.44u 9 Gs () Reduced first-order model ( s+ ) s + 3s+ 9 x.9686x+.46u.8955s G() s +.44 y.46x.37u s s+.353 G( s) /( s+ )
6 Continued Fraction Approximation Routh array For an asymptotically stable system Bs () Gs () R () s R () s As () c + c3 + c5 + c c4 c6 + R () s + R () s + R3() s s s s With m coefficients, the approximating transfer function is m-th order. Bm () s Gm( s) (with Rm ) Am () s Calculation of coefficients ( ) ( ) As () a + as + P () s Bs () b+ bs + P () s a a a c a / b b c b / p b b () () () p p p c 3 p / p () i ( i ) ( i ) ( i ) ( i ) p p p p / p + + () i () i () i i i p c i+ p / p p p Interpretation of continued fraction approximation he model reduction presuppose that the inner most F bloc may be eliminated. Rm () s he condition for good approximation: c / s R () s m m s For discrete-time system Y ( z) c z U ( z ) + Y ( z) i i i+ i+ U ( z) c Y( z) + U ( z) i+ i i i Yi+ ( z) + ci ciz ciz Yi() z U ( z) c U () z i+ i i Yi( z) c iz Yi+ () z Ui( z) c U i () z i + c i c iz + (Used in lattice algorithms) Unimodular (det)
7 Example Gs () s + 3s + 59s s +.354s +.44s 3 3 Routh array c c c c.9 c.44 c Reduced-order models G () s.5 + s/ s G() s.5+ /(.7683/ s+ /( s/.473)) /(s+).65s +.985s +.986s Moment Matching ransfer function as infinite series H( z ) h + hz + hz + Matching of the reduced-order model, B m /A m H m m m m m Am( z ) + az + + amz M nh n ( z ) Moment: n Moment matching M M M H m z dh dz m B ( z ) b + bz + + bz z dhm dhm + dz dz z z (,,3, )
8 Example.4 (Moment matching).z 3 H( z ).z.54z.54z.7z.8z Reduced-order model Bm ( z ) b bz H m( z ) Am ( z ) + + az Moment matching b + b Hm z + a m z ( + a ) m 3 ( + a ) z Reduced first-order model dh ab + b 4.99 dz dh ( ab ab) dz z H( z ).7993z Pade Approximation runcated aylor series expansion Bs () Gs () Gm() s g + gs + + gm s As () Pade approximation B () s G () sa () s m m m Example.5 First order: m Ys () 458 Gs ().6s.876s.898s U() s s + 3s + 59s B ()/ s A() s b /( + as) b /( + as ).6 s B()/ s A() s /(+.3 s) Second order: B ()/ s A () s ( b + bs)/( + as+ as) ( b + bs)/( + as+ as ).6s+.876s.898s 3 B s A s s + s+ s ()/ () (.645 )/( )
9 Example.6 (Unstable Pade approximation) B () s (+ ) s A () s ( + s+ s ) b ( ) 3 + s s + s + ( g gs ) b g ( ag g) s ags +as B() s b b g B() s A() s ( + as ) a g/ g A() s ( s) Unstable he Pade approximation can be very poor. Describing Function Analysis Analysis of Linear systems Laplace and Fourier transforms Describing function analysis Extension of frequency analysis to nonlinear systems Based on harmonic analysis Assumptions Assumed periodic solution is sufficiently close to sinusoidal oscillation y(t)csin(ωt). he nonlinearity can be represented by un(x,x ). Fourier series expansion f ( x i x ) cos( ) sin( ) a + a x + a x ce / / where a ()cos( ), ()sin( ) / f t t dt b f t t dt / ( +ϕ ) c a + b, ϕ arctan( a / b )
10 Description of periodic oscillation with input amplitude C c( C) iϕ b + ia N( C) e Describing function C C / where a nc ( sin tc, cos t )cos( t) dt, ω ω ω / / b nc ( sin ωtc, ωcos ωt)sin( t) dt / Describing function N(C) is obtained as the amplitude dependent gain N(C) and its phase shift ϕ (ω) for the nonlinear elements. Harmonics caused by nonlinearities are ignored in describing function analysis Characteristic equation Bs ( ) G () s Apxt ( ) () + Bpnxx ( ) (, ), p d/ dt As ( ) ( + G ( pnc ) ( )) xt () G ( p) / N( C) If the Nyquist curves of G (iω) and /N(C) cross, it indicates the possible existence of a limit cycle with amplitude C and frequency ω. Example.7 (A rate-limited servo) (tan S S / ), + S C C> S NC ( ) π C C, C S For small C, /N(C) (max.) Phase angle of /N(C) is π G (.44 i ).667 he Nyquist contour and the describing function cross when K>6. his implies that the limit cycle appears when K>6. Saturatingamplifier K K6 K
11 Condition of application Assumption:the input to nonlinearity is a sinusoid Linear element is usually low-pass: higher harmonics will be attenuated he conditions G ( iω) G ( iω) for,3, and G ( iω ) as G (s) must not have any imaginary poles si. he function n(x,x ) should have finite partial derivatives w.r.t. x and x and must not be an explicit function of time. he zero-order coefficient c of n(x,x ). If the prerequisites are not satisfied, the describing function analysis may predict oscillations that do not exist and may fail to predict periodic solutions that indeed do exist. Example.8 (Pulse-width modulation, PWM) PWM leads to asymmetric inputs It requires more complicated analysis with a nonzero constant term c of the describing function. Consider a PWM system with switching nonlinearity he pulse-width modulating signal z is often a sawtooth-shaped signal of high frequency. he pulse-width modulated signal is a square wave with nonzero mean. ut () nvv (,) (/ z)() xt + nz If the modulating frequency ω z is high compare to G (s), the high frequency components deriving from z are efficiently absorbed by the low-pass lin G (s). ut () (/ z )() xt (linearized form) (n is the ordinary describing function terms)
12 Balanced Model Reduction in Identification Usage of balanced model reduction Further reduction from alinearized model Identify a reduced-order linear model Identification assisted with model reduction Estimation of high-order model (with sufficient excitation) Model reduction with elimination of less important states Example.9 (Model reduction in identification) S: y +.9y+. u + d ( d ) Identification as first-order model: biased y+.9838y +.638u Estimate -th order model, then reduce to first-order model x+.9466x +.9u.9z +.9 Y( z) U() z y.9x.9u z.9466 Example. (An impulse response test) Impulse response data t t gt (). (.5) +. (.75) State-space model for impulse response, g x + x + u Balanced model reduction with n t t gt ().38 (.53) +.38 (.73) ( ) y g g g g x n n Impulse response coefficients
Lecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
More informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationIntro. Computer Control Systems: F8
Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2
More informationNeed for transformation?
Z-TRANSFORM In today s class Z-transform Unilateral Z-transform Bilateral Z-transform Region of Convergence Inverse Z-transform Power Series method Partial Fraction method Solution of difference equations
More informationModel structure. Lecture Note #3 (Chap.6) Identification of time series model. ARMAX Models and Difference Equations
System Modeling and Identification Lecture ote #3 (Chap.6) CHBE 70 Korea University Prof. Dae Ryoo Yang Model structure ime series Multivariable time series x [ ] x x xm Multidimensional time series (temporal+spatial)
More informationSystem Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang
System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang 1-1 Course Description Emphases Delivering concepts and Practice Programming Identification Methods using Matlab Class
More informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationEE 380. Linear Control Systems. Lecture 10
EE 380 Linear Control Systems Lecture 10 Professor Jeffrey Schiano Department of Electrical Engineering Lecture 10. 1 Lecture 10 Topics Stability Definitions Methods for Determining Stability Lecture 10.
More informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationBasic Procedures for Common Problems
Basic Procedures for Common Problems ECHE 550, Fall 2002 Steady State Multivariable Modeling and Control 1 Determine what variables are available to manipulate (inputs, u) and what variables are available
More informationKalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and
Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1
More informationLecture Note #7 (Chap.11)
System Modeling and Identification Lecture Note #7 (Chap.) CBE 702 Korea University Prof. Dae Ryoo Yang Chap. Real-time Identification Real-time identification Supervision and tracing of time varying parameters
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 informationCHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 5-1 Road Map of the Lecture V Laplace Transform and Transfer
More informationFourier Series. Spectral Analysis of Periodic Signals
Fourier Series. Spectral Analysis of Periodic Signals he response of continuous-time linear invariant systems to the complex exponential with unitary magnitude response of a continuous-time LI system at
More information( ) ( = ) = ( ) ( ) ( )
( ) Vρ C st s T t 0 wc Ti s T s Q s (8) K T ( s) Q ( s) + Ti ( s) (0) τs+ τs+ V ρ K and τ wc w T (s)g (s)q (s) + G (s)t(s) i G and G are transfer functions and independent of the inputs, Q and T i. Note
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
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 informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More 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 informationChemical Process Dynamics and Control. Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University
Chemical Process Dynamics and Control Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University 1 Chapter 4 System Stability 2 Chapter Objectives End of this
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationMEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction
MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214 Outline Solving State Equations Variation
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 informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
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 informationPrüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid
More informationControl Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationECEEN 5448 Fall 2011 Homework #5 Solutions
ECEEN 5448 Fall 211 Homework #5 Solutions Professor David G. Meyer December 8, 211 1. Consider the 1-dimensional time-varying linear system ẋ t (u x) (a) Find the state-transition matrix, Φ(t, τ). Here
More informationOutput Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)
Output Feedback and State Feedback EL2620 Nonlinear Control Lecture 10 Exact feedback linearization Input-output linearization Lyapunov-based control design methods ẋ = f(x,u) y = h(x) Output feedback:
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 22: Balanced Realization Readings: DDV, Chapter 26 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology April 27, 2011 E. Frazzoli
More informationControl Systems. Frequency domain analysis. L. Lanari
Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic
More informationTransfer Functions. Chapter Introduction. 6.2 The Transfer Function
Chapter 6 Transfer Functions As a matter of idle curiosity, I once counted to find out what the order of the set of equations in an amplifier I had just designed would have been, if I had worked with the
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 informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
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 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 informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationExam in Systems Engineering/Process Control
Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
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 informationCourse roadmap. ME451: Control Systems. Example of Laplace transform. Lecture 2 Laplace transform. Laplace transform
ME45: Control Systems Lecture 2 Prof. Jongeun Choi Department of Mechanical Engineering Michigan State University Modeling Transfer function Models for systems electrical mechanical electromechanical Block
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time
More informationUnit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace
Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,
More informationControl Systems. Dynamic response in the time domain. L. Lanari
Control Systems Dynamic response in the time domain L. Lanari outline A diagonalizable - real eigenvalues (aperiodic natural modes) - complex conjugate eigenvalues (pseudoperiodic natural modes) - phase
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationControl Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationECE557 Systems Control
ECE557 Systems Control Bruce Francis Course notes, Version.0, September 008 Preface This is the second Engineering Science course on control. It assumes ECE56 as a prerequisite. If you didn t take ECE56,
More informationLecture 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster.
Lecture 8 Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture
More informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationChapter 3 Data Acquisition and Manipulation
1 Chapter 3 Data Acquisition and Manipulation In this chapter we introduce z transf orm, or the discrete Laplace Transform, to solve linear recursions. Section 3.1 z-transform Given a data stream x {x
More informationEE Control Systems LECTURE 9
Updated: Sunday, February, 999 EE - Control Systems LECTURE 9 Copyright FL Lewis 998 All rights reserved STABILITY OF LINEAR SYSTEMS We discuss the stability of input/output systems and of state-space
More informationECE504: Lecture 9. D. Richard Brown III. Worcester Polytechnic Institute. 04-Nov-2008
ECE504: Lecture 9 D. Richard Brown III Worcester Polytechnic Institute 04-Nov-2008 Worcester Polytechnic Institute D. Richard Brown III 04-Nov-2008 1 / 38 Lecture 9 Major Topics ECE504: Lecture 9 We are
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More information2.004 Dynamics and Control II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.004 Dynamics and Control II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8 I * * Massachusetts
More informationAssignment 11 Assigned Mon Sept 27
Assignment Assigned Mon Sept 7 Section 7., Problem 4. x sin x dx = x cos x + x cos x dx ( = x cos x + x sin x ) sin x dx u = x dv = sin x dx du = x dx v = cos x u = x dv = cos x dx du = dx v = sin x =
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 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 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 informationSolutions to Exercises 6.1
34 Chapter 6 Conformal Mappings Solutions to Exercises 6.. An analytic function fz is conformal where f z. If fz = z + e z, then f z =e z z + z. We have f z = z z += z =. Thus f is conformal at all z.
More informationMAE143a: Signals & Systems (& Control) Final Exam (2011) solutions
MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions Question 1. SIGNALS: Design of a noise-cancelling headphone system. 1a. Based on the low-pass filter given, design a high-pass filter,
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationSTABILITY OF CLOSED-LOOP CONTOL SYSTEMS
CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 10-1 Road Map of the Lecture X Stability of closed-loop control
More informationMechatronics Engineering. Li Wen
Mechatronics Engineering Li Wen Bio-inspired robot-dc motor drive Unstable system Mirko Kovac,EPFL Modeling and simulation of the control system Problems 1. Why we establish mathematical model of the control
More informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
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 informationCBE507 LECTURE III Controller Design Using State-space Methods. Professor Dae Ryook Yang
CBE507 LECTURE III Controller Design Using State-space Methods Professor Dae Ryook Yang Fall 2013 Dept. of Chemical and Biological Engineering Korea University Korea University III -1 Overview States What
More information2006 Fall. G(s) y = Cx + Du
1 Class Handout: Chapter 7 Frequency Domain Analysis of Feedback Systems 2006 Fall Frequency domain analysis of a dynamic system is very useful because it provides much physical insight, has graphical
More informationPoles and Zeros in z-plane
M58 Mixed Signal Processors page of 6 Poles and Zeros in z-plane z-plane Response of discrete-time system (i.e. digital filter at a particular frequency ω is determined by the distance between its poles
More information9.5 The Transfer Function
Lecture Notes on Control Systems/D. Ghose/2012 0 9.5 The Transfer Function Consider the n-th order linear, time-invariant dynamical system. dy a 0 y + a 1 dt + a d 2 y 2 dt + + a d n y 2 n dt b du 0u +
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
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 informationECE317 : Feedback and Control
ECE317 : Feedback and Control Lecture : Routh-Hurwitz stability criterion Examples Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationModule 4. Related web links and videos. 1. FT and ZT
Module 4 Laplace transforms, ROC, rational systems, Z transform, properties of LT and ZT, rational functions, system properties from ROC, inverse transforms Related web links and videos Sl no Web link
More informationMATH 452. SAMPLE 3 SOLUTIONS May 3, (10 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic.
MATH 45 SAMPLE 3 SOLUTIONS May 3, 06. (0 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic. Because f is holomorphic, u and v satisfy the Cauchy-Riemann equations:
More informationChapter 6: The Laplace Transform. Chih-Wei Liu
Chapter 6: The Laplace Transform Chih-Wei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace
More informationLOWELL WEEKLY JOURNAL. ^Jberxy and (Jmott Oao M d Ccmsparftble. %m >ai ruv GEEAT INDUSTRIES
? (») /»» 9 F ( ) / ) /»F»»»»»# F??»»» Q ( ( »»» < 3»» /» > > } > Q ( Q > Z F 5
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationAero III/IV Conformal Mapping
Aero III/IV Conformal Mapping View complex function as a mapping Unlike a real function, a complex function w = f(z) cannot be represented by a curve. Instead it is useful to view it as a mapping. Write
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationStability. X(s) Y(s) = (s + 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = +2 and s = -2. Y(s) 8X(s) G 1 G 2
Stability 8X(s) X(s) Y(s) = (s 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = 2 and s = -2 If all poles are in region where s < 0, system is stable in Fourier language s = jω G 0 - x3 x7 Y(s)
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationDiscrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture
Discrete Systems Mark Cannon Hilary Term 22 - Lecture 4 Step response and pole locations 4 - Review Definition of -transform: U() = Z{u k } = u k k k= Discrete transfer function: Y () U() = G() = Z{g k},
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationMath 211. Substitute Lecture. November 20, 2000
1 Math 211 Substitute Lecture November 20, 2000 2 Solutions to y + py + qy =0. Look for exponential solutions y(t) =e λt. Characteristic equation: λ 2 + pλ + q =0. Characteristic polynomial: λ 2 + pλ +
More informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More informationCHAPTER 3. Analytic Functions. Dr. Pulak Sahoo
CHAPTER 3 Analytic Functions BY Dr. Pulak Sahoo Assistant Professor Department of Mathematics University Of Kalyani West Bengal, India E-mail : sahoopulak1@gmail.com 1 Module-4: Harmonic Functions 1 Introduction
More informationLinear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output
Linear Systems Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output Our interest is in dynamic systems Dynamic system means a system with memory of course including
More information