{"id":2085,"date":"2019-12-15T13:52:30","date_gmt":"2019-12-15T13:52:30","guid":{"rendered":"https:\/\/pressbooks.library.ryerson.ca\/controlsystems\/?post_type=chapter&p=2085"},"modified":"2021-01-14T16:18:48","modified_gmt":"2021-01-14T16:18:48","slug":"14-3-solved-examples-for-polar-plots","status":"publish","type":"chapter","link":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/chapter\/14-3-solved-examples-for-polar-plots\/","title":{"raw":"14.3 Solved Examples for Polar Plots","rendered":"14.3 Solved Examples for Polar Plots"},"content":{"raw":"14.3.1 Example<\/b>\r\n\r\n<\/b>Consider the following transfer function:\r\n

[latex]G(s) = \\frac{200}{s^3+11s^2+38s+4}[\/latex]<\/p>\r\nConsider its frequency response, [latex]G(j \\omega )[\/latex], at a specific frequency of [latex] \\omega =1[\/latex] rad\/sec. Show its rectangular and polar forms.\r\n

[latex]G(j \\omega )= \\frac{200}{(j \\omega )^3 +11(j \\omega )^2 +38(j \\omega )+4}[\/latex]<\/p>\r\n

[latex]= \\frac{200}{-j( \\omega )^3 -11( \\omega )^2 +38(j \\omega )+4}=[\/latex] [latex] \\frac{200}{(4-11 \\omega^2 )+j \\omega (38- \\omega^2 )}[\/latex]<\/p>\r\n

[latex]G(j1)= \\frac{200}{-7+j37} =[\/latex] [latex] \\frac{200(-7-j37)}{49+1369} =-0.9873-j5.2186[\/latex]<\/p>\r\nThe polar form of this function is:\r\n

[latex] \\mid G(j1) \\mid = \\sqrt{(-0.9873)^2+(-5.2186)^2}=5.3112[\/latex]<\/p>\r\n

[latex] \\angle G(j1) =-1.7578rad=-100.71 ^\\circ [\/latex]<\/p>\r\n

[latex]G(j \\omega )= \\mid G(j \\omega ) \\mid \\cdot e^{(j \\angle G( \\omega )}[\/latex]<\/p>\r\n

[latex]G(j1)=5.3112 \\cdot e^{-j1.7578}[\/latex]<\/p>\r\n

[latex]G(j1) = 5.3112 \\cdot e^{-j1.7578}=-0.9873-j5.2186[\/latex]<\/p>\r\n14.3.2 Example<\/b>\r\n\r\nConsider a simple first order system, with one real pole:\r\n

[latex]G(s) = \\frac{1}{10s+1}[\/latex]<\/p>\r\n

[latex]G(j \\omega )=[\/latex] [latex] \\frac{1}{10j \\omega +1} =[\/latex] [latex] \\frac{1}{ \\sqrt{(10 \\omega )^2+1} } \\angle -tan^{-1}(10 \\omega )[\/latex]<\/p>\r\n

[latex]M(j \\omega )= \\frac{1}{ \\sqrt{(10 \\omega )^2 +1} }[\/latex]<\/p>\r\n

[latex] \\Phi (j \\omega )= -tan^{-1}(10 \\omega )[\/latex]<\/p>\r\nNow consider the rectangular representation of the same frequency response function [latex]G(j \\omega ) [\/latex]:\r\n

[latex]G(j \\omega )= \\frac{1}{1+10j \\omega )} = \\frac{1(1-10j \\omega )}{ (1+10j \\omega)(1-10j \\omega)} = \\frac{1}{1+100 \\omega^2 }-j \\frac{10 \\omega }{1+100 \\omega^2 }[\/latex]<\/p>\r\n

[latex]Re( \\omega ) = \\frac{1}{1+100 \\omega^2 }[\/latex]<\/p>\r\n

[latex]Im( \\omega )=- \\frac{10 \\omega }{1+100 \\omega^2 }[\/latex]<\/p>\r\nThe standard frequency response plot (Bode Plot) with magnitude in decibels and phase in degrees is shown below. For the Polar Plot, crossovers with the Imaginary and Real axes can be calculated analytically by setting first the Real, then the Imaginary part to zero, and solving for frequency. In this example:\r\n

[latex]Re( \\omega )=0 \\implies [\/latex] [latex] \\omega = \\infty [\/latex]<\/p>\r\n

[latex]Im ( \\infty )=0[\/latex]<\/span><\/p>\r\n

[latex]Im ( \\infty )=0 \\implies [\/latex]\u00a0 [latex] \\omega = 0,[\/latex] [latex] \\omega = \\infty [\/latex]<\/p>\r\n

[latex] Re(0) =1[\/latex]<\/p>\r\n

[latex] Re(\\infty)=0[\/latex]<\/p>\r\n \r\n

\"\"<\/p>\r\n

\"\"<\/p>\r\n \r\n\r\nThis indicates that the polar plot starts at (1, j0) location for [latex] \\omega\u00a0 =\u00a0 0[\/latex] (DC condition), and ends at (0, j0) for [latex] \\omega =\u00a0 \\infty[\/latex] . The sense of increasing frequency [latex] \\omega [\/latex] should always be shown on the polar plot. The polar plot of the system [latex]G(s)[\/latex] is shown.\r\n\r\nTo do plot polar plots in MATLAB, use subroutine Nyquist - see below<\/strong>. The second plot (on the following page) shows a so-called Nyquist contour<\/strong>, which will be discussed in detail later. The Nyquist contour consists of the polar plot for positive frequencies,[latex]0 < \\omega < + \\infty [\/latex], and its mirror image for negative frequencies, [latex]- \\infty < \\omega < 0[\/latex].\r\n\r\n \r\n

\"\"<\/p>\r\n

\"\"<\/p>\r\n \r\n\r\n14.3.3 Example<\/b>\r\n\r\nA process transfer function is described as follows: [latex]G(s) = \\frac{10}{s^3+4s^2+6s+8}.[\/latex] Frequency plots of [latex]G(s)[\/latex] are shown. Sketch a polar plot for [latex]G(s)[\/latex].\r\n\r\nSolution:\u00a0\u00a0 <\/strong>It is helpful to construct a table with the important coordinates:\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
Frequency<\/td>\r\nPhase<\/td>\r\nMagnitude<\/td>\r\n<\/tr>\r\n
[latex] \\omega = 0 [\/latex] rad\/s<\/td>\r\n[latex] \\phi = 0 ^ \\circ[\/latex]<\/td>\r\n[latex]1.25[\/latex]<\/td>\r\n<\/tr>\r\n
[latex] \\omega = 1.4[\/latex] rad\/s<\/td>\r\n[latex] \\phi = -90 ^ \\circ[\/latex]<\/td>\r\n[latex]1.77[\/latex]<\/td>\r\n<\/tr>\r\n
[latex] \\omega = 2.4[\/latex] rad\/s<\/td>\r\n[latex] \\phi = -180 ^ \\circ[\/latex]<\/td>\r\n[latex]0.625[\/latex]<\/td>\r\n<\/tr>\r\n
[latex] \\omega = + \\infty [\/latex] rad\/s<\/td>\r\n[latex] \\phi = -270 ^ \\circ[\/latex]<\/td>\r\n[latex]0[\/latex]<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nThe resulting polar plot can be also plotted using MATLAB subroutine Nyquist<\/strong>.\r\n\r\n \r\n

\"\"<\/p>\r\n

\"\"<\/p>\r\n \r\n\r\n14.3.4 Example<\/b>\r\n\r\nConsider a unity feedback control system under Proportional Control. The process transfer function is described as follows:\r\n

[latex]G(s) = \\frac{5}{s(s+1)(s+5)}[\/latex]<\/p>\r\nFrequency plots of [latex]G(s)[\/latex] are shown. It is helpful to construct a table with the important coordinates read off the plot. Note that this is a Type I system, with an integrator, and therefore its polar plot will begin with an infinite gain at the DC level.\r\n\r\n\r\n\r\n\r\n\r\n\r\n
Frequency<\/td>\r\nPhase<\/td>\r\nMagnitude in dB<\/td>\r\nMagnitude in Volt\/Volt<\/td>\r\n<\/tr>\r\n
[latex] \\omega = 0[\/latex]rad\/s<\/td>\r\n[latex] \\phi = -90 ^ \\circ[\/latex]<\/td>\r\n[latex] + \\infty [\/latex]<\/td>\r\n[latex] + \\infty [\/latex]<\/td>\r\n<\/tr>\r\n
[latex] \\omega = 224[\/latex]rad\/s<\/td>\r\n[latex] \\phi = -180 ^ \\circ[\/latex]<\/td>\r\n[latex]-15.56[\/latex] dB<\/td>\r\n[latex] 0.1667 [\/latex] Volt\/Volt<\/td>\r\n<\/tr>\r\n
[latex] \\omega = + \\infty [\/latex]rad\/s<\/td>\r\n[latex] \\phi = -270 ^ \\circ[\/latex]<\/td>\r\n[latex] - \\infty [\/latex]<\/td>\r\n[latex]0[\/latex]<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n \r\n

\"\"<\/p>\r\n \r\n\r\nThe resulting polar plot is shown.\r\n\r\n \r\n

\"\"<\/p>","rendered":"

14.3.1 Example<\/b><\/p>\n

<\/b>Consider the following transfer function:<\/p>\n

[latex]G(s) = \\frac{200}{s^3+11s^2+38s+4}[\/latex]<\/p>\n

Consider its frequency response, [latex]G(j \\omega )[\/latex], at a specific frequency of [latex]\\omega =1[\/latex] rad\/sec. Show its rectangular and polar forms.<\/p>\n

[latex]G(j \\omega )= \\frac{200}{(j \\omega )^3 +11(j \\omega )^2 +38(j \\omega )+4}[\/latex]<\/p>\n

[latex]= \\frac{200}{-j( \\omega )^3 -11( \\omega )^2 +38(j \\omega )+4}=[\/latex] [latex]\\frac{200}{(4-11 \\omega^2 )+j \\omega (38- \\omega^2 )}[\/latex]<\/p>\n

[latex]G(j1)= \\frac{200}{-7+j37} =[\/latex] [latex]\\frac{200(-7-j37)}{49+1369} =-0.9873-j5.2186[\/latex]<\/p>\n

The polar form of this function is:<\/p>\n

[latex]\\mid G(j1) \\mid = \\sqrt{(-0.9873)^2+(-5.2186)^2}=5.3112[\/latex]<\/p>\n

[latex]\\angle G(j1) =-1.7578rad=-100.71 ^\\circ[\/latex]<\/p>\n

[latex]G(j \\omega )= \\mid G(j \\omega ) \\mid \\cdot e^{(j \\angle G( \\omega )}[\/latex]<\/p>\n

[latex]G(j1)=5.3112 \\cdot e^{-j1.7578}[\/latex]<\/p>\n

[latex]G(j1) = 5.3112 \\cdot e^{-j1.7578}=-0.9873-j5.2186[\/latex]<\/p>\n

14.3.2 Example<\/b><\/p>\n

Consider a simple first order system, with one real pole:<\/p>\n

[latex]G(s) = \\frac{1}{10s+1}[\/latex]<\/p>\n

[latex]G(j \\omega )=[\/latex] [latex]\\frac{1}{10j \\omega +1} =[\/latex] [latex]\\frac{1}{ \\sqrt{(10 \\omega )^2+1} } \\angle -tan^{-1}(10 \\omega )[\/latex]<\/p>\n

[latex]M(j \\omega )= \\frac{1}{ \\sqrt{(10 \\omega )^2 +1} }[\/latex]<\/p>\n

[latex]\\Phi (j \\omega )= -tan^{-1}(10 \\omega )[\/latex]<\/p>\n

Now consider the rectangular representation of the same frequency response function [latex]G(j \\omega )[\/latex]:<\/p>\n

[latex]G(j \\omega )= \\frac{1}{1+10j \\omega )} = \\frac{1(1-10j \\omega )}{ (1+10j \\omega)(1-10j \\omega)} = \\frac{1}{1+100 \\omega^2 }-j \\frac{10 \\omega }{1+100 \\omega^2 }[\/latex]<\/p>\n

[latex]Re( \\omega ) = \\frac{1}{1+100 \\omega^2 }[\/latex]<\/p>\n

[latex]Im( \\omega )=- \\frac{10 \\omega }{1+100 \\omega^2 }[\/latex]<\/p>\n

The standard frequency response plot (Bode Plot) with magnitude in decibels and phase in degrees is shown below. For the Polar Plot, crossovers with the Imaginary and Real axes can be calculated analytically by setting first the Real, then the Imaginary part to zero, and solving for frequency. In this example:<\/p>\n

[latex]Re( \\omega )=0 \\implies[\/latex] [latex]\\omega = \\infty[\/latex]<\/p>\n

[latex]Im ( \\infty )=0[\/latex]<\/span><\/p>\n

[latex]Im ( \\infty )=0 \\implies[\/latex]\u00a0 [latex]\\omega = 0,[\/latex] [latex]\\omega = \\infty[\/latex]<\/p>\n

[latex]Re(0) =1[\/latex]<\/p>\n

[latex]Re(\\infty)=0[\/latex]<\/p>\n

 <\/p>\n

\"\"<\/p>\n

\"\"<\/p>\n

 <\/p>\n

This indicates that the polar plot starts at (1, j0) location for [latex]\\omega\u00a0 =\u00a0 0[\/latex] (DC condition), and ends at (0, j0) for [latex]\\omega =\u00a0 \\infty[\/latex] . The sense of increasing frequency [latex]\\omega[\/latex] should always be shown on the polar plot. The polar plot of the system [latex]G(s)[\/latex] is shown.<\/p>\n

To do plot polar plots in MATLAB, use subroutine Nyquist – see below<\/strong>. The second plot (on the following page) shows a so-called Nyquist contour<\/strong>, which will be discussed in detail later. The Nyquist contour consists of the polar plot for positive frequencies,[latex]0 < \\omega < + \\infty[\/latex], and its mirror image for negative frequencies, [latex]- \\infty < \\omega < 0[\/latex].\n\n \n\n\n

\"\"<\/p>\n

\"\"<\/p>\n

 <\/p>\n

14.3.3 Example<\/b><\/p>\n

A process transfer function is described as follows: [latex]G(s) = \\frac{10}{s^3+4s^2+6s+8}.[\/latex] Frequency plots of [latex]G(s)[\/latex] are shown. Sketch a polar plot for [latex]G(s)[\/latex].<\/p>\n

Solution:\u00a0\u00a0 <\/strong>It is helpful to construct a table with the important coordinates:<\/p>\n\n\n\n\n\n\n\n
Frequency<\/td>\nPhase<\/td>\nMagnitude<\/td>\n<\/tr>\n
[latex]\\omega = 0[\/latex] rad\/s<\/td>\n[latex]\\phi = 0 ^ \\circ[\/latex]<\/td>\n[latex]1.25[\/latex]<\/td>\n<\/tr>\n
[latex]\\omega = 1.4[\/latex] rad\/s<\/td>\n[latex]\\phi = -90 ^ \\circ[\/latex]<\/td>\n[latex]1.77[\/latex]<\/td>\n<\/tr>\n
[latex]\\omega = 2.4[\/latex] rad\/s<\/td>\n[latex]\\phi = -180 ^ \\circ[\/latex]<\/td>\n[latex]0.625[\/latex]<\/td>\n<\/tr>\n
[latex]\\omega = + \\infty[\/latex] rad\/s<\/td>\n[latex]\\phi = -270 ^ \\circ[\/latex]<\/td>\n[latex]0[\/latex]<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

The resulting polar plot can be also plotted using MATLAB subroutine Nyquist<\/strong>.<\/p>\n

 <\/p>\n

\"\"<\/p>\n

\"\"<\/p>\n

 <\/p>\n

14.3.4 Example<\/b><\/p>\n

Consider a unity feedback control system under Proportional Control. The process transfer function is described as follows:<\/p>\n

[latex]G(s) = \\frac{5}{s(s+1)(s+5)}[\/latex]<\/p>\n

Frequency plots of [latex]G(s)[\/latex] are shown. It is helpful to construct a table with the important coordinates read off the plot. Note that this is a Type I system, with an integrator, and therefore its polar plot will begin with an infinite gain at the DC level.<\/p>\n\n\n\n\n\n\n
Frequency<\/td>\nPhase<\/td>\nMagnitude in dB<\/td>\nMagnitude in Volt\/Volt<\/td>\n<\/tr>\n
[latex]\\omega = 0[\/latex]rad\/s<\/td>\n[latex]\\phi = -90 ^ \\circ[\/latex]<\/td>\n[latex]+ \\infty[\/latex]<\/td>\n[latex]+ \\infty[\/latex]<\/td>\n<\/tr>\n
[latex]\\omega = 224[\/latex]rad\/s<\/td>\n[latex]\\phi = -180 ^ \\circ[\/latex]<\/td>\n[latex]-15.56[\/latex] dB<\/td>\n[latex]0.1667[\/latex] Volt\/Volt<\/td>\n<\/tr>\n
[latex]\\omega = + \\infty[\/latex]rad\/s<\/td>\n[latex]\\phi = -270 ^ \\circ[\/latex]<\/td>\n[latex]- \\infty[\/latex]<\/td>\n[latex]0[\/latex]<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

 <\/p>\n

\"\"<\/p>\n

 <\/p>\n

The resulting polar plot is shown.<\/p>\n

 <\/p>\n

\"\"<\/p>\n","protected":false},"author":162,"menu_order":3,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-2085","chapter","type-chapter","status-publish","hentry"],"part":2068,"_links":{"self":[{"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/pressbooks\/v2\/chapters\/2085","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/wp\/v2\/users\/162"}],"version-history":[{"count":23,"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/pressbooks\/v2\/chapters\/2085\/revisions"}],"predecessor-version":[{"id":2733,"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/pressbooks\/v2\/chapters\/2085\/revisions\/2733"}],"part":[{"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/pressbooks\/v2\/parts\/2068"}],"metadata":[{"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/pressbooks\/v2\/chapters\/2085\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/wp\/v2\/media?parent=2085"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/pressbooks\/v2\/chapter-type?post=2085"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/wp\/v2\/contributor?post=2085"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/pressbooks.library.torontomu.ca\/controlsystems\/wp-json\/wp\/v2\/license?post=2085"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}