Week11
Week11
(KON 313E)
Week 11
• The Nyquist criterion relates stability of a Closed Loop (CL) system to the Open Loop (OL) Frequency
Response (FR) and the OL poles and zeros
• The Nyquist criterion provides information on the transient and steady-state error.
Introduction to Nyquist Criterion
• OL system
𝑂𝐿 → 𝐺 𝑠 𝐻 𝑠
• CL system
𝐶𝐿 → 𝑇 𝑠 = 𝐺 𝑠 / 1 + 𝐺 𝑠 𝐻 𝑠
We obtain:
𝑁𝐺 𝑁𝐻
𝐺 𝑠 𝐻 𝑠 =
𝐷𝐺 𝐷𝐻
and thus
𝐷𝐺 𝐷𝐻 + 𝑁𝐺 𝑁𝐻
𝐹 𝑠 =1+𝐺 𝑠 𝐻 𝑠 =
𝐷𝐺 𝐷𝐻
𝐺 𝑠 𝑁𝐺 𝐷𝐻
𝐶𝐿𝑇𝐹 = =
1+𝐺 𝑠 𝐻 𝑠 𝐷𝐺 𝐷𝐻 + 𝑁𝐺 𝑁𝐻
If we take a complex number on the s-plane and substitute it into a function, F(s), another complex number results.
Let say the complex number to be mapped is s=4+3j and the mapping function 𝐹 𝑠 = 𝑠 2 + 2𝑠 + 1
Slide 9
Introduction to Nyquist Criterion
Some examples of contour mapping:
Slide 10
Introduction to Nyquist Criterion
𝐷𝐺 𝐷𝐻 +𝑁𝐺 𝑁𝐻
• Now, assume that the mapping function is 𝐹 𝑠 = 1 + 𝐺 𝑠 𝐻 𝑠 =
𝐷𝐺 𝐷𝐻
Introduction to Nyquist Criterion
➢ Each pole or zero of 1 + 𝐺 𝑠 𝐻 𝑠 whose vector undergoes a complete rotation of contour A (which
means inside contour A) must yield a change of 360° in the resultant, R, or a complete rotation of
contour B
➢ A zero inside a Clockwise (CW) contour A yields a CW rotation of contour B
➢ A pole inside a CW contour A yields a CCW rotation of contour B
➢ 𝑁 = 𝑃 − 𝑍
𝐷𝐺 𝐷𝐻 + 𝑁𝐺 𝑁𝐻
➢ N, # CCW rotations of contour B about the origin 𝐹 𝑠 =1+𝐺 𝑠 𝐻 𝑠 =
𝐷𝐺 𝐷𝐻
➢ P, # poles of 1 + G(s)H(s) inside contour A
➢ Z, # zeros of 1 + G(s)H(s) inside contour A
𝐷𝐺 𝐷𝐻 +𝑁𝐺 𝑁𝐻
➢ Now, lets extend the contour A to include the entire RHP for 𝐹 𝑠 = 1 + 𝐺 𝑠 𝐻 𝑠 = 𝐷𝐺 𝐷𝐻
➢ Z, # RHP CL poles
➢ CL stability!
➢ P, # RHP OL poles
➢ Easy ∞
➢ N, # CCW rotations of contour B about origin
➢Adjustment
➢Lets map 𝐺 𝑠 𝐻 𝑠 instead of 1 + 𝐺 𝑠 𝐻 𝑠
➢N, # CCW rotations of contour B about -1, instead of
origin.
Introduction to Nyquist Criterion
Nyquist stability criterion:
Z=P−N
# of CCW rev.
# of RHP # of RHP around −1
CL Pole OL Pole
• The contour that encloses the right half-plane can be mapped through the
function G(s)H(s) by substituting points along the contour into G(s)H(s).
• The points along the positive extension of the imaginary axis yield the polar
frequency response of G(s)H(s).
Z=P−N
# of CCW rev.
# of RHP # of RHP around −1
CL Pole OL Pole
Sketching the Nyquist Diagram
• Let us handle the following control system to show how to sketch the Nyquist diagram
Sketching the Nyquist Diagram
• We need to first find the open loop transfer function (G(s)H(s))
𝟓𝟎𝟎
𝑮 𝒔 =
(𝒔 + 𝟏)(𝒔 + 𝟑)(𝒔 + 𝟏𝟎)
• The sketch of the Nyquist diagram is then performed by complex arithmetic with vector
evaluations.
Sketching the Nyquist Diagram
• We need to first find the open loop transfer function (G(s)H(s))
500
𝐺 𝑠 𝐻(𝑠) =
(𝑠 + 1)(𝑠 + 3)(𝑠 + 10)
• Then the jw form is as follows
500
𝐺 𝑗𝑤 𝐻 𝑗𝑤 = 𝐺 𝑠 𝐻 𝑠 ቚ =
𝑠→𝑗𝑤 (𝑗𝑤 + 1)(𝑗𝑤 + 3)(𝑗𝑤 + 10)
• In polar coordinates:
500
𝐺 𝑗𝑤 𝐻 𝑗𝑤 =
𝑤 2 + 1 𝑤 2 + 9 𝑤 2 + 100
Note: For inverse tangent calculation, use always four-quadrant inverse tangent functions (For example
«atan2d» function in matlab)!
𝑤 𝑤
i.e. Define tan−1 as atan2d(w,3) or Define tan−1 as atan2d(w,-3) . This function results in four
3 −3
quadrant as expected!
Sketching the Nyquist Diagram
500
𝐺 𝑗𝑤 𝐻 𝑗𝑤 =
𝑤 2 + 1 𝑤 2 + 9 𝑤 2 + 100
𝑤→∞ 𝐺 𝑗∞ 𝐻 𝑗∞ = 0 ∠𝐺 𝑗∞ 𝐻 𝑗∞ = −270°
1
•𝐺 𝑗𝑤 𝐻 𝑗𝑤 = 𝐺 𝑠 𝐻 𝑠 ȁ𝑠→𝑗𝑤 =
𝑗𝑤+1
• In polar coordinates:
1 𝑤
•𝐺 𝑗𝑤 𝐻 𝑗𝑤 = ∠ − tan−1
𝑤 2 +1 1
𝑤 = 0−
• What to do now?
• We detour around the poles.
• Can detour to the right or left.
• If we detour to the left, then the poles count as
unstable open loop poles.
a. poles on contour; b. detour right; c. detour left
• P=3
• Assume we detour to the right.
• P=0
• Both can be done but We will detour right, in
our examples
Sketching the Nyquist Diagram
Nyquist diagram of the given unity feedback system , • At point A (w=0+), the two open-loop poles at the origin contribute
𝑆+2 2X90=180deg , and the zero contributes 0deg. The total angle at
where 𝐺 𝑠 = 𝑠2
point A is thus 0-180=-180deg
• We can use the Nyquist diagram to determine a system's stability, using the simple equation
𝑍 = 𝑃— 𝑁
• If the closed-loop system has a variable gain in the loop, one question we would like to ask is, "For
what range of gain is the system stable?“
• We have already two answers: (1) Root Locus and (2) Routh Hurwitz
• The system is marginally stable if the Nyquist diagram intersects the real axis at —1. Thus,
➢ If the open-loop system contains a variable gain, K, set K = 1 and sketch the Nyquist diagram.
➢ Adjust the value of K to yield stability, based upon the Nyquist criterion (scale the Nyquist plot by
1/𝐾 to touch the critical point -1.).
Stability via the Nyquist Diagram
𝑤 𝑤 𝑤
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = − tan−1+ tan −1
+ tan −1
0 3 5
𝑤 𝑤
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = −[90 + tan−1 + tan−1 ]
3 5
Stability via the Nyquist Diagram
𝐾
𝐺 𝑗𝑤 𝐻 𝑗𝑤 =
𝑤 𝑤 2 + 9 𝑤 2 + 25
−1
𝑤 −1
𝑤
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = −[90 + tan + tan ]
3 5
• Set 𝐾 = 1 and draw the Nyquist diagram using the contour A
??
Next, we need to find the point where the Nyquist diagram intersects the negative real axis.
Stability via the Nyquist Diagram
1
𝐺 𝑗𝑤 𝐻 𝑗𝑤 ቚ =
𝐾=1 (𝑗𝑤)(𝑗𝑤 + 3)(𝑗𝑤 + 5)
−8𝑤 2 − 𝑗(15𝑤 − 𝑤 3 )
𝐺 𝑗𝑤 𝐻 𝑗𝑤 ቚ =
𝐾=1 64𝑤 4 + 𝑤 2 15 − 𝑤 2 2
𝐺 𝑗𝑤 𝐻 𝑗𝑤 ቚ
𝐾=1
−8𝑤 2 −𝑗(15𝑤 − 𝑤 3 )
= +
64𝑤 4 + 𝑤 2 15 − 𝑤 2 2 64𝑤 4 + 𝑤 2 15 − 𝑤 2 2
Hence K can be increased by 1∕0.0083 (𝑲 = 𝟏Τ𝟎. 𝟎𝟎𝟖𝟑) before the Nyquist encircles -1. Thus,
• for stability 𝐾 < 120.5
• for marginal stability 𝐾 = 120.5
Stability via the Nyquist Diagram
𝑤 𝑤 𝑤
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = tan−1
− tan −1
+ (− tan −1
)
0.5 0 −2
−1
𝑤 −1
𝑤 Reminder: For inverse tangent calculation,
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = tan − 90 − tan ) use always four-quadrant inverse tangent
0.5 −2 functions (For example «atan2d» function in
𝑤 𝑤 matlab)!
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = tan−1 − 90 + tan−1
0.5 −2
𝑤 𝑤
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = −90 + tan−1 + tan−1
0.5 −2
Stability via the Nyquist Diagram
𝐾 𝑤 2 + 0.25 𝑤 = 0+
𝐺 𝑗𝑤 𝐻 𝑗𝑤 =
𝑤 𝑤2 + 4
𝑤 𝑤
∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = −90 + tan−1 + tan−1
0.5 −2
• Set 𝐾 = 1 and draw the Nyquist diagram
𝑤 = ⋯ → 𝐺 𝑗𝑤 𝐻 𝑗𝑤 = ⋯ ;
𝑤 = ⋯ → ∠𝐺 𝑗𝑤 𝐻 𝑗𝑤 = ⋯
𝑤 = ∞ → 𝐺 𝑗∞ 𝐻 𝑗∞ = 0
𝑤 = 0−
∠𝐺 𝑗∞ 𝐻 𝑗∞ = −90°
Stability via the Nyquist Diagram
• Next, we need to find the point where the Nyquist diagram intersects the negative real axis
−2.5𝑤 2 𝑤(1 − 𝑤 2 )
𝐺 𝑗𝑤 𝐻 𝑗𝑤 ቚ = 2 +𝑗 2
𝐾=1 (𝑤 ) (𝑤 2 + 4) (𝑤 ) (𝑤 2 + 4)
• If we set the imaginary part to zero
𝑤(1 − 𝑤 2 ) = 0
• We obtain 𝑤 = 1. Now if substitute this value to real part of 𝐺 𝑗𝑤 𝐻 𝑗𝑤 ȁ𝐾=1 , we obtain:
−2.5𝑤 2
อ = −0.5
(𝑤 2 ) (𝑤 2 + 4)
𝑤=1
Stability via the Nyquist Diagram
We know that P=1 thus for stability, N=1. Therefore, the system is stable if the critical point lies inside the
contour (N=1) so that Z=P-N=0
𝑤 = 0+
𝐾(𝑠 + 0.5)
𝐺 𝑠 𝐻(𝑠) =
𝑠(𝑠 − 2)
-0.5
𝑤 = 0−
Stability via the Nyquist Diagram
𝑤 = 0+ 𝐾(𝑠 + 0.5)
𝐺 𝑠 𝐻(𝑠) =
𝑠(𝑠 − 2)
-0.5
𝑤 = 0−
Stability via the Nyquist Diagram
• Hence K must increased by 𝟏Τ𝟎 . 𝟓 before the Nyquist encircles -1. Thus,
• for stability 𝐾 > 2
• for marginal stability 𝐾 = 2
𝐾(𝑠 + 0.5)
𝑤 = 0+ 𝐺 𝑠 𝐻(𝑠) =
𝑠(𝑠 − 2)
𝑤 = 0−
Stability via the Nyquist Diagram
𝑤 = 0+
𝐾(𝑠 + 0.5)
𝐺 𝑠 𝐻(𝑠) =
𝑠(𝑠 − 2)
𝑤 = 0−
Stability via the Nyquist Diagram
𝑤 = 0+
𝐾(𝑠 + 0.5)
𝐺 𝑠 𝐻(𝑠) =
𝑠(𝑠 − 2)
𝑤 = 0−