AA/EE/ME 548: Problem Session Notes #5

Size: px
Start display at page:

Download "AA/EE/ME 548: Problem Session Notes #5"

Transcription

1 AA/EE/ME 548: Problem Session Notes #5 Review of Nyquist and Bode Plots. Nyquist Stability Criterion. LQG/LTR Method Tuesday, March 2, 203 Outline:. A review of Bode plots. 2. A review of Nyquist plots and Nyquist stability criterion. 3. LQG/LTR method - example. A Review of Bode Plots The frequency response of a linear system can be computed from its transfer function, by setting the input u(t) to: u(t) = exp{jωt} = cos(ωt)+jsin(ωt) () The resulting output can now be found as: y(t) = G(jω)exp{jωt} = M exp{j(ωt+ϕ)} = M cos(ωt+ϕ)+jm sin(ωt+φ) (2) where G(jω) denotes the transfer function, and M and ϕ are its gain and phase: M = G(jω) ϕ = arctan Im{G(jω)} Re{G(jω)} (3) The frequency response G(jω) can thus be represented by two curves: the gain curve and the phase curve. The gain curve represents G(jω) as a function of frequency ω, and the phase curve the ϕ as a function of frequency. 2 Nyquist Stability Criterion The Nyquist stability criterion allows us to determine if a system is stable or unstable. In addition, it also provides us with a measure of the degree of stability through the definition of stability margins. The given stability criterion also indicates how an unstable system should be changed to make it stable. Before stating the Nyquist stability criterion, let s first recall the definition of loop-transfer function L(s). 2. The Loop-Transfer Function Consider a block diagram, depicted in Figure, which represents the stochastic dynamical system controlled by a dynamical LQG controller. The loop-transfer function of the given closed-loop system can be found as: L(s) := F(s)G(s) (4) and it is obtained by breaking the feedback loop between the controller and the system. Thus, the loop transfer function is simply the transfer function from input u(t) to output y(t).

2 r(t) = LQG Controller F(s) u(t) System G(s) = C(sI - A) - B y(t) Fig.. Block diagram of a dynamic regulator, consisting of an optimal Kalman filter and an optimal LQR statefeedback controller. 2.2 The Nyquist Plot As stated in [2], the Nyquist plot of the loop-transfer function L(s) is formed by tracing s C around the Nyquist D contour,, denotes as Γ C. This contour consists of the imaginary axis combined with an arc at infinity connecting the endpoints of the imaginary axis. The image of L(s), when s traverses Γ, gives a closed curve in the complex plane an that curve is referred to as the Nyquist plot for L(s). There is a bit of subtlety in the Nyquist plot when the loop-transfer function has poles on the imaginary axis because the gain is infinite at the poles. To solve this problem, contour Γ C is modified to include small deviations that avoid any poles on the imaginary axis. The deviation consists of a small semicircles to the right of the imaginary axis pole location. Fig. 2. (a) The Nyquist D contour and (b) The Nyquist plot. The Nyquist contour encloses the right half-plane, with a small semicircle around any poles of L(s) on the imaginary axis (here illustrated at the origin) an an arc at infinity, represented by R. The Nyquist plot is the image of the loop-transfer function L(s) when s traverses the Nyquist contour Γ in the clockwise direction. The solid line corresponds to ω > 0 and the dashed line to ω < 0. The plots are reproduced from [2], and the curve is generated for L(s) =.4 exp{ s} (s+) The Nyquist Stability Criterion We will state two versions of the Nyquist stability criteria. The first one represents the special case when loop-transfer function L(s) is assumed to be stable. Theorem. (Simplified Nyquist criterion) Let L(s) be the loop-transfer function for a negative feedback system (as shown in Figure ) and assume that L has no unstable poles, except for single poles on the imaginary axis. Then the closed-loop system is stable if and only if the closed contour given by Ω = {L(jω) : < ω < } C has no net encirclements of the critical point s =. Theorem requires that loop-transfer function L(s) has no unstable poles. In some cases, however, this requirement is not satisfied, and a more general result is needed. Nyquist originally considered this general case, which can be summarized as Theorem 2. 2

3 Theorem 2. (Nyquist s Stability Theorem) Consider a closed-loop system with loop-transfer function L(s) that has P poles in the region enclosed by the Nyquist contour. Let N be the net number of clockwise encirclements of critical point - by L(s) when s encircles the Nyquist contour Γ C in the clockwise direction. The closed-loop system then has Z = N +P poles in the right half-plane. Note : Theorem 2 states that if L(s) has P unstable poles, then the Nyquist curve for L(s) should have P counterclockwise encirclements of the critical point - (so that N = P). In particular, this requires that the magnitude L(jω c ) > for some frequency ω c, corresponding to the crossing of the negative real axis. Note 2: As pointed out in [2], care has to be taken to get the right sign of the encirclements. The Nyquist contour has to be traversed clockwise, which means that frequency ω moves from to and N is positive if the Nyquist curve winds clockwise. If the Nyquist curve winds counterclockwise, then N will be negative (the desired case if P 0). 2.4 Stability Margins In real applications, it may not be enough to simply determine whether a system is stable or not. Very often we will also want to have some margins of stability, to describe how stable the system is and how robust it is to perturbations. One of the common ways to express these stability margins is the use of gain and phase margins, inspired by Nyquists stability criterion. It is known that an increase in controller gain simply expands the Nyquist plot radially. Similarly, an increase in the phase of the controller twists the Nyquist plot. It therefore becomes easy to read the amount of gain and phase that can be added without causing the system to become unstable from the Nyquist plot. Definition : The gain margin g m of a system is defined as the smallest amount that the open-loop gain can be increased before the closed-loop system goes unstable. For a system whose phase decreases monotonically as a function of frequency starting at 0 degrees, the gain margin can be computed based on the smallest frequency where the phase of the loop-transfer function L(s) equals to -80 degrees. Let ω pc represent the frequency at which that happens. We refer to ω pc the phase crossover frequency. Using ω pc, the gain margin can be computed as: g m = L(jω pc ) (5) Definition 2: The phase margin is defined as the amount of phase lag required to reach the stability limit. Let ω gc denote the gain crossover frequency, the smallest frequency where the loop-transferfunction L(s) has unit magnitude. Then for a system with monotonically decreasing gain, the phase margin is defined as: ϕ m = π +argl(jω gc ) (6) Note 3: The gain and phase margins have simple geometric interpretations on the Nyquist diagram of the loop-transfer function. The gain margin is given by the inverse of the distance to the nearest point between - and 0, where the loop transfer function crosses the negative real axis. The phase margin is given by the smallest angle on the unit circle between - and the loop-transfer function. When the gain or phase is monotonic, this geometric interpretation agrees with the formulas above. 2.5 Stability Margins from Bode Plot Formany systems,the gainand phaseand marginscanbe determined fromthe Bode plot ofthe loop-transfer function. 3

4 To find the gain margin, we first find the phase crossover frequency ω pc where the phase is -80 degrees. The gain margin is the inverse of the gain at that frequency. To determine the phase margin, we first determine the gain crossover frequency ω gc, i.e., the frequency where the gain of the loop-transfer function is. The phase margin is the phase of the loop-transfer function at that frequency plus 80 degrees. Fig.3. Stability margins. The gain margin g m and phase margin ϕ m are shown on the Nyquist plot (a) and the Bode plot (b). The gain margin corresponds to the smallest increase in gain that creates an encirclement, and the phase margin is the smallest change in phase that creates an encirclement. Figures reproduced from [2]. 2.6 Example Consider the following stochastic dynamical system: [ ] ẋ = x+ 0 [ 0 ] u+ [ ] 0 w 0 y = [ 0 ] x+v (7) with noise covariance matrices Q n and R n given as: [ ] Q n = ρ, R n = (8) If our goal is to solve the LQG problem with cost matrices Q and R: [ ] Q = σ, R = (9) we could compute the optimal steady-state control gain K as: K = α [ ], where α := 2+ 4+ρ (0) and the optimal steady-state Kalman gain L as: [ L = β, where β := 2+ ] 4+σ () Please find the Nyquist plot of the given closed-loop system if α = 8 and β = 27. Use Matlab function nyquist(), if necessary. What are the gain and the phase margins in this case? 4

5 Solution: To solve the given problem, we observe the gains are given as: K = 8 [ ] ] L = 27 ] (2) Using equations: α := 2+ 4+ρ β := 2+ 4+σ (3) we can find ρ and σ as: ρ = 252 σ = 62 (4) The loop-transfer function of the closed-loop system using the optimal steady-state LQG can now be found as: H LQG (s) = F(s)G(s) (5) where F(s) denotes the transfer function of the LQG controller and G(s) the transfer function of the openloop system: F(s) = K(sI A+BK +LC) L G(s) = C(sI A) B (6) The gain and phase margins can now be found using Matlab functions nyquist() and margin(). The complete code of this example is given below. %Example %%System dynamics A = [, ; 0 ]; B = [0; ]; 6 C = [, 0]; D = [0]; %%Cost matrices rho = 252 Q = rho [;] [,]; R = []; %%Noise covariance matrices sigma = 62; 6 Qn = sigma [;] [,]; Rn = []; %%Controller and observer gains alpha = 8; 2 beta = 27; Kc = alpha [,]; Kf = beta [;]; %%Open loop transfer function 26 [num P, den P] = ss2tf (A,B,C,D) ; %%Dynamical controller W = [Q, zeros (2, ) ; zeros (, 2), R]; 5

6 V = [Qn, zeros (2, ) ; zeros (, 2), Rn]; 3 [Af,Bf,Cf,Df] = lqg(a,b,c,d,w,v) ; %%Controller s transfer function [num F, den F] = ss2tf (Af, Bf, Cf, Df) ; G = tf (num P, den P) 36 F = tf (num F, den F) 4 %%Nyquist plot and margins figure () nyquist (F G) ; [Gm,Pm, Wcg, Wcp] = margin(f G) figure () margin(f G) ; 3 LQG/LTR Method - Matlab Software and Example 3. Stability Margins and Minimum Singular Values Let s recall the minimum singular value of I +H LQG (s) can be found as: inf ω σ(i +H LQG(jω)) := α (7) It can be shown the following inequality relations between gain margin, g m, phase margins, ϕ m, and the minimum singular value, α, are satisfied: 3.2 Matlab Software +α < g m < α ( ) ) cos α2 < ϕ m < cos ( α2 2 2 As shown in [3], Matlab functions ltru() and ltry() from the Robust Control Toolbox can be used to do LQG/LTR design. The first function (ltru()) does loop recovery at the input, with a fixed state-feedback matrix K, whereas the latter function does loop recovery at the output, with fixed Kalman gain L. The command for ltru() is given as: [af, bf, cf, df, svl] = ltru(sys, dim, Kc, Xi, Th, r, w, svk) where: [ ] A B Input sys defines the system to be controlled, sys :=, C D Input dim defines the size of matrix A, Input Kc is the full-state feedback matrix, Input Xi is the nominal process noise covariance matrix Q n, Input Th is the measurement noise covariance matrix R n, Input r is a row vector of increasing fictitious-noise coefficients Input w is a row vector of frequencies for Bode plots of the singular values of the loop gain Input svk is the singular-value Bode plot of H LQR (s) and it may be omitted The function ltru() producesbodeplots ofthe singularvaluesofthe loop-transferfunction H LQG (s) foreach of the values in vector r and returns a controller F(s), specified by a state-space realization (af,bf,cf,df), for the last value of r specified in the vector r. The singular values of H LQG (s) are stored in the output data svl. (8) 6

7 3.3 Example Let s consider the following helicopter model: ẋ = x u [ y = where the incremental outputs denote: y - the vertical velocity (knots/hr) y 2 - the pitch altitude (radians) and the inputs are: u - the collective rotor thrust u 2 - the differential collective rotor thrust Given cost matrices Q and R: and noise covariance matrices: w ] x (9) Q = C T C, R = I 2 2 (20) Q n = BB T, R n = I 2 2 (2) the minimum cost matrix S and the minimum error covariance matrix P were found as: P = S = (22) and the state-space representation of the optimal LQG controller is: ] ˆx = ˆx u [ y [ ] u opt = ˆx (23) The optimal loop-transfer function H LQG (s) can now be computed as: H LQG (s) = F(s)C(sI A) B (24) and the closed-loop state-space representation can be found as: [ẋ ] [ [ = x ˆx u x] 0 0 y] [ ] y = (25) ][ xˆx 7

8 The state-space representation is then used to plot the singular values of I + H LOG (s) versus frequency (Matlab function sigma()). The plot of the singular values if depicted in Figure 4, from which the minimum singular value is determined to be α = From this, we can obtain the bounds on the gain margin in each input channel using the following inequality: +α < g m < α (26) For the given system, we thus obtain: Consider now the problem of recovering gain margins: < l i < (27) < α < (28) Using inequality (26), that corresponds to α = Using Matlab function ltru(), the value of fictitious noise that corresponds to α = 0.65 was found to be r = 00. Fig.4. Singular values of I+H LQG(S) as functions of frequency. The x-axis represents frequency and y-axis singular values. Figure taken from [3]. Sources for Today s Lecture:. Karl J. Aestrom, Richard M. Murray Feedback Systems - An Introduction for Scientists and Engineers. Princeton University Press, 2009, Chapter 8, p Karl J. Aestrom, Richard M. Murray Feedback Systems - An Introduction for Scientists and Engineers. Princeton University Press, 2009, Chapter 9, p Peter Dorato, Chaouki T. Abdullah, Vito Cerone Linear Quadratic Control: An Introduction. Krieger Publishing, 2000 Chapter 7, p

Linear Control Systems Lecture #3 - Frequency Domain Analysis. Guillaume Drion Academic year

Linear Control Systems Lecture #3 - Frequency Domain Analysis. Guillaume Drion Academic year Linear Control Systems Lecture #3 - Frequency Domain Analysis Guillaume Drion Academic year 2018-2019 1 Goal and Outline Goal: To be able to analyze the stability and robustness of a closed-loop system

More information

x(t) = x(t h), x(t) 2 R ), where is the time delay, the transfer function for such a e s Figure 1: Simple Time Delay Block Diagram e i! =1 \e i!t =!

x(t) = x(t h), x(t) 2 R ), where is the time delay, the transfer function for such a e s Figure 1: Simple Time Delay Block Diagram e i! =1 \e i!t =! 1 Time-Delay Systems 1.1 Introduction Recitation Notes: Time Delays and Nyquist Plots Review In control systems a challenging area is operating in the presence of delays. Delays can be attributed to acquiring

More information

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods

More information

The Nyquist criterion relates the stability of a closed system to the open-loop frequency response and open loop pole location.

The Nyquist criterion relates the stability of a closed system to the open-loop frequency response and open loop pole location. Introduction to the Nyquist criterion The Nyquist criterion relates the stability of a closed system to the open-loop frequency response and open loop pole location. Mapping. If we take a complex number

More information

Control Systems I. Lecture 9: The Nyquist condition

Control Systems I. Lecture 9: The Nyquist condition Control Systems I Lecture 9: The Nyquist condition Readings: Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Jacopo Tani Institute for Dynamic Systems and Control

More information

Control Systems. Frequency Method Nyquist Analysis.

Control Systems. Frequency Method Nyquist Analysis. Frequency Method Nyquist Analysis chibum@seoultech.ac.kr Outline Polar plots Nyquist plots Factors of polar plots PolarNyquist Plots Polar plot: he locus of the magnitude of ω vs. the phase of ω on polar

More information

Control Systems I. Lecture 9: The Nyquist condition

Control Systems I. Lecture 9: The Nyquist condition Control Systems I Lecture 9: The Nyquist condition adings: Guzzella, Chapter 9.4 6 Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Emilio Frazzoli Institute

More information

Return Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems

Return Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems Spectral Properties of Linear- Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018! Stability margins of single-input/singleoutput (SISO) systems! Characterizations

More information

Systems Analysis and Control

Systems 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 information

Topic # Feedback Control Systems

Topic # Feedback Control Systems Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the

More information

Classify a transfer function to see which order or ramp it can follow and with which expected error.

Classify a transfer function to see which order or ramp it can follow and with which expected error. Dr. J. Tani, Prof. Dr. E. Frazzoli 5-059-00 Control Systems I (Autumn 208) Exercise Set 0 Topic: Specifications for Feedback Systems Discussion: 30.. 208 Learning objectives: The student can grizzi@ethz.ch,

More information

Systems Analysis and Control

Systems 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 information

Course Outline. Closed Loop Stability. Stability. Amme 3500 : System Dynamics & Control. Nyquist Stability. Dr. Dunant Halim

Course Outline. Closed Loop Stability. Stability. Amme 3500 : System Dynamics & Control. Nyquist Stability. Dr. Dunant Halim Amme 3 : System Dynamics & Control Nyquist Stability Dr. Dunant Halim Course Outline Week Date Content Assignment Notes 1 5 Mar Introduction 2 12 Mar Frequency Domain Modelling 3 19 Mar System Response

More information

Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31

Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31 Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured

More information

1 (s + 3)(s + 2)(s + a) G(s) = C(s) = K P + K I

1 (s + 3)(s + 2)(s + a) G(s) = C(s) = K P + K I MAE 43B Linear Control Prof. M. Krstic FINAL June 9, Problem. ( points) Consider a plant in feedback with the PI controller G(s) = (s + 3)(s + )(s + a) C(s) = K P + K I s. (a) (4 points) For a given constant

More information

MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Intro Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /5/27 Outline Closed Loop Transfer

More information

Autonomous Mobile Robot Design

Autonomous Mobile Robot Design Autonomous Mobile Robot Design Topic: Guidance and Control Introduction and PID Loops Dr. Kostas Alexis (CSE) Autonomous Robot Challenges How do I control where to go? Autonomous Mobile Robot Design Topic:

More information

Topic # Feedback Control Systems

Topic # Feedback Control Systems Topic #19 16.31 Feedback Control Systems Stengel Chapter 6 Question: how well do the large gain and phase margins discussed for LQR map over to DOFB using LQR and LQE (called LQG)? Fall 2010 16.30/31 19

More information

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Lecture 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 information

Intro to Frequency Domain Design

Intro to Frequency Domain Design Intro to Frequency Domain Design MEM 355 Performance Enhancement of Dynamical Systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Closed Loop Transfer Functions

More information

Controls Problems for Qualifying Exam - Spring 2014

Controls 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 information

H(s) = s. a 2. H eq (z) = z z. G(s) a 2. G(s) A B. s 2 s(s + a) 2 s(s a) G(s) 1 a 1 a. } = (z s 1)( z. e ) ) (z. (z 1)(z e at )(z e at )

H(s) = s. a 2. H eq (z) = z z. G(s) a 2. G(s) A B. s 2 s(s + a) 2 s(s a) G(s) 1 a 1 a. } = (z s 1)( z. e ) ) (z. (z 1)(z e at )(z e at ) .7 Quiz Solutions Problem : a H(s) = s a a) Calculate the zero order hold equivalent H eq (z). H eq (z) = z z G(s) Z{ } s G(s) a Z{ } = Z{ s s(s a ) } G(s) A B Z{ } = Z{ + } s s(s + a) s(s a) G(s) a a

More information

Analysis of Discrete-Time Systems

Analysis 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 information

ROOT LOCUS. Consider the system. Root locus presents the poles of the closed-loop system when the gain K changes from 0 to. H(s) H ( s) = ( s)

ROOT LOCUS. Consider the system. Root locus presents the poles of the closed-loop system when the gain K changes from 0 to. H(s) H ( s) = ( s) C1 ROOT LOCUS Consider the system R(s) E(s) C(s) + K G(s) - H(s) C(s) R(s) = K G(s) 1 + K G(s) H(s) Root locus presents the poles of the closed-loop system when the gain K changes from 0 to 1+ K G ( s)

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect

More information

LQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin

LQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system

More information

Systems Analysis and Control

Systems 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 information

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities. 19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that

More information

Topic # Feedback Control

Topic # Feedback Control Topic #4 16.31 Feedback Control Stability in the Frequency Domain Nyquist Stability Theorem Examples Appendix (details) This is the basis of future robustness tests. Fall 2007 16.31 4 2 Frequency Stability

More information

Digital Control Systems

Digital Control Systems Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist

More information

Analysis of Discrete-Time Systems

Analysis 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 information

Control Systems I Lecture 10: System Specifications

Control Systems I Lecture 10: System Specifications Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture

More information

MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/8/25 Outline Closed Loop Transfer Functions

More information

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design

More information

FEL3210 Multivariable Feedback Control

FEL3210 Multivariable Feedback Control FEL3210 Multivariable Feedback Control Lecture 5: Uncertainty and Robustness in SISO Systems [Ch.7-(8)] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 5:Uncertainty and Robustness () FEL3210 MIMO

More information

CDS 101/110a: Lecture 8-1 Frequency Domain Design

CDS 101/110a: Lecture 8-1 Frequency Domain Design CDS 11/11a: Lecture 8-1 Frequency Domain Design Richard M. Murray 17 November 28 Goals: Describe canonical control design problem and standard performance measures Show how to use loop shaping to achieve

More information

Design Methods for Control Systems

Design Methods for Control Systems Design Methods for Control Systems Maarten Steinbuch TU/e Gjerrit Meinsma UT Dutch Institute of Systems and Control Winter term 2002-2003 Schedule November 25 MSt December 2 MSt Homework # 1 December 9

More information

OPTIMAL CONTROL AND ESTIMATION

OPTIMAL CONTROL AND ESTIMATION OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION

More information

MAS107 Control Theory Exam Solutions 2008

MAS107 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 information

Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

3 Stabilization of MIMO Feedback Systems

3 Stabilization of MIMO Feedback Systems 3 Stabilization of MIMO Feedback Systems 3.1 Notation The sets R and S are as before. We will use the notation M (R) to denote the set of matrices with elements in R. The dimensions are not explicitly

More information

Lecture 6. Chapter 8: Robust Stability and Performance Analysis for MIMO Systems. Eugenio Schuster.

Lecture 6. Chapter 8: Robust Stability and Performance Analysis for MIMO Systems. Eugenio Schuster. Lecture 6 Chapter 8: Robust Stability and Performance Analysis for MIMO Systems Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 6 p. 1/73 6.1 General

More information

Problem Set 4 Solution 1

Problem Set 4 Solution 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 4 Solution Problem 4. For the SISO feedback

More information

CDS 101/110a: Lecture 10-1 Robust Performance

CDS 101/110a: Lecture 10-1 Robust Performance CDS 11/11a: Lecture 1-1 Robust Performance Richard M. Murray 1 December 28 Goals: Describe how to represent uncertainty in process dynamics Describe how to analyze a system in the presence of uncertainty

More information

CONTROL SYSTEM STABILITY. CHARACTERISTIC EQUATION: The overall transfer function for a. where A B X Y are polynomials. Substitution into the TF gives:

CONTROL SYSTEM STABILITY. CHARACTERISTIC EQUATION: The overall transfer function for a. where A B X Y are polynomials. Substitution into the TF gives: CONTROL SYSTEM STABILITY CHARACTERISTIC EQUATION: The overall transfer function for a feedback control system is: TF = G / [1+GH]. The G and H functions can be put into the form: G(S) = A(S) / B(S) H(S)

More information

Subject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)

Subject: Optimal Control Assignment-1 (Related to Lecture notes 1-10) Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must

More information

Richiami di Controlli Automatici

Richiami di Controlli Automatici Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici

More information

K(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

K(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 information

Robust fixed-order H Controller Design for Spectral Models by Convex Optimization

Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Alireza Karimi, Gorka Galdos and Roland Longchamp Abstract A new approach for robust fixed-order H controller design by

More information

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state

More information

MAE 143B - Homework 9

MAE 143B - Homework 9 MAE 143B - Homework 9 7.1 a) We have stable first-order poles at p 1 = 1 and p 2 = 1. For small values of ω, we recover the DC gain K = lim ω G(jω) = 1 1 = 2dB. Having this finite limit, our straight-line

More information

Class 13 Frequency domain analysis

Class 13 Frequency domain analysis Class 13 Frequency domain analysis The frequency response is the output of the system in steady state when the input of the system is sinusoidal Methods of system analysis by the frequency response, as

More information

Lecture 1: Feedback Control Loop

Lecture 1: Feedback Control Loop Lecture : Feedback Control Loop Loop Transfer function The standard feedback control system structure is depicted in Figure. This represend(t) n(t) r(t) e(t) u(t) v(t) η(t) y(t) F (s) C(s) P (s) Figure

More information

Uncertainty and Robustness for SISO Systems

Uncertainty and Robustness for SISO Systems Uncertainty and Robustness for SISO Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Outline Nature of uncertainty (models and signals). Physical sources of model uncertainty. Mathematical

More information

Lecture 11. Frequency Response in Discrete Time Control Systems

Lecture 11. Frequency Response in Discrete Time Control Systems EE42 - Discrete Time Systems Spring 28 Lecturer: Asst. Prof. M. Mert Ankarali Lecture.. Frequency Response in Discrete Time Control Systems Let s assume u[k], y[k], and G(z) represents the input, output,

More information

Stability and Robustness 1

Stability and Robustness 1 Lecture 2 Stability and Robustness This lecture discusses the role of stability in feedback design. The emphasis is notonyes/notestsforstability,butratheronhowtomeasurethedistanceto instability. The small

More information

EL2520 Control Theory and Practice

EL2520 Control Theory and Practice EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller

More information

Lecture 4: Analysis of MIMO Systems

Lecture 4: Analysis of MIMO Systems Lecture 4: Analysis of MIMO Systems Norms The concept of norm will be extremely useful for evaluating signals and systems quantitatively during this course In the following, we will present vector norms

More information

Robust Control 3 The Closed Loop

Robust Control 3 The Closed Loop Robust Control 3 The Closed Loop Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /2/2002 Outline Closed Loop Transfer Functions Traditional Performance Measures Time

More information

Theory of Robust Control

Theory of Robust Control Theory of Robust Control Carsten Scherer Mathematical Systems Theory Department of Mathematics University of Stuttgart Germany Contents 1 Introduction to Basic Concepts 6 1.1 Systems and Signals..............................

More information

Outline. Classical Control. Lecture 1

Outline. Classical Control. Lecture 1 Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction

More information

MAE 143B - Homework 9

MAE 143B - Homework 9 MAE 43B - Homework 9 7.2 2 2 3.8.6.4.2.2 9 8 2 2 3 a) G(s) = (s+)(s+).4.6.8.2.2.4.6.8. Polar plot; red for negative ; no encirclements of, a.s. under unit feedback... 2 2 3. 4 9 2 2 3 h) G(s) = s+ s(s+)..2.4.6.8.2.4

More information

r + - FINAL June 12, 2012 MAE 143B Linear Control Prof. M. Krstic

r + - FINAL June 12, 2012 MAE 143B Linear Control Prof. M. Krstic MAE 43B Linear Control Prof. M. Krstic FINAL June, One sheet of hand-written notes (two pages). Present your reasoning and calculations clearly. Inconsistent etchings will not be graded. Write answers

More information

16.30/31, Fall 2010 Recitation # 2

16.30/31, Fall 2010 Recitation # 2 16.30/31, Fall 2010 Recitation # 2 September 22, 2010 In this recitation, we will consider two problems from Chapter 8 of the Van de Vegte book. R + - E G c (s) G(s) C Figure 1: The standard block diagram

More information

Topic # Feedback Control Systems

Topic # 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 information

Homework 7 - Solutions

Homework 7 - Solutions Homework 7 - Solutions Note: This homework is worth a total of 48 points. 1. Compensators (9 points) For a unity feedback system given below, with G(s) = K s(s + 5)(s + 11) do the following: (c) Find the

More information

LQG/LTR CONTROLLER DESIGN FOR AN AIRCRAFT MODEL

LQG/LTR CONTROLLER DESIGN FOR AN AIRCRAFT MODEL PERIODICA POLYTECHNICA SER. TRANSP. ENG. VOL. 8, NO., PP. 3 4 () LQG/LTR CONTROLLER DESIGN FOR AN AIRCRAFT MODEL Balázs KULCSÁR Department of Control and Transport Automation Budapest University of Technology

More information

STABILITY OF CLOSED-LOOP CONTOL SYSTEMS

STABILITY 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 information

Qualitative Graphical Representation of Nyquist Plots

Qualitative Graphical Representation of Nyquist Plots Qualitative Graphical presentation of Nyquist Plots R. Zanasi a, F. Grossi a,, L. Biagiotti a a Department of Engineering Enzo Ferrari, University of Modena and ggio Emilia, via Pietro Vivarelli 0, 425

More information

Lecture 9. Introduction to Kalman Filtering. Linear Quadratic Gaussian Control (LQG) G. Hovland 2004

Lecture 9. Introduction to Kalman Filtering. Linear Quadratic Gaussian Control (LQG) G. Hovland 2004 MER42 Advanced Control Lecture 9 Introduction to Kalman Filtering Linear Quadratic Gaussian Control (LQG) G. Hovland 24 Announcement No tutorials on hursday mornings 8-9am I will be present in all practical

More information

Robust Multivariable Control

Robust 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 information

EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.

EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators. Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3

More information

Nyquist Stability Criteria

Nyquist Stability Criteria Nyquist Stability Criteria Dr. Bishakh Bhattacharya h Professor, Department of Mechanical Engineering IIT Kanpur Joint Initiative of IITs and IISc - Funded by MHRD This Lecture Contains Introduction to

More information

EE3CL4: Introduction to Linear Control Systems

EE3CL4: Introduction to Linear Control Systems 1 / 30 EE3CL4: Introduction to Linear Control Systems Section 9: of and using Techniques McMaster University Winter 2017 2 / 30 Outline 1 2 3 4 / 30 domain analysis Analyze closed loop using open loop

More information

Nyquist Criterion For Stability of Closed Loop System

Nyquist Criterion For Stability of Closed Loop System Nyquist Criterion For Stability of Closed Loop System Prof. N. Puri ECE Department, Rutgers University Nyquist Theorem Given a closed loop system: r(t) + KG(s) = K N(s) c(t) H(s) = KG(s) +KG(s) = KN(s)

More information

Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming

Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming Gorka Galdos, Alireza Karimi and Roland Longchamp Abstract A quadratic programming approach is proposed to tune fixed-order

More information

Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year

Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year Automatic Control 2 Loop shaping Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 21-211 1 / 39 Feedback

More information

Nyquist Plots / Nyquist Stability Criterion

Nyquist Plots / Nyquist Stability Criterion Nyquist Plots / Nyquist Stability Criterion Given Nyquist plot is a polar plot for vs using the Nyquist contour (K=1 is assumed) Applying the Nyquist criterion to the Nyquist plot we can determine the

More information

Analysis of SISO Control Loops

Analysis of SISO Control Loops Chapter 5 Analysis of SISO Control Loops Topics to be covered For a given controller and plant connected in feedback we ask and answer the following questions: Is the loop stable? What are the sensitivities

More information

ECE382/ME482 Spring 2005 Homework 7 Solution April 17, K(s + 0.2) s 2 (s + 2)(s + 5) G(s) =

ECE382/ME482 Spring 2005 Homework 7 Solution April 17, K(s + 0.2) s 2 (s + 2)(s + 5) G(s) = ECE382/ME482 Spring 25 Homework 7 Solution April 17, 25 1 Solution to HW7 AP9.5 We are given a system with open loop transfer function G(s) = K(s +.2) s 2 (s + 2)(s + 5) (1) and unity negative feedback.

More information

NPTEL Online Course: Control Engineering

NPTEL Online Course: Control Engineering NPTEL Online Course: Control Engineering Ramkrishna Pasumarthy Assignment-11 : s 1. Consider a system described by state space model [ ] [ 0 1 1 x + u 5 1 2] y = [ 1 2 ] x What is the transfer function

More information

1 (20 pts) Nyquist Exercise

1 (20 pts) Nyquist Exercise EE C128 / ME134 Problem Set 6 Solution Fall 2011 1 (20 pts) Nyquist Exercise Consider a close loop system with unity feedback. For each G(s), hand sketch the Nyquist diagram, determine Z = P N, algebraically

More information

Optimal control and estimation

Optimal control and estimation Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011

More information

(Refer Slide Time: 2:11)

(Refer Slide Time: 2:11) Control Engineering Prof. Madan Gopal Department of Electrical Engineering Indian institute of Technology, Delhi Lecture - 40 Feedback System Performance based on the Frequency Response (Contd.) The summary

More information

ECE317 : Feedback and Control

ECE317 : Feedback and Control ECE317 : Feedback and Control Lecture : Steady-state error Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling Analysis Design Laplace

More information

Closed Loop Identification Of A First Order Plus Dead Time Process Model Under PI Control

Closed Loop Identification Of A First Order Plus Dead Time Process Model Under PI Control Dublin Institute of Technology RROW@DIT Conference papers School of Electrical and Electronic Engineering -6- Closed Loop Identification Of First Order Plus Dead Time Process Model Under PI Control Tony

More information

Lecture 15 Nyquist Criterion and Diagram

Lecture 15 Nyquist Criterion and Diagram Lecture Notes of Control Systems I - ME 41/Analysis and Synthesis of Linear Control System - ME86 Lecture 15 Nyquist Criterion and Diagram Department of Mechanical Engineering, University Of Saskatchewan,

More information

6 OUTPUT FEEDBACK DESIGN

6 OUTPUT FEEDBACK DESIGN 6 OUTPUT FEEDBACK DESIGN When the whole sate vector is not available for feedback, i.e, we can measure only y = Cx. 6.1 Review of observer design Recall from the first class in linear systems that a simple

More information

Topic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback

Topic # 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 information

Frequency domain analysis

Frequency domain analysis Automatic Control 2 Frequency domain analysis Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011

More information

D(s) G(s) A control system design definition

D(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 information

IC6501 CONTROL SYSTEMS

IC6501 CONTROL SYSTEMS DHANALAKSHMI COLLEGE OF ENGINEERING CHENNAI DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING YEAR/SEMESTER: II/IV IC6501 CONTROL SYSTEMS UNIT I SYSTEMS AND THEIR REPRESENTATION 1. What is the mathematical

More information

Control Systems Design

Control 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 information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3. 8. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid

More information

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Leonid Freidovich Department of Mathematics Michigan State University MI 48824, USA e-mail:leonid@math.msu.edu http://www.math.msu.edu/

More information

Introduction to Nonlinear Control Lecture # 4 Passivity

Introduction 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 information

Frequency Response Techniques

Frequency Response Techniques 4th Edition T E N Frequency Response Techniques SOLUTION TO CASE STUDY CHALLENGE Antenna Control: Stability Design and Transient Performance First find the forward transfer function, G(s). Pot: K 1 = 10

More information

CDS 101/110: Lecture 3.1 Linear Systems

CDS 101/110: Lecture 3.1 Linear Systems CDS /: Lecture 3. Linear Systems Goals for Today: Describe and motivate linear system models: Summarize properties, examples, and tools Joel Burdick (substituting for Richard Murray) jwb@robotics.caltech.edu,

More information

MIMO Toolbox for Matlab

MIMO Toolbox for Matlab MIMO Toolbox for Matlab Oskar Vivero School of Electric and Electronic Engineering University of Manchester, Manchester, UK M60 1QD Email: oskar.vivero@postgrad.manchester.ac.uk Jesús Liceaga-Castro Departamento

More information