Control Systems Notes
Control Systems Notes
IV Semester: ECE
2 2 - 3 30 70 100
OBJECTIVES:
Control systems: Introduction, open loop and closed loop systems, examples, mathematical models and
differential equations of physical systems, concept of transfer function, translational and rotational
mechanical systems, electrical systems, force voltage and force current analogy. Block Diagrams: Block
diagram representation of various systems, block diagram algebra, characteristics of feedback systems
DC servomotors, signal flow graph, Mason’s gain formula.
Time response analysis: Standard test signals, shifted unit step, ramp and impulse signals, shifting
theorem, convolution integral, impulse response, unit step response of first and second order systems,
time domain specifications, steady state errors and error constants effects of proportional, derivative and
proportional derivative, proportional integral and PID controllers.
Concept of stability: Necessary and sufficient conditions for stability, Routh’s and Routh Hurwitz stability
criterions and limitations.
Root locus technique: Introduction, root locus concept, construction of root loci, graphical determination
of ‘k’ for specified damping ratio, relative stability, effect of adding zeros and poles on stability.
Frequency domain analysis: Introduction, frequency domain specifications, stability analysis from Bode
plot, polar plot, Nyquist plot, calculation of gain margin and phase margin, determination of transfer
function from bode plot, correlation between time and frequency responses. Compensators: Lag, lead, lag
lead networks.
Text Books:
Reference Books:
Web References:
1. https://www.researchgate.net
2. https://www.aar.faculty.asu.edu/classes
3. https://www.facstaff.bucknell.edu/
4. https://www.electrical4u.com
5. https://www.crectirupati.com
E-Text Books:
1. https://www.jntubook.com/
2. https://www.freeengineeringbooks.com
Course Outcome:
UNIT I
Control Systems -Introduction
A control system is a system, which provides the desired response by controlling the output.
The following figure shows the simple block diagram of a control system.
Here, the control system is represented by a single block. Since, the output is controlled by
varying input, the control system got this name. We will vary this input with some mechanism.
In the next section on open loop and closed loop control systems, we will study in detail about
the blocks inside the control system and how to vary this input in order to get the desired
response.
Examples − Traffic lights control system, washing machine
Traffic lights control system is an example of control system. Here, a sequence of input signal
is applied to this control system and the output is one of the three lights that will be on for some
duration of time. During this time, the other two lights will be off. Based on the traffic study at a
particular junction, the on and off times of the lights can be determined. Accordingly, the input
signal controls the output. So, the traffic lights control system operates on time basis.
∙ In continuous time control systems, all the signals are continuous in time. But, in discrete time control systems, there
exists one or more discrete time signals.
∙ SISO (Single Input and Single Output) control systems have one input and one output. Whereas, MIMO (Multiple
Inputs and Multiple Outputs) control systems have more than one input and more than one output.
The error detector produces an error signal, which is the difference between the input and the
feedback signal. This feedback signal is obtained from the block (feedback elements) by
considering the output of the overall system as an input to this block. Instead of the direct input,
the error signal is applied as an input to a controller.
So, the controller produces an actuating signal which controls the plant. In this combination, the
output of the control system is adjusted automatically till we get the desired response. Hence,
the closed loop control systems are also called the automatic control systems. Traffic lights
control system having sensor at the input is an example of a closed loop control system.
The differences between the open loop and the closed loop control systems are mentioned in the
following table.
Open Loop Control Systems Closed Loop Control Systems
These are also called as non-feedback These are also called as feedback
control systems. control systems.
Inaccurate. Accurate.
Control Systems - Feedback
If either the output or some part of the output is returned to the input side and utilized as part of
the system input, then it is known as feedback. Feedback plays an important role in order to
improve the performance of the control systems. In this chapter, let us discuss the types of
feedback & effects of feedback.
Types of Feedback
There are two types of feedback −
∙ Positive feedback
∙ Negative feedback
Positive Feedback
The positive feedback adds the reference input, R(s) and feedback output. The following figure
shows the block diagram of positive feedback control system.
The concept of transfer function will be discussed in later chapters. For the time being, consider
the transfer function of positive feedback control system is,
Where,
∙ T is the transfer function or overall gain of positive feedback control system.
Negative Feedback
Negative feedback reduces the error between the reference input, R(s) and system output. The
following figure shows the block diagram of the negative feedback control system.
Where,
∙ T is the transfer function or overall gain of negative feedback control system.
Effects of Feedback
Let us now understand the effects of feedback.
∙ If the value of (1+GH) is less than 1, then the overall gain increases. In this case, 'GH' value is negative because the gain
of the feedback path is negative.
∙ If the value of (1+GH) is greater than 1, then the overall gain decreases. In this case, 'GH' value is positive because the
gain of the feedback path is positive.
In general, 'G' and 'H' are functions of frequency. So, the feedback will increase the overall gain
of the system in one frequency range and decrease in the other frequency range.
So, we got the sensitivity of the overall gain of closed loop control system as the reciprocal of
(1+GH). So, Sensitivity may increase or decrease depending on the value of (1+GH).
∙ If the value of (1+GH) is less than 1, then sensitivity increases. In this case, 'GH' value is negative because the gain of
feedback path is negative.
∙ If the value of (1+GH) is greater than 1, then sensitivity decreases. In this case, 'GH' value is positive because the gain
of feedback path is positive.
In general, 'G' and 'H' are functions of frequency. So, feedback will increase the sensitivity of
the system gain in one frequency range and decrease in the other frequency range. Therefore,
we have to choose the values of 'GH' in such a way that the system is insensitive or less
sensitive to parameter variations.
∙ In Equation 2, if the denominator value is zero (i.e., GH = -1), then the output of the control system will be infinite. So,
the control system becomes unstable.
Therefore, we have to properly choose the feedback in order to make the control system
stable. Effect of Feedback on Noise
To know the effect of feedback on noise, let us compare the transfer function relations with and
without feedback due to noise signal alone.
Consider an open loop control system with noise signal as shown below.
∙ Get the differential equation in terms of input and output by eliminating the intermediate variable(s).
Example
Consider the following electrical system as shown in the following figure. This circuit consists
of resistor, inductor and capacitor. All these electrical elements are connected in series. The
input voltage applied to this circuit is vi and the voltage across the capacitor is the output voltage
vo.
Here, we represented an LTI system with a block having transfer function inside it. And this
block has an input X(s) & output Y(s).
Mass
Mass is the property of a body, which stores kinetic energy. If a force is applied on a body
having mass M, then it is opposed by an opposing force due to mass. This opposing force is
proportional to the acceleration of the body. Assume elasticity and friction are negligible.
Fm∝a
⇒Fm=Ma=Md2x/dt2
F=Fm=Md2/xdt2
Where,
∙ F is the applied force
∙ M is mass
∙ a is acceleration
∙ x is displacement
Spring
Spring is an element, which stores potential energy. If a force is applied on spring K, then it is
opposed by an opposing force due to elasticity of spring. This opposing force is proportional to
the displacement of the spring. Assume mass and friction are negligible.
F∝x
⇒Fk=Kx
F=Fk=Kx
Where,
∙ F is the applied force
∙ K is spring constant
∙ x is displacement
Dashpot
If a force is applied on dashpot B, then it is opposed by an opposing force due to friction of the
dashpot. This opposing force is proportional to the velocity of the body. Assume mass and
elasticity are negligible.
Fb∝ν
⇒Fb=Bν=Bdx/dt
F=Fb=Bdx/dt
Where,
∙ Fb is the opposing force due to friction of dashpot
∙ v is velocity
∙ x is displacement
Moment of Inertia
In translational mechanical system, mass stores kinetic energy. Similarly, in rotational
mechanical system, moment of inertia stores kinetic energy.
If a torque is applied on a body having moment of inertia J, then it is opposed by an opposing
torque due to the moment of inertia. This opposing torque is proportional to angular
acceleration of the body. Assume elasticity and friction are negligible.
T=Tj=Jd2/θdt2
Where,
∙ T is the applied torque
Tj∝α
⇒Tj=Jα=Jd2/θdt2
∙ Tj is the opposing torque due to moment of inertia
∙ J is moment of inertia
∙ α is angular acceleration
∙ θ is angular displacement
Torsional Spring
In translational mechanical system, spring stores potential energy. Similarly, in rotational
mechanical system, torsional spring stores potential energy.
If a torque is applied on torsional spring K, then it is opposed by an opposing torque due to the
elasticity of torsional spring. This opposing torque is proportional to the angular displacement
of the torsional spring. Assume that the moment of inertia and friction are negligible.
Tk∝θ
⇒Tk=Kθ
T=Tk=Kθ
Where,
∙ T is the applied torque
Dashpot
If a torque is applied on dashpot B, then it is opposed by an opposing torque due to the
rotational friction of the dashpot. This opposing torque is proportional to the angular velocity
of the body. Assume the moment of inertia and elasticity are negligible.
Tb∝ω
⇒Tb=Bω=Bdθdt
T=Tb=Bdθdt
Where,
∙ Tb is the opposing torque due to the rotational friction of the dashpot
()
⇒V=Ld2qdt2+Rdqdt+ 1c q (Equation 3)
By comparing Equation 1 and Equation 3, we will get the analogous quantities of the
translational mechanical system and electrical system. The following table shows these
analogous quantities.
Translational Mechanical System Electrical System
Force(F) Voltage(V)
Mass(M) Inductance(L)
Displacement(x) Charge(q)
Velocity(v) Current(i)
Similarly, there is torque voltage analogy for rotational mechanical systems. Let us now discuss
about this analogy.
Torque(T) Voltage(V)
Force(F) Current(i)
Mass(M) Capacitance(C)
Velocity(v) Voltage(V)
Similarly, there is a torque current analogy for rotational mechanical systems. Let us now
discuss this analogy.
Torque(T) Current(i)
Moment of inertia(J) Capacitance(C)
In this chapter, we discussed the electrical analogies of the mechanical systems. These analogies
are helpful to study and analyze the non-electrical system like mechanical system from
analogous electrical system.
The above block diagram consists of two blocks having transfer functions G(s) and H(s). It is
also having one summing point and one take-off point. Arrows indicate the direction of the flow
of signals. Let us now discuss these elements one by one.
Block
The transfer function of a component is represented by a block. Block has single input and
single output.
The following figure shows a block having input X(s), output Y(s) and the transfer function
G(s).
Transfer Function,G(s)=Y(s)X(s)
⇒Y(s)=G(s)X(s)
Output of the block is obtained by multiplying transfer function of the block with
The following figure shows the summing point with two inputs (A, B) and one output (Y). Here,
the inputs A and B are having opposite signs, i.e., A is having positive sign and B is having
negative sign. So, the summing point produces the output Y as the difference of A and B.
Y = A + (-B) = A - B.
The following figure shows the summing point with three inputs (A, B, C) and one output (Y).
Here, the inputs A and B are having positive signs and C is having a negative sign. So, the
summing point produces the output Y as
Y = A + B + (−C) = A + B − C.
Take-off Point
The take-off point is a point from which the same input signal can be passed through more than
one branch. That means with the help of take-off point, we can apply the same input to one or
more blocks, summing points.
In the following figure, the take-off point is used to connect the same input, R(s) to two more
blocks.
In the following figure, the take-off point is used to connect the output C(s), as one of the inputs
to the summing point.
Series Connection
Series connection is also called cascade connection. In the following figure, two blocks having
transfer functions G1(s) and G2(s) are connected in series.
Similarly, you can represent series connection of ‘n’ blocks with a single block. The transfer
function of this single block is the product of the transfer functions of all those ‘n’ blocks.
Parallel Connection
The blocks which are connected in parallel will have the same input. In the following figure,
two blocks having transfer functions G1(s) and G2(s) are connected in parallel. The outputs of
these two blocks are connected to the summing point.
Similarly, you can represent parallel connection of ‘n’ blocks with a single block. The transfer
function of this single block is the algebraic sum of the transfer functions of all those ‘n’ blocks.
Feedback Connection
As we discussed in previous chapters, there are two types of feedback — positive feedback and
negative feedback. The following figure shows negative feedback control system. Here, two
blocks having transfer functions G(s) and H(s) form a closed loop.
Similarly, you can represent the positive feedback connection of two blocks with a single block.
The transfer function of this single block is the closed loop transfer function of the positive
feedback, i.e., G(s)1−G(s)H(s)
Summing point has two inputs R(s) and X(s). The output of it is {R(s)+X(s)}.
So, the input to the block G(s) is {R(s)+X(s)} and the output of it is –
Y(s)=G(s){R(s)+X(s)}
⇒Y(s)=G(s)R(s)+G(s)X(s) (Equation 1)
Now, shift the summing point after the block. This block diagram is shown in the following
figure.
Output of the block G(s) is G(s)R(s).
The output of the summing point is
Y(s)=G(s)R(s)+X(s) (Equation 2)
Compare Equation 1 and Equation 2.
The first term ‘G(s)R(s)′ is same in both the equations. But, there is difference in the second
term. In order to get the second term also same, we require one more block G(s). It is having the
input X(s) and the output of this block is given as input to summing point instead of X(s). This
block diagram is shown in the following figure.
Here, X(s)=Y(s)=G(s)R(s)
When you shift the take-off point before the block, the output Y(s) will be same. But, there is
difference in X(s) value. So, in order to get same X(s) value, we require one more block G(s). It
is having the input R(s) and the output is X(s). This block diagram is shown in the following
figure.
∙ Rule 3 − Check for the blocks connected in feedback loop and simplify.
∙ Rule 4 − If there is difficulty with take-off point while simplifying, shift it towards right.
∙ Rule 5 − If there is difficulty with summing point while simplifying, shift it towards left.
∙ Rule 6 − Repeat the above steps till you get the simplified form, i.e., single block.
Note − The transfer function present in this single block is the transfer function of the overall
block diagram.
Example
Consider the block diagram shown in the following figure. Let us simplify (reduce) this block
diagram using the block diagram reduction rules.
Step 1 − Use Rule 1 for blocks G1 and G2. Use Rule 2 for blocks G3 and G4. The modified
block diagram is shown in the following figure.
Step 2 − Use Rule 3 for blocks G1G2 and H1. Use Rule 4 for shifting take-off point after the
block G5. The modified block diagram is shown in the following figure.
Step 3 − Use Rule 1 for blocks (G3+G4) and G5. The modified block diagram is shown in the
following figure.
Step 4 − Use Rule 3 for blocks (G3+G4)G5 and H3. The modified block diagram is shown in the
following figure.
Step 5 − Use Rule 1 for blocks connected in series. The modified block diagram is shown in the
following figure.
Step 6 − Use Rule 3 for blocks connected in feedback loop. The modified block diagram is
shown in the following figure. This is the simplified block diagram.
Therefore, the transfer function of the system is
Y(s)R(s)=G1G2G25(G3+G4)(1+G1G2H1){1+(G3+G4)G5H3}G5−G1G2G5(G3+G4)H2 Note −
Follow these steps in order to calculate the transfer function of the block diagram having
multiple inputs.
∙ Step 1 − Find the transfer function of block diagram by considering one input at a time and make the remaining inputs
as zero.
∙ Step 3 − Get the overall transfer function by adding all those transfer functions.
The block diagram reduction process takes more time for complicated systems. Because, we
have to draw the (partially simplified) block diagram after each step. So, to overcome this
drawback, use signal flow graphs (representation).
In the next two chapters, we will discuss about the concepts related to signal flow graphs, i.e.,
how to represent signal flow graph from a given block diagram and calculation of transfer
function just by using a gain formula without doing any reduction process.
Node
Node is a point which represents either a variable or a signal. There are three types of nodes —
input node, output node and mixed node.
∙ Input Node − It is a node, which has only outgoing branches.
∙ Mixed Node − It is a node, which has both incoming and outgoing branches.
Example
Let us consider the following signal flow graph to identify these nodes.
∙ The nodes present in this signal flow graph are y1, y2, y3 and y4.
∙ y1 and y4 are the input node and output node respectively.
Branch
Branch is a line segment which joins two nodes. It has both gain and direction. For example,
there are four branches in the above signal flow graph. These branches have gains of a, b, c and
-d.
Step 2 − Signal flow graph for y3=a23y2+a53y5 is shown in the following figure.
Step 3 − Signal flow graph for y4=a34y3 is shown in the following figure.
Step 4 − Signal flow graph for y5=a45y4+a35y3 is shown in the following figure.
Step 5 − Signal flow graph for y6=a56y5 is shown in the following figure.
Step 6 − Signal flow graph of overall system is shown in the following figure.
∙ Represent the transfer functions inside the blocks of block diagram as gains of the branches in signal flow graph.
∙ Connect the nodes as per the block diagram. If there is connection between two nodes (but there is no block in
between), then represent the gain of the branch as one. For example, between summing points, between summing
point and takeoff point, between input and summing point, between take-off point and output.
Example
Let us convert the following block diagram into its equivalent signal flow graph.
Represent the input signal R(s) and output signal C(s) of block diagram as input node R(s) and
output node C(s) of signal flow graph.
Just for reference, the remaining nodes (y1 to y9) are labelled in the block diagram. There are
nine nodes other than input and output nodes. That is four nodes for four summing points, four
nodes for four take-off points and one node for the variable between blocks G1 and G2.
The following figure shows the equivalent signal flow graph.
With the help of Mason’s gain formula (discussed in the next chapter), you can calculate the
transfer function of this signal flow graph. This is the advantage of signal flow graphs. Here, we
no need to simplify (reduce) the signal flow graphs for calculating the transfer function.
Δ=1−(sumofallindividualloopgains)
+(sumofgainproductsofallpossibletwonontouchingloops)
−(sumofgainproductsofallpossiblethreenontouchingloops)+...
Δi is obtained from Δ by removing the loops which are touching the ith forward path.
Consider the following signal flow graph in order to understand the basic terminology involved
here.
Path
It is a traversal of branches from one node to any other node in the direction of branch arrows. It
should not traverse any node more than once.
Examples − y2→y3→y4→y5 and y5→y3→y2
Forward Path
The path that exists from the input node to the output node is known as forward path.
Examples − y1→y2→y3→y4→y5→y6 and y1→y2→y3→y5→y6.
Forward Path Gain
It is obtained by calculating the product of all branch gains of the forward path.
Examples − abcde is the forward path gain of y1→y2→y3→y4→y5→y6 and abge is the forward
path gain of y1→y2→y3→y5→y6.
Loop
The path that starts from one node and ends at the same node is known as loop. Hence, it is a
closed path.
Examples − y2→y3→y2 and y3→y5→y3.
Loop Gain
It is obtained by calculating the product of all branch gains of a loop.
Examples − bj is the loop gain of y2→y3→y2 and gh is the loop gain of y3→y5→y3.
Non-touching Loops
These are the loops, which should not have any common node.
Examples − The loops, y2→y3→y2 and y4→y5→y4 are non-touching.
Calculation of Transfer Function using Mason’s Gain Formula
Let us consider the same signal flow graph for finding transfer function.
∙ Number of forward paths, N = 2.
Δ=1−(sumofallindividualloopgains)
+(sumofgainproductsofallpossibletwonontouchingloops)
−(sumofgainproductsofallpossiblethreenontouchingloops)+...
Substitute the values in the above equation,
Δ=1−(bj+gh+cdh+di+f)+(bjdi+bjf)−(0)
⇒Δ=1−(bj+gh+cdh+di+f)+bjdi+bjf
There is no loop which is non-touching to the first forward path.
So, Δ1=1.
Similarly, Δ2=1. Since, no loop which is non-touching to the second forward path.
Substitute, N = 2 in Mason’s gain formula
T=C(s)R(s)=Σ2i=1PiΔiΔ
T=C(s)R(s)=P1Δ1+P2Δ2Δ
Substitute all the necessary values in the above equation.
T=C(s)R(s)=(abcde)1+(abge)11−(bj+gh+cdh+di+f)+bjdi+bjf
⇒T=C(s)R(s)=(abcde)+(abge)1−(bj+gh+cdh+di+f)+bjdi+bjf
Therefore, the transfer function is -
T=C(s)R(s)=(abcde)+(abge)1−(bj+gh+cdh+di+f)+bjdi+bjf
UNIT 2
∙ Transient response
∙ Steady state response
The response of control system in time domain is shown in the following figure.
Here, both the transient and the steady states are indicated in the figure. The responses corresponding to these states are known
as transient and steady state responses.
c(t)=ctr(t)+css(t)c(t)=ctr(t)+css(t)
Where,
Transient Response
After applying input to the control system, output takes certain time to reach steady state. So, the output will be in transient
state till it goes to a steady state. Therefore, the response of the control system during the transient state is known as transient
response.
The transient response will be zero for large values of ‘t’. Ideally, this value of ‘t’ is infinity and practically, it is five times
constant.
limt→∞ctr(t)=0limt→∞ctr(t)=0
So, the unit impulse signal exists only at ‘t’ is equal to zero. The area of this signal under small interval of time around ‘t’ is
equal to zero is one. The value of unit impulse signal is zero for all other values of ‘t’.
u(t)=1;t≥0u(t)=1;t≥0
=0;t<0=0;t<0
Following figure shows unit step signal.
So, the unit step signal exists for all positive values of ‘t’ including zero. And its value is one during this interval. The value of
the unit step signal is zero for all negative values of ‘t’.
r(t)=t;t≥0r(t)=t;t≥0
=0;t<0=0;t<0
We can write unit ramp signal, r(t)r(t) in terms of unit step signal, u(t)u(t) as
r(t)=tu(t)r(t)=tu(t)
p(t)=t22;t≥0p(t)=t22;t≥0
=0;t<0=0;t<0
We can write unit parabolic signal, p(t)p(t) in terms of the unit step signal, u(t)u(t)as,
p(t)=t22u(t)p(t)=t22u(t)
So, the unit parabolic signal exists for all the positive values of ‘t’ including zero. And its value increases non-linearly with
respect to ‘t’ during this interval. The value of the unit parabolic signal is zero for all the negative values of ‘t’.
We know that the transfer function of the closed loop control system has unity negative feedback
as, C(s)R(s)=G(s)1+G(s)C(s)R(s)=G(s)1+G(s)
The power of s is one in the denominator term. Hence, the above transfer function is of the first order and the system is said to
be the first order system.
C(s)=(1sT+1)R(s)C(s)=(1sT+1)R(s)
Where,
Follow these steps to get the response (output) of the first order system in the time domain.
signal. For example, the response of the system for an impulse input is called as impulse response. Step Response
of First Order System
Consider the unit step signal as an input to first order system.
So, r(t)=u(t)r(t)=u(t)
Apply Laplace transform on both the sides.
R(s)=1sR(s)=1s
C(s)=1s(sT+1)=As+BsT+1C(s)=1s(sT+1)=As+BsT+1
⇒1s(sT+1)=A(sT+1)+Bss(sT+1)⇒1s(sT+1)=A(sT+1)+Bss(sT+1)
On both the sides, the denominator term is the same. So, they will get cancelled by each other. Hence, equate the numerator
terms.
1=A(sT+1)+Bs1=A(sT+1)+Bs
By equating the constant terms on both the sides, you will get A = 1.
Substitute, A = 1 and equate the coefficient of the s terms on both the sides.
0=T+B⇒B=−T0=T+B⇒B=−T
⇒C(s)=1s−1s+1T⇒C(s)=1s−1s+1T
The unit step response, c(t) has both the transient and the steady state terms.
ctr(t)=−e−(tT)u(t)ctr(t)=−e−(tT)u(t)
css(t)=u(t)css(t)=u(t)
The following figure shows the unit step response.
The value of the unit step response, c(t) is zero at t = 0 and for all negative values of t. It is gradually increasing from zero
value and finally reaches to one in steady state. So, the steady state value depends on the magnitude of the input.
We know that the transfer function of the closed loop control system having unity negative feedback
as
The power of ‘s’ is two in the denominator term. Hence, the above transfer function is of the second order and the system is
said to be the second order system.
All the time domain specifications are represented in this figure. The response up to the settling time is known as transient
response and the response after the settling time is known as steady state response.
Delay Time
It is the time required for the response to reach half of its final value from the zero instant. It is denoted by
tdtd. Consider the step response of the second order system for t ≥ 0, when ‘δ’ lies between zero and one.
Therefore, at t=tdt=td, the value of the step response will be 0.5. Substitute, these values in the above equation.
By using linear approximation, you will get the delay time td as
Rise Time
It is the time required for the response to rise from 0% to 100% of its final value. This is applicable for the under-damped
systems. For the over-damped systems, consider the duration from 10% to 90% of the final value. Rise time is denoted by tr.
From above equation, we can conclude that the rise time trtr and the damped frequency ωdωd are inversely proportional to
each other.
Peak Time
It is the time required for the response to reach the peak value for the first time. It is denoted by tptp. At t=tpt=tp, the first
derivate of the response is zero.
We know the step response of second order system for under-damped case is
From the above equation, we can conclude that the peak time tptp and the damped frequency ωdωd are inversely proportional
to each other.
Peak Overshoot
Peak overshoot Mp is defined as the deviation of the response at peak time from the final value of response. It is also called the
maximum overshoot.
Mp=c(tp)−c(∞)Mp=c(tp)−c(∞)
Where,
Where,
ess=ess1+ess2+ess3ess=ess1+ess2+ess3
⇒ess=0+0+1=1⇒ess=0+0+1=1
UNIT 3
Stability
Stability is an important concept. In this chapter, let us discuss the stability of system and types of systems based on
This is the response of first order control system for unit step input. This response has the values between 0 and 1. So, it is
bounded output. We know that the unit step signal has the value of one for all positive values of tincluding zero. So, it is
bounded input. Therefore, the first order control system is stable since both the input and the output are bounded. Types
of Systems based on Stability
We can classify the systems based on stability as follows.
a0sn+a1sn−1+a2sn−2+...+an−1s1+ans0=0a0sn+a1sn−1+a2sn−2+...+an−1s1+ans0=0
Note that, there should not be any term missing in the nth order characteristic equation. This means that the nth order
characteristic equation should not have any coefficient that is of zero value.
∙ Fill the first two rows of the Routh array with the coefficients of the characteristic polynomial as mentioned in the table
below. Start with the coefficient of snsn and continue up to the coefficient of s0s0.
∙ Fill the remaining rows of the Routh array with the elements as mentioned in the table below. Continue this process till
you get the first column element of row s0s0 is anan. Here, anan is the coefficient of s0s0in the characteristic
polynomial.
Note − If any row elements of the Routh table have some common factor, then you can divide the row elements with that factor
for the simplification will be easy.
The following table shows the Routh array of the nth order characteristic polynomial.
a0sn+a1sn−1+a2sn−2+...+an−1s1+ans0a0sn+a1sn−1+a2sn−2+...+an−1s1+ans0
snsn a0a0 a2a2 a4a4 a6a .. ..
6 . .
⋮⋮ ⋮⋮ ⋮⋮ ⋮⋮
s1s1 ⋮⋮ ⋮⋮
s0s0 anan
Example
Let us find the stability of the control system having characteristic equation,
s4+3s3+3s2+2s+1=0s4+3s3+3s2+2s+1=0
All the coefficients of the characteristic polynomial, s4+3s3+3s2+2s+1s4+3s3+3s2+2s+1are positive. So, the control
system satisfies the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
s4s4 11 33 11
s3s3 33 22
s1s1 (73×2)−(1×3)73=57(73×2)−(1×3)73=57
s0s0 11
All the elements of the first column of the Routh array are positive. There is no sign change in the first column of the Routh
array. So, the control system is stable.
Let us find the stability of the control system having characteristic equation,
s4+2s3+s2+2s+1=0s4+2s3+s2+2s+1=0
All the coefficients of the characteristic polynomial, s4+2s3+s2+2s+1s4+2s3+s2+2s+1are positive. So, the control system
satisfied the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
s4s4 11 11 11
s3s3 21 21
s1s1
s0s0
The row s3s3 elements have 2 as the common factor. So, all these elements are divided by 2.
Special case (i) − Only the first element of row s2s2 is zero. So, replace it by ϵϵand continue the process of completing the
Routh table.
s4s4 1 1 1
s3s3 1 1
s2s2 ϵϵ 1
s1s1 (ϵ×1)−(1×1)ϵ=ϵ−1ϵ(ϵ×1)−(1×1)ϵ=ϵ−1ϵ
s0s0 1
s4s4 1 1 1
s3s3 1 1
s2s2 0 1
s1s1 -∞
s0s0 1
There are two sign changes in the first column of Routh table. Hence, the control system is
unstable. All the Elements of any row of the Routh array are zero
∙ Write the auxilary equation, A(s) of the row, which is just above the row of zeros.
∙ Differentiate the auxiliary equation, A(s) with respect to s. Fill the row of zeros with these coefficients.
Example
Let us find the stability of the control system having characteristic equation,
s5+3s4+s3+3s2+s+3=0s5+3s4+s3+3s2+s+3=0
Step 1 − Verify the necessary condition for the Routh-Hurwitz stability.
All the coefficients of the given characteristic polynomial are positive. So, the control system satisfied the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
s5s5 1 1 1
s4s4 31 31 31
s2s2
s1s1
s0s0
The row s4s4 elements have the common factor of 3. So, all these elements are divided by 3.
Special case (ii) − All the elements of row s3s3 are zero. So, write the auxiliary equation, A(s) of the row
s4s4. A(s)=s4+s2+1A(s)=s4+s2+1
s5s5 1 1 1
s4s4 1 1 1
s3s3 42 21
s1s1 (0.5×1)−(1×2)0.5=−1.50.5=−3(0.5×1)−(1×2)0.5=−1.50.5=−3
s0s0 1
There are two sign changes in the first column of Routh table. Hence, the control system is unstable.
In the Routh-Hurwitz stability criterion, we can know whether the closed loop poles are in on left half of the ‘s’ plane or on the
right half of the ‘s’ plane or on an imaginary axis. So, we can’t find the nature of the control system. To overcome this
limitation, there is a technique known as the root locus. We will discuss this technique in the next two chapters.
ROOT LOCUS
In the root locus diagram, we can observe the path of the closed loop poles. Hence, we can identify the nature of the control
system. In this technique, we will use an open loop transfer function to know the stability of the closed loop control system.
We know that, the characteristic equation of the closed loop control system is
1+G(s)H(s)=01+G(s)H(s)=0
Where,
∙ N(s) represents the numerator term having (factored) nth order polynomial of ‘s’.
∙ D(s) represents the denominator term having (factored) mth order polynomial of ‘s’.
⇒D(s)+KN(s)=0⇒D(s)+KN(s)=0
Case 1 − K = 0
If K=0K=0, then D(s)=0D(s)=0.
That means, the closed loop poles are equal to open loop poles when K is zero.
Case 2 − K = ∞
K(1K+N(s)D(s))=0⇒1K+N(s)D(s)=0K(1K+N(s)D(s))=0⇒1K+N(s)D(s)=0
If K=∞K=∞, then N(s)=0N(s)=0. It means the closed loop poles are equal to the open loop zeros when K is
infinity. From above two cases, we can conclude that the root locus branches start at open loop poles and end at open
loop zeros.
⇒G(s)H(s)=−1+j0⇒G(s)H(s)=−1+j0
The angle condition is the point at which the angle of the open loop transfer function is an odd multiple of 1800.
Magnitude of G(s)H(s)G(s)H(s) is -
|G(s)H(s)|=(−1)2+02−−−−−−−−−√=1|G(s)H(s)|=(−1)2+02=1
The magnitude condition is that the point (which satisfied the angle condition) at which the magnitude of the open loop transfer
function is one.
Rule 1 − Locate the open loop poles and zeros in the ‘s’ plane.
We know that the root locus branches start at the open loop poles and end at open loop zeros. So, the number of root locus
branches N is equal to the number of finite open loop poles P or the number of finite open loop zeros Z, whichever is greater.
If the angle of the open loop transfer function at a point is an odd multiple of 1800, then that point is on the root locus. If odd
number of the open loop poles and zeros exist to the left side of a point on the real axis, then that point is on the root locus
branch. Therefore, the branch of points which satisfies this condition is the real axis of the root locus branch.
∙ If P=ZP=Z, then all the root locus branches start at finite open loop poles and end at finite open loop zeros. ∙ If
P>ZP>Z , then ZZ number of root locus branches start at finite open loop poles and end at finite open loop zeros and
P−ZP−Z number of root locus branches start at finite open loop poles and end at infinite open loop zeros. ∙ If P<ZP<Z , then
P number of root locus branches start at finite open loop poles and end at finite open loop zeros and Z−PZ−P number of root
locus branches start at infinite open loop poles and end at finite open loop zeros. So, some of the root locus branches approach
infinity, when P≠ZP≠Z. Asymptotes give the direction of these root locus branches. The intersection point of asymptotes on
the real axis is known as centroid.
We can calculate the centroid α by using this formula,
α=∑Realpartoffiniteopenlooppoles−∑RealpartoffiniteopenloopzerosP−Zα=∑Realpartoffiniteopenlooppoles−∑Realpartoffiniteo
penloopzerosP−Z
The formula for the angle of asymptotes θ is
θ=(2q+1)1800P−Zθ=(2q+1)1800P−Z
Where,
q=0,1,2,....,(P−Z)−1q=0,1,2,....,(P−Z)−1
Rule 5 − Find the intersection points of root locus branches with an imaginary axis.
We can calculate the point at which the root locus branch intersects the imaginary axis and the value of K at that point by using
the Routh array method and special case (ii).
∙ If all elements of any row of the Routh array are zero, then the root locus branch intersects the imaginary axis and
vice-versa.
∙ Identify the row in such a way that if we make the first element as zero, then the elements of the entire row are zero.
Find the value of Kfor this combination.
∙ Substitute this K value in the auxiliary equation. You will get the intersection point of the root locus branch with an
imaginary axis.
∙ If there exists a real axis root locus branch between two open loop poles, then there will be a break-away point in
between these two open loop poles.
∙ If there exists a real axis root locus branch between two open loop zeros, then there will be a break-in point in between
these two open loop zeros.
Note − Break-away and break-in points exist only on the real axis root locus branches.
The Angle of departure and the angle of arrival can be calculated at complex conjugate open loop poles and complex conjugate
open loop zeros respectively.
ϕ=∑ϕP−∑ϕZϕ=∑ϕP−∑ϕZ
Example
Let us now draw the root locus of the control system having open loop transfer function,
G(s)H(s)=Ks(s+1)(s+5)G(s)H(s)=Ks(s+1)(s+5)
Step 1 − The given open loop transfer function has three poles at s=0,s=−1s=0,s=−1 and s=−5s=−5. It doesn’t have any
zero. Therefore, the number of root locus branches is equal to the number of poles of the open loop transfer
function. N=P=3N=P=3
The three poles are located are shown in the above figure. The line segment between s=−1s=−1 and s=0s=0 is one branch of
root locus on real axis. And the other branch of the root locus on the real axis is the line segment to the left ofs=−5s=−5. Step
2 − We will get the values of the centroid and the angle of asymptotes by using the given formulae.
Centroid α=−2α=−2
The angle of asymptotes are θ=600,1800θ=600,1800 and 30003000.
The centroid and three asymptotes are shown in the following figure.
Step 3 − Since two asymptotes have the angles of 600600 and 30003000, two root locus branches intersect the imaginary
axis. By using the Routh array method and special case(ii), the root locus branches intersects the imaginary axis at j5– √j5
and −j5–√−j5.
There will be one break-away point on the real axis root locus branch between the poles s=−1s=−1 and s=0s=0. By
following the procedure given for the calculation of break-away point, we will get it as s=−0.473s=−0.473. The root locus
diagram for the given control system is shown in the following figure.
In this way, you can draw the root locus diagram of any control system and observe the movement of poles of the closed loop
transfer function.
From the root locus diagrams, we can know the range of K values for different types of damping.
∙ If we include a pole in the open loop transfer function, then some of root locus branches will move towards right half of
‘s’ plane. Because of this, the damping ratio δδ decreases. Which implies, damped frequency ωdωd increases and
the time domain specifications like delay time tdtd, rise time trtr and peak time tptp decrease. But, it effects the
system stability.
∙ If we include a zero in the open loop transfer function, then some of root locus branches will move towards left half of
‘s’ plane. So, it will increase the control system stability. In this case, the damping ratio δδincreases. Which implies,
damped frequency ωdωd decreases and the time domain specifications like delay time tdtd, rise time trtr and peak
time tptp increase.
So, based on the requirement, we can include (add) the open loop poles or zeros to the transfer
If a sinusoidal signal is applied as an input to a Linear Time-Invariant (LTI) system, then it produces the steady state output,
which is also a sinusoidal signal. The input and output sinusoidal signals have the same frequency, but different amplitudes and
phase angles.
r(t)=Asin(ω0t)r(t)=Asin(ω0t)
G(s)=G(jω)G(s)=G(jω)
c(t)=A|G(jω0)|sin(ω0t+∠G(jω0))c(t)=A|G(jω0)|sin(ω0t+∠ G(jω0))
∙ The amplitude of the output sinusoidal signal is obtained by multiplying the amplitude of the input sinusoidal signal
and the magnitude of G(jω)G(jω) at ω=ω0ω=ω0.
∙ The phase of the output sinusoidal signal is obtained by adding the phase of the input sinusoidal signal and the phase of
G(jω)G(jω) at ω=ω0ω=ω0.
Where,
Here, f0f0 is the frequency of the input sinusoidal signal. Similarly, you can follow the same procedure for closed loop control
system.
Consider the transfer function of the second order closed loop control system as,
T(s)=C(s)R(s)=ω2ns2+2δωns+ω2nT(s)=C(s)R(s)=ωn2s2+2δωns+ωn2
Substitute, s=jωs=jω in the above equation.
T(jω)=ω2n(jω)2+2δωn(jω)+ω2nT(jω)=ωn2(jω)2+2δωn(jω)+ωn2
⇒T(jω)=ω2n−ω2+2jδωωn+ω2n=ω2nω2n(1−ω2ω2n+2jδωωn)⇒T(jω)=ωn2−ω2+2jδωωn+ωn2=ωn2ωn2(1−ω2ωn2+2jδ
ω ωn)
⇒T(jω)=1(1−ω2ω2n)+j(2δωωn)⇒T(jω)=1(1−ω2ωn2)+j(2δωωn)
Magnitude of T(jω)T(jω) is -
M=|T(jω)|=1(1−u2)2+(2δu)2−−−−−−−−−−−−−−√M=|T(jω)|=1(1−u2)2+(2δu)2
Phase of T(jω)T(jω) is -
∠T(jω)=−tan−1(2δu1−u2)∠ T(jω)=−tan−1(2δu1−u2)
Resonant Frequency
It is the frequency at which the magnitude of the frequency response has peak value for the first time. It is denoted by ωrωr.
At ω=ωrω=ωr, the first derivate of the magnitude of T(jω)T(jω) is zero.
Differentiate MM with respect to uu.
dMdu=−12[(1−u2)2+(2δu)2]−32[2(1−u2)(−2u)+2(2δu)(2δ)]dMdu=−12[(1−u2)2+(2δu)2]−32[2(1−u2)(−2u)+2(2δu)
( 2δ)]
⇒dMdu=−12[(1−u2)2+(2δu)2]−32[4u(u2−1+2δ2)]⇒dMdu=−12[(1−u2)2+(2δu)2]−32[4u(u2−1+2δ2)]
] ⇒4ur(u2r−1+2δ2)=0⇒4ur(ur2−1+2δ2)=0
⇒u2r−1+2δ2=0⇒ur2−1+2δ2=0
⇒u2r=1−2δ2⇒ur2=1−2δ2
⇒ur=1−2δ2−−−−−−√⇒ur=1−2δ2
Substitute, ur=ωrωnur=ωrωn in the above equation.
ωrωn=1−2δ2−−−−−−√ωrωn=1−2δ2
⇒ωr=ωn1−2δ2−−−−−−√⇒ωr=ωn1−2δ2
Resonant Peak
It is the peak (maximum) value of the magnitude of T(jω)T(jω). It is denoted by MrMr.
At u=uru=ur, the Magnitude of T(jω)T(jω) is -
Mr=1(1−u2r)2+(2δur)2−−−−−−−−−−−−−−−√Mr=1(1−ur2)2+(2δur)2
2 ⇒Mr=12δ1−δ2−−−−−√⇒Mr=12δ1−δ2
Resonant peak in frequency response corresponds to the peak overshoot in the time domain transient response for certain
values of damping ratio δδ. So, the resonant peak and peak overshoot are correlated to each other.
Bandwidth
It is the range of frequencies over which, the magnitude of T(jω)T(jω) drops to 70.7% from its zero frequency
value. At ω=0ω=0, the value of uu will be zero.
Substitute, u=0u=0 in M.
M=1(1−02)2+(2δ(0))2−−−−−−−−−−−−−−−−√=1M=1(1−02)2+(2δ(0))2=1
⇒2=(1−u2b)2+(2δ)2u2b⇒2=(1−ub2)2+(2δ)2ub2
Let, u2b=xub2=x
⇒2=(1−x)2+(2δ)2x⇒2=(1−x)2+(2δ)2x
⇒x2+(4δ2−2)x−1=0⇒x2+(4δ2−2)x−1=0
⇒x=−(4δ2−2)±(4δ2−2)2+4−−−−−−−−−−−√2⇒x=−(4δ2−2)±(4δ2−2)2+42
x=1−2δ2+(2δ2−1)2+1−−−−−−−−−−−√x=1−2δ2+(2δ2−1)2+1
⇒x=1−2δ2+(2−4δ2+4δ4)−−−−−−−−−−−−√⇒x=1−2δ2+(2−4δ2+4δ4)
Substitute, x=u2b=ω2bω2nx=ub2=ωb2ωn2
ω2bω2n=1−2δ2+(2−4δ2+4δ4)−−−−−−−−−−−−√ωb2ωn2=1−2δ2+(2−4δ2+4δ4)
⇒ωb=ωn1−2δ2+(2−4δ2+4δ4)−−−−−−−−−−−−√−−−−−−−−−−−−−−−−−−−−−−√⇒ωb=ωn1−2δ2+(2−4δ2+4δ4)
Bandwidth ωbωb in the frequency response is inversely proportional to the rise time trtr in the time domain transient
∙ Magnitude plot
∙ Phase plot
In both the plots, x-axis represents angular frequency (logarithmic scale). Whereas, yaxis represents the magnitude (linear
scale) of open loop transfer function in the magnitude plot and the phase angle (linear scale) of the open loop transfer function
in the phase plot.
The magnitude of the open loop transfer function in dB is -
M=20log|G(jω)H(jω)|M=20log|G(jω)H(jω)|
ϕ=∠G(jω)H(jω)ϕ=∠ G(jω)H(jω)
Const KK 00 20logK20logK 00
ant
The magnitude plot is having magnitude of 0 dB upto ω=1τω=1τ rad/sec. From ω=1τω=1τ rad/sec, it is having a slope of 20
dB/dec. In this case, the phase plot is having phase angle of 0 degrees up to ω=1τω=1τ rad/sec and from here, it is having
0
phase angle of 90 . This Bode plot is called the asymptotic Bode plot.
As the magnitude and the phase plots are represented with straight lines, the Exact Bode plots resemble the asymptotic Bode
plots. The only difference is that the Exact Bode plots will have simple curves instead of straight lines.
Similarly, you can draw the Bode plots for other terms of the open loop transfer function which are given in the
∙ Represent the open loop transfer function in the standard time constant form.
∙ Consider the starting frequency of the Bode plot as 1/10th of the minimum corner frequency or 0.1 rad/sec whichever is
smaller value and draw the Bode plot upto 10 times maximum corner frequency.
∙ Draw the magnitude plots for each term and combine these plots properly.
∙ Draw the phase plots for each term and combine these plots properly.
Note − The corner frequency is the frequency at which there is a change in the slope of the magnitude
plot. Example
Consider the open loop transfer function of a closed loop control system
G(s)H(s)=10s(s+2)(s+5)G(s)H(s)=10s(s+2)(s+5)
Let us convert this open loop transfer function into standard time constant form.
G(s)H(s)=10s2(s2+1)5(s5+1)G(s)H(s)=10s2(s2+1)5(s5+1)
⇒G(s)H(s)=s(1+s2)(1+s5)⇒G(s)H(s)=s(1+s2)(1+s5)
So, we can draw the Bode plot in semi log sheet using the rules mentioned earlier.
∙ If the phase cross over frequency ωpcωpc is greater than the gain cross over frequency ωgcωgc, then the control
system is stable.
∙ If the phase cross over frequency ωpcωpc is equal to the gain cross over frequency ωgcωgc, then the control system
is marginally stable.
∙ If the phase cross over frequency ωpcωpc is less than the gain cross over frequency ωgcωgc, then the control system is
unstable.
Gain Margin
Gain margin GMGM is equal to negative of the magnitude in dB at phase cross over frequency.
GM=20log(1Mpc)=20logMpcGM=20log(1Mpc)=20logMpc
Where, MpcMpc is the magnitude at phase cross over frequency. The unit of gain margin (GM) is
Where, ϕgcϕgc is the phase angle at gain cross over frequency. The unit of phase margin is degrees. The
stability of the control system based on the relation between gain margin and phase margin is listed below.
∙ If both the gain margin GMGM and the phase margin PMPM are positive, then the control system is stable. ∙ If both
the gain margin GMGM and the phase margin PMPM are equal to zero, then the control system
is marginally stable.
∙ If the gain margin GMGM and / or the phase margin PMPM are/is negative, then the control system is unstable.
Control Systems - Polar Plots
In the previous chapters, we discussed the Bode plots. There, we have two separate plots for both magnitude and phase as the
function of frequency. Let us now discuss about polar plots. Polar plot is a plot which can be drawn between magnitude and
phase. Here, the magnitudes are represented by normal values only.
The Polar plot is a plot, which can be drawn between the magnitude and the phase angle of G(jω)H(jω)G(jω)H(jω) by
varying ωω from zero to ∞. The polar graph sheet is shown in the following figure.
This graph sheet consists of concentric circles and radial lines. The concentric circles and the radial lines represent the
magnitudes and phase angles respectively. These angles are represented by positive values in anti-clock wise direction.
Similarly, we can represent angles with negative values in clockwise direction. For example, the angle 2700in anti-clock wise
direction is equal to the angle −900in clockwise direction.
Example
Consider the open loop transfer function of a closed loop control system.
G(s)H(s)=5s(s+1)(s+2)G(s)H(s)=5s(s+1)(s+2)
Let us draw the polar plot for this control system using the above rules.
M=5ω(ω2+1−−−−−√)(ω2+4−−−−−√)M=5ω(ω2+1)(ω2+4)
ϕ=−900−tan−1ω−tan−1ω2ϕ=−900−tan−1ω−tan−1ω2
Step 2 − The following table shows the magnitude and the phase angle of the open loop transfer function at ω=0ω=0 rad/sec
and ω=∞ω=∞ rad/sec.
0 ∞ -90 or 270
∞ 0 -270 or 90
So, the polar plot starts at (∞,−900) and ends at (0,−2700). The first and the second terms within the brackets indicate the
magnitude and phase angle respectively.
Step 3 − Based on the starting and the ending polar co-ordinates, this polar plot will intersect the negative real axis. The phase
angle corresponding to the negative real axis is −1800 or 1800. So, by equating the phase angle of the open loop transfer
function to either −1800 or 1800, we will get the ωω value as 2–√2.
By substituting ω=2–√ω=2 in the magnitude of the open loop transfer function, we will get M=0.83M=0.83. Therefore,
the polar plot intersects the negative real axis when ω=2–√ω=2 and the polar coordinate is (0.83,−1800). So, we can draw the
polar plot with the above information on the polar graph sheet.
∙ If the enclosed ‘s’ plane closed path contains only poles, then the direction of the encirclement in the
G(s)H(s)G(s)H(s) plane will be opposite to the direction of the enclosed closed path in the ‘s’ plane. ∙ If the enclosed
‘s’ plane closed path contains only zeros, then the direction of the encirclement in the G(s)H(s)G(s)H(s) plane will be
in the same direction as that of the enclosed closed path in the ‘s’ plane.
Let us now apply the principle of argument to the entire right half of the ‘s’ plane by selecting it as a closed path. This selected
path is called the Nyquistcontour.
We know that the closed loop control system is stable if all the poles of the closed loop transfer function are in the left half of
the ‘s’ plane. So, the poles of the closed loop transfer function are nothing but the roots of the characteristic equation. As the
order of the characteristic equation increases, it is difficult to find the roots. So, let us correlate these roots of the characteristic
equation as follows.
∙ The Poles of the characteristic equation are same as that of the poles of the open loop transfer function. ∙
The zeros of the characteristic equation are same as that of the poles of the closed loop transfer function.
We know that the open loop control system is stable if there is no open loop pole in the the right half of the ‘s’
plane. i.e.,P=0⇒N=−ZP=0⇒N=−Z
We know that the closed loop control system is stable if there is no closed loop pole in the right half of the ‘s’
plane. i.e.,Z=0⇒N=PZ=0⇒N=P
Nyquist stability criterion states the number of encirclements about the critical point (1+j0) must be equal to the poles of
characteristic equation, which is nothing but the poles of the open loop transfer function in the right half of the ‘s’ plane. The
shift in origin to (1+j0) gives the characteristic equation plane.
∙ Locate the poles and zeros of open loop transfer function G(s)H(s)G(s)H(s) in ‘s’ plane.
∙ Draw the polar plot by varying ωω from zero to infinity. If pole or zero present at s = 0, then varying ωω from 0+ to
infinity for drawing polar plot.
∙ Draw the mirror image of above polar plot for values of ωω ranging from −∞ to zero (0−if any pole or zero present at
s=0).
∙ The number of infinite radius half circles will be equal to the number of poles or zeros at origin. The infinite radius half
circle will start at the point where the mirror image of the polar plot ends. And this infinite radius half circle will end
at the point where the polar plot starts.
After drawing the Nyquist plot, we can find the stability of the closed loop control system using the Nyquist stability criterion.
If the critical point (-1+j0) lies outside the encirclement, then the closed loop control system is absolutely stable. Stability
Analysis using Nyquist Plots
From the Nyquist plots, we can identify whether the control system is stable, marginally stable or unstable based on the values
of these parameters.
∙ If the phase cross over frequency ωpcωpc is greater than the gain cross over frequency ωgcωgc, then the control
system is stable.
∙ If the phase cross over frequencyωpcωpc is equal to the gain cross over frequency ωgcωgc, then the control system
is marginally stable.
∙ If phase cross over frequency ωpcωpc is less than gain cross over frequency ωgcωgc, then the control system is
unstable.
Gain Margin
The gain margin GMGM is equal to the reciprocal of the magnitude of the Nyquist plot at the phase cross over
frequency. GM=1MpcGM=1Mpc
Where, MpcMpc is the magnitude in normal scale at the phase cross over frequency.
Phase Margin
The phase margin PMPM is equal to the sum of 1800and the phase angle at the gain cross over
frequency. PM=1800+ϕgcPM=1800+ϕgc
Where, ϕgcϕgc is the phase angle at the gain cross over frequency.
The stability of the control system based on the relation between the gain margin and the phase margin is listed below.
∙ If the gain margin GMGM is greater than one and the phase margin PMPMis positive, then the control system
is stable.
∙ If the gain margin GMGM is equal to one and the phase margin PMPM is zero degrees, then the control system
is marginally stable.
∙ If the gain margin GMGM is less than one and / or the phase margin PMPMis negative, then the control system
is unstable.
Vo(s)Vi(s)=1α(s+1τs+1ατ)Vo(s)Vi(s)=1α(s+1τs+1ατ)
Where,
τ=R2Cτ=R2C
α=R1+R2R2α=R1+R2R2
So, in order to produce the phase lag at the output of this compensator, the phase angle of the transfer function should be
negative. This will happen when α>1α>1.
Lead Compensator
The lead compensator is an electrical network which produces a sinusoidal output having phase lead when a sinusoidal input is
applied. The lead compensator circuit in the ‘s’ domain is shown in the following figure.
Here, the capacitor is parallel to the resistor R1R1 and the output is measured across resistor
$R_2. The transfer function of this lead compensator is -
Vo(s)Vi(s)=β(sτ+1βsτ+1)Vo(s)Vi(s)=β(sτ+1βsτ+1)
Where,
τ=R1Cτ=R1C
β=R2R1+R2β=R2R1+R2
From the transfer function, we can conclude that the lead compensator has pole at s=−1βs=−1β and zero
ats=−1βτs=−1βτ. Substitute, s=jωs=jω in the transfer function.
Vo(jω)Vi(jω)=β(jωτ+1βjωτ+1)Vo(jω)Vi(jω)=β(jωτ+1βjωτ+1)
Lag-Lead Compensator
Lag-Lead compensator is an electrical network which produces phase lag at one frequency region and phase lead at other
frequency region. It is a combination of both the lag and the lead compensators. The lag-lead compensator circuit in the ‘s’
domain is shown in the following figure.
This circuit looks like both the compensators are cascaded. So, the transfer function of this circuit will be the product of
transfer functions of the lead and the lag compensators.
Vo(s)Vi(s)=β(sτ1+1βsτ1+1)1α(s+1τ2s+1ατ2)Vo(s)Vi(s)=β(sτ1+1βsτ1+1)1α(s+1τ2s+1ατ2)
We know αβ=1αβ=1.
⇒Vo(s)Vi(s)=(s+1τ1s+1βτ1)(s+1τ2s+1ατ2)⇒Vo(s)Vi(s)=(s+1τ1s+1βτ1)(s+1τ2s+1ατ2)
Where,
τ1=R1C1τ1=R1C1
τ2=R2C2τ2=R2C2
Proportional Controller
The proportional controller produces an output, which is proportional to error signal.
u(t)∝e(t)u(t)∝e(t)
⇒u(t)=KPe(t)⇒u(t)=KPe(t)
U(s)=KPE(s)U(s)=KPE(s)
U(s)E(s)=KPU(s)E(s)=KP
The block diagram of the unity negative feedback closed loop control system along with the proportional controller is shown in
the following figure.
The proportional controller is used to change the transient response as per the requirement.
Derivative Controller
The derivative controller produces an output, which is derivative of the error signal.
u(t)=KDde(t)dtu(t)=KDde(t)dt
U(s)=KDsE(s)U(s)=KDsE(s)
U(s)E(s)=KDsU(s)E(s)=KDs
The derivative controller is used to make the unstable control system into a stable one.
Integral Controller
The integral controller produces an output, which is integral of the error signal.
u(t)=KI∫e(t)dtu(t)=KI∫e(t)dt
Apply Laplace transform on both the sides -
U(s)=KIE(s)sU(s)=KIE(s)s
U(s)E(s)=KIsU(s)E(s)=KIs