EE Control Systems

Similar documents
Bode Diagrams School of Mechanical Engineering ME375 Frequency Response - 29 Purdue University Example Ex:

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

CDS 101: Lecture 8.2 Tools for PID & Loop Shaping

732 Appendix E: Previous EEE480 Exams. Rules: One sheet permitted, calculators permitted. GWC 352,

1the 1it is said to be overdamped. When 1, the roots of

Lab(8) controller design using root locus

School of Mechanical Engineering Purdue University. ME375 Frequency Response - 1

Frequency Response Methods

CHAPTER NINE. Frequency Response Methods

MEM 255 Introduction to Control Systems: Analyzing Dynamic Response

CONTROL SYSTEMS. Chapter 7 : Bode Plot. 40dB/dec 1.0. db/dec so resultant slope will be 20 db/dec and this is due to the factor s

Chapter 2 Feedback Control Theory Continued

2.004 Dynamics and Control II Spring 2008

Answer: 1(A); 2(C); 3(A); 4(D); 5(B); 6(A); 7(C); 8(C); 9(A); 10(A); 11(A); 12(C); 13(C)

Dynamic Response of Linear Systems

Chapter 7: The z-transform. Chih-Wei Liu

Exponential Moving Average Pieter P

EE 205 Dr. A. Zidouri. Electric Circuits II. Frequency Selective Circuits (Filters) Bode Plots: Complex Poles and Zeros.

Class 07 Time domain analysis Part II 2 nd order systems

Lecture 13. Graphical representation of the frequency response. Luca Ferrarini - Basic Automatic Control 1

COMM 602: Digital Signal Processing

University of California at Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences

Course Outline. Problem Identification. Engineering as Design. Amme 3500 : System Dynamics and Control. System Response. Dr. Stefan B.

ME 375 FINAL EXAM Friday, May 6, 2005

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Dr. Seeler Department of Mechanical Engineering Fall 2009 Lafayette College ME 479: Control Systems and Mechatronics Design and Analysis

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Chapter 7 z-transform

EE 4343 Lab#4 PID Control Design of Rigid Bodies

Chapter 4 : Laplace Transform

Crash course part 2. Frequency compensation

Dynamic System Response

The z-transform can be used to obtain compact transform-domain representations of signals and systems. It

Ch3 Discrete Time Fourier Transform

Sinusoidal Steady-state Analysis

Jitter Transfer Functions For The Reference Clock Jitter In A Serial Link: Theory And Applications

Voltage controlled oscillator (VCO)

Contents Kreatryx. All Rights Reserved.

MM4 System s Poles and Feedback Characteristics

Introduction to Control Systems

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 6 - Small Signal Stability

Introduction to Signals and Systems, Part V: Lecture Summary

EE Midterm Test 1 - Solutions

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

6.003: Signals and Systems. Feedback, Poles, and Fundamental Modes

UNIVERSITY OF BOLTON SCHOOL OF ENGINEERING BENG (HONS) IN MECHANICAL ENGINEERING SEMESTER 1 EXAMINATION 2016/2017

Automatic Control Systems

6.003 Homework #12 Solutions

Appendix: The Laplace Transform

Part I: Covers Sequence through Series Comparison Tests

A New Accurate Analytical Expression for Rise Time Intended for Mechatronics Systems Performance Evaluation and Validation

6.003 Homework #12 Solutions

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

6.003 Homework #3 Solutions

2C09 Design for seismic and climate changes

EE422G Homework #13 (12 points)

Run-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE

Time-Domain Representations of LTI Systems

Curve Sketching Handout #5 Topic Interpretation Rational Functions

Lecture 3. Digital Signal Processing. Chapter 3. z-transforms. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev. 2016

Solutions of Chapter 5 Part 1/2

Exam. Notes: A single A4 sheet of paper (double sided; hand-written or computer typed)

Analysis of the No-Load Characteristic of the Moving Coil Linear Compressor

Definition of z-transform.

Systems Analysis and Control

P.3 Polynomials and Special products

(Figure 2.9), we observe x. and we write. (b) as x, x 1. and we write. We say that the line y 0 is a horizontal asymptote of the graph of f.

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t,

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Series: Infinite Sums

MATH 31B: MIDTERM 2 REVIEW

CDS 101: Lecture 5.1 Controllability and State Space Feedback

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.

State space systems analysis

FIR Filter Design: Part II

Student s Printed Name:

The z transform is the discrete-time counterpart of the Laplace transform. Other description: see page 553, textbook.

Lesson 10: Limits and Continuity

Linear time invariant systems

Active Damping Control using Optimal Integral Force Feedback

Chapter 7: Numerical Series

Chapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:

Lecture 30: Frequency Response of Second-Order Systems

Butterworth LC Filter Designer

Chapter 3. z-transform

2D DSP Basics: 2D Systems

Fig. 2. Block Diagram of a DCS

Signals & Systems Chapter3

Orthogonal Gaussian Filters for Signal Processing

lecture 3: Interpolation Error Bounds

Numerical Methods in Fourier Series Applications

ECEN620: Network Theory Broadband Circuit Design Fall 2014

Chapter 6 Principles of Data Reduction

PID SELF-TUNING CONTROLLERS USED IN DIFFERENT STRUCTURES OF THE CONTROL LOOPS

Transcription:

Copyright FL Lewis 7 All rights reserved Updated: Moday, November 1, 7 EE 4314 - Cotrol Systems Bode Plot Performace Specificatios The Bode Plot was developed by Hedrik Wade Bode i 1938 while he worked at Bell Labs Here we shall show how performace specificatios i terms of Bode plots i the frequecy domai are related to time domai performace Badwidth ad Rise Time The Bode plot of the trasfer fuctio α 1 H() s = = s+ α s+ 1 is show The break frequecy occurs at 1 rad/sec, the magitude of the pole -3db -5-1 System: sys Frequecy (rad/sec): 998 Magitude (db): -3-15 - -5-3 -35-4 -3 Magitude (db) B -6-9 1-1 1 1 1 1 1 3 Frequecy (rad/sec) The 3dB cutoff frequecy, or badwidth, B is the frequecy at which the frequecy magitude respose has decreased by 3dB from its low frequecy value I this example B = α = 1 rad / s αt t/ τ The impulse respose of this system is ht () = e = e, where the time costat is τ = 1/ α = 1/ B 1

The step respose rise time is give by tr = τ The settlig time is ts = 5τ The time costat is iversely related to the badwidth Therefore, as badwidth icreases, the system respose becomes faster COMPLEX POLE PAIR A trasfer fuctio with a complex pair of poles ad o fiite zeros ca be writte as H ( = = s + αs + s s Δ( + ζ + The umerator is chose to scale the trasfer fuctio so that the DC gai (eg set s=) is equal to oe The deomiator is the Characteristic polyomial which ca be writte i several atural or caoical forms, icludig Δ ( = s + α s + = s + ζ s + Oe may also write Δ( = s + α s + = ( s + α) + β where β + α = These variables mea somethig i terms of time domai performace as we have see They also mea somethig i the frequecy domai, particularly the dampig ratio ζ ad the atural frequecy 5 The Bode plot for H() s = is show s + s+ 5 4 M P 3-3db Magitude (db) 1-1 - System: sys Frequecy (rad/sec): 778 Magitude (db): -3-3 -4-5 -6 r B -45-9 -135-18 1 1 1 1 Frequecy (rad/sec)

Recall that for complex poles, the step respose is faster tha τ due to the oscillatory compoets However, the settlig time is ts = 5τ ad is closely related to the badwidth; it decreases as badwidth icreases The resoat frequecy is give for ζ 77 by The maximum value of the Bode plot at resoace is give by 1 M p = ζ 1 ζ r = 1 ζ These fuctios are show i the figure From either of these, oe may compute the dampig ratio ad hece the percet overshoot i the time domai Figure 811 from Dorf ad Bishop editio 1 The quality factor 1 Q = = ζ α measures the sharpess of the resoat peak i the Bode plot Note that this is effectively determied solely by the dampig ratio The poles are complex if Q> 1/ I terms of the quality factor oe may write the characteristic polyomial i the odimesioal form 3

s Δ( = 1 s + + 1 Q Figure from Dorf ad Bishop, Moder Cotrol Systems, editio 8 Bode Desig i terms of the Ope-Loop Gai Cosider the trackig cotroller give i the figure The plat is H( ad the compesator K(; the feedback gai is k The fuctio of the tracker is to make the output y(t) follow the commad or referece iput r(t) by makig the trackig error e(t)=r(t)-y(t) small The disturbace is d(t) r(t) e(t) d(t) H( y(t) The closed-loop trasfer fuctio is 4

Y ( H ( T ( = = R( 1+ H ( The deomiator is Δ ( s ) = 1+ H ( 1+ kg( where G(=K(H( is the ope-loop gai Note that we use the same symbol for the deomiator of T( as for the state-variable characteristic polyomial Δ ( = si A However, 1+kG( is actually a polyomial fractio, whose umerator is the system characteristic polyomial May desig techiques rely o tryig to determie closed-loop properties from opeloop properties I root locus desig, oe uses the ope-loop gai G( to estimate the locatios of the closed-loop poles, which are the roots of the umerator of Δ ( s ) = 1+ kg( The key poit of RL desig is that it is easier to plot the locatios of the closed-loop poles versus the feedback gai parameter k tha it is to fid the actual closed-loop poles themselves This was extremely importat i days before digital computers whe fidig roots of high-order polyomials was difficult, ad it also gives great isight ito the properties of the closed-loop system Similarly, Bode desig uses the Bode plots of the ope-loop trasfer fuctios H( ad K(H( to select the compesator K( to give desirable closed-loop properties icludig stability, good POV, ad fast trasiet respose Steady-State Error Recall that the system is type N if there are N itegrators (ie N poles at s=) i the feedforward path K(H(=G( For zero steady-state error i respose to a uit step commad, or a uit step disturbace, oe requires the system to be of type 1 The Bode plot of the itegrator compesator K(= 1/s is give i the figure It has a costat slope of =1, or - db/decade, ad a agle of -9 o Therefore, a system of type oe has a slope of =-1 at low frequecies To get zero steady-state error i respose to a uit step, oe must add a itegrator to obtai such a slope, uless the Bode plot of H( already has this slope at low frequecies 5 Magitude (db) -5-1 -15 - -89-89 -894-896 -898-9 -9-94 -96-98 -91 1 1 1 Frequecy (rad/sec) 5

Crossover Frequecy The crossover frequecy c is where the loop gai G(= K(H( has a gai of uity, ie kg( j c ) = 1 The closed-loop trasfer fuctio is Y ( H ( T ( = = R( 1+ H ( Whe << c oe has kg( j ) >> 1 so that T( j) 1 ad the closed-loop gai is uity Whe >> c oe has kg( j ) << 1 so that T( j) kg( j) ad T( j) falls off like kg( j ) Therefore, the (closed-loop) badwidth B is about equal to the (ope loop gai) crossover frequecy To get faster resposes, oe eeds to icrease the crossover frequecy c Gai Margi, Phase Margi The gai margi ad phase margi deped o both of the ope-loop gai Bode plots, magitude ad phase By extractig iformatio from both plots, GM ad PM provide closed-loop stability ad performace criteria The closed-loop deomiator is Δ ( s ) = 1+ H ( 1+ kg( The closed-loop poles are give by Δ ( = 1 + H( 1 + kg( = or kg( s ) = 1 Therefore, stability may be studied i terms of whe kg( j ) has a magitude of oe ad a phase of 18 o The gai margi is the gai icrease required to make kg( j ) = 1 whe its phase is -18 o The phase margi is the phase shift required to make phase( kg( j )) = 18 o whe kg( j ) = 1 If the magitude of kg( j ) is α whe its phase is -18 deg, the the gai margi is 1/α The logarithm of 1/α is egative the logarithm of α Therefore, i terms of db, if α is d db, the 1/α is simply d db Note that kg( j ) = 1 meas kg( j ) = db 6

The Bode plot of kg() s = ss ( + )( s+ 3) is show Lookig at the crossover frequecy c, where the magitude is equal to oe (ie db), the phase margi = -155 (-18)= 745 deg Lookig at the frequecy where the phase is -18 deg, the gai margi is -35 db 5 Gm = 35 db (at 45 rad/sec), Pm = 745 deg (at 37 rad/sec) 35 db Magitude (db) -5-1 c 1/GM -155 o -15-9 -135-18 -5 PM -7 1-1 -1 1 1 1 1 The dampig ratio icreases with phase margi accordig to the figure Therefore, to icrease dampig ratio, we eed to icrease phase margi Relatio betwee PM ad dampig ratio for a secod order system Figure 91 from Dorf ad Bishop ed 1 7