second order systems – Steady state errors – P, PI, PD and PID compensation – Analysis

using MATLAB.
UNIT III FREQUENCY RESPONSE ANALYSIS 9
Frequency response – Bode plot – Polar plot – Nyquist plot – Frequency domain
specifications from the plots – Constant M and N circles – ζichol‘s chart – Use of
ζichol‘s chart in control system analysis – Series – Parallel – Series-parallel
compensators – Lead – Lag – Lead lag compensators – Analysis using MATLAB.
UNIT IV STABILITY ANALYSIS 9
Stability – Routh-hurwitz criterion – Root locus technique – Construction of root locus –
Stability – Dominant poles – Application of root locus diagram – Nyquist stability
criterion – Relative stability – Analysis using MATLAB.
UNIT V STATE VARIABLE ANALYSIS & DIGITAL CONTROL
SYSTEMS 9
State space representation of continuous time systems – State equations – Transfer
function from state variable representation – Solutions of the state equations - Concepts
of controllability and observability – State space representation for discrete time systems
– Sampled data control systems – Sampling theorem – Sample and hold – Open loop and closed loop sampled data
systems.
Total : 45
TEXTBOOKS
1. Nagrath, J., and Gopal, ε., ―Control System Engineering‖, 5th Edition, ζew Age
International Publishers, 2007.
2. Gopal, ε., ―Control System Principles and Design‖, 2nd Edition, TεH, 2002.
REFERENCES
1. Benjamin C. Kuo., ―Automatic Control Systems‖, 7th Edition, PHI, 1995.
2. Schaum‘s Outline Series., ―Feedback and Control Systems‖, 2nd Edition, TεH,
2007.
3. John J. D‘azzo and Constantine H. Houpis., ―δinear Control System Analysis and
Design‖, 8th Edition, TεH, Inc., 1995.
4. Richard C. Dorf, and Robert H. Bishop., ―εodern Control Systems‖, 2nd Edition,
Addison Wesley, 1999.

CONTROL SYSTEMS

UNIT I

CONTROL SYSTEM MODELING

Basic elements of control system – Open loop and closed loop systems – Differential equation
– Transfer function – Modeling of electric systems – Translational and rotational mechanical Basic elements of control
systems – Block diagram reduction techniques – Signal flow graph.
system
In recent years, control systems have gained an increasingly importance in the development and advancement of
the modern civilization and technology. Figure shows the basic components of a control system. Disregard the
complexity of the system; it consists of an input (objective), the control system and its output (result). Practically our
day-to-day activities are affected by some type of control systems. There are two main branches of control systems:
1) Open-loop systems and
2) Closed-loop systems.

Open-loop systems:
The open-loop system is also called the non-feedback system. This is the simpler of the two systems. A simple
example is illustrated by the speed control of an automobile as shown in Figure 1-2. In this open-loop system, there is
no way to ensure the actual speed is close to the desired speed automatically. The actual speed might be way off the
desired speed because of the wind speed and/or road conditions, such as uphill or downhill etc.

Basic open-loop system

Closed-loop systems:
The closed-loop system is also called the feedback system. A simple closed-system is shown in Figure
1-3. It has a mechanism to ensure the actual speed is close to the desired speed automatically.

Transfer Function
A simpler system or element maybe governed by first order or second order differential equation

When several elements are connected in sequence, say ―n‖ elements, each one with first order, the total order of the
system will be nth order

In general, a collection of components or system shall be represented by nth order differential equation
 In control systems, transfer function characterizes the input output relationship of components or systems that can
be described by Liner Time Invariant Differential Equation

In the earlier period, the input output relationship of a device was represented graphically
In a system having two or more components in sequence, it is very difficult to find graphical relation
between the input of the first element and the output of the last element. This problem is solved by transfer function
Definition of Transfer Function
Transfer function of a LTIV system is defined as the ratio of the Laplace Transform of the output variable to the
Laplace Transform of the input variable assuming all the initial condition as zero
Properties of Transfer Function
The transfer function of a system is the mathematical model expressing the differential equation that relates the
output to input of the system
The transfer function is the property of a system independent of magnitude and the nature of the input
The transfer function includes the transfer functions of the individual elements. But at the same time, it does not
provide any information regarding physical structure of the system
The transfer functions of many physically different systems shall be identical
If the transfer function of the system is known, the output response can be studied for various types of inputs to
understand the nature of the system
If the transfer function is unknown, it may be found out experimentally by applying known inputs to the device and
studying the output of the system
How you can obtain the transfer function (T. F.)
Write the differential equation of the system
Take the L. T. of the differential equation, assuming all initial condition to be zero.
Take the ratio of the output to the input. This ratio is the T. F.
Mathematical Model of control systems
A control system is a collection of physical object connected together to serve an objective. The mathematical
model of a control system constitutes a set of differential equation.

1. Mechanical Translational systems

The model of mechanical translational systems can obtain by using three basic elements mass, spring and dash-
pot.When a force is applied to a translational mechanical system, it is opposed by opposing forces due to mass, friction and
elasticity of the system. The force acting on a mechanical body is governed by ζewton‘s second law of motion. For
translational systems it states that the sum of forces acting on a body is zero.
Force balance equations of idealized elements

Consider an ideal mass element shown in fig. which has negligible friction and elasticity. Let a force be applied on it.
The mass will offer an opposing force which is proportional to acceleration of a body.
 Let f = applied force
fm =opposing force due to mass

d 2x
Here fm M
dt 2

d 2x
By ζewton‘s second law, f = f m= M 2
dt
Consider an ideal frictional element dash-pot shown in fig. which has negligible mass and elasticity. Let a force be applied
on it. The dashpot will be offer an opposing force which is proportional to velocity of the body.

Let f = applied force

f b = opposing force due to friction
dx
Here, f b B
dt
dx
By ζewton‘s second law, f = f b B
dt
Consider an ideal elastic element spring shown in fig. which has negligible mass and friction.

Let f = applied force

f k = opposing force due to elasticity
Here, f k x
By ζewton‘s second law, f = f k x
Mechanical Rotational Systems

The model of rotational mechanical systems can be obtained by using three elements, moment of inertia [J] of mass,
dash pot with rotational frictional coefficient [B] and torsional spring with stiffness[k].
 When a torque is applied to a rotational mechanical system, it is opposed by opposing torques due to moment of inertia,
friction and elasticity of the system. The torque acting on rotational mechanical bodies is governed by ζewton‘s second
law of motion for rotational systems.

Torque balance equations of idealized elements

Consider an ideal mass element shown in fig. which has negligible friction and elasticity. The opposing torque due to
moment of inertia is proportional to the angular acceleration.

Let T = applied torque

Tj =opposing torque due to moment of inertia of the body
d2
Here Tj J
dt 2
By ζewton‘s law
d2
T= Tj = J
dt 2
Consider an ideal frictional element dash pot shown in fig. which has negligible moment of inertia and elasticity. Let a
torque be applied on it. The dash pot will offer an opposing torque is proportional to angular velocity of the body.

Let T = applied torque

Tb =opposing torque due to friction
d
Here Tb B ( 1 2 )
dt
By ζewton‘s law
d
T= Tb = B ( 1 2 )
dt
Consider an ideal elastic element, torsional spring as shown in fig. which has negligible moment of inertia and friction. Let
a torque be applied on it. The torsional spring will offer an opposing torque which is proportional to angular displacement
of the body.

Let T = applied torque

 Tk =opposing torque due to friction
Here Tk K( 1 2 )
By ζewton‘s law
T = Tk K( 1 2 )
Modeling of electrical system

Electrical circuits involving resistors, capacitors and inductors are considered. The behaviour of such systems is
governed by Ohm‘s law and Kirchhoff‘s laws
Resistor: Consider a resistance of ‗R‘ carrying current ‗i‘ Amps as shown in Fig (a), then the voltage drop across
it is v = R I

Inductor: Consider an inductor ―δ‘ H carrying current ‗i‘ Amps as shown in Fig (a), then the voltage drop across it
can be written as v = L di/dt

Capacitor: Consider a capacitor ―C‘ F carrying current ‗i‘ Amps as shown in Fig (a), then the voltage drop across it
can be written as v = (1/C) i dt

Steps for modeling of electrical system

Apply Kirchhoff‘s voltage law or Kirchhoff‘s current law to form the differential equations describing electrical
circuits comprising of resistors, capacitors, and inductors
Form Transfer Functions from the describing differential equations
Then simulate the model

Example
 t
1
R1i(t ) R2 i (t ) i (t )dt v1 (t )
c0
t
1
R2 i (t ) i (t )dt v2 (t )
c0

Electrical systems

LRC circuit. Applying Kirchhoff‘s voltage law to the system shown. We obtain the following equation;

Resistance circuit

Equation (1) & (2) give a mathematical model of the circuit. Taking the L.T. of equations (1)&(2), assuming zero initial
conditions, we obtain

Armature-Controlled dc motors

The dc motors have separately excited fields. They are either armature-controlled with fixed field or field-controlled with
fixed armature current. For example, dc motors used in instruments employ a fixed permanent-magnet field, and the
controlled signal is applied to the armature terminals.
Consider the armature-controlled dc motor shown in the following figure.
 Ra = armature-winding resistance, ohms
La = armature-winding inductance, henrys
ia = armature-winding current, amperes
if = field current, a-pares
ea = applied armature voltage, volt
eb = back emf, volts
θ = angular displacement of the motor shaft, radians
T = torque delivered by the motor, Newton*meter
J = equivalent moment of inertia of the motor and load referred to the motor shaft

kg.m2
f = equivalent viscous-friction coefficient of the motor and load referred to the motor shaft. Newton*m/rad/s
T = k1 i a ψ where ψ is the air gap flux, ψ = kf if , k1 is constant
For the constant flux

Where Kb is a back emf constant -------------- (1)

The differential equation for the armature circuit

The armature current produces the torque which is applied to the inertia and friction; hence

Assuming that all initial conditions are condition are zero/and taking the L.T. of equations (1), (2) & (3), we
obtain

Kpsθ(s) = Eb(s)

(Las+Ra) Ia(s) +Eb(s) =Ea(s) (Js2+fs)

θ (s) = T(s) = KIa(s)
The T.F can be obtained is

Analogous Systems
Let us consider a mechanical (both translational and rotational) and electrical system as shown in the fig.
From the fig (a)
d 2x dx
We get M D Kx f ----------------- (1)
dt 2 dt
From the fig (b)
d2 d
We get J 2 D K T -------------- (2)
dt dt
From the fig (c)
d 2q dq 1
We get L 2 R q v(t ) ------------ (3)
dt dt C

Where q i dt

They are two methods to get analogous system. These are (i) force- voltage (f-v) analogy and (ii) force-current (f-c)
analogy

Force –Voltage Analogy

Force – Current Analog

Problem
Find the system equation for system shown in the fig. And also determine f-v and f-i analogies.

For free body diagram M1

For free body diagram M2

Force –voltage analogy

From eq (1) we get

From eq (2) we get

From eq (3) and (4) we can draw f-v analogy

Force–current analogy

From eq (1) we get

From eq (2) we get

From eq (5) and (6) we can draw force-current analogy

The system can be represented in two forms:

Block diagram representation

Signal flow graph

Block diagram
A pictorial representation of the functions performed by each component and of the flow of signals.
Basic elements of a block diagram
Blocks
Transfer functions of elements inside the blocks
Summing points
Take off points
Arrow
Block diagram
A control system may consist of a number of components. A block diagram of a system is a pictorial representation of
the functions performed by each component and of the flow of signals. The elements of a block diagram are block, branch
point and summing point.
 Summing point Branch point

+ G(s)
―
-

Block
In a block diagram all system variables are linked to each other through functional blocks. The functional block or
simply block is a symbol for the mathematical operation on the input signal to the block that produces the output.

Summing point

Although blocks are used to identify many types of mathematical operations, operations of addition and
subtraction are represented by a circle, called a summing point. As shown in Figure a summing point may have one or
several inputs. Each input has its own appropriate plus or minus sign.
A summing point has only one output and is equal to the algebraic sum of the inputs.

A takeoff point is used to allow a signal to be used by more than one block or summing point.
The transfer function is given inside the block
• The input in this case is E(s)
• The output in this case is C(s)
• C(s) = G(s) E(s)
 Functional block – each element of the practical system represented by block with its T.F.
Branches – lines showing the connection between the blocks
Arrow – associated with each branch to indicate the direction of flow of signal
Closed loop system
Summing point – comparing the different signals
Take off point – point from which signal is taken for feed back
Advantages of Block Diagram Representation
Very simple to construct block diagram for a complicated system
Function of individual element can be visualized
Individual & Overall performance can be studied
Over all transfer function can be calculated easily
Disadvantages of Block Diagram Representation
No information about the physical construction
Source of energy is not shown

Simple or Canonical form of closed loop system

R(s) – Laplace of reference input r(t)

C(s) – Laplace of controlled output c(t)
E(s) – Laplace of error signal e(t)
B(s) – Laplace of feed back signal b(t)
G(s) – Forward path transfer function
H(s) – Feed back path transfer function

Block diagram reduction technique

Because of their simplicity and versatility, block diagrams are often used by control engineers to describe all
types of systems. A block diagram can be used simply to represent the composition and interconnection of a system.
Also, it can be used, together with transfer functions, to represent the cause-and-effect relationships throughout the
system. Transfer Function is defined as the relationship between an input signal and an output signal to a device.

Block diagram rules

Cascaded blocks

Moving a summer beyond the block

moving a summer ahead of block

Moving a pick-off ahead of block

Moving a pick-off behind a block

Eliminating a feedback loop

 Cascaded Subsystems

Parallel Subsystems

Feedback Control System

Procedure to solve Block Diagram Reduction Problems

Step 1: Reduce the blocks connected in series
Step 2: Reduce the blocks connected in parallel
Step 3: Reduce the minor feedback loops
Step 4: Try to shift take off points towards right and Summing point towards left
Step 5: Repeat steps 1 to 4 till simple form is obtained
Step 6: Obtain the Transfer Function of Overall System
Problem 1

Obtain the Transfer function of the given block diagram

Combine G1, G2 which are in series

Combine G3, G4 which are in Parallel

Reduce minor feedback loop of G1, G2 and H1

Transfer function

2. Obtain the transfer function for the system shown in the fig

3. Obtain the transfer function C/R for the block diagram shown in the fig

Solution
The take-off point is shifted after the block G2

Reducing the cascade block and parallel block

Replacing the internal feedback loop

Equivalent block diagram

Transfer function

Signal Flow Graph Representation

Signal Flow Graph Representation of a system obtained from the equations, which shows the flow of the signal

Signal flow graph

A signal flow graph is a diagram that represents a set of simultaneous linear algebraic equations. By taking Laplace
transfer, the time domain differential equations governing a control system can be transferred to a set of algebraic equation
in s-domain. A signal-flow graph consists of a network in which nodes are connected by directed branches. It depicts the
flow of signals from one point of a system to another and gives the relationships among the signals.

Basic Elements of a Signal flow graph

Node - a point representing a signal or variable.

Branch – unidirectional line segment joining two nodes.
Path – a branch or a continuous sequence of branches that can be traversed from one node to another node.
Loop – a closed path that originates and terminates on the same node and along the path no node is met twice.
Nontouching loops – two loops are said to be nontouching if they do not have a common node.

Mason’s gain formula

The relationship between an input variable and an output variable of signal flow graph is given by the net gain
between the input and the output nodes is known as overall gain of the system. εason‘s gain rule for the determination of
the overall system gain is given below.

Where M= gain between Xin and Xout

Xout =output node variable
Xin= input node variable
N = total number of forward paths
Pk= path gain of the kth forward path
=1-(sum of loop gains of all individual loop) + (sum of gain product of all possible combinations of two nontouching
loops) – (sum of gain products of all possible combination of three nontouching loops)
Problem
 UNIT I
CONTROL SYSTEM MODELLING
PART – A

1. What is control system?

2. Define open loop control system.
3. Define closed loop control system.
4. Define transfer function.
5. What are the basic elements used for modeling mechanical rotational system?
6. Name two types of electrical analogous for mechanical system.
7. What is block diagram?
8. What is the basis for framing the rules of block diagram reduction technique?
9. What is a signal flow graph?
10. What is transmittance?
11. What is sink and source?
12. Define non- touching loop.
13. Write Masons Gain formula.
14. Write the analogous electrical elements in force voltage analogy for the elements of
mechanical translational system.
15. Write the analogous electrical elements in force current analogy for the elements of
mechanical translational system.
16. Write the force balance equation of m ideal mass element.
17. Write the force balance equation of ideal dashpot element.
18. Write the force balance equation of ideal spring element.
19. What is servomechanism?
PART – B
1. Write the differential equations governing the Mechanical system shown in fig.
and determine the transfer function. (16)

2. Determine the transfer function Y2(S)/F(S) of the system shown in fig. (16)

3. Write the differential equations governing the Mechanical rotational system shown in
fig. Draw the Torque-voltage and Torque-current electrical analogous circuits. (16)

4. Determine the overall transfer function C(S)/R(S) for the system shown in fig. (16)
5. For the system represented by the block diagram shown in fig. Determine C1/R1 and
C2/R1. (16)

6. Find the overall gain of the system whose signal flow graph is shown in fig. (16)

7. Draw a signal flow graph and evaluate the closed loop transfer function of a system
Whose block is shown in fig?
8. Derive the transfer function for Armature controlled DC servo motor. (16)
9. Derive the transfer function for Field controlled DC servo motor. (16)

UNIT II
TIME RESPONSE ANALYSIS

Time response analysis – First order systems – Impulse and step response analysis of second order systems – Steady
state errors – P, PI, PD and PID compensation – Analysis using MATLAB

Time response analysis

Introduction

After deriving a mathematical model of a system, the system performance analysis can be done in various methods.
In analyzing and designing control systems, a basis of comparison of performance of various control systems should
be made. This basis may be set up by specifying particular test input signals and by comparing the responses of
various systems to these signals.
The system stability, system accuracy and complete evaluation are always based on the time response analysis and
the corresponding results

Next important step after a mathematical model of a system is obtained.

To analyze the system‘s performance.
Normally use the standard input signals to identify the characteristics of system‘s response
o Step function
 o Ramp function
o Impulse function
o Parabolic function
o Sinusoidal function
Time response analysis

It is an equation or a plot that describes the behavior of a system and contains much information about it with respect
to time response specification as overshooting, settling time, peak time, rise time and steady state error. Time response is
formed by the transient response and the
steady Time response = Transient response + Steady state response state response.

Transient time response (Natural response) describes the behavior of the system in its first short time until arrives the
steady state value and this response will be our study focus. If the input is step function then the output or the response
is called step time response and if the input is ramp, the response is called ramp time response ... etc.
Classification of Time Response
Transient response
Steady state response
y(t) = yt(t) + yss(t)
Transient Response
The transient response is defined as the part of the time response that goes to zero as time becomes very large. Thus
yt(t) has the property
Lim yt(t) = 0
t -->
The time required to achieve the final value is called transient period. The transient response may be exponential or
oscillatory in nature. Output response consists of the sum of forced response (form the input) and natural response (from
the nature of the system).The transient response is the change in output response from the beginning of the response to the
final state of the response and the steady state response is the output response as time is approaching infinity (or no more
changes at the output).

Steady State Response

The steady state response is the part of the total response that remains after the transient has died out. For a position
control system, the steady state response when compared to with the desired reference position gives an indication of the
final accuracy of the system. If the steady state response of the output does not agree with the desired reference exactly, the
system is said to have steady state error.
Typical Input Signals
Impulse Signal
Step Signal
Ramp Signal
Parabolic Signal

Time Response Analysis & Design

Two types of inputs can be applied to a control system
Command Input or Reference Input yr(t)
Disturbance Input w(t) (External disturbances w(t) are typically uncontrolled variations in the load on a control
system)
In systems controlling mechanical motions, load disturbances may represent forces.
In voltage regulating systems, variations in electrical load area major source of disturbances.
Test Signals
Input r(t) R(s)

Step Input A A/s

2
Ramp Input At
A/s
 2 3
Parabolic Input
At / 2 A/s

Impulse Input (t) 1

Transfer Function
One of the types of Modeling a system
Using first principle, differential equation is obtained
Laplace Transform is applied to the equation assuming zero initial conditions
Ratio of LT (output) to LT (input) is expressed as a ratio of polynomial in s in the transfer function
Order of a system
The Order of a system is given by the order of the differential equation governing the system
Alternatively, order can be obtained from the transfer function
In the transfer function, the maximum power of s in the denominator polynomial gives the order of the system
Dynamic Order of Systems
Order of the system is the order of the differential equation that governs the dynamic behaviour
Working interpretation: Number of the dynamic elements / capacitances or holdup elements between a
manipulated variable and a controlled variable
Higher order system responses are usually very difficult to resolve from one another
The response generally becomes sluggish as the order increases

System Response
First-order system time response
 Transient
 Steady-state

Second-order system time response

 Transient
 Steady-state
First Order System
Y ( s) K K
R( s) 1 K sT 1 sT

Step Response of First Order System

Evolution of the transient response is determined by the pole of the transfer function at s=-1/t where t is the time
constant
Also, the step response can be found:
 (s 1 ) C (s) K
s
K K K
C (s)
s (s 1 ) s s 1
t/
c(t ) K (1 e )u (t )

Impulse response Exponential

K
1 sT

Step response K K Step, exponential

-
s s 1/ T

Ramp response K KT KT Ramp, step, exponential

2
- -
s s s 1/ T

Second-order systems

LTI second-order system

Second-Order Systems
Second order system responses

Overdamped response:
Poles: Two real at

1 , 2
Natural response: Two exponentials with time constants equal to the reciprocal of the pole location

1t 2t
c(t ) k1e k2e
Poles: Two complex at

Underdamped response:

d j d

ζatural response: Damped sinusoid with an exponential envelope whose time constant is equal to the reciprocal of the pole‘s
radian frequency of the sinusoid, the damped frequency of oscillation, is equal to the imaginary part of the poles

dt
c(t ) Ae cos( d t )
Undamped response:

Poles: Two imaginary at

j 1
Natural response: Undamped sinusoid with radian frequency equal to the imaginary part of the poles

c(t ) A cos( 1t )
Critically damped responses:

Poles: Two real at

Natural response: One term is an exponential whose time constant is equal to the reciprocal of the pole location. Another ter
product of time and an exponential with time constant equal to the reciprocal of the pole location

1t 1t
c(t ) k1e k2te

Second order system responses damping cases

Second- order step response

complex poles

Steady State Error

Consider a unity feedback system

Transfer function between e(t) and r(t)

E(s) 1 R(s)
or E(s)
R(s) 1 G(s) 1 G(s)

Steady state error is

sR(s)
ess lime(t) limsE(s) lim1
t s 0 s 0 G(s)

Steady State Error

Output Feedback Control Systems

Feedback only the output signal

– Easy access
– Obtainable in practice
PID Controllers

Proportional controllers
– pure gain or attenuation
Integral controllers
– integrate error
Derivative controllers
– differentiate error
Proportional Controller

Controller input is error (reference output)

Controller output is control signal
P controller involves only a proportional gain (or attenuation)

Integral Controller
Integral of error with a constant gain
Increase system type by 1
Infinity steady-state gain
Eliminate steady-state error for a unit step input
Integral Controller

Derivative Control

Differentiation of error with a constant gain

Reduce overshoot and oscillation
Do not affect steady-state response
Sensitive to noise
Controller Structure
Single controller
P controller, I controller, D controller
Combination of controllers
PI controller, PD controller
PID controller
Controller Performance
P controller
PI controller
PD controller
PID controller
Design of PID Controllers
Based on the knowledge of P, I and D
– trial and error
– manual tuning
– simulation
Design of PID Controllers
Time response measurements are particularly simple.
A step input to a system is simply a suddenly applied input - often just a constant voltage applied through a switch.
The system output is usually a voltage, or a voltage output from a transducer measuring the output.
A voltage output can usually be captured in a file using a C program or a Visual Basic program.
 You can use responses in the time domain to help you determine the transfer function of a system.
First we will examine a simple situation. Here is the step response of a system. This is an example of really "clean"
data, better than you might have from measurements. The input to the system is a step of height 0.4. The goal is to
determine the transfer function of the system.
Impulse Response of A First Order System
The impulse response of a system is an important response. The impulse response is the response to a unit impulse.
The unit impulse has a Laplace transform of unity (1).That gives the unit impulse a unique stature. If a system has a unit
impulse input, the output transform is G(s), where G(s) is the transfer function of the system. The unit impulse response is
therefore the inverse transform of G(s), i.e. g(t), the time function you get by inverse transforming G(s). If you haven't
begun to study Laplace transforms yet, you can just file these last statements away until you begin to learn about Laplace
transforms. Still there is an important fact buried in all of this.
Knowing that the impulse response is the inverse transform of the transfer function of a system can be useful in
identifying systems (getting system parameters from measured responses).
In this section we will examine the shapes/forms of several impulse responses. We will start with simple first
order systems, and give you links to modules that discuss other, higher order responses.
A general first order system satisfies a differential equation with this general form.

If the input, u(t), is a unit impulse, then for a short instant around t = 0 the input is infinite.
Let us assume that the state, x(t), is initially zero, i.e. x(0) = 0. We will integrate both sides of the differential equation
from a small time, ฀, before t = 0, to a small time, after t = 0. We are just taking advantage of one of the properties of the
unit impulse.
The right hand side of the equation is just Gdc since the impulse is assumed to be a unit impulse - one with unit area.
Thus, we have:

We can also note that x(0) = 0, so the second integral on the right hand side is zero. In other words, what the impulse
does is it produces a calculable change in the state, x(t), and this change occurs in a negligibly short time (the duration of
the impulse) after t = 0 That leads us to a simple strategy for getting the impulse response. Calculate the new initial
condition after the impulse passes. Solve the differential equation - with zero input - starting from the newly calculated
initial condition.
UNIT II
TIME RESPONSE ANALYSIS
PART-A
1. What is Proportional controller and what are its advantages?
2. What is the drawback in P-controller?
3. What is integral control action?
4. What is the advantage and disadvantage in integral controller?
5. What is PI controller?
6. What is PD controller?
7. What is PID controller?
8. What is time response?
9. What is transient and steady state response?
10.Name the test signals used in control system.
11. Define Step signal:
12. Define Ramp signal:
13. Define parabolic signal:
14. What is an impulse signal?
15. What is the order of a system?
16. Define Damping ratio.
17. Give the expression for damping ratio of mechanical and electrical system.
18. How the system is classified depending on the value of damping?
29. What will be the nature of response of a second order system with different types of
damping?
20. Sketch the response of a second order under damped system.
21. What is damped frequency of oscillation?
22. List the time domain specifications:
24. Define Delay time.
25. Define rise time.
26. Define Peak time.
27. Define Peak overshoot.
28. Define settling time.
29. What is type number of a system? What is its significance?
30. Distinguish between type and order of a system:
31. What is steady state error?
32. Define acceleration error constant:
33. What are generalized error coefficients?
34. Give the relation between generalized and static error coefficients:
35. Mention two advantages of generalized error constants over static error constants
PART-B
1. (a) Derive the expressions and draw the response of first order system for unit step input. (8)
(b) Draw the response of second order system for critically damp case and when input is unit step. (8)
2. Derive the expressions for Rise time, Peak time, and Peak overshoot. (16)
3. A potential control system with velocity feedback is shown in fig. What is response of the system for unit step input?

4. Measurements conducted on a Servomechanism show the system response to be c(t)=1+0.2 ê 60t-1.2 ê –10 t. when
subjected to a unit step. Obtain an expression for closed loop transfer function. (16)
5. A positional control system with velocity feedback is shown in fig. What is the response c(t) to the unit step input. Given
that =0.5.and also calculate rise time, peak time, Maximum overshoot and settling time. (16)
6. A unity feedback control system has an open loop transfer functionG(S) = 10/S(S+2). Find the rise time, percentage
over shoot, peak time and settling time(16)
7. A closed loop servo is represented by the differential equation, where c is the displacement of the output shaft, r is the
displacement of the input shaft and e= r-c. Determine undamped natural frequency, damping ratio and percentage
maximum overshoot for unit step input. (16)
8. For a unity feedback control system the open loop transfer function G(S) = 10(S+2)/ S2 (S+1). Find
(a) Position, velocity and acceleration error constants.
(b) The steady state error when the input is R(S) where R(S) =3/S –2/S2 +1/3S3 (16)
9. The open loop transfer function of a servo system with unity feedback system is G(S) = 10/ S(0.1S+1).Evaluate the static
error constants of the system. Obtain the steady state error of the system when subjected to an input given by the
polynomial r(t) = a0+a1t +a2 /2 t2 . (16)

Unit III
Frequency response analysis

Frequency Response - Bode Plot, Polar Plot, Nyquist Plot - Frequency Domain Frequency Response
specifications from the plots - Constant M and N Circles - ζichol‘s Chart - Use of
The frequency
ζichol‘s Chart in Control System Analysis. Series, Parallel, series-parallel Compensators
response of a system is a
- Lead, Lag, and Lead Lag Compensators, Analysis using MATLAB.
frequency dependent
function which expresses
how a sinusoidal signal of a given frequency on the system input is transferred through the system. Time-varying signals at
least periodical signals —which excite systems, as the reference (set point) signal or a disturbance in a control system or
measurement signals which are inputs signals to signal filters, can be regarded as consisting of a sum of frequency
components. Each frequency component is a sinusoidal signal having certain amplitude and a certain frequency. (The
Fourier series expansion or the Fourier transform can be used to express these frequency components quantitatively.) The
frequency response expresses how each of these frequency components is transferred through the system. Some
components may be amplified, others may be attenuated, and there will be some phase lag through the system.
 The frequency response is an important tool for analysis and design of signal filters (as low pass filters and high
pass filters), and for analysis, and to some extent, design, of control systems. Both signal filtering and control systems
applications are described (briefly) later in this chapter. The definition of the frequency response — which will be given in
the next section — applies only to linear models, but this linear model may very well be the local linear model about some
operating point of a non-linear model. The frequency response can found experimentally or from a transfer function model.
It can be presented graphically or as a mathematical function.

Bode plot

• Plots of the magnitude and phase characteristics are used to fully describe the frequency response

• A Bode plot is a (semilog) plot of the transfer function magnitude and phase angle as a function of frequency

• The gain magnitude is many times expressed in terms of decibels (dB)

dB = 20 log10 A

Bode plot procedure

There are 4 basic forms in an open-loop transfer function G(jω)H(jω)

Gain factor, K
(jω)±p factor: pole and zero at origin
(1+jωT)±q factor
Quadratic factor

r
2
1 j2 2
n n

Gain margin and Phase margin

Gain margin:

The gain margin is the number of dB that is below 0 dB at the phase crossover frequency (ø=-180º). It can also be
increased before the closed-
Term Corner Frequency Slope Change
loop system db/dB in slope becomes
unstable
20/jw --- -20
Phase margin: 1/(1+j4w) wc1=1/4=0.25 -20 -20-20=-40

The phase margin

is the number of degrees the
phase of that is above -180º at the gain crossover frequency

Gain margin and Phase margin

Bode Plot – Example

For the following T.F draw the Bode plot and obtain Gain cross over frequency (w gc) , Phase cross over frequency ,
Gain Margin and Phase Margin.

G(s) = 20 / [s (1+3s) (1+4s)]

Solution:

The sinusoidal T.F of G(s) is obtained by replacing s by jw in the given T.F

G(jw) = 20 / [jw (1+j3w) (1+j4w)]

Corner frequencies: wc1= 1/4 = 0.25 rad /sec ;

wc2 = 1/3 = 0.33 rad /sec

 1/(1+j3w) wc2=1/3=0.33 -20 -40-20=-60

Choose a lower corner frequency and a higher Corner frequency

wl= 0.025 rad/sec ; wh = 3.3 rad / sec

Calculation of Gain (A) (MAGNITUDE PLOT)

A @ wl ; A= 20 log [ 20 / 0.025 ] = 58 .06 dB

A @ wc1 ; A = [Slope from wl to wc1 x log (wc1 / wl ] + Gain (A)@wl

= - 20 log [ 0.25 / 0.025 ] + 58.06

= 38.06 dB

A @ wc2 ; A = [Slope from wc1 to wc2 x log (wc2 / wc1 ] + Gain (A)@ wc1

= - 40 log [ 0.33 / 0.25 ] + 38

= 33 dB

A @ wh ; A = [Slope from wc2 to wh x log (wh / wc2 ] + Gain (A) @ wc2

= - 60 log [ 3.3 / 0.33 ] + 33

= - 27 dB

Calculation of Phase angle for different values of frequencies [PHASE PLOT]

Ø = -90O- tan -1 3w – tan -1 4w

When

Frequency in rad / Phase Angle in degrees

sec
o
w=0 Ø= -90
o
w=0.025 Ø= - 99
o
w=0.25 Ø= -172
o
w=0.33 Ø= -188
o
w=3.3 Ø= - 259
w=∞ Ø= - 270
o
 Calculations of Gain cross over frequency
The frequency at which the dB magnitude is Zero wgc = 1.1 rad / sec
Calculations of Phase cross over frequency
The frequency at which the Phase of the system is - 180o
wpc = 0.3 rad / sec
Gain Margin
The gain margin in dB is given by the negative of dB magnitude of G(jw) at phase cross over frequency
GM = - { 20 log [G( jwpc )] = - { 32 } = -32 dB
Phase Margin
Ґ = 180o+ Øgc= 180o + (- 240o) = -60o
Conclusion
For this system GM and PM are negative in Values. Therefore the system is unstable in nature.
Polar plot

To sketch the polar plot of G(jω) for the entire range of frequency ω, i.e., from 0 to infinity, there are four
key points that usually need to be known:
(1) the start of plot where ω = 0,
(2) the end of plot where ω = ∞,
(3) where the plot crosses the real axis, i.e., Im(G(jω)) = 0, and
(4) where the plot crosses the imaginary axis, i.e., Re(G(jω)) = 0.

BASICS OF POLAR PLOT

The polar plot of a sinusoidal transfer function G(jω) is a plot of the magnitude of G(jω) Vs the phase of

G(jω) on polar co-ordinates as ω is varied from 0 to ∞.

(ie) |G(jω)| Vs G(jω) as ω → 0 to ∞.

Polar graph sheet has concentric circles and radial lines.

Concentric circles represents the magnitude.

 Radial lines represents the phase angles.

In polar sheet

+ve phase angle is measured in ACW from 0

-ve phase angle is measured in CW from 0

PROCEDURE

Express the given expression of OLTF in (1+sT) form.

Substitute s = jω in the expression for G(s)H(s) and get G(jω)H(jω).

Get the expressions for | G(jω)H(jω)| & G(jω)H(jω).

Tabulate various values of magnitude and phase angles for different values of ω ranging from 0 to ∞.

Usually the choice of frequencies will be the corner frequency and around corner frequencies.

Choose proper scale for the magnitude circles.

Fix all the points in the polar graph sheet and join the points by a smooth curve.

Write the frequency corresponding to each of the point of the plot.

MINIMUM PHASE SYSTEMS

Systems with all poles & zeros in the Left half of the s-plane – Minimum Phase Systems.

For Minimum Phase Systems with only poles

Type No. determines at what quadrant the polar plot starts.

Order determines at what quadrant the polar plot ends.

Type No. → No. of poles lying at the origin

Order → Max power of s in the denominator polynomial of the transfer function.

GAIN MARGIN

Gain Margin is defined as the factor by which the system gain can be increased to drive the system to the

verge of instability .

For stable systems,

ωgc ωpc

G(j )H(j ) at ω=ωpc 1

GM = in positive dB
More positive the GM, more stable is the system.

For marginally stable systems,

ωgc = ωpc

G(j )H(j ) at ω=ωpc = 1

GM = 0 dB

For Unstable systems,

ωgc ωpc

G(j )H(j ) at ω=ωpc > 1

GM = in negative dB

Gain is to be reduced to make the system stable

Note:

If the gain is high, the GM is low and the system s step response shows high overshoots and long settling

time.

On the contrary, very low gains give high GM and PM, but also causes higher ess, higher values of rise

time and settling time and in general give sluggish response.

Thus we should keep the gain as high as possible to reduce ess and obtain acceptable response speed

and yet maintain adequate GM and PM.

An adequate GM of 2 ( ie 6 dB ) and a PM of about 30 is generally considered good enough as a thumb

rule.

At = pc , G(j )H(j ) = -180

Let G(j )H(j ) at = pc be taken as B

If the gain of the system is increased by a factor 1/B, then the G(j )H(j ) at = pc becomes B(1/B)

= 1 and hence the G(j )H(j ) locus pass through -1+j0 point driving the system to the verge of instability.

GM is defined as the reciprocal of the magnitude of the OLTF evaluated at the phase cross over

frequency.

GM = Kg = 1 / G(j )H(j ) = pc

GM in dB = 20 log (1/B) = - 20 log B

PHASE MARGIN

Phase Margin is defined as the additional phase lag that can be introduced before the system becomes

unstable .

Let A be the point of intersection of G(j )H(j ) plot and a unit circle centered at the origin.

Draw a line connecting the points O & A and measure the phase angle between the line OA and

+ve real axis.

This angle is the phase angle of the system at the gain cross over frequency.

G(j gc)H(j gc) = gc

If an additional phase lag of PM is introduced at this frequency, then the phase angle G(j gc)H(j gc) will
become 180 and the point ‗A‘ coincides with (-1+j0) driving the system to the verge of instability.
This additional phase lag is known as the Phase Margin.
= 180 + G(j gc)H(j gc)

= 180 + gc

[Since gc is measured in CW direction, it is taken as negative]

For a stable system, the phase margin is positive.
A Phase margin close to zero corresponds to highly oscillatory system.

A polar plot may be constructed from experimental data or from a system transfer function
If values of w are marked along the contour, a polar plot has the same information as a Bode plot
Usually, the shape of a polar plot is of most interest

Nyquist Plot

The Nyquist plot is a polar plot of the function

The Nyquist stability criterion relates the location of the roots of the characteristic equation to the open-loop frequency
response of the system. In this, the computation of closed-loop poles is not necessary to determine the stability of the
system and the stability study can be carried out graphically from the open-loop frequency response. Therefore.
experimentally determined open-loop frequency response can be used directly for the study of stability. When the feedback
path is closed. The Nyquist criterion has the following features that make it an alternative method that is attractive for the
analysis and design of control systems. 1. In addition to providing information on absolute and relative.

Nyquist Plot Example

Consider the following transfer function

Cha ge it fro s do ai to j do ai :

Find the magnitude and phase angle equations:

E aluate ag itude a d phase a gle at ω= 0+ a d ω = +∞

Draw the nyquist plot:

Frequency domain specifications

The resonant peak Mr is the maximum value of jM(jw)j.

_ The resonant frequency !r is the frequency at which the peak resonance Mr occurs.
_ The bandwidth BW is the frequency at which(jw) drops to 70:7% (3 dB) of its zero-frequency value.

Mr indicates the relative stability of a stable closed loop system.

A largeMr corresponds to larger maximum overshoot of the step response.
Desirable value: 1.1 to 1.5
BW gives an indication of the transient response properties of a control system.
A large bandwidth corresponds to a faster rise time. BW and rise time tr are inversely proportional.
BW also indicates the noise-filtering characteristics and robustness of the system.
Increasing wn increases BW.
BW and Mr are proportional to each other
Constant M and N circles
Consider a candidate design of a loop transfer function L( jω) shown on the RHS.
Evaluate T( jω) from L( jω) in the manner of frequency point by frequency point.
Alternatively, the Bode plot of L( jω) can also be show on the complex plane to form its Nyquist plot.

M circles (constant magnitude of T)

In order to precisely evaluate |T( jω)| from the Nyquist plot of L( jω), a tool called M circle is developed as followed.
Let L( jω)=X+jY, where X is the real and Y the imaginary part. Then

Rearranging the above equation, it gives

X2(1−M2)−2M2X−M2+(1−M2)Y2=0
That is, all (X, Y) pair corresponding to a constant value of M for a circle on the complex plane. Therefore, we have the
following (constant) M circles on the complex plane as shown below.
N circles (constant phase of T)
Similarly, it can be shown that the phase of T( jω) be

It can be shown that all (X, Y) pair which corresponds to the same constant phase of T (i.e., constant N) forms a circle on
the complex plane as shown below.

Example
Nyquist plot of L( jω), and M-N circles of T( jω)
Nichols Chart

The Nyquist plot of L( jω) can also be represented by its polar form using dB as magnitude and degree as phase .

All L( jω) which corresponds to a constant |T( jω)| can be draw as a locus of M circle on this plane as shown below
And ll L( jω) which corresponds to a constant α( jω) can be draw as a locus of M circle on this plane as shown below

Combining the above two groups of M circles and N circles, we have the Nichols chart as shown below
Use of Nichols chat in control system analysis
And the resultant close-loop system has the following Bode plot
 TYPES OF COMPENSATION
Series Compensation or Cascade Compensation
This is the most commonly used system where the controller is placed in series with the controlled process.
Figure shows the series compensation.

Feedback Compensation or Parallel Compensation

This is the system where the controller is placed in the senor feedback path as shown in Fig.

State Feedback Compensation

This is a system which generates the control signal by feeding back the state variables through constant real gains.
The scheme is termed state feedback. It is shown in Fig.
The compensation schemes shown in Figs above have one degree of freedom, since there is only one controller in each
system. The demerit with one degree of freedom controllers is that the performance criteria that can be realized are limited.
That is why there are compensation schemes which have two degree freedoms, such as:
(a) Series-feedback compensation
(b) Feedforward compensation
Series-Feedback Compensation
Series-feedback compensation is the scheme for which a series controller and a feedback controller are used. Figure 9.6
shows the series-feedback compensation scheme.

Feedforward Compensation
The feedforward controller is placed in series with the closed-loop system which has a controller in the forward path
Orig. 9.71. In Fig. 9.8, Feedforward the is placed in parallel with the controller in the forward path. The commonly used
controllers in the above-mentioned compensation schemes are now described in the section below.

Lead Compensator
It has a zero and a pole with zero closer to the origin. The general form of the transfer function of the load compensator is

Subsisting

Transfer function
 Lag Compensator
It has a zero and a pole with the zero situated on the left of the pole on the negative real axis. The general form of the
transfer function of the lag compensator is

Therefore, the frequency response of the above transfer function will be

Now comparing with

 Therefore

Lag-Lead Compensator
The lag-lead compensator is the combination of a lag compensator and a lead compensator. The lag-section is
provided with one real pule and one real zero, the pole being to the right of zero, whereas the lead section has one real polo
and one real cam with the zero being to the right of the pole.
The transfer function of the lag-lead compensator will be

The figure shows lag lead compensator

The above transfer functions are comparing with

Then
 Therefore

UNIT III
FREQUENCY RESPONSE ANALYSIS

PART -A
1. What is frequency response?
2. What are advantages of frequency response analysis?
3. What are frequency domain specifications?
4. Define Resonant Peak.
5. What is resonant frequency?
6. Define Bandwidth.
7. What is cut-off rate?
8. Define gain margin.
9. Define phase margin.
10. What is phase and Gain cross-over frequency?
11. What is Bode plot?
12. Define corner frequency.
13. What are the advantages of Bode Plot?
14. What is a Nichols plot?
15. What are M and N circles?
16. What is Nichols chart?
17. What are the advantages of Nichols chart?
18. What is polar plot?
19. What is minimum phase system?
20. What are All-Pass systems?
PART-B
1. Plot 1. Plot the Bode diagram for the following transfer function and obtain the gain and phase cross over frequencies:
G(S) = 10/ S(1+0.4S) (1+0.1S) (16)
2. The open loop transfer function of a unity feedback system is G(S) = 1/ S (1+S) (1+2S) Sketch the Polar plot and
determine the Gain margin and Phase margin (16)
3. Sketch the Bode plot and hence find Gain cross over frequency, Phase cross over frequency, Gain margin and Phase
margin. G(S) = 0.75(1+0.2S)/ S(1+0.5S) (1+0.1S)(16)
4. Sketch the Bode plot and hence find Gain cross over frequency, Phase cross over frequency, Gain margin and Phase
margin. G(S) = 10(S+3) / S(S+2) (S2+4S+100) (16)
5. Sketch the polar plot for the following transfer function and find Gain cross over frequencies, Phase cross over
frequency, Gain margin and Phase margi G(S) = 10(S+2) (S+4)/ S (S2 -3S+10) (16)
6. Construct the polar plot for the function GH(S) =2(S+1)/ S2. Find Gain cross over frequency, Phase cross over
frequency, Gain margin and Phase margin. (16)
7. Plot the Bode diagram for the following transfer function and obtain the gain and phase cross over frequencies. G(S)
=KS2 / (1+0.2S) (1+0.02S). Determine the value of K for a gain cross over frequency of 20 rad/sec. (16)
8. Sketch the polar plot for the following transfer function and find Gain cross over frequency, Phase cross over frequency,
Gain margin and Phase margin.G(S) = 400/ S (S+2) (S+10) (16)
9. A unity feedback system has open loop transfer function G(S) = 20/ S(S+2) (S+5).Using ζichol‘s chart determine the
closed loop frequency Response and estimate all thefrequency domain specifications. (16)
10. Sketch the Bode plot and hence find Gain cross over frequency, Phase cross over frequency, Gain margin and Phase
margin: G(S) = 10(1+0.1S)/ S (1+0.01S) (1+S). (16)
11. Write short notes on correlation between the time and frequency response? (16)
12. What is compensation? Why it is needed for control system? Explain the types of compensation (16)
13. Realize the basic compensators using electrical network and obtain the transfer function. (16)
14. Design a suitable lead compensators for a system with unity feedback and having open loop transfer function G(S)= K/
S(S+1) (S+4) to meet the specifications.
(i)Damping ratio=0.5
(ii) Undamped natural frequency _n =2 rad/sec. (16)
15. A unity feedback system has an open loop transfer function G(S) = K/ S(S+1) (0.2S+1).Design a suitable phase lag
compensators to achieve the following specifications Kv= 8 and Phase margin 40 deg with usual notation. (16)
16. Explain the procedure for lead compensation and lag compensation. (16)
17. Explain the design procedure for lag-lead compensation. (16)
18. Consider a type 1 unity feedback system with an OLTF G(S) =K/S (S+1) (S+4).
The system is to be compensated to meet the following specifications Kv > 5sec
and PM>43 deg. Design suitable lag compensators.
UNIT IV
STABILITY ANALYSIS

Stability – Routh-hurwitz criterion – Root locus technique – Construction of root locus –Stability – Dominant poles –
Application of root locus diagram – Nyquist stability Criterion – Relative stability – Analysis using MATLAB.

Stability
A system is stable if any bounded input produces a bounded output for all bounded initial conditions.

Basic concept of stability

Stability of the system and roots of characteristic equations

Characteristic Equation
Consider an nth-order system whose the characteristic equation (which is also the denominator of the transfer function) is

a( s) sn a1s n 1
a2 s n 2
 a n 1 s 1 a0 s 0

Routh Hurwitz Criterion

Goal: Determining whether the system is stable or unstable from a characteristic equation in polynomial form
without actually solving for the roots
Routh’s stability criterion is useful for determining the ranges of coefficients of polynomials for stability, especially
when the coefficients are in symbolic (non numerical) form
To find Kmar & ω

A necessary condition for Routh’s Stability

A necessary condition for stability of the system is that all of the roots of its characteristic equation have negative
real parts, which in turn requires that all the coefficients be positive.
A necessary (but not sufficient) condition for stability is that all the coefficients of the polynomial characteristic
equation are positive & none of the co-efficient vanishes.
Routh’s formulation requires the computation of a triangular array that is a function of the coefficients of the
polynomial characteristic equation.
A system is stable if and only if all the elements ofthe first column of the Routh array are positive

Method for determining the Routh array

Consider the characteristic equation

a( s) 1 s n a1s n 1
a2 s n 2
a3 s n 3
 a n 1 s 1 a0 s 0

Routh array method

Then add subsequent rows to complete the Routh array

Compute elements for the 3rd row:
1 a3 a2 a1
b1 ,
a1
1 a5 a4 a1
b2 ,
a1
1 a7 a6 a1
b3
a1

Given the characteristic equation,

a( s) s6 4s 5 3s 4 2s 3 s2 4s 4
Is the system described by this characteristic equation stable?

Answer:

All the coefficients are positive and nonzero

Therefore, the system satisfies the necessary condition for stability
We should determine whether any of the coefficients of the first column of the Routh array are negative

s6 : 1 3 1 4
5
s : 4 2 4 0
4
s : 52 0 4
3
s : 2 12 5 0
2
s : ? ?
1
s : ? ?
0
s : ?

s6 : 1 3 1 4
5
s : 4 2 4 0
4
s : 52 0 4
3
s : 2 12 5 0
2
s : 3 4
1
s : 76 15 0
0
s : 4

The elements of the 1st column are not all positive. Then the system is unstable.

Special cases of Routh’s criteria

Case 1: All the elements of a row in a RA are zero

Form Auxiliary equation by using the co-efficient of the row which is just above the row of zeros

Find derivative of the A.E.

Replace the row of zeros by the co-efficient of dA(s)/ds

 complete the array in terms of these coefficients

analyze for any sign change, if so, unstable

no sign change, find the nature of roots of AE

non-repeated imaginary roots - marginally stable

repeated imaginary roots – unstable

Case 2: First element of any of the rows of RA is

Zero and the same remaining row contains atleast one non-zero element

Substitute a small positive no. ‗ ‘ in place of zero and complete the array.

Examine the sign change by taking δt = 0

Root Locus Technique

Introduced by W. R. Evans in 1948

Graphical method, in which movement of poles in the s-plane is sketched when some parameter is varied
The path taken by the roots of the characteristic equation when open loop gain K is varied from 0 to ∞ are called
root loci
Direct Root δocus = 0 < k < ∞
Inverse Root Locus = - ∞ < k < 0

Root Locus Analysis

The roots of the closed-loop characteristic equation define the system characteristic responses
Their location in the complex s-plane lead to prediction of the characteristics of the time domain responses in terms
of:
o damping ratio,
o natural frequency, wn
o damping constant, first-order modes
o Consider how these roots change as the loop gain is varied from 0 to

Basics of Root Locus

Symmetrical about real axis

RL branch starts from OL poles and terminates at OL zeroes
No. of RL branches = No. of poles of OLTF
Centroid is common intersection point of all the asymptotes on the real axis
Asymptotes are straight lines which are parallel to Rδ going to ∞ and meet the Rδ at ∞
ζo. of asymptotes = ζo. of branches going to ∞
At Break Away point , the RL breaks from real axis to enter into the complex plane
At BI point, the RL enters the real axis from the complex plane

Constructing Root Locus

Locate the OL poles & zeros in the plot

Find the branches on the real axis

Find angle of asymptotes & centroid

 Φa= ±180º(2q+1) / (n-m)

 a = (Σpoles - Σzeroes) / (n-m)
 Find BA and BI points

Find Angle Of departure (AOD) and Angle Of Arrival (AOA)

AOD = 180º- (sum of angles of vectors to the complex pole from all other poles) +

(Sum of angles of vectors to the complex pole from all zero)

AOA = 180º- (sum of angles of vectors to the complex zero from all other zeros) + (sum of angles of vectors to the
complex zero from poles)

Find the point of intersection of RL with the imaginary axis.

Application of the Root Locus Procedure

Step 1: Write the characteristic equation as

1 F ( s) 0

Step 2: Rewrite preceding equation into the form of poles and zeros as follows
m
(s zj)
j 1
1 K n
0
(s pi )
i 1

Step 3:

Locate the poles and zeros with specific symbols, the root locus begins at the open-loop poles and ends at the open-
loop zeros as K increases from 0 to infinity
If open-loop system has n-m zeros at infinity, there will be n-m branches of the root locus approaching the n-m
zeros at infinity

Step 4:

The root locus on the real axis lies in a section of the real axis to the left of an odd number of real poles and zeros

Step 5:

The number of separate loci is equal to the number of open-loop poles

Step 6:

The root loci must be continuous and symmetrical with respect to the horizontal real axis

Step 7:

The loci proceed to zeros at infinity along asymptotes centered at centroid and with angles

n m
pi zj
i 1 j 1
a
n m
(2k 1)
a (k 0,1, 2, n m 1)
n m
Step 8:

The actual point at which the root locus crosses the imaginary axis is readily evaluated by using Routh‘s criterion
Step 9:

Determine the breakaway point d (usually on the real axis)

Step 10:

Plot the root locus that satisfy the phase criterion

P(s) (2k 1) k 1, 2,

Step 11:

Determine the parameter value K1 at a specific root using the magnitude criterion
n
(s pi )
i 1
K1 m

Nyquist Stability Criteria (s zj)

j 1 s s1

The Routh-Hurwitz criterion is a method for determining whether a linear system is stable or not by examining the
locations of the roots of the characteristic equation of the system. In fact, the method determines only if there are roots that
lie outside of the left half plane; it does not actually compute the roots. Consider the characteristic equation.
To determine whether this system is stable or not, check the following conditions

1. Two necessary but not sufficient conditions that all the roots have negative real parts are
a) All the polynomial coefficients must have the same sign.
b) All the polynomial coefficients must be nonzero.
2. If condition (1) is satisfied, then compute the Routh-Hurwitz array as follows

where the ai's are the polynomial coefficients, and the coefficients in the rest of the table are computed using the
following pattern
3. The necessary condition that all roots have negative real parts is that all the elements of the first column of the array have
the same sign. The number of changes of sign equals the number of roots with positive real parts.
4. Special Case 1: The first element of a row is zero, but some other elements in that row are nonzero. In this case, simply
replace the zero elements by " ", complete the table development, and then interpret the results assuming that " " is a small
number of the same sign as the element above it. The results must be interpreted in the limit as 0.
5. Special Case 2: All the elements of a particular row are zero. In this case, some of the roots of the polynomial are located
symmetrically about the origin of the s-plane, e.g., a pair of purely imaginary roots. The zero rows will always occur in a
row associated with an odd power of s. The row just above the zero rows holds the coefficients of the auxiliary polynomial.
The roots of the auxiliary polynomial are the symmetrically placed
roots. Be careful to remember that the coefficients in the array skip powers of s from one coefficient to the next.
Let P = no. of poles of q(s)-plane lying on Right Half of s-plane and encircled by s-plane contour.
Let Z = no. of zeros of q(s)-plane lying on Right Half of s-plane and encircled by s-plane contour.

For the CL system to be stable, the no. of zeros of q(s) which are the CL poles that lie in the right half of s-plane
should be zero. That is Z = 0, which gives N = -P.
Therefore, for a stable system the no. of ACW encirclements of the origin in the q(s)-plane by the contour Cq must
be equal to P.

Nyquist modified stability criteria

We know that q(s) = 1+G(s)H(s)

Therefore G(s)H(s) = [1+G(s)H(s)] – 1
The contour Cq, which has obtained due to mapping of Nyquist contour from s-plane to q(s)-plane (ie)[1+G(s)H(s)]
-plane, will encircle about the origin.
The contour CGH, which has obtained due to mapping of Nyquist contour from s-plane to G(s)H(s) -plane, will
encircle about the point (-1+j0).
 Therefore encircling the origin in the q(s)-plane is equivalent to encircling the point -1+j0 in the G(s)H(s)-plane.

St
ab
Un
sta
Problem

Sketch the Nyquist stability plot for a feedback system with the following open-loop transfer function

1
G( s) H ( s) 2
s(s s 1)

Solution

For section ab, s = j , :0

1
G( j ) H ( j ) 2
j (1 j )

(i) 0 : G(j )H(j ) –1 – j

(ii) = 1 : G(j )H(j ) –1+ j0

(iii) : G(j )H(j ) 0 –270o

On section bcd, s = Rej R ; therefore i.e. section bcd maps onto the origin of the G(s)H(s)-plane

1
G( s) H ( s) 0
R3
Section de maps as the complex image of the polar plot as before
Relative stability

The main disadvantage of a Bode plot is that we have to draw and consider two different curves at a time, namely,
magnitude plot and phase plot. Information contained in these two plots can be combined into one named polar plot.
The polar plot is for a frequency range of 0 . while the Nyquist plot is in the frequency range of The
information on the negative frequency is redundant because the magnitude and real part of G ( jw) an are even functions.
. In this section. We consider how to evaluate the system performance in terms of relative stability using a Nyquist plot.
The open-loop system represented by this plot will become unstable beyond a certain value. As shown in the Nyquist plot of
Fig. the intercept of magnitude 'a on the negative real axis corresponds lost phase shift of - 180° and - 1 represents the
amount of increase in gain that can be tolerated before closed-loop system tends toward instability. As 'a' approaches
( 1 j 0) point the relative stability is reduced, The gain and phase margins are represented as follows in the Nyquist plot.
Gain margin
1
As system gain is increased by a factor 1/a, the open loop G( jw) H ( jw) will increase by a factor a ( ) 1and the
a
system would be driven to instability.Thus, the gain margin is the reciprocal of the gain at the frequency at which the phase
angle of the Nyquist plot is - 1800. The gain rnargin, usually measured in dB, is a positive quantity given by
GM = 20 log adB

Phase Margin m

Importance of the phase margin has already in the content of Bode. Phase margin is defined as the change in open-loop phase
shift required al unity gain to make a closed loop system unstable. A closed-loop system will be unstable if the Nyquist plot
encircles 1 j 0 point. Therefore, the angle required to make this system marginally stable in a closed loop is the phase
margin .In order to measure this angle, we draw a circle with a radius of 1, and find the point of intersection of the Nyquist
plot with this circle, and measure the phase shift needed for this point to be at an angle of 1800. If may be appreciated that the
system having plot of Fig with larger PM is more stable than the one with plot of Fig.
 UNIT IV

STABILITY ANALYSIS
PART-A

1. Define BIBO Stability.

2. What is impulse response?
3. What is characteristic equation?
4. How the roots of characteristic equation are related to stability?
5. What is the necessary condition for stability?
6. What is the relation between stability and coefficient of characteristic polynomial?
7. What will be the nature of impulse response when the roots of characteristic equation
8. What will be the nature of impulse response if the roots of characteristic equation are
9. What is principle of argument?
10. What is the necessary and sufficient condition for stability?
11. What is routh stability condition?
12. What is auxiliary polynomial?
13. What is quadrantal symmetry?
14. In routh array what conclusion can you make when there is a row of all zeros?
15. What is limitedly stable system?
16. What is Nyquist stability criterion?
17. What is root locus?
18. How will you find root locus on real axis?
19. What are asymptotes?
20. What is centroid, how it is calculated?
21. What is breakaway point and break in point?
22. What is dominant pole?
PART-B
1. Using Routh criterion determine the stability of the system whose characteristics
equation is S4+8S3+18S2+16S+5 =0 . (16)
2. F(S) =S6 +S5-2S4-3S3-7S2-4S1-4 =0.Find the number of roots falling in the RHS plane and LHS plane. (16)
3. Draw the Nyquist plot for the system whose open loop transfer function is
G(S) H(S) =K/S (S+2) (S+10).Determine the range of K for which closed loop system is stable. (16)
4. Construct Nyquist plot for a feedback control system whose open loop transfer
function is given by G(S)H(S) =5/ S(1-S).comment on the stability of open loop and
closed loop transfer function. (16)
5. Sketch the Nyquist plot for a system with the open loop transfer function
G(S)H(S) =K (1+0.5S) (1+S) / (1+10S) (S-1).determine the range of values of K for
which the system is stable. (16)
6. Sketch the root locus for the open loop transfer function of unity feedback control system given below: G(S)
H(S)=K/S(S+2)(S+4). (16)