Lab Manual 6EE7 Cs
Lab Manual 6EE7 Cs
Lab Manual 6EE7 Cs
TECHNOLOGY, JAIPUR
LAB-MANUAL
VI SEM EE
6EE07
CONTROL SYSTEM LAB
(b) Ramp
(c) Impulse
5. To design 2nd order electrical network and study its trarient response for step input and
following cases. (a) Under damped system (b) Over damped System. (c) Critically damped
system
6. To Study the frequency response of following compensating Networks, plot the graph and
final out corner frequencies. (a) Lag Network (b) Lead Network (c) Log-lead Network.
7. To draw characteristics of a.c servomotor
8. To perform experiment on Potentiometer error detector.
9. Check for the stability of a given closed loop system.
10. Plot bode plot for a 2nd order system and find GM and PM.
ROTOR # 1
1.
2.
3.
4.
5.
Find UJT static emitter characteristics and study the variation in peak point and valley
point
6.
Study and test firing circuits for SCR-R, RC and UJT firing circuits.
ROTOR # 2
7.
Study and test 3-phase diode bridge rectifier with R and RL loads. Study the effect of
filters.
8.
Study and obtain waveforms of single-phase half wave controlled rectifier with and
without filters. Study the variation of output voltage with respect to firing angle.
9.
10.
Study and obtain waveforms of single-phase full controlled bridge converter with R and
RL loads. Study and show rectification and inversion operations with and without
freewheeling diode.
11.
Control the speed of a dc motor using single-phase half controlled bridge rectifier and
full controlled bridge rectifier. Plot armature voltage versus speed characteristics.
LAB ETHICS
DOs
1. Enter the lab on time and leave at proper time.
2. Keep the bags outside in the racks.
3. Utilize lab hours in the corresponding experiment.
4. Shut down the computers before leaving the lab.
5. Maintain the decorum of the lab.
Donts
1. Dont bring any external material in the lab.
2. Dont make noise in the lab.
3. Dont bring the mobile in the lab. If extremely necessary then keep ringers off.
4. Dont enter in server room without permission of lab incharge.
5. Dont litter in the lab.
6. Dont delete or make any modification in system files.
7. Dont carry any lab equipments outside the lab
We need your full support and cooperation for smooth functioning of the lab.
INSTRUCTIONS
WHEN WORKING AT COMPUTER LAB
BEFORE ENTERING IN THE LAB
1. All the students are supposed to prepare the theory regarding the next program.
2. Students are supposed to bring the practical file and the lab copy.
3. Previous program should be written in the practical file.
4. Algorithm & Program of the current program should be written in the lab copy.
5. Any student not following these instructions will be denied entry in the lab and
Sessional Marks will be affected.
INSTRUCTIONS
WHEN WORKING AT CONTROL SYSTEM LAB
BEFORE ENTERING IN THE LAB
6. All the students are supposed to prepare the theory regarding the present
Experiment.
7. Students are supposed to bring the practical file and the lab copy.
8. Previous experiment should be written in the practical file.
9. Object, Apparatus Table & Brief Theory of the current practical should be
written in the lab copy.
10. Any student not following these instructions will be denied entry in the lab and
Sessional Marks will be affected.
EXPERIMENT # 1
OBJECT: INTRODUCTION TO MATLAB COMPUTING SOFTWARE.
1.1 What is MATLAB?
MATLAB is a software package for high-performance numerical computation and
visualization. It provides an interactive environment with hundreds of built-in functions for technical
computation, graphics, and animation. Best of all, it also provides easy extensibility with its own
high-level programming language. The name MATLAB stands for MATrix LABoratory. MATLAB is
the product of the Math Works Inc.
1.2 Functions of MATLAB
MAT Labs built-in functions provide excellent tools for linear algebra computations, data analysis,
signal processing, optimization, numerical solutions of ordinary differential equations (ODEs),
quadrature, and many other types of scientific computations. Most of these functions use state-of-the
art algorithms. There are numerous functions for 2-D and 3-D graphics as well as for animation. Also,
for those who cannot do without their FORTRAN or C codes, MATLAB even provides an external
interface to run those programs from within MATLAB. The user, however, is not limited to the builtin functions; he can write his own functions in the MATLAB language. Once written, these functions
behave just like the built-in functions. MAT Labs language is very easy to learn and to use.
There are also several optional Toolboxes available from the developers of MATLAB. These
Toolboxes are collections of functions written for special applications such as Symbolic
Computation, Image Processing, Statistics, Control System Design, Neural Networks, etc.
1.3 General Ideas of MATLAB
Here we discuss some basic feature and commands. To begin, let us look at the general structure of
the MATLAB environment.
1.3.1
MATLAB Windows
On almost all systems, MATLAB works through three basic windows, which are discussed below
1. Command window: This is main window. It is characterized by the MATLAB command prompt
>>. When you launch the application program, MATLAB puts you in this window. All commands,
including those for running user-written programs, are typed in this window at the MATLAB prompt.
In MATLAB 6, this window is a part of the MATLAB window that contains four other smaller
windows.
6
Launch Pad: This sub window lists all MATLAB related applications and toolboxes that are
installed on your machine. You can launch any of the listed applications by double clicking on
them.
Workspace: This sub window lists all variables that you have generated so far and shows
their type and size. You can do various things with these variables, such as plotting, by
clicking on a variable and these variables, such as plotting, by clicking on a variable and then
using the right button on the mouse to select your option.
Command History: All commands typed on the MATLAB prompt in the command window
get recorded, even across multiple sessions (you worked on Monday, then on Thursday, and
then on next Wednesday, and so on), in this window. You can select a command from this
window with the mouse and execute it in the command window by double clicking on it. You
can also select a set of commands from this window and create an M-file with the right click
of the mouse (and selecting the appropriate option from the menu).
Current Directory: This is where all files from the current directory are listed. You can do
file navigation here. You also have several options of what you can do with a file once you
select it (with a mouse click). To see the options, click the right button of the mouse after
selecting a file. You can run M-files, rename them, delete them, etc.
2. Graphics window: The output of all graphics commands typed in the command window are
flushed to the graphics or Figure window, a separate gray window with (default) white background
color. The user can create as many figure windows as the system memory will allow.
3. Edit window: This is where you write, edit, create, and save your own programs in files called
M-files. You can use any text editor to carry out these tasks. On most systems, MATLAB provides
its own built-in editor. However, you can use your own editor by typing standard file-editing
command that you normally use on your system. From within MATLAB, the command is typed at
the MATLAB prompt following the special character !. The exclamation character prompts
MATLAB to return the control temporarily to the local operating system, which executes the
command the following the ! character. After the editing is completed, the control is returned to
MATLAB. For example, on UNIX systems, typing! VI myprogram.m at the MATLAB prompt (and
hitting the return key at the end) invokes the VI editor on the file myprogram.m Typing! Emacs
myprogram.m invokes the emacs editor.
7
1.3.2 Input-Output: - MATLAB supports interactive computation, taking the input from the screen,
and flushing the output to the screen. In addition, it can read input files and write output files. The
following features hold for all forms of input-output:
Data type: The fundamental data type in MATLAB is the array. It encompasses several
distinct data objects-integers, doubles (real numbers), matrices character strings, structures
and cells. In most cases, however, you never have to worry about the data type or the data
object declarations. For example, there is no need to declare variables as real or complex.
When a real number is entered as the value of a variable, MATLAB automatically sets the
variable to be real (double).
Case sensitivity: MATLAB is case-sensitive; that is, it differentiates betn the lowercase and
uppercase letters. Thus a and A are different variables. Most MATLAB commands and builtin functions calls are typed in lowercase letters. You can turn case sensitivity on and off with
the casesen command. However we do not recommend it.
Output display: The output of every\ command is displayed on the screen unless MATLAB
is directed otherwise. A semicolon at the end of a command suppresses the screen output,
except for graphic and on line help commands.
1.3.3 File types: - MATLAB has three types of files for storing information:
M-files are standard ASCII text files, with a .m extension to the filename. There are two types
of these files: script files and function files. Most programs you write in MATLAB are saves
as M-files. All built-in functions in MATLAB are M-files, most of which reside on your
8
computer in precompiled format. Some built-in functions are provided with source code in
readable M-files so that can be copied and modified.
Mat-files are binary data files, with a .mat extension to the filename. Mat-files are created by
MATLAB when you save data with the save command. The data is written in a special format
that only MATLAB can read. Mat-files can be loaded into MATLAB with the load command.
Mex-files are MATLAB-callable FORTRAN and C programs, with a .mex extension to the
filename. Use of these files requires some experience with MATLAB and a lot of patience.
whos
what
clear
clear x y z
clear all
clc
home
same as clc
clf
cd
dir
ls
path
editpath
copyfile
copies a file
mkdir
creates a directory
clock
date
more
ver
gives the license and the version information about MATLAB installed on
your computer
bench
1.5.4 Termination
control-c
quit
quits MATLAB
exit
same as quit
RESULT:- We have successfully studied the basics of the MATLAB programing Software.
VIVA-VOCA QUESTIONS/ANSWERS:
Q. 1 What is MATLAB?
Ans. MATLAB is a numerical computing environment and programming language.
MATLAB is
a software package for high-performance numerical computation and visualization. MATLAB means
MATrix+LABoratory.
Q.2 Who invented the MATLAB software?
Ans. MATLAB was invented in the late 1970s by Cleve Moler, then the chairman of the computer
science department at the University of New Mexico. He designed it to give it to his students access
to LINPACK and EISPACK without having to learn Fortran. It soon spreads to other universities and
found a strong audience within the applied mathematics community. The MATLAB was rewritten in
C and founded the MathWorks Incorporated in 1984 to continue its development.
10
EXPERIMENT # 2
11
OBJECT: - To design 1st order R-C circuits and observes its response with the following inputs
and traces the curve.
(a) Step input
(b) Ramp input
(c) Impulse input
APPARATUS REQUIRED: PC with MATLAB 7.0.4 software.
THEORY:First order system:
For a linear variant system, when the power of s in denominator is one then it is called First order
system.
Block diagram representation is shown in fig 2.1:
(t ) I (t ) R
1
i (t )d (t )
C
( s ) R[ I ( s)]
1
I (s)
Cs
12
( s ) I ( s ) R
Cs
(1)
V0 (s)
V (s)
V
V
1
Cs
R+
(2)
1
Cs
1
1 RCs
1
1 RCs
Or transfer function =
Put RC=T
C ( s)
1
R ( s) Ts 1
We conclude that this is a Ist order transfer function. So R-C circuit is Ist order system. With transfer
function.
C ( s)
1
R ( s) Ts 1
Case (1):
For unit step input response
1
R( s)
s
1
1
C (s)
*
Ts 1 s
1
T
C (s)
s Ts 1
1
1
s s 1
T
Taking inverse Laplace
C (t ) 1 e
e(t ) 1 1 e
e (t ) e
for t 0
T
t
lim e(t ) 0
t
Graph of C(t) V/S t is shown in fig 3.3.
ss
13
T
T
S
Ts 1
C (t ) t T T e
for t0
e(t ) r (t ) c(t )
t
e(t ) T 1 e T
e(t ) T
ess tlim
e (t ) e
t0
ess= e() = 0.
B) For unit ramp input
C (t ) t T T e
e(t ) T 1 e
es s= T.
t0
t0
17
EXPERIMENT # 3
OBJECT: - DESIGN THE SECOND ORDER ELECTRICAL NETWORK AND STUDY IT
TRANSIENT RESPONSE FOR STEP INPUT UNDER FOLLOWING CASES.
CASE 1: Undamped
(=0)
CASE 2: Underdamped (=0.5)
CASE 3: Critical damped (=1)
CASE 4: Over damped
(=1.5)
APPARATUS REQUIRED: - PC with MATLAB 7.0.4 software.
THEORY:6.1 TIME RESPONSE OF SECOND ORDER SYSTEM: - The block diagram of second order
control system is in fig.6.1.
Fig.3.1
. (3.1)
From equation (3.1)
. (3.2)
18
. (3.3)
Break the equation (6.3) by partial fraction, put
. (3.4)
Multiply equation (3.4) by s and put s=0
. (3.5)
Laplace inverse pf equation (6.5)
. (3.6)
Put
19
. (3.7)
The error signal for the system
. (3.8)
The steady state value of c (t)
20
From the expression (3.7), it is clear that the time constant is 1/ and the response having
damped oscillation with overshoot and undershoot. Such response is known as underdamped
response.
. (3.9)
Thus at n the system will oscillate (with =0). The damped frequency always less then the undamped
frequency (n) because of factor . If the system having certain value of then it is not possible to
measure undapmed natural frequency experimentally. The observed frequency is the damped
frequency (d) which is equal to n 1- 2. The response is shown in Fig.
or
. (3.10)
. (3.11)
Multiply both the sides by (s +n) 2 and put s = - n
. (3.12)
. (3.13)
Inverse Laplace of equation (3.13)
. (3.14)
. (3.15)
Inverse Laplace of
22
. (3.18)
Break the equation (3.18) by partial fraction
Equation (3.18) can be written as
.(3.19)
23
Put
. (3.20)
Inverse Laplace of equation (3.20)
24
. (3.21)
From equation (3.21) we get two times constant
From equation (3.21) we observe that when is greater than one there are two exponential term, the
first term has a time constant T1 which is smaller than the time constant of other exponential term
(having time constant T2 ), in other word we can say that the first exponential term decaying much
faster than other exponential term. So, for time response we can neglect it, then
.(3.22)
and time constant
.(3.23)
For different values of the curves of C (t) is shown in Fig.3.5
From the curves it is clear that the overdamped systems are sluggish.
n=6
num= [wn*wn]
den= [1 2*zeta*wn wn*wn]
g=tf (num, den)
step (g)
CASE 2:- Underdamped
zeta=0.5
n=6
num= [wn*wn]
den= [1 2*zeta*wn wn*wn]
g=tf (num, den)
step (g)
CASE 3:- Critical damped
zeta=1
n=6
num= [wn*wn]
den= [1 2*zeta*wn wn*wn]
g=tf (num, den)
step (g)
CASE 4:- Overdamped
zeta=1.5
n=6
num= [wn*wn]
den= [1 2*zeta*wn wn*wn]
g=tf (num, den)
step (g)
Case 1:- Output for the undamped case.
zeta = 0
n =6
num =36
den =1 0 36
36
Transfer function= -------s2 + 36
26
Step Response
2
1.8
1.6
A m p litu d e
1.4
1.2
1
0.8
0.6
0.4
0.2
0
10
12
Time (sec)
Fig.3.6:- Undamped
Case 2:- Output for the under damped case.
zeta=0.5000
n = 6
num =36
den =1 6 36
36
Transfer function = -------------s2 + 6 s + 36
Step Response
1.4
1.2
A m p litu d e
0.8
0.6
0.4
0.2
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
Time (sec)
Fig.3.7:-Underdamped
Case 3:-Output for the critically damped case.
Zeta =1
n =6
num =36
27
den =1
12
36
36
Transfer function= ----------------s2 + 12 s + 36
Step Response
1
0.9
0.8
A m p litu d e
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
Time (sec)
1
0.9
0.8
A m p litu d e
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0.5
1.5
2.5
Time (sec)
28
Peak overshoot in step input response reduces as damping is increased from 0.2 to 0.6.
Response peak in frequency response reduces as damping is increased from 0.2 to 0.6.
Non of the above statement is true
Statement (i) is true but statement (ii) is false.
Statement (ii) is true but statement (i) is false.
Both the statement are true.
(d)
Q.8 Undamped natural frequency n and response frequency r of a unity feedback system
with open-loop transfer function.
; < 1/2
are related as
(a) n =r (b) n >r
(c)n <r
(d) None
29
Ans:
(b)
Q.10 For the system c(s)/r(s) =16/ s2+8s+16. The nature of the time response will be.
(a)Overdamped
(b)Underdamped
(c) Critically damped
(d)None of the
Ans:
(c)
EXPERIMENT # 4
OBJECT: - TO DRAW THE CHARACTERISTIC OF A.C. SERVO MOTOR.
APPRATUS REQUIRED: - P.C. With MAT LAB 7.0.4 Software, 2-Phase A.C. Servo Motor etc
THEORY: - Servo-Motor is a closed loop control system in which a small input power controls a
larger output power in a strictly proportionate manner. The control variable is mechanical position,
displacement, velocity or torque etc. its used in control system where very high acceleration is
desired.
A.C. Servo-Motor is 2-phase A.C. Induction motor. Its called 2-phase A.C. servo motor as there are
two windings, one winding called fixed or reference winding is supplied with a fixed voltage &
frequency from a constant voltage source. Second winding is called control winding with variable
supply voltage of same frequency.
The stator windings for the fixed and control winding are placed 90
apart in space. So the torque is produced most efficiently on the shaft when the phase winding axes
are in space & volt is in the two phases are in qardrature.
The instantaneous control volt is given by:Ec(t) = Ec(t) Sin wt , for
Ec(t)>0(1)
Ec(t)= |Ec(t)| Sin(180+wt) , for Ec(t)<0 (2)
The change in (-Ec) reverse the direction of rotation of motor. Due to constant reference voltage the
torque (T) & angular speed (0) are also function of the control volts Ec (t).
The Servomotor provides a large torque at zero speed. This torque is necessary for rapid acceleration.
30
(9)
T = J. d20/dt2 + B.d0/dt(10)
J = moment of inertia of motor & load referred to motor shaft.
B = viscous friction coefficient of the motor & load referred to motor shaft
31
Fig:-4.1
32
Fig 4.3
Where, X/R=Small (linear)
[V4>V3>V2>V1]
X/R=Large (Non-linear)
RESULT: - We have successfully studied to the A.C. servomotor and in that control .voltage
Result in the development of the motor torque (Tm). The servomotor provides a large
Torque at zero speed. This torque is necessary for rapid acceleration.
THE ANALYSIS OF RESULT:
IN THE 2-phase A.C. Servo-motor is also an Induction Motor having drag-up type rotor
construction, The control voltage Ec(t) is applied to the to the control winding & a fixed voltage
having a phase difference of 90 w.r.t. control winding voltage is applied to the reference winding. The
control voltage result in the development of the motor Torque (Tm).
The Servo-motor provides a large at zero speed. This torque is necessary for rapid Acceleration.
SERVO-MOTORS are widely used in Radars, Electro-mechanical, Actuators,
Computers, Machine tools, tracking and Quittance system, process controllers
and ROBOTS.
VIVA-VOCA QUESTIONS/ANSWERS:
33
34
EXPERIMENT # 5
OBJECT: - TO PERFORM EXPERIMENT ON POTENTIOMETER ERROR DETECTOR.
APPARATUS REQUIRED: 1) Experimental Board of Potentiometer.
2) Power supply (230 V).
3) Connecting leads.
THEORY: INTRODUCTION
A potentiometer is an electromechanical transducer which converts the mechanical energy
(Displacement) into electrical energy (Voltage). It is also called error detecting device.
Potentiometric transducers are relatively inexpensive and easy to apply.
Precision potentiometers are simple rotary devices for obtaining shaft position information. The most
straight forward application is the conversion of mechanical position to a voltage. Basically a
precision potentiometer consists of a resistive element with a movable arm or slider in contact with
the element. As the arm (slider) rotates, resistance varies between the end of the resistive element and
the slider, indicating shaft position. The resistive element can be made of wire, conductive film or
cermet element.
Potentiometers used for servomechanisms are generally above 7/8 to 3 and half inches in diameter.
The early models were mostly of wire wound type. Current technology provides other choices such
as stability, longer life and lower sensitivity to environment. Potentiometer can be excited with
alternating and direct current. Single turn potentiometers have a rotation i.e. usually limited to 350
degrees.
POTENTIOMETER CHARACTERISTICS.
A linear potentiometer produces a resistance change i.e. linearly related to the shaft position. A
position of rotation will produces 50% of maximum resistance and a position of rotations will
produce 75% of maximum resistance. Linearity is specified as the deviation (in percentage of total
resistance) of the actual resistance. This is called normal or independent linearity.
35
Where
yi
yt
36
Where
i
t
The applied voltage is normally D.C. and polarity of output voltage describes the relative
position of the shafts. In case of A.C., the phase difference will find the relative position of the two
shafts.
Resolution for Potentiometer: -It is defined as the ratio of change in the output
voltage in
one step to the supply voltage.
OPERATING INSTRUCTIONS
1. Connect the mains cord to the supply line.
2. You can study the potentiometer as transducer by measuring voltage across the output points
of the potentiometer and other the readings in the table.
3. Plot the graph for above readings.
4. Now you can study the system as null detector by keeping left hand potentiometer at zero
degrees and treating it as reference, rotate the right hand potentiometer to read the error
directly on the DPM. Note that these are ten turn potentiometers.
5. You may also plot input position verses output position when meter is indicating null.
SOME MORE IMPORTANT INSTRUCTION.
1. D.P.M. can be connected across the O/P and GND point of either potentiometer or it can be
connected across the variable terminals of both pots for studying error detector function.
37
OBESRVATION:SR. NO.
1.
2.
Angular position
Output Voltage
38
3.
4.
5.
6.
7.
8.
9.
10.
RESULT: - After performing this experiment we have obtain the relation between angular
Position and output voltage which can be shown in graphical manner.
ANALYSIS OF RESULT: - we are obtained the graph between angular position and output voltage
and this graph shows that the output voltage varies linearly with the angular position of the shaft.
APPLICATION OF POTENTIOMETER:A. SELF-BALANCING POTENTIOMETER
The merits of dc potentiometer for accurate measurement of dc voltage have prompted the
development of self-balancing (automatic) potentiometer for the measurement of dc voltages. They
are widely used for indication and record of voltage developed by various transducers. The automatic
balancing is provided be an ac servomotor. The servomotor drives the tapping point of the
potentiometer, the writing mechanism, and the pointer for indication. The motor comes to rest when
the unbalance voltage becomes zero. The block diagram is shown in figure No. 9.6. The output
displacement and voltage feedback are related by the calibration constant of the potentiometer wire.
The amplifier is tuned to 50 Hz and has high gain. The unbalance dc voltages are converted into
proportional ac voltage of 50 Hz be means of a synchronous chopper. The forward path consist of the
servomotor and the inertial load of the writing pen. By virtue of the negative feedback, the bandwidth
of the system is increased, and the system is used to record faithfully voltages of frequency range dc
to 5 Hz. The system is considered to function as voltage to displacement converter. The output shaft
can be coupled to a shaft angle encoder it, in addition to analog indication, digital display is desired.
39
Q.4
Ans.
Q.5
Ans.
Q.6
Ans.
ANGULAR POSITION
IN DEGREE
OUTPUT VOLTAGE
IN MV.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
10
30
60
90
120
150
180
210
240
270
300
330
360
10
29
59
90
120
150
181
211
240
269
299
329
359
ANGULAR POSITION
IN DEGREE
RESISTANCE IN
K-Ohm
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
00
30
60
90
120
150
180
210
240
270
300
.001
.080
.162
.247
.331
.413
.496
.577
.659
.742
.825
41
12.
13.
330
360
.907
.990
I/P ANGULAR
POSITION IN
DEGREE
20
60
90
140
180
240
240
320
340
350
O/P ANGULAR
POSITION IN
DEGREE
30
80
120
170
220
265
320
350
360
360
ERROR
IN DEGREE
10
19
28
30
39
27
40
30
20
9
EXPERIMENT # 6
OBJECT: - CHECK FOR THE STABILITY OF A GIVEN CLOSED LOOP SYSTEM.
APPARATUS REQUIRED: - P. C. with MAT LAB 7.0.4 software.
THEORY:CONCEPT OF STABILITY
The concept of stability is very important to analyze and design the system. A system is said to
be stable if its response cannot be made to increase indefinitely by the application of a bounded
input excitation. If the output approaches towards infinite value for sufficiently large time, the
system is said to be unstable.
A linear time invariant (LTI) system is stable if
(1) The system is excited by a bounded input, the output is bounded. (BIBO stability criteria)
(2) In the absence of the input, the output tends to zero(the equilibrium state of the system).
This is known as asymptotic stable.
NECESSARY BUT NOT SUFFICIENT CONDITIONS FOR STABILITY:Consider a system with characteristic equation
a0 sm+a1sm-1+.+am=0
(a)All the coefficients of the equation should have same sign,
(b)there should be no missing term.
If above two conditions are not satisfied the system will be unstable. But if all the coefficients
have same sign and there is no missing term we have no guarantee that the system will be
stable.
There are two methods for check the stability of a given closed loop.
<1>ROUTH-HURWITZ CRITERION:Consider the following characteristics polynomial
a0s n+a1sn-1+.+an=0. ................... (1)
42
Where the coefficients a0, a1.an are all of the same sign and none is zero.
STEP<1>Arrange all the coefficients of above equ. (1) in two rows.
ROW 1
a0
a2
a4
ROW 2
a1
a3
a5
STEP<2> From these two rows form a third row
ROW1
a0
a2
a4
ROW2
a1
a3
a5
ROW3
b1
b3
b5
Where,
b1=-1a1 aoa2a1a3
b3=-1a1 a0a4a1a5
STEP<3>From second and third row, form a fourth row
ROW1
a0
a2
a4
ROW2
a1
a3
a5
ROW3
b1
b3
b5
ROW4
c1
c3
c5
c1=-1b1 a1a3b1b3
c3=-1b1 a1a5b1b5
STEP<4> Continue this procedure of forming a new rows
<1.1>STATEMENT OF ROUTH-HURWITZ CRITERION
Routh-Hurwitz criterion states that the system is stable if and only if all the element in first
column has the same algebraic sign. If all elements are not of the same sign then the number
Of sign changes of the elements in first column equals the number of roots of the characteristic
equation in right half of the s-plane (or equals to the number of roots with positive real
Parts)
INPUT:- A closed loop control systems has the characteristic equation given by
S3+4.5S2+3.5S+1.5=0
Investigate the stability using Routh-Hurwitz criterion.
OUTPUT:S3
1
3.5
2
S
4.5
1.5
S1
3.5
S0
1.5
No. of sign changes in first column=0
No. of roots in right half of S- plane=0
Hence, System is stable.
<2.>ROOT LOCUS METHOD:Root locus method is a graphical method in which roots of the characteristic equation are plotted
in S-plane for the different values of parameter. The locus of the roots of the characteristic
Equation when gain is varied from zero to infinity is called Root-locus.
<2.1>RULES FOR CONSTRUCTION OF ROOT LOCII:Following are the rules to sketch the root locus plot.
Rule1: The root locus is symmetrical about the real axis.
Rule2: The root locii starts from an open loop pole with K=0 e.g. For the system having
G(s) H(s)=K(s+3)/(s+2)..(2)
Find the starting point of root locii.
Solution:-According to the rule the root locii starts from s=-2
43
Rule3: The root locii will terminate either on an open loop zeros or on infinity with K= e.g.,
Find the ending point of the root locii given in equation (2). According to the rule the root locii
will terminate at s=-3
Rule4:-If
N=No. of separate locii
P=No. of finite poles
Z=No. Of finite zeros then
Number of root locii will be equal to the no. Of poles if number of poles are more than number
of zeros i.e.P>Z
N=P if p>z
If z>p, then number of root locii will be equal to the number of zeros.
If p=z, then No. of root locii=poles =zeros.
E.g.find the number of separate root locii for the system given by the equation (2).
Solution:
p=1
Z=1
N=1
Rule5: ROOT LOCII ON THE REAL AXIS.
Any point on the real axis is a part of the root locus if and only if the number of poles and zeros
to its right is odd.
Rule6:ASYMPTOTES
The branches of the root locus tend to infinity along a set of straight line called a asymptotes.
These asymptotes making an angle with real axis and is given by
Q= (2K+1)1800/p-z where K=0,1,2,
The total number of asymptotes=p-z
e.g.
If G(s) H(s)=K/s(s2+6s+10)(3)
p=3
z=0
No. of asymptotes=p-z=3-0=3
K=0 Q1= (2*0+1)*1800/3=600
K=1
Q2= (2*1+1)*1800/3=1800
K=2
Q3= (2*2+1)*1800/3=3000
Rule7: CENTROID OF ASYMPTOTES
The point of intersection of asymptotes with real axis is called centroid of asymptotes (SA)
and is given by
SA= (sum of poles-sum of zeros)/p-z
E.g. Find the centroid of asymptotes of the system given by equation(3)
Solution: There are three poles at s1=0,s2=-3+j1,s3=-3-j1
No. of zeros=0
So, centroid
SA=0-3+J1-3-J1-0/3=-2
Rule8: ANGLE OF DEPARTURE & ANGLE OF ARRIVAL OF THE ROOT LOCII
The angle of departure of the root locus from a complex pole is given by
SD=1800-(sum of the angles of vectors drawn to this pole from other poles) + (sum of angles of
vectors drawn to this pole from the zeros).
The angle of arrival at a complex zero is given by
SA==1800-(sum of angles of vectors drawn to this zero from other zeros) + (sum of angles of
vectors drawn to this zero from poles).
E.g. For
G(s) H(s) =k/s(s+6)(s2+4s+13)
Determine the angle of departure from complex poles.
SD=1800-(P1+P2+P2)
44
=1800-(1230+370+900) =-700
So, angle of departure at (-2+j3) =-700
So, angle of departure at (-2-j3)=700
Rule9: BREAKAWAY POINT ON REAL AXIS
If the root locus lies between two adjacent open loop poles on the real axis then there will be at
least one breakaway point , because the roots move towards each other as K is increased and
meet at a point. At this point K is maximum. If we increase the value of K between two poles the
root locus breaks in two parts.
Similarly if root locus lies between two adjacent zeros on real axis, then there will be at least one
break in point. If the root locus lies between an open loop pole and zero, then there will be no
breakaway or break in point or may be both occur.
The breakaway or break in points can be determined from the roots of
Dk/ds=0
e.g. if
G(s) H(s)=K/s(s2+6s+10)determine the breakaway point.
1+G(s) H(s)=1+k/s(s2+6s+10)
S (s2+6s+10)+k=0
Or, k=-s3-6s2-10s
dk/ds=-32-12s-10=0
s1=-1.1835 & s2=-2.815 are the breakaway points.
Rule10: The intersection of root locus branches with jw-axis can be determined through Routh
Hurwitz criterion.
e.g. If G(s)H(s)=k/s(s2+6s+10).Find the intersection of root
locii with the imaginary axix.
Solution: The characteristics equation s3+6s+10s+k=0
S3
1
10
S2
6
k
1
S
60-k/6
S0
k
Hence, we get a zero row if k=60
The auxiliary equation A(s)=6s2+k
6s2+k=0
6s2+60=0
S=+-j3.16
The root locus branches cross the imaginary axis at s=+-j3.16 for k=60.
INPUT: Plot the root locii for the closed loop control system and check the stalibity with
K
G(s)= ---------------------s(s+1)(s2+4s+5)
, H(s)=1
break
end
OUTPUT:
poles =
-10.0156
0.0078 + 9.5998i
0.0078 - 9.5998i
k = 71
1.I have studied the concept of stability is very important to analyse and design the system .
2.Any system is said to be stable if its response cannot be made to increase indefinitely by the
application of a bounded input excited .
3.A linear time invariant (LTI) sytem is stable if
(a) The system is excited by a bounded input , the output is bounded(BIBO).
(b)In the absence of the input , output tends towards zero.
Some standard criteria must be study to check the stability.i.e.
(a)all the coefficient of the equestion should have same sign.
(b) There should be no missing term.
Otherwise system will be unstable.
VIVA-VOCA QUESTIONS/ANSWERS:
Q1 Why concept of stability is nessary?
Ans: Concept of stability is nessary to analyse and design the system.
Q.2 what is the meaning of stability?
Ans: A system is said to be stable if its response cannot be made to increase indefinite by the
application of a bounded input excitation.
Q.3 What is linear time invariant(LTI)system.
Ans: A system which output is not vary with
time called LTI.
Q.4 What is the rule for stability of linear invareiant system?
Ans (a) The system is excitedby a bounded input, the output is bounded.
(b) In the absence of the input, the output tends towards zero.
Q.5 What is an asymptotic stability?
Ans: In the absence of input, the output tends towards zero(the equilibrium state of a system),
is known as asymptotic stable.
Q.6 What are the necessary conditions for stability?
Ans: There are some necessary conditions for stability(a) All the coefficients of the should have same sign.
(b) There should be no missing term.
Q.7 What is the statement of Routh-Hurwitz criterion?
Ans: Routh-Hurwitz criterion states that the system is stable if and only if all the elements in
the first column have the same algebraic sign. If all elements are not of the same sign then the
number of sign changes of the elements in first column equals the characteristics equation in
the right half of the s-plan.
Q.8 What is root locus?
Ans: The locus of the roots of the characteristics equation when gain varied from zero to
infinity is called root locus.
Q.9 What is the meaning of equilibrium state ?
Ans: When input is absence then output tends to zero is called equilibrium state
48
EXPERIMENT # 7
OBJECT: -PLOT BODE PLOT OF GIVEN TRANSFER FUNCTION IN SECOND ORDER
CONTROL SYSTEM.
Gs=25s2+4s+25
APPARATUS REQUIRED : - P C with MAT lab 7.0.4 software
THEORY: Simplified and desirable form but actual solution of the problem must still be accomplished .
solution of the problem may be either in the form of analysis or design for which following
alternatives are possible.
Transformed differential equation keeping all initial condition zero may be obtained after
manipulating the transfer function and the transient solution can also be obtained by inserting
the proper boundary conditions and then using laplace transform inversion.
The transfer function can be graphically represented on the s-plane to evaluate the roots of the
characteristic equation, open loop frequency response or closed loop frequency response
through graphical manipulations.
The transfer function G(j) just by replacing s by j and which can then be used to obtain
open loop and closed loop frequency response curve directly.
Advantages :
One of the advantage of frequency response method of control system analysis and design is
that in this method sinusoidal signal is used as standard test signal which is easily available.
Secondly, with the use of experimental data unknown transfer function of complicated
components and system can be determine
In frequency response method many powerful graphical tools are available such as Bode plot,
magnitude versus phase plot and nyquist plot etc.
given transfer function just by replacing s by j hence sinusodial transfer function and
frequency transfer function is nothing but transfer function in which s is replaced by j.
For example: s transfer function is given as
Ts=KS+1ss+5s2+2s+1
BODE PLOT
Bode plot developed by H.W BODE can be simplified because they can be approximated
as a sequence of straight line. Straight line approximation simplifies the evaluation of
magnitude and phase frequency response. Bode plot consist of two plots. One is logarithm
of magnitude versus frequency and other is phase angle versus frequency. Both the plots
are drawn on single semi log graph paper. This plot can be sketched from their asymptotic
properties as no detailed plotting is required. Bode plot are also known as asymptotic plots
and bode plot can be sketched by approximating the magnitude and phase width
asymptotic straight line which requires very less time. Some of the important parameters
of frequency response like gain margin, phase margin, gain crossover frequency and phase
crossover frequency can easily be determined from bode plot.
LIMITATION
Bode plot can be used to determined stability of only minimum phase system (minimum
phase system is the system which has all its poles and zeros on the left half of s-plane.
Bode plot in comparison to polar plot can be more easily and quickly constructed.
For example
Gs=25s2+4s+25
Change the transfer function in s+1 form
=1s225+4s25+1
Putting an equation s=j
=11-225+0.15
We determine the values of magnitude in db and phase angles in degrees at different
values of frequency
S.NO
FACTOR
CORNER
MAGNITUDE IN DB
FREQUENCY
1.
11-225+0.15
=5
Straight line of
constant slope
-20db/decade passing
through 0 db line at
50
=5
1.
0.1
-0.039
2.
0.2
-0.06
3.
0.3
-0.17
4.
-0.355
5.
-0.80
6.
-1.603
7.
-0.000872
8.
9.
4.57
10.
2.50
11.
17.06
51
Bode Diagram
20
10
Magnitude (dB)
0
-10
-20
-30
-40
-50
-60
0
Phase (deg)
-45
-90
-135
-180
-1
10
10
10
Frequency (rad/sec)
Fig 7.1
SOURCE CODE: num=[0 0 25];
den=[1 4 25];
g=tf(num,den)
bode(g)
grid on
[mag,phase,w]=bode(g)
OUTPUT: - The bode plot of the given transfer function is shown in figure.
52
10
Magnitude (dB)
20
0
-20
-40
Phase (deg)
-60
0
-45
-90
-135
-180
-1
10
10
10
10
Frequency (rad/sec)
Fig 7.2
ANALYSIS OF RESULT:
Successfully done the bode plot of given transfer function and successfully plot the
graph.
VIVA-VOCA QUESTIONS/ANSWERS:
Q.1 which software is used in MAT LAB?
Ans: MAT LAB 7.0.4 software are used.
Q.2 why bode plot are used?
Ans: Bode plot can be used to determined stability of only minimum phase system (minimum phase
system is the system which has all its poles and zeros on the left half of s-plane.
Q.3 which parameters are used in bode plot?
Ans: Some of the important parameters of frequency response like gain margin, phase margin, gain
crossover frequency and phase crossover frequency can easily be determined from bode plot.
Q.4 what is the unit of frequency?
Ans: Unit of frequency is rad/sec.
Q.5 what is bode plot?
Ans: Bode plot can be simplified because they can be approximated as a sequence of straight lineis
known as bode plot.
Q.6 How much plot consist of bode plote?
Ans: Bode plot consist of two plots. One is logarithm of magnitude versus frequency and other is
phase angle versus frequency. Both the plots are drawn on single semi log graph paper.
Q.7 In transfer function s is replaced in which form?
Ans: In transfer function s is replaced by j.
53
54