STABILITY ANALYSIS TECHNIQUES
|
|
- Polly Johnston
- 6 years ago
- Views:
Transcription
1 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 Here, we look at determining system stability using various methods DEFINITION: AsystemisBIBOstableiffaboundedinputproducesa bounded output Check by first writing system input output relationship as G(z) Y (z) = 1 + GH(z) R(z) = K m (z z i ) n (z p i ) R(z) Assume for now that all the poles {pi } are distinct and different from the poles in R(z) Then, Y (z) = k 1z + + k nz } z p 1 {{ z p n} Response to initial conditions + Y R (z) }{{} Response to R(z) If the system is stable, the response to initial conditions must decay to zero as time progresses Z 1 [ ki z z p i So, the system is stable if p i < 1 ] = k i (p i ) k 1[k]
2 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 2 {pi } are the roots of 1 + GH(z) = 0 So,therootsof1 + GH(z) = 0 must lie within the unit circle of the z-plane Same result even if poles are repeated, but harder to show If the magnitude of a pole pi =1, thenthesystemismarginally stable The unforced response does not decay to zero but also does not increase to However, itispossibletodrivethesystemwitha bounded input and have the output go to Therefore,amarginally stable system is unstable Bilinear transformation The stability criteria for a discrete-time system is that all itspoleslie within the unit circle on the z-plane Stability criteria for cts-time systems is that the poles be inthelhp Simple tool to test for continuous-time stability Routh test Can we use the Routh test to determine stability of a discrete-time system (either directly or indirectly)? To use the Routh test, we need to do a z-plane to s-plane conversion that retains stability information The s-plane version of the z-plane system does NOT need to correspond in any other way That is, The frequency responses may be different The step responses may be different z-plane w-plane
3 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 3 Since only stability properties are maintained by the transform, it is not accurate to label the destination plane the s-plane It is often called the w-plane, and the transformation between the z-plane and the w-plane is called the w-transform Atransformthatsatisfiestheserequirementsisthebilineartransform Recall: H(w) = H(z) z= 1+(T/2)w 1 (T/2)w and H(z) = H(w) w= 2 z 1 T z+1 Three things to check: 1 Unit circle in z-plane jω-axis in w-plane 2 Inside unit circle in z-plane LHP in w-plane 3 Outside unit circle in z-plane RHP in w-plane If true, 1 Take H(z) H(w) via the bilinear transform 2 Perform Routh test on H(w) CHECK: Let z = re jωt Then,zis on the unit circle if r = 1, z is inside the unit circle if r < 1 and z is outside the unit circle if r > 1 z = re jωt w = 2 z 1 T z + 1 z=re jωt = 2 T re jωt 1 re jωt + 1 Expand e jωt = cos(ωt ) + j sin(ωt ) and use the shorthand c = cos(ωt ) and s = sin(ωt ) Alsonotethats 2 + c 2 = 1 w = 2 [ ] rc + jrs 1 T rc + jrs + 1
4 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 4 ] +1 = 2 T = 2 T = 2 T [ ][ (rc 1) + jrs (rc + 1) jrs (rc + 1) + jrs (rc + 1) jrs [ (r 2 c 2 1) + j (rs)(rc + 1) j (rs)(rc 1) + r 2 s 2 ] (rc + 1) 2 + (rs) 2 [ r 2 ] 1 + j 2 [ ] 2rs r 2 + 2rc + 1 T r 2 + 2rc + 1 Notice that the real part of w is 0 when r = 1 (w is on the imaginary axis), the real part of w is negative when r < 1 (w in LHP), and that the real part of w is positive when r > 1 (w in RHP) Therefore, the bilinear transformation does exactly what we want When r = 1, w = j 2 T which will be useful to know 2sin(ωT ) 2 + 2cos(ωT ) = j 2 ( ) ωt T tan, 2 The following diagram summarizes the relationship between the s-plane, z-plane, and w-plane: ➃ ➄ s-plane z-plane w-plane j ω s 2 j ω s 4 j ω s j ω 4 s 2 ➂ ➁ ➀ ➆ ➅ ➂ ➅ ➃ ➄ ➁ ➆ R = 1 ➀ ➂ ➅ ➃ R = 2 ➄ T ➁ j 2 T ➀ ➆
5 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES : Discrete-time stability via Routh Hurwitz test Review of Routh test Let H(w) = b(w) a(w) is the characteristic polynomial a(w) a(w) = a n w n + a n 1 w n 1 + +a 1 w + a 0 Case 0: Case 1: If any of the a n are negative then the system is unstable (unless ALL are negative) Form Routh array: w n a n a n 2 a n 4 w n 1 a n 1 a n 3 a n 5 w n 2 b 1 b 2 w n 3 c 1 c 2 w 1 j 1 b 1 = 1 a n 1 a n a n 1 w 0 k 1 a n 2 a n 3 b 2 = 1 a n 1 a n a n 1 a n 4 a n 5 c 1 = 1 b 1 a n 1 a n 3 b 1 b 2 c 2 = 1 b 1 a n 1 a n 5 b 1 b 3 TEST: Number of RHP roots = number of sign changes in left column Case 2: Case 3: If one of the left column entries is zero, replace it with ϵ as ϵ 0 Suppose an entire row of the Routh array is zero, the w i 1 th row The w i th row, right above it, has coefficients α 1,α 2,
6 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 6 Then, form the auxiliary equation: α 1 w i + α 2 w i 2 + α 3 w i 4 + =0 This equation is a factor of the characteristic equation and must be tested for RHP roots (it WILL have non-lhp roots we might want to know how many are RHP) EXAMPLE: Consider: r(t) T K 1 e st s 1 s(s + 1) y(t) G(s) = ( 1 e Ts)( ) 1 s s(s + 1) From z-transform tables: ( ) [ ] z 1 1 G(z) = Z z s 2 (s + 1) ( )( z 1 (e T + T 1)z 2 + (1 e T Te T ) )z = z (z 1) 2 (z e T ) Let T = 01s = z (z 1)(z 0905) Perform the bilinear transform G(w) = G(z) z= 1+(T/2)w 1 (T/2)w = G(z) z= 1+005w 1 005w The characteristic equation is: = w w w w
7 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 7 (system 0 = 1 + KG(w) = ( K )w 2 + ( K )w + 381K w 2 ( K ) 381K K < 23, 813 w 1 ( K ) K < 203 w 0 381K K > 0 So, for stability, 0 < K < 203 NOTE: The equivalent continuous-timesystemis: r(t) K 1 s(s + 1) y(t) T (s) = KG(s) 1 + KG(s) Characteristic equation: s(s + 1) + K = 0 s 2 1 K s 1 1 s 0 K Stable for all K > 0 sample and hold destabilizes the system EXAMPLE: Let sdothesameexample,butwitht = 1s(not01s) (math happens) 0 = 1 + KG(w) = ( K )w 2 + ( K )w K w 2 ( K ) 0924K K < 262 w 1 ( K ) K < 239 w K K > 0
8 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 8 So, for stability, 0 < K < 239 This is a much more restrictive range than when T = 01s slow sampling really destabilizes a system
9 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES : Jury s stability test H(z) H(w) Routh is complicated and error-prone Jury made a direct test on H(z) for stability Disadvantage (?) another test to learn Let T (z) = b(z) a(z), a(z) = characteristic polynomial a(z) = an z n + a n 1 z n 1 + +a 1 z + a 0 = 0, a n > 0 Form Jury array: z 0 z 1 z 2 z n k z n 1 z n a 0 a 1 a 2 a n k a n 1 a n a n a n 1 a n 2 a k a 1 a 0 b 0 b 1 b 2 b n k b n 1 b n 1 b n 2 b n 3 b k 1 b 0 c 0 c 1 c 2 c n k c n 2 c n 3 c n 4 c k 2 l 0 l 1 l 2 l 3 l 3 l 2 l 1 l 0 m 0 m 1 m 2 Quite different from Routh array Every row is duplicated in reverse order Final row in table has three entries (always) Elements are calculated differently a 0 a n k b 0 b k = c k = a n a k b n 1 b n 1 k b k d k = c 0 c n 2 c n 2 k c k
10 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 10 Stability criteria is different a(z) z=1 > 0 ( 1) n a(z) z= 1 > 0 n = order of a(z) a 0 < a n b 0 > b n 1 c 0 > c n 2 d 0 > d n 3 m 0 > m 2 First, check that a(1) >0, ( 1) n a( 1) >0 and a 0 < a n (relatively few calculations) If not satisfied, stop Next, construct array Stop if any condition not satisfied EXAMPLE: r(t) T = 1s 1 e st s K s(s + 1) y(t) r[k] K 0368z z z y[k] Characteristic equation: 0 = 1 + KG(z) = 1 + K (0368z ) z z
11 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 11 The Jury array is: = z 2 + (0368K 1368)z + ( K ) z 0 z 1 z K 0368K The constraint a(1) >0 yields K K = 0632K > 0 K > 0 The constraint ( 1) 2 a( 1) >0 yields K K = 0104K > 0 K < 263 The constraint a0 < a 2 yields K < 1 K < = 239 So, 0 < K < 239 (Same result as on pg 4 8 using bilinear rule) EXAMPLE: Supposethatthecharacteristicequationforaclosed-loop discrete-time system is given by the expression: a(z) = z 3 18z z 020 = 0 a(1) = = 005 > 0 ( 1) 3 a( 1) = [ ] > 0 a0 =02 < a 3 = 1 Jury array: z 0 z 1 z 2 z
12 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 12 b 0 = = 096 b 1 = = 159 b 2 = = 069 b0 =096 > b 2 =069 The system is stable
13 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 13 ŷ 44: Root-locus and Nyquist tests For cts-time control, we examined the locations of the roots ofthe closed-loop system as a function of the loop gain K Root locus r(t) K D(s) G(s) y(t) H(s) T (s) = KD(s)G(s) 1 + KD(s)G(s)H(s) Let L(s) = D(s)G(s)H(s) (The looptransferfunction ) Developed rules for plotting the roots of the equation Root Locus Drawing Rules Applied them to plotting roots of 1 + K b(s) a(s) = KL(s) = 0 Now, we have the digital system: r(t) K D(z) T G(s) {}}{ 1 e st G p (s) s H(s) y(t) T (z) = So, we let L(z) = D(z)GH(z) Poles are roots of 1 + KL(z) = 0 KD(z)G(z) 1 + KD(z)GH(z),
14 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 14 This is exactly the same form as the Laplace-transform root locus Plot roots in exactly the same way EXAMPLE: GH(z) = 0368(z ) (z 1)(z 0368) D(z) = 1 K = 239 numd=0368*[1 0717]; dend=conv([1-1],[1-0368]); d=tf(numd,dend,-1); rlocus(d); The Nyquist test In continuous-time control we also used the Nyquist test to assess stability I(s) I(s) Zoom II III ρ 0 I R(s) θ R(s) IV The Nyquist D path encircles the entire (unstable) RHP The Nyquist plot is a polar plot of L(s) evaluated on the D path Adjustments to D shape are made if pole on the jω-axis
15 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 15 The Nyquist test evaluated stability by looking at the Nyquist plot N=No of CW encirclements of 1 in Nyquist plot P=No of open-loop unstable poles (poles inside D shape) Z=No of closed-loop unstable poles Z = N + P, EXAMPLE: r(t) K s(s + 1) Z = 0 for stable closed-loop system y(t) This gives: L(s) = 1 s(s + 1) Pole at origin: Need detour s = ρe jθ, ρ 1 Resulting Nyquist map has infinite radius Cannot draw to scale No poles inside modified- D curve: P = 0 Z = N + P = 0 Stable system Note that increasing the gain K onlymagnifiestheentireplotthe 1 point is not encircled for K > 0 (infinite gain margin) Nyquist test for discrete systems Three different ways to do the Nyquist test for discrete systems Based on three different representations of the characteristic eqn L (s) = 0 L = DGH L(z) = L(w) = 0
16 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES L (s) = 0 We know that L (s) is periodic in jω s Therefore, the D curve does not need to encircle the entire RHP to encircle all unstable poles [If there were any, there would be an infinite number] Modify D curve to be: j ω s 2 j ω s 2 I(s) R(s) Evaluate L (s) on new contour and plot polar plot Same Nyquist test as before L(z) = 0 We can do the Nyquist test directly using z-transforms The stable region is the unit circle The z-domain Nyquist plot is done using a Nyquist curve which is the unit circle Nyquist test changes because we are now encircling the STABLE region (albeit CCW) Z = # closed-loop unstable poles I(z) P = # open-loop unstable poles N = # CCW encirclements of 1 in Nyquist plot R(z) Z = P N Probably difficult to evaluate L(z) z=e jθ for π θ π unless using adigitalcomputer
17 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES L(w) = 0 Now, we convert L(z) L(w) I(w) L(w) = L(z) z= 1+(T/2)w 1 (T/2)w Bilinear transform maps unit circle to jω-axis in w-plane R(w) Use standard continuous-time test in w-plane Summary: Open-loop fn Range of variable Rule GH (s) s = jω, ω s /2 ω ω s /2 Z = P + N cw GH(z) z = e jωt, π ωt π Z = P N ccw = P + N cw GH(w) w = jω w, ω w Z = P + N cw All three methods produce identical Nyquist plots Nyquist Diagrams Note that the sampled system does not have gain margin (a = 0418, GM = 239) andhassmaller PM than cts-time system Imaginary Axis PM a Real Axis
18 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES : Bode methods Bode plots are an extremely important tool for analyzing and designing control systems They provide a critical link between continuous-time and discrete-time control design methods Recall: Bode plots are plots of frequency response of a system: Magnitude and Phase In s-plane, H(s) s= jω is frequency response for 0 ω< In z-plane, H(z) z=e jωt is frequency response for 0 ω ω s /2 Straight-line tools of s-plane analysis DON T WORK! They are based on geometry and geometry has changed jω-axis to z-unit circle BUT in w-plane, H(w) w= jωw is the frequency response for 0 ω w < Straight-linetoolswork,butfrequencyaxisiswarped PROCEDURE: 1 Convert H(z) to H(w) by H(w) = H(z) z= 1+(T/2)w 1 (T/2)w 2 Simplify expression to rational-polynomial in w 3 Factor into zeros and poles in standard Bode Form (Refer to review notes) 4 Plot the response exactly the same way as an s-plane Bode plot Note: Plots are versus log 10 ω w ω w = 2 ( ) ωt T tan Can 2 re-scale axis in terms if ω if we want EXAMPLE: Example seen before with T = 1 second
19 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 19 (1,2) (3) (4) G(w) = = = [ 1+05w 1 05w Let G(z) = 0368z z z [ 1+05w 1 05w] ] 2 [ w 1 05w] (1 + 05w)(1 05w) (1 05w) 2 (1 + 05w) (1 + 05w)(1 05w) (1 05w) (w 2)(w ) w(w ) Magnitude (db) G( jω w ) = ( j ω w 2 1 )( j ω w 1214 ( + 1) jω w j ω w ) Bode Plots Phase (deg) Frequency (warped rads/sec) Gain margin and phase margin work the SAME way we expect
20 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 20 WAIT! We have discussed frequency-response methods without verifying that discrete-time frequency response means the same thing as continuous-time frequency response Verify X (z) G(z) Y (z) Let x[k] =sin(ωkt) X (z) = z sin ωt (z e jωt )(z e jωt ) Y (z) = G(z)X (z) = G(z)z sin ωt (z e jωt )(z e jωt ) Do partial-fraction expansion Y (z) k 1 = z z e + k 2 jωt z e + Y g(z) jωt Yg (z) is the response due to the poles of G(z) IF the system is stable, the response due to Y g (z) 0 as t So, as t we say Y ss(z) z = k 1 = k 1 z e + k 2 jωt z e jωt G(z) sin ωt z e jωt z=e jωt = G(e jωt ) sin ωt e jωt e jωt = G(e jωt ) 2 j = G(e jωt ) e j G(e jωt ) 2 j
21 ECE4540/5540, STABILITY ANALYSIS TECHNIQUES 4 21 Similarly, k 2 = G(e jωt ) e j G(e jωt ) 2( j) = G(e jωt ) e j G(e jωt ) 2 j Combining and solving for yss [k] y ss [k] =k 1 (e jωt ) k + k 2 (e jωt ) k = G(e jωt ) e jωkt+ j G(e jωt ) e jωkt j G(e jωt ) 2 j = G(e jωt ) sin(ωkt + G(e jωt )) Sure enough, G(e jωt ) is magnitude response to sinusoid, and G(e jωt ) is phase response to sinusoid Closed-loop frequency response We have looked at open-loop concepts and how they apply to closed loop systems our end product Closed-loop frequency response usually calculated by computer: G(z) 1 + G(z),forexample In general, if G(e jωt ) large, T (e jωt ) 1 If G(e jωt ) small, T (e jωt ) G(e jωt ) Closed-loop bandwidth similar to open-loop bandwidth If PM = 90,thenCLBW= OL BW If PM = 45,thenCLBW= 2 OL BW
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 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 informationDESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD
206 Spring Semester ELEC733 Digital Control System LECTURE 7: DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD For a unit ramp input Tz Ez ( ) 2 ( z ) D( z) G( z) Tz e( ) lim( z) z 2 ( z ) D( z)
More informationFREQUENCY-RESPONSE ANALYSIS
ECE450/550: Feedback Control Systems. 8 FREQUENCY-RESPONSE ANALYSIS 8.: Motivation to study frequency-response methods Advantages and disadvantages to root-locus design approach: ADVANTAGES: Good indicator
More informationControl 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 information1 (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 informationSTABILITY ANALYSIS. Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated using cones: Stable Neutral Unstable
ECE4510/5510: Feedback Control Systems. 5 1 STABILITY ANALYSIS 5.1: Bounded-input bounded-output (BIBO) stability Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated
More informationSystems 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 informationControl 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 informationDIGITAL CONTROLLER DESIGN
ECE4540/5540: Digital Control Systems 5 DIGITAL CONTROLLER DESIGN 5.: Direct digital design: Steady-state accuracy We have spent quite a bit of time discussing digital hybrid system analysis, and some
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 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 informationTopic # 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 informationROOT 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 informationControl 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 informationSolutions to Skill-Assessment Exercises
Solutions to Skill-Assessment Exercises To Accompany Control Systems Engineering 4 th Edition By Norman S. Nise John Wiley & Sons Copyright 2004 by John Wiley & Sons, Inc. All rights reserved. No part
More informationECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27
1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING
KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK SUB.NAME : CONTROL SYSTEMS BRANCH : ECE YEAR : II SEMESTER: IV 1. What is control system? 2. Define open
More informationFREQUENCY-RESPONSE DESIGN
ECE45/55: Feedback Control Systems. 9 FREQUENCY-RESPONSE DESIGN 9.: PD and lead compensation networks The frequency-response methods we have seen so far largely tell us about stability and stability margins
More informationr + - 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 informationRobust 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 informationLecture 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 informationPlan of the Lecture. Goal: wrap up lead and lag control; start looking at frequency response as an alternative methodology for control systems design.
Plan of the Lecture Review: design using Root Locus; dynamic compensation; PD and lead control Today s topic: PI and lag control; introduction to frequency-response design method Goal: wrap up lead and
More informationFrequency 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 informationIntro 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 informationECE317 : Feedback and Control
ECE317 : Feedback and Control Lecture : Stability Routh-Hurwitz stability criterion Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling
More informationLecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types
Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types Venkata Sonti Department of Mechanical Engineering Indian Institute of Science Bangalore, India, 562 This
More informationFrequency 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 informationStability Analysis Techniques
Stability Analysis Techniques In this section the stability analysis techniques for the Linear Time-Invarient (LTI) discrete system are emphasized. In general the stability techniques applicable to LTI
More informationFrequency 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 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 informationData Based Design of 3Term Controllers. Data Based Design of 3 Term Controllers p. 1/10
Data Based Design of 3Term Controllers Data Based Design of 3 Term Controllers p. 1/10 Data Based Design of 3 Term Controllers p. 2/10 History Classical Control - single controller (PID, lead/lag) is designed
More informationModule 3F2: Systems and Control EXAMPLES PAPER 2 ROOT-LOCUS. Solutions
Cambridge University Engineering Dept. Third Year Module 3F: Systems and Control EXAMPLES PAPER ROOT-LOCUS Solutions. (a) For the system L(s) = (s + a)(s + b) (a, b both real) show that the root-locus
More information1 (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 informationChapter 7. Digital Control Systems
Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,
More informationMEM 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 informationControl Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:
More informationSoftware Engineering/Mechatronics 3DX4. Slides 6: Stability
Software Engineering/Mechatronics 3DX4 Slides 6: Stability Dr. Ryan Leduc Department of Computing and Software McMaster University Material based on lecture notes by P. Taylor and M. Lawford, and Control
More information6.302 Feedback Systems Recitation 27: Final Recitation and Review Prof. Joel L. Dawson
6.302 Feedbac Systems We re going to spend some time in this recitation revisiting some of the very early examples of recitation #1. To aid us in this review, the final class exercise involves exploring
More informationAnalysis 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 informationDelhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:
Serial : 0. LS_D_ECIN_Control Systems_30078 Delhi Noida Bhopal Hyderabad Jaipur Lucnow Indore Pune Bhubaneswar Kolata Patna Web: E-mail: info@madeeasy.in Ph: 0-4546 CLASS TEST 08-9 ELECTRONICS ENGINEERING
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 informationClass 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 informationDigital 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 informationINTRODUCTION TO DIGITAL CONTROL
ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant
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 informationRoot Locus Methods. The root locus procedure
Root Locus Methods Design of a position control system using the root locus method Design of a phase lag compensator using the root locus method The root locus procedure To determine the value of the gain
More informationECE 486 Control Systems
ECE 486 Control Systems Spring 208 Midterm #2 Information Issued: April 5, 208 Updated: April 8, 208 ˆ This document is an info sheet about the second exam of ECE 486, Spring 208. ˆ Please read the following
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 24: Compensation in the Frequency Domain Overview In this Lecture, you will learn: Lead Compensators Performance Specs Altering
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 informationThe Frequency-response Design Method
Chapter 6 The Frequency-response Design Method Problems and Solutions for Section 6.. (a) Show that α 0 in Eq. (6.2) is given by α 0 = G(s) U 0ω = U 0 G( jω) s jω s= jω 2j and α 0 = G(s) U 0ω = U 0 G(jω)
More information(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications:
1. (a) The open loop transfer function of a unity feedback control system is given by G(S) = K/S(1+0.1S)(1+S) (i) Determine the value of K so that the resonance peak M r of the system is equal to 1.4.
More informationMAE 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 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 informationEE402 - Discrete Time Systems Spring Lecture 10
EE402 - Discrete Time Systems Spring 208 Lecturer: Asst. Prof. M. Mert Ankarali Lecture 0.. Root Locus For continuous time systems the root locus diagram illustrates the location of roots/poles of a closed
More informationVALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur 603 203. DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING SUBJECT QUESTION BANK : EC6405 CONTROL SYSTEM ENGINEERING SEM / YEAR: IV / II year
More informationExplanations and Discussion of Some Laplace Methods: Transfer Functions and Frequency Response. Y(s) = b 1
Engs 22 p. 1 Explanations Discussion of Some Laplace Methods: Transfer Functions Frequency Response I. Anatomy of Differential Equations in the S-Domain Parts of the s-domain solution. We will consider
More informationControl Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho Tel: Fax:
Control Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 02-760-4253 Fax:02-760-4435 Introduction In this lesson, you will learn the following : How to determine
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 informationNyquist 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 informationECE317 : 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 informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (HS 25) Solution Exercise Set Loop Shaping Noele Norris, 9th December 26 Exercise (A Non-minimum Phase System) To increase the rise time of the system, we
More informationSoftware Engineering 3DX3. Slides 8: Root Locus Techniques
Software Engineering 3DX3 Slides 8: Root Locus Techniques Dr. Ryan Leduc Department of Computing and Software McMaster University Material based on Control Systems Engineering by N. Nise. c 2006, 2007
More informationControl 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 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 informationECE 345 / ME 380 Introduction to Control Systems Lecture Notes 8
Learning Objectives ECE 345 / ME 380 Introduction to Control Systems Lecture Notes 8 Dr. Oishi oishi@unm.edu November 2, 203 State the phase and gain properties of a root locus Sketch a root locus, by
More informationLecture 8 - IIR Filters (II)
Lecture 8 - IIR Filters (II) James Barnes (James.Barnes@colostate.edu) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423 1 / 27 Lecture 8 Outline Introduction Digital
More informationIf you need more room, use the backs of the pages and indicate that you have done so.
EE 343 Exam II Ahmad F. Taha Spring 206 Your Name: Your Signature: Exam duration: hour and 30 minutes. This exam is closed book, closed notes, closed laptops, closed phones, closed tablets, closed pretty
More information7.4 STEP BY STEP PROCEDURE TO DRAW THE ROOT LOCUS DIAGRAM
ROOT LOCUS TECHNIQUE. Values of on the root loci The value of at any point s on the root loci is determined from the following equation G( s) H( s) Product of lengths of vectors from poles of G( s)h( s)
More informationELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 2010/2011 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques
CAIRO UNIVERSITY FACULTY OF ENGINEERING ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 00/0 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques [] For the following system, Design a compensator such
More informationIntroduction. 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 informationLaplace Transform Analysis of Signals and Systems
Laplace Transform Analysis of Signals and Systems Transfer Functions Transfer functions of CT systems can be found from analysis of Differential Equations Block Diagrams Circuit Diagrams 5/10/04 M. J.
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 informationH(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 informationCONTROL 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 informationTransient Response of a Second-Order System
Transient Response of a Second-Order System ECEN 830 Spring 01 1. Introduction In connection with this experiment, you are selecting the gains in your feedback loop to obtain a well-behaved closed-loop
More informationECE382/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 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 informationBoise State University Department of Electrical Engineering ECE461 Control Systems. Control System Design in the Frequency Domain
Boise State University Department of Electrical Engineering ECE6 Control Systems Control System Design in the Frequency Domain Situation: Consider the following block diagram of a type- servomechanism:
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 informationCHAPTER # 9 ROOT LOCUS ANALYSES
F K א CHAPTER # 9 ROOT LOCUS ANALYSES 1. Introduction The basic characteristic of the transient response of a closed-loop system is closely related to the location of the closed-loop poles. If the system
More informationCourse 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 informationCourse Summary. The course cannot be summarized in one lecture.
Course Summary Unit 1: Introduction Unit 2: Modeling in the Frequency Domain Unit 3: Time Response Unit 4: Block Diagram Reduction Unit 5: Stability Unit 6: Steady-State Error Unit 7: Root Locus Techniques
More informationMEM 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 informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 05 IIR Design 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationDamped Oscillators (revisited)
Damped Oscillators (revisited) We saw that damped oscillators can be modeled using a recursive filter with two coefficients and no feedforward components: Y(k) = - a(1)*y(k-1) - a(2)*y(k-2) We derived
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 12: Overview In this Lecture, you will learn: Review of Feedback Closing the Loop Pole Locations Changing the Gain
More informationControl Systems. EC / EE / IN. For
Control Systems For EC / EE / IN By www.thegateacademy.com Syllabus Syllabus for Control Systems Basic Control System Components; Block Diagrammatic Description, Reduction of Block Diagrams. Open Loop
More informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (Autumn 27) Solution Exercise Set 2 Loop Shaping clruch@ethz.ch, 8th December 27 Exercise (A Non-minimum Phase System) To decrease the rise time of the system,
More informationReturn 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 informationThe stability of linear time-invariant feedback systems
The stability of linear time-invariant feedbac systems A. Theory The system is atrictly stable if The degree of the numerator of H(s) (H(z)) the degree of the denominator of H(s) (H(z)) and/or The poles
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 informationControl Systems. University Questions
University Questions UNIT-1 1. Distinguish between open loop and closed loop control system. Describe two examples for each. (10 Marks), Jan 2009, June 12, Dec 11,July 08, July 2009, Dec 2010 2. Write
More informationLecture 3 - Design of Digital Filters
Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can
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 informationUse: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)}
1 VI. Z Transform Ch 24 Use: Analysis of systems, simple convolution, shorthand for e jw, stability. A. Definition: X(z) = x(n) z z - transforms Motivation easier to write Or Note if X(z) = Z {x(n)} z
More informationStep Response Analysis. Frequency Response, Relation Between Model Descriptions
Step Response Analysis. Frequency Response, Relation Between Model Descriptions Automatic Control, Basic Course, Lecture 3 November 9, 27 Lund University, Department of Automatic Control Content. Step
More informationEC6405 - CONTROL SYSTEM ENGINEERING Questions and Answers Unit - I Control System Modeling Two marks 1. What is control system? A system consists of a number of components connected together to perform
More informationToday. ESE 531: Digital Signal Processing. IIR Filter Design. Impulse Invariance. Impulse Invariance. Impulse Invariance. ω < π.
Today ESE 53: Digital Signal Processing! IIR Filter Design " Lec 8: March 30, 207 IIR Filters and Adaptive Filters " Bilinear Transformation! Transformation of DT Filters! Adaptive Filters! LMS Algorithm
More information