NPR COLLEGE OF ENGINEERING & TECHNOLOGY

DEPARTMENT OF ELECTRIAL & ELECTRONICS ENGINEERING

ACADEMIC YEAR 2012-13

SUB NAME: CONTROL SYSTEMS

SUB CODE: EE2253
YEAR/SEM: II/IV
LECTURE NOTES

Prepared By:
A.R.SALINIDEVI
AP/EEE
DEPARTMENT OF EEE

EE2253

CONTROL SYSTEMS

L T PC
3104

(Common to EEE, EIE & ICE)

UNIT I SYSTEMS AND THEIR REPRESENTATION
Basic elements in control systems Open and closed loop systems Electrical analogy of
Mechanical and thermal systems Transfer function Synchros AC and DC
servomotors Block diagram reduction techniques Signal flow graphs.
UNIT II TIME RESPONSE
Time response Time domain specifications Types of test input I and II order
system response Error coefficients Generalized error series Steady state error P,
PI, PID modes of feedback control.
UNIT III FREQUENCY RESPONSE
Frequency response Bode plot Polar plot Determination of closed loop response
from open loop response Correlation between frequency domain and time domain
specifications.
UNIT IV STABILITY OF CONTROL SYSTEM
Characteristics equation Location of roots in S plane for stability Routh Hurwitz
criterion Root locus construction Effect of pole, zero addition Gain margin and
phase margin Nyquist stability criterion.
UNIT V COMPENSATOR DESIGN
Performance criteria Lag, lead and lag-lead networks Compensator design using
bodePlots.
L: 45 T: 15 TOTAL: 60
PERIODS
TEXT BOOKS
1. Nagrath.I.J and Gopal.M, Control Systems Engineering, New Age International
Publishers, 2003.
2. Benjamin C. Kuo, Automatic Control systems, Pearson Education, New Delhi,
2003.
REFERENCE BOOKS
1. Ogata.K, Modern Control Engineering, 4th edition, PHI, New Delhi, 2002.
2. Norman S. Nise, Control Systems Engineering, 4th Edition, John Wiley, New Delhi,
2007.
3. Samarajit Ghosh, Control systems, Pearson Education, New Delhi, 2004
4. Gopal.M, Control Systems, Principles and Design, Tata McGraw Hill, New Delhi,
2002.

UNIT I
SYSTEMS AND THEIR REPRESENTATION

Basic elements in control systems Open and closed loop systems Electrical
analogy of Mechanical and thermal systems Transfer function Synchros AC
and DC servomotors Block diagram reduction techniques Signal flow graphs.
Control System

Control System means any quantity of interest in a machine or mechanism is

maintained or altered in accordance with desired manner

A system which controls the output quantity is called a control system

Definitions:

Controlled Variable: It is the quantity or condition that is measured & controlled

Controller: Controller means measuring the value of the controlled variable of the system &
applying the manipulated variable to the system to correct or to limit the deviation of
the measured value to the desired value.

Plant: A plant is a piece of equipment, which is a set of machine parts functioning

together. The purpose of which is to perform a particular operation. Example:
Furnace, Space craft etc.,

System: A system is a combination of components that works together & performs

certain objective.

Disturbance: A disturbance is a signal that tends to affect the value of the output of a
system. If disturbances created inside the system, it is called internal While an
external disturbance is generated outside the system.

Feedback Control: It is an operation that, in the presence of disturbance tends to

reduce the difference between the output of a system & some reference input.

Servo Mechanism: servo mechanism is a feedback controlled system in which the

output is some mechanical position, velocity or acceleration.

Open loop System: In an Open loop System, the control action is independent of the
desired output. OR When the output quantity of the control system is not fed back to
the input quantity, the control system is called an Open loop System.

Closed loop System: In the Closed loop Control System the control action is
dependent on the desired output, where the output quantity is considerably controlled
by sending a command signal to input quantity
3

Classification of Control Systems

The Control System can be classified mainly depending upon,
(a) Method of analysis & design, as Linear & Non- Linear Systems.
(b)The type of the signal, as Time Varying, Time Invariant, Continuous data, Discrete data
systems etc.
(c)The type of system components, as Electro Mechanical, Hydraulic, Thermal, Pneumatic
Control systems etc.
(d)The main purpose, as Position control & Velocity control Systems.
Linear & Non -Linear Systems:
In a linear system, the principle of superposition can be applied. In non- linear system,
principle & homogeneity.
Time Varying & Time Invariant Systems:

While operating a control system, if the parameters are unaffected by the time, then
the system is called Time Invariant Control System.

Most physical systems have parameters changing with time. If this variation is
measurable during the system operation then the system is called Time Varying
System.

If there is no non-linearity in the time varying system, then the system may be called
as Linear Time varying System.

Discrete Data Systems:

If the signal is not continuously varying with time but it is in the form of pulses.
Thenthe control system is called Discrete Data Control System.

If the signal is in the form of pulse data, then the system is called Sampled Data
Control System.

Here the information supplied intermittently at specific instants of time. This has the
advantage of Time sharing system.

On the other hand, if the signal is in the form of digital code, the system is called
Digital Coded System.

Here use of Digital computers, p, c is made use of such systems are analyzed by
the Z- transform theory.

Continuous Data Systems:

If the signal obtained at various parts of the system are varying continuously with time,then
the system is called Continuous Data Control Systems.
4

Adaptive Control systems:

In some control systems, certain parameters are either not constant or vary in an
unknown manner.

If the parameter variations are large or rapid, it may be desirable to design for the
capability of continuously measuring them & changing the compensation, so that the
system performance criteria can always satisfied.

This is called Adaptive Control Systems.

Optimal Control System:

Optimal Control System is obtained by minimizing and/or maximizing the performance
index. This index depends upon the physical system & skill.
Single Variable Control System:
In simple control system there will be One input & One output such systems are calledSingle
variable System (SISO Single Input & Single Output).
Multi Variable Control System:
In Multivariable control system there will be more than one input & correspondingly more
output - Multiple Inputs & Multiple Outputs)
Basic elements of control 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.

CONTROL
SYSTEM

Input/
Objective

Output/
Result

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. 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.

Desired
speed

Controller

Engine

Actual
speed

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

+
Desired
Speed

Engine

Controller

Actual
Speed

_
Speed
Sensor

Comparison between Open loop & Closed loop Gain

Open loop System

Closed loop System

An open loop system has the ability

to perform accurately, if its
calibration is good. If the calibration
is not perfect its performance will go
down

A closed loop system has got the

ability to perform accurately because
of the feed back

It is easier to build

It is difficult to build

In general it is more stable as the

feedback is absent

Less Stable Comparatively

If non- present; the system operation

is not good

Even under the presence of non loop

system

Feed back is absent.

Feed back is present.

Example :(i)Traffic Control System.

Example :( i) Pressure Control

System.

(ii)Control of furnace for coal

heating.

(ii) Speed Control System

(iii)An Electric Washing Machine.

(iii) Robot Control System.

(iv) Temperature Control System

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

The input- output relationship in a linear time invariant system is defined by the
transfer function.

The features of the transfer functions are,

7

 It is applicable to Linear Time Invariant system.

It is the ratio between the Laplace Transform of the o/p variable to the
Laplace Transform of the i/p variable.
It is assumed that initial conditions are zero.
It is independent of i/p excitation.

It is used to obtain systems o/p response.

An equation describing the physical system has integrals & differentials, the step involved in
obtaining the transfer function are;
(1) Write the differential equation of the system.
(2) Replace the terms

d
by S &
dt

dt by 1/S.

(3) Eliminate all the variables except the desired variables.

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.
8

Mathematical Model of control systems

A physical system is a collection of physical objects connected together to serve an

objective. An idealized physical system is called a Physical model.

Once a physical model is obtained, the next step is to obtain Mathematical model.
When a mathematical model is solved for various i/p conditions, the result represents
the dynamic behavior of the system.

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.

Analogous System:
The concept of analogous system is very useful in practice. Since one type of system may be
easier to handle experimentally than another. A given electrical system consisting of
resistance, inductance & capacitances may be analogous to the mechanical system consisting
of suitable combination of Dash pot, Mass & Spring.
The advantages of electrical systems are,

Many circuit theorems, impedance concepts can be applicable.

An Electrical engineer familiar with electrical systems can easily analyze the system
under study & can predict the behavior of the system.

The electrical analog system is easy to handle experimentally.

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 Newtons second law of motion. For
translational systems it states that the sum of forces acting on a body is zero.
Translational System:
It has 3 types of forces due to elements.
1. Inertial Force: Due to inertial mass,

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
Here fm M

d 2x
dt 2

By Newtons second law, f = f m= M

d 2x
dt 2

2. Damping Force [Viscous Damping]: Due to viscous damping, it is proportional to

velocity & is given by,

Damping force is denoted by either D or B or F

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
10

Let f = applied force

f b = opposing force due to friction
Here, f b B

dx
dt

By Newtons second law, f = f b B

dx
dt

3. Spring Force: Spring force is proportional to displacement.

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 Newtons second law, f = f k x
Mechanical Rotational Systems
The model of rotational mechanical systems can be obtained by using three elements,
11

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 Newtons 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
Here Tj J

d 2
dt 2

By Newtons law
d 2
T= Tj = J 2
dt

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
Here Tb B

d
(1 2 )
dt

By Newtons law

12

T= Tb = B

d
(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 Newtons law
T = Tk K (1 2 )
Analogous quantities in translational & Rotational system:
The electrical analog of the mechanical system can be obtained by,
(i) Force Voltage analogy: (F.V)
(ii) Force Current analogy: (F.I)

Mechanical
Translational System

Mechanical
Rotational System

F.VAnalogy

Force (F)

Torque (T)

Voltage (V)

Current (I)

Mass (M)

Moment of Inertia
(M)

Inductance
(L)

Capacitance (C)

Viscous friction (D or B
or F)

Viscous friction (D
or B or F)

Resistance
(R)

Conductance (G)

Spring stiffness (k)

Torsional spring
stiffness

Reciprocal of
Capacitance
(1/C)

Reciprocal of
Inductance (1/L)

Linear displacement ( x ) Angular

Charge (q)

Flux ( )

Linear velocity ( u )

Current (i)

Voltage (v)

Sl.No.

F.IAnalogy

displacement ( )

Angular Velocity
(w)

13

Modeling of electrical system

Electrical circuits involving resistors, capacitors and inductors are considered. The
behaviour of such systems is governed by Ohms law and Kirchhoffs 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 L 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 Kirchhoffs voltage law or Kirchhoffs 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

14

R1i(t ) R2 i (t )

1
i (t )dt v1 (t )
c 0

1
R2 i (t ) i (t )dt v2 (t )
c0
Electrical systems
LRC circuit. Applying Kirchhoffs 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.
15

 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 ia

where is the air gap flux, = kf if,

For the constant flux

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

16

k1 is constant

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)

17

The T.F can be obtained is

Field-Controlled dc motors

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
18

T = k1 ia

where is the air gap flux, = kf if ,

k1 is constant

In the field controlled DC motor, the armature current is fed from a constant current source.

The KVL equation for the field circuit is,

On Laplace Transform,

The torque equation is,

Where, KT is Motor torque constant.

On Laplace Transform,

19

The complete block diagram is as shown below,

Substituting the value of I f (s ) from equation (2) in equation (1), we get

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 2 D Kx f ----------------- (1)
dt
dt
20

From the fig (b)

We get J

d 2
d
D
K T -------------- (2)
2
dt
dt

From the fig (c)

We get

d 2q
dq 1
R q v(t ) ------------ (3)
2
dt C
dt

Where q i dt
They are two methods to get analogous system.
These are
(i)

force- voltage (f-v) analogy and

(ii)

(ii) force-current (f-c) analogy

Force Voltage Analogy

Force Current Analog

21

PROBLEM

1. 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

22

From eq (2) we get

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

Forcecurrent analogy

From eq (1) we get

From eq (2) we get

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

23

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.
It is a representation of the control system giving the inter-relation between the
transfer function of various components.
The block diagram is obtained after obtaining the differential equation & Transfer
function of all components of a control system.
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.
The arrow head pointing towards the block indicates the i/p
& pointing away from the block indicates the o/p.

Transfer function G(s)

After obtaining the block diagram for each & every component, all blocks are combined to
obtain a complete representation. It is then reduced to a simple form with the help of block
diagram algebra.
24

Basic elements of a block diagram

Blocks

Transfer functions of elements inside the blocks

Summing points

Take off points

Arrow

Block diagram

Branch point

Summing 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.
25

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

26

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) Feedback 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-andeffect 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

27

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

28

Eliminating a feedback loop

Cascaded Subsystems

Parallel Subsystems

Feedback Control System

29

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. Reduce the Block Diagrams shown below:

Solution:
By eliminating the feed-back paths, we get

30

Combining the blocks in series, we get

Eliminating the feedback path, we get

Transfer function

2. Obtain the Transfer function of the given block diagram

31

Solution:
Shifting the take-off T1 beyond the block G3 we get

Combining G2&G3 and eliminating H3 we get

Eliminating the feedback path

32

Transfer function

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

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

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
34

Signal Flow Graph Representation

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

For complicated systems, Block diagram reduction method becomes tedious &
time consuming. An alternate method is that signal flow graphs developed by S.J.
Mason.

In these graphs, each node represents a system variable & each branch connected
between two nodes acts as Signal Multiplier.

The direction of signal flow is indicated by an arrow.

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.

Definitions

Node: A node is a point representing a variable.

Transmittance: A transmittance is a gain between two nodes.
Branch: A branch is a line joining two nodes. The signal travels along a branch.
35

 Input node [Source]: It is a node which has only out going signals.
Output node [Sink]: It is a node which is having only incoming signals.
Mixed node: It is a node which has both incoming & outgoing branches (signals).
Path: It is the traversal of connected branches in the direction of branch arrows. Such
that no node is traversed more than once.
Loop: It is a closed path.
Loop Gain: It is the product of the branch transmittances of a loop.
Non-Touching Loops: Loops are Non-Touching, if they do not possess any common
node.
Forward Path: It is a path from i/p node to the o/p node w hich doesnt cross any
node more than once.
Forward Path Gain: It is the product of branch transmittances of a forward path.
Masons 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. Masons gain rule for the determination of the overall system gain is given below.

Where M= gain between Xin and Xout

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 non touching loops)
(sum of gain products of all possible combination of three nontouching loops)

PROBLEM

1. Obtain the closed loop TF C(S)/R(S) by using Masons gain formula.

36

Solution:
Masons gain formula is,

No of forward paths

No. of individual loops:

Gain Products of all possible combinations of two non-touching loops:

No. of three non-touching loops = 0.

37

(2) Obtain the closed loop TF, by using Masons gain formula.

Solution:
Masons gain formula is,

Forward Paths:

No. of individual loops:

Two Non-touching loops:

38

3. Obtain C(s)/R(s) by Block Diagram Reduction method & verify the result by signal flow
graph

Solution:
Re-arranging the summing points,

39

Signal flow graphs:

No. of forward paths:

No. of individual loops:

Masons gain formula is,

Modeling of Elements of Control Systems

A feedback control system usually consists of several components in addition to the actual
process. These are error detectors, power amplifiers, actuators, sensors etc. Let us now
discuss the physical characteristics of some of these and obtain their mathematical models.
DC Servo Motor

A DC servo motor is used as an actuator to drive a load. It is usually a DC motor of

low power rating.

40

DC servo motors have a high ratio of starting torque to inertia and therefore they have
a faster dynamic response.

DC motors are constructed using rare earth permanent magnets which have high
residual flux density and high coercively.

As no field winding is used, the field copper losses am zero and hence, the overall
efficiency of the motor is high.

The speed torque characteristic of this motor is flat over a wide range, as the armature
reaction is negligible.

Moreover speed in directly proportional to the armature voltage for a given torque.
Armature of a DC servo motor is specially designed to have low inertia.

In some application DC servo motors are used with magnetic flux produced by field
windings.

The speed of PMDC motors can be controlled by applying variable armature voltage.
These are called armature voltage controlled DC servo motors.

Wound field DC motors can be controlled by either controlling the armature voltage
or controlling rho field current. Let us now consider modelling of these two types or
DC servo motors.

(a) Armature controlled DC servo motor

The physical model of an armature controlled DC servo motor is given in

The armature winding has a resistance R a and inductance La.

The field is produced either by a permanent magnet or the field winding is separately
excited and supplied with constant voltage so that the field current If is a constant.

When the armature is supplied with a DC voltage of e a volts, the armature rotates and
produces a back e.m.f eb.
41

The armature current ia depends on the difference of eb and en. The armature has a
remnant of inertia J, frictional coefficient B0

The angular displacement of the motor is 8.

The torque produced by the motor is given by

Where K T is the motor torque constant.

The back emf is proportional to the speed of the motor and hence

The differential equation representing the electrical system is given by

Taking Laplace transform of equation from above equation

The mathematical model of the mechanical system is given by

Taking Laplace transform

Solving for (s) ,we get

42

The block diagram representation of the armature controlled DC servo motor is developed in
steps

Combining these blocks we have

Usually the inductance of the armature winding is small and hence neglected

Where

Field Controlled Dc Servo Motor

The field servo motor

43

The electrical circuit is modeled as

Where
Motor gain constant

Motor time constant

Field time constant

The block diagram is as shown as

44

AC Servo Motors

An AC servo motor is essentially a two phase induction motor with modified

constructional features to suit servo applications.

The schematic of a two phase or servo motor is shown

It has two windings displaced by 90 on the stator One winding, called as reference
winding, is supplied with a constant sinusoidal voltage.

The second winding, called control winding, is supplied with a variable control
voltage which is displaced by -- 90 out of phase from the reference voltage.

The major differences between the normal induction motor and an AC servo motor are

The rotor winding of an ac servo motor has high resistance (R) compared to its
X
inductive reactance (X) so that its R ratio is very low.
X
For a normal induction motor, R ratio is high so that the maximum torque is obtained
in normal operating region which is around 5% of slip.

The torque speed characteristics of a normal induction motor and an ac servo motor
are shown in fig

The Torque speed characteristic of a normal induction motor is highly nonlinear and
has a positive slope for some portion of the curve.

This is not desirable for control applications. as the positive slope makes the systems
unstable. The torque speed characteristic of an ac servo motor is fairly linear and has
45

negative slope throughout.

The rotor construction is usually squirrel cage or drag cup type for an ac servo motor.
The diameter is small compared to the length of the rotor which reduces inertia of the
moving parts.

Thus it has good accelerating characteristic and good dynamic response.

The supplies to the two windings of ac servo motor are not balanced as in the case of
a normal induction motor.

The control voltage varies both in magnitude and phase with respect to the constant
reference vulture applied to the reference winding.

The direction of rotation of the motor depends on the phase ( 90) of the control
voltage with respect to the reference voltage.

For different rms values of control voltage the torque speed characteristics are shown
in Fig.

The torque varies approximately linearly with respect to speed and also controls
voltage.

The torque speed characteristics can be linearised at the operating point and the
transfer function of the motor can be obtained.

Synchros
A commonly used error detector of mechanical positions of rotating shafts in AC control
systems is the Synchro.
It consists of two electro mechanical devices.

Synchro transmitter

Synchro receiver or control transformer.

46

The principle of operation of these two devices is sarne but they differ slightly in their
construction.

The construction of a Synchro transmitter is similar to a phase alternator.

The stator consists of a balanced three phase winding and is star connected.

The rotor is of dumbbell type construction and is wound with a coil to produce a
magnetic field.

When a no voltage is applied to the winding of the rotor, a magnetic field is produced.

The coils in the stator link with this sinusoidal distributed magnetic flux and voltages
are induced in the three coils due to transformer action.

Than the three voltages are in time phase with each other and the rotor voltage.

The magnitudes of the voltages are proportional to the cosine of the angle between the
rotor position and the respective coil axis.

The position of the rotor and the coils are shown in Fig.

When 90 the axis of the magnetic field coincides with the axis of coil S2 and
maximum voltage is induced in it as seen.
47

For this position of the rotor, the voltage c, is zero, this position of the rotor is known
as the 'Electrical Zero' of die transmitter and is taken as reference for specifying the
rotor position.

In summary, it can be seen that the input to the transmitter is the angular position of
the rotor and the set of three single phase voltages is the output.

The magnitudes of these voltages depend on the angular position of the rotor as given

Hence

Now consider these three voltages to he applied to the stator of a similar device called
control transformer or synchro receiver.

The construction of a control transformer is similar to that of the transmitter except

that the rotor is made cylindrical in shape whereas the rotor of transmitter is dumbell
in shape.

Since the rotor is cylindrical, the air gap is uniform and the reluctance of the magnetic
path is constant.

This makes the output impedance of rotor to be a constant.

Usually the rotor winding of control transformer is connected teas amplifier which
requires signal with constant impedance for better performance.

A synchro transmitter is usually required to supply several control transformers and

hence the stator winding of control transformer is wound with higher impedance per
phase.

Since the some currents flow through the stators of the synchro transmitter and
receiver, the same pattern of flux distribution will be produced in the air gap of the
control transformer.
48

The control transformer flux axis is in the same position as that of the synchro
transmitter.

Thus the voltage induced in the rotor coil of control transformer is proportional to the
cosine of the angle between the two rotors.

**************************

UNIT I
CONTROL SYSTEM MODELLING
49

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)

50

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)
51

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

Time response Time domain specifications Types of test input I and II order
system response Error coefficients Generalized error series Steady state error
P, PI, PID modes of feedback control

Time response analysis

52

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 systems performance.

The first step in the analysis of a control system is, describing the system in terms of a
mathematical model.

The system is modelled by defining its transfer function.

The next step would be, to obtain its response, both transient and steady state, to a
specific input.

The input can be a time varying function which may be described by known
mathematical functions or it may be a random signal.

Moreover these input signals may not be known apriority.

Thus it is customary to subject the control system to some standard input test signals
which strain the system very severely.

These standard input signals arc: an impulse. a step, a ramp and a parabolic input.
Analysis and design of control systems are carried out, defining cream performance
measures for the system, using these standard test signals.

It is also pertinent to mention that any arbitrary time function can be expressed in
terms of linear combinations of these test signals and hence, if the system is linear, the
output of the system can be obtained easily by using superstition principle.

Further, convolution integral can also be used to determine the response of a linear
system for any given input, if the response is known for a step or an impulse input.

Normally use the standard input signals to identify the characteristics of systems
response
Step function
Ramp function
Impulse function
53

 Parabolic function
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 state response.
Time response = Transient response + Steady 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).

54

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

Standard Test Signals

Impulse Signal
An impulse signal is shown in Fig.

The impulse function is zero for all t 0 and its infinity at t 0 .It rises to infinity at
t 0 and comes uback to zero at t 0 enclosing a finite area. If this area is A it is called as
impulse function of strength A. If A 1 it is called a unit impulse function. Thus an impluse
signal is denoted by f (t ) A (t ) .

step Signal
A Step signal is shown in Fig.
55

It is zero for t 0 and suddenly rises to a value A at t=0 and remains at this value for t>0: it is
denoted by f (t ) A u(t ) . If A=1, it is called a unit step function.

Ramp Signal
A Ramp signal is shown in Fig.

It is zero for t<0 and uniformly increases with a slope equal to A. It is denoted by
f (t ) A t .If the slope is unity then it is called a unit ramp signal.

Parabolic Signal
A Parabolic Signal is shown in Fig.

A parabolic signal is denoted by f (t )

At 2
If A is equal to unity it is known as a unit
2

parabolic signal.

Impulse response

Exponential

K
1 sT

56

Step response

Ramp response

K
K
s s 1/ T
K KT
KT
2
s
s s 1/ T

Step, exponential

Ramp, step, exponential

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

First order system

Response to a unit step input
Consider a feedback system with G( s)

1
as shown in fig
s

The closed loop transfer function of the system is given by

57

---------------------------(1)
For a unit step input R ( s )

1
and the output is given by
s

---------------------------(2)
Inverse laplace transformation yields
---------------------- (3)
The plot of c(t) Vs t is shown in the fig

The response is an exponentially increasing function and it approaches a value of

unity t .At t the response reaches a value

This is 63.2 percent of the steady value. This time is known as the here constant of
the system. One of the characteristics which we would like to know about the system
is its speed of response or how fast the response is approaching the final value. The
time constant is indicative of this measure and the speed of response is inversely
proportional to the time constant of the system.
Another important characteristic of the system is the error between the desired value
and the actual value under steady state conditions. This quantity is known as the
steady state error of the system and is denoted by e ss
The error E(s) for a unity feedback system is given by
E(s) =R(s)-C(s)

For the system under consideration

58

At t e(t ) 0 . Thus the output of the first order system approaches the reference
input. Which is the desired output,without any error.
Response to a Unit Ramp Input or Unit Velocity Input
The response of the system in Fig.for a unit ramp input, for which.

is given by,

The time response is obtained by taking inverse Laplace transform of eqn.

If above eqn. is differentiated we get

which is the response of the system to a step input. Thus no additional infommtion
about
the speed of response is obtained by considering a ramp input. But let us see
the effect on the steadystate error As before,

Thus the steady state error is equal to the time constant of the system. The first order system,
therefore, cannot track the ramp input without a finite steady state error. If the time constant
is reduced not only the speed of response increases but also the steady state error for ramp
59

input decreases. Hence the ramp input is important to the extent that it produces a finite
steady state error. Instead of finding the entire response, it is sufficient to estimate the steady
state value by using the final value theorem. Thus

The response of a first order system for unit ramp input is plotted in Fig.

Response to a Unit Parabolic or Acceleration Input

The response of a first order system to a unit parabolic input, for which

Differentiating above equation

which is the response of the first order system to unit velocity input. Thus subjecting the first
order system to a unit parabolic input does not give any additional information regarding
transient behavior of the system. But, the steady state error, for a parabolic input is given by.
60

Thus a first order system has infinite state error for a parabolic input .The steady state error
can be easily obtained by using the final value theorem as:

Summarizing the analysis of first order system. we can say that the step input yields the
desired information about the speed of transient response. Isis observed that the speed of
response is inversely proportional lathe time constant r of the system. The ramp and parabolic
inputs do not give any additional information regarding the speed of response. However, the
steady statesman is different for these three different inputs. For a step input, the steady state
error e0 is zero, for a velocity input there is a finite error equal to the time constant r of the
system and for an acceleration input the steady state error is infinity.
It is clear from the discussion above, that it is sufficient to study the behavior of any system
to a unit step input for understanding its transient response and use the velocity input and
acceleration input fur understanding the steady state behavior of the system.
Response to a unit ramp input
For a unit ramp input

And the output is given by

Taking inverse Laplace transform, we get the time response c(t)

The time response for a unit ramp input is plotted

61

Second order system

Response to Unit Step Input
Consider a Type second order system as shown in Fig. Since G(s) has Pole at the origin,

The closed loop transfer function is given by

The transient response of any system depends on the poles of the transfer function T(s).the
roots of the denominator polynomial in s of T(s) are the poles of the transfer function. Thus
the denominator polynomial of T(s) is given by

Is known as the characteristic polynomial of the system and D(S)=0is known as the
characteristic equation of the system.

--------------(1)
Where

62

= natural frequency

2 K

= damping factor

The poles of T(s), or the roots of the characteristic equation

s 2 2n s n 0
2

2n s 4 2 n 4 n
2

s1, 2

n j n 1 2

Are given by

n j d
Where d n 1 2 is known as the damped natural frequency of the system. If >1, the

roots s1, s2 are real and we have an over damped system. If =1 the system is known as a
critically damped system. The more common case of <1 is known as the under damped
system.
For a unit step input R(s) =1/s and equation (1) can be written as
C(s) =T(s).R(s)

n 2
1
. --------------- (2)
2
2
s 2n s n s

Splitting eqn.(2) into partial fractions assuming to be less than 1, we have

K 2 s K3
K1
2
s s 2n s n 2

Evaluating K1, K2 and K3 by the usual procedure we have

C ( s)

s 2n
1

s ( s n ) 2 n 2 (1 2 )

n (1 2 )
s n
1

s ( s n ) 2 n 2 (1 2 )
( s n ) 2

63

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

c(t ) k1e 1t k2e 2t

Poles: Two complex at

Underdamped response:

d j d
Natural response: Damped sinusoid with an exponential envelope whose time constant is equal to the
reciprocal of the poles real part. The radian frequency of the sinusoid, the damped frequency of oscillation,
is equal to the imaginary part of the poles

c(t ) Ae d t cos(d t )
Undamped response:
Poles: Two imaginary at

j1
64

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
1

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

c(t ) k1e 1t k2te1t

Second order system responses damping cases

Second- order step response

complex poles

Time Domain Specifications of a Second Order System

The performance of a system is usually evaluated in terms of the following qualities. .

How fast it is able to respond to the input.

How fast it is reaching the desired output
65

What is the error between the desired output and the actual output, once the transients
die down and steady slate is achieved
Does it oscillate around the desired value, and
Is the output continuously increasing with time or is it bounded.
The last aspect is concerned with the stability of the system and we would require the system
to be stable. This aspect will be considered later. The first four questions will be answered in
terms of time domain specifications of the system based on its response to a unit step input.
These are the specifications to be given for the design of a controller for a given system.
We have obtained the response of a type 1 second order system to a unit step input. The step
response of a typical underdamped second order system is plotted in Fig.

It is observed that, for an underdamped system, there are two complex conjugate poles.
Usually, even if a system is of higher order, the two complex conjugate poles nearest to the
j - axis (called dominant poles) are considered and the system is approximated by a second
order system. Thus, in designing any system, certain design specifications are given based on
the typical underdamped step response shown as Fig.
The design specifications are
Delay time t d : It is the time required for the response to reach 50% of the steady state value
for the first time.
Rise time t r : It is the time required for the response to mach 100% of the steady state value
for under damped systems. However, for over damped systems, it is taken as the time
required for the response to rise from 10% to 90% of the steady state value.
Peak time t p : It is the time required for the response to reach the maximum or Peak value of
the response.
Peak overshoot M p : It is defined as the difference between the peak value of the response
and the steady state value. Iris usually expressed in percent of the steady state value. If the
time for the peak is t p , percent peak overshoot is given by,

66

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

ess lime(t) limsE(s) lim
t

s 0

s 0

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
67

sR(s)
1 G(s)

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
68

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.
69

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.

************************

70

UNITII
IIII
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:
71

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)

72

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

Frequency response Bode plot Polar plot Determination of closed loop

response from open loop response Correlation between frequency domain
and time domain specifications.

Frequency response
Introduction

There are basically two methods of analysis of systems to determine certain

properties, so that design procedures can be developed, using these properties as
performance measures.

If a time signal like step and ramp are used to excite the system and its time response
is obtained, we call it a time response analysis or time domain analysis.

On the other hand, if a sinusoidal signal of variable frequency is used to excite the
system and the magnitude and phase of the steady state output from the system is
measured, we call it frequency response analysis or frequency domain analysis.

Both the methods have their own advantages and disadvantages.

Frequency Response of a Control System

73

The time domain response of second order system to a unit step unit. In this section
the frequency response of a second order system is obtained and a correlation between
time domain at frequency domain response.
Consider a second order system with the transfer function

The steady state sinusoidal response is obtained by substituting s j in eqn

Normalising the frequency , with respect to natural frequency n by defining a variable

We have

The magnitude and angle of the frequency response is obtained as

The time response for a unit sinusoidal input with frequency is given by

The magnitude and phase of steadystate sinusoidal response for variable frequency can be
plotted. It is to be noted that when

u0

M 1

u 1

1
2
M 0 and

90

180

Let us examine whether the magnitude response given by

74

Monotonically decrease s from 1 to 0 or it attains a maximum value and then decreases to

zero. If it attains the maximum value at some frequency ,its derivative should be zero at that
frequency.

This frequency where the magnitudes become maximum is known as the resonance
frequency. Substituting u u r ,we get the maximum value of the magnitude response. This
value of M M r is known as the resonance peak. Hence

Similarly the phase angle at resonance frequency is given by

75

Fig.(a) Frequency response :Magnitude Vs Normalised frequency

Fig. (b) Frequency response: Phase angle Vs Normalised frequency
One important observation can be made about the resonant peak. It is clear that the resonant
peak depends only on the damping factor . For a given , the resonance frequency is
indicative of the natural frequency of the system. In other words, resonant frequency is a
4
measure of the speed of response of the system since settling time t s
n
Just as peak overshoot and settling time are used as performance measures of a control
system in time domain, the resonant peak and resonance frequency can be used as
performance measures for a control system in frequency domain.
Another important performance measure, in frequency domain, is the bandwidth of the
system. From Fig. we observe that for < 0. 707 and u > ur, the magnitude decreases
1
monotonically. The frequency ub where the magnitude becomes
is known as the cut off
2
1
frequency. At this frequency, the magnitude shill be 20 log
= -3db.
2
As M=1 at u=0 control system are considered as lowpass filters and the frequency at which
the magnitude falls to -3db is known as the bandwidth of the system. This shown in the fig .
frequencies beyond u=ub are greatly attenuated.

An expression for bandwidth can be obtained by equating M=0.707 at u=ub in equation

76

Squaring on both sides of above equation and simplifying, we have

Solving for ub we get

Only positive sign is considered in above equation because ub must be positive and real

The band width is indicative of the noise characteristics of the system. If the bandwidth is
more the system is more susceptible to noise signals. Also for a given , the bandwidth b is
a measured n and hence the speed of response. If the bandwidth is more, the speed of
response is high.
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+jT)q factor

77

Quadratic factor

1 j 2

n n 2

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-loop system becomes unstable
Phase margin:
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:

78

Term

Corner Frequency

Slope
db/dB

Change
in slope

20/jw

---

1/(1+j4w)

wc1=1/4=0.25

-20

-20-20=-40

1/(1+j3w)

wc2=1/3=0.33

-20

-40-20=-60

-20

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

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

79

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

= -90O- tan -1 3w tan -1 4w
When

Frequency in rad /
sec

Phase Angle in degrees

w=0

= -90

w=0.025

= - 99

w=0.25

= -172

w=0.33

= -188

w=3.3

= - 259

w=

= - 270

o
o
o

o
o
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
80

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.

81

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 systems 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.
82

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

Let

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

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

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(jgc)H(jgc) = gc

If an additional phase lag of PM is introduced at this frequency, then the phase

angle G(jgc)H(jgc) 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(jgc)H(jgc)
= 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.

83

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

Time Domain Analysis Vs Frequency Domain Analysis

In the following a comparison of time domain and frequency domain analysis is given.

Variable frequency, sinusoidal signal generators are readily available and precision
measuring instruments are available for measurement of magnitude and phase angle.
The time response for a step input is more difficult to measure with accuracy.

It is easier to obtain the transfer function of a system by a simple frequency domain

test. Obtaining transfer function from the step response is more tedious.

If the system has large time constants, it makes more time to reach steady state at each
frequency of the sinusoidal input. Hence time domain method is preferred over
frequency domain method in such systems.

In order to do a frequency response test on a system, the system has to be isolated and
the sinusoidal signal has to be applied to the system. This may not be possible in
systems which cannot be interrupted. In such cases, a step signal or an impulse signal
may be given to the system to find its transfer function. Hence for systems which
cannot be interrupted, time domain method is more suitable.

The design of a controller is easily done in the frequency domain method than in time
domain method. For a given set of performance measures in frequency domain.

The effect of noise signals can be assessed easily in frequency domain rather than
time domain.

The most important advantage of frequency domain analysis is the ability to obtain
the relative stability of feedback control systems. The Routh Hurwitz criterion is
essentially a time domain method which determines the absolute stability of a system.
The determination of relative stability by Routh Hurwitz criterion is cumbersome.
Nyquist criterion, will not only give stability but also relative stability of the system
without actually finding the roots of the characteristic equation.

Since the time response and frequency response of a system are related through Fourier
transform, the time response can be easily obtained from the frequency response. The
correlation between time and frequency response can be easily established so that the time
domain performance measures can be obtained from the frequency domain specifications and
vice versa.
Correlation between frequency domain and time domain specifications
The time domain specifications are obtained by subjecting the second order system to a unit
step input.
The important time domain specifications are,
84

Peak overshoot

M p e /

1 2

For 0 1

Damped frequency of oscillation

d n 1 2

Settling time
ts

In the frequency domain, the second order system is subjected to constant amplitude, variable
frequency; sinusoidal input and the magnitude and phase response are obtained. The
important frequency domain specifications are

Resonant peak
Mr

For <0.707

2 1 2

Resonance frequency

r n 1 2 2

Band width

b n [1 2 2 4 4 ]
2

1
2

Correlation between M p and M r

Both M p and M r are dependent on the damping factor only and hence they are both
inductive of damping in the system.
Given M p , the resonant peak M r can be evaluated provided is less than 0.707. This
condition is usually satisfied by many practical control system as is seldom greater than
0.707. Thus the resonant peak M r and peak overshoot are well correlated. A plot of M p and
M r with respect to damping factor is given in the fig.

85

The correlation between the time domain specifications of damped natural frequency domain

specification of resonant frequency. For a given damping factor , the ratio r is fixed and
d
given the frequency domain specification, the corresponding time domain specification and
vice versa, can be easily obtained.
Settling Time and Bandwidth
The speed of response is indicated by settling time in time domain. The bandwidth a
frequency domain concept is also indicative of speed of response. Thus we can see that there
is a perfect correlation between time domain and frequency domain performance measures
given one, the other can be obtained easily. Of course one should keep in mind that the
correlation is valid only for < 0.707, which is usually satisfied in many practical control
systems.

*************************

UNIT III
IVIV IV
FREQUENCY RESPONSE
PART -A
1. What is frequency response?
2. What are advantages of frequency response analysis?
3. What are frequency domain specifications?
86

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)

87

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
Nichols 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

Time response Characteristics equation Location of roots in S plane for

stability Routh Hurwitz criterion Root locus construction Effect of pole,
zero addition Gain margin and phase margin Nyquist stability criterion.
88

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 ) s n a1 s n 1 a 2 s n 2 a n 1 s 1 a 0 s 0

Location of roots in S plane for stability

S.no

Nature of closed

Locations of closed

loop poles

loop poles in s-plane

89

Step response

Stability
condition

1.

Real negative i.e in

Absolutely

L.H.S of s-plane

stable

Complex conjugate
2.

with negative real

Absolutely

part i.e in L.H.S of

stable

s-plane
Real positive i.e in

3.

unstable

R.H.S of s-plane
Complex conjugate
with positive real

4.

unstable

part i.e in R.H.S of

S-plane
Non repeated pair

5.

on imaginary axis

Marginally or

without any pole in

critically stable

R.H.S of s-plane
Repeated pair on
imaginary axis
without any pole
in R.H.S of s- plane

6.

Unstable

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

Rouths 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 Rouths 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.
90

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.

Rouths 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 of the 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 an 1s1 a0 s 0

Routh array method

Then add subsequent rows to complete the Routh array
sn :

a2

a4

s n 1 :

a1

a3

a5

n2

b1

b2

b3

s n 3 :

c1

c2

c3

s2 :

s1 :

s0 :

Compute elements for the 3rd row:

b1

1 a3 a2 a1
,
a1

b2

1 a5 a4 a1
,
a1

b3

1 a7 a6 a1
a1

Given the characteristic equation,

a( s) s 6 4s 5 3s 4 2s 3 s 2 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
91

s6 :

1 4

4 0

s :
4

s : 52

12 5 0

s :

s :

s :

s :

s6 :

1 4

s :

4 0

s4 :

52

s :

s2 :

s : 76 15
1

s0 :

12 5 0
4
0

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

Special cases of Rouths 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
92

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

Examine the sign change by taking Lt = 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 Locus = 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 RL going to and meet the RL
at
No. of asymptotes = No. 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)
93

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

1 K

Step 3:

(s z
j 1
n

( s pi )

i 1

Locate the poles and zeros with specific symbols, the root locus begins at the openloop 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

pi z j
i 1

j 1

n94 m

(2k 1)
nm

(k 0,1, 2, n m 1)

Step 8:
The actual point at which the root locus crosses the imaginary axis is readily evaluated
by using Rouths criterion
m

d z
j 1

i 1

1
d pi

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

K1

Nyquist Stability Criteria

(s p )
i

i 1
m

(s z
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

95

where the ai's are the polynomial coefficients, and the coefficients in the rest of the
table are computed using the following pattern

96

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 splane 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
le
Un
sta

97

Problem
Sketch the Nyquist stability plot for a feedback system with the following open-loop transfer
function
G( s) H ( s)

1
s ( s s 1)
2

Solution
For section ab, s = j, : 0
G ( j ) H ( j )

1
j (1 2 j )

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

(ii) = 1 : G(j )H(j ) 1+ j0
(iii) : G(j )H(j ) 0270o

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

G( s) H ( s)

1
0
R3

Section de maps as the complex image of the polar plot as before

98

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
As system gain is increased by a factor 1/a, the open loop G( jw) H ( jw) will increase by
1
a factor a ( ) 1and the system would be driven to instability.Thus, the gain margin is the
a

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
99

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.

100

UNITIV
IVIV 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)
101

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)
7. Sketch the root locus for the open loop transfer function of unity feedback control system
given below: G(S) H(S) =K/S(S+1) (S+2).Also find K of breakaway point. (16)

UNIT V
COMPENSATOR DESIGN
102

Performance criteria Lag, lead and lag-lead networks Compensator design using
bode Plots.

Realisation of Compensators and their Characteristics

The three types of compensators, viz, lead, lag and lag-lead compensators can be
realized by electrical, mechanical.

Pneumatic, hydraulic and other components. We will discuss the realization of these
compensators by electrical RC networks.

We will also obtain their frequency response characteristics which are useful in their
design.

Lead Compensator
We have seen in section that the pole zero configuration of a lead compensator is as
given in fig

The transfer function of a lead compensator is given by,

This transfer function can be realised by an RC network shown in Fig

103

The transfer function of the network is given by

Comparing the above equation

Since 1 it is a lead compensator. Given the time constant and the d.c gain of the
network the three components of RC network namely R1, R2 and C have to determined.
To obtain the frequency response of the network, we have

This means that the network produces an attenuation of , and this d.c attenuation has to be
cancelled by using an amplifier of gain A

.thus the compensating network takes the

form as shown in fig.

Thus, the sinusoidal transfer function of the lead network is given by

104

Under the steady conditions the output of this network leads the input and hence this network
is known as a lead network.
The bode plot of equation is shown in the fig.

G( j ) is given by

The angle is a function of and it will be a maximum when tan is a maximum with
respect to

The above equation can be written as

105

 m is the geometric mean of the two corner frequencies

Substituting the value of m we get

From the right angled triangle

We have

Solving for ,we get

Also, the gain at m is obtained as

Substituting m

we have

106

1
1
and

The lead compensator is required to provide the necessary angle at the gain crossover
frequency, to obtain the required phase margin. Knowledge the value of required angle m to
be provided by the network parameter can be obtained .the value of required to get a
maximum phase lead m is tabulated

20

0.49

40

0.217

60

0.072

80

0.0076

From the table it can be seen that for obtaining larger phase angles, has to be reduced to
very low value. For very low values of . the pole of the lead network is located for away
from the j ,axis and from the Bode magnitude diagram it can be seen that the network has a
1
db at high frequencies. Than the high frequency noise signals are

amplified and signal to noise ratio becomes very poor. Hence it snot desirable to have low
values . A value of = 0.1 is considered to be suitable. When large phase angles are
desired, two lead networks in cascade may be used.

large gain of 20 log

Lag Compensator
The transfer function f a lag compensator is given by,

1 ,Ensures that the pole is to the right of the zero. The pole zero configurations is given
in fig

The RC network realization is shown in the fig

107

Where

The sinusoidal transfer function is given by

The d.c gain of this network is unity and hence no d.c amplifier is required as in the case of a
lead network.

108

Under steady state conditions the output of the network lags the input and therefore this
network is known as a lag network. The maximum phase lag m is again obtained from the
equation and frequency at which it occurs, replacing by .At high frequencies the network
produces an attenuation of 20 log and hence signal to noise ratio is improved by this
network.
Lag-Lead Compensator
The general form of the transfer function of lag lead compensator is

Transfer function given by eqn. can be realised by an RC network as shown in Fig

The transfer function of this network can be derived as

109

Comparing the equation

And

This means that and cannot be chosen independently.

Therefore we can write Gc(s) as

The Bode plot is shown in fig

110

Cascade Compensation in Frequency Domain Using Bode Plots

In order to design a compensator using Bode plot, the specifications are to be given in
frequency domain.

If specifications are given in time domain, they can be easily converted into frequency
domain using time and frequency domain correlations.

The frequency domain Specifications are usually given by:

Phase margin m or peak resonance M r which indicates relative stability.
Bandwidth b or resonance frequency r , which indicates rise time or
settling time
Error constants which indicates steady state errors.

The design is carried out in the frequency domain using Bode plots; the time response
of the compensated system must be obtained to check the time response
specifications.

This is necessary because the time and frequency domain specifications arc correlated
under the assumption that the compensated system is of second order, or it has a pair
of dominant closed loop poles.

Based on the time response characteristics. The design may be altered suitably.

The frequency domain method of design is easy to apply but the time response
specifications cannot be directly controlled in this method.

Design of Lead Compensator

The procedure for the design of a lead compensator is developed in the following.
1. The system gain K is adjusted to satisfy the steady stale error criterion. as specified by
the appropriate error constants K p, K v or K a or the steady state error ess.

111

2. The Bode plot of the system with the desired K is plotted. The phase margin of the
uncompensated system is read from the graph ( 1 ). If this is satisfactory no
compensator is necessary. If the phase margin falls short of the desired phase margin,
additional phase lead m must be provided by a lead compensator at the gain cross
over frequency.
3. Since maximum phase lead. m , of a RC lead network occurs at a frequency m ,
which is the geometric mean of the two corner frequencies to, and et, of the network,
it is desirable to have this maximum phase lead occur at the gain cross over frequency
of the uncompensated system. Thus the two corner frequencies of the lead
compensator must be located on either side of the gain crossover frequency gc1 .But
the gain curve of the compensator also effects the gain curve of the uncompensated
system. It lifts up the gain curve around the gain cross over frequency and hence the
compensated gain curve crosses 0 db line at a slightly higher frequency gc 2 . This is
shown in Fig. for a type one, second order system.

The frequency m ,at which maximum phase lead occurs in a phase lead network is
made to coincide with the gain cross over frequency gc1 , of the uncompensated
system. The effect of this is twofold.
The gain cross over frequency moves to gc 2 which is higher than gc1 . Hence the
phase margin obtained after the compensation is not equal to 1 m as desired
The new phase margin has to be read at the new gain cross over frequency gc 2 .and is
equal to 2 3 as shown in the fig. Obviously 2 1 and 3 m and hence the
phase margin obtained will fall short of the desired value.
It is clear from that if m of lead network is made to coincide with gc1 ,desired phase
margin is not obtained for the compensated network. Thus it is desirable to make
m coincide with gc 2 so that maximum phase lead occurs at the new gain crossover
frequency. But this new gain crossover frequency is not known before hand.
The above difficulty is overcome as follows. The phase lead network provides a gain
1
of 10 log at m .if m is to be coincide with the new gain cross over frequency

gc 2 , gc 2 must be chosen such that the uncompensated system has a gain of

10 log

at this frequency

112

 The maximum phase angle to be provided by the lead network is obtained as

follows:
At the gain crossover frequency of the uncompensated system. The system has a
phase margin of 1 which is less than the desired phase margin d . Hence an
additional phase lead of d 1 , has to be provided by the lead network. But
thesystem will not be providing a phase lead of 1 after compensation because
the gain crossover is going to shift to a new higher value gc 2 . at which the
phase angle provided by the system will be 2 . as shown in Fig. which is less
than 1 . Also gc 2 and hence 2 is not known to start with. In order to
compensate this unknown shortfall, an additional phase angle ranging from
5to 20

is to be provided. Thus the maximum phase lead m to be provided by

the lead network is taken to be,

may be around 5 if the slope of the uncompensated system at gc1 is -40

db/dec or less and may go to 15 to 20 if it is 60 db/dec.
The RC phase lead network is designed as follows:
Maximum phase lead network provided is given by,

Then is calculated from,

We have to make m gc 2 , where to gc 2 is the frequency at which.

113

Knowing , we can locate gc 2 from the magnitude plot.

From these values, the corner frequencies of lead network are obtained as

The transfer function of the lead network is

It is observed that the gain cross over frequency is increased after compensation.
Gain cross over frequency can be taken as a measure of the bandwidth of the
system. Since the bandwidth is increased by the lead compensation, Speed of
response is improved. The actual bandwidth can be obtained by transferring the
data from Bode plot to Nichols chart.

1
.the noise

frequencies' are amplified and hence the signal to noise ratio deteorates.

Since the gain at higher frequencies is increased by 20 log

The values of R, R2 and Cane be calculated using,

Since one of the three elements R1, R2and C can be chosen arbitrarily, as there are only two
equations, impedance level required or the cost of the compensator maybe used to select this
element.
The design procedure may be summarized as follows.
(i) Select the gain K to satisfy steady state error requirements.
(ii) Using this value of K draw the Bode magnitude and phase plots of the uncompensated
system. Obtain the phase margin of the uncompensated system. Let this be 1
(iii)Determine the phase lead required.

Where
114

 d - Desired phase margin

1 -Available phase margin
-Additional phase angle required to compensate the reduction in
phase angle due to increase in gain crossover frequency after
compensation.

(iv)Determine of the RC network as

If m required is more than 60o it is recommended to use two RC networks in

cascade each providing an angle of

m
2

(v)Find the gain of the phase lead network at m , which is 10 log

, Locate the

frequency at which the gain of the uncompensated system is 10 log

.Let this

frequency be gc 2 .Make gc 2 = m .Read the phase margin 2 , provided by the

system at this frequency m .If 1 2 is ,increase and recalculate

m , and m until 1 2 .
(vi) Compute the two corner frequencies of the lead network

The transfer function of lead network

(vii) Draw the bode plot of the compensated system and check the phase margin.
If it falls short of the desired phase margin increase and redesign.
(viii) Obtain the elements of RC phase lead network
Design of Lag Compensator
The procedure for design of a lag compensator is developed in the following:
115

1. The gain constant is set as per the steady state error requirements.
2. Bode plot is drawn. If the phase margin obtained is adequate no other compensation is
required. If phase margin is not adequate, locate a frequency gc 2 , at which the desired phase
margin is available.
3. The gain of the uncompensated system at this frequency is measured. This gain will be
positive for type 1 or type 0 systems.
4. The attenuation characteristic of phase lag network at high frequencies is used to reduce
the gain of the uncompensated system to zero at gc 2 , so that this frequency becomes the gain
crossover frequency of the compensated system, as shown in Fig. The phase lag introduced
by the lag network is detrimental to the system. Hence the frequency at which maximum
phase lag occurs for the phase lag network must be located for away from the gain crossover

5. The corner frequencies of the lag network are located far to the left of the new gain
cross over frequency gc 2 ' so that the phase lag provided by the lag network does not
appreciably affect the phase of the system at this frequency. However a small negative
angle will be added to the system phase at gc 2 , as shown in Fig. To compensate this, the
phase margin to be provided is taken to be slightly higher than the desired phase margin.
Thus the phase margin is taken to be

m d
is taken to be in the range of 50o to 15o. If the corner frequencies of lag network are
closer to the new gain cross over frequency. higher value of is chosen. The upper

116

corner frequency 2 , of the lag network is placed, usually, one octave to 1 decade
below the new gain cross over frequency gc 2 .
The design steps are summerised as follows
1. Choose the gain K to satisfy steady state error specifications.
2. Draw the Bode plot with this value of K. If the phase margin of this system is
inadequate, design a lag compensator.
3. The phase margin to be provided is given by

Where , ranges between 5 to 15.

4. Locate the frequency gc 2 , at which the phase margin is available on the gain plot
of G(s). Find the magnitude of G(s) at this frequency.
5. The attenuation to be provided by the lag network is given by,

Find the value of from this equation.

6. The upper corner frequency, 2 of the lag compensator is chosen to be between

gc 2
2

gc 2
10

to

.Thus

7. With and known, the phase lag network can be designed. The transfer function
of the Lag network is,

8. Draw the Bode plot of the compensated system to check the specifications.
If necessary, choose a different value of and redesign to satisfy the specifications.
Design of Lag-Lead Compensator

117

A lead compensator increases the gain cross over frequency and hence the bandwidth
is also increased. For higher order systems with large error constants, the phase lead
required may be large, which in turn results in larger bandwidth.

This may be undesirable due to noise considerations. On the other hand a lag
compensator results in a lower bandwidth and generally a sluggish system.

To satisfy- the additional requirement on bandwidth, a lag or lead compensator alone

may not be satisfactory.

Hence a lag-lead compensator is used. The procedure for the design of lag-lead
compensator is as follows:
I. The required error constant is satisfied by choosing a proper gain. The phase margin
and the bandwidth are obtained from the Nichols chart. If the phase margin falls short and
the bandwidth is smaller than the desired value, a lead compensator is designed.
2. If the bandwidth is larger than the desired value, lag compensator is attempted, if the
phase margin desired is available at any frequency. If this results in a lower bandwidth
than the desired bandwidth, a lag lead compensator is to be designed.
3. To design a lag-lead compensator, we start with the design of a lag compensator. A lag
compensator is designed to partially satisfy the phase margin requirement. It means that
the gain cross over frequency gc 2 is chosen to be higher than that to be used, if full lag
compensation is designed. This ensures that the bandwidth is not reduced excessively by
the lag compensator.
4. The and gc 2 are known and the comer frequencies of the lag compensator can be
obtained. Since

the frequency at which 20 log is available on the magnitude

plot is the new gain crossover frequency gc 3 . The maximum phase lead available is,

With these parameters, a lead compensator is designed. Here, for the lead network,

The log magnitude vs phase angle curve of the lag-lead compensated network is drawn
on the Nichols chart, from which the bandwidth can be obtained.
118

***********************

UNIT V

COMPENSATOR DESIGN
PART-A
1. What are the two methods of designing a control system?
2. What is the time domain specifications needed to design a control system?
3. What are the frequency domain specifications needed to design a control system?
4. What is compensation?
5. What is compensator?
6. What are the different types of compensator available?
7. What are the two types of compensation schemes?
8. What is feedback compensation?
9. Why compensation is necessary in feedback control system?
10. What is lag compensation?
PART -B
1. What is compensation? Why it is need for control system? Explain the types of
Compensation? What is an importance of compensation?
2.Realise the basic compensators using electrical network and obtain the transfer function.
3. Design suitable lead compensators for a system unity feedback and having open loop
transfer function G(S)= K/ S(S+1) to meet the specifications.(i) The phase margin of the
system 45, (ii) Steady state error for a unit ramp input 1/15, (iii) The gain cross over
frequency of the system must be less than 7.5 rad/sec.
4. A unity feed back system has an open loop transfer function G(S) = K/ S(S+1) (0.2S+1).
Design a suitable phase lag compensators to achieve following specifications Kv= 8
and Phase margin 40 deg with usual notation.
5. Explain the procedure for lead compensation and lag compensation?
6. Explain the design procedure for lag- lead compensation
7. Consider a type 1 unity feed back 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.
8. Using Electrical lead network derive the transfer function.
9. Using Electrical lag network derive the transfer function
119

10. Using Electrical lag-lead network derive the transfer function

11. Design a lead compensator for a unity feedback system with open loop transfer function
G(S) = K/ S(S+1) (S+5) to satisfy the following specifications (i) Kv > 50 (ii) Phase
Margin
is 20 .
12. Design a lead compensator for G(S) =K / S2 (0.2S+1) to meet the following
Specifications (i) Acceleration ka=10; (ii) P.M=35.
13. Design a Lag compensator for the unity feedback system whose closed loop transfer
function C(s) / R(s) = K / (s (s+4) (s+80) + K) is to meet the following specifications
P.M 33 . And Kv 30.
14. Consider the unity feed back system whose OLTF is G(s) = K / s (s+3)(s+6). Design a
lag-lead compensator to meet the following specifications.
(i) Velocity error constant, Kv = 80, (ii) P.M _ 35 (16)
15. A unity feedback system has an OLTF G(s) = K / s(s+2)(s+60). Design a Lead-Lag
compensator is to meet the following specifications.
(i) P.M is atleast 40 , (ii) Steady state error for ramp input 0.04 rad.

*********************************END*********************************

120

121

122

123

124