Rules to draw Polar Plot | Control Systems - Electrical Engineering (EE) PDF Download

Introduction

• The polar plot of a sinusoidal transfer function is the plot of the magnitude G(jω) versus the phase angle of G(jω) on the polar coordinates. The frequency in the polar plot is varied from zero to infinity. The polar plot is drawn on the polar sheet, which is the form of a graph, and the graph consists of concentric circles and radial lines.
• The concentric circles on the polar sheet graph represent the magnitude, and the radial lines represent phase angles. Each point on the graph displays information about the magnitude and the phase angle.
• The example of a polar graph is shown below:
  
• The positive angle in a polar graph is measured in the anti-clockwise direction, while the negative angle is measured clockwise direction. Both the angles are measured with respect to the reference point, i.e., 0 degree axis.
• The transfer function G(jω) in the rectangular form can be written as:
  G(jω) = GR(jω) + GI(jω)
  Where,
  GR(jω) is the real part of the transfer function G(jω)
  GI(jω) is the imaginary part of the transfer function G(jω)
• As discussed, the angular frequency in the polar plot varies from zero to infinity. We should not get confused with the Nyquist plot, which is the extension of the polar plot. The angular frequency (ω) varies from zero to infinity, while the Nyquist plot varies from a negative value of infinity to positive infinity.
• The primary advantage of the polar plot is that it depicts the frequency response characteristics of a system over the entire frequency range in a single plot. Since everything seems to be a single block, it fails to show the contributions of each factor of the open-loop transfer function.
• The integral factor of the polar plot is given by:
  G(s) = 1/s
  Where,
  S = jω
  s is the transfer function
  G(jω) = 1/jω
  It is the negative imaginary axis.
• The derivative factor of the polar plot is given by:
  G(jω) = jω
  It is the positive imaginary axis.
• Effect of addition of pole and zero
  • The addition of pole in the polar plot will shift its end by -90 degrees.
  • The addition of zero in the polar plot will shift its end by +90 degrees.

Type and order of the system

• The type of the system in the polar plot determines the quadrant at which the polar plot starts. The start of the polar plot of all poles is shown in the below diagram.

• The order of the system in the polar plot determines the quadrant at which the polar plot ends. The end of the polar plot of all poles is shown in the below diagram.
  

Polar plot of standard transfer function

The polar plots of various standard transfer functions.
Type 0

Order: 1
Let, G(s) = 1/(1 + sT)
Put, s =  jω
G(jω) = 1/(1 + jωT)
The above transfer function in the form of magnitude and angle can be represented as:
G(jω) = 1/[(1 + ω² T²)^1/2∠tan^-1ωT]

If we consider the angle part in the numerator, we need to insert a negative sign due to the transition from the denominator to the numerator or vice-versa.
G(jω) = 1∠-tan-1ωT /(1 + ω2 T2)1/2
Let's find the value of the above function at zero and infinity.
When, ω = 0
G(jω) = 1∠0/1
G(jω) = 1∠0
It is because tan-10 = 0
When, ω = infinity
G(jω) = 0∠-90
It is because tan-1 ∞ = 90 degrees
The polar plot at value 0 and infinity will appear as:

Order 2:
Since, it is a 2 order system; the function includes the highest derivative variable (s) with the power 2.
Let, G(s) = 1/(1 + sT₁) (1 + sT₂)
Put, s= jω
G(jω) = 1/(1 + jωT₁) (1 + jωT₂)
The above transfer function in the form of magnitude and angle can be represented as:
G(jω) = 1/[(1 + ω² T₁ ²)^1/2∠tan^-1ωT₁][(1 + ω² T₂ ²)^1/2∠tan_-1ωT₂]
If we consider the angle part in the numerator, we need to insert a negative sign due to the transition from the denominator to the numerator, as shown below:
G(jω) = 1∠(-tan^-1ωT₁ - tan^-1ωT₂) /(1 + ω² T₁ ²)^1/2(1 + ω² T₂ ²)^1/2
Let's find the value of the above function at zero and infinity.
When, ω = 0
G(jω) = 1∠(-0 - 0)/1x1
G(jω) = 1∠0
It is because tan-10 = 0
When, ω = infinity
G(jω) = 0∠(-90 - 90)/1
G(jω) = 0∠-180
It is because tan^-1 ∞ = 90 degrees
The polar plot at value 0 and infinity will appear as:

Type 1

Order: 1
Let, G(s) = 1/s
Put, s= jω
G(jω) = 1/jω
G(jω) = 1/(ω∠90)
If we consider the angle part in the numerator, we need to insert a negative sign due to the transition from the denominator to the numerator, as shown below:
G(jω) =1∠-90 /ω
Let's find the value of the above function at zero and infinity.
When, ω = 0
G(jω) =∞ ∠-90
It is because the function is directly divided by ω. 1/0 = infinity.
tan^-10 = 0
When, ω = infinity
G(jω) = 0∠-90
It is because the function is directly divided by ω. 1/infinity = zero.
tan^-1 ∞ = 90 degrees
The polar plot at value 0 and infinity will appear as:

Order: 2
Let, G(s) = 1/s(1 + sT)
Put, s= jω
G(jω) = 1/ jω(1 + jωT)
The above transfer function in the form of magnitude and angle can be represented as:
G(jω) = 1/ω∠90 [(1 + ω² T²)^1/2∠tan^-1ωT]
If we consider the angle part in the numerator, we need to insert a negative sign due to the transition from the denominator to the numerator or vice-versa.
G(jω) = 1∠(-90 - tan^-1ωT) / ω(1 + ω² T²)^1/2
Let's find the value of the above function at zero and infinity.
When, ω = 0
G(jω) =∞ ∠-90
It is because the function is directly divided by ω. 1/0 = infinity.
tan^-10 = 0
When, ω = infinity
G(jω) = 0∠(-90-90)
G(jω) = 0∠-180
It is because the function is directly divided by ω. 1/infinity = zero.
tan^-1 ∞ = 90 degrees
The polar plot at value 0 and infinity will appear as:

Order: 3
Let, G(s) = 1/s(1 + sT₁) (1 + sT₂)
Put, s = jω
G(jω) = 1/ jω (1 + jωT₁) (1 + jωT₂)
The above transfer function in the form of magnitude and angle can be represented as:
G(jω) = 1/ ω∠90 [(1 + ω² T₁ ²)^1/2∠tan^-1ωT₁][( 1 + ω² T₂ ²)^1/2∠tan^-1ωT₂]
If we consider the angle part in the numerator, we need to insert a negative sign due to the transition from the denominator to the numerator, as shown below:
G(jω) = 1∠(-90 -tan^-1ωT₁ - tan^-1ωT₂) /ω(1 + ω² T₁²)^1/2(1 + ω² T₂²)^1/2
Let's find the value of the above function at zero and infinity.
When, ω = 0
G(jω) =∞ ∠-90
It is because the function is directly divided by ω. 1/0 = infinity.
tan^-10 = 0
When, ω = infinity
G(jω) = 0∠(-90 -90 -90)
G(jω) = 0∠-270
It is because the function is directly divided by ω. 1/infinity = zero.
tan^-1 ∞ = 90 degrees
The polar plot at value 0 and infinity will appear as:

Type 2

Order: 4
Let, G(s) = 1/s²(1 + sT₁) (1 + sT₂)
Put, s= jω
G(jω) = 1/ j² ω² (1 + jωT₁) (1 + jωT₂)
The above transfer function in the form of magnitude and angle can be represented as:
G(jω) = 1/ ω²∠-180 [(1 + ω² T₁ ²)^1/2∠tan^-1ωT₁][( 1 + ω² T₂ ²)^1/2∠tan^-1ωT₂]

If we consider the angle part in the numerator, we need to insert a negative sign due to the transition from the denominator to the numerator, as shown below:
G(jω) = 1∠(-180 -tan^-1ωT₁ - tan^-1ωT₂) /(1 + ω² T₁²)^1/2(1 + ω² T₂²)^1/2
Let's find the value of the above function at zero and infinity.
When, ω = 0
G(jω) = ∞ ∠(-180 - 0 - 0)/1x1
G(jω) = ∞ ∠-180
It is because the function is directly divided by ω. 1/0 is equal to infinity.
tan-10 = 0
When, ω = infinity
G(jω) = 0∠(-180 - 90 - 90)/1
G(jω) = 0∠-360
It is because the function is directly divided by ω. 1/infinity is equal to zero.
tan^-1 ∞ = 90 degrees
The polar plot at value 0 and infinity will appear as:

Gain margin and phase margin of polar plot

• The gain margin is Kg. It is given by:
  Kg = 1/GB
• GB is the value at a point (B) on the magnitude circle that cuts the 180 degree axis. The gain margin is positive if the point B lies within unit circle. Otherwise, it is negative.
• The phase margin is given by:
  Y = 180 + theta
  Where,
• Theta is the phase angle of G(jω) at the gain crossover frequency. The phase angle is calculated when the magnitude curve line intersects with the external unity circle. The line drawn from that intersection point to the end of the graph determines the angle theta. It can be positive or negative.
• Both the phase margin and the gain margin can be better understood with the help of an example.

Example: The open loop transfer function of a unity feedback system is given by G(s) = 1/s(s + 1)(2s + 1). Sketch the polar plot and also determine the gain margin and the phase margin.

The transfer function is given by:
G(s) = 1/s(s + 1)(2s + 1)
The above function clearly depicts that the system is of type 1 and order 3. It is in the form:
G(s) = 1/s(1 + sT₁₎ (1 + sT₂)
Put, s = jω
G(s) = 1/ jω (jω + 1)(2jω + 1)
The above transfer function in the form of magnitude and angle can be represented as:
G(jω) = 1/ ω∠90 [(1 + ω²)^1/2∠tan^-1ω][( 1 + ω²4)^1/2∠tan^-1ω2]

If we consider the angle part in the numerator, we need to insert a negative sign due to the transition from the denominator to the numerator, as shown below:
G(jω) = 1∠(-90 -tan^-1ω - tan^-1 2ω) / ω(1 + ω²)^1/2(1 + 4ω²)^1/2
Now, let us separate the magnitude and angle terms from the above equation.
|G(jω)| = 1/ ω((1 + ω²) (1 + 4ω²))^1/2
|G(jω)| = 1/ ω(1 + 4ω² + ω²+ 4ω4)^1/2
|G(jω)| = 1/ ω(1 + 5ω² + 4ω4)1/2
Angle G(jω) = -90 -tan^-1ω - tan^-1 2ω
We know the value of the above function at zero and infinity.
When, ω = 0
G(jω) =∞ ∠-90
When, ω = infinity
G(jω) = 0∠-270
Let's find the magnitude and phase of G(jω) at different frequencies.

The polar plot is shown below:

• Here, every two lines have a gap of 15 degrees. We have specified both the positive and negative angle value at a point. It is based on the concept that the positive angles are measured anti-clockwise and the negative angles are measured clockwise. As we start from the 0 angle in the clockwise direction, we can notice the increase in the negative values of the angle. Similarly, as we begin from the 0 angle in the clockwise direction, we can notice the rise in the positive values of the angle.
• Let's calculate the gain margin and the phase margin.
• We can see in the polar plot that the magnitude circle cuts the 180 degree axis at point 0.7. Hence, it will be the value of GB.
  The gain margin is Kg = 1/GB
  Kg = 1/0.7
  Kg = 1.428
• The phase margin is given by:
  Y = 180 + theta
• We can clearly see the point marked with the red color. It is the intersection point of the magnitude curve with the unity circle. The line drawn from the intersection point (marked in red) determines the theta angle, which is equal to (-168) degrees.
  So, phase angle is equal to
  Y = 180 - 168
  Y = 12 degrees
• Thus, the gain margin is 1.428 and the phase angle is 12 degrees.