Control Systems Lab Manual 18EEL66 Student
Control Systems Lab Manual 18EEL66 Student
Prepared By
Mrs. Ruma Sinha
GLOBAL ACADEMY OF TECHNOLOGY
VISION
Become a premier institution imparting quality education in engineering and
management to meet the changing needs of society
MISSION
M1. Create environment conducive for continuous learning through quality
teaching and learning processes supported by modern infrastructure
i
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
VISION
To be a foremost education program in the field of Electrical & Electronics
Engineering to train the students to face global challenges and succeed in
their careers.
MISSION
M1: Strengthen learning environment that facilitates quality education in the
field of Electrical & Electronics Engineering.
M2: Enhance the Industry Institute interaction continuously to enable
students work on real time engineering problems.
M3: Improve the quality of value added programs to enhance the technical
and intellectual capabilities of students ensuring their success in competitive
examinations for higher studies, employment, and research.
ii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
iii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
iv
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
v
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
Regulations Governing
THE DEGREE OF BACHELOR OF ENGINEERING
ATTENDANCE REQUIREMENT
Each semester is considered as a unit and the candidate has to put in a minimum
attendance of 85% in each subject with a provision of condonation of 10% of the
attendance by the Vice-Chancellor on the specific recommendation of the Principal
of the college where the candidate is studying, showing some reasonable cause
such as medical grounds, participation in University level sports, cultural activities,
seminars, workshops and paper presentation, etc.
The basis for the calculation of the attendance shall be the period prescribed by
the University by its calendar of events. For the first semester students, the same is
reckoned from the date of admission to the course as per CET allotment.
The students shall be informed about their attendance position periodically by the
colleges so that the students shall be cautioned to make up the shortage.
A Candidate having shortage of attendance in one or more subjects shall have to
repeat the whole semester and such candidates shall not be permitted to take
admission to next higher semester. Such students shall take readmission to the
same semester in the subsequent academic year.
Course objectives
Upon completion of this course, students are expected to:
1. Study speed torque characteristics of AC, DC servomotor and Synchro pair.
2. Learn to simulate system in time domain and analyze the effect of addition of
poles, zeros, and loop gain and pole location on stability.
3. Understand frequency response characteristics of second order system to
determine frequency domain specification.
4. Gain practical Knowledge to design lag, lead, lag–lead compensating network and
PID controller.
5. Compare Bode plot, Nyquist plot and Root locus with reference to stability.
vii
Dept. of EEE 8 Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
Course outcome
Cognitive No. of
COs COURSE OUTCOMES
Level sessions
viii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
SEMESTER – VI B.E ELECTRICAL AND ELECTRONICS ENGINEERING (EEE)
CONTROL SYSTEMS LABORATORY
Subject Code 18EEL66 IA Marks 40
Number of Practical Hours/Week 03 Exam Hours 03
Total Number of Practical Hours 42 Exam Marks 60
Credits - 02
SL No Experiments Page No
xii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
Experiments 8 to 11 must be done using standard simulation package
(a) To simulate a typical second order system and determine step response
and evaluate time response specifications.
8.
(b) To evaluate the effect of additional poles and zeros on time response
of second order system.
(c) To evaluate the effect of pole location on stability.
(a) To simulate a D.C. Position control system and obtain its step response.
(b) To verify the effect of input waveform, loop gain and system type on steady
9.
state errors.
(c) To perform trade-off study for lead compensator.
(d) To design PI controller and study its effect on steady state error
(a) To examine the relationship between open-loop frequency response and
10. stability, open loop frequency and closed loop transient response
(b) To study the effect of open loop gain on transient response of closed
loop system using root locus
(a) To study the effect of open loop poles and zeros on root locus contour
11. (b) Comparative study of Bode, Nyquist and root locus with respect to stability
Revised Bloom’s
L3 – Applying, L4 – Analysing, L5 – Evaluating, L2 –
Taxonomy level Understanding
Graduates Attribute (As per NBA)
Engineering Knowledge Problem Analysis Individual and Team work
Communication
xi
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
Conduct of Practical Examination
1. All laboratory experiments are to be included for practical examination.
2. Breakup of marks and the instructions printed on the cover page of
answer script to be strictly adhered by the examiners.
3. Students can pick one experiment from the questions lot prepared by the
examiners.
4. Change of experiment is allowed only once and 20% Marks allotted to the
procedure part to be made zero.
1 Write-Up 15
2 Conduction 70
3 Viva Voce 15
TOTAL 100
2017 SCHEME
(6th SEMESTER)
1 Observation Book 10
2 Record 15
3 Viva 5
TOTAL 30
xii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
Average of Weekly
1 30
Entries
Internal Assessment
2 10
Reduced To
TOTAL 40
xi
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
Cycles of Experiments
Cycle - 1
1.a) Experiment to draw the speed torque characteristic of a two phase AC servomotor.
1.b) Experiment to draw the speed torque characteristic of a DC servomotor.
2.a) To design a RC lag compensating network for the given specifications., viz., the
maximum phase lag and the frequency at which it occurs, and to obtain its frequency
response.
2.b) To determine experimentally the transfer function of the lag compensating network.
3.a) To design a passive RC lead compensating network for the given specifications., viz., the
maximum phase lead and the frequency at which it occurs, and to obtain its frequency
response.
3.b) To determine experimentally the transfer function of the lead compensating network.
4. Simulation of a typical second order system and determination of step response
and evaluation of time domain specifications.
5. a) Simulate a DC Position control system using MATLAB and obtain its step response.
5. b) To verify the effect of input waveform, loop gain and type of system on steady state
error.
6. To perform a trade-off study for lead compensation.
xii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
Cycle 2
1. Experiment to draw and study of Synchro pair characteristics.
2. Experiment to draw the frequency response characteristic of a given lag-
lead compensating network.
3. To determine the frequency response of a second order system and evaluation of
frequency domain specifications experimentally.
4.a) To determine the effect of open loop poles and zeroes on root locus contour.
4.b) To estimate the effect of open loop gain on the transient response of closed
loop system by using Root locus
5. Comparative study of Bode, Nyquist and Root locus with respect to stability.
6. To study the effect of P, PI, PD and PID controller on the step response of a
xiii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
List of Experiments
SL No Experiments
Using MATLAB/SCILAB a) Simulation of a typical second order system and
determination of step response and evaluation of time- domain specifications, b)
1. Evaluation of the effect of additional poles and zeroes on time response of second
order system, c) Evaluation of effect of pole location on stability, d) Effect of loop
gain of a negative feedback system on stability
(a) To design a passive RC lead compensating network for the given specifications,
viz., the maximum phase lead and the frequency at which it occurs and to obtain its
2. frequency response.
(b) To determine experimentally the transfer function of the lead compensating
network.
(a) To design RC lag compensating network for the given specifications., viz., the
maximum phase lag and the frequency at which it occurs, and to obtain its
3. frequency response.
(b) To determine experimentally the transfer function of the lag compensating
network.
Experiment to draw the frequency response characteristic of a given lag- lead
4.
compensating network.
To study the effect of P, PI, PD and PID controller on the step response of a
5. feedback control system (using control engineering trainer/process control
simulator). Verify the same by simulation.
a) Experiment to draw the speed – torque characteristic of a two - phase A.C.
servomotor.
6.
b) Experiment to draw speed torque characteristic of a D.C. servomotor.
xii
Dept. of EEE Control Systems Lab Manual / 17EEL67
GLOBAL ACADEMY OF TECHNOLOGY
its step response, b) To verify the effect of the input wave form, loop gain system
type on steady state errors, c) To perform a trade-off study for lead compensation,
d) To design a PI controller and study its effect on steady state error.
Using MATLAB/SCILAB a) To examine the relationships between open-loop
frequency response and stability , open loop frequency and closed loop transient
9.
response, b) To study the effect of addition closed loop poles and zeroes on the
closed loop transient response
Using MATLAB/SCILAB a) Effect of open loop and zeroes on root locus contour,
b) To estimate the effect of open loop gain on the transient response of closed loop
10.
system by using Root locus c) Comparative study of Bode, Nyquist and Root locus
with respect to Stability.
11. Experiment to draw to syncro pair characteristics.
12. Using MATLAB, represent the given system in state space.
13. Using MATLAB, check the Controllability and Observability of a given system.
xv
Dept. of EEE Control Systems Lab Manual / 17EEL67
INTRODUCTION
Control Engineering is of significant interest in most areas of industry - new and established.
Control System applications are assuming an increasingly important role in the developing
world, making it crucial for the students to be exposed to contemporary control system
equipment in a realistic manner, in order to connect theoretical material taught in lecture
courses with the realities of physical hardware and simulation through ‘higher-level’ technical
computing language. Control engineering is based on the foundations of feedback theory and
linear system analysis. Accordingly, control engineering is not limited to any engineering
discipline but is applicable to aeronautical, chemical, mechanical, environmental, civil, and
electrical engineering. A control system is an interconnection of components forming a system
configuration that will provide a desired system response. The basis for analysis of a system is
the foundation provided by linear system, which assumes a cause-effect relationship for the
components of a system. A system is a collection of elements or components that are organized
for a common purpose and a control system is a system of devices or set of devices, that manages,
commands, directs or regulates the behavior of other devices or systems to achieve desired
results.
Introduction to MATLAB:
The name MATLAB stands for MATrix LABoratory. MATLAB is a high-performance language for
technical computing. It integrates computation, visualization, and programming environment.
Furthermore, MATLAB is a modern programming language environment: it has sophisticated
data structures, contains built-in editing and debugging tools, and supports object-oriented
programming. These factors make MATLAB an excellent tool for teaching and research.
MATLAB has many advantages compared to conventional computer languages (e.g., C,
FORTRAN) for solving technical problems. MATLAB is an interactive system whose basic data
element is an array that does not require dimensioning. The software package has been
commercially available since 1984 and is now considered as a standard tool at most universities
and industries worldwide. It has powerful built-in routines that enable a very wide variety of
computations. It also has easy to use graphics commands that make the visualization of results
immediately available. Specific applications are collected in packages referred to as toolbox.
GLOBAL ACADEMY OF TECHNOLOGY
There are toolboxes for signal processing, symbolic computation, control theory, simulation,
optimization, and several other fields of applied science and engineering.
Control Systems Lab manual 18EEL66
EXPERIMENT: 01 (a)
INTRODUCTION:
A servo motor is a rotary actuator that is designed for precise position control for closed loop
application. It consists of an electric motor, a feedback device, and a controller. They are able
to accommodate complex motion patterns and profiles better than any other type of motor.
Although they are small in size they have higher power capability and are extremely energy
efficient. Unlike large industrial motors, a servo motor is not used for continuous energy
conversion.
There are two types of servo motors, AC servos and DC servos. The main difference between
the two motors is their source of power. AC servo motors rely on an electric outlet, rather
than batteries like DC servo motors. While DC servo motor performance is dependent only
on voltage, AC servo motors are dependent on both frequency and voltage. Because of the
complexity of the power supply, AC servo motors can handle high surges, which is why they
are often used in industrial machinery.
AIM: To conduct a load test on AC servo motor and to obtain its speed torque
characteristic.
APPARATUS REQUIRED:
1. AC servomotor module with loading arrangement,
2. Digital multimeter
3. Connecting wires.
CIRCUIT DIAGRAM:
Reference
Winding
+ +
2
3
- 0
DC AC servo
Eb Motor V
Motor
Control
A
Winding
- Control C
Voltage
THEORY:
AC servo motors are basically two-phase squirrel cage induction motors and are used for low
power applications. Nowadays, three phase squirrel cage induction motors have been
modified such that they can be used in high power servo systems.
PROCEDURE:
1. Connections are made as shown in the circuit diagram.
2. Ensure that the potentiometer and control voltage autotransformer are in min position
and load switch is in OFF position.
3. Switch on the power supply.
4. Control voltage switch is put on and the control voltage is adjusted to the rated value 220
V.
5. The load switch is put on and the readings are noted.
6. Loading is done by varying the potentiometer till the rated armature current of PMDC
machine is reached (250 mA)
7. Readings at different steps are noted down.
8. Adjust the potentiometer to its original position and switch off the load switch.
9. Bring autotransformer to zero position and switch off the supply.
10. Repeat the above procedure by setting the control voltage to 200V&175 V.
11. Measure the armature resistance of the PMDC machine using multimeter.
TABULAR COLUMN
Speed Torque
Vc (volts) Ia (Amp) Va (volts) P (watts)
(rpm) (Gm-cm)
CALCULATION:
P P 10 5 60
Power P VaIa and Torque T
2 N 9.81
With Vc = 220 V,
P P 10 5 60
i. P V a I a = 0.205 watt, T =112.8 gm-cm
2 N 9.81
Expected graph
N
(rpm)
Torque
(gm-cm)
INFERENCE:
For a two phase servomotor with variation of current, variation of speed, torque are
obtained for different level of control winding voltage. Speed torque characteristics of two phase
servomotor is drawn.
EXPERIMENT: 01 (b)
INTRODUCTION:
Generally, the DC servo motor includes a DC source separately in the field of the armature
winding. The motor can be controlled either by managing the field current otherwise the armature
current. The armature control has some benefits compare with field control. Similarly, field control
has come benefits compare with armature control. The controlling of this motor can be done
based on the application used. This motor offers a quick and accurate response to begin or end
command signals because of the small armature inductive reactance. These motors are utilized in
several devices and numerically controlled equipment.
AIM: To conduct a load test on DC servo motor and to obtain its speed torque characteristic.
APPARATUS REQUIRED:
1. DC servomotor module with loading arrangement,
2. Digital multimeter
3. Connecting wires.
CIRCUIT DIAGRAM
+ Mechanical
A
Load
THEORY:
A DC servo motor consists of a small DC motor, feedback potentiometer, gearbox, motor drive
electronic circuit and electronic feedback control loop. The stator of the motor consists of a
cylindrical frame and the magnet is attached to the inner periphery of the frame. A brush is built
with an armature coil that supplies the current to the commutator. At the back of the shaft, a
detector is built into the rotor in order to detect the rotation speed. With this construction, it is
simple to design a controller using simple circuitry because the torque is proportional to the
amount of current flow through the armature. And also the instantaneous polarity of the control
voltage decides the direction of torque developed by the motor. DC servo motor provides very
accurate and also fast respond to start or stop command signals due to the low armature
inductive reactance. DC servo motors are used in similar equipment and computerized
numerically controlled machines.
DC Servomotor AC Servomotor
It has more stability problems. It has less stability problems
It requires frequent maintenance due to the It requires less maintenance due to the
presence of commutator. absence of commutator.
The life of DC servo motor depends on the The life of AC servo motor depends on
life on brush life bearing life
These motors are used for high power These motors are used for low power
applications. applications.
PROCEDURE:
TABULAR COLUMN
CALCULATION:
2 N
T= (F1 – F2) x r, r=radius of brake drum, P = T x x10-5
60
With Va = 23.2 V,
2 N
i. T= (F1 – F2) x r = 29.595 gm-cm, P T 10 5 = 0.0425 watt
60
EXPECTED GRAPH:
N
(rpm)
Torque
INFERENCE:
For a DC servomotor with variation of current, variation of speed and torque are obtained for
different level of control winding voltage. Speed torque characteristics of two phase servomotor
is drawn.
EXPERIMENT: 02
Synchro pair characteristics
INTRODUCTION:
The Synchro is a type of transducer which transforms the angular position of the shaft into an electric
signal. It is used as an error detector and as a rotary position sensor. The error occurs in the system
because of the misalignment of the shaft. The transmitter and the control transformer are the two main
parts of the synchro.
The synchro system is of two types. They are
1. Torque Transmission type Synchro.
2. Control Type Synchro.
1. Torque Transmission type Synchros:
This type of synchros has small output torque, and hence they are used for running the very light load like a
pointer. The control type Synchro is used for driving the large loads.
2. Control Type Synchros System:
The control type synchros is used for error detection in positional control systems. Their systems consist two
units. They are
(i) Synchro Transmitter
(ii) Synchro receiver
APPARATUS REQUIRED:
1. Synchro transmitter and receiver module
2. Connecting wires
diagram:
DPM
Meter
input
Transmitter Receiver
0 S1 S1 0
S2 S2
90 90
270 270
S3 S3
180 180
R1 SW1 R2 R1 SW1 R2
THEORY:
A Synchro Transmitter Receiver Trainer is an electromagnetic transducer commonly used to convert an
angular position of a shaft into an electric signal. It is usually called a synchro transmitter. Its construction
is similar to that of a three phase alternator. The stator of the synchros is made of steel for reducing the
iron losses. The stator is slotted for housing the three phase windings. The axis of the stator winding is
kept 120º apart from each other and is Y connected as shown in fig(a).
Fig(a) Fig(b)
The rotor is a dumb bell construction and wound with a concentric coil. The system set up is made up of
synchro transmitter and synchro receiver on a single rigid base provided with suitable switches. The system
also contains a step down transformer for providing excitation to the rotors. Suitable test points for rotor (R1
and R2) and stator (S1, S2 and S3) as shown in fig(b). The circuit diagram of synchro transmitter and receiver
as shown in fig(c),
Fig(c)
Advantages of synchro
The operating speed of synchro is very high, they have full 360 degree rotation without
electrical break and they are highly reliable.
Disadvantages of synchro
They have full range in order of 30 degrees linearity is not much better than modern
potentiometers.
PROCEDURE:
TABULAR COLUMN:
SL# Stator terminal Voltage
Rotor position
VS1S2 VS2S3 VS3S1
(degrees)
(volts) (volts) (volts)
1. 0 38.3 70.1 31.5
2. 20 15.5 67 51
3. 40 8 55.3 65
4 60 31.1 38.3 70
5 80 50.1 16.3 67.2
6 100 63.7 7.3 56
7 120 69.5 29.5 39.5
8 140 67.2 50.8 16.4
9 160 56.2 63.9 7.2
10 180 38.2 69.9 30.4
11 200 16.1 67.0 50.5
12 220 7.9 55.7 64.2
13 240 23.8 38 70.2
14 260 31.2 15.9 67.3
Connection Diagram:
DPM
Met
er
inp
ut
Transmitter Receiver
0 S1 S1 0
S2 S2
90 90
270 270
S3 S3
180 180
R1 SW1 R2 R1 SW1 R2
PROCEDURE:
1. Connect mains supply cable
2. Connect S1, S2 & S3 terminals of transmitter to the S1, S2 and S3 terminals of receiver
respectively with patch cords.
3. Switch on SW1 and SW2 and also the main supply.
4. Move the pointer of rotor position of synchro transmitter and observe the
receiver’s rotor position.
5. Note down the transmitter’s angular position and receiver’s angular position and plot
a graph.
TABULAR COLUMN:
Angular position(degrees)
SL#
Synchro transmitter Synchro receiver
1 0 -2
2 20 17
3 40 36
4 60 56
5 80 78
6 100 98
7 120 120
8 140 139
9 160 159
10 180 179
11 200 198
12 220 218
13 240 238
14 260 258
15 280 275
16 300 299
17 320 319
18 340 340
19 360 360
EXPECTED GRAPH:
VS3S1
VS2S3 VS1S2
R
M
S
V
O
L
T
A
G
E
0
Angular position in degree
Receiver’s
Angular
position
ACTUAL GRAPH:
The properties of synchro are studied and graphs are plotted. Because of
transformer action, voltages are induced in each of the stator coil.
EXPERIMENT: 03
FREQUENCY RESPONSE OF SECOND ORDER SYSTEMS
INTRODUCTION:
The response of a system can be partitioned into both the transient response and the
steady state response. We can find the transient response by using Fourier integrals.
The steady state response of a system for an input sinusoidal signal is known as
the frequency response.
The order of the system is defined by the number of independent energy storage
elements in the system, and intuitively by the highest order of the linear differential
equation that describes the system. In a transfer function representation, the order is
the highest exponent in the transfer function. The order of the characteristic equation
can also be referred as the order of the system. For a second order system the order of
characteristic equation is 2.
AIM:
a. To determine the frequency response characteristics of a second order system to a
sinusoidal input and determine its frequency response specifications: resonant
peak M r , resonant frequency r and the bandwidth (BW or b ) using hardware
component.
b. Verify the result using simulation package
APPARATUS REQUIRED:
1. Decade resistance box (1)
2. Decade inductance box (1)
3. Decade capacitance box (1)
4. Signal generator (1)
5. Oscilloscope. (1)
THEORY:
A SECOND ORDER SYSTEM SET UP: A series RLC circuit can be used to study the
frequency response of a second order system. The voltage response across the
capacitor to a step input in a series RLC circuit is analogous to step response of any
second order system. Therefore, while using a series RLC circuit to study the frequency
response characteristics, the voltage across the capacitor is to be used as the output
and the input to the circuit can be a step voltage or a square wave of suitable frequency
from a signal generator.
Vin Vout
Signal
Generator
The transfer function of the series RLC network with the capacitor voltage as output can
be obtained as,
1
V o
( s) LC
V i
( s)
s
2 R
s
1
L LC
It is similar to the standard form of closed loop transfer function of a second order system given
by:
C (s) 2 1
n , where ωn = natural frequency = and
R( s) s2 2 n
2
n LC
R C
ξ = damping ratio = .
2 L
With proper selection of L & C (in order to have suitable n) and by varying the value
of resistance R different damping ratios can be obtained.
The frequency domain specifications are resonant peak, resonant frequency and
bandwidth, cut off frequency, phase margin, gain margin, gain crossover frequency,
phase crossover frequency.
Resonant frequency - It is the frequency at which the magnitude of the frequency
response has peak value for the first time.
Resonant peak - It is the peak (maximum) value of the magnitude. A system, with larger
resonant peak will exhibit a larger overshoot.
Bandwidth – In feedback control system, the range of frequencies over which, the
Bode Diagram
10
Mr
0
-20 BW
-30
ωr ωc
-40
0
-45
Phase (deg)
-90
-135
-180
0 1 2
10 10 10
Frequency (rad/sec)
PROCEDURE:
1) The RLC series circuit is connected as shown in the figure using decade boxes.
2) The input to the circuit is to be taken from a signal generator so that a variable
frequency sinusoidal input can be given and the output (voltage across the
capacitor) is to be observed in a CRO. These devices are properly connected.
3) Set L=10 mH and C= 0.1 F so that n=31662.8 radians / sec.
4) Now set R in order to have an under damped response with =0.2
5) Switch on the signal generator and give a sine wave input with a suitable
magnitude, say 1V.
6) Now slowly vary the frequency from a low value to a large value well beyond n
8) Draw the frequency response graph to a proper scale on the graph sheet and from the
graph determine the frequency response specifications resonant peak M p , resonant
9) Calculate the theoretical values using the expressions given in the previous section and
compare with the experimental results.
10) Repeat the experiments for two other values of damping ratios by changing the
value of resistance R.
TABULAR COLUMN:
Input magnitude Vi = 3 Volts.
Frequency Output Phase
Frequency in Gain
f magnitude Angle
(Radians/sec) (dB)
(Hz (V0max) (Degrees)
)
235 1476.55 3 0 0
CALCULATION:
C
0.5R 0.5 R 3.162 10 3
L
For ξ = 0.2, R = 126.5 ohm
RESULT:
Performance Value from Value from
parameters experiment simulation
Resonant peak 4.0823 dB
Resonant frequency 25132.7 rad/sec
Band width 44000 dB
APPARATUS REQUIRED:
1. CRO
2. Signal generator
3. decade resistance and capacitance box
4. MATLAB Software package
CIRCUIT DIAGRAM:
Signal
Generator E1 E2
THEORY:
The lead compensator is an electrical network which produces a sinusoidal output having phase lead
DESIGN:
G(s) = E (s)
2 R 2
1R Cs 1 1 s ,
E (
1
s ) R R 1 2
1
R1R Cs 2
1 s
R R 1 2
where R1 C
α=
R 2
(0<α<1)
R1 R 2
This shows that ωm is the geometric mean of the two corner frequencies.
1 sin
m
The value of α for maximum phase lead Φm is given by α=
1 sin
m
0
Selecting maximum phase lead Φm = 40
0
We get α=
1 sin 40 0.2174
0
1 sin 40
R2 = R1
0.2174
1590
441.68 ,
1 1 0.2174
1
ω2= =28930 rad/s, f2 = 4604Hz.
m= 1 2
13482rad / s 2145 Hz
PROCEDURE:
1. The circuit is connected after selecting suitable values of resistors and capacitor as per design
for the given maximum phase lead.
2. Give the input to the network from a signal generator having variable frequency sine wave.
Adjust the magnitude of the input to a suitable value and note down this value.
TABULAR COLUMN:
(Note: The phase angle should be noted down with a positive sign for lead network)
Bode Diagram
0
Magnitude (dB)
-5
-10
-15
60
Phase (deg)
30
0
2 3 4 5 6
10 10 10 10 10
Frequency (rad/sec)
C = 0.1e-6;
sys = tf(num,den)
bode(sys);
grid;
RESULT:
INFERENCE:
LEAD compensator has been designed to obtain a specific value of maximum phase lead. It is observed through
experiment that that output leads the input and the same phase lead has been obtained. Transfer function of
the system is determined experimentally.
INTRODUCTION:
Compensation is the modification of dynamic and/or static characteristics of the control system by
adding some subsystem with the original system. The cascade combination of original system together
with the subsystem is called compensated system. This subsystem which is added for improved static
and/or dynamic characteristics is called compensator. Modification of system dynamics (equations by
which system is described) is required in order to obtain performance specification, like in case of Lead
and Lag compensator which are often designed to satisfy phase and gain margins, these bode plot of
the system itself changes after the insertion of compensator, bode plot indirectly represents system
dynamics.
The Lag Compensator is an electrical network which produces a sinusoidal output having the phase lag
when a sinusoidal input is applied.
Necessity of compensators
1. In order to obtain the desired performance of the system, we use compensating networks.
2. To compensate a unstable system to make it stable.
3. A compensating network is used to minimize overshoot.
4. Compensating networks increase the steady state accuracy of the system
AIM:
i) To design a lag network for a maximum phase angle of 54° choosing a suitable value of time constant and
obtain the frequency response.
ii) To determine the transfer function of the lag compensating network.
iii) Verify the result using simulation package.
APPARATUS REQUIRED:
1. CRO
2. Signal generator
3. decade resistance and capacitance box
4. MATLAB simulation Package
CIRCUIT DIAGRAM:
Signal
Generator
E1 E2
THEORY:
Effect of Phase Lag Compensation
1. Gain crossover frequency increases.
2. Bandwidth decreases.
3. Phase margin will be increase.
4. Response will be slower before due to decreasing bandwidth, the rise time and the settling
time become larger.
5. Improves the steady state performance in the expense of slower settling time.
Advantages of Phase Lag Compensation
1. Phase lag network allows low frequencies and high frequencies are attenuated.
2. Due to the presence of phase lag compensation the steady state accuracy increases.
Disadvantages of Phase Lag Compensation
Due to the presence of phase lag compensation the speed of the system decreases.
DESIGN:
G(s)=
E (s)
2
1R Cs 2 1 s
where R2 C
E (s)
1 1 R1 R Cs 2
1 s
β= R R
1 2
(β>1)
R 2
This shows that ωm is the geometric mean of the two corner frequencies.
1 sin
m
The value of β for maximum phase lag Φm is given by β=
1 sin
m
We get β=
1 sin 54 10
0
1 sin 54
R2 = 1k and
C
R1 = R2 (β-1) = 10x103(10-1) = 9 kΩ.
1
Higher corner frequency ω2= =10000rad/s, f2=1591Hz,
So ωm= 1 2
=3162 radians/sec, fm=503Hz.
Transfer function
1+𝑠𝜏 1+0.0001𝑠
= =
𝑎+𝑠𝜏𝛽 1+0.001𝑠
PROCEDURE:
1. The circuit is connected after selecting suitable values of resistors and capacitor as per design
for the given maximum phase lag.
2. Give the input to the network from a signal generator having variable frequency sine wave.
Adjust the magnitude of the input to a suitable value and note down this value
TABULAR COLUMN:
1 20 126.66 5 0 0
2 40 266.32 5 0 0
(Note: The phase angle should be noted down with a negative sign for lag network)
Bode Diagram
0
-5
Magnitude (dB)
-10
-15
-20
0
Phase (deg)
-30
-60
1 2 3 4 5 6
10 10 10 10 10 10
Frequency (rad/sec)
1. Calculate lower corner frequency (ω1) which is the frequency at which magnitude is 3 dB less
than the highest value.
2. Calculate higher corner frequency (ω2) which is the frequency at which magnitude is 3 dB more
than the lowest frequency.
3. Starting magnitude=0; so no need to calculate K.
1
1 s
2
4. Transfer function G(s)= K
1
1 s
1
CALCULATIONS:
From graph w1 = 930 rad/sec , w2 = 7100 rad/sec
1+0.00014𝑠
=
1+0.00107𝑠
C = 0.1e-6;
sys = tf(num,den)
bode(sys);
grid;
RESULT:
Transfer function
1 + 0.0001𝑠 1 + 0.0001𝑠 1 + 0.00014𝑠
1 + 0.001𝑠 1 + 0.001𝑠 1 + 0.00107𝑠
INFERENCE:
LAG compensator has been designed to obtain a specific value of maximum phase lag. Same phase lag has
been obtained experimentally and it is observed that output lags the input. Transfer function of the system is
determined experimentally.
EXPERIMENT: 06
INTRODUCTION:
A lead–lag compensator is a component in a control system that improves an undesirable frequency
response in a feedback and control system. It is a fundamental building block in classical control theory.
Lead–lag compensators influence disciplines as varied as robotics, satellite control, automobile
diagnostics, LCD displays and laser frequency stabilization. They are an important building block in
analog control systems, and can also be used in digital control.
AIM: To design a lag-lead network for a maximum phase lag of 55° and a maximum phase lead of 55°
choosing a suitable value of time constant and obtain the frequency response.
APPARATUS REQUIRED:
1. CRO
2. Signal generator
3. decade resistance and capacitance box
CIRCUIT DIAGRAM:
Signal
Generator
E1 E2
THEORY:
Lag-Lead compensator is an electrical network which produces phase lag at one frequency region and
DESIGN:
G(s)=
E ( s)
2 RC s 1 R C s 1 1 1 2 2
1 s 1
1 s 2
E ( s)
1 R C s 1 R C s 1 R C 2s
1 1 2 2 1 1
1 s 1 s 2
1
1 j 1
1 j 2
G(s) =
2
1 j 1
1 j
1
The frequency at which G (jω) = 0 is given by 0
1 2
1 1
The pole located at s together with zero located at s constitutes the lag element
1 1
1
while the zero located at s 1 2
together with pole located at s constitute the
2 2
1 1 1
lead element. Corner frequencies are respectively: , , ,
1 1 2 2
1 sin
m 2
Let, Φm=250, =2.5, 2.5
1 sin 1
m
Let f1=50Hz, ω2=314.16 radians /second and f2=2.5 x 50 Hz = 125 Hz, ω2=785.4 radians /second
1
2
, 1
=1.273 x 10-3
1
1
0
= 1241 rad/sec
1 2
Frequency at which phase difference is zero. F0 = 197 Hz.
R1 = 9 k R2 = 1 K
C1 = 0.14 µF C2 = 0.5 µF
EXPECTED GRAPH:
PROCEDURE:
TABULAR COLUMN:
Input voltage E1= 5 Volts
Sl No Frequency f =2 f Magnitude Phase angle Gain in dB
(Hz) (radians/sec) Of output Of E2- Φ 20log(E2/E1)
E2(volts) (degrees)
1. 1 23 144.5 3.5 -20 -3.09
.
2. 35 219.9 2.5 -25 -6.02
RESULT:
Transfer function
RESULT:
Experimental Value Designed Value
INFERENCE:
LAG-LEAD compensator has been designed to obtain a specific value of maximum phase lag and phase lead.
It is observed that output lags the input for lower frequency range and becomes in phase with the input at a
particular frequency and then leads.
Study the effect of P, PI, PD and PID controller on the step response
of a feedback Control System using simulation
INTRODUCTION:
PID controllers are commercially successfully and widely used as controllers in industries.
For example, in a typical paper mill 90% of the controllers used are PID controller. The PID
controller uses a proportional mode, an Integral mode and a Derivative mode. The first
letters of these modes make up the name PID controller. Depending upon the applications
one or more of these modes are used. For example, in a liquid control system, a zero steady
state error is required; thus a PI controller is used. Whereas in a temperature control system
where zero steady state error is not essential, a P controller can be used.
P Control. Output power is directly proportional to control error. The higher the proportion
coefficient, the less the output power at the same control error. Proportional control can be
recommended for fast-response systems with a large transmission coefficient.
PI Control. Output power equals to the sum of proportional and integration coefficients. The
higher the proportion coefficient, the less the output power at the same control error. The
higher the integration coefficient, the slower the accumulated integration coefficient. PI
control provides zero control error and is insensitive to interference of the measurement
channel. The PI control disadvantage is slow reaction to disturbances.
P-D Control - Output power equals to the sum of proportional and differential coefficients.
P-D control for lower order system apparently is not very useful as it cannot reduce the
steady state error to zero. But for higher order processes, it can be shown that the stability
of the closed loop system can be improved using P-D controller.
PID Control - Output power equals to the sum of three coefficients: proportional, integral
and differential. The higher the proportion coefficient, the less the output power at the same
control error. The higher the integration coefficient, the slower the accumulated integration
coefficient. The higher the differentiation coefficient, the greater the response of the
system to the disturbance. The PID controller is used in inertial systems with relatively low
AIM: To study the effect of P, PI, PD and PID controller on the step response of a feedback
control system using simulation.
BLOCK DIAGRAM:
THEORY:
Guidelines for selection of controller mode:
Proportional controller - It is simple regulating type; tuning is easy. But it normally introduces
steady state error. It is recommended for process transfer functions having a pole at origin,
or for transfer functions having a single dominating pole.
Integral Control: It does not exhibit steady state error, but is relatively slow responding. It is
particularly effective for: (i) very fast process, with high noise level (ii) process dominated by
dead time (iii) high order system with all time constants of the same magnitude.
Proportional plus Integral (P-I) Control: It does not cause offset associated with proportional
control. It also yields much faster response than integral action alone. It is widely used for
process industries for controlling variables like level, flow, pressure, etc., those do not have
large time constants.
Proportional plus Derivative (P-D) Control: It is effective for systems having large number of
time constants. It results in a more rapid response and less offset than is possible by pure
proportional control. But one must be careful while using derivative action in control of very
fast processes, or if the measurement is noisy (e.g. flow measurement).
In the early history of automatic process control the PID controller was implemented
as a mechanical device in steering system of Ships.
Electronic analog PID control loops are often found within more complex electronic
systems, for example, the head positioning of a disk drive, the power conditioning of
a power supply, or even the movement-detection circuit of a modern seismometer.
SIMULINK MODEL:
PROCEDURE:
1. System transfer function is taken as
GRAPH:
Controller: P [Kp=10]
Without Controller
Ki 500
Ki 10 Controller: PID [Kp+ + sKd = 500 + + 100s]
Controller: PI [Kp+ = 500 + ] s s
s s
TABULAR COLUMN:
Kp = 10 0.35 0
2. P
Kp = 500 0.037 50.57%
3. PI 0.03 48.45%
4. PD 0.037 24.6%
5. PID 0 0
INFERENCE
For the system given without controller, steady state error is 0.95 but overshoot = 0
By adding proportional controller, error can be reduced, but with the increase of Kp
as the error reduces overshoot increases. Error can not be eliminated by proportional
controller.
By adding PD controller, error does not change but, overshoot reduces.
PI controller reduces steady state error, has very less effect on overshoot.
When PID controller is connected, the system steady state as well as transient
response are improved.
EXPERIMENT - 8
INTRODUCTION:
The time response of control system consists of two parts. Transient response and
its independent variable. Analysis of response means to see the variation of output with
respect to time. The output of the system takes some finite time to reach to its final
value. Every system has a tendency to oppose the oscillatory behavior of the system
which is called damping. The damping is measured by a factor called damping ratio of
the system. If the damping is very high then there will not be any oscillations in the
output. The output is purely exponential. Such system is called an over damped system.
If the damping is less compared to over damped case then the system is called a critically
damped system. If the damping is very less then, the system is called under damped
system. With no damping system is un-damped.
Denominator of the transfer function equated to zero is the characteristic equation and
the order of the characteristic equation is called the order of a system. So if the order
of the characteristic equation is 2 then the system is called as second order system.
There are a number of factors that make second order systems important. They are
simple and exhibit oscillations and overshoot. Higher order systems are based on
second order systems. In case of mechanical second order systems, energy is stored in
the form of inertia whereas in case of electrical systems, energy can be stored in a
capacitor or inductor.
AIM:
THEORY:
C (s) 2
n , where ωn = natural frequency and ξ = damping ratio.
R( s) s2 2 n
2
n
For,
1< ξ < ∞ Over-damped System
0<ξ<1 Under-damped system
ξ=1 critically damped system
ξ=0 undamped system
1 2
tr where tan - 1
2
n 1
Theoretically, M p e 2
1
Steady-state error (ess ) is the difference between actual output and desired output
at the infinite range of time.
PROCEDURE:
1. In MATLAB m file type the program
PROGRAM:
% Determination of time response specifications of second order system
clc;
num=4; den=[1 1.2 4];
sys=tf(num, den);
step(sys)
2. Run the program, the time response curve will be obtained. Right click on the
figure characteristics settling time, peak response, Rise time. Bringing
the cursor at the point the values for rise time, peak time, overshoot and
settling time can be noted down.
PROGRAM:
% Time response specifications of second order system with addition of pole
clc;
num = 4; den=[1 1.2 4];
den_P = conv([1 1], den)
sys=tf(num, den_P);
step(sys)
PROGRAM:
% Time response specifications of second order system with addition of zero
clc;
num = 4;
den=[1 1.2 4];
num_Z = conv([1 1], num)
sys=tf(num_Z, den);
step(sys)
PROGRAM:
% Time response specifications of second order system with changed pole
location
clc;
num=4; den=[1 1.2 4];
den_P1 = conv([1 2], den)
sys=tf(num, den_P1)
step(sys)
7. Each time tabulate the data from the figure by right click characteristics.
Tabular Column:
CALCULATION:
Transfer function of the given system is
4
s 2 1.2s 4
C (s) 2
Comparing with standard second order system, n
R( s) s2 2 n
2
n
ξ = 0.3 and ωn = 2
1 2
tan - 1 = 1.266 rad
tr = 0.983 sec
2
n 1
Mp e 2 100 % = 40.8%
1
RESULT: For different ωn and ζ transfer function is calculated and time response
specifications are determined.
INFERENCE:
An additional zero in the left half of s plane makes the system faster and more
oscillatory. This can be seen from the simulations. As the zero moves along the
negative real axis toward the origin, the time to the first peak of the step response
decreases monotonically while the percent overshoot increases monotonically.
Also, it takes longer for the system to settle to the final value of the response.
The addition of left half pole tends to slow down the system response. As the pole
moves towards the origin in s plane, the rise time increases and the maximum
overshoot decreases.
EXPERIMENT - 9
INTRODUCTION:
It is required to keep the position of the load constant. Such a system where output
position is kept constant is called as position control system. The output position is sensed
and feedback to the potentiometer used as an error detector. For any change in the output
position θL, the potentiometer generates an error signal proportional to the difference θr
and θL. The θr is the reference position corresponding to the ideal output position. The
error signal is given to the amplifier and the output of the amplifier is given to the armature
of a d.c. servo motor. The d.c. servomotor maintains the output shaft position constant.
The entire scheme is called d.c. position control system.
AIM : TO Simulate the D.C position control system using Matlab and to observe the step
response when the system is (a) under damped (b) critically damped / over damped
THEORY:
Applications:
i. Electromechanical Actuators
ii. Air Craft Control Systems
iii. Process controllers
iv. Robotics
A dc-motor position control system is shown in figure. The following parameters and
variables are defined: e-error voltage , er=the reference input, θL=the load position, KA-
the amplifier gain=200 , ea-the motor input voltage, eb-the back emf, ia- the motor current
, Tm- the motor torque , Jm-the motor inertia=0.03, Bm- the motor viscous-friction co-
efficient=10, JL- the load inertia=0.05 , Kt-motor torque constant=21, Kb-the back emf
constant=15.5V/rpm, Ke-the error-detector gain=E/2π, where E=2π,Ra-the motor
resistance=1.15, θe= θr-θL
Now represent the transfer functions of different blocks, using Simulink and
interconnect them as shown in the following figure.
PROCEDURE:
1) Click on step input block and set final value as 1 and other parameters to 0.
2) Click on amplifier gain and set gain to 200. Similarly set gains of flexible shaft and error
detector to 1 each.
3) Click on transfer function block and set numerator value as 18.62 and
4) denominator to [0.08 293].
5) Save the interconnected block diagram.
6) Click on simulation tool bar and go to simulation parameters and set start time to 0.0 and
stop time to 0.5sec.
7) Click on the start simulation button and observe the waveforms of error, speed and
output θL. The output θL corresponds to critically/over damped second order step
response.
8) Change the denominator value of transfer function to [0.08 20] and repeat step 6 after
setting stop time to 1 sec. The output θL in this case corresponds to under damped
second order step response.
GRAPH:
INFERENCE:
DC position control system is simulated using MATLAB Simulink and response of the system
with variation of damping is observed. With reduced damping oscillation increases.
INTRODUCTION:
Steady-state error is defined as the difference between the input (command) and the
output of a system in the limit as time goes to infinity (i.e. when the response has reached
steady state). The steady-state error will depend on the type of input (step, ramp, etc.) as
well as the system type (0, 1, 2 …). In addition to these, imperfections in the system
components, such as static friction, backlash, and amplifier drift, as well as aging or
deterioration, will also cause errors at steady state.
Steady-state error analysis is only useful for stable systems. The stability of the system must
be checked before performing a steady-state error analysis. Many of the techniques that
we present will give an answer even if the error does not reach a finite steady-state value.
AIM: To verify the effect of input waveform, loop gain and system type on steady state
error.
THEORY:
The deviation of the output of control system from desired response during steady state is
known as steady state error. It is represented as ess. We can find steady state error using
the final value theorem as follows.
e ss
Lim e t
t
Lim sE s
s 0
1 1
Unit ramp, t Velocity ∞ 0
k v k v
Acceleration
1
Unit parabolic, t2 1 ∞ ∞
k a
k a
PROCEDURE:
1. System transfer function is taken as
( s 1)(s 2)
G( s) for type 0 system
( s 3)(s 5)
( s 1)(s 2)
G( s) for type 1 system
s( s 3)(s 5)
( s 1)(s 2)
G(s) for type 2 system
s 2 ( s 3)(s 5)
2. Connect the system using SIMULINK as shown in Simulink model. The block
ramp input can be obtained from source, sum and gain block from math
operations, transfer function and integrator block from continuous, mux from
commonly used block and scope from sink. A ramp input and integrator
together can produce parabolic input.
WAVEFORM:
With step input With ramp input
Type 0 Type 0
Type 1 Type 1
Type 2 Type 2
With parabolic input
Type 0 Type 1
Type 2
Step 0.882 0 0
INFERENCE:
Steady state error values are obtained for type 0, type 1 and type 2 system with step, ramp
and parabolic input. The result satisfies the theoretical values as shown in table 11.1.
INTRODUCTION:
A lead compensator can increase the stability or speed of response of a system. a lag compensator
can reduce (but not eliminate) the steady-state error. Depending on the effect desired, one or
more lead and lag compensators may be used in various combinations. A system which has one
pole and one dominating zero (the zero which is closer to the origin than all over zeros is known
as dominating zero.) is known as lead network. If we want to add a dominating zero for
compensation in control system then we have to select lead compensation network.
AIM: To perform a trade off study for lead compensation using MATLAB
THEORY:
Steps to design Lead Compensator using Bode Plot:
Step 1: Obtain the Bode plot of the given system.
Step 2: Find the phase lead angle to be provided by the compensator.
Is the required Phase Margin of compensated system and ε = 5 to 12°, this is to compensate
the shift in gain cross over frequency.
1 sin m 1
Step 3: Calculate and mag = -20log .
1 sin m
Step 4: Obtain the frequency (ωc) corresponding to this magnitude from the bode plot of
uncompensated system.
1 1
Step 5: But this new gain cross over frequency c . Thus T
T c
1 1
Step 6: Corner frequencies of the lead compensator are 1
and 2
T T
1 s 1
Step 7: Transfer function of the compensator is G c s
s 2
Step 8: Obtain the PM and GM of the system with compensator. If they do not meet the
required specifications, continue from step 2 with updated ε.
1. In order to obtain the desired performance of the system, we use compensating networks.
Compensating networks are applied to the system in the form of feed forward path gain
adjustment.
2. Compensate an unstable system to make it stable.
3. A compensating network is used to minimize overshoot.
4. These compensating networks increase the steady state accuracy of the system. An important
point to be noted here is that the increase in the steady state accuracy brings instability to the
system.
5. Compensating networks also introduces poles and zeros in the system thereby causes changes
in the transfer function of the system. Due to this, performance specifications of the system
change.
20
PROBLEM: Design a lead compensator for the system given by phase margin is
s ( s 0.5s )
BLOCK DIAGRAM:
Gc(s)
PROCEDURE:
1. Write the program in Matlab m file
Program1:
% Input the value of phi required for compensation
req_phi = 55;
Program 2:
req_phi = 55;
phi1=req_phi - pm + 6;
wc=9.49;
T=1/(sqrt(alpha)*wc);
w1=-1/T;
w2=-1/(alpha*T);
Note: Take the reading of phase margin and gain margin to check if the compensator
satisfies the given requirement. Vary the tolerance for further adjustment.
Bode Diagram
Gm = Inf dB (at Inf rad/sec) , Pm = 18 deg (at 6.17 rad/sec)
50
40
30
20
Magnitude (dB)
10
-10
System: sys
Frequency (rad/sec): 9.5
-20 Magnitude (dB): -7.24
-30
-40
-50
-90
Phase (deg)
-135
-180
-1 0 1 2
10 10 10 10
Frequency (rad/sec)
Bode Diagram
Gm = Inf dB (at Inf rad/sec) , Pm = 54.9 deg (at 9.49 rad/sec)
50
sys_overall
Magnitude (dB)
-50
-100
-90
Phase (deg)
-135
-180
-1 0 1 2 3
10 10 10 10 10
Frequency (rad/sec)
INFERENCE:
Phase margin of the uncompensated system = 18 ° and gain margin = infinity.
Required phase margin = 55°
502893( s 4.123)
Transfer function of the compensator Gc( s)
(s 21.87)
502893( s 4.123) 20
Transfer function of the overall system = Goverall ( s)
( s 21.87) s( s 0.5s)
Phase margin of the overall system after compensation = 54.9° and gain margin = infinity.
By adding the Lead compensator the phase margin of the system is improved to meet the required
specification.
EXPERIMENT: 10
ROOT LOCUS PLOT USING MATLAB
INTRODUCTION:
The root locus is a plot of the roots of the characteristic equation of the closed-loop system
as a function of one system parameter varies, such as the gain of the open-loop transfer
function. It is a method that determines how the poles move around the S-plane as we
change one control parameter. This plot was introduced by Evans in 1948 and has been
developed and used extensively in control engineering.
AIM:
1. To obtain the root locus plot of a given open loop transfer function using
MATLAB and determine
a. The range of open loop gain K for which the closed loop system is
stable.
b. The breakaway and/or break-in points and the corresponding gain.
2. Verify the effect of open loop poles and zeros on root locus
3. Verify the effect of open loop gain on root locus.
THEORY:
The Root Locus Method tells us the position of the poles in the s-plane for each value of a
control parameter. By using root locus method the designer can predict the effects of location
of closed loop poles by varying the gain value or adding open loop poles and/or open loop
zeroes. The closed loop poles are the roots of the characteristic equation. This method brings
in to focus the complete dynamic response of the system. By using root locus method the
designer can predict the effects of location of closed loop poles by varying the gain value or
adding open loop poles and/or open loop zeroes. The closed loop poles are the roots of the
characteristic equation. This method gives the information about absolute and relative
stability of the system.
PROCEDURE:
1. Assign the numerator (num) and denominator (den) polynomial.
2. Define the system by using the command
sys=tf(num,den)
3. The root locus plot is obtained by using the command
rlocus(sys)
4. To determine the range of K for which the system is stable: In order to find
the range of K for stability of closed loop system the point of intersection of
the root locus with the imaginary axis (if it intersects) is required. It can be
found using the function
rlocfind(sys).
5. Use the above function from the command window, select the point of
intersection with the imaginary axis on the root locus by moving the crosshair
cursor and click on it. The gain and the corresponding locations of all poles will
be displayed. Now decide the range of K for stability of closed loop system.
Example:
k
1.
s(s 4)(s 5)
Program 1:
clc
clear all
k=1
sys = zpk([], [0 -4 -5],1);
rlocus(sys)
title('Root Locus Plot for k/s(s+4)(s+5)');
rlocfind(sys)
Root Locus
15
System: sys
Gain: 274
Pole: 0.409 + 5.27i
10
Damping: -0.0774
Overshoot (%): 128
Frequency (rad/sec): 5.29
5
Imaginary Axis
0
System: sys
Gain: 13.1
Pole: -1.47 - 3.62e-008i
-5
Damping: 1
Overshoot (%): 0
Frequency (rad/sec): 1.47
-10
-15
-14 -12 -10 -8 -6 -4 -2 0 2 4
Real Axis
𝐾
2. G(S)H(S)=𝑆(𝑆 2 +4𝑆+13)
Program 2.
clc
clear all
k=1
num = k;
den =[1 4 13 0];
sys = tf(num/den)
rlocus(sys);
𝐾
title('Root Locus Plot for G(S)H(S)=𝑆(𝑆 2+4𝑆+13)');
rlocfind(sys)
RESULT:
𝐾 (ifany)
1. 1.472 ±j4.472 ------ 1.42 ±j4.47 ----
𝑆(𝑆+4)(𝑆+5)
𝐾
2. ---- ±j3.6 -33.69° --- ±j3.6 -33.69°
𝑆(𝑆 2 + 4𝑆 + 13)
3. Addition of pole:
𝑲
G(S)H(S) =
𝑺(𝑺+𝟒)(𝑺+𝟓)(𝒔+𝟏)
Program 3
clc
clear all
k=1
sys = zpk([], [0 -1 -4 -5], 1);
rlocus(sys)
title('Root Locus Plot with addition of pole’);
rlocus(sys)
4.Addition of zero:
𝑲(𝒔+𝟏)
G(S)H(S) =𝑺(𝑺+𝟒)(𝑺+𝟓)
PROGRAM
clc
clear all
k=1
sys = zpk([-1], [0 -4 -5]
rlocus(sys)
title('Root Locus Plot with addition of zero');
rlocus(sys)
Theoretical verification.
𝑲
1. G(S)H(S) =
𝑺(𝑺+𝟒)(𝑺+𝟓)
2*0 1
θ 1= *180 =60o
3 0
2 *1 1 2*2 1
θ 2= *180 =180o θ3= *180 =300o
3 0 3 0
(𝑠𝑢𝑚 𝑜𝑓 𝑟𝑒𝑎𝑙 𝑝𝑎𝑟𝑡𝑠 𝑜𝑓 𝑝𝑜𝑙𝑒)−(𝑠𝑢𝑚 𝑜𝑓 𝑟𝑒𝑎𝑙 𝑝𝑎𝑟𝑡𝑠 𝑜𝑓 𝑧𝑒𝑟𝑜𝑠)
6. centroid=
𝑝−𝑧
0−4−5−0
= = -3.
3−0
7. Break away point:
The characteristic equation is
1+ G(S)H(S)=0
𝐾
1+ 𝑆(𝑆+4)(𝑆+5)
K=-(s3+ s2+20s+k)
dk/ds=0
3s+2s+20=0
s=-1.472,-4.52
s=-1.472 is valid break away point.
S1 180 - k/9
s0 k
𝑲
2. G(S)H(S)= 𝟐
𝑺(𝑺 +𝟒𝑺+𝟏𝟑)
2*0 1
θ 1= *180 =60o
3 0
2 *1 1 2*2 1
θ 2= *180 =180o θ 3= *180 =300o
3 0 3 0
(𝑠𝑢𝑚 𝑜𝑓 𝑟𝑒𝑎𝑙 𝑝𝑎𝑟𝑡𝑠 𝑜𝑓 𝑝𝑜𝑙𝑒)−(𝑠𝑢𝑚 𝑜𝑓 𝑟𝑒𝑎𝑙 𝑝𝑎𝑟𝑡𝑠 𝑜𝑓 𝑧𝑒𝑟𝑜𝑠)
6. centrod=
𝑝−𝑧
s3+ 4s2+13s+k=0
Routh array
s3 1 13
s2 4 K
S1 (52-k)/4 0
s0 k
EXPERIMENT: 11
INTRODUCTION: Magnitude and phase shift of the frequency response can be plotted directly
on linear scale but to see the effect in a wide range of frequency the linear plot becomes
cumbersome. Thus The plot of the magnitude and phase shift with respect to frequency is
plotted in logarithmic scale . This is known as Bode Plot after Hendrick Wade Bode.
Hendrick W. Bode was an American engineer, researcher, inventor,
author and scientist, of Dutch ancestry. As a pioneer of
modern control theory and electronic telecommunications he
revolutionized both the content and methodology of his chosen fields
of research. He made important contributions to the design, guidance
and control of anti-aircraft systems during World War II. He helped
develop the automatic artillery weapons that defended London from
the V-1 flying bombs during world war II. He also made important
contributions to control systems theory and mathematical tools for
the analysis of stability of linear systems, inventing Bode plots, gain Hendrick Wade Bode
margin and phase margin. [ 1905-1982]
AIM: To obtain comparative study of Bode, Nyquist and Root locus with respect to stability
using Matlab.
THEORY:
Condition for determination of Stability from Bode Plot: For a system to be stable, the gain
crossover frequency is to be less than phase cross over frequency or both the gain margin
and phase margin should be positive.
Nyquist Stability Criterion: Nyquist stability criterion states that For the closed loop system
to be stable the number of encirclements of the nyquist plot about the critical point (-1+j0)
must be equal to the poles of the open loop transfer function in the right half of the 's' plane.
Nyquist plot has advantages over the root locus and Routh-Horwitz because it easily handles
time delays.
Gain Margin: Greater will the gain margin greater will be the stability of the system. It refers
to the amount of gain, which can be increased or decreased without making the system
unstable. It is usually expressed in dB.
Phase Margin: Greater will the phase margin greater will be the stability of the system. It
refers to the phase which can be increased or decreased without making the system unstable.
It is usually expressed in terms of phase angle.
Gain Cross Over Frequency: It refers to the frequency at which magnitude curve cuts the zero
dB axis in the bode plot.
Phase Cross Over Frequency: It refers to the frequency at which phase curve cuts the
negative times the 180o axis in this plot.
Corner Frequency: The frequency at which the two asymptotes cuts or meet each other is
known as break frequency or corner frequency.
The margins are undefined when there is no crossover at 0dB and if there is no crossover at
180°.
K
G( s)
(s 1)(s 2)(s 5)
PROCEDURE:
A. Root locus: Plot the root locus for the same system using MATLAB and obtain the
range of K for stability.
Program:
sys=zpk([], [-1 -2 -5], 1);
rlocus(sys);
rlocfind(sys);
B. Bode Plot:
The Bode diagrams can be obtained using MATLAB either by executing MATLAB
statements one by one through command window or by creating a matlab file (.m file)
of program statements and executing it as follows.
1. In the matlab editor enter the program to obtain the Bode diagrams for the given
transfer functions. A sample program is given below:
%bode diagrams
%plot of system1
sys=zpk([], [-1 -2 -5], 60);
% One value less than Kmargin (in place of 60) for stable system,
% one value greater than Kmargin for unstable system
% one value = Kmargin for marginaaly stable system
margin(sys);
iv. Find the phase angle by clicking on the phase curve at = g. The
phase margin is calculated as PM=180 + phase angle at = g.
4. The GM, PM and the cross over frequencies can also be obtained using the
following function.
[GM PM WCF WCG]= margin (sys);
Where GM=Gain margin is in abs units ( 20logGM is the GM in db) PM=Phase margin
WCF=Phase cross over frequency WCG=Gain cross over frequency.
These parameters will be displayed in the command window.
Note: Verify the result with a gain less than Kmargin, equal to Kmargin and greater
than Kmargin and conclude.
C. Nyquist plot
Dept. of Electrical and Electronics Engineering - 96
Global Academy of Technology
Control Systems Lab manual 18EEL66
1. The Nyquist plot can be obtained using MATLAB either by executing MATLAB
statements one by one through command window or by creating a matlab file (.m
file) of program statements and executing it as follows:
Program:
2. Check the number of open loop poles (P) in the right half of s plane.
Find out the number of encirclements (N) about the point -1+j0 in
anticlockwise direction.
3. I f N=P the system is stable. Otherwise the system is not stable.
4. Find out the point where magnitude = 1, corresponding frequency is gain
cross over frequency.
5. Find out the point when phase angle = -1800, corresponding frequency is
phase crossover frequency.
6. Alternatively right click on the plot -> Characteristics -> All stability Margin, Now if
cursor is brought on the dots the GM, corresponding frequency i.e phase
crossover frequency, PM and the corresponding frequency i.e gain crossover
frequency can be noted down.
Note: Verify the result with a gain less than Kmargin, equal to Kmargin and greater than
Kmargin and conclude
TABULAR COLUMN
1. Define system
• A system is a collection of elements or components that are organized for a common purpose.
2. Define control system
• A control system is a system of devices or set of devices, that manages, commands, directs or
regulates the behavior of other devices or systems to achieve desired results.
3. Define transfer function
• For a linear, time invariant, single input single output system Transfer function is defined as the ratio
of laplace transform of output to laplace transform of input with all initial conditions set to zero.
4. Define poles and zeros.
• Values of s which makes transfer function equal to infinity is called poles. Poles can be obtained by
the solving the equation denominator of transfer function to zero.
• Values of s which makes transfer function equal to zero is called zeros. zeros can be obtained by the
solving the equation numerator of transfer function to zero.
5. Define order of the system.
• Servo motors are intended and designed to use in motion control applications which require high
accuracy positioning, quick reversing and exceptional performance. A servo motor is a linear or
rotary actuator that provides fast precision position control for closed-loop position control
applications. Unlike large industrial motors, a servo motor is not used for continuous energy
conversion.
Control Systems Lab manual 18EEL66
9. What is ac servo motor?
• AC servo motors are basically two-phase squirrel cage induction motors and are used for low power
applications. Nowadays, three phase squirrel cage induction motors have been modified such that
they can be used in high power servo systems.
• The main difference between a standard split-phase induction motor and AC motor is that the
squirrel cage rotor of a servo motor has made with thinner conducting bars, so that the motor
resistance is higher. Based on the construction there are two distinct types of AC servo motors, they
are synchronous type AC servo motor and induction type AC servo motor.
10. What is dc servo motor?
• A DC servo motor consists of a small DC motor, feedback potentiometer, gearbox, motor drive
electronic circuit and electronic feedback control loop. It is more or less similar to the normal DC
motor. The rotor consists of brush and shaft. A commutator and a rotor metal supporting frame are
attached to the outside of the shaft and the armature winding is coiled in the rotor metal supporting
frame. A brush is built with an armature coil that supplies the current to the commutator. At the
back of the shaft, a detector is built into the rotor in order to detect the rotation speed.
• With this construction, it is simple to design a controller using simple circuitry because the torque is
proportional to the amount of current flow through the armature.
• DC servo motor provides very accurate and also fast respond to start or stop command signals due
to the low armature inductive reactance.
11. List the differences between AC servo motor and DC servo motor
AC servo motor DC Servo motor
It delivers high power output Delivers low output of about 0.5 W to 100 W
about 5 to 20%
The life of DC servo motor depends on the life on The life of AC servo motor depends on bearing life
brush life
• Applications of AC servomotor
• automation, robotics, CNC machinery, and other applications involving a high level of
precision
• DC servo motors are used in similar equipments and computerized numerically controlled
machines.
19. What is damping ratio? Specify the range of damping ratio for un-damped, under damped,
critically damped and over damped system.
• Damping ratio basically indicates the amount of damping present in the overall system
denoted by zeta, where damping is a counter force. If the ratio is zero, that indicates there is no
damping present and as such the system will oscillate forever (oscillations are the natural
frequency response). A system with ratio between 0 & 1 has some amount of damping and
accordingly the oscillations will die down slowly or quickly as damping ratio changes from 0 to 1.
• A system with damping ratio 1 has no overshoots and the quickest rise time.
• Systems with damping ratio more than 1 are overdamped i.e. such systems have more than
desired damping (friction) that causes the system response to be sluggish.
• ζ= 1 is called critically damped, ζ >1 is called overdamped
• ζ= 0 is undamped and the oscillations never decay to zero
20. What are the time domain specifications?
• Time response specifications are delay time, rise time, maximum overshoot, peak time and
settling time.
21. Define delay time, rise time, settling time, peak time, overshoot and steady state error.
• Settling Time Ts is the time required for the response to remain within a certain percent of its
final value, typically 2% to 5%.
• Delay time (td) is the time required to reach at 50% of its final value by a time response signal
during its first cycle of oscillation.
• Rise time (tr) is the time required to reach at final value by a under damped time response
signal during its first cycle of oscillation. If the signal is over damped, then rise time is counted as
the time required by the response to rise from 10% to 90% of its final value.
• Peak time (tp) is simply the time required by response to reach its first peak i.e. the peak of
first cycle of oscillation, or first overshoot.
Mr
M is greater than or equal to 2 is defined as bandwidth.
1
• Corner Frequency: The frequency at which the magnitude is equal to 2 times of its
maximum value is known as cut off frequency or corner frequency. In other words, frequency at
which the magnitude is 3 dB less than its maximum value is called the cut-off frequency.
• Gain Margin: Greater will the gain margin greater will be the stability of the system. It refers
to the amount of gain, which can be increased or decreased without making the system unstable.
It is usually expressed in dB.
• Phase Margin: Greater will the phase margin greater will be the stability of the system. It refers
to the phase which can be increased or decreased without making the system unstable. It is
usually expressed in phase.
• Gain Cross Over Frequency: It refers to the frequency at which magnitude curve cuts the zero
dB axis in the bode plot.
Control Systems Lab manual 18EEL66
• Phase Cross Over Frequency: It refers to the frequency at which phase curve cuts the negative
times the 180o axis in this plot.
• Corner Frequency: The frequency at which the two asymptotes cuts or meet each other is
known as break frequency or corner frequency.
25. Difference between time domain and frequency domain analysis
Time Domain analysis Frequency Domain Analysis
A time domain graph shows how a signal The frequency domain graph shows how
changes over time much of the signal lies within each given
frequency band over a range of
frequencies.
Time domain is the domain for analysis of Frequency domain is the domain for
mathematical functions or signals with analysis of mathematical functions or
respect to time. signals with respect to frequency.
Transfer function calculation is not possible Transfer function of a system can be
experimentally using time domain analysis. experimentally calculated using frequency
domain analysis