Stability analysis of feedback systems
|
|
- Gloria Hopkins
- 6 years ago
- Views:
Transcription
1 Chapter 1 Stability analysis of feedback systems 1.1 Introduction A control system must be asymptotically stable. A method for determination of the stability property of a control system will be described in the following. The method is based on the frequency response ofthelooptransferfunction 1,L(jω),anditisdenotedNyquist sstability criterion. This is a graphical analysis method. There are also algebraic analysis methods, as the Routh s stability criterion which is based on the coefficients of the characteristic polynomial of the control system. Such stability analysis methods are not described in this book since they have limited practical importance(and the mathematical operations become quite complicated except for the simplest models). 1.2 Stability analysis of dynamic systems in general This section defines the different stability properties that a dynamic systemcanhaveintermsofimpulseresponseofthesystem. Then,the corresponding transfer function pole locations in the complex plane are derived. In Section 1.3, and subsequent sections, these results are applied loop. 1 The loop transfer function is the product of all the transfer functions in the closed 1
2 2 to feedback(control) systems, which are a special case of dynamic systems Stability properties and impulse response The different stability properties can be defined in several ways. I have chosentousetheimpulseresponse ofthesystemasthebasisfordefinition of the stability properties. The impulse response is the time-response in outputvariableofthesystemduetoanimpulseontheinput. Usingthe impulse response makes it relatively simple to relate stability properties to thepoles ofthesystem(thisisbecausetheimpulseresponseandthepoles are closely connected), as you will see soon. Some words about the impulse signal: It is a time-signal which in principle has infinitely short duration and infinite amplitude, but so that the integral ofthesignal theintegralisequaltotheareaunderthetime-functionof thesignal isfinite. Thisareaisalsocalledthestrength oftheimpulse. Animpulseofstrengthoneiscalledaunitimpulse,δ(t). Thesquarepulse infigure1.1approachesanimpulsefunctionas goestozero. A t Figure1.1: Thesquarepulseapproachesanimpulsefunctionas goestozero. Herearethestabilitydefinitions: Adynamicsystemhasoneofthe following stability properties: Asymptotically stable system: The stationary impulse response,
3 3 h(t), is zero: lim h(t)=0 (1.1) t Marginally stable system: The stationary impulse response is different from zero, but limited: 0< lim t h(t)< (1.2) Unstable system: The stationary impulse response is unlimited: lim h(t)= (1.3) t A system is stable if it is either asymptotically stable or marginally stable. Figure 1.2 depicts the different stability properties. Asymtotically stable system h(t) 0 t Impulse, System Impulse response y(t)=h(t) Marginally stable system 0 t Unstable system Figure 1.2: Different stability properties Inthedefinitionaboveitisassumedthattheimpulseintheinputustarts fromzeroandthattheimpulseresponseintheoutputyhaszeroinitial value. Inpracticetheseinitialvaluesmaybedifferentfromzeroifthe system initially is some other operating point than the zero operating point.
4 4 Oneproblemwiththeidealimpulsefunctionisitcannotbegenerated fullyinpractice,butinpracticethereishardlyeveraneedtoperform impulse response experiments to determine stability properties. It is more useful as a conceptual definition of stability, cf. the next section Stability properties and poles Inmanycasesitwouldbequiteimpracticaliftheonlywaytodetermine thestabilitypropertyofasystemwastodoexperimentsortorun simulations, to obtain the impulse response(or the step response). Fortunately, we can conclude about the stationary part of the impulse response, and hence conclude about the stability property, by just analyzing the mathematical model of the system. This is because the stationary impulse response is a function of the poles of the system. The connection between the impulse response and the poles is derived in the following. Letusassumethatthemathematicalmodelofthesystemisatransfer function, H(s), from input signal u to output signal y: y(s) = H(s)u(s) (1.4) Theinputuwillbeaunitimpulse,thatis,u(t)=δ(t). Itcanbeshown thatthelaplacetransformofδ(t)is1. Letusdenotetheimpulseresponse h(t). The Laplace transform of h(t) is h(s)=h(s)l{δ(t)}=h(s) 1=H(s) (1.5) Thus, the Laplace transform of the impulse response equals the transfer function of the system. We need the impulse response time-function h(t) since it is the basis for the definitions of the different stability properties. With results from Laplace transform theory, h(t) can be calculated using the following formula: h(t) = 1 d m 1 lim s p i (m 1)! i ds m 1 (s p i ) m H(s) e st (1.6) }{{} h(s) = lim s p i i (s p i)h(s) e st ifm=1 (1.7) }{{} {p i }isthesetofpolesinh(s),andhencetherootsofthedenominatorof H(s). m is the multiplicity of the poles(so-called simple poles have h(s)
5 5 m = 1). The denominator polynomial of H(s) is called the characteristic polynomial,, a(s). The poles are the roots of a(s). Consequently, the poles are the solutions of the characteristic equation An example: Transfer function H(s)= a(s) = 0 (1.8) 1 (s+2)(s+3) hasthepolesp 1 = 2andp 2 = 3,andthecharacteristicpolynomialis (1.9) a(s)=(s+2)(s+3)=s 2 +5s+6 (1.10) We will now use(1.6) and(1.7) to connect the different stability properties tothepoleplacementinthecomplexplane. Letusfirstassumethatthe polesofh(s)aresimple. Thenm=1,andh(t)isgivenby(1.7). Apoleis generally a complex number: p i =a i +jb i (1.11) wherea i istherealpartandb i istheimaginarypartofthepole. (1.7) impliesthath(t)isequaltothesumoftheipartialresponsesofthetotal impulse response: h(t)= h i (t) (1.12) i where h i (t)=k i e p it =k i e (a i+jb i )t =k i e a it e jb it (1.13) Herek i issomeconstant. Theterme jb it isacomplexnumberontheunity circleandthereforeithasabsolutevalueequalto1. 2 Thus,itistheterm e a it whichdeterminesthesteady-state(t )absolutevalueofthe partialresponseh i (t)accordingtothefollowinganalysis: Supposethattherealpart,a i,ofthepoleisstrictlynegative,thatis, a i <0,whichmeansthatthepoleliesinthelefthalfplane. This impliese a it 0,andthereforeh i (t) 0ast. Supposethattherealpart,a i,ofthepoleiszero,thatisa i =0, whichmeansthatthepoleliesontheimaginaryaxis. Thisimplies e a it =1,andthereforeh i (t)goestowardsaconstantvaluedifferent fromzeroast. 2 Ifapolehasaimaginarypartbdifferenfromzero,theremustbeacomplexconjugate pole with imaginary part b. However, this fact does not influate the conclusions of the analysis we are performing.
6 6 Supposethattherealpart,a i,ofthepoleisstrictlypositive,thatis, a i >0,whichmeansthatthepoleliesintheleftrightplane. This impliese a it,andthereforeh i (t) ast. From the above analysis we can conclude as follows for transfer functions havingpolemultiplicityone: (1)Ifeachofthepolesliesinthelefthalf plane, the system is asymptotically stable, because then each of the partial impulseresponseterms,h i (t),goestowards0ast. (2)Ifonepolelies ontheimaginaryaxiswhiletherestofthepolesliesonthelefthalfplane, thesystemismarginallystable,becausethenoneoftheh i (t)-termsgoes towardsaconstantvaluedifferentfrom0ast. (3)Ifatleastoneof thepolesliesintherighthalfplane,thesystemisunstable,becausethen atleastonetermh i (t)goesto ast. Multiple poles: It would have been nice to conclude about stability and poleplacementnow,butwehavetolookcloseratthecaseofmultiple polesofh(s). Theimpulseresponseh(t)isgivenby(1.6). Supposethat thatthemultiplicityofthepolep i ism=2. Thecorrespondingpartial impulse response becomes { d [ h i (t)= lim (s p i ) 2 H(s)e st]} (1.14) s pi ds Here,theterm d ds (est )isequaltote st,whichmeansthath i (t)willcontain termsaste p it. Byperformingthesameanalysisasforsimplepoles,wewill findthefollowing: (1)h i (t) 0forapolewithnegativerealpart(since te p it goestowardszerobecausee p it decreasesfasterthantincreases). (2) h i (t) forapoleontheimaginaryaxis(te p it equalst). (3)h i (t) forapolehavingpositiverealpart. Wewillgetthesameresultsifthe multiplicity m is greater than two. Now we can conclude by stating the following correspondence between stability and pole placement: Asymptotically stable system: Each of the poles of the transfer function lies strictly in the left half plane(has strictly negative real part). Marginallystable system: Oneormorepolesliesonthe imaginaryaxis(haverealpartequaltozero),andallthesepolesare distinct. Besides, no poles lie in the right half plane. Unstable system: Atleastonepoleliesintherighthalfplane(has realpartgreaterthanzero). Or: Therearemultiplepolesonthe imaginary axis.
7 7 Figure 1.3 gives a illustration of the relation between stability property and pole placement. Left half plane Im Right half plane Asymptotically stable pole area Unstable pole area Re Figure 1.3: The relation between stability property and pole placement Example 1.1 Stability property of some simple dynamic systems The first order transfer function H 1 (s)= 1 s+1 haspolep= 1whichliesinthelefthalfplane. Thus,thesystemis asymptotically stable. (1.15) The transfer function H 2 (s)= 1 s (1.16) (whichisthetransferfunctionofanintegrator)haspolep=0,whichlies ontheimaginaryaxisandhasmultiplicityone. So,thesystemis marginally stable. The transfer function H 3 (s)= 1 s 2 (1.17) havepolesp=0,whichisontheimaginaryaxiswithmultiplicitytwo. The system is therefore unstable. The transfer function H 4 (s)= 1 s 1 (1.18)
8 8 haspolep=+1,whichliesintherighthalfplane. Thesystemistherefore unstable. [End of Example 1.1] Example 1.2 Stability property of a mass-spring-damper Figure 1.4 shows a mass-spring-damper-system.y is position. F is applied F [N] m K [N/m] D [N/(m/s)] 0 y [m] Figure 1.4: Mass-spring-damper force. D is damping constant. K is spring constant. Newton s 2. Law gives the following mathematical model: mÿ(t)=f(t) Dẏ(t) Ky(t) (1.19) ThetransferfunctionfromtheforceF topositionyis H(s)= y(s) F(s) = 1 ms 2 +Ds+K (1.20) Assumethatm=20kg,D=4N/(m/s),andK=2N/m. Whatisthe stability property of the system? The characteristic polynomial becomes which has roots a(s)=ms 2 +Ds+K=20s 2 +4s+2 (1.21) p 1, p 1 = 4± = 0.1 ± j0.3 (1.22) which are the poles of H(s). Both these poles have strictly negative real parts( 0.1). The system is therefore asymptotically stable. Figure 1.5 shows the poles(marked as crosses) in the complex plane. Figure 1.6 shows the impulse response. AssumenowthatthedamperisremovedsothatD=0. Thenthe characteristic polynomial is a(s)=ms 2 +K (1.23)
9 9 Figure 1.5: The poles of a mass-spring-damper plotted in the complex plane. Thepolesarep 1,2 = 0,1±j0,3. Figure 1.6: The impulse response for a mass-spring-damper with m = 20 kg, D=4N/(m/s)andK f =2N/m. Thesystemetisasymptoticallystable. andthepolesare K p 1, p 1 =±j =±j0.32 (1.24) m which lies on the imaginary axis, and they have multiplicity one. The system is then marginally stable. Figure 1.7 shows the impulse response. [End of Example 1.2]
10 10 Figure 1.7: The impulse response of the mass-spring-damper with m = 20 kg, D=0N/(m/s)ogK f =2N/m. Thesystemismarginallystabile. 1.3 Transfer functions used for stability analysis of feedback systems Which is the transfer function to be used to determine the stability analysis of feedback systems, e.g. control systems? Figure 1.8 shows a block diagram(transfer function based block diagram) of a control system. Wemustselectatransferfunctionfromoneoftheinputsignalstothe closedlooptooneoftheoutputsignalsfromtheloop. Letusselectthe transferfunctionfromthesetpointy msp totheprocessmeasurementy m. This transfer function is y m (s) y msp (s) = H c(s)h u (s)h m (s) 1+H c (s)h u (s)h m (s) = L(s) =T(s) (1.25) 1+L(s) which is the tracking transfer function of the control system. (If we had selected some other transfer function, for example the transfer function from the disturbance to the process output variable, the result of the analysis would have been the same.) L(s) is the loop transfer function of the control system: L(s)=H c (s)h u (s)h m (s) (1.26) Figure1.9showsacompactblockdiagramofacontrolsystem. The transferfunctionfromy msp toy m isthetrackingfunction: s T(s)= L(s) 1+L(s) = n L (s) d L (s) 1+ n L(s) d L (s) = n L (s) d L (s)+n L (s) (1.27) wheren L (s)andd L (s)arethenumeratoranddenominatorpolynomialsof L(s), respectively. The characteristic polynomial of the tracking function is c(s)=d L (s)+n L (s) (1.28)
11 11 Process disturbance Process d(s) Reference or setpoint y r (s) H sm (s) Reference in measurement unit y mr (s) Model of sensor with scaling Manipulating Control error variable e m (s) Hc (s) u(s) Controller Measurement y m (s) H u (s) H d (s) Actuator transfer function H s (s) Sensor (measurement) with scaling Disturbance transfer function y(s) Process output variable Figure 1.8: Block diagram(transfer function based) of a control system The stability of the control system is determined by the placement of the roots of(1.28) in the complex plane. Example 1.3 Stability of a feedback control system See Figure 1.8 which shows a block diagram(transfer function based) of a control system. We assume that the individual transfer functions with parameter values are as follows: H sm (s)=k sm =1 (1.29) y msp L(s) y m Figure1.9: Compactblockdiagramofacontrolsystemwithsetpointy msp as inputvariableandprocessmeasurementy m asoutputvariable
12 K p =1(asympt. stable) K p =2(marg. stable) K p =4(unstable) p 1 = 1.75 p 1 = 2 p 1 = 2.31 p 2 = 0.12+j0.74 p 2 =j p 2 =0.16+j1.31 p 3 = 0.12 j0.74 p 3 = j p 3 =0.16 j1.31 Table1.1: PolesofthetrackingtransferfunctionforvariousK p -values 12 H s (s)=k s =1 (1.30) H c (s)=k p (proportionalcontroller) (1.31) 1 H u (s)= (s+1) 2 (1.32) s 1 H d (s)= (s+1) 2 (1.33) s The stability property of the control system can be determined from the placementofthepolesofthetrackingtransferfunction,t(s),whichisthe transferfunctionfromthereferencey SP totheprocessoutputvariabley. The tracking transfer function is Inserting(1.29) (1.33) gives T(s)= y(s) y SP (s) = H sm(s)h c (s)h u (s) 1+H s (s)h c (s)h u (s) H yr,y(s)= The characteristic polynomial of T(s) is (1.34) K p s 3 +2s 2 +s+k p (1.35) a(s)=s 3 +2s 2 +s+k p (1.36) Table1.1showsthepolesforthreedifferentK p -values. 3 Figure1.10showsthestepresponseinyforthethreeK p -values(itisa unitstepinr). [End of Example 1.3] 1.4 Nyquist s stability criterion The Nyquist s stability criterion will now be derived. We start with a little rewriting: Therootsof(1.28)arethesameastherootsof d L (s)+n L (s) =1+ n L(s) =1+L(s)=0 (1.37) d L (s) d L (s) 3 The poles can be calculated using the MATLAB-function roots or pzmap or pole or using the LabVIEW-function Complex Polynomial Roots.
13 13 Figure 1.10: Example 1.3: Step response in the process output variable y for threedifferentk p -values which, therefore, too can be denoted the characteristic equation of the control system. (1.37) is the equation from which the Nyquist s stability criterion will be derived. In the derivation we will use the Argument variation principle: Argument variationprinciple: Givenafunctionf(s)wheresisa complexnumber. Thenf(s)isacomplexnumber,too. Aswithall complexnumbers,f(s)hasanangleorargument. Ifsfollowsa closed contour Γ(gamma) in the complex s-plane which encircles a numberofpolesandanumberofzerosoff(s),seefigure1.11,then the following applies: arg Γ f(s)=360 (numberofzerosminusnumberofpolesoff(s)insideγ) (1.38) wherearg Γ f(s)meansthechangeoftheangleoff(s)whenshas followed Γ once in positive direction of circulation(i.e. clockwise).
14 14 Γ contour Im(s) Right half plane Re(s) P CL poles of closed loop system in right half plane P OL poles of open loop system in right half plane Positive direction of circulation Figure 1.11: s shall follow the Γ contour once in positive direction (counter clockwise). Forourpurpose,weletthefunctionf(s)intheArgumentvariation principle be f(s)=1+l(s) (1.39) TheΓcontourmustencircletheentirerighthalfs-plane,sothatweare certainthatallpolesandzerosof1+l(s)areencircled. Fromthe Argument Variation Principle we have: arg Γ d L (s)+n L (s) [1+L(s)] = arg Γ d L (s) = 360 (numberofrootsof(d L +n L )inrhp (1.40) minusnumberrootsofd L inrhp) (1.41) = 360 (numberpolesofclosedloopsysteminrhp minusnumberpolesofopensysteminrhp) = 360 (P CL P OL ) (1.42) whererhpmeansrighthalfplane. By opensystem wemeanthe (imaginary)systemhavingtransferfunctionl(s)=n L (s)/d L (s),i.e.,the
15 15 original feedback system with the feedback broken. The poles of the open systemaretherootsofd L (s)=0. Finally, we can formulate the Nyquist s stability criterion. But before we dothat,weshouldremindourselveswhatweareafter,namelytobeable todeterminethenumberpolesp CL oftheclosedloopsysteminrhp.it those poles which determines whether the closed loop system(the control system)isasymptoticallystableornot. If P CL =0theclosedloopsystem is asymptotically stable. Nyquist sstabilitycriterion: LetP OL bethenumberofpolesofthe opensystemintherighthalfplane,andletarg Γ L(s)betheangular changeofthevectorl(s)asshavefollowedtheγcontouroncein positivedirectionofcirculation. Then,thenumberpolesP CL ofthe closedloopsystemintherighthalfplane,is P CL = arg ΓL(s) 360 +P OL (1.43) IfP CL =0,theclosedloopsystemisasymptoticallystable. Letustakeacloserlookatthetermsontherightsideof(1.43): P OL are therootsofd L (s),andthereshouldnotbeanyproblemcalculatingthese roots. Todeterminetheangularchangeofthevector1+L(s). Figure1.12 showshowthevector(orcomplexnumber)1+l(s)appearsinanyquist diagram foratypicalplotofl(s). ANyquistdiagramissimplya CartesiandiagramofthecomplexplaneinwhichLisplotted. 1+L(s)is thevectorfromthepoint( 1, 0j),whichisdenotedthecriticalpoint,to the Nyquist curve of L(s). More about the Nyquist curve of L(jω) LetustakeamoredetailedlookattheNyquistcurveofLassfollowsthe Γ contour in the s-plane, see Figure In practice, the denominator polynomial of L(s) has higher order than the numerator polynomial. This impliesthatl(s)ismappedtotheoriginofthenyquistdiagramwhen s =. Thus,thewholesemicircularpartoftheΓcontourismappedto the origin. TheimaginaryaxisconstitutestherestoftheΓcontour. Howisthe mappingofl(s)assrunsalongtheimaginaryaxis? Ontheimaginary axiss=jω,whichimpliesthatl(s)=l(jω),whichisthefrequency responseofl(s). Aconsequenceofthisisthatwecaninprinciple
16 16 The critical point 1 Im L(s) Negative ω Infinite ω 0 Re L(s) 1 + L(s) Nyquist curve of L(s) Positive ω Decreasing ω Figure1.12: TypicalNyquistcurveof L(s). Thevector1+L(s)isdrawn. determine the stability property of a feedback system by just looking at the frequency response of the open system, L(jω). ωhasnegativevalueswhens=jωisonthenegativeimaginaryaxis. For ω < 0 the frequency response has a mathematical meaning. From general properties of complex functions, L( jω) = L(jω) (1.44) and L( jω) = L(jω) (1.45) ThereforetheNyquistcurveofL(s)forω<0willbeidenticaltothe Nyquistcurveofω>0,butmirroredabouttherealaxis. Thus,weonly needtoknowhowl(jω)ismappedforω 0. TherestoftheNyquist curvethencomesbyitself! ActuallyweneednotdrawmoreoftheNyquist curve(forω>0)thanwhatissufficientfordeterminingifthecritical point is encircled or not. WemustdosomeextraconsiderationsifsomeofthepolesinL(s),which arethepolesoftheopenloopsystem,lieintheorigin. Thiscorrespondsto pure integrators in control loop, which is a common situation in feedback control systems because the controller usually has integral action, as in a PI or PID controller. If L(s) contains integrators, the Γ contour must go outside theorigo. Buttotheleftortotheright? Wechoosetotheright, seefigure1.13. (Wehavetherebydecidedthattheoriginbelongstothe
17 17 Im(s) Infinetely small radius Im L(s) Infinitely large radius from origin 1 0 Re(s) 0 Re L(s) 1 + L(s) Figure 1.13: Left diagram: If L(s) has a pole in origin, the Γ contour must pass the origin along an arbitrarily small semicircle to the right. Right diagram: AtypicalNyquistcurveof L. lefthalfplane. ThisimpliesthatP OL doesnotcountthesepoles.) The radius of the semicircle around origin is arbitrarily small. The Nyquist curvethenbecomesasshowninthediagramtotherightinthesame figure. The arbitrarily small semicircle in the s-plane is mapped to an infinitelylargesemicircleinthel-plane. Theisbecauseass 0,theloop transfer function is approximately L(s) K s (ifweassumeonepoleintheorigin). Onthesmallsemicircle, s=re jθ (1.46) which gives L(s) K r e jθ (1.47) Whenr 0andwhensimultaneouslyθgoesfrom+90 via0 to 90, the Nyquist curve becomes an infinitely large semicircle, as shown. The Nyquist s stability criterion for non-rational transfer functions The Nyquist s stability criterion gives information about the poles of feedback systems. So far it has been assumed that the loop transfer
18 18 function L(s) is a rational transfer function. What if L(s) is irrational? Here is one example: L(s)= 1 s e τs (1.48) wheree τs representstimedelay. InsuchcasesthetrackingratioT(s)will alsobeirrational,andthedefinitionofpolesdoesnotapplytosuch irrational transfer functions. Actually, the Nyquist s stability criterion can be used as a graphical method for determining the stability property on basis of the frequency response L(jω). Nyquist s special stability criterion Inmostcasestheopensystemisstable,thatis,P OL =0. (1.43)then becomes P CL = arg Γ[L(s)] 360 (1.49) This implies that the feedback system is asymptotically stable if the Nyquist curve does not encircle the critical point. This is the Nyquist s special stability criterion or the Nyquist s stability criterion for open stable systems. The Nyquist s special stability criterion can also be formulated as follows: The feedback system is asymptotically stable if the Nyquist curve of L has thecriticalpointonitsleftsideforincreasing ω. Another way to formulate Nyquist s special stability criterion involves the amplitudecrossoverfrequency ω c andthephasecrossoverfrequency ω 180. ω c isthefrequencyatwhichthel(jω)curvecrossestheunitcircle,while ω 180 isthefrequencyatwhichthel(jω)curvecrossesthenegativereal axis. In other words: L(jω c ) =1 (1.50) and argl(jω 180 )= 180 (1.51) See Figure Note: The Nyquist diagram contains no explicit frequency axis. We can now determine the stability properties from the relation between these two crossover frequencies: Asymptoticallystableclosedloopsystem:ω c <ω 180 Marginallystableclosedloopsystem:ω c =ω 180 Unstableclosedloopsystem:ω c >ω 180
19 19 Im L(s) L(jω 180 ) j Unit circle 1 0 Re L(s) L(jω c ) Decreasing ω Positive ω Figure1.14: Definitionofamplitudecrossoverfrequencyω c andphasecrossover frequencyω 180 The frequency of the sustained oscillations There are sustained oscillations in a marginally stable system. The frequencyoftheseoscillationsis ω c =ω 180.Thiscanbeexplainedas follows: Inamarginallystablesystem,L(±jω 180 )=L(±jω c )= 1. Therefore,d L (±jω 180 )+n L (±jω 180 )=0,whichisthecharacteristic equationoftheclosedloopsystemwith±jω 180 insertedfors. Therefore, thesystemhas±jω 180 amongitspoles. Thesystemusuallyhave additionalpoles,buttheylieinthelefthalfplane. Thepoles±jω 180 leads tosustainedsinusoidaloscillations. Thus,ω 180 (orω c )isthefrequencyof the sustained oscillations in a marginally stable system. 1.5 Stability margins StabilitymarginsintermsofgainmarginGM and phasemarginpm An asymptotically stable feedback system may become marginally stable if the loop transfer function changes. The gain margin GM and the phase margin PM [radiansordegrees]arestabilitymargins whichintheirown
20 20 ways expresses how large parameter changes can be tolerated before an asymptotically stable system becomes marginally stable. Figure 1.15 shows the stability margins defined in the Nyquist diagram. GM is the Im L(s) L(jω 180 ) 1/GM j Unity circle 1 PM 0 Re L(s) L(jω c ) Figure1.15: GainmarginGM andphasemarginpm definedinthenyquist diagram (multiplicative, not additive) increase of the gain that L can tolerate at ω 180 beforethelcurve(inthenyquistdiagram)passesthroughthe critical point. Thus, L(jω 180 ) GM =1 (1.52) which gives GM= 1 L(jω 180 ) = 1 ReL(jω 180 ) (Thelatterexpressionin(1.53)isbecauseatω 180,ImL=0sothatthe amplitudeisequaltotheabsolutevalueoftherealpart.) (1.53) Ifweusedecibelastheunit(likeintheBodediagramwhichwewillsoon encounter), then GM [db]= L(jω 180 ) [db] (1.54) ThephasemarginPM isthephasereductionthatthelcurvecantolerate atω c beforethelcurvepassesthroughthecriticalpoint. Thus, argl(jω c ) PM= 180 (1.55) which gives PM =180 +argl(jω c ) (1.56)
21 21 Wecannowstateasfollows: Thefeedback(closed)systemis asymptotically stable if GM>0dB =1andPM >0 (1.57) This criterion is often denoted the Bode-Nyquist stability criterion. Reasonable ranges of the stability margins are and 2 6dB GM 4 12dB (1.58) 30 PM 60 (1.59) Thelargervalues,thebetterstability,butatthesametimethesystem becomes more sluggish, dynamically. If you are to use the stability margins as design criterias, you can use the following values(unless you have reasons for specifying other values): GM 2.5 8dBandPM 45 (1.60) Forexample,thecontrollergain,K p,canbeadjusteduntiloneofthe inequalitiesbecomesanequality. 4 Itcanbeshown 5 thatforpm 70,thedampingofthefeedbacksystem approximately corresponds to that of a second order system with relative damping factor ζ PM 100 (1.61) Forexample,PM =50 ζ= Stability margins in terms of maximum sensitivity amplitude, S(jω) max An alternative quantity of a stability margin, is the minimum distance fromthel(jω)curvetothecriticalpoint. Thisdistanceis 1+L(jω),see Figure1.16. So,wecanusetheminimalvalueof 1+L(jω) asastability margin. However,itismorecommontotaketheinverseofthedistance: Thus,astabilitymarginisthemaximum valueof1/ 1+L(jω). And 4 Butyoushoulddefinitelycheckthebehaviourofthecontrolsystembysimulation,if possible. 5 The result is based on the assumption that the loop transfer function is L(s) = ω 2 0/[(s+2ζω 0)s] which gives tracking transfer function T(s) = L(s)/[1 + L(s)] = ω 2 0/ [ s 2 +2ζω 0 s+ω0] 2. ThephasemarginPM canbecalculatedfroml(s).
22 22 Im L(s) 1 0 Re L(s) 1+L(jω c ) min = S(jω c ) max L(jω) Figure 1.16: The minimum distance between the L(jω) curve and the critical point can be interpreted as a stability margin. This distance is 1+L min = S max. since1/[1+l(s)]isthesensitivityfunctions(s),then S(jω) max represents a stability margin. Reasonable values are in the range dB S(jω) max dB (1.62) Ifyouuse S(jω) max asacriterionforadjustingcontrollerparameters,you can use the following criterion(unless you have reasons for some other specification): S(jω) max =2.0 6dB (1.63) 1.6 Stability analysis in a Bode diagram ItismostcommontouseaBodediagramforfrequencyresponsebased stability analysis of closed loop systems. The Nyquist s Stability Criterion says: The closed loop system is marginally stable if the Nyquist curve(of L)goesthroughthecriticalpoint,whichisthepoint( 1,0). Butwhereis the critical point in the Bode diagram? The critical point has phase (angle) 180 andamplitude1=0db.thecriticalpointtherefore constitutestwolinesinabodediagram: The0dBlineintheamplitude diagramandthe 180 lineinthephasediagram. Figure1.17shows typical L curves for an asymptotically stable closed loop system. In the figure,gm,pm,ω c andω 180 areindicated. Example 1.4 Stability analysis of a feedback control system
23 23 [db] L(jω) 0 db ω c GM ω (logarithmic) [degrees] arg L(jω) -180 PM ω 180 Figure 1.17: Typical L curves of an asymptotically stable closed loop system withgm,pm,ω c andω 180 indicated Given a feedback control system with structure as shown in Figure The loop transfer function is 1 L(s)=H c (s)h p (s)= K p }{{}(s+1) 2 = s }{{} H c(s) H p(s) K p (s+1) 2 s =n L(s) d L (s) (1.64) We will determine the stability property of the control system for different valuesofthecontrollergaink p inthreeways: Poleplacement,Nyquist s Stability Criterion, and simulation. The tracking transfer function is T(s)= y m(s) y msp (s) = L(s) 1+L(s) = n L (s) d L (s)+n L (s) = The characteristic polynomial is K p s 3 +2s 2 +s+k p (1.65) c(s)=s 3 +2s 2 +s+k p (1.66) Figures showthestepresponseafterastepinthesetpoint,the poles,thebodediagramandnyquistdiagramforthreek p valueswhich result in different stability properties. The detailed results are shown below.
24 24 y msp Controller H c (s) Process with measurement and scaling H p (s) y m Figure 1.18: Example 1.4: Block diagram of feedback control system K p =1: Asymptoticallystablesystem,seeFigure1.19. Fromthe BodediagramwereadoffstabilitymarginsGM=6.0dB=2.0and PM =21. weseealsothat S(jω) max =11dB=3.5(alargevalue, butitcorrespondswiththesmallthephasemarginofpm =20 ). K p =2: Marginallystablesystem,seeFigure1.20. FromtheBode diagram,ω c =ω 180. TheLcurvegoesthroughthecriticalpointin thenyquistdiagram. S max hasinfinitelylargevalue(sincethe minimumdistance,1/ S max,between L andthecriticalpointis zero). LetuscalculatetheperiodT p oftheundampedoscillations: Since ω 180 =1.0rad/s,theperiodisT p =2π/ω 180 =6.28s,whichfitswell with the simulation shown in Figure K p =4: Unstablesystem,seeFigure1.21. FromtheBodediagram, ω c >ω 180. FromtheNyquistdiagramweseethattheLcurvepasses outside the critical point. (The frequency response curves of M and N havenophysicalmeaninginthisthecase.) [End of Example 1.4] 1.7 Robustness in term of stability margins Per definition the stability margins expresses the robustness of the feedback control system against certain parameter changes in the loop transfer function: Thegainmargin GM ishowmuchtheloopgain,k,canincrease beforethesystembecomesunstable. Forexample,isGM =2when
25 25 K=1.5,thecontrolsystembecomesunstableforK largerthan 1.5 2=3.0. Thephasemargin PM ishowmuchthephaselagfunctionofthe loop can be reduced before the loop becomes unstable. One reason of reducedphaseisthatthetimedelayincontrolloopisincreased. A changeofthetimedelayby τ introducesthefactore τs inl(s) andcontributestoarglwith τ ω[rad]or τ ω 180 π [deg]. L ishowevernotinfluencedbecausetheamplitudefunctionofe τs is1, independentofthevalueofτ. Thesystembecomesunstableifthe timedelayhaveincreasedby τ max suchthat 6 PM= τ max ω c 180 π [deg] (1.67) which gives the following maximum change of the time delay: τ max = PM ω c π 180 (1.68) IfyouwanttocalculatehowmuchthephasemarginPM isreduced ifthetimedelayisincreasedby τ,youcanusethefollowing formula which stems from(1.67): PM = τ ω c 180 π [deg] (1.69) For example, assume that a given control system has ω c =0.2rad/minandPM=50. Ifthetimedelayincreasesby1min, thephasemarginisreducedby PM= π =11.4,i.e. from 50 to RememberthatPM isfoundatω c.
26 Figure 1.19: Example 1.4: Step response(step in setpoint), poles, Bode diagramandnyquistdiagramwithk p =1. Thecontrolsystemisasymptotically stable. 26
27 Figure 1.20: Example 1.4: Step response (step in setpoint), poles, Bode diagram and Nyquist diagram with K p = 2. The control system is marginally stable. 27
28 Figure 1.21: Example 1.4: Step response(step in setpoint), poles, Bode diagramandnyquistdiagramwithk p =4. Thecontrolsystemisunstable. 28
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 informationMAS107 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 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 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 informationOutline. 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 information7.2 Controller tuning from specified characteristic polynomial
192 Finn Haugen: PID Control 7.2 Controller tuning from specified characteristic polynomial 7.2.1 Introduction The subsequent sections explain controller tuning based on specifications of the characteristic
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 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 informationLinear 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 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 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 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 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 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 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 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 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 informationChapter 2. Classical Control System Design. Dutch Institute of Systems and Control
Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral
More informationControl of Electromechanical Systems
Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance
More informationDr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review
Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More 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 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 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 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 informationx(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 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 informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationSystems 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 informationClassify 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 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 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 informationLecture 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 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 informationControl Systems. Root Locus & Pole Assignment. L. Lanari
Control Systems Root Locus & Pole Assignment L. Lanari Outline root-locus definition main rules for hand plotting root locus as a design tool other use of the root locus pole assignment Lanari: CS - Root
More informationAdditional Closed-Loop Frequency Response Material (Second edition, Chapter 14)
Appendix J Additional Closed-Loop Frequency Response Material (Second edition, Chapter 4) APPENDIX CONTENTS J. Closed-Loop Behavior J.2 Bode Stability Criterion J.3 Nyquist Stability Criterion J.4 Gain
More informationPrüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid
More informationIC6501 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 informationME 304 CONTROL SYSTEMS Spring 2016 MIDTERM EXAMINATION II
ME 30 CONTROL SYSTEMS Spring 06 Course Instructors Dr. Tuna Balkan, Dr. Kıvanç Azgın, Dr. Ali Emre Turgut, Dr. Yiğit Yazıcıoğlu MIDTERM EXAMINATION II May, 06 Time Allowed: 00 minutes Closed Notes and
More informationECSE 4962 Control Systems Design. A Brief Tutorial on Control Design
ECSE 4962 Control Systems Design A Brief Tutorial on Control Design Instructor: Professor John T. Wen TA: Ben Potsaid http://www.cat.rpi.edu/~wen/ecse4962s04/ Don t Wait Until The Last Minute! You got
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 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 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 informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 5. 2. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationSTABILITY ANALYSIS TECHNIQUES
ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response
More informationAPPLICATIONS FOR ROBOTICS
Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table
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 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 informationLABORATORY INSTRUCTION MANUAL CONTROL SYSTEM I LAB EE 593
LABORATORY INSTRUCTION MANUAL CONTROL SYSTEM I LAB EE 593 ELECTRICAL ENGINEERING DEPARTMENT JIS COLLEGE OF ENGINEERING (AN AUTONOMOUS INSTITUTE) KALYANI, NADIA CONTROL SYSTEM I LAB. MANUAL EE 593 EXPERIMENT
More informationHomework 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 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 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 informationEECE 460 : Control System Design
EECE 460 : Control System Design SISO Pole Placement Guy A. Dumont UBC EECE January 2011 Guy A. Dumont (UBC EECE) EECE 460: Pole Placement January 2011 1 / 29 Contents 1 Preview 2 Polynomial Pole Placement
More informationChapter 9 Robust Stability in SISO Systems 9. Introduction There are many reasons to use feedback control. As we have seen earlier, with the help of a
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Robust
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 informationThe 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 informationAutomatic 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 informationThe loop shaping paradigm. Lecture 7. Loop analysis of feedback systems (2) Essential specifications (2)
Lecture 7. Loop analysis of feedback systems (2). Loop shaping 2. Performance limitations The loop shaping paradigm. Estimate performance and robustness of the feedback system from the loop transfer L(jω)
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 informationNotes for ECE-320. Winter by R. Throne
Notes for ECE-3 Winter 4-5 by R. Throne Contents Table of Laplace Transforms 5 Laplace Transform Review 6. Poles and Zeros.................................... 6. Proper and Strictly Proper Transfer Functions...................
More informationControl for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e
Control for Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e Motion Systems m F Introduction Timedomain tuning Frequency domain & stability Filters Feedforward Servo-oriented
More informationDesign 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 informationAlireza Mousavi Brunel University
Alireza Mousavi Brunel University 1 » Control Process» Control Systems Design & Analysis 2 Open-Loop Control: Is normally a simple switch on and switch off process, for example a light in a room is switched
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 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 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 informationCDS 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 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 informationECE382/ME482 Spring 2005 Homework 6 Solution April 17, (s/2 + 1) s(2s + 1)[(s/8) 2 + (s/20) + 1]
ECE382/ME482 Spring 25 Homework 6 Solution April 17, 25 1 Solution to HW6 P8.17 We are given a system with open loop transfer function G(s) = 4(s/2 + 1) s(2s + 1)[(s/8) 2 + (s/2) + 1] (1) and unity negative
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 informationFeedback Control of Linear SISO systems. Process Dynamics and Control
Feedback Control of Linear SISO systems Process Dynamics and Control 1 Open-Loop Process The study of dynamics was limited to open-loop systems Observe process behavior as a result of specific input signals
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 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 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 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 informationQualitative 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 informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More 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 informationTable of Laplacetransform
Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014
Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can
More informationAndrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain)
1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Linear systems (frequency domain) 2 Motivations Consider an LTI system Thanks to the Lagrange s formula we can compute the motion of
More informationAN INTRODUCTION TO THE CONTROL THEORY
Open-Loop controller An Open-Loop (OL) controller is characterized by no direct connection between the output of the system and its input; therefore external disturbance, non-linear dynamics and parameter
More informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More 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 informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
More informationControl Systems II. ETH, MAVT, IDSC, Lecture 4 17/03/2017. G. Ducard
Control Systems II ETH, MAVT, IDSC, Lecture 4 17/03/2017 Lecture plan: Control Systems II, IDSC, 2017 SISO Control Design 24.02 Lecture 1 Recalls, Introductory case study 03.03 Lecture 2 Cascaded 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 informationCompensator Design to Improve Transient Performance Using Root Locus
1 Compensator Design to Improve Transient Performance Using Root Locus Prof. Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning
More informationDigital Control: Summary # 7
Digital Control: Summary # 7 Proportional, integral and derivative control where K i is controller parameter (gain). It defines the ratio of the control change to the control error. Note that e(k) 0 u(k)
More informationDr Ian R. Manchester
Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus 7 Root Locus 2 Assign
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 informationCDS 101/110: Lecture 10.3 Final Exam Review
CDS 11/11: Lecture 1.3 Final Exam Review December 2, 216 Schedule: (1) Posted on the web Monday, Dec. 5 by noon. (2) Due Friday, Dec. 9, at 5: pm. (3) Determines 3% of your grade Instructions on Front
More informationStability of Feedback Control Systems: Absolute and Relative
Stability of Feedback Control Systems: Absolute and Relative Dr. Kevin Craig Greenheck Chair in Engineering Design & Professor of Mechanical Engineering Marquette University Stability: Absolute and Relative
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Introduction to Describing Functions Hanz Richter Mechanical Engineering Department Cleveland State University Introduction Frequency domain methods
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time
More informationDefinition of Stability
Definition of Stability Transfer function of a linear time-invariant (LTI) system Fs () = b 2 1 0+ b1s+ b2s + + b m m m 1s - - + bms a0 + a1s+ a2s2 + + an-1sn- 1+ ansn Characteristic equation and poles
More informationToday (10/23/01) Today. Reading Assignment: 6.3. Gain/phase margin lead/lag compensator Ref. 6.4, 6.7, 6.10
Today Today (10/23/01) Gain/phase margin lead/lag compensator Ref. 6.4, 6.7, 6.10 Reading Assignment: 6.3 Last Time In the last lecture, we discussed control design through shaping of the loop gain GK:
More information