CHP 4
CHP 4
1
CHAPTER 4
FREQUENCY RESPONSE AND STABILITY ANALYSIS
4.1. Introduction
2
Cont. …
➢ If stable, linear, time-invariant system subjected sinusoidal input signal, the
steady-state output will also be a sinusoidal signal of the same frequency,
but with possibly different magnitude and phase angle of the input signal.
p(s) p(s)
➢ Let G s = = , and
q(s) s+s1 s+s2 …(s+sn )
3
Cont. …
𝐱 𝐭 = 𝐀 𝐬𝐢𝐧 𝛚𝐭 −is sinusoidal input signal and y t − is sinusoidal output
signal.
Y(s) Aω
= G s ⇒ Y s = G s X(s), X s =
X(s) s2 +ω2
Aω
Y s =G s X s = G(s)
s2 +ω2
k ഥ
k
Y s = + + partial fraction of G s term
s+jω s−jω
due to transfer function
due to input ↑
↑ transient response
steady state response (yss ) (it will die out as time increases
so we can neglect it)
4
Cont. …
➢ For a stable system, s1 , s2 , … , sn have negative real parts. Therefore, as time
approaches infinity, the terms e−s1 t , e−s2 t , … , e−sn t approach zero.
k ഥ
k
Yss (s) = + where, kത − complex conjugate of k.
s+jω s−jω
A G(jω) A G(jω)
Yss t = −e−jϕ 𝑒 −𝑗𝜔𝑡 + ejϕ 𝑒 𝑗𝜔𝑡 = −𝑒 −𝑗(𝜔𝑡+ϕ) + 𝑒 𝑗(𝜔𝑡+ϕ)
2j 2j
𝑒 𝑗(𝜔𝑡+ϕ) −𝑒 −𝑗(𝜔𝑡+ϕ)
Yss t = A G(jω) = A G(jω) sin 𝜔𝑡 + ϕ ,
2𝑗
𝐞𝐣𝛉 − 𝐞−𝐣𝛉
= 𝐣𝐬𝐢𝐧𝛉
𝟐
6
Cont. …
➢ A positive phase angle is called phase lead while negative phase angle is
called phase lag. 7
Cont. …
➢ Consider the network given by
1
s+
T1
G s =
1
s+
T2
For the sinusoidal input, x t = k sin ωt
a) Find the steady state output.
b) Determine whether this network is a lead network or lag network.
1
jω+T T2 (1+jωT1 ) T2 1+T2 1 ω2
Solution: G jω = 1
1
= , G jω = ,
jω+T T1 (1+jωT2 ) T1 1+T2 1 ω2
2
b) IF T1 > T2 , then tan−1 ωT1 > tan−1 ωT2 , so the network is lead network. Lag
network if T1 < T2 . 8
Cont. …
❖ Advantages of frequency response analysis
➢ Can infer performance and stability from same plot. It is very easy to
visualize the effects of noise disturbance and parameter variations.
11
4.3. Frequency domain specifications
➢ Consider the negative feedback system shown in Figure 4.2.
Fig.4.2.
C(s) G(s)
Gcl s = =
R(s) 1 + H s G(s)
12
Cont. …
G(jω)
Gcl jω =
1 + H jω G(jω)
❖ Resonant peak (𝐌𝐫 ): The maximum value of Gcl (jω) as ω is varied is termed
as resonance peak (Mr ). It is the peak (maximum) value of Gcl jω .
Mr = Gcl jω |ω=ωr
13
Cont. …
❖ Resonant frequency (𝝎𝒓 ): It is the frequency at which the magnitude of the
frequency response has peak value for the first time. It can obtain from
𝑑
Gcl jω 𝜔=𝜔𝑟 =0
𝑑𝜔
1
Gcl jω 𝜔=𝜔𝑏 = = 0.707
2
➢ Small bandwidth means lower frequencies are passed and hence response
is sluggish.
14
Cont. …
➢ If bandwidth is large, higher frequencies will pass through the system. This
indicates a faster rise in time. Large bandwidth indicates susceptibility to
noise.
❖ Cut-off rate: It is the slop of the log-magnitude curve near the cut-off
frequency.
➢ The factor by which the gain H(jω)G(jω) dB increases without making the system
unstable (the gain required to lift magnitude curve to 0dB).
1
ii. Determine, Gm = ; In decibels, Gm │dB = 20 log Gm = −20 log H(jω)G(jω) │ω=ωpc
H(jω)G(jω) │ω=ωpc
➢ The factor by which the phase in degree can be increased without making the
system unstable. It is a measure of relative stability.
Pm = 180° + ϕ jω |ω=ωgc
➢ Drive the resonant frequency (𝜔𝑟 ), resonant peak (𝑀𝑟 ) and bandwidth (𝜔𝑏 ).
2 −1Τ2
𝑑 𝑑 ω2 2ξω 2
G(jω) 𝜔=𝜔𝑟 =0 ⇒ 1− + │𝜔=𝜔𝑟 = 0
𝑑𝜔 𝑑𝜔 ω2 n ωn
18
Cont. …
ω
For simplicity let, u =
ωn
d 2 −1Τ2 │ ωr
1 − u2 2 + 2ξu u=ur = 0, where ur =
du ωn
1 2 −3Τ2
− 1 − u2 2
+ 2ξu 2 1 − u2 −2u + 2 2ξu 2ξ │u=ur = 0
2
1 2 −3Τ2
− 1 − u2 2 + 2ξu 4u u2 − 1 + 2ξ2 │u=ur = 0
2
−4u u2 −1+2ξ2
2 3Τ2 |u=ur = 0 ⇒ 4u u2 − 1 + 2ξ2 |u=ur = 0
2 1−u2 2 + 2ξu
ωr
ur = → ωr = ωn ur = ωn 1 − 2ξ2 19
ωn
Cont. …
1 1
Mr = G(jω) |ω=ωr = |ω=ωr = ,
2 2
ω2 2ξω 2 ω2 r 2ξωr 2
1− 2 + ω 1− 2 + ω
ω n n ω n n
ii. Mr is greater than unity and ωr is less than ωn for 0 < ξ < 0.707
iii. There is no resonant peak and the maximum value of Mr is unity and ωr = 0 and
ξ > 0.707
20
Cont. …
1 1 1
G(jω) |ω=ωb = = 0.707 ⇒ |ω=ωb =
2 2 2
ω2 2ξω 2
1− 2 + ω
ω n n
2
1 1 𝜔2 𝑏 2ξωb 2 ωb
= ⇒ 1− + = 2, let 𝑢𝑏 =
2 2 ω2 n ωn ωn
𝜔2 2ξωb 2
1− 2 𝑏 + ωn
ω n
21
Cont. …
ωb
ub = ⇒ 𝜔𝑏 = 𝜔𝑛 𝑢𝑏 = ωn 1 − 2ξ2 ± 4ξ4 − 4ξ2 + 2
ωn
i. ωb ∝ ωn
iii. For 0 < ξ < 0.707, bandwidth and Mr are directly proportional.
23
4.5. Bode plot and Bode stability criterion
Bode plot
➢ Bode diagram consists of two graphs:
available.
𝐺 𝑗𝜔 = 𝐺 𝑠 |𝑠=𝑗𝜔
25
Cont. …
i. Constant gain k
±
ω2 j2ζω
iv. Quadratic factor (i.e., quadratic pole or zero) 1 − 2 +
ωn ωn
26
Cont. …
➢ Once we become familiar with the logarithmic plots of these basic factors, it
is possible to utilize them to construct a composite logarithmic plot for any
general form of H jω G(jω) by sketching the curves for each factor and
adding individual curves graphically, because adding the logarithms of the
gains corresponds to multiplying them together.
k T1 s + 1 T2 s + 1 … Tn s + 1
G s =
p1 s + 1 p2 s + 1 … pn s + 1
27
Cont. …
➢ The sinusoidal transfer function is given by
tan−1 p2 ω − ⋯ − tan pn ω
28
Cont. …
❖ Bode plot of constant gain 𝐤: Assume unity feedback system
Fig.4.3: Bode plot of constant gain k, (a) Magnitude plot and (b) phase plot.
29
Cont. …
❖ Bode plot of Derivative or Integral factors (zero or pole at origin) 𝒋𝝎 ±:
If 𝜔 = 1, 𝐺 𝑗𝜔 |𝑑𝐵 = 0
If 𝜔 = 10, 𝐺 𝑗𝜔 |𝑑𝐵 = 20 𝑑𝐵 Τ𝑑𝑒𝑐𝑎𝑑𝑒
1 1
➢ Consider open-loop transfer function G s = ⇒ G jω =
𝑠 𝑗𝜔
1
𝐺 𝑗𝜔 |𝑑𝐵 = 20 log = −20 log ω and ϕ 𝑗𝜔 = −90°, as 𝜔 varies from 0 to ∞.
𝜔2
If 𝜔 = 1, 𝐺 𝑗𝜔 |𝑑𝐵 = 0
If 𝜔 = 10, 𝐺 𝑗𝜔 |𝑑𝐵 = −20 𝑑𝐵 Τ𝑑𝑒𝑐𝑎𝑑𝑒
30
Cont. …
1
Fig.4.4: Bode diagram of (a) 𝐺 𝑗𝜔 = jω, (b) 𝐺 𝑗𝜔 =
𝑗𝜔 31
Cont. …
➢ In Bode diagrams, frequency ratios are expressed in terms of octaves or decades. An
octave is a frequency band from ω1 to 2ω1 and decade is a frequency band from ω1 to
10ω1 , where ω1 , is any frequency value.
1
➢ The slopes of the log-magnitude curves for the factors jω and are 20 𝑑𝐵 Τ𝑑𝑒𝑐𝑎𝑑𝑒
jω
and −20 𝑑𝐵 Τ𝑑𝑒𝑐𝑎𝑑𝑒 respectively.
1 𝑛
➢ Note: If the transfer function contains the factor jω 𝑛 and , the log magnitude
jω
becomes, respectively,
1
20 log jω 𝑛
= 20n log jω = 20𝑛 log 𝜔 𝑑𝐵 and 20 log = −20n log jω = −20𝑛 log 𝜔 𝑑𝐵
jω 𝑛
For ωT ≪ 1, ϕ 𝑗𝜔 = tan−1 0 = 0°
1 1
➢ Consider open-loop transfer function G s = ⇒ 𝐺 jω =
Ts++1 1+jωT
1 1
G jω |dB = 20 log = 20 log = −20 log 1 + T 2 ω2 𝑑𝐵
1+jωT 1+T2 ω2
➢ Cut off frequency (corner frequency or break frequency) is the frequency at which the
For ωT ≪ 1, ϕ 𝑗𝜔 = − tan−1 0 = 0°
𝑗𝜔 𝑗𝜔 2 𝜔2 𝑗2𝜁𝜔
➢ Consider second order zero 𝐺 𝑗𝜔 = 1 + 2𝜁 + =1− +
𝜔𝑛 𝜔𝑛 𝜔𝑛 2 𝜔𝑛
• The magnitude and phase of the quadratic factor depend on both the corner
frequency and the damping ratio.
• Damping ratio determine magnitude resonant peak. It is large for small value of ζ.
37
Cont. …
2
𝜔2 𝑗2𝜁𝜔 𝜔 2 2𝜁𝜔 2
G jω |dB = 20 log 1 − + = 20 log 1− + 𝑑𝐵
𝜔𝑛 2 𝜔𝑛 𝜔𝑛 𝜔𝑛
𝜔 2
For 𝜔 → 0 𝜔 ≪ 𝜔𝑛 ⇒ ≪ 1 , G jω |dB = 20 log 1 = 20 log 1 = 0 𝑑𝐵
𝜔𝑛
𝜔 2 𝜔 4
For high frequency 𝜔 ≫ 𝜔𝑛 , 4𝜁 2 ≪ because 𝜁 is very low.
𝜔𝑛 𝜔𝑛
𝜔 4 𝜔
G jω |dB = 20 log = 40 log dB (this the straight line of slope 40 𝑑𝐵Τ𝑑𝑒𝑐𝑎𝑑𝑒.
𝜔𝑛 𝜔𝑛
corner frequency is 𝜔 = 𝜔𝑛 , G jω |dB = 0 dB).
2𝜁𝜔
−1 𝜔𝑛
ϕ 𝑗𝜔 = tan
𝜔 2
1−
𝜔𝑛
0
For 𝜔 → 0, ϕ 𝑗𝜔 = tan−1 = 0°
1
38
Cont. …
𝜔 2𝜁 2𝜁
For 𝜔 = 𝜔𝑛 ( = 1), ϕ 𝑗𝜔 = tan−1 = tan−1 = tan−1 ∞ = 90°
𝜔𝑛 1−1 0
𝜔
For → ∞, ϕ 𝑗𝜔 = 180°
𝜔𝑛
𝜔
For > 1, ϕ 𝑗𝜔 − is becomes negative. Thus after this, the angle contribution
𝜔𝑛
1 1
❖ Consider second order zero 𝐺 𝑗𝜔 = 𝑗𝜔 𝑗𝜔 2
= 𝜔2 𝑗2𝜁𝜔
1+2𝜁 + 1− 2 +
𝜔𝑛 𝜔𝑛 𝜔𝑛 𝜔𝑛
39
Cont. …
2
1 𝜔 2 2𝜁𝜔 2
G jω |dB = 20 log 𝜔2 𝑗2𝜁𝜔
= −20 log 1− + 𝑑𝐵
1− 2 + 𝜔 𝜔𝑛 𝜔𝑛
𝜔𝑛 𝑛
𝜔
For ≪ 1, G jω |dB = −20 log 1 = −20 log 1 = 0 𝑑𝐵
𝜔𝑛
𝜔 𝜔 2 𝜔 4
For high frequency 𝜔 ≫ 𝜔𝑛 ( ≫ 1), 4𝜁 2 ≪ because 𝜁 is very low.
𝜔𝑛 𝜔𝑛 𝜔𝑛
𝜔 4 𝜔
G jω |dB = −20 log = − 40 log (this the straight line of slope
𝜔𝑛 𝜔𝑛
− 40 𝑑𝐵Τ𝑑𝑒𝑐𝑎𝑑𝑒. Corner frequency is 𝜔 = 𝜔𝑛 , G jω |dB = 0 dB.
2𝜁𝜔
𝜔𝑛
ϕ 𝑗𝜔 = − tan−1 𝜔 2
1−
𝜔𝑛
0
For 𝜔 → 0, ϕ 𝑗𝜔 = tan −1 = 0°
1
40
Cont. …
𝜔 2𝜁 2𝜁
For 𝜔 = 𝜔𝑛 ( = 1), ϕ 𝑗𝜔 = − tan−1 = − tan−1 = − tan−1 ∞ = −90°
𝜔𝑛 1−1 0
𝜔
For → ∞, ϕ 𝑗𝜔 = −180°
𝜔𝑛
𝜔
For > 1, ϕ 𝑗𝜔 − is becomes positive. Thus in such case the angle contribution
𝜔𝑛
𝜔
Example: = 2, 𝜁 = 0.3, ϕ 𝑗𝜔 = − tan−1 −0.4 = 21.8°. Thus the actual angle is
𝜔𝑛
41
Cont. …
➢ Then identify the corner frequencies associated with these basic factors.
➢ The exact curve, which lies close to the asymptotic curve, can be obtained
by adding proper corrections.
For 𝜔 → 0, 𝐺 𝑗𝜔 𝑑𝐵 = 0 dB
For 𝜔 = 4, 𝐺 𝑗𝜔 𝑑𝐵 = 20 log 1 + 0.25 × 4 2 = 20 log 2 = 3 dB
𝜔 = 4 − is corner frequency
44
Cont. …
1 1
➢ Pole: G jω = and G jω = , G jω dB = −20 log 1 + 0.5ω 2 and
j0.5ω+1 j0.1ω+1
45
Cont. …
46
Cont. …
ϕ 𝑗𝜔 = tan−1 0 + tan−1 0.25𝜔 − tan−1 0.5𝜔 − tan−1 0.1𝜔
47
Cont. …
2) Draw the Bode diagram for the unity feedback system with forward transfer
100(𝑠+1)
function given by: 𝐺 𝑠 =
𝑠 2 +110𝑠+1000
100(𝑠+1) 100(𝑠+1)
Solution:H s 𝐺 𝑠 = 𝐺(𝑠) = =
𝑠 2 +110𝑠+1000 (𝑠+10)(𝑠+100)
➢ Zero: 𝐺 𝑗𝜔 = 𝑗𝜔 + 1, 𝐺 𝑗𝜔 𝑑𝐵 = 20 log 1 + 𝜔 2
For 𝜔 → 0, 𝐺 𝑗𝜔 𝑑𝐵 = 0 dB
For 𝜔 = 1, 𝐺 𝑗𝜔 𝑑𝐵 = 20 log 1 + 1 2 = 20 log 2 = 3 dB
𝜔 = 1 − is corner frequency 48
Cont. …
1 1
➢ Pole: 𝐺 𝑗𝜔 = and 𝐺 𝑗𝜔 = , 𝐺 𝑗𝜔 𝑑𝐵 = −20 log 1 + 0.1𝜔 2
(𝑗0.1𝜔+1) (𝑗0.01𝜔+1)
49
Cont. …
ϕ 𝑗𝜔 = tan−1 0 + tan−1 𝜔 − tan−1 0.1𝜔 − tan−1 0.01𝜔 (deg)
50
Cont. …
3) Draw the Bode Diagram for the unity feedback system with forward transfer
120𝑆
function given by: 𝐺 𝑠 =
(𝑆+3)(𝑆+6)
120𝑆
Solution: H s 𝐺 𝑠 = 𝐺(𝑠) =
(𝑆+3)(𝑆+6)
120 𝑗𝜔 120(𝑗𝜔) 𝑗𝜔
𝐺 𝑗𝜔 = 𝐺 𝑠 |𝑠=𝑗𝜔 = = 𝑗𝜔 𝑗𝜔 = 20Τ3 𝑗𝜔 𝑗𝜔
(𝑗𝜔+3)(𝑗𝜔+6) 3×6( +1)( +1) ( +1)( +1)
3 6 3 6
1 1 𝜔2
➢ Pole: 𝐺 𝑗𝜔 = 𝑗𝜔 and 𝐺 𝑗𝜔 = 𝑗𝜔 , 𝐺 𝑗𝜔 𝑑𝐵 = −20 log 1 + and
( 3 +1) ( 6 +1) 9
𝜔2
𝐺 𝑗𝜔 𝑑𝐵 = −20 log 1 +
36
−1 𝜔 −1 𝜔
• ϕ 𝑗𝜔 = 90° − tan − tan
3 6
52
Cont. …
53
Cont. …
4) Draw the Bode Diagram for the unity feedback system with the forward transfer
120𝑠 2
function given by: 𝐺 𝑠 =
(𝑆+3)(𝑆+6)
100𝑠 2
Solution: H(s)𝐺 𝑠 =
(𝑆+3)(𝑆+6)
1 1 𝜔2
➢ Pole: 𝐺 𝑗𝜔 = 𝑗𝜔 and 𝐺 𝑗𝜔 = 𝑗𝜔 , 𝐺 𝑗𝜔 𝑑𝐵 = −20 log 1 + and
( 3 +1) ( 6 +1) 9
𝜔2
𝐺 𝑗𝜔 𝑑𝐵 = −20 log 1 +
36
−1 𝜔 −1 𝜔
• ϕ 𝑗𝜔 = 180° − tan − tan
3 6
55
Cont. …
56
Cont. …
Exercise 4.2: Draw the bode diagram of the following unity feedback system with
10(s+3) 20
a) G s = b) 𝐺 𝑠 =
s s+2 s2 +s+2 𝑠2 𝑠+3 𝑠+4
s(s+2)
e) G s =
s+2 s2 +4s+3
Exercise 4.3: Give your own examples and show clearly how the transfer function
is obtained from bode plot of the system.
57
Cont. …
➢ If the phase margin is positive then the system is stable otherwise unstable.
➢ For the system to be stable both gain margin and phase margin needs to be
greater than zero.
58
Cont. …
Bode stability requirement
❖ For a given open-loop transfer function (𝐻 𝑠 𝐺(𝑠)), the closed-loop system is stable
if the following conditions are met.
➢ Minimum Phase System: has no poles located on or to the right of the imaginary
axis, except of a single pole at the origin.
➢ Assume that the open-loop frequency response has only a single phase crossover
frequency and a single gain crossover frequency.
❖ Consider the bode diagram of a unity feedback system shown in Figure 4.7 below
59
Cont. …
Exercise 4.4: For a unity feedback system with forward transfer function
50
𝐺 𝑠 =
𝑠(0.4𝑠+1)
a) Find the gain crossover frequency (𝜔𝑔𝑐 ) and phase margin (𝑃𝑚 ),
b) Find ωn , ξ, ωr , Mr & ωb
61
4.6. Log-magnitude-versus-phase plots (Nichols Plots) and Stability
❖ Log-magnitude-versus-phase plots
➢ It is a plot of the logarithmic magnitude in decibels versus the phase angle or phase
margin for a frequency range of interest.
➢ The phase margin is the difference between the actual phase angle ϕ jω and −180°, that
is, ϕ jω − −180° = 180° + ϕ jω .
➢ Advantages of the log-magnitude vs phase plot are the relative stability of the closed-loop
system can be determined quickly and compensation can be worked out easily. 62
Cont. …
Example: For the following sinusoidal transfer function sketch the log-magnitude
in decibels versus phase plot in degree.
1 1
a) G jω = b) G jω = c) G jω = 1 + jωT
jω 1+jωT
1 1
d) G jω = e−jωL e) G jω = f) G jω = ω2 2ξω
for 𝛏 > 0
jω(1+jωT) 1− 2 +j
ω n ωn
Solution:
1 1
a) G jω dB = 20 log = 20 log = −20 log ω and ϕ jω = −90°
ω2 ω
If ω → 0, G jω dB → +ve
If ω = 1, G jω dB =0
If ω → ∞, G jω dB → −ve
63
Cont. …
1
b) G jω dB = 20 log = −20 log 1 + (ωT)2 and ϕ jω = − tan−1 ωT
1+(ωT)2
If ω = 0, G jω dB = 0 dB and ϕ jω = 0°
If ωT = 1, G jω dB = −3 dB and ϕ jω = −45°
If ω → ∞, G jω dB → −∞ and ϕ jω = −90°
64
Cont. …
d) G jω dB = 20 log 1 = 0 dB and ϕ jω = −ω𝐿
If ω = 0, G jω dB = 0 dB and ϕ jω = 0°
If ω → ∞, G jω dB = 0 dB and ϕ jω = −∞
65
Cont. …
2ξω
1 ω2 2 2ξω 2 ωn
f) G jω dB = 20 log = −20 log 1− +( ) and ϕ jω = − tan−1 ω2
2 ω2 n ωn 1− 2
ω2 2ξω 2 ω n
1− 2 +( )
ω n ωn
If ω = 0, G jω dB = 0 dB and ϕ jω = 0°
If ω = ωn , G jω dB = −20 log 2𝜉 𝑑𝐵 and ϕ jω = −90°
𝜔
If → ∞, G jω dB = −∞ and ϕ jω = −180°
𝜔𝑛
If ω = ωr , G jω dB = Mr
66
Cont. …
➢ If both gain margin and phase margin are positive the system is stable otherwise
unstable.
67
Cont. …
Exercise 4.5: For the following open-loop transfer functions sketch log-magnitude
versus phase plot.
2 50
a) G s = b) G s =
s(s+1)(s+2) s(s+5)
68
4.6. Polar Plot and stability
❖ Polar plot
➢ Thus, the polar plot is the locus of vectors G(jω) ∠G(jω) as ω is varied from
zero to infinity.
➢ One disadvantage is that the plot does not clearly indicate the contributions
of each individual factor of the open-loop transfer function.
➢ Given 𝐺(𝑠)
i. Obtain 𝐺 𝑗𝜔 = 𝐺(𝑠)│𝑠=𝑗𝜔
70
Cont. …
iii. Determine 𝐺 𝑗𝜔 as well as ∠𝐺 𝑗𝜔 at frequency 𝜔 = 0 and 𝜔 = ∞ using limit.
𝐺 𝑗𝜔 𝜔→0 𝐺 𝑗𝜔 𝜔→∞
ቊ and ቊ
∠𝐺 𝑗𝜔 𝜔→0 ∠𝐺 𝑗𝜔 𝜔→∞
1
a) G jω = b) G jω = jω
jω
1
c) G jω = d) G jω = 1 + jωT
1+jωT
e−jωL 1
e) G jω = f) G jω = jω 2
, for ξ > 0
1+jωT jω
1+2ξ +
ωn ωn
Solution:
1 1 −jω 1
a) G jω = = × = −j
jω jω −jω ω
1 1
G jω = −j = ∠ − 90°
ω ω
lim G jω = −𝑗∞ = ∞∠ − 90°
𝜔→0
lim G jω = −𝑗0 = 0∠ − 90°
𝜔→∞
72
Cont. …
b) G jω = jω = ω∠90°
lim G jω = 𝑗0 = 0∠90°
𝜔→0
lim G jω = 𝑗∞ = ∞∠90°
𝜔→∞
1
c) G jω =
1+jωT
1 1−jωT
G jω = ×
1+jωT 1−jωT
1 ωT 1 −1 ωT
G jω = − j = ∠ − tan
1+ ωT 2 1+ ωT 2 1+ ωT 2
lim G jω = 1 − j0 = 1∠0°
ω→0
lim G jω = 0 − j0 = 0∠ − 90°
ω→∞
1
lim G jω = 1Τ2 − j 1Τ2 = ∠ − 45° 73
ωT=1 2
Cont. …
d) G jω = 1 + jωT
G jω = 1 + jωT = 1 + ωT 2 ∠ tan−1 ωT
lim G jω = 1 + j0 = 1∠0°
ω→0
lim G jω = 1 + j∞ = ∞∠90°
ω→∞
e−jωL
e) G jω =
1+jωT
1 1−jωT
G jω = e−jωL ×
1+jωT 1−jωT
1 ωT 1
G jω = −j e−jωL = ∠ − ωL − tan −1 ωT
1+ ωT 2 1+ ωT 2 1+ ωT 2
lim G jω = 1 − j0 = 1∠0°
ω→0
lim G jω = 0 − j0 = 0∠ − ∞°
ω→∞
74
Cont. …
1
f) G jω = jω jω 2
, for ξ > 0
1+2ξ +
ωn ωn
ω 2 ω
1 1 1− −j2ξ
ωn ωn
G jω = ω 2 ω
= ω 2 ω
× ω 2 ω
1− ωn
+j2ξωn
1− ωn
+j2ξωn
1− ωn
−j2ξ ωn
ω 2 ω
1− ω 2ξω
n n
G jω = 2 −j 2
ω 2 ω 2 ω 2 ω 2
1− + 2ξω 1− + 2ξω
ωn n ωn n
ω
1 2ξω
G jω = ∠ − tan−1 n
ω 2
2 1−
ω 2 ω 2 ωn
1− ωn
+ 2ξω
n
1 1
a) G s = b) G s =
s(Ts+1) T1 s+1 T2 s+1
1 k
c) G s = d) G s =
T1 s+1 T2 s+1 (T3 s+1) s(Ts+1)
1 1
e) G s = f) G s =
s2 T1 s+1 T2 s+1 s+4 s+2
jω jω 2 20
g) G jω = 1 + 2ξ + , for ξ > 0 h) 𝐺 𝑠 =
ωn ωn 𝑠(𝑠+1)(𝑠+2)
20
i) G s =
s(s+1)
76
Cont. …
❖ Polar plot stability criterion
➢ Gain margin and phase margin can be determined from, the polar plot by drawing
unit circle.
➢ In the polar plot, a line may be drawn from the origin to the point at which the unit
circle crosses the 𝐺(𝑗𝜔)1ocus. The angle from the negative real axis to this line is
the phase margin. If the phase margin is positive the system is stable and if phase
margin is negative the system is unstable.
Pm = 180° + ϕ jω │ω=ωgc
➢ The gain margin expressed in decibels is positive if 𝐺𝑚 > 1 or G(jω) │𝜔=𝜔𝑝𝑐 < 1
➢ Thus, a positive gain margin (in decibels) means that the system is stable, and a
negative gain margin (in decibels) means that the system is unstable.
Example: Consider the polar plot shown below with unit circle.
78
Cont. …
➢ The point ‘A’ corresponds to phase crossover frequency (i.e., ϕ jω = −180°), while
point ‘B’ corresponds to gain crossover frequency (i.e., G(jω) = 1).
i. The plot may not cross real axis (point ‘A’ doesn’t exists).
➢ If these small lags are accounted for, the so-called first or second-order systems
may become unstable.
𝐺𝑚 │𝑑𝐵 = −20 log G jω = −20 log 1Τ2 = 20 log 2, Since 𝐺𝑚 > 0, the system is stable.
80
Cont. …
iii. The plot may cross 𝑅𝑒 −axis at −1 + 𝑗0. Then G(jω) = 1
𝐺𝑚 │𝑑𝐵 = −20 log G jω = −20 log 1 = 0 dB , the system is marginally stable.
iv. The plot may cross 𝑅𝑒 −axis beyond −1 + 𝑗0.
Example: The plot cross 𝑅𝑒 −axis at −2 + 𝑗0 , then G(jω) = 2
𝐺𝑚 │𝑑𝐵 = −20 log G jω = −20 log 2 = −6 dB, Since 𝐺𝑚 < 0, the system is unstable.
Exercise 4.7: Identify which of the following system is stable and unstable.
81
Cont. …
Exercise 4.8: For the system shown below by block diagram, sketch the
polar plot for k = 30, 60 & 120 and from polar plot determine whether the
closed-loop system is stable for each value of this gain.
82
4.8. Nyquist plot and Nyquist stability criterion
❖ Nyquist plot
➢ Nyquist plot is the mirror image of polar plot. It is based on polar plot.
k
Example: Draw the Nyquist plot of G s =
s(Ts+1)
k
G jω =
jω jωT+1
k −jω 1−jωT
G jω = ×
jω 1+jωT −jω 1−jωT
kT k k −1 ωT
G jω = − −j = ∠ − 90° − tan
1+ ωT 2 ω 1+ ωT 2 ω 1+ ωT 2
Y(s) G(s)
=
R(s) 1 + H s G(s)
➢ For stability of the closed-loop system all the roots of the characteristic
equation must lie to the left half of s-plane.
85
Cont. …
➢ The characteristic equation is
k s + z1 s + z2 … s + zm
𝑄 𝑠 =1+H s G s =1+ =0
s + p1 s + p2 … s + pn
s + p1 s + p2 … s + pn + 𝑘 s + z1 s + z2 … s + zm
𝑄 𝑠 = =0
s + p1 s + p2 … s + pn
➢ From equation (4.1) and (4.2), we find that the poles of Q(s) are the same as the
poles of open-loop transfer function and zeros of Q(s) are the root of characteristic
equation. 86
Cont. …
➢ The system is stable if all the roots of characteristic equation are lie in the
➢ Therefore for the system to be stable the zeros of 𝑄(𝑠) lie in the left half of s-
plane.
➢ Segment 𝐶2 and 𝐶3 cover the whole imaginary axis from −𝑗∞ to 𝑗∞.
➢ Segment 𝐶1 is semicircle of radius infinity that encloses the entire right half of s-
plane.
88
Cont. …
➢ Mapping of C into Q(s) is obtained for segment C2 and C3 by substituting s = jω in
Q(s) and varying ω from 0 to j∞ and −j∞ to 0 respectively.
➢ Segment C1 can be mapped into Q(s) plane by substituting s = Rejθ and varying θ
from + πΤ2 to − πΤ2 and radius approach infinity.
❖ With mapping Q(s) of closed path C, we can find the following on the basis of
principle of argument.
1) If the total number of zeros of Q(s) in the right half of s-plane is equal to z, then the
number of clockwise encirclements of the origin of Q s −plane is z, hence the
system is unstable (i.e. there is no open-loop poles (𝑝 = 0) in the right half -plane).
2) If the total number of poles of Q(s) in the right half of s-plane is equal to p, then the
number of counterclockwise encirclements of the origin of Q s −plane is p, hence
the system is stable (i.e. there is no zeros of Q(s) (𝑧 = 0) in the RH s-plane ).
89
Cont. …
3) If the difference between number of poles and zeros of Q(s) in the right half of s-
plane is N ( i.e., N = z − p ), then the number of encirclements of the origin of
Q s − plane is N.
➢ For closed-loop system to be stable the number of zeros of Q(s) that lie to the right
of s-plane is zero. That is z = 0, which gives N = −p. Therefore the system is stable,
the number of counterclockwise encirclements of origin Q s −plane must be
equal to p.
⇒ H s G s = −1 (4.3)
Example: when a system has one open-loop pole (p = 1) in the right half of
s-plane and has one clockwise encirclement of −1 + j0, i.e., N = 1, then
number of closed-loop pole is two (z = 2) in the right half of s-plane. N = z −
p ⇒ z = N + P = 1 + 1 = 2, hence the system is unstable.
93
Cont. …
➢ If open-loop pole is one ( 𝑝 = 1 ) in the right s-plane and has one
counterclockwise encirclement of −1 + 𝑗0, i.e., 𝑁 = −1, 𝑁 = 𝑧 − 𝑝 ⇒ 𝑧 =
𝑁 + 𝑝 = −1 + 1 = 0. Therefore there is no closed-loop poles in the right half
of s-plane , hence the system is stable but open-loop is unstable.
➢ In modified closed contour/path only a very small portion of the plane is not
encircled. 94
Cont. …
➢ This small portion approaches zero as the radius of semicircle 𝜖 → 0, thus
ensuring induction of all poles and zeros in the right half of s-plane.
2) Varies ω from −∞ to ∞ and map all points on HG −plane. This can be done by
plotting polar plot for a given transfer function by varying ω from 0 to ∞. Then plot
for ω = −∞ to 0 is the mirror image of the polar plot with respect to real axis.
3) The mapping of infinite semicircle arc can be done by substituting s = Rejθ and
varying θ from πΤ2 to −πΤ2 and tending R to infinity.
4) If the open-loop transfer function has any poles on jω −axis, then map the segment
by letting s = ϵejθ , (i. e. ϵ → 0) and varying θ from −πΤ2 to πΤ2.
96
Cont. …
Example: For the system shown in figure 4.12(a) below, sketch the
Nyquist plot for k = 2 and use Nyquist criterion to determine whether
the closed-loop system is stable for this gain. Find the range of k for
the system to be stable.
k 1−s 2 1−s
H s G s = = , for k = 2
s+2 s+3 s+2 s+3
𝟐𝛚 𝛚𝟐 −𝟏𝟏
R e −intercept, Im G jω =0 ⇒ = 0 ⇒ ω = ωx = 11
𝛚𝟐 +𝟒 𝛚𝟐 +𝟗
99
Cont. …
➢ Mapping of section C2 : This section in s-plane is mapped to HG −plane by
jθ 2 1−Rejθ
G Re = G s │s=Rejθ =
Rejθ +2 Rejθ +3
jθ 2 Rejθ 1ΤRejθ −1
lim G Re = lim = lim ϵe−jθ
R→∞ R→∞ R2 ej2θ 1+2ΤRejθ 1+3ΤRejθ ∈→0
when θ = πΤ2, G Rejθ = lim ϵe−jπΤ2 and when θ = − πΤ2, G Rejθ = lim ϵejπΤ2
∈→0 ∈→0
➢ From the above value, section C2 in s-plane is mapped into a circular arc of radius
zero (0) at origin in the HG −plane.
100
Cont. …
➢ The complete Nyquist plot is
➢ The intersection of Nyquist plot with real axis is obtained by equating imaginary
part of H(jω)G(jω) to zero as solved above.
𝟐
𝟐𝛚 𝛚𝟐 −𝟏𝟏 𝟐 𝟏+ 𝟏𝟏
= 0 ⇒ ω = ωx = 11 radΤs, ∴ H(jω)G(jω) │ω=ωx = 11 = = 0.4
𝛚𝟐 +𝟒 𝛚𝟐 +𝟗 𝟐 𝟐
𝟒+ 𝟏𝟏 𝟗+ 𝟏𝟏
101
Cont. …
➢ Therefore the Nyquist plot cuts the real axis at 𝑅𝑒 𝐻𝐺 = −0.4 and the number of
encirclement of −1 + 𝑗0 point is zero. That is 𝑁 = 0, the number of open-loop poles in
the right of s-plane is zero. Z= 𝑁 + 𝑃 = 0 + 0 = 0
𝟐
𝒌 𝟏+ 𝟏𝟏
H(jω)G(jω) │ω= 11 = = k(0.2)
𝟐 𝟐
𝟒+ 𝟏𝟏 𝟗+ 𝟏𝟏
➢ For the system to be stable, the number of encirclements of Nyquist plot with −1 + 𝑗𝑜
point must be zero.
∴ k 0.2 < 1 ⇒ 𝑘 < 5, so that the range of k is, 0 < 𝑘 < 5. 102
Cont. …
Exercise 4.9:
5
H s G s =
s(s + 1)(s + 2)
Draw the Nyquist plot and determine whether the system is stable or
unstable.
103
End
104