Nyquist Plots
Nyquist Plots
Introduction
• Nyquist plots are the continuation of polar
  plots for finding the stability of the closed loop
  control systems by varying ω from −∞ to ∞.
• Nyquist plots are used to draw the complete
  frequency response of the open loop transfer
  function.
• The Nyquist stability criterion determines the
  stability of a closed-loop system from its open-
  loop frequency response and open-loop poles.
The Nyquist Criterion can be expressed as,
            Z=P+N
  where Z = number of zeros of 1 + G(s)H(s) on the
             right-half s-plane
        N = net encirclements around the point (-1+j0).
             (clockwise encirclements are taken as
              positive and anticlockwise encirclements
              are negative)
         P = number of poles of G(s)H(s) in the
             right-half of s-plane
The stability of linear control systems using the Nyquist
stability criterion, three possibilities can occur:
1. There is no encirclement of the (–1+j0) point. This
    implies that the system is stable if there are no poles
    of G(s)H(s) in the right-half of s- plane; otherwise, the
    system is unstable.
2. There are one or more counterclockwise
    encirclements of the (–1+j0) point. In this case the
    system is stable if the number of counterclockwise
    encirclements is the same as the number of poles of
    G(s)H(s) in the right-half of s- plane; otherwise, the
    system is unstable.
3. There are one or more clockwise encirclements of the
    (–1+j0) point. In this case the system is unstable.
Consider the closed-loop system shown in Fig.
       put s = jω
                            10
       G(jω)H(jω) =
                       (jω+1)(jω+2)
                                       10
              =
                    ω2+1∠tan     −1ω    { ω2+4∠tan   −1(ω2)}
                        10                                      ω
              =                        − ∠tan −1ω − ∠tan −1    ( )
                    ω2+1 { ω2+4}                                2
                 10
       M=                        ;
             ω2+1 { ω2+4}
                                             ω
       Ø = − ∠tan −1ω − ∠tan −1             ( )
                                             2
lim 𝑀∠Ø = 5∠0      (point a’)
ω⟶0
lim 𝑀∠Ø = 0∠-180   (point b’)
ω⟶∞
Section II : To find the image of path ‘bcd’
      put s = lim 𝑅𝑒𝑗θ in G(s)H(s)
             𝑅⟶∞
     Here , θ changes from +90 ⟶ 0 ⟶ -90
                         𝑗θ           10
     Then, lim 𝐺𝐻( 𝑅𝑒 ) = lim      θ            𝑗θ
             𝑅⟶∞                 𝑅⟶∞ (𝑅𝑒𝑗 +1)(𝑅𝑒 +2)
                                          10
                             =    lim   𝑗θ   𝑗θ
                                 𝑅⟶∞ (𝑅𝑒 )(𝑅𝑒 )
                                       10
                             =    lim 2 𝑗2θ
                                 𝑅⟶∞ (𝑅 𝑒 )
                             = 0∠-2θ
                             = 0∠-180⟶ 0⟶ 180
                                   ⤒b’ ⤒c’ ⤒d’
Hence, the infinite semicircle ‘bcd’ on the s-
plane is mapped to the origin of the G(s)H(s)-
plane.
Section III: To find the image of path ‘da’
Path d’a’ is the mirror image of the path a’b’
with respect to real axis.
The complete Nyquist plot is shown below
Since, (-1+j0) point is not encircled by the plot
a’b’c’d’a’ in the GH(s) plane , N= 0
Step 3: The Nyquist stability criterion is Z = P + N
Hence, Z = 0 + 0 = 0
      ⟹No roots of the system lie to the right-
half of s-plane.
Hence the closed-loop control system is
stable.
Matlab
                  Problem
Sketch the Nyquist diagram for the system
shown in the following figure, and then
determine the system stability using the Nyquist
criterion
Matlab
                     Problem
Using Nyquist stability criterion, Investigate the
stability of a closed-loop system whose open-loop
transfer function is given by,
                                  (𝑠+2)
                    G(s)H(s) =
                               (𝑠+1)(𝑠−1)
Solution:
Step 1: Plot the poles of G(s)H(s) on the s-plane
       put s = jω
                               (jω+2)
        G(jω)H(jω) =
                            (jω+1)(jω−1)
                                                 ω
                               { ω2+4 ∠ tan−1 ( 2 )}
                =                                          ω
                    ω2+1 ∠ tan−1 ω { ω2+1 ∠ tan−1 (−1)}
                                                 ω
                                { ω2+4 ∠ tan−1 ( 2 )}
                =
                        ω2+1 2∠ tan−1 ω {∠180        − tan−1 ω}
              { ω2+4}
       M=               ;
               ω2+1
                      ω
       Ø = ∠ tan−1 ( ) − ∠ tan−1 ω − 180 + ∠ tan−1 ω
                      2
                   −1 ω
          =   ∠ tan ( )       − 180
                      2
lim 𝑀∠Ø = 2∠-180   (point a’)
ω⟶0
lim 𝑀∠Ø = 0∠-90    (point b’)
ω⟶∞
Section II : To find the image of path ‘bcd’
      put s = lim 𝑅𝑒𝑗θ in G(s)H(s)
              𝑅⟶∞
      Here , θ changes from +90 ⟶ 0 ⟶ -90θ
                          𝑗θ        (𝑅𝑒𝑗  +2)
      Then, lim 𝐺𝐻( 𝑅𝑒 ) = lim      θ       𝑗θ
             𝑅⟶∞                𝑅⟶∞ (𝑅𝑒𝑗 +1)(𝑅𝑒 −1)
                                         𝑅𝑒𝑗θ
                               = lim    𝑗θ    𝑗θ
                                 𝑅⟶∞ (𝑅𝑒 )(𝑅𝑒 )
                                       1
                               = lim     θ
                                 𝑅⟶∞ (𝑅𝑒𝑗 )
                               = 0∠-θ
                               = 0∠-90⟶ 0⟶ 90
                                    ⤒b’ ⤒c’ ⤒d’
Hence, the infinite semicircle ‘bcd’ on the s-
plane is mapped to the origin of the G(s)H(s)-
plane.
Section III: To find the image of path ‘da’
Path d’a’ is the mirror image of the path a’b’
with respect to real axis.
The complete Nyquist plot is shown below
Since, (-1+j0) point is encircled in anticlockwise
direction by the plot a’b’c’d’a’ in the GH(s) plane,
N = -1
Step 3:
The Nyquist stability criterion is Z = P + N
      Hence, Z = 1 - 1 = 0
      ⟹No roots of the system lie to the right-
half of s-plane. Hence the closed-loop control
system is stable.
Using Matlab
                Special Case
If the poles of G(s)H(s) lie at the origin of the
s- plane, then they are taken to the left-side of
the s- plane (right-side of the s- plane) by
drawing an indent ‘efa’ of radius, r⟶0 as
shown. Then find the number of encirclements
made by the image of the contour ‘abcdefa’
about (- 1+j0) point on the GH(s)- plane.
                      Problem
A negative feedback control system is characterized by
                                            5
an open-loop transfer function, G(s)H(s) =      .
                                             𝑠(𝑠+1)
Investigate the closed-loop stability of the system using
Nyquist stability criterion.
Solution
Step 1: Plot the poles of G(s)H(s) on the s-plane.
       The pole at the origin is taken to the left-side of
        the s-plane by drawing an indent of zero radius
        around this pole.
        Since the pole at the origin is taken to the left-
        side of the s-plane, P=0.
The contour ‘abcdefa’ that includes at the origin
to the left side of the s-plane
Step 2: To find N:
Section I : To find the image of path ab.
                            5
            G(s)H(s)   =
                         𝑠(𝑠+1)
                Put s = jω
                               5
         G(jω)H(jω) =
                           𝑗ω(𝑗ω+1)
                                      5
                       =
                           ω∠90 (ω2+1)∠ tan−1 ω
                                   5
                       =
                           ω (ω2+1)∠90+tan−1 ω
               5
      M=
           ω (ω2+1)
      Ø = −90 − tan−1 ω
lim 𝑀∠Ø = ∞∠-90    (point a’)
ω⟶0
lim 𝑀∠Ø = 0∠-180   (point b’)
ω⟶∞
Section II : To find the image of path ‘bcd’
      put s = lim 𝑅𝑒𝑗θ in G(s)H(s)
              𝑅⟶∞
      Here , θ changes from +90 ⟶ 0 ⟶ -90
                          𝑗θ         5
      Then, lim 𝐺𝐻( 𝑅𝑒 ) = lim 𝑗θ 𝑗θ
             𝑅⟶∞                𝑅⟶∞ 𝑅𝑒 (𝑅𝑒 +1)
                                           5
                               = lim    𝑗θ   𝑗θ
                                 𝑅⟶∞ (𝑅𝑒 )(𝑅𝑒 )
                                        5
                               = lim 2 𝑗2θ
                                 𝑅⟶∞ (𝑅 𝑒 )
                               = 0∠-2θ
                               = 0∠-180⟶ 0⟶ 180
                                      ⤒b’ ⤒c’ ⤒d’
Hence, the infinite semicircle ‘bcd’ on the s-
plane is mapped to the origin of the G(s)H(s)-
plane.
Section III: To find the image of path ‘de’
Path d’e’ is the mirror image of the path a’b’
with respect to real axis.
Section IV : To find the image of path efa
             put s = lim 𝑟𝑒𝑗∅ in G(s)H(s)
                    𝑟⟶0
      Here , ∅ changes from -90 ⟶ 0 ⟶ +90
                        𝑗∅          5
      Then, lim 𝐺𝐻( 𝑟𝑒 ) = lim 𝑗∅ 𝑗∅
             𝑟⟶0              𝑟⟶0 𝑟𝑒 (𝑟𝑒 +1)
                                        5
                                = lim 𝑗∅
                                  𝑟⟶0 (𝑟𝑒 )
                                 = ∞∠-∅
                                 = ∞∠90⟶ 0⟶ -90
                                      ⤒e’ ⤒f’ ⤒a’
The complete Nyquist plot is shown
Since -1+j0 point is not encircled by the plot ,
N=0
Step 3 : Z = P + N = 0 + 0 = 0
Hence the closed loop system is stable.
Matlab
Let us re work the problem by including the pole
at origin to the right-half of the s-plane and the
contour ‘abcdefa’
Step 1: P = 1
The mapping of sections ‘ab’, ’bcd’, and ‘de’ on
to the G(s)H(s)-plane remains same
Section IV :
Mapping of section ‘efa’:
put s = lim 𝑟𝑒𝑗∅ in G(s)H(s)
       𝑟⟶0
      Here , ∅ changes from -90 ⟶ -180 ⟶ -270
                             𝑗∅          5
            Then, lim 𝐺𝐻( 𝑟𝑒 ) = lim 𝑗∅ ∅
                   𝑟⟶0                𝑟⟶0 𝑟𝑒 (𝑟𝑒𝑗 +1)
                                         5
                               =   lim     ∅
                                   𝑟⟶0 (𝑟𝑒𝑗 )
                               = ∞∠-∅
                               = ∞∠90⟶ 180⟶ 270
                                   ⤒e’ ⤒f’  ⤒a’
The complete Nyquist plot is shown
Since -1+j0 point is encircled by the plot in
anticlockwise direction, N = -1
Step 3 : Z = 1 -1 = 0
Hence the closed loop system is stable.
                          Problem
                                                      𝑘
Sketch the Nyquist plot for G(s)H(s) =
                                                 𝑠(𝑠+1)(𝑠+2)
Find the range of k for closed-loop stability
Solution
Step 1: Plot the poles of GH(s) on the s-plane
                                    𝑘
      =                                              ω
            ω∠90 { (ω2+1)∠ tan−1 ω}{ (ω2+4)∠ tan−1 ( 2 )}
                                𝑘
       =                                         ω
            ω (ω2+1) (ω2+4)∠90+tan−1 ω+tan−1 ( 2 )
                 𝑘
    M=
           ω (ω2+1) (ω2+4)
                     −1            −1 ω
    Ø = −90 − tan        ω−     tan ( )
                                      2
lim 𝑀∠Ø = ∞∠-90    (point a’)
ω⟶0
lim 𝑀∠Ø = 0∠-270   (point b’)
ω⟶∞
Section II : To find the image of path ‘bcd’
      put s = lim 𝑅𝑒𝑗θ in G(s)H(s)
              𝑅⟶∞
      Here , θ changes from +90 ⟶ 0 ⟶ -90
                          𝑗θ             𝑘
      Then, lim 𝐺𝐻( 𝑅𝑒 ) = lim 𝑗θ 𝑗θ               θ
             𝑅⟶∞                𝑅⟶∞ 𝑅𝑒 (𝑅𝑒 +1)(𝑅𝑒𝑗 +2)
                                            𝑘
                               = lim    𝑗θ   𝑗θ   𝑗θ
                                 𝑅⟶∞ (𝑅𝑒 )(𝑅𝑒 )(𝑅𝑒 )
                                        𝑘
                               = lim 3 𝑗3θ
                                 𝑅⟶∞ (𝑅 𝑒 )
                               = 0∠-3θ
                               = 0∠-270⟶ 0⟶ 270
                                   ⤒b’   ⤒c’ ⤒d’
Hence, the infinite semicircle ‘bcd’ on the s-
plane is mapped to the origin of the G(s)H(s)-
plane.
Section III: To find the image of path ‘de’
Path d’e’ is the mirror image of the path a’b’
with respect to real axis.
Section IV : To find the image of path efa
             put s = lim 𝑟𝑒𝑗∅ in G(s)H(s)
                    𝑟⟶0
      Here , ∅ changes from -90 ⟶ 0 ⟶ +90
                        𝑗∅            𝑘
      Then, lim 𝐺𝐻( 𝑟𝑒 ) = lim 𝑗∅ 𝑗∅            ∅
             𝑟⟶0              𝑟⟶0 𝑟𝑒 (𝑟𝑒 +1)(𝑟𝑒𝑗 +2)
                                        𝑘
                                = lim 𝑗∅
                                  𝑟⟶0 (𝑟𝑒 )
                                 = ∞∠-∅
                                 = ∞∠90⟶ 0⟶ -90
                                      ⤒e’ ⤒f’ ⤒a’
The complete Nyquist plot is shown
To find M1 :
At point x’, phase = -180
                    −1         −1 ω
     ⟹ −90    − tan ω − tan ( )          = -180
                                  2
                    −1 3ω
                tan       2 = 0
                      2−ω
                 2 − ω2 = 0
                     ω = √2 rad/sec
                 M1 = │GH(jω)│ ω = √2
                                 𝑘             𝑘
                     =       2       2     =
                         ω (ω +1) (ω +4)       6
Since P is zero , N must be zero for Z to be zero.
N will be zero if and only if -1+j0 is not encircled
by the Nyquist plot
For N to be zero, M1 < 1
         𝑘
Hence,       <1
         6
      ⟹k < 6
Since k is always positive, for closed-loop
stability : 0 < k < 6
                     Problem
The open loop transfer function of a negative unity
                               𝑘(𝑠+3)(𝑠+5)
feed back system is given by                 .
                                (𝑠−2)(𝑠−4)
Find the range of k for closed − loop stability
Solution
Step 1 : Plot the poles of GH(s) on the s-plane
X X X X
                   Put s = jω
                                  𝑘(𝑗ω+3)(𝑗ω+5)
              G(jω)H(jω) =
                                  (𝑗ω−2)(𝑗ω−4)
                                          ω                       ω
                     𝑘{    𝜔2+9 ∠ tan−1 ( 3 )} { (ω2+25)∠ tan−1 ( 5 )}
              =                          ω                        ω
                   { (ω2+4)∠ tan−1 (−2)} { (ω2+16)∠ tan−1 (−4)}
                                                        ω                ω
                     𝑘{    𝜔2+9     ω2+25 } {∠ tan−1 3 + ∠ tan−1 ( 5 )}
               =                                    ω              ω
                   {      𝜔2+4    𝜔2+16 } {∠ tan−1 −2 + ∠ tan−1 (−4)}
     𝑘{    𝜔2+9       ω2+25 }
M=
      {    𝜔2+4       𝜔2+16 }
                  ω                 ω                 ω               ω
Ø = ∠ tan−1           + ∠ tan−1 ( )} - {∠ tan−1             + ∠ tan−1 ( )}
                  3                 5                 −2              −4
            ω                 ω                   ω                        ω
 = tan−1          + tan−1 ( )} - {180 − tan−1              + 180 − tan−1 ( )}
            3                 5                   2                        4
                          ω             ω              ω             ω
Ø = - 360 + tan−1               + tan−1 ( ) + tan−1         + tan−1 ( )
                          3             5              2             4
            15𝑘
lim 𝑀∠Ø =         ∠-360   (point a’)
ω⟶0          8
lim 𝑀∠Ø = k ∠0            (point b’)
ω⟶∞
Section II : To find the image of path ‘bcd’
                   𝑘(𝑗ω+3)(𝑗ω+5)
      G(s)H(s) =
                    (𝑗ω−2)(𝑗ω−4)
      put s = lim 𝑅𝑒𝑗θ in G(s)H(s)
              𝑅⟶∞
      Here , θ changes from +90 ⟶ 0 ⟶𝑗θ -90 θ
                          𝑗θ     𝑘(𝑅𝑒 +3)(𝑅𝑒𝑗 +5)
      Then, lim 𝐺𝐻( 𝑅𝑒 ) = lim       𝑗θ     θ
             𝑅⟶∞                𝑅⟶∞ (𝑅𝑒 −2)(𝑅𝑒𝑗 −4)
                                     𝑘(𝑅𝑒𝑗θ)(𝑅𝑒𝑗θ)
                               = lim     𝑗θ   𝑗θ
                                 𝑅⟶∞ (𝑅𝑒 )(𝑅𝑒 )
                               = lim k
                                   𝑅⟶∞
                               = k∠0
Section III: To find the image of path ‘de’
Path d’e’ is the mirror image of the path a’b’
with respect to real axis.
The Complete Nyquist plot is shown
To find M1 :
At point x’, phase = -180
                    ω            ω            ω             ω
⟹ - 360 +   tan−1
                    3
                        + tan−1 ( ) + tan−1
                                 5            2
                                                  + tan−1
                                                            4
                                                                = -180