Chapter-5 Int. To Control System 2023
Chapter-5 Int. To Control System 2023
(EPCE3204)
Credit: 3 Hrs (Lecture: 2 Hr, Tut: 3Hr & Lab: 3 Hr)
_______________________________________________________________
Chapter-Five
Outline: Frequency Domain Analysis
1. Introduction
• Methods of frequency response
• Correlation b/n Frequency domain & Time domain analysis
• Advantage & Disadvantage of frequency domain analysis
2. Bode plot
• Step to sketch bode plot
• Stability analysis using bode plot
3. polar plot
• Polar plot of standard function
4. Nyquist plot
• Nyquist analysis
• Nyquist stability criterion Course Instructor: Tamiru Getahun G. (Lecturer @ASTU, EPCE)
1. Objectives
After Completing this chapter, student will able to:
• Analyze the system behavior/ Dynamics.
• Analysis Steady state response of different systems.
• Plot the magnitude and phase curve of the system using frequency Domain
methods.
• Check the Stability of the system.
5.1 Introduction
• Frequency Response, we mean the steady-state response of a
system to a sinusoidal input.
Where :
5/29/2023 6
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 6
Frequency domain Analysis
• Then, the steady-state response becomes
• where
5/29/2023 7
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 7
Frequency domain Analysis
• Similarly, we obtain the following expression for G(–jω):
where
Then ,
5/29/2023 8
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 8
Frequency domain Analysis
5/29/2023 9
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 9
Frequency domain Analysis
• Generally,
• Example: For the system given below, the transfer function G(s) is
5/29/2023 10
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 10
Frequency domain Analysis
5/29/2023 11
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 11
Frequency domain Analysis
Cont…
• Basic Factors of G( jω)H(jω)
The basic factors that very frequently occur in an arbitrary
transfer function
5/29/2023 13
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 13
Frequency domain Analysis
The Gain K
• A number greater than unity has a positive value in
decibels, while a number smaller than unity has a
negative value.
• The log-magnitude curve for a constant gain K is a
horizontal straight line at the magnitude of 20 log K
decibels.
• The phase angle of the gain K is zero.
• The effect of varying the gain K in the transfer function
is that it raises or lowers the log-magnitude curve of the
transfer function by the corresponding constant amount,
but it has no effect on the phase curve.
5/29/2023 14
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 14
Frequency domain Analysis
Cont…
• A number–decibel conversion line
•As a number increases by a
factor of 10, the
corresponding decibel value
increases by a factor of 20.
This may be seen from the
following:
5/29/2023 15
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 15
Frequency domain Analysis
5/29/2023 16
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 16
Frequency domain Analysis
Cont…
5/29/2023 19
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 19
Frequency domain Analysis
Cont…
• At zero frequency, the phase angle is 0°.
• At the corner frequency, the phase angle is
w=1/T
5/29/2023 20
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 20
Frequency domain Analysis
• The Log-magnitude curve, together with the asymptotes,
and phase-angle curve for 1+jωT as follows.
5/29/2023 21
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 21
Frequency domain Analysis
• For high frequencies such that ω >> ωn, the log magnitude
becomes
• The phase angle of the quadratic factor [1+ 2ξ(jω/ ωn) + (jω/ ωn) 2 ]-1 is
5/29/2023 25
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 25
Frequency domain Analysis
The Resonant Frequency ωr and the Resonant Peak Value Mr
• The magnitude of
is a minimum.
5/29/2023 26
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 26
Frequency domain Analysis
• The above equation can be written
5/29/2023 27
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 27
Frequency domain Analysis
• For ξ ≥ 0.707, Mr = 1
• As ξ approaches zero, Mr approaches infinity.
• This means that if the undamped system is excited at its
natural frequency, the magnitude of G(jω) becomes
infinity.
• The relationship between Mr and ξ
• At the resonant frequency ωr ,
5/29/2023 28
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 28
Frequency domain Analysis
5/29/2023 29
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 29
Frequency domain Analysis
EXAMPLE: Draw the Bode diagram for the following transfer function
Solution:
1) rewrite the sinusoidal transfer function G(jω)H(jω)
as a product of basic factors
5/29/2023 30
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 30
Frequency domain Analysis
• Note that the last term has the damping ratio of 0.3536.
5/29/2023 31
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 31
Frequency domain Analysis
5/29/2023 32
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 32
Frequency domain Analysis
5/29/2023 33
ASTU – EPCE Dpt Introduction to Control Systems By Tamiru G. G.(Lecturer) 33
Frequency domain Analysis
• There are several ways of displaying frequency response data, including Bode' plots
and Nyquist plots.
• Bode' plots use frequency as the horizontal axis and use two separate plots to
display amplitude and phase of the frequency response.
• Nyquist plots display both amplitude and phase angle on a single plot, using
frequency as a parameter in the plot.
• Nyquist plots have properties that allow you to see whether a system is stable or
unstable.
Nyquist Plot
• A Nyquist plot is a polar plot of the frequency response function of a linear system.
• That means a Nyquist plot is a plot of the transfer function, G(s) with s = jω. That
means you want to plot G(j ω).
• G(j ω) is a complex number for any angular frequency, ω, so the plot is a plot of
complex numbers.
• The complex number, G(j ω), depends upon frequency, so frequency will be a
parameter if you plot the imaginary part of G(j ω) against the real part of G(j ω).
To sketch the polar plot of G(jω) for the entire range of frequency ω, i.e., from 0 to
infinity, there are four key points that usually need to be known:
1) The start of plot where ω = 0,
2) The end of plot where ω = ∞,
3) Where the plot crosses the real axis, i.e., Im(G(jω)) = 0, and
4) Where the plot crosses the imaginary axis, i.e., Re(G(jω)) = 0.
Magnitude
Phase
Magnitude
Phase
Polar Plot
Problem-4: Sketch the Nyquist diagram for the system shown in the following figure,
and then determine the system stability using the Nyquist criterion.
(a)
Now that we have expressions for the magnitude and phase of the frequency response, we
can sketch the polar plot using the 4 key points.
500 ∞ 3
𝐺 ∞ 𝐻(∞) = =0 ∠𝐺 ∞ 𝐻 ∞ = tan−1 = −3 × 90° = −270°
∞ 30
Point 3: Where the plot crosses the real axis, i.e., Im(G(jω)) = 0
Take the imaginary part of equation (a), and put equal to zero, to get the value of
frequency ω at the interception of real axis.
Take the real part of equation (a), and put equal to zero, to get the value of frequency
ω at the interception of imaginary axis.
Problem-5: Sketch the polar plot for the following transfer function.
Representing G(s) in the frequency response form G( jω ) by replacing s = jω, and then
Multiply both numerator and denominator by the conjugate of denominator.
(a)
At frequency ω = 0, we only observe the most significant terms that take the effect.
Magnitude at ω = 0:
Phase at ω = 0:
Magnitude at ω = ∞:
Phase at ω = ∞:
Point 3: Where the plot crosses the real axis, i.e., Im(G(jω)) = 0
Take the imaginary part of equation (a), and put equal to zero, to get the value of
frequency ω at the interception of real axis.
Therefore, the intersection point between the polar plot and the real axis, when 𝝎 = 𝟓,
is located at;
Take the real part of equation (a), and put equal to zero, to get the value of frequency
ω at the interception of imaginary axis.
Therefore, the intersection point between the polar plot and the imaginary axis is when
𝜔 = ∞ is located at;
Polar Plot
|G(jω)| ∠G(jω)
ω=0 ∞ -90o
ω=∞ 0 -270o
Cross Re: 1
−
𝜔= 5 3
Cross Img:
𝜔=∞ 0
1
𝐺𝐻 𝑗𝜔 =
𝜔 4 𝜔 2 + 𝑝2
The phase of GH(jω) denoted by, 𝜑, is obtained as;
𝜔
𝜑 = ∠𝐺𝐻 0 = − tan−1
𝑝
1
𝐺𝐻(∞) = =0
∞ ∞+ 𝑝2
It is conclude that;
1) The poles of 1 + G(s)H(s) are the same as the poles of G(s)H(s), the open-loop system.
2) The zeros of 1 + G(s)H(s) are the same as the poles of T(s), the closed-loop system.
• Equivalently, the poles of the closed-loop transfer function, or the zeros of the
denominator, 1 + GH(s), of the closed-loop transfer function, must lie in the left-
half plane.
• The Nyquist Stability Criterion establishes the number of zeros of 1 + GH(s) in the
right-half plane directly from the Nyquist Stability Plot of GH(s).
Where,
If N ≤ 0 and Po = 0, then the system is absolutely stable if and only if N = 0; that is, if
and only if the (-1,0) point does not lie in the shaded region.
Problem-7: Sketch the polar plot for the following transfer function and check
the stability of the system.
𝜔
𝜑 = ∠𝐺𝐻 0 = −90° − tan−1
1
Polar Plot
|GH(jω)| ∠GH(jω)
ω=0 ∞ -90o
ω=∞ 0 -180o
• Clearly, the (-1,0) point is not in the shaded region; therefore it is not enclosed by
the contour and so N ≤ 0.
• The poles of GH(s) are at s = 0 and s = -1, neither of which are in the right-hand-
plane RHP; hence Po = 0. Thus N = -Po = 0, and the system is absolutely stable.
Problem-8: The Nyquist Stability Plot for GH(s) = 1/s(s-1) is given in the figure below.
• The region to the right of the contour has been shaded and the (-1,0) point is enclosed; then
N > 0. (It is clear that N = 1).
• The poles of GH are at s = 0 and s = +1, the latter pole being in the RHP. Hence Po = 1.