Difference between revisions of "Frequency Domain Design"
(→MATLAB code) 

(13 intermediate revisions by the same user not shown)  
Line 5:  Line 5:  
{{chaptertable left}}  {{chaptertable left}}  
== Textbook Contents ==  == Textbook Contents ==  
−  {{am05pdfam08loopsyn  +  {{am05pdfam08loopsyn28Sep12Frequency Domain Design}} 
* 1. A Basic Feedback Loop  * 1. A Basic Feedback Loop  
* 2. Performance Specifications  * 2. Performance Specifications  
Line 16:  Line 16:  
== Lecture Materials ==  == Lecture Materials ==  
−  * [[Lecture:  +  * Lectures: [[Lecture: Loop ShapingLoop Shaping]], [[Lecture: Limits of PerformanceLimits of Performance]] 
* [[#ExercisesAdditional Exercises]]  * [[#ExercisesAdditional Exercises]]  
Line 22:  Line 22:  
* [[#Frequently Asked QuestionsFrequently Asked Questions]]  * [[#Frequently Asked QuestionsFrequently Asked Questions]]  
* [[#ErrataErrata]]  * [[#ErrataErrata]]  
−  * Wikipedia entries:  +  * Wikipedia entries: [http://en.wikipedia.org/wiki/Leadlag_compensator leadlag compensator] 
* [[#Additional InformationAdditional Information]]  * [[#Additional InformationAdditional Information]]  
{{chaptertable end}}  {{chaptertable end}}  
== Chapter Summary ==  == Chapter Summary ==  
+  This chapter describes the use of the ''loop shaping'' to design a control law to achieve a desired level of performance:  
+  <ol>  
+  <li><p>Control design in the frequency domain consists of designing a ''feedback controller'' <math>C(s)</math> and a ''feedforward controller'' <math>F(s)</math> as illustrated in the basic feedback diagram below:  
+  <center>[[Image:loopsyn_fbksys.png]]</center>  
+  The external signals are the reference signal <math>r</math>, the load disturbance <math>d</math> and the measurement noise <math>n</math>. The process output is <math>\eta</math>, and the control signal is <math>u</math>. This control structure is also known as a ''two degree of freedom controller''.  
+  </p></li>  
−  == Exercises ==  +  <li><p> 
+  The primary transfer functions that define the input/output characteristics of the system are called the ''Gang of Six'':  
+  <center><amsmath>  
+  \begin{aligned}  
+  & TF=\frac{PCF}{1+PC}, \quad & \quad  
+  & T=\frac{PC}{1+PC}, \quad & \quad  
+  & PS=\frac{P}{1+PC}, \\[1ex]  
+  & CFS=\frac{CF}{1+PC}, \quad & \quad  
+  & CS=\frac{C}{1+PC}, \quad & \quad  
+  & S=\frac{1}{1+PC}.  
+  \end{aligned}  
+  </amsmath></center>  
+  The transfer functions in the first column give the response of the process output and control signal to the reference signal. The second column gives the response of the control variable to the load disturbance and the noise, and the final column gives the response of the process output to those two inputs. When <math>F(s) = 1</math>, the system is said to have pure error feedback and the relevant input/output transfer functions are given by the ''Gang of Four'', given by the transfer functions in the right two columns.  
+  </p></li>  
+  
+  <li><p> Feedforward compensation is used to provided improved response to reference signals and measured disturbances (2 DOF system). The figure below illustrates one possible structure for designing a feedforward compensator:  
+  <center>[[Image:loopsyn_twodof.png]]</center>  
+  Three feedforward elements are present: <math>F_m(s)</math> sets the desired output value, <math>F_u(s)</math> generates the feedforward command <math>u_{ff}</math> and <math>F_d(s)</math> attempts to cancel disturbances. The transfer function from reference input to process output is given by  
+  <center><amsmath>  
+  G_{yr}(s)=\frac{P(CF_m+F_u)}{1+PC}=F_m+\frac{PF_uF_m}{1+PC},  
+  </amsmath></center>  
+  The first term represents the desired transfer function and the second term should be made small. Perfect feedforward compensation is obtained by choosing <math>F_u =F_m/P</math>.  
+  </p></li>  
+  
+  <li><p>  
+  The performance of a system can be given in terms of the characteristics of the frequency response between an input and output. A ''resonant peak'' is a maximum of the gain, and the peak frequency is the corresponding  
+  frequency. The ''gain crossover frequency'' is the frequency where the open loop gain is equal one. The ''bandwidth'' is defined as the frequency range where the closed loop gain is <amsmath>1/\sqrt{2}</amsmath> of the lowfrequency gain  
+  (lowpass), midfrequency gain (bandpass) or highfrequency gain (highpass).  
+  </p></li>  
+  
+  <li><p> The ''sensitivity function'' <math>S = 1/(1 + PC)</math> describes how disturbances are attenuated by closing the feedback loop. Disturbances with frequencies such that <math>S(i\omega)<1</math> are attenuated, but disturbances with frequencies such that <math>S(i\omega)>1</math> are amplified by feedback. The maximum sensitivity <math>M_s</math>, which occurs at the frequency <math>\omega_{ms}</math>, is a measure of the largest amplification of the disturbances. The ''complementary sensitivity function'' <math>T = PC / (1 + PC)</math> describes how well the controller tracks a references signal. The ''maximum complementary sensitivity'', <math>M_t</math>, which occurs at the frequency <math>\omega_{mt}</math>, is the peak value of the magnitude of the complementary sensitivity function. It provides the maximum amplification from the reference signal to the output signal.  
+  </p></li>  
+  
+  <li><p> There are approximate relations between specifications in the time and frequency domains. In the frequency  
+  domain the response time can be characterized by the closed loop ''bandwidth'' <math>\omega_b</math>, the ''gain crossover frequency'' <math>\omega_{gc}</math> and the sensitivity frequency <math>\omega_{ms}</math>. The product of bandwidth and rise time is approximately constant, <amsmath>T_r\omega_b\approx 2</amsmath>, so decreasing the rise time corresponds to increasing the closed loop bandwidth. The overshoot of the step response <math>M_p</math> is related to the resonant peak <math>M_r</math> of the frequency response in the sense that a larger peak normally implies a larger overshoot.  
+  </p></li>  
+  
+  <li><p> Using frequency domain specifications, controllers can often be designed by shaping the loop transfer function, <math>L = PC</math>. The gain curve and sensitivity functions for a typical loop transfer function are shown below:  
+  <center>[[Image:loopsyn_typicalbode.png]]</center>  
+  The plot on the left shows the gain curve and the plots on the right show the sensitivity function and complementary sensitivity function. The gain crossover frequency <math>\omega_{gc}</math> and the slope <math>n_{gc}</math> of the gain curve at crossover are important parameters that determine the robustness of closed loop systems. At low frequency, a large magnitude for <math>L</math> provides good load disturbance rejection and reference tracking, while at high frequency a small loop gain is used to avoid amplifying measurement noise.  
+  </p></li>  
+  
+  <li><p> A common control element to shape the loop transfer function is to use a pole/zero pair:  
+  <center><amsmath>  
+  C(s) = k \frac{s + a}{s + b}.  
+  </amsmath></center>  
+  The compensator is called a ''lead compensator'' if <math>a<b</math> the controller provides phase lead between frequencies <math>\omega = a</math> and <math>\omega = b</math>. The compensator is called a ''lag compensator''  
+  if <math>a>b</math> and this form is often used to provide low frequency gain. Lead and lag compensation are closely related to PD and PI controllers.  
+  </p></li>  
+  
+  <li><p> Feedback control systems have a number of fundamental limits, usually exacerbated by the presence of right half plane poles and zeros. For systems with right half plane poles or zeros, we can decompose the process dynamics into a minimum phase transfer function (no right half plane poles or zeros) and an all pass transfer function (gain = 1):  
+  <center><amsmath>  
+  P(s)=P_{mp}(s)P_{ap}(s),  
+  </amsmath></center>  
+  The ''gain crossover inequality''  
+  <center><amsmath>  
+  \arg P_{ap}(i\omega_\text{gc})  
+  \leq \pi\varphi_m+n_\text{gc}\frac{\pi}{2} =: \varphi_l.  
+  </amsmath></center>  
+  provides a relationship between the phase margin <amsmath>\varphi_m</amsmath>, the slope of the gain curve <math>n_{gc}</math>. For processes with near pole/zero cancellations in the right half plane, the gain crossover inequality limits the maximum amount of achievable phase margin.  
+  </p></li>  
+  
+  <li><p> Another fundamental limit is given by ''Bode's integral formula'', which states that for systems with a loop transfer function that goes to zero faster than <math>1/s</math> as <amsmath>s\rightarrow\infty</amsmath>, the sensitivity function must satisfy  
+  <center><amsmath>  
+  \int_0^\infty\log{S(i\omega)}\, d\omega  
+  = \int_0^\infty \log{\frac{1}{1+L(i\omega)}}\, d\omega  
+  = \pi\sum p_k,  
+  </amsmath></center>  
+  where <math>p_k</math> are the poles in the right halfplane. This conservation law shows that to get lower sensitivity in one frequency range, we must get higher sensitivity in some other region. An analogous formula exists for the complementary sensitivity function in the presence of right half plane zeros.  
+  </p></li>  
+  </ol>  
+  
+  {{chaptertable begin}}  
+  {{chaptertable left}}  
+  
+  <!  
+  == Additional Exercises ==  
+  The following exercises cover some of the topics introduced in this chapter. Exercises marked with a * appear in the printed text.  
<ncl>Frequency Domain Design Exercises</ncl>  <ncl>Frequency Domain Design Exercises</ncl>  
+  >  
== Frequently Asked Questions ==  == Frequently Asked Questions ==  
<ncl>Frequency Domain Design FAQ</ncl>  <ncl>Frequency Domain Design FAQ</ncl>  
== Errata ==  == Errata ==  
−  <ncl>Frequency Domain Design errata</ncl>  +  <ncl>Frequency Domain Design errata v2.11b</ncl> 
+  * [[:Category:Frequency Domain Design errataFull list of errata starting from first printing]]  
* {{submitbug}}  * {{submitbug}}  
+  {{chaptertable right}}  
+  
+  == MATLAB code ==  
+  The following MATLAB scripts are available for producing figures that appear in this chapter.  
+  * Figure 11.4: {{matlabfileloopsynsteering_feedforward.m}}  
+  * Figure 11.5: {{matlabfileloopsynrefresp.m}}  
+  * Figure 11.7: {{matlabfileloopsyndistresp.m}}  
+  * Figure 11.10: {{matlabfileloopsynafmtap_ltfbode.m}}  
+  * Figure 11.12, 11.1820: {{matlabfileloopsynpvtol_nested.m}}, {{matlabdirpvtol}}  
+  ** See also: [[Vectored thrust aircraft]] example  
+  See the [[softwaresoftware page]] for more information on how to run these scripts.  
+  
+  <!  
== Additional Information ==  == Additional Information ==  
+  >  
+  {{chaptertable end}} 
Latest revision as of 21:27, 23 November 2012
Prev: PID Control  Chapter 11  Frequency Domain Design  Next: Robust Performance 
In this chapter we continue to explore the use of frequency domain techniques for design of feedback systems. We begin with a more thorough description of the performance specifications for controls systems, and then introduce the concept of "loop shaping" as a mechanism for designing controllers in the frequency domain. We also introduce some fundamental limitations to performance for systems with right half plane poles and zeros.
Textbook ContentsFrequency Domain Design (pdf, 28Sep12)

Lecture MaterialsSupplemental Information

Chapter Summary
This chapter describes the use of the loop shaping to design a control law to achieve a desired level of performance:
Control design in the frequency domain consists of designing a feedback controller and a feedforward controller as illustrated in the basic feedback diagram below:
The external signals are the reference signal , the load disturbance and the measurement noise . The process output is , and the control signal is . This control structure is also known as a two degree of freedom controller.
The primary transfer functions that define the input/output characteristics of the system are called the Gang of Six:
The transfer functions in the first column give the response of the process output and control signal to the reference signal. The second column gives the response of the control variable to the load disturbance and the noise, and the final column gives the response of the process output to those two inputs. When , the system is said to have pure error feedback and the relevant input/output transfer functions are given by the Gang of Four, given by the transfer functions in the right two columns.
Feedforward compensation is used to provided improved response to reference signals and measured disturbances (2 DOF system). The figure below illustrates one possible structure for designing a feedforward compensator:
Three feedforward elements are present: sets the desired output value, generates the feedforward command and attempts to cancel disturbances. The transfer function from reference input to process output is given by
The first term represents the desired transfer function and the second term should be made small. Perfect feedforward compensation is obtained by choosing .
The performance of a system can be given in terms of the characteristics of the frequency response between an input and output. A resonant peak is a maximum of the gain, and the peak frequency is the corresponding frequency. The gain crossover frequency is the frequency where the open loop gain is equal one. The bandwidth is defined as the frequency range where the closed loop gain is of the lowfrequency gain (lowpass), midfrequency gain (bandpass) or highfrequency gain (highpass).
The sensitivity function describes how disturbances are attenuated by closing the feedback loop. Disturbances with frequencies such that are attenuated, but disturbances with frequencies such that are amplified by feedback. The maximum sensitivity , which occurs at the frequency , is a measure of the largest amplification of the disturbances. The complementary sensitivity function describes how well the controller tracks a references signal. The maximum complementary sensitivity, , which occurs at the frequency , is the peak value of the magnitude of the complementary sensitivity function. It provides the maximum amplification from the reference signal to the output signal.
There are approximate relations between specifications in the time and frequency domains. In the frequency domain the response time can be characterized by the closed loop bandwidth , the gain crossover frequency and the sensitivity frequency . The product of bandwidth and rise time is approximately constant, , so decreasing the rise time corresponds to increasing the closed loop bandwidth. The overshoot of the step response is related to the resonant peak of the frequency response in the sense that a larger peak normally implies a larger overshoot.
Using frequency domain specifications, controllers can often be designed by shaping the loop transfer function, . The gain curve and sensitivity functions for a typical loop transfer function are shown below:
The plot on the left shows the gain curve and the plots on the right show the sensitivity function and complementary sensitivity function. The gain crossover frequency and the slope of the gain curve at crossover are important parameters that determine the robustness of closed loop systems. At low frequency, a large magnitude for provides good load disturbance rejection and reference tracking, while at high frequency a small loop gain is used to avoid amplifying measurement noise.
A common control element to shape the loop transfer function is to use a pole/zero pair:
The compensator is called a lead compensator if the controller provides phase lead between frequencies and . The compensator is called a lag compensator if and this form is often used to provide low frequency gain. Lead and lag compensation are closely related to PD and PI controllers.
Feedback control systems have a number of fundamental limits, usually exacerbated by the presence of right half plane poles and zeros. For systems with right half plane poles or zeros, we can decompose the process dynamics into a minimum phase transfer function (no right half plane poles or zeros) and an all pass transfer function (gain = 1):
The gain crossover inequality
provides a relationship between the phase margin , the slope of the gain curve . For processes with near pole/zero cancellations in the right half plane, the gain crossover inequality limits the maximum amount of achievable phase margin.
Another fundamental limit is given by Bode's integral formula, which states that for systems with a loop transfer function that goes to zero faster than as , the sensitivity function must satisfy
where are the poles in the right halfplane. This conservation law shows that to get lower sensitivity in one frequency range, we must get higher sensitivity in some other region. An analogous formula exists for the complementary sensitivity function in the presence of right half plane zeros.
Frequently Asked Questions
Errata

MATLAB codeThe following MATLAB scripts are available for producing figures that appear in this chapter.
See the software page for more information on how to run these scripts. 