Automation & Control-Study Material (Unit 1,2,3,4)
Automation & Control-Study Material (Unit 1,2,3,4)
OF ME
UNIT - I:
Introduction: Concept of control system, Classification of control systems - Open loop and
closed loop control systems, Differences, Examples of control systems- Effects of feedback,
Feedback Characteristics.
Transfer Function Representation: Block diagram algebra, Determining the Transfer function
from Block Diagrams, Signal flow graphs (SFG) - Reduction using Mason’s gain formula-
Transfer function of SFG’s.
UNIT - II:
Time Response Analysis: Standard test signals, Time response of first order systems,
Characteristic Equation of Feedback control systems, Transient response of second order
systems - Time domain specifications, Steady state response, Steady state errors and error
constants. PID controllers: Effects of proportional derivative, proportional integral systems on
steady state error.
UNIT - III:
State Space Analysis of Continuous Systems: Concepts of state, state variables and state
model, Derivation of state models from block diagrams, Diagonalization, Solving the time
invariant state equations, State Transition Matrix and it’s properties, Concepts of
Controllability and Observability.
UNIT - IV:
Stability Analysis in S-Domain: The concept of stability – Routh-Hurwitz’s stability criterion –
qualitative stability and conditional stability – Limitations of Routh-Hurwitz’s stability.
Root Locus Technique: Concept of root locus - Construction of root locus, Effects of adding
poles and zeros to G(s) H(s) on the root loci.
Frequency Response Analysis: Introduction, Frequency domain specifications, Bode plot
diagrams-Determination of Phase margin and Gain margin, Stability analysis from Bode
plots, Polar plots.
DEPT. OF ME
UNIT-I
INTRODUCTION
A control system manages commands, directs or regulates the behavior of other devices or
systems using control loops. It can range from a single home heating controller using
a thermostat controlling a domestic boiler to large Industrial control systems which are used
for controlling processes or machines. 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.
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.
DEPT. OF ME
In open loop control systems, output is not fed-back to the input. So, the control action is
independent of the desired output.
The following figure shows the block diagram of the open loop control system.
In closed loop control systems, output is fed back to the input. So, the control action is
dependent on the desired output.
The following figure shows the block diagram of negative feedback closed loop control
system.
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.
DEPT. OF ME
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.
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)R(s) and feedback output. The following
figure shows the block diagram of positive feedback control system
DEPT. OF ME
he 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,
Negative Feedback
Negative feedback reduces the error between the reference input, R(s)R(s) and system
output. The following figure shows the block diagram of the negative feedback control
system.
Where,
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.
Consider an open loop control system with noise signal as shown below.
DEPT. OF ME
The control systems can be represented with a set of mathematical equations known
as mathematical model. These models are useful for analysis and design of control systems.
Analysis of control system means finding the output when we know the input and
DEPT. OF ME
mathematical model. Design of control system means finding the mathematical model when
we know the input and the output.
Block Diagrams
Block diagrams consist of a single block or a combination of blocks. These are used to
represent the control systems in pictorial form.
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).
DEPT. OF ME
Summing Point
The summing point is represented with a circle having cross (X) inside it. It has two or more
inputs and single output. It produces the algebraic sum of the inputs. It also performs the
summation or subtraction or combination of summation and subtraction of the inputs based
on the polarity of the inputs. Let us see these three operations one by one.
The following figure shows the summing point with two inputs (A, B) and one output (Y).
Here, the inputs A and B have a positive sign. So, the summing point produces the output, Y
as sum of A and B i.e. = A + B.
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 i.e
Y = A + (-B) = A - B.
DEPT. OF ME
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.
DEPT. OF ME
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.
Block diagram algebra is nothing but the algebra involved with the basic elements of the
block diagram. This algebra deals with the pictorial representation of algebraic equations.
Series Connection
Series connection is also called cascade connection. In the following figure, two blocks
having transfer functions G1(s)G1(s) and G2(s)G2(s) are connected in series.
DEPT. OF ME
That means we can represent the series connection of two blocks with a single block. The
transfer function of this single block is the product of the transfer functions of those two
blocks. The equivalent block diagram is shown below.
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)G1(s) and G2(s)G2(s) are connected in parallel.
The outputs of these two blocks are connected to the summing point.
That means we can represent the parallel connection of two blocks with a single block. The
transfer function of this single block is the sum of the transfer functions of those two blocks.
The equivalent block diagram is shown below.
DEPT. OF ME
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)G(s) and H(s)H(s) form a closed loop.
This means we can represent the negative feedback connection of two blocks with a single
block. The transfer function of this single block is the closed loop transfer function of the
negative feedback. The equivalent block diagram is shown below.
DEPT. OF ME
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.,
There are two possibilities of shifting summing points with respect to blocks −
The first term ‘G(s)R(s)′‘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)G(s).
It is having the input X(s)X(s) and the output of this block is given as input to summing point
instead of X(s)X(s). This block diagram is shown in the following figure.
DEPT. OF ME
The first term ‘G(s)R(s)′ is same in both equations. But, there is difference in the second term.
In order to get the second term also same, we require one more block 1/G(s). It is having the
DEPT. OF ME
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.
Shifting a Take-off Point form a Position before a Block to a position after the Block
Consider the block diagram shown in the following figure. In this case, the take-off point is
present before the block.
When you shift the take-off point after the block, the output Y(s) will be same. But, there is
difference in X(s) value. So, in order to get the same X(s) value, we require one more
DEPT. OF ME
block 1/G(s). It is having the input Y(s) and the output is X(s) this block diagram is shown in
the following figure.
Shifting Take-off Point from a Position after a Block to a position before the Block
Consider the block diagram shown in the following figure. Here, the take-off point is present
after the block.
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.
DEPT. OF ME
The concepts discussed in the previous chapter are helpful for reducing (simplifying) the
block diagrams.
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 2 − Repeat step 1 for remaining inputs.
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).
Examples:
1. Consider the block diagram shown in the following figure. Let us simplify (reduce) this
block diagram using the block diagram reduction rules.
DEPT. OF ME
DEPT. OF ME
DEPT. OF ME
DEPT. OF ME
Signal flow graph is a graphical representation of algebraic equations. In this chapter, let us
discuss the basic concepts related signal flow graph and also learn how to draw signal flow
graphs.
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.
Example
Let us consider the following signal flow graph to identify these nodes.
DEPT. OF ME
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.
Represent all the signals, variables, summing points and take-off points of block
diagram as nodes in signal flow graph.
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.
DEPT. OF ME
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 labeled 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 G1and G2.
Let us now discuss the Mason’s Gain Formula. Suppose there are ‘N’ forward paths in a signal
flow graph. The gain between the input and the output nodes of a signal flow graph is
DEPT. OF ME
nothing but the transfer function of the system. It can be calculated by using Mason’s gain
formula.
Where,
Δ=1−(sum of all individual loop gains) +(sum of gain products of all possible two
nontouching loops)−(sum of gain products of all possible three nontouching loops)
+….
Δ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.
DEPT. OF ME
Loop
The path that starts from one node and ends at the same node is known as a loop. Hence, it
is a closed path.
y2→y3→y2, y4→y5→y4.
First non-touching loops pair is -
Gain product of first non-touching loops pair l1l4=bjdi
Second non-touching loops pair is - y2→y3→y2, y5→y5.
Gain product of second non-touching loops pair is l1l5=bjf
Higher number of (more than two) non-touching loops are not present in this signal flow
graph.We know,
DEPT. OF ME
DEPT. OF ME
DEPT. OF ME
UNIT-II
TIME RESPONSE ANALYSIS
We can analyze the response of the control systems in both the time domain and the
frequency domain. We will discuss frequency response analysis of control systems in later
chapters. Let us now discuss about the time response analysis of control systems.
Transient response
Steady state response
The response of control system in time domain is shown in the following figure.
Where,
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.
Example
Let us find the transient and steady state terms of the time response of the control system
Here, the second term will be zero as t denotes infinity. So, this is the transient term.
And the first term 10 remains even as t approaches infinity. So, this is the steady state term.
Standard Test Signals
The standard test signals are impulse, step, ramp and parabolic. These signals are used to
know the performance of the control systems using time response of the output.
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’.
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’.
DEPT. OF ME
So, the unit ramp signal exists for all positive values of‘t’ including zero. And its value
increases linearly with respect to‘t’ during this interval. The value of unit ramp signal is zero
for all negative values of‘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’.
In this chapter, let us discuss the time response of the first order system. Consider the
following block diagram of the closed loop control system. Here, an open loop transfer
function, 1/sT is connected with a unity negative feedback.
DEPT. OF ME
So, r(t)=δ(t)
Apply Laplace transform on both the sides.
R(s) =1
The unit impulse response, c(t) is an exponential decaying signal for positive values of ‘t’ and
it is zero for negative values of ‘t’.
So, r(t)=u(t)
DEPT. OF ME
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)+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=−T
The unit step response, c(t) has both the transient and the steady state terms.
DEPT. OF ME
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.
On both the sides, the denominator term is the same. So, they will get cancelled by each
other. Hence, equate the numerator terms.
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=−T
Similarly, substitute B = −T and equate the coefficient of s2 terms on both the sides. You will
get C=T2
Substitute A = 1, B = −T and C=T2 in the partial fraction expansion of C(s).
The unit ramp response, c(t) has both the transient and the steady state terms.
The unit ramp response, c(t) follows the unit ramp input signal for all positive values of t.
But, there is a deviation of T units from the input signal.
\
DEPT. OF ME
The unit parabolic response, c(t) has both the transient and the steady state terms.
From these responses, we can conclude that the first order control systems are not stable
with the ramp and parabolic inputs because these responses go on increasing even at infinite
amount of time. The first order control systems are stable with impulse and step inputs
because these responses have bounded output. But, the impulse response doesn’t have
steady state term. So, the step signal is widely used in the time domain for analyzing the
control systems from their responses.
In this chapter, let us discuss the time response of second order system. Consider the
following block diagram of closed loop control system. Here, an open loop transfer
function, ωn2 / s(s+2δωn) is connected with a unity negative feedback.
DEPT. OF ME
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.
Where,
Follow these steps to get the response (output) of the second order system in the time
domain.
So, the unit step response of the second order system is having damped oscillations
(decreasing amplitude) when ‘δ’ lies between zero and one.
Case 4: δ > 1
We can modify the denominator term of the transfer function as follows −
DEPT. OF ME
Since it is over damped, the unit step response of the second order system when δ > 1 will
never reach step input in the steady state.
Follow the procedure involved while deriving step response by considering the value
of R(s) as 1 instead of 1/s.
Do the differentiation of the step response.
The following table shows the impulse response of the second order system for 4 cases of
the damping ratio.
DEPT. OF ME
In this chapter, let us discuss the time domain specifications of the second order system. The
step response of the second order system for the underdamped case is shown in the
following figure.
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
DEPT. OF ME
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.
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.
At t = t1 = 0, c(t) = 0.
We know that the final value of the step response is one. Therefore, at t=t2, the value of step
response is one. Substitute, these values in the following equation.
DEPT. OF ME
From above equation, we can conclude that the rise time tr and the damped
frequency ω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 tp. At t=tp the first derivate of the response is zero.
We know the step response of second order system for under-damped case is
DEPT. OF ME
From the above equation, we can conclude that the peak time tp and the damped
frequency ω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.
Where,c(tp) is the peak value of the response, c(∞) is the final (steady state) value of the
response.
From the above equation, we can conclude that the percentage of peak overshoot %Mp will
decrease if the damping ratio δ increases.
Settling time
It is the time required for the response to reach the steady state and stay within the specified
tolerance bands around the final value. In general, the tolerance bands are 2% and 5%. The
settling time is denoted by ts.
DEPT. OF ME
By equating these two transfer functions, we will get the un-damped natural frequency ωn as
2 rad/sec and the damping ratio δ as 0.5.
We know the formula for damped frequency ωd as
DEPT. OF ME
Substitute the above necessary values in the formula of each time domain specification and
simplify in order to get the values of time domain specifications for given transfer function.
The following table shows the formulae of time domain specifications, substitution of
necessary values and the final values
DEPT. OF ME
The deviation of the output of control system from desired response during steady state is
known as steady state error. It is represented as ess. We can find steady state error using
the final value theorem as follows.
Where,
The following table shows the steady state errors and the error constants for standard input
signals like unit step, unit ramp & unit parabolic signals.
Where, Kp, Kv and Ka are position error constant, velocity error constant and acceleration
error constant respectively.
Note − If any of the above input signals has the amplitude other than unity, then multiply
corresponding steady state error with that amplitude.
DEPT. OF ME
Note − We can’t define the steady state error for the unit impulse signal because, it exists
only at origin. So, we can’t compare the impulse response with the unit impulse input
as t denotes infinity
We will get the overall steady state error, by adding the above three steady state errors.
ess = ess1+ess2+ess3
⇒ess=0+0+1=1⇒ess=0+0+1=1
Therefore, we got the steady state error ess as 1 for this example.
Steady State Errorsfor Non-Unity Feedback Systems
Consider the following block diagram of closed loop control system, which is having non unity
negative feedback.
DEPT. OF ME
We can find the steady state errors only for the unity feedback systems. So, we have to
convert the non-unity feedback system into unity feedback system. For this, include one
unity positive feedback path and one unity negative feedback path in the above block
diagram. The new block diagram looks like as shown below.
Simplify the above block diagram by keeping the unity negative feedback as it is. The
following is the simplified block diagram
DEPT. OF ME
This block diagram resembles the block diagram of the unity negative feedback closed loop
control system. Here, the single block is having the transfer function G(s) / [
1+G(s)H(s)−G(s)] instead of G(s).You can now calculate the steady state errors by using
steady state error formula given for the unity negative feedback systems.
Note − It is meaningless to find the steady state errors for unstable closed loop systems. So,
we have to calculate the steady state errors only for closed loop stable systems. This means
we need to check whether the control system is stable or not before finding the steady state
errors. In the next chapter, we will discuss the concepts-related stability.
The various types of controllers are used to improve the performance of control systems. In
this chapter, we will discuss the basic controllers such as the proportional, the derivative and
the integral controllers.
Proportional Controller
The proportional controller produces an output, which is proportional to error signal.
The block diagram of the unity negative feedback closed loop control system along with the
proportional controller is shown in the following figure.
DEPT. OF ME
Derivative Controller
The derivative controller produces an output, which is derivative of the error signal.
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.
DEPT. OF ME
The proportional derivative controller is used to improve the stability of control system
without affecting the steady state error.
The block diagram of the unity negative feedback closed loop control system along with the
proportional integral controller is shown in the following figure.
DEPT. OF ME
The proportional integral controller is used to decrease the steady state error without
affecting the stability of the control system.
The block diagram of the unity negative feedback closed loop control system along with the
proportional integral derivative controller is shown in the following figure.
DEPT. OF ME
UNIT - IV
STABILITY ANALYSIS IN S-DOMAIN
Stability is an important concept. In this chapter, let us discuss the stability of system and
types of systems based on stability.
Whatis Stability?
A system is said to be stable, if its output is under control. Otherwise, it is said to be unstable.
A stable system produces a bounded output for a given bounded input.
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 t including zero. So, it is bounded input. Therefore, the
first order control system is stable since both the input and the output are bounded.
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.
If all the roots of the characteristic equation exist to the left half of the ‘s’ plane, then the
control system is stable. If at least one root of the characteristic equation exists to the right
half of the ‘s’ plane, then the control system is unstable. So, we have to find the roots of the
characteristic equation to know whether the control system is stable or unstable. But, it is
difficult to find the roots of the characteristic equation as order increases.
So, to overcome this problem there we have the Routh array method. In this method, there
is no need to calculate the roots of the characteristic equation. First formulate the Routh
table and find the number of the sign changes in the first column of the Routh table. The
number of sign changes in the first column of the Routh table gives the number of roots of
characteristic equation that exist in the right half of the ‘s’ plane and the control system is
unstable.
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 sn and
continue up to the coefficient of s0.
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 an.
Here, an is the coefficient of s0 in 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.
DEPT. OF ME
Example:
Let us find the stability of the control system having characteristic equation,
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.
Example
DEPT. OF ME
Let us find the stability of the control system having characteristic equation,
The row s3 elements have 2 as the common factor. So, all these elements are divided by 2.
Special case (i) − Only the first element of row s2 is zero. So, replace it by ϵ and continue the
process of completing the Routh table.
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’s array are zero
In this case, follow these two steps −
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,
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.
DEPT. OF ME
DEPT. OF ME
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.
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
DEPT. OF ME
From above two cases, we can conclude that the root locus branches start at open loop poles
and end at open loop zeros.
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 –
The magnitude condition is that the point (which satisfied the angle condition) at which the
magnitude of the open loop transfers function is one.
The root locus is a graphical representation in s-domain and it is symmetrical about the real
axis. Because the open loop poles and zeros exist in the s-domain having the values either as
real or as complex conjugate pairs. In this chapter, let us discuss how to construct (draw) the
root locus.
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.
DEPT. OF ME
N=P if P≥Z
N=Z if P<Z
Rule 3 − Identify and draw the real axis root locus branches.
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=Z, then all the root locus branches start at finite open loop poles and end at finite
open loop zeros.
If P>Z, then Z number of root locus branches start at finite open loop poles and end
at finite open loop zeros and P−Z number of root locus branches start at finite open
loop poles and end at infinite open loop zeros.
If P<Z , then P number of root locus branches start at finite open loop poles and end
at finite open loop zeros and Z−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≠Z. Asymptotes give the
direction of these root locus branches. The intersection point of asymptotes on the real axis
is known as centroid.
Rule 5 − Find the intersection points of root locus branches with an imaginary axis.
DEPT. OF ME
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 K for 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.
Example
Let us now draw the root locus of the control system having open loop transfer
function,
Step 1 − The given open loop transfer function has three poles at s = 0,
s = -1, s = -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=3
The three poles are located are shown in the above figure. The line segment between s=−1,
and s=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 of s=−5.
Step 2 − We will get the values of the centroid and the angle of asymptotes by using the
given formulae.
Centroid
DEPT. OF ME
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 and
There will be one break-away point on the real axis root locus branch between the poles s
=−1 and s=0. By following the procedure given for the calculation of break-away point, we
will get it as s =−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.
DEPT. OF ME
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 increases and the time
domain specifications like delay time td, rise time tr and peak time tp 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 decreases and the time domain specifications like delay time td, rise
time tr and peak time tp increase.
So, based on the requirement, we can include (add) the open loop poles or zeros to the
transfer function.
DEPT. OF ME
DEPT. OF ME
DEPT. OF ME
Effects of adding a pole or a zero to the root locus of a second- order system
We discussed how we could change the value of gain K to change the position of the
closed-loop poles. This corresponds to placing a proportional gain, K, in cascade with the
system G(s) and finding the closed-loop poles for different values of gain, K. However,
proportional control is a simple form of control; it does not provide us with zero steady
example, in some control design problems, to produce the performance required in the
design specifications we need to move the poles to some positions on the s-plane, which
may not lie on a root locus defined by the simple proportional gain K. To be able to move
the poles to any position on the s-plane, we need to use a more complicated controller.
For example, we may need to add a zero or a pole to the controller and see how this will
affect the root locus and hence the position of the closed-loop poles. Examples of
controllers with poles or zeros are:
Thus, we need to know how the root locus will change if we add a pole or a zero.
To investigate this, we will use a simple example.
The poles are given by s = –p1 and s = –p2 and the simple root locus plot for this
system is shown in Figure 13.13(a). When we add a zero at s = –z1 to the controller,
the open-loop transfer function will change to:
DEPT. OF ME
DEPT. OF ME
UNIT-IV
Where,
Here, f0 is the frequency of the input sinusoidal signal. Similarly, you can follow the same
procedure for closed loop control system.
Frequency Domain Specifications
The frequency domain specifications are
Resonant peak
Resonant frequency
Bandwidth.
Consider the transfer function of the second order closed control system as
DEPT. OF ME
Resonant Peak
It is the peak (maximum) value of the magnitude of T(jω). It is denoted by Mr.
At u=ur, the Magnitude of T(jω) is -
DEPT. OF ME
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ω) drops to 70.7% from its zero
frequency value.
At ω=0, the value of u will be zero.
Substitute, u=0 in M.
Bandwidth ωb in the frequency response is inversely proportional to the rise time tr in the
time domain transient response.
Bode plots
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 plot is a horizontal line, which is independent of frequency. The 0 dB line
itself is the magnitude plot when the value of K is one. For the positive values of K, the
horizontal line will shift 20logK dB above the 0 dB line. For the negative values of K, the
horizontal line will shift 20logK dB below the 0 dB line. The Zero degrees line itself is the
phase plot for all the positive values of K.
DEPT. OF ME
The magnitude plot is a line, which is having a slope of 20 dB/dec. This line started
at ω=0.1rad/sec having a magnitude of -20 dB and it continues on the same slope. It is
touching 0 dB line at ω=1 rad/sec. In this case, the phase plot is 900 line.
Consider the open loop transfer function G(s)H(s)=1+sτ.
Magnitude
Phase angle
The magnitude plot is having magnitude of 0 dB upto ω=1τω=1τ rad/sec. From ω=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τ rad/sec and from here, it is having phase angle of 900. 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 table.
Represent the open loop transfer function in the standard time constant form.
DEPT. OF ME
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 syste
If the phase cross over frequency ωpc is greater than the gain cross over
frequency ωgc, then the control system is stable.
If the phase cross over frequency ωpc is equal to the gain cross over frequency ωgc,
then the control system is marginally stable.
If the phase cross over frequency ωpcis less than the gain crosses over frequency ω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)=20logMpc
Where, MpcMpc is the magnitude at phase cross over frequency. The unit of gain margin
(GM) is dB.
Phase Margin
The formula for phase margin PMPM is
PM=1800+ϕgc
Where, ϕgc is the phase angle at gain cross over frequency. The unit of phase margin
is degrees.
NOTE:
The stability of the control system based on the relation between gain margin and phase
margin is listed below.
If both the gain margin GM and the phase margin PM are positive, then the control
system is stable.
If both the gain margin GM and the phase margin PM are equal to zero, then the
control system is marginally stable.
If the gain margin GM and / or the phase margin PM are/is negative, then the control
system is unstable.
Polar plots
DEPT. OF ME
Polar plot is a plot which can be drawn between magnitude and phase. Here, the
magnitudes are represented by normal values only.
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 2700 in anti-clock wise
direction is equal to the angle −900 in clockwise direction.
Find the ending magnitude and the phase of G(jω)H(jω) by substituting ω=∞ So, the
polar plot ends with this magnitude and the phase angle.
Check whether the polar plot intersects the real axis, by making the imaginary term
of G(jω)H(jω) equal to zero and find the value(s) of ω.
Check whether the polar plot intersects the imaginary axis, by making real term
of G(jω)H(jω) equal to zero and find the value(s) of ω.
For drawing polar plot more clearly, find the magnitude and phase of G(jω)H(jω) by
considering the other value(s) of ω.
Example
Consider the open loop transfer function of a closed loop control system.
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.
DEPT. OF ME
By substituting ω=√2 in the magnitude of the open loop transfer function, we will
get M=0.83. Therefore, the polar plot intersects the negative real axis when ω=√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.
Nyquist Plots
Nyquist plots are the continuation of polar plots for finding the stability of the closed loop
control systems by varying ω from −∞ to ∞. That means, Nyquist plots are used to draw the
complete frequency response of the open loop transfer function.
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 Nyquist contour.
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
DEPT. OF ME
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) 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.
If the phase cross over frequency ωpc is greater than the gain cross over
frequency ωgc, then the control system is stable.
If the phase cross over frequency ωpc is equal to the gain cross over frequency ωgc,
then the control system is marginally stable.
If phase cross over frequency ωpc is less than gain cross over frequency ωgc, then the
control system is unstable.
Gain Margin
The gain margin GM is equal to the reciprocal of the magnitude of the Nyquist plot at the
phase cross over frequency.
Where, Mpc is the magnitude in normal scale at the phase cross over frequency.
Phase Margin
The phase margin PM is equal to the sum of 1800 and the phase angle at the gain cross over
frequency.
PM=1800+ϕgc
Where, ϕ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 GM is greater than one and the phase margin PM is positive, then
the control system is stable.
If the gain margin GMs equal to one and the phase margin PM is zero degrees, then
the control system is marginally stable.
If the gain margin GM is less than one and / or the phase margin PM is negative, then
the control system is unstable.
DEPT. OF ME
UNIT - III
STATE SPACE ANALYSIS OF CONTINUOUS SYSTEMS
The state space model of Linear Time-Invariant (LTI) system can be represented as,
X˙=AX+BU
Y=CX+DU
The first and the second equations are known as state equation and output equation respectively.
Where,
X and X˙ are the state vector and the differential state vector respectively.
U and Y are input vector and output vector respectively.
A is the system matrix.
B and C are the input and the output matrices.
D is the feed-forward matrix.
There are two storage elements (inductor and capacitor) in this circuit. So, the number of
the state variables is equal to two and these state variables are the current flowing through
the inductor, i(t) and the voltage across capacitor, vc(t).
From the circuit, the output voltage, v0(t) is equal to the voltage across capacitor, vc(t).
DEPT. OF ME
And u(t)=u
DEPT. OF ME
Then,
Here, D=[0].
DEPT. OF ME
Example:
Find the state space model for the system having transfer function.
and u(t)=u
Then, the state equation is
DEPT. OF ME
DEPT. OF ME
Y=CX+DU
⇒ (sI−A)X(s)=BU(s)
⇒ X(s) = (sI−A)−1BU(s)
⇒Y(s) = [C (sI−A)−1B+D]U(s)
The above equation represents the transfer function of the system. So, we can calculate the
transfer function of the system by using this formula for the system represented in the state
space model.
Example:
Let us calculate the transfer function of the system represented in the state space model as,
Therefore, the transfer function of the system for the given state space model is
From the above relation, we can write the state transition matrix ϕ(t) as
So, the zero input response can be obtained by multiplying the state transition
matrix ϕ(t) with the initial conditions matrix.
Properties of the state transition matrix
If t=0, then state transition matrix will be equal to an Identity matrix.
ϕ(0)=I
Inverse of state transition matrix will be same as that of state transition matrix just by
replacing‘t’ by ‘-t’.
ϕ(t1+t2)=ϕ(t1)ϕ(t2)
Controllability and Observability
Let us now discuss controllability and observability of control system one by one.
Controllability
A control system is said to be controllable if the initial states of the control system are
transferred (changed) to some other desired states by a controlled input in finite duration of
time.
Find the determinant of matrix QcQc and if it is not equal to zero, then the control
system is controllable.
DEPT. OF ME
Observability
A control system is said to be observable if it is able to determine the initial states of the
control system by observing the outputs in finite duration of time.
Find the determinant of matrix QoQo and if it is not equal to zero, then the control
system is observable.
Example:
Let us verify the controllability and observability of a control system which is represented in
the state space model as,
Since the determinant of matrix Qc is not equal to zero, the given control system is
controllable.
For n=2, the matrix Qo will be –
DEPT. OF ME
Since, the determinant of matrix Qo is not equal to zero, the given control system is
observable. Therefore, the given control system is both controllable and observable.