CE Unit4
CE Unit4
RV College of world
Engineering
CONTROL ENGINEERING
(18AS62)
Department of Aerospace Engineering
R V College of Engineering
Bangalore
Go, change the
RV College of world
Engineering
UNIT-4
FREQUENCY RESPONSE
COURSE COORDINATOR
CONTENTS
1. Magnitude and phase,
2. Frequency response of higher order systems,
3. Bode Plot
4. Polar Plot
5. Nichols Plot,
6. Bandwidth,
7. Nyquist stability criterion,
8. Gain and phase margins.
9. Standard control actions: Proportional control, steady state error constants, system type, tracking control and
integral control, lag compensator, transient response improvement and derivative control, lead compensators.
Gain Margin
The greater the Gain Margin (GM), the greater the stability of the system. The
gain margin refers to the amount of gain, which can be increased or decreased
without making the system unstable. It is usually expressed as a magnitude in dB.
We can usually read the gain margin directly from the Bode plot (as shown in the
diagram above). This is done by calculating the vertical distance between the
magnitude curve (on the Bode magnitude plot) and the x-axis at the frequency
where the Bode phase plot = 180°. This point is known as the phase crossover
frequency.
It is important to realize that the Gain and the Gain Margin are not the same
things. In fact, the Gain Margin is the negative of the gain (in decibels, dB).
This will make sense when we look at the Gain margin formula.
Department of Aerospace Engineering, Bangalore
R V College of Engineering
Where G is the gain. This is the magnitude (in dB) as read from the vertical axis
of the magnitude plot at the phase crossover frequency.
In our example shown in the graph above, the Gain (G) is 20. Hence using our
formula for gain margin, the gain margin is equal to 0 – 20 dB = -20 dB
(unstable).
Phase Margin
The greater the Phase Margin (PM), the greater will be the stability of the system.
The phase margin refers to the amount of phase, which can be increased or
decreased without making the system unstable. It is usually expressed as a phase
in degrees.
We can usually read the phase margin directly from the Bode plot (as shown in
the diagram above). This is done by calculating the vertical distance between the
phase curve (on the Bode phase plot) and the x-axis at the frequency where the
Bode magnitude plot = 0 dB. This point is known as the gain crossover frequency.
It is important to realize that the phase lag and the Phase Margin are not the
same things. This will make sense when we look at the phase margin
formula.
Department of Aerospace Engineering, Bangalore
R V College of Engineering
Where is the phase lag (a number less than 0). This is the phase as read from
the vertical axis of the phase plot at the gain crossover frequency.
In our example shown in the graph above, the phase lag is -189°. Hence using our
formula for phase margin, the phase margin is equal to -189° – (-180°) = -9°
(unstable).
1. Gain Margin: Greater will the gain margin greater will be the stability of the
system. It refers to the amount of gain, which can be increased or decreased
without making the system unstable. It is usually expressed in dB.
2. Phase Margin: Greater will the phase margin greater will be the stability of
the system. It refers to the phase which can be increased or decreased without
making the system unstable. It is usually expressed in phase.
7. Factors: Every loop transfer function {i.e. G(s) × H(s)} product of various
factors like constant term K, Integral factors (jω), first-order factors
(1+jωT)(± n) where n is an integer, second-order or quadratic factors.
8. Slope: There is a slope corresponding to each factor and slope for each
factor is expressed in the dB per decade.
9. Angle: There is an angle corresponding to each factor and angle for each
factor is expressed in the degrees.
Now there are some results that one should remember in order to plot the Bode curve.
These results are written below:
1. Constant term K: This factor has a slope of zero dB per decade. There is no corner
frequency corresponding to this constant term. The phase angle associated with this
constant term is also zero.
2. Integral factor 1/(jω)n: This factor has a slope of -20 × n (where n is an integer)dB per
decade. There is no corner frequency corresponding to this integral factor. The phase
angle associated with this integral factor is -90 × n. Here n is also an integer.
3. First order factor 1/ (1+jωT): This factor has a slope of -20 dB per decade. The corner
frequency corresponding to this factor is 1/T radian per second. The phase angle
associated with this first factor is -tan– 1(ωT).
4. First order factor (1+jωT): This factor has a slope of 20 dB per decade. The corner
frequency corresponding to this factor is 1/T radian per second. The phase angle
associated with this first factor is tan– 1(ωT) .
5. Second order or quadratic factor : [{1/(1+(2ζ/ω)} × (jω) + {(1/ω2)} × (jω)2)]: This factor
has a slope of -40 dB per decade. The corner frequency corresponding to this factor is
ωn radian per second. The phase angle associated with this first factor is
Keeping all the above points in mind, we are able to draw a Bode plot for any kind of control system. Now let
us discuss the procedure of drawing a Bode plot:
Note: Bode plot is a frequency response plot, where the response depends on the
frequency. Hence, both magnitude and phase depends on frequency, as shown in the
generic plot above.
Important: Once the question is given, first check whether it is in the time constant form,
that exists in the form of
Step 1: Convert the given system in terms of (j) – Frequency domain. Thus, replace ‘s’ by (j)
Here, our focus is to find the Magnitude of the given
system. For this we have identify the Corner
frequencies.
We know that the frequency and time period are
inversely proportional. Thus we can identify the time
period from Time constant form and then we can arrive
in getting the corner frequencies
With the obtained corner frequencies and Amplitude, we have to represent the same in Semi-log plot.
X- Axis (Frequency)
Here all the ones should be replaced
with the powers of 10 (In semi log plot).
Since we have selected the first
frequency as 0.1, the plot will start from
the value 0.1 followed by 1,10,100, …..
Y- Axis (Magnitude)
Choose and indicate the appropriate
values in the graph as per the Amplitude
values obtained.
Since there are three terms and the first term corresponds to (1/s), then we have the phase angle for the
first term to be -90. and if ‘s’ is there in the numerator, then phase angle will be+90 degree.
Similarly, (1/s2) corresponds to -180 and if s2 is in the numerator, it will have the values +180 degree.
Therefore, Phase angle = = First term (1) +Second Term(2) + Third Term (3)
But the second term has the complex number (a+ib). Hence the phase angle for complex term can be
identified by tan-1(b/a). If the complex term is in the denominator, the corresponding phase angle
will have negative value and if it is there on the numerator, it is considered to have positive angle.
In order to have the phase angle plot more precise, addition to the corner frequencies considered,
few more frequency values are incorporated to get the smooth phase angle curve.
Hence, using the expression the corresponding values of phase angle is calculated and plotted.
Other side of the Y-axis in the Semi-log chart is used as Phase angle information. Also for magnitude
plot, scale is used to connect the points, whereas, free hand drawing is used for drawing phase plot.
In order to draw the graph, we should understand, that the graph will have starting and ending point:
Starting point is determined by the Type of system and the ending point is determined by the Order of the
system
Effect of adding pole to the system: Only the ending point will be affected.
If we are adding a Pole to the system, then our ending will shift by angle -90
If we are adding a Zero to the system, then our ending will shift by an angle +90
P
O
L
A
R
G
R
A
P
H
Increase the number of corner frequencies just below the upper limit
Polar Graph
Converted
Rectangular Graph
Use Calculator to
convert
Department of Aerospace Engineering, Bangalore
R V College of Engineering
From the given problem, it is clear that it is Type-1 system and having order value of 3 (If gives rise to s3 term)
The generic model graph for both Type and Order is shown below:
Therefore the model graph will start from Type-1 and end in order-3.
Nyquist Plot is used to determine the Stability of the closed loop system from the given Open
loop system
Objectives:
To Draw the Nyquist Plot and to determine the stability of the system
Example 1: Draw the Nyquist plot and determine the Stability of the system given below
For the given open loop system, the poles are identified from the denominator. i.e. s=-1 and +1.
Plotting this in the S-plane, shows that the system is unstable. However, we have to check for its stability
in the closed loop system.
Department of Aerospace Engineering, Bangalore
R V College of Engineering
For the Open loop system we mark the poles in the S-plane and identify whether the system is stable or
not. However to check for Closed loop system, the characteristic equation is obtained and to check the
stability of the closed loop system, it is recommended to have the plot in GH-plane.
Hence, there is a need to transform the S-plane to GH-plane, which is called as Mapping. For this
purpose, we will be using Nyquist Contour.
Department of Aerospace Engineering, Bangalore
R V College of Engineering
Since the given problem in open loop is unstable, we are going to consider the entire right side of the S-
plane, dividing it into 3 regions (C1, C2 and C3)
C1 = varied from 0 to infinity.
C2 = Right half Circle = Rejθ (R = Radius of the Circle and θ is the angle) – From the polar plot it moves from +90
to -90
C3 = varied - ∞ (Minus Infinity) to 0.
=2
=0
= -180
=0
=∞
Plotting Region 1 (C1) in GH plane
= -90
Note: In S-plane, the origin is 0, but while plotting the Nyquist plot, we can observe that the origin
will be (-1+j0)
Therefore Z = P – N
Z=1–1=0
If there are no Zeros on the RHS, then the Closed loop system is stable.
Though the Open Loop System is UNSTABLE, we can see that the Closed Loop System is STABLE.
Constant magnitude loci that are M-circles and constant phase angle loci that are N-circles are the
fundamental components in designing the Nichols chart.
The constant M and constant N circles in G (jω) plane can be used for the analysis and design of
control systems.
However the constant M and constant N circles in gain phase plane are prepared for system design
and analysis as these plots supply information with fewer manipulations.
Gain phase plane is the graph having gain in decibel along the ordinate (vertical axis) and phase
angle along the abscissa (horizontal axis).
The M and N circles of G (jω) in the gain phase plane are transformed into M and N contours in
rectangular co-ordinates.
A point on the constant M loci in G (jω) plane is transferred to gain phase plane by drawing the
vector directed from the origin of G (jω) plane to a particular point on M circle and then
measuring the length in db and angle in degree.
Department of Aerospace Engineering, Bangalore
R V College of Engineering
The critical point in G (jω), plane corresponds to the point of zero decibel and -180o in the gain
phase plane. Plot of M and N circles in gain phase plane is known as Nichols chart /plot.
The Nichols plot is named after the American engineer N.B Nichols who formulated this plot.
Compensators can be designed using Nichols plot. Nichols plot technique is however also used in
designing of dc motor. This is used in signal processing and control design.
Nyquist plot in complex plane shows how phase of transfer function and frequency variation of
magnitude are related.We can find out the gain and phase for a given frequency.
Angle of positive real axis determines the phase and distance from origin of complex plane
determines the gain. There are some advantages of Nichols plot in control system engineering. They
are:
Gain and phase margin can be determined easily and also graphically.
Closed loop frequency response is obtained from open loop frequency response.
1 • Introduction.
6 • Logarithmic coordinate.
Polar Plots
R V College of Engineering
Viewpoints of analyzing
control system behavior
• Routh-Hurwitz ( s j )
• Root locus
( s j )
• Bode diagram (plots)
• Nyquist plots ( s j )
• Nicols plots ( s j )
• Time domain
( s j )
R V College of Engineering
r (t ) A sin t y (t ) B sin( t )
L.T.I system
Magnitude: B Phase:
A
r (t ) y (t )
G(s)
+
- Steady state response
H(s)
Y ( s) G( s)
R( s ) 1 G ( s ) H ( s )
s j s j
Magnitude:
G ( j ) Phase: G( j )
1 G ( j ) H ( j ) [1 G( j ) H ( j )]
R V College of Engineering
2 2
Decade : dec log 10 Octave : oct log 2
1 1
dB
1 2 3 4 10 20 100
• The gain magnitude is many times expressed in terms of
decibels (dB)
dB = 20 log10 A
where A is the amplitude or gain
– a decade is defined as any 10-to-1 frequency range
– an octave is any 2-to-1 frequency range
R V College of Engineering
Y (s) k ( s z1 )( s z 2 )
R ( s ) ( s p1 )( s p2 )( s 2 as b)
Case I : k GH (dB)
Magnitude:
0.1 1 10
k dB 20 log k (dB)
GH
Phase: 180 0
0 900
o
,k 0
k o
180 ,k 0
R V College of Engineering
1
Case II :
sp
GH (dB)
Magnitude: p2
1 p 1
( j ) p
20 p log (dB )
dB 0.1 1 10
GH
Phase:
1 900
( j ) p
(90o ) p
p 1
900
p2
180 0
R V College of Engineering
Case III : p
s
p2
GH (dB) p 1
Magnitude:
( j ) p
dB
20 p log (dB)
0.1 1 10
GH
Phase: 0
p2
180
0
p 1
90
( j ) (90 ) p
p o
900
180 0
a 1
CaseofIVEngineering
R V College : or ( s 1) 1
( s a) a
a 1
Magnitude:
GH (dB)
(1 j ) 1 20 log 1 ( ) 2
a dB a
10 log[ 1 ( ) 2 ]
a 0.1 1 10
a 0 dB 10 log 1 0
a
a 1 j dB 20 log GH
a a a a 450
dB [20 log 20 log a] 180 0
a 1 j1 dB 10 log 2 3.01 900
Phase:
1 900
(1 j ) 0 tan
0
a a 180 0
a 0 GH tan 1 0 0o
a
a GH tan 1 90o
a
( s a)
Case Vof: Engineering
R V College 1
or ( s 1)
a a a 1
Magnitude:
GH (dB)
(1 j ) 20 log 1 ( ) 2
a dB a
10 log[ 1 ( ) 2 ]
a 0.1 1 10
a 0 dB 10 log 1 0
a
a 1 j dB 20 log GH
a a a a 450
dB 20 log 20 log a 180 0
a 1 j1 dB 10 log 2 3.01 900
Phase:
1 900
(1 j ) tan
a a 180 0
a 0 GH tan 1 0 0o
a
a GH tan 1 90o
a
R V College of Engineering
2
CaseofVIEngineering
R V College : T ( s) 2 n
s 2n s n2
n2 2n
T ( j ) T ( j ) tan 1
( n 2 ) 2 jn ( n 2 )
2 2
2
1 n
T ( j ) T ( j ) tan 1
2 2
(1 ( ) ) j 2 1 ( )
n n n
1 , 1
0 ,
n
n 0 0
T ( j ) 20 log( 2 ) , 1 T ( j ) 90 , 1
0
n 180 o n
40 log( ) , 1 , 1
n n n
R V College of Engineering
n
R V College of Engineering
Example : 50(s 2)
T (s )
s (s 10)
1 s 2 10
T ( s ) 10( )( )( )
s 2 s 10
R V College of Engineering
R V College of Engineering k ( s z1 )
Minimum phase system T ( s) n , zi 0, pi 0
s ( s p1 )
Type 0 : (i.e. n=0)
0dB/dec
GH (dB)
k p p1
T (s)
( s p1 ) A
20 log K p A 0.1 p1 p1 10 p1
R V College of Engineering
Type I : (i.e. n=1)
k v p1 -20dB/dec
T (s) GH (dB)
s ( s p1 )
A
20 log K v A 0
10 p1
0.1 p1 1 p1
Kv
20 log 0dB -40dB/dec
j 0
0 k v
R V College of Engineering
Type 2 : (i.e. n=2)
-40dB/dec
k a p1
T ( s) 2 GH (dB)
s ( s p1 )
A
20 log K a A
0
0.1 p1 1 p1 10 p1
Ka
20 log 0dB -60dB/dec
( j 0 ) 2
0 ka
2
R V College of Engineering
Relative Stability
A transfer function is called minimum phase when all the poles and zeroes are
LHP and non-minimum-phase when there are RHP poles or zeroes.
The gain margin (GM) is the distance on the bode magnitude plot from
the amplitude at the phase crossover frequency up to the 0 dB point.
GM=-(dB of GH measured at the phase crossover frequency)
The phase margin (PM) is the distance from -180 up to the phase at the
gain crossover frequency. PM=180+phase of GH measured at the gain
crossover frequency
R V College of Engineering
Open loop transfer function : G ( s) H ( s)
Closed-loop transfer function : 1 G (s) H (s)
GH
180 0 Stable system
900
900
p
180 0
P.M.<0
Unstable system
R V College of Engineering
R V College of Engineering
H ( j )
K j (1 j 1 ) 1 2 2 ( j 2 ) ( j 2 ) 2
N
(1 j a ) 1 2 b ( j b ) ( j b ) 2
( j ) N
0 dB
Zero : K ( j ) N 0 dB (1 / K )1 / N
0 dB +20 dB
–20 dB 0 dB
ω ω
Phase One Phase One
Decade Decade
0° +90°
–45° +45°
–90° 0°
ω ω
• We can also take the Bode plot and extract the transfer
function from it (although in reality there will be error
associated with our extracting information from the
graph)
• First, determine the constant gain factor, K
• Next, move from lowest to highest frequency noting the
appearance and order of the poles and zeros
R V College of Engineering
Frequency Response Plots
where K is a variable (constant) gain and G(s) is the plant under consideration.
The gain margin is defined as the change in open loop gain required to
make the system unstable. Systems with greater gain margins can withstand
greater changes in system parameters before becoming unstable in closed
loop. Keep in mind that unity gain in magnitude is equal to a gain of zero in
dB
The phase margin is defined as the change in open loop phase shift required
to make a closed loop system unstable.
The phase margin is the difference in phase between the phase curve and -
180 deg at the point corresponding to the frequency that gives us a gain of
0dB (the gain cross over frequency, Wgc).
Likewise, the gain margin is the difference between the magnitude curve
and 0dB at the point corresponding to the frequency that gives us a phase of
-180 deg (the phase cross over frequency, Wpc).
R V College of Engineering
Gain and Phase Margin
-180
R V College of Engineering
Examples - Bode
R V College of Engineering
Examples - Bode
R V College of Engineering
Examples – Bode