Lead/Lag Compensator Frequency Domain Properties and Design Methods

Similar documents
Dead-beat controller design

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

, which yields. where z1. and z2

Differentiation Applications 1: Related Rates

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >

NUMBERS, MATHEMATICS AND EQUATIONS

Lecture 7: Damped and Driven Oscillations

We can see from the graph above that the intersection is, i.e., [ ).

Section I5: Feedback in Operational Amplifiers

Thermodynamics Partial Outline of Topics

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b

Preparation work for A2 Mathematics [2017]

Determining the Accuracy of Modal Parameter Estimation Methods

LHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers

Review Problems 3. Four FIR Filter Types

Cambridge Assessment International Education Cambridge Ordinary Level. Published

COMP 551 Applied Machine Learning Lecture 11: Support Vector Machines

Lecture 6: Phase Space and Damped Oscillations

Pattern Recognition 2014 Support Vector Machines

ENSC Discrete Time Systems. Project Outline. Semester

Math 105: Review for Exam I - Solutions

B. Definition of an exponential

This section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving.

COMP 551 Applied Machine Learning Lecture 5: Generative models for linear classification

AP Statistics Notes Unit Two: The Normal Distributions

Medium Scale Integrated (MSI) devices [Sections 2.9 and 2.10]

CHM112 Lab Graphing with Excel Grading Rubric

Experiment #3. Graphing with Excel

Trigonometric Ratios Unit 5 Tentative TEST date

Chapter 3 Kinematics in Two Dimensions; Vectors

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture

Root locus ( )( ) The given TFs are: 1. Using Matlab: >> rlocus(g) >> Gp1=tf(1,poly([0-1 -2])) Transfer function: s^3 + 3 s^2 + 2 s

Corrections for the textbook answers: Sec 6.1 #8h)covert angle to a positive by adding period #9b) # rad/sec

Hypothesis Tests for One Population Mean

Kinetic Model Completeness

2004 AP CHEMISTRY FREE-RESPONSE QUESTIONS

PHYS 314 HOMEWORK #3

5 th grade Common Core Standards

ECE 2100 Circuit Analysis

Computational modeling techniques

Preparation work for A2 Mathematics [2018]

Part a: Writing the nodal equations and solving for v o gives the magnitude and phase response: tan ( 0.25 )

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018

MATHEMATICS SYLLABUS SECONDARY 5th YEAR

Lab 11 LRC Circuits, Damped Forced Harmonic Motion

7 TH GRADE MATH STANDARDS

1 The limitations of Hartree Fock approximation

ECE 5318/6352 Antenna Engineering. Spring 2006 Dr. Stuart Long. Chapter 6. Part 7 Schelkunoff s Polynomial

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.

MODULE FOUR. This module addresses functions. SC Academic Elementary Algebra Standards:

Dispersion Ref Feynman Vol-I, Ch-31

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

3.4 Shrinkage Methods Prostate Cancer Data Example (Continued) Ridge Regression

ANSWER KEY FOR MATH 10 SAMPLE EXAMINATION. Instructions: If asked to label the axes please use real world (contextual) labels

Biplots in Practice MICHAEL GREENACRE. Professor of Statistics at the Pompeu Fabra University. Chapter 13 Offprint

Frequency Response Analysis: A Review

READING STATECHART DIAGRAMS

Lecture 20a. Circuit Topologies and Techniques: Opamps

[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )

2.161 Signal Processing: Continuous and Discrete Fall 2008

T(s) 1+ T(s) 2. Phase Margin Test for T(s) a. Unconditionally Stable φ m = 90 o for 1 pole T(s) b. Conditionally Stable Case 1.

CHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came.

Thermodynamics and Equilibrium

Section 6-2: Simplex Method: Maximization with Problem Constraints of the Form ~

A solution of certain Diophantine problems

A New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation

Physics 2010 Motion with Constant Acceleration Experiment 1

ENGI 4430 Parametric Vector Functions Page 2-01

Hubble s Law PHYS 1301

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

4th Indian Institute of Astrophysics - PennState Astrostatistics School July, 2013 Vainu Bappu Observatory, Kavalur. Correlation and Regression

Lecture 24: Flory-Huggins Theory

제어이론복습 강의보조자료. 박상혁

Algebra2/Trig: Trig Unit 2 Packet

Phys. 344 Ch 7 Lecture 8 Fri., April. 10 th,

Kinematic transformation of mechanical behavior Neville Hogan

ChE 471: LECTURE 4 Fall 2003

Calculus Placement Review. x x. =. Find each of the following. 9 = 4 ( )

Relationship Between Amplifier Settling Time and Pole-Zero Placements for Second-Order Systems *

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

Math Foundations 20 Work Plan

Building to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.

37 Maxwell s Equations

Lecture 12: Chemical reaction equilibria

Tree Structured Classifier

A proposition is a statement that can be either true (T) or false (F), (but not both).

Five Whys How To Do It Better

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India

BASD HIGH SCHOOL FORMAL LAB REPORT

Computational modeling techniques

A Correlation of. to the. South Carolina Academic Standards for Mathematics Precalculus

The standards are taught in the following sequence.

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

COMP 551 Applied Machine Learning Lecture 9: Support Vector Machines (cont d)

and the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Activity Guide Loops and Random Numbers

SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical model for microarray data analysis

Chapter 8. Root Locus Techniques

Transcription:

Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin When using such cmpensatrs in carrying ut ple-placement rt lcus-based design, we ften used the frm s K In the same way, mst tables f Laplace transfrm pairs express a first rder system as G s We s a chse the ntatin in Definitin because the units f thse parameters are, indeed, radians/sec Since in engineering the symbl is cmmnly used t represent frequency in radians/secnd, we have simply chsen t highlight this fact, when addressing these cmpensatrs in the frequency dmain Example Cnsider s K Because, this is a lead cmpensatr In relatin t the rt s lcus, this cmpensatr adds angle (r phase t the rt lcus angle criterin Furthermre, the zer lead s (r ccurs befre the ple n the negative real axis (a Fr K, this is a unity static gain cmpensatr Give its Bde plt Then verlay a straight-line apprximatin The Straight-Line Apprximatin Bde Plt The symmetry in bth plts is n accident The magnitude plt straightline apprximatin has turn-n frequencies At the zer break frequency the slpe changes frm zer t /decade At the ple break frequency the - /decade slpe assciated with this ple turns n, and cancels that f the zer The straight-line phase apprximatin is trickier By itself, the phase f the numeratr f (s ges frm t 9 as : Hence, it has nt nly a turn-n frequency, but als a turn-ff frequency The frmer is at, and the latter is at Because this turn-n/turn-ff interval spans tw decades, the slpe is 45 / decade ver the interval [, ]The denminatr term will clearly have a negative slpe ver the interval [, ] Because these intervals verlap, the slpes will cancel each ther in that regin

(b At what frequency will the phase cmpensatr cntribute the mst phase t the pen lp system At 3rad / s In fact, frm the symmetry f the magnitude plt, we can speculate that the frequency at which the phase is imum is the gemetric mean f [which is equivalent t lg( [lg( lg( ]/ ] (, : (c Estimate the amunt f amplificatin the cmpensatr gives at Then infer a general expressin fr it At 3rad / s the cmpensatr magnitude is abve the static gain f It is als belw the cmpensatr high frequency gain Hence we can infer that at the cmpensatr gain is higher that then static gain, where lg( / / (d Summarize the gain and phase relatins fr the lead cmpensatr lead is imum s K at the frequency where its phase s SUMMARY: The frequency f imum phase lead is, and the amplificatin re: static gain is where lg( / / (e Infer a similar summary fr the lag cmpensatr SUMMARY: The frequency f imum phase lag (ie negative phase is, and the attenuatin re: static gain is / where lg( / (ie this will be a negative value Example Suppse that we have experimentally arrived at the FRF at the right fr a given plant that is t be cntrlled (afr unity feedback prprtinal cntrl with K, estimate the CL GM and PM gc r / s PM 5 and pc r / s PM is undefined (b By incrprating the unity static gain cmpensatr in Example, determine hw will the OL Bde plt be affected at 3rad / s Fr s, the cmpensatr Bde plt gives M ( 3 and ( 3 55 Hence, at this frequency, the OL s magnitude will becme - + = -, and the OL phase will becme -7 + 55 = -5 (c Determine hw much we wuld need t gain up the cmpensatr in (b s that the pen lp gain crssver frequency is at 3 r/s [Nte: this wuld result in a ~75 CL PM] Answer: We wuld need K / K 3 6

3 (darrive at a transfer functin mdel fr the plant, using the plant Bde plt At lw frequencies, the magnitude slpe is -/dec This can be mdeled by a /s term This lpe breaks t a -4/dec value at the frequency 5 r/s Hence, we will mdel the plant as G p s ( s 5 T find the value fr, we will use the lw frequency gain 46/ 46 M p ( M p ( (5 5 Hence, ur plant mdel is G p s ( s 5 (efr the cmpensatr achieved by (b-c, and the plant mdel in (d, use the cmmand margin t verify the design The plt at right verifies the design (fplt the CL step respnse, and cmment n it The plt is at lwer right The CL has unity static gain, achieves steady state in ~ secnd, and has n versht Remark It shuld be nted that we did nt use the mdel in (d t design the cntrl system We used the experimental FRF fr the plant We nly used the mdel as a way f analytically verifying the design, since we d nt have a real plant t implement it n

4 Derivatin f Phase Cmpensatr Design Equatins fr Changing Clsed Lp Phase Margin Here we will derive the mathematical frmulas needed t design a unity static gain phase-lead cmpensatr: s with s ( First, we write the cmpensatr FRF: i ( i i Hence: ( M G ( i and ( tan ( / tan ( / c (afind the frequency at which ( is imum Slutin: OK Here s the first f the tw tricks that I use t slve this prblem: Recall (ha! ha! the trig identity: tan( tan( tan( tan( tan( T use this, let tan ( / and tan ( / We then have tan[ ( ] ( / ( /, and hence: ( / ( / ( tan where d( T find the frequency ( achieves its imum value, we need t slve Nw, n the surface it d may seem like this is a real pain, due t the arctan functins Hwever recall that the arctan is a mntnically increasing functin n the interval ( /, / T see why this matters, recall the chain rule fr differentiatin (the secnd trick: g(x dg[ f ( ] d dg df df d If is strictly mntnic, that means that it derivative is strictly greater than (r less than zer Hence, if in this case df we set the abve equatin t zer, then it must be that Applying this secnd trick, we get d d( d d d ( ( ( ( Fr the rightmst quantity t equal zer, its numeratr must equal zer; that is: ( ( ( ( ( ( ( The slutin t this is = the frequency where ( achieves its imum (b Find the expressin fr the cntrller magnitude at Slutin: This is relatively easy:

5 M ( ( ( Hence: M ( Nw, the static gain f the cntrller is D ( and its high frequency gain is D ( If we cnvert these gains t decibels, then we have D ( lg, D( and Magnitude ( 5 5 6 Bde Diagram D ( i lg In wrds, the cntrller achieves half f its ttal gain at T visualize this, let The Bde plt fr this cntrller is and Phase (deg 3 - - 3 Frequency (rad/sec shwn at right We see that the phase is imum at is, and the assciated gain is half way frm t (cderive the expressin fr ( Slutin: T this end, define / [NOTE: On p349 the term α used by the authrs is /! Hence, my α is the inverse f theirs] We then have ( tan ( tan (/ Applying the abve trig identity gives tan( Nw cnsider the fllwing right triangle at right We see that 4 ( Hence sin(, r sin( sin( We can nw summarize the design prcedure fr using a lead cntrller t increase the pen lp phase by a specified amunt : Having specified sin(, cmpute Since /, by knwing α we knw the rati f the sin( cntrller break frequencies Next, decide at what frequency this phase gain shuld ccur This decisin will ften be made in cnjunctin with ther clsed lp specificatins (eg errr cnstants, etc, and s n ne prcedure can sin( be described here Since and, we can slve these fr the unknwns & : sin( Frm the first, we have / Substituting this int the secnd gives: sin( Having this, we can nw slve fr / sin(

6 Remark The abve develpment is nt given in the bk (it is left as a hmewrk prblem Furthermre, expressins (639 and (64 n p365 invlve the quantities While I d briefly use the variable α, it is nt the same as the authrs On p366 the authrs use the quantities z & p While ne might presume these are the cmpensatr zer and ple, respectively, they are nt As evidenced by (64 they are the negatives f the zer and ple Finally, I chse t use cmpensatr break frequencies design These break frequencies are directly read frm an pen lp Bde plt T d & & fr the simple reasn that we are engaging in frequency dmain-based cntrller Implementatin f the Design Equatins- There are tw standard methds f implementatin Methd : When the PM is specified but n value fr is specified In this case, fr a specified PM ne wuld first identify the current OL gain crssver frequency, current PM The difference PM-PM wuld be raise the OL gain by / And s, except fr the fact that if we center the cmpensatr in rder t read the, it will will n lnger be the gain crssver frequency If we have the freedm t use a nn-unity static gain cmpensatr, then we can simply make its static gain equal t / Then stature as the gain crssver frequency If we have n such freedm (and ften we d nt- since by adjusting the cmpensatr static gain, we wuld alter steady state errr, then we have t resrt t the fllwing apprach: will retain its Identify the frequency at which the OL magnitude is currently /, and hpe that the OL phase has nt changed much in relatin t what it is at Typically, it will be less than what it was at Fr this reasn, we wuld add a cushin t PM PM at the utset The steps fr this methd are summarized belw Cmpute the current PM, and then cmpute the needed PM PM cushin A reasnable first chice fr the cushin is 5 sin( Cmpute an initial value fr sin( 3 Lcate the frequency, O L, where the current pen lp has G ( i / [This is because, if yu were t nw re-center the cmpensatr at, this wuld becme ( new ( 4 Cmpute the pen lp phase ( ld L O L ( O ] 5 Cmpute 8 ( new O L ( If this value is clse t the required PM, missin accmplished If this value is deemed t be unacceptably lw, then return t and add a larger cushin [If this value is deemed t be unacceptably large (which rarely happens!, then return t and reduce the cushin ] 6 Cmpute /, fllwed by (new gc

7 Example 3 [See bk Example 65 n p 367] Cnsider a unity feedback system with pen lp transfer functin G / s( s with FRF at right Suppse that we desire PM 45 withut changing the OL lw frequency gain Using margin(g we get PM 8 We will chse 45 8 5 3 (ie a cushin Then a=(+sind(3/(-sind(3 = 3546, r a=*lg(a =5 3 Fr GO L ( i / 5 5 Using the data cursr, we find 4r / s ( ld 4 Using the data cursr, we then find O L ( 67, hence ( new O ( 67 3 35 L ( new 5 8 O L ( 8 35 45 PM Missin accmplished (in a single try! 6 / 33r / s, fllwed by 758r / s 5 Hence, s s 33 35 The clsed lp step s s 758 respnses fr the cntrllers G C and s 33 G C 35 are shwn at right Clearly, the s 758 inclusin f the PM specificatin resulted in a huge imprvement Befre leaving this example, in rder t better tie this design methd t the rt lcus ple placement methd, we ffer the rt lcus crrespnding t the pen lp system s 33 G K 35 s 758 s( s [Nte: K crrespnds t ur final OLTF] The rt lcus plt at the right reveals that fr K (required t satisfy the steady state errr specificatin, the cmplex clsed lp ples are imally damped [As an aside: Whereas the data cursr states that thse ples crrespnd t % versht, we see frm the abve step respnse, that the actual versht is ~% This increase is due t the clsed lp zer that the phase cmpensatr intrduced] QUESTION: In view f the abve, why didn t we simply design the cntrller via the rt lcus ple-placement methd? It seems like we are just learning a new methd that gives the same results we wuld have gtten using a methd we already knw? ANSWER: Recall that the frequency dmain methd never actually used direct knwledge f the plant transfer functin Even thugh we btained the plant FRF by assuming it was knwn, very ften a plant FRF is btained frm experimental

8 data Hence, the answer is: Because ften we d nt have a mathematical mdel fr the plant transfer functin, the frequency dmain methd can be mre apprpriate As reasnable as the abve answer is, there is yet anther additinal answer that is ften as, if nt mre imprtant T this end, suppse that there is a sec time delay between when the errr is cmputed and when it is received by the cntrller T incrprate this delay int the pen lp, st recall that fr e( t E( s, then e( t t e E( s And s nw the final pen lp transfer functin is: G( s s 33 35 e s 758 s( s s The clsed lp unit step respnses withut and with this time delay included are shwn at right Clearly, the time delay had a destabilizing effect The term in the pen lp transfer functin means that it is n lnger a rati f plynmials, frm which ne simply btains ples and zers Hence, nw, the rt lcus methd des nt even apply i t Even s, this term is trivial t incrprate int the pen lp FRF This is because e has a magnitude f (and s des nthing t the pen lp magnitude, and cntributes angle t And s, at the new gain e st t ( crssver frequency 4r / s this time delay will reduce the pen lp phase by t ( (4(8 / 4 Hence, when the delay is included, the clsed lp PM 45 4 OUCH! Basically, we are back t where we were prir t incrprating the phase cmpensatr Yes, that IS true Hwever, had that cmpensatr NOT been incrprated, the clsed lp PW wuld be PM 8 4 6 In layman s terms, we wuld have been SOL Methd : When the PM is specified and als a value fr is specified In sme ways this case is easier that Methd It typically requires the cmbined use f a lead and a lag cmpensatr The lead cmpensatr is use t adjust the phase at Then the lag cmpensatr is used as a high frequency attenuatr t frce the OL gain t be at cmpensatr needs t be used as a high frequency amplifier] [Nte: It can als happen that, as ppsed t a lag cmpensatr, a secnd lead Example 4 Cnsider again the plant in Example 5: / s( s The Bde plt is given again here fr cnvenience G p

9 Suppse that we require that the clsed lp have PM 6 at an pen lp gain crssver frequency 58 r gc / will first design a unity static gain lead cmpensatr with and r / s The abve design equatins give: & 8, and s G lead (8/ 5( s 5 /( s 8 The resulting OL Bde plt is shwn belw 5 c s We We see that the phase is what we want at r / s, but the OL gain is abut - S we will use a secnd lead cmpensatr as a high frequency amplifier G lead c ( / ( s /( s