Control System
Control System
in
1
www.aktutor.in
QUANTUM SERIES
For
B.Tech Students of Third Year
of All Engineering Colleges Affiliated to
Dr. A.P.J. Abdul Kalam Technical University,
Uttar Pradesh, Lucknow
(Formerly Uttar Pradesh Technical University)
CONTROL SYSTEM
By
Varun Gupta
TM
Obtain transfer functions to predict the correct operation of open loop and closed
CO 1 loop control systems and identify the basic elements, structures and the K3
characteristics of feedback control systems.
Measure and evaluate the performance of basic control systems in time domain.
CO 2 K4
Design specification for different control action.
Analyze the stability of linear time-invariant systems in time domain using Routh-
CO 3 K4
Hurwitz criterion and root locus technique.
Determine the stability of linear time-invariant systems in frequency domain using
CO 4 K4
Nyquist criterion and Bode plot.
Design different type of compensators to achieve the desired performance of
CO 5 control System by root locus and Bode plot method. Develop and analyze the K5
intermediate states of the system using state space analysis.
KL- Bloom’s Knowledge Level (K1, K2, K3, K4, K5, K6)
K1 – Remember K2 – Understand K3 – Apply K4 – Analyze K5 – Evaluate K6 – Create
Detailed Syllabus:
Unit-I:
Control System Concepts: Elements of control systems, concept of open loop and closed loop systems,
Examples and application of open loop and closed loop systems, Mathematical Modelling of Physical Systems
(Electro Mechanical), Determination of transfer function by block diagram reduction techniques and signal
flow method using Mason’s gain formula, Basic Characteristics of negative feedback control systems.
Control System Components: Constructional and working concept of AC & DC servomotor, synchro’s,
stepper motor and tachometer.
Unit-II:
Time Response Analysis: Standard test signals, time response analysis of first and second order systems, time
response specifications of second order system for unit step input, location of roots of characteristics equation
and corresponding time response, steady state errors and error constants.
Basic modes of feedback control: Proportional, Derivative, Integral and PID controllers.
Unit-III:
Stability and Algebraic Criteria: Concept of stability and its necessary conditions, Routh-Hurwitz criteria
and its limitations.
Root Locus Technique: Salient features of root locus plot, Procedure for plotting root locus, root contours.
Unit-V
Introduction to Design: The design problems and preliminary considerations of lead, lag and lead-lag
compensation networks, design of closed loop systems using compensation techniques in time and frequency
domains.
State Space Technique: The concept of state & space, State-space model of physical system, conversion of
state-space to transfer function model and vice-versa, State transition matrix, Concept of controllability and
observability and their testing.
Text Book:
1. I. J. Nagrath & M. Gopal, “Control System Engineering”, 6th Ed. New Age International
Publishers, 2018. 4th Edition
2. M. Gopla, “Control System Principles and Design”, McGraw Hill 4th Edition
3. Ogata, "Modern Control Engineering, 5th Edition", Pearson Education, 2015
4. B.C. Kuo & Farid Golnaraghi, “Automatic Control Systems”, 10th Edition, McGraw Hill
5. D. Roy Choudhary, “Modern Control Engineering”, Prentice Hall of India.
1. Salivahanan, "Control Systems Engineering", Pearson Education, 2015.
6. Ambikapathy, “Control Systems”, Khanna Publishers
Reference Books:
1. (Schaums Outlines Series) Joseph J. Distefano III, Allen R. Stubberud, Ivan J. Williams, “Control
Systems”, 3rd Edition, McGraw Hill, Special Indian Edition, 2010.
2. Norman S. Mise, Control System Engineering, Wiley Publishing Co.
3. Ajit K Mandal, “Introduction to Control Engineering” New Age International.
4. R.T. Stefani, B.Shahian, C.J.Savant and G.H. Hostetter, “Design of Feedback Control Systems”
Oxford University Press.
5. Samarjit Ghosh, “Control Systems theory and Applications”, Pearson Education.
1 Control System
Concepts
• Open Loop and Closed Loop Systems and their Sensitivity Analysis
A. Concept Outline : Part-4 ........................................................... 1–30C
B. Long and Medium Answer Type Questions ............................. 1–31C
Control System Concepts
www.aktutor.in 1–2 C (EN-Sem-5)
PART-1
Concept of Control System,
Physical Systems and their Mathematical Modeling.
Questions-Answers
Que 1.1. Explain open loop and closed loop control system with
the help of suitable examples.
Answer
A. Open loop control system :
1. The open loop control system is also known as control system without
feedback or non-feedback control system.
2. In open loop systems the control action is independent of the desired
output.
3. In this system the output is not compared with the reference input.
4. The component of the open loop systems are controller and controlled
process.
Controlled
Input Controller Output
process
Fig. 1.1.1.
Examples :
1. Automatic washing machine
2. Immersion rod
3. A field control DC motor.
B. Closed loop control system :
1. Closed loop control system is also known as feedback control system.
Control System
www.aktutor.in 1–3 C (EN-Sem-5)
Fig. 1.1.2.
Examples :
1. Air conditioners
2. Autopilot aeroplane
3. Electric iron.
Que 1.2. Compare the open loop control system and closed loop
control system, also give few examples for each system.
AKTU 2013-14, Marks 05
Answer
Answer
Translation mechanical systems :
1. A translational spring-mass-damper system is shown in Fig. 1.3.1.
f(t) Applied force
x Displacement
K
M
d2 x dx
M =–D – Kx + f(t) ...(1.3.1)
dt2 dt
4. Rearranging eq. (1.3.1), we get
d2 x dx
M D Kx = f(t) ...(1.3.2)
dt 2 dt
Force-voltage analogy :
1. The analogy of eq. (1.3.2) with the voltage equation of an electrical
circuit can be established by considering electrical circuit shown in
Fig. 1.3.2.
R L
e(t) i C
Answer
Translational mechanical system : Refer Q. 1.3, Page 1–4C, Unit-1.
Equation of spring-mass-damper system
d2 x dx
M D Kx = f(t) ...(1.4.1)
dt 2 dt
Force-Current analogy :
1. The nodal equation as obtained for the circuit shown in Fig. 1.4.1 is
analogous to spring-mass-damper system described by the eq. (1.4.1).
iL iR iC
i(t) L e
R C
1
iR = e ...(1.4.4)
R
de
iC = C ...(1.4.5)
dt
4. Substituting eq. (1.4.3), (1.4.4) and (1.4.5) in eq. (1.4.2),
1 1 de
edt e C
L
= i(t) ...(1.4.6)
R dt
5. The voltage e is related to flux linkages associated with inductance L as
follows :
d
e= ...(1.4.7)
dt
d
6. Putting e = in eq. (1.4.6), we get
dt
1 1 d d2
C 2 = i(t)
L R dt dt
d 2 1 d 1
C = i(t) ...(1.4.8)
dt2 R dt L
7. Comparing eq. (1.4.1) and (1.4.8), we get
Table 1.4.1. Force-Current Analogy
Mechanical system Electrical system
Force, f Current, i
Velocity, v Voltage, e
Displacement, x Flux linkage,
Mass, M Capacitance, C
Damping coefficient, D Conductance, G
Stiffness, K (compliance, 1/K) Inductance, L
Control System
www.aktutor.in 1–7 C (EN-Sem-5)
PART-2
Construction and Working of AC and DC Servomotor,
Synchros, Stepper Motor and Tachometer.
Questions-Answers
Answer
Construction :
1. Fig. 1.5.1 shows the schematic diagram of a 2 AC servomotor.
2. The stator has two distributed windings which are displaced from each
other by 90 electrical degrees.
3. One winding is called the reference or fixed phase and other winding is
called control phase.
Working principle :
1. Reference phase is supplied from a constant voltage source Vr 0°.
The other winding i.e., control phase is supplied with a variable voltage
of the same frequency as the reference phase but its phase is displaced
by 90° (electrically).
Control System Concepts
www.aktutor.in 1–8 C (EN-Sem-5)
Reference winding
Vr
Applications :
1. In feedback control systems
2. In aircraft and spacecraft
3. In tracking and guidance systems
4. In robotics
5. In radar and machine tools
6. In process controller
7. In radio controlled cars and airplanes.
Answer
1. Consider the armature controlled DC motor (DC servo motor) and
assume that the demagnetizing effect of armature reaction is neglected,
field voltage is constant and magnetic circuit is linear. Armature
controlled DC servo motor is shown in Fig. 1.6.1.
Control System
www.aktutor.in 1–9 C (EN-Sem-5)
Ra La
B
ia J
V E T
if
Lf
Constant
Fig. 1.6.1.
2. Let Ra = Armature resistance
La = Armature inductance
ia = Armature current
E = Induced emf in armature
V = Applied armature voltage
= Angular displacement of the motor shaft
T = Torque developed by motor
J = Equivalent moment of inertia of motor shaft
and load referred to the motor
B = Equivalent viscous friction coefficient.
3. Apply KVL in armature circuit
La dia
V= Ra ia E ...(1.6.1)
dt
4. In the armature-controlled DC motor, the field current is held constant.
For a constant field current, the flux becomes constant, and the torque
becomes directly proportional to the armature current so that
T ia
T = Kia ...(1.6.2)
5. When armature is rotating, an emf is induced
K b d
E= ...(1.6.3)
dt
6. The armature current produces the torque which is applied to the inertial
mass and friction hence the force balance equation is
Jd 2 Bd
= T = KIa ...(1.6.4)
dt 2 dt
7. Taking the Laplace transform on both sides of eq. (1.6.4) and (1.6.1)
(sLa + Ra) Ia(s) + E(s) = V(s) [Initial condition is zero]
(Js2 + Bs) (s) = T(s) = KIa(s) ...(1.6.5)
Control System Concepts
www.aktutor.in 1–10 C (EN-Sem-5)
+ 1 Ia(s) K (s)
V(s)
– sLa + Ra s(Js + B)
E(s)
sK b
Fig. 1.6.2.
Answer
1. A schematic diagram of a field controlled DC motor (DC servo motor)
shown in Fig. 1.7.1.
Rf Ra
ia = Constant
ef if Lf J ea
T
B
Fig. 1.7.1.
T = K2 i f ...(1.7.2)
where K2 is a constant.
5. The equations for this system are
dif
Lf R f if = ef ...(1.7.3)
dt
Jd 2 Bd
and = T = K2 i f ...(1.7.4)
dt 2 dt
6. By taking the Laplace transform on both sides of eq. (1.7.3) and (1.7.4)
where all initial conditions are zero, we get
(Lfs + Rf) if(s) = Ef(s) ...(1.7.5)
(Js2 + Bs) (s) = K2if(s) ...(1.7.6)
7. From the above equations, the transfer function of this system is obtained
as
(s) K2
= ...(1.7.7)
E f (s ) s( Lf s Rf )( Js B)
8. Block diagram is shown in Fig. 1.7.2.
1 If(s) K2
Ef(s) (s)
L fs + R f s(Js + B)
Fig. 1.7.2.
Answer
A. Principle of synchros :
1. It is a rotary transducer that converts angular displacement into an AC
voltage or an AC voltage into an angular displacement.
2. A synchros system consists of
i. A control transmitter (CX) and
ii. A control transformer (CT).
3. The control transmitter consists of a stator and a rotor. The rotor is a
dumb-bell shaped magnetic structure. The supply is given to the rotor
by means of slip rings, which are actually mounted on the stator housing.
4. The secondaries are in the skewed slot all along the periphery of the
stator and are 120° apart because of their mechanical displacement.
5. The induced secondary voltage will depend upon the angle of the rotor
shaft. For reference the zero degree position of the shaft is defined
when the rotor is in alignment with the coil S2.
Control System Concepts
www.aktutor.in 1–12 C (EN-Sem-5)
S2
Stator
Slip rings
R1
AC supply
R2 S1
S3
R1 Synchro Control R1
transmitter (CX) transformer (CT)
S2 S2 Output voltage
AC source proportional to
sin ( r – L)
Stator Stator
R2 R2
S3 S3
S1 S1
r L
Fig. 1.8.2. Synchro error detector.
Answer
A. Stepper motor :
1. A stepper motor is a form of AC motor. The input given to this motor is
in the form of electric pulses.
2. For every input pulse, the motor shaft turns through a specified number
of degrees, called a step.
3. Shaft of stepper motors moves through one angular step for each
input pulse. The range of step size may vary from 0.72° to 90°.
B. Principle of operation :
1. Stators of stepper motors have salient poles on which concentrated
windings are wound. These windings may be appropriately connected
so as to result in two, three or four phase windings, on the stator.
2. The rotors of stepper motors carries no winding. The rotors are made
from either permanent-magnet or ferromagnetic material.
3. A stepper motor is fed through an external drive circuit. The function
of this drive circuit is to receive input voltage pulses and to deliver
appropriate currents to the stator windings of the stepper motor.
4. These currents set up air-gap field which moves through one angular
step for each input pulse. The rotor follows the axis of air-gap magnetic
field on account of the development of reluctance torque and/or
permanent-magnet torque.
5. Stepper motors vary widely in their designs and configurations, these
can be classified into three main types; variables-reluctance motors,
permanent-magnet motors and hybrid motors.
C. Operation of variable reluctance stepper motor :
Rotor
SA A
SB
SD D B
C
Stator
Vd SC
Answer
Tachometer :
1. Tachometer is electromechanical device that convert mechanical energy
into electrical energy.
2. The device works as a generator with the output voltage proportional to
the magnitude of the angular velocity.
Types of tachometer :
A. DC tachometer :
1. The most common type of DC tachometer contains an iron-core rotor.
Control System Concepts
www.aktutor.in 1–16 C (EN-Sem-5)
er = K r
Fig. 1.10.1.
B. AC tachometer :
Reference
winding Vr cos ct
1
2 2
1 Quadrature
winding
Rotor
Vout
Fig. 1.10.2. AC tachometer.
1. It consists of two stator field coils or windings placed 90° electrical apart
or mounted at right angles to each other. These windings are reference
winding and quadrature winding.
2. A rotating element, i.e., a rotor, is placed in the air gap between the
magnetic structure. Generally the rotor is a thin aluminium cup to
minimize the losses and low inertia.
3. A sinusoidal voltage of rated value is applied to the primary winding. A
secondary winding is placed at a 90 degree angle mechanically with
respect to the primary winding.
4. When the rotor shaft is rotated, the magnitude of the sinusoidal output
voltage is proportional to the rotor speed.
Mathematical modeling of tachometers :
1. The dynamics of the tachometer can be represented by the equation
Control System
www.aktutor.in 1–17 C (EN-Sem-5)
d(t)
et(t) = K t = kt(t) ...(1.10.1)
dt
where, et(t) = Output voltage
(t) = Rotor displacement
(t) = Rotor velocity
Kt = Tachometer constant
2. The transfer function of the tachometer is obtained by taking the Laplace
transform of eq. (1.10.1),
Et (s)
= Kts
(s )
er e Power
+ Controller M Load
– amplifier
et T Tachometer
PART-3
Transfer Function Models, Block Diagram Algebra,
Signal Flow Graph, Mason’s Gain Formula.
Questions-Answers
Answer
A. Transfer function : Transfer function is defined as the Laplace
transform of output to the Laplace transform of input when all initial
conditions are set to zero.
L c(t) C (s)
T(s) =
L r t R(s) All initial condition are zero
B. Advantages :
1. Stability analysis of the system can easily be carried out.
2. The use of Laplace transform approach allows converting integral-
differential time domain equations to simple algebraic equations.
C. Disadvantages :
1. It is applied only to linear and time invariant systems.
2. Initial conditions do affect the system performance but these are
neglected while determining transfer function and hence lose their
importance.
R1
C
vi(t) R2 vo (t)
i(t)
Fig. 1.12.1.
Answer
1. Taking Laplace transform of circuit components,
R1
+
1/Cs R2
Vi(s) Vo(s)
I(s)
–
Fig. 1.12.2.
Control System
www.aktutor.in 1–19 C (EN-Sem-5)
Answer
A. Branch point : A branch point or a take off point is a point from which
the signal from a block goes concurrently to other blocks or running
points.
B. Summing point : The output of control system is feedback to the
summing point where it is compared with the reference input.
C. Blocks :
1. A block diagram of a system is a pictorial representation of the functions
performed by each component and of the flow of signals.
2. A system is represented by using block diagrams. These are easy to
construct even for the complicated system. It is easy to visualize the
function of individual element in the block diagram of a system.
Take-off point
R(s)
+ G(s) C(s)
–
Summing point H(s)
Fig. 1.13.1. Schematic diagram of closed loop system.
D. Block diagram reduction rules :
1. Series connection :
R(s) Y(s) C(s) R(s) C(s)
G1(s) G2(s) G1(s)G2(s)
Control System Concepts
www.aktutor.in 1–20 C (EN-Sem-5)
2. Parallel connection :
+ C(s)
R(s) G1(s)
± R(s) C(s)
G1(s) ± G 2(s)
G2(s)
3. Positive feedback connection :
R(s) + C(s)
G1(s)
+ R(s) G1(s) C(s)
1 – G1(s)H1(s)
H 1(s)
a + (a + b)G a + (a + b)G
G G
+ +
b b
G
7. Interchanging input :
+ +
a + + d a d
+ + + +
c
b c b
8. Multiplication by (–1) to the input :
+ a +
a c c
– +
–1 b
b
Control System
www.aktutor.in 1–21 C (EN-Sem-5)
R G C R G C
R
R 1/G
R G C R G C
R G
R
G3
R – +
+ + G1 G2 G4 G6
– –
+
H1 G5
H2
Fig. 1.14.1.
Answer
G1 G2
1. R + G3 – G4 + G5 G6
– 1+ G1 G2 H 1
H2
G1 G2 (G 3 – G 4 + G 5)
2. R + G6
– 1+ G1 G2 H 1
H2
G1 G2 (G3 – G 4 + G5)
R G6
3. [ 1+ G1 G2 H 1 + G1 G2 H 2(G3 – G4 + G 5)]
Que 1.15. Reduce the block diagram shown in Fig. 1.15.1 to a single
block representation.
G2
+ + + + +
R G1 G3 G4 C
– – – –
H1 H2
H3
Fig. 1.15.1.
Answer
Step 1 :
G 2/G 1
+
G1 + + a C
G3 G4
R 1 + G 1H 1 – – +
H2
H3
H 3/G 3
H2
+
G1 + G 3 G4
R C
1 + G1H1 – 1 + G 3H 2 +
H3/G4
G2(1 + G3H2)
G 1G 3G 4
+ +
G1 + G 3G 4
C
1 + G 1H 1 – 1 + G 3H 2
H 3/G 4
G1 G3 G4 + G2 + G2 G3 H2
R C
1 + G 1H 1 + G 3 H2 + G 4 H 3 + G1 G3 H1 H2 + G1 G4 H 1 H 3
Que 1.16. Determine the transfer function C(s)/R(s) for the block
diagram shown in Fig. 1.16.1.
H3
S1 + – T1 + S
3
R(s) G1 G2 G3 G4 C(s)
– – T2
S2
H1 H2
Fig. 1.16.1.
Answer
Step 1 : Shift S2 before G1, we get
H3 / G1
S1 – + S3
R(s)
G1 G2 G3 G4
– + – T2 C(s)
S2
H1 H2
Step 2 : G1 and G2 are in cascade.
H3 / G1
S1 – S3
R(s) + C(s)
G1 G2 G3 G4
– + – T2
S2
H1 H2
H3 / G1
R(s) – + + S3 C(s)
+ G1 G2 G3 G4
– S1 – T2
S2
H1 H2
H3 / G1
R(s) – G1 G2 + S3 C(s)
+ G3 G4
1+G 1G 2H1 – T2
S2
H2
H 3/G1 G4
– S3
+ G1 G2 + C
G3 G4
R(s) 1+G1 G2 H1 – T2
S2 H2
H3/G1 G4
–
G1 G2 G3 G4
C(s)
R(s) + 1+G1 G2 H 1 1+G3 G4 H2
C (s ) G1G2G3G4
=
R( s) (1 G1G2 H1 ) (1 G3G4 H2 ) G2G3 H3
G4
R(s) C(s)
+ + G1 + G2 G3 ++
– – –
H2
H1
Fig. 1.17.1.
Answer
Step 1:
G4/G2
R(s) + C(s)
+ + G1 +
–
G2 G3 +
– –
H2
H1
Step 2:
R(s) + G4 C(s)
+ G1 + G2 G3 +
– – – G2
H2
H1
Step 3 :
L1
R(s) G4
+ + G1 –
G2 G3 + C(s)
– – G2
H2
G
H1 G3 + 4
G2
R(s) G2G3 + G4
+ +– G1 C(s)
– 1 + (G2G3 + G4)H2
H 1G2
G3G2 + G4
Control System Concepts
www.aktutor.in 1–26 C (EN-Sem-5)
Answer
1. The overall gain can be determined by Mason’s gain formula given
below
k
Pk k
T=
k1
Answer
1. Source node : The node having only outgoing branches is known as
source or input node. e.g., x0 is source node.
2. Sink node : The node having only incoming branches is known as sink
or output node. e.g., x5 is sink node.
3. Chain node : A node having incoming and outgoing branches is known
as chain node. e.g. x1, x2, x3 and x4.
G42
G33
G13 G35
Fig. 1.19.1.
x1 x2 1 x1 x2 1
Source Sink
(a) ( b)
Fig. 1.19.2.
L1 L2 L3
10. Loop gain : The product of all the gains of the branches forming a loop
is called loop gain. For a self loop, gain indicated along it is its gain.
Generally such loop gain are denoted by ‘L’ e.g. L1, L2 etc.
G1 G2 G3 x4
x1 x2 L1 x3
– H1
Fig. 1.19.4.
Loop gain, L1 = – G2H1.
Que 1.20. Find the transfer function of the signal flow graph
shown in Fig. 1.20.1, using Mason’s gain formula.
G10
G9 G7
X2 G8
X3 X4 X5 X6
X1 X7
G1 G2 G3 G4 G5 G6
–H 1 – H2 –H3
H4
Fig. 1.20.1.
Answer
1. Forward path and gains are :
P1 = G1 G2 G3 G4 G5 G6, 1 = 1
P2 = G1 G9 G4 G5 G6, 2 = 1 – G7
P3 = G1 G9 G10 G6, 3 = 1 – G7 – G8 + G7 G8
P4 = G1 G2 G3 G10 G6, 4 = 1 – G8
Control System
www.aktutor.in 1–29 C (EN-Sem-5)
X7
Pk k
k 1
3. Using Mason’s gain formula, T(s) = =
X1
G1G2G3G4G5G6 G1G9G4 G5G6 (1 G7 ) G1G9G10G6 (1 G7 G8 G7G8 )
=
1 (G2 H1 G4 H2 G5 H3 G3G4 G5 H4 G7 G8 H1G9G10 H4
H1G9 G5 H4G4 H4 G3G10 H2G10 H3 ) (G2 H1G4 H2 G2 H1G5 H3
G2 H1G8 G2 H1 H2G10 H3 H1G9G4 G5 H4 G8 G7G8 G7G4 H2
G7G5 H3 G7 H2G10 H3 H4G3G10G7 H4 G3G10G5 H3 )
Que 1.21. Draw the signal flow graph and determine the overall
transfer function of the block diagram shown in Fig. 1.21.1.
G4
R + + + + – C
G1 G2 G3
– – –
H1
Fig. 1.21.1.
Answer
1. Signal flow graph of given block diagram :
Control System Concepts
www.aktutor.in 1–30 C (EN-Sem-5)
– G4
1 1 2 G1 1 3 G2 4 G 5
3
R C
–H 1
–1
–H 1
Fig. 1.21.2.
PART-4
Open Loop and Closed Loop System and their Sensitivity Analysis.
Questions-Answers
Answer
i. Overall Gain :
For open loop :
C ( s)
T(s) = = G(s)
R( s)
R(s) G(s) C(s)
Fig. 1.22.1.
For closed loop :
C(s) G(s)
T(s) =
R(s) 1 G(s)H (s)
+
R(s) G(s) C(s)
–
H(s)
Fig. 1.22.2.
So, for a closed loop control system the gain G(s) reduces by a factor
1/[1 + G(s)H(s)] as compared to open loop control system.
ii. Stability :
1. Consider an open loop system with overall transfer function as
K
R(s) C(s)
s+1
Fig. 1.22.3.
K
G(s) =
sT
The open loop pole is located at s = – T.
2. If unity negative feedback is introduced in the system, the overall
transfer function of a closed loop system becomes
K
C ( s) sT K
=
R( s) K s (K T)
1
sT
Control System Concepts
www.aktutor.in 1–32 C (EN-Sem-5)
× ×
–T Re – (K + T) Re
(a) (b )
Fig. 1.22.4.
4. Thus, the feedback controls the time response i.e., dynamics of the
system by adjusting location of its poles. The stability of the system
depends on the location of poles in s-plane. Thus it can be concluded
that the feedback affects the stability of the system.
iii. Noise and Disturbances :
1. Every control system has some non-linearities present in it. The
dominant non-linearities like friction, dead zone, saturation etc., affect
the output of the system adversely.
2. Some external disturbance signals also make the system output
inaccurate.
3. The examples of such external disturbances are high frequency noise
in electronic applications, thermal noise in amplifier tubes, wind gusts
on antenna of radar system etc.
4. The disturbance may be in the forward path, feedback path or output
of a system.
a. Disturbance in the forward path :
Td(s)
+
R(s) + G1 (s) + G2 (s) C(s)
–
H(s)
Fig. 1.22.5.
1. G(s) = G2(s)
H(s) = – G1(s) H(s)
C (s) G2 (s)
=
Td (s) 1 [G2 (s) ( G1 (s) H(s))]
C (s) G2 (s)
=
Td (s) 1 G1 (s) G2 (s) H (s)
Td (s) G2 (s)
C(s) =
1 G1 (s) G2 (s) H (s)
2. In the denominator assume that 1 << G1(s) G2(s) H(s), hence we get
Td (s)
C(s) =
G1 (s) H (s)
Thus to make the effect of disturbance on the output as small as
possible the G1(s) must be selected as large as possible.
Control System
www.aktutor.in 1–33 C (EN-Sem-5)
R(s)
+ G1(s) G2(s) C(s)
–
+
H2(s) H 1(s)
+
T d(s)
Fig. 1.22.6.
H(s)
Fig. 1.22.7.
C (s) 1 1
=
Td (s) 1 [ G(s) H (s)] 1 G(s) H (s)
2. For large value of G(s) H(s), 1 in denominator can be neglected
[ 1 << G(s) H(s)]
Td (s)
C(s) =
G(s) H (s)
3. Thus if disturbance is affecting the output directly then by changing
the values of G(s), H(s) or both, the effect of disturbance can be
minimized.
The feedback minimizes the effect of disturbance signals occurring in
the control system.
Control System Concepts
www.aktutor.in 1–34 C (EN-Sem-5)
R(s) s+1 K
+ C(s)
– s+5 s(s + 1)
Fig. 1.23.1.
Answer
Sensitivity : Sensitivity is the change in variable due to variation in
parameters of control system.
e.g.
T
% Change in T 100
SGT T
% Change in G = G
100
G
Numerical :
K
Given : G(s) = , H(s) = 1
s(s 5)
4. Putting s = j, K = 1
j( j 5)
= 2
( j)2 5 j 1
5. Again, putting = 5
j5 ( j5 5) j 25 ( j 1)
=
[ 25 j 25 1]2 [ j25 24]2
( 25 j 25)
=
( j 25 24)2
252 252
| SKT | = = 1.02
252 24 2
H(s)
Fig. 1.24.1.
Answer
i. Sensitivity :
T
100
% Change in T
SGT = T
% Change in G G
100
G
G dT
SGT =
T dG
C G
We know T=
R 1 GH
G (1 GH ) GH
SGT = 2
G (1 GH )
1 GH
1 1
SGT = (1 GH )
(1 GH )2 1 GH
So sensitivity of feedback system (with respect to G) is reduced by factor
of 1 + GH.
ii. Error : E(s) = R(s) – H(s) C(s) ...(1.24.1)
Control System Concepts
www.aktutor.in 1–36 C (EN-Sem-5)
G(s)
But, C(s) = R(s)
1 G (s ) H ( s )
H (s) G(s)
E(s) = R(s) R(s)
1 G(s) H (s)
1 G (s) H ( s) H ( s)G (s)
= R( s)
1 G ( s ) H (s )
1
E(s) = R(s) ...(1.24.2)
1 G ( s) H ( s )
1
Hence, error E(s) is reduced by
1 G( s) H ( s)
G
iii. Gain : Gain = 1 GH
Hence, Gain reduces by factor of 1 + GH.
0.50
Fig. 1.25.1.
Answer
R1
C
vi(t) R2 vo (t)
i(t)
Fig. 1.
Ans. Refer Q. 1.12, Unit-1.
Fig. 2.
Ans. Refer Q. 1.17, Unit-1.
G10
G9 G7
X2 G8
X3 X4 X5 X6
X1 X7
G1 G2 G3 G4 G5 G6
–H 1 – H2 –H3
H4
Fig. 3.
Ans. Refer Q. 1.20, Unit-1.
R + + + + – C
G1 G2 G3
– – –
H1
Fig. 4.
Control System
www.aktutor.in 1–39 C (EN-Sem-5)
H(s)
Fig. 5.
Ans. Refer Q. 1.24, Unit-1.
0.50
Fig. 6.
Ans. Refer Q. 1.25, Unit-1.
Control System
www.aktutor.in 2–1 C (EN-Sem-5)
PART-1
Standard Test Signals, Time Response of First and
Second Order Systems, Time Response Specifications.
Questions-Answers
Que 2.1. Explain various standard test signals, and also find
Answer
1. Unit step : Signals which start at time t = 0 and have magnitude of
unity are called unit step signals.
They are represented by a unit step function u(t).
u(t)
t
0
Fig. 2.1.1. Unit step.
Control System
www.aktutor.in 2–3 C (EN-Sem-5)
r(t) = mt
m = tan
t
0
Fig. 2.1.2.
3. Unit impulse : Signals which act for very small time but have large
amplitude are called unit impulse functions.
They are represented by (t).
They are defined mathematically as,
0; t0
(t) =
1 ; t0
(t) dt =1
x(t)
(t)
–t t
Fig. 2.1.3.
4. Unit Parabolic Signal :The continuous-time unit parabolic function
p(t), also called acceleration signal starts at t = 0, and is defined as :
t2
; for t 0
p(t) = 2
0 ; for t 0
t2
or p(t) = u(t)
2
Time Response Analysis
www.aktutor.in 2–4 C (EN-Sem-5)
p(t)
0 t
Fig. 2.1.4. Unit parabolic signal.
Relation :
1. Relation between impulse and step signal :
d
(t) = u(t)
dt
2. Relation between step and ramp signal :
d
u(t) = r (t )
dt
3. Relation between ramp and parabolic signal :
d t2
r(t) = .
dt 2
Que 2.2. Discuss the time response of first order system with
unit step, unit impulse and unit ramp inputs.
Answer
1. Consider a first order system with unity feedback.
R(s)
+ 1/sT C(s)
–
Fig. 2.2.1.
1
C ( s) sT 1
2. =
R( s) 1 (1 sT )
1
sT
1
C(s) = R( s) ...(2.2.1)
(1 sT )
3. Response to unit step input :
r(t) = u(t)
1
R(s) =
s
Putting in eq. (2.2.1), we get
1 1
C(s) =
s (1 sT )
Control System
www.aktutor.in 2–5 C (EN-Sem-5)
= 1 1
s 1
s
T
Taking inverse Laplace transform, we have
1 1
L–1 [C(s)] = c(t) = L1 L1
s 1
s
T
c(t) = 1 – e–t/T
4. Response to unit impulse input :
r(t) = (t)
R(s) = 1
Putting in eq. (2.2.1),
1 1 1
C(s) = 1
(sT 1) T 1
s
T
Taking inverse Laplace transform, we have
1 1 1
L–1[C(s)] = L
T 1
s
T
1 t/T
e
c(t) =
T
5. Response to unit ramp input :
r(t) = t
1
R(s) =
s2
Putting in eq. (2.2.1), we get
1 1
C(s) = R(s)
(sT 1) s 2 ( sT 1)
1 T T
= 2
s s 1
s T
Taking inverse Laplace transform, we have
1 1 1
L–1[C(s)] = L1 2 TL1 TL1
s s 1
s
T
c(t) = t – T + Te–t/T
Que 2.3. Give comparison between open loop and closed loop
system. The impulse response of unity feedback closed loop system
is, c(t) = – te– t + 2e– t, find the open loop transfer function.
AKTU 2014-15, Marks 10
Time Response Analysis
www.aktutor.in 2–6 C (EN-Sem-5)
Answer
Comparison between open loop and closed loop System : Refer
Q. 1.2, Page 1 –3C, Unit-1.
Numerical :
Given : r(t) = δ(t) , c(t) = – te– t + 2e– t
To Find : Transfer function.
Answer
1. Consider the second order system with unity feedback.
R(s) + 2n
s(s + 2 n) C(s)
–
Fig. 2.4.1.
Put d = n 1 2
1 s 2 n
C (s) = ...(2.4.4.)
s (s n )2 d2
4. Rewrite eq. (2.4.4)
1 s n n d
C(s) = ...(2.4.5)
s (s n )2 d2 d (s n )2 2d
5. Taking inverse Laplace transform of eq. (2.4.5),
t n nt
c(t) = 1 e n cos d t e sin d t
d
Put d = n 1 2
Time Response Analysis
www.aktutor.in 2–8 C (EN-Sem-5)
ent
c(t) = 1 [( 1 2 )cos d t sin d t ] ...(2.4.6)
1 2
6. Put sin = 1 2 ,
cos =
e n t
c(t) = 1 (sin cos d t cos sin d t)
1 2
e n t
c(t) = 1 sin( d t ) ...(2.4.7)
1 2
where d = n 1 2
1 2
and = tan 1
7. Eq. (2.4.7) is rewritten as
1 2
e n t
c(t) = 1 sin ( n 1 2 )t tan 1 ...(2.4.8)
1 2
c(t)
<1
t
Fig. 2.4.2. Underdamped oscillations.
b. Undamped case (= 0) :
Put = 0
Eq. (2.4.8) becomes,
c(t) = 1 – sin (nt + /2)
c(t) = 1 – cos nt
Control System
www.aktutor.in 2–9 C (EN-Sem-5)
c(t)
=0
t
Fig. 2.4.3. Undamped oscillations.
c(t)
2%
1
98% =1
0 t
4T
Fig. 2.4.4. Critically damped.
d. Overdamped case ( > 1) :
2n
C (s) =
s n
n 1 s n – n 2 1 s
2
Taking inverse Laplace both side,
2 1 n t 2 1 n t
e e
At > 1 c(t) = 1
2 2 1( 2 1) 2 2 1.( 2 1)
neglect
2 1 n t
e
c(t) = 1
2 2 1( 2 1)
c(t)
>1
t
O
Fig. 2.4.5. Overdamped.
Time Response Analysis
www.aktutor.in 2–10 C (EN-Sem-5)
2
AKTU 2013-14, Marks 6
3
Answer
1.
[r(t) – (t)] T
Transducer
Fig. 2.5.1.
2. The rotation of motor can be described by the differential equation
.. .
J (t) + B (t) = Tm + Tw ...(2.5.1)
where, Tm = Torque developed by motor
Tw = Wind torque (disturbance)
B= Coefficient of viscous friction
J = Moment of inertia of motor and load
3. The torque developed by the motor is assumed to be proportional to
u(t), the input voltage to the motor, so that
Tm(t) = Ku(t) = K[r(t)– (t)]
4. Putting Tm(t) in eq. (2.5.1)
J
(t) B (t) = K[r(t) – (t)] [Tw = 0]
J
(t) B (t) + K(t) =Kr(t) ...(2.5.2)
5. Taking Laplace transform on both sides
Js 2(s) + Bs (s) + K(s) = Kr(s)
[Js 2 + Bs + K] (s) = Kr(s)
(s) K K/J
=
r (s) [ Js 2 Bs K ] 2 B K
s s
J J
6. Putting K/J = n2
Control System
www.aktutor.in 2–11 C (EN-Sem-5)
B
= 2n
J
( s) C( s) n2
= 2
r (s) R( s) s 2 n s n 2
Thus, the servomechanism system is a second order system.
PART-2
Steady State Errors and Errors Constants.
Questions-Answers
Answer
1. The difference between the steady-state response and desired reference
give steady state error.
2. If the actual output of control system during steady state error deviates
from the reference input, then system possess steady state error.
3. Steady error helps in determining the accuracy, so the steady state
error should be minimum.
4. The steady state performance of a control system is assessed by the
magnitude of the steady state error possessed by the system and the
system input can be step or ramp or parabolic.
5. The magnitude of steady state error in a closed loop control system
depends on its open loop control function. i.e., G(s)H(s) of the system.
Time Response Analysis
www.aktutor.in 2–12 C (EN-Sem-5)
R(s) + E(s)
G(s) C(s)
B(s) –
H(s)
Answer
i. Unit step input :
1; t 0
Here r(t) =
0; t 0
Control System
www.aktutor.in 2–13 C (EN-Sem-5)
1
R(s) =
s
R(s)
E(s) = ...(2.7.1)
1 G( s) H ( s)
1
E(s) =
s sG (s) H (s)
1
ess = lim sE(s) = lim s
s0 s0 s sG (s) H ( s)
1 1
=
1 lim G (s) H ( s) 1 K p
s 0
where, Kp = lim
s0
G( s) H ( s)
Kp = Position error constant.
1
For step input, ess =
1 Kp
Case 1 : For type ‘0’
Kp = Constant
ess = Constant
Case 2 : For type ‘1’
Kp =
1
ess = = 0
1
Case 3 : For type ‘2’
Kp =
1
ess = =0
1
ii. Ramp input :
t ; t 0
Here r(t) =
0 ; t 0
1
R(s) =
s2
1
E(s) =
s 2 s 2G (s) H ( s)
1
ess = lim sE( s) = lim s
s0 s0 s2 s2G( s) H (s)
1
=
lim sG(s) H (s)
s0
1
ess =
Kv
Time Response Analysis
www.aktutor.in 2–14 C (EN-Sem-5)
Kv = lim sG( s) H ( s)
s0
Kv = Velocity error constant.
Case 1 : For type ‘0’
Kv = lim sG( s) H ( s) = 0
s0
1
ess = =
Kv
Case 2 : For type ‘1’
Kv = constant
1
ess = = constant
Kv
Case 3 : For type ‘2’
Kv =
ess = 0
iii. Parabolic input :
(1 / 2) t 2 ; t 0
Here r(t) =
0 ; t0
1
R(s) =
s3
1
E(s) =
s3 s3G( s) H ( s)
1
ess = lim sE(s) lim
s0 s 0 s 2 s 2 G ( s ) H (s )
1 1
ess =
lim s2G( s) H ( s) K a
s 0
10
i. G(s) = and H (s) = 0.7
s2
5 2
ii. G(s) = , H(s) = 0.6 AKTU 2013-14, Marks 6
2
s + 3s + 5 3
Answer
Answer
K
Given : G(s) = , H(s) = 1
s(s 2) (s2 2s 5)
t2
r(t) = 2 4 t
2
To Find : Steady state error, ess.
K
1. G(s) H(s) =
s(s 2)( s2 2s 5)
KP = lim
s0
G( s) H ( s)
=
Kv = lim s G( s) H ( s)
s0
sK K
= lim =
s( s 2)( s2 2s 5) 10
s0
Ka = lim s 2 G( s) H ( s) = 0
s0
t2
2. Now, r(t) = 2 + 4t +
2
A1 = 2 (Constant)
A2 = 4 (Coeff. of t)
A3 = 1 (Coeff. of t2/2)
A1 A A 2 4 1
ess = 2 3 =
1 K p Kv Ka K 0
10
= 0 + 40/K + =
Que 2.10. Discuss different type of test signal used for analysis of
control system in time domain.
The reference input to a unity feedback system is shown in
Fig. 2.10.1. The open loop transfer function of the system is
400(s 1)
G(s) =
( s 2)( s 8)
Calculate the steady state error.
r(t)
2
t
10 20 30 40
Fig. 2.10.1.
Answer
Types of test signal : Refer Q. 2.1, Page 2–2C, Unit-2.
Numerical :
400( s 1)
Given : G(s) = , H(s) = 1
( s 2)( s 8)
To Find : Steady state error.
1. Laplace of periodic waveform using formula :
R1 (s)
R(s) =
1 e TS
2. First we have to find Laplace of r(t),
r(t)
t
10
Fig. 2.10.2.
By integration method
10
R1(s) = 2 e st dt
0
2 10
= e st
( s) 0
2
R1(s) = (1 e 10 s )
s
R1 (s)
4. Here, T = 20, So R(s) =
1 eTS
2
(1 e10 s )
R(s) = s
(1 e20 s )
2
(1 e 10 s )
R(s) = s
(1 e10 s )(1 e10 s )
2
R(s) =
s(1 e10 s )
5. Using formula
sR( s)
ess = lim
s0 1 G( s) H (s)
2
s
s (1 e10 s )
= lim
s0 400( s 1)
1
( s 2)( s 8)
Time Response Analysis
www.aktutor.in 2–18 C (EN-Sem-5)
2
= lim (1 e10 s )
s 0 ( s 2)( s 8) 400( s 1)
( s 2)( s 8)
2(s 2)( s 8)
= lim
s 0 [( s 2)( s 8) 400(s 1)](1 e10 s )
2(0 2)(0 8)
=
[(0 2)(0 8) 400(0 1)](1 e0 )
2 28 32
=
2 8 400 416
ess = 0.0769
Que 2.11. For the system shown in figure, determine the type of
system, error coefficient and the error for the following inputs :
i. r(t) = 6
ii. r(t) = 8t
iii. r(t) = 10 + 4t + 15t2
+ + 20
R(s) 20 2 C(s)
– – (s + 2)(s + 5)
2(s + 1)
Answer
Step 1 :
+ 20
R(s) 20 2 C(s)
– s 2 +47s + 50
Step 2 :
+ 800
R(s) C(s)
– s 2 +47s + 50
800
So, G(s) = and H(s) = 1
s 2 47 s 50
800
G(s) H(s) =
s0 (s 2 47 s 50)
Control System 2–19 C (EN-Sem-5)
www.aktutor.in
Type of system is zero.
i. r(t) = 6
6 1
So, R(s) = A = 6, which is coefficient of
s s
800
Now, Kp = lim G( s) H ( s) =
s0 50
A 6 6 50
Now, ess = = 0.353
1 K p 1 800 850
50
ii. r(t) = 8 t
8 1
So, R(s) = A = 8, which is coefficient of 2
s2 s
Now, Kv = lim sG( s) H ( s)
s0
s 800
= lim =0
s0 s 2 47 s 50
A 8
Now, ess = =
Kv 0
iii. r(t) = 10 + 4t + 15t2
t2
= 10 + 4t + 30
2
Now, So, A1 = 10, A2 = 4, A3 = 30
K = lim G (s) H (s)
p s0
800
= = 0.941
850
Kv = lim
s0
sG( s) H (s)
=0
Ka = lim s2 G(s) H (s)
s 0
=0
A1 A A
Now, ess = 2 3
1 K p Kv Ka
10 4 30
= =
1 0.941 0 0
Que 2.12. For a unity feedback system having
35( s 4)
G(s) =
s(s + 2)( s + 5)
find (i) the type of the system, (ii) all error coefficients and (iii)
errors for ramp input with magnitude 5. AKTU 2017-18, Marks 10
Time Response Analysis
www.aktutor.in 2–20 C (EN-Sem-5)
Answer
35( s 4)
Given : G(s) = , H(s) = 1
s(s 2)(s 5)
To Find : Type of the system, All error coefficients and Errors for
ramp input with magnitude 5.
i. Type : System is type 1 system. As we have s 1 in common in
denominator.
ii. Error coefficients :
Kp = lim G (s) H (s)
s0
35(s 4)
= lim 1=
s0 s( s 2)(s 5)
Kv = lim
s0
s G(s) H ( s)
s 35 ( s 4)
= lim 1
s0 s (s 2) (s 5)
35 4
= = 14
10
Ka = lim s 2 G( s) H ( s)
s0
35( s 4)
= lim s2 1 = 0
s0 s( s 2)(s 5)
iii. Error : Magnitude, A = 5 and Kv = 14
A 5
ess =
K v 14
ess = 0.357.
K ( s + 2)
G(s)= . Find the static error constant and ess due to
s( s 3 + 7 s 2 + 42s)
Answer
K ( s 2)
Given : G(s) = , r(t) = t2u(t)
s(s3 7 s 2 42 s)
To Find : Ka and ess.
t2
1. r(t) = t2u(t) = 2 u( t)
2
A = 2. [Coeff. of t2/2]
Control System
www.aktutor.in 2–21 C (EN-Sem-5)
K ( s 2)
2. Now, G(s)H(s) =
s(s3 7 s 2 42 s)
K ( s 2)
= 2 2
s ( s 7 s 42)
3. Now, Ka = lim s 2 G( s) H(s)
s0
s2 K (s 2)
= lim
s 0 s ( s2 7 s 42)
2
2K K
=
42 21
A 2 42
4. ess =
Ka K K
21
PART-3
Time Response Specifications, Design Specifications
of Second Order Systems.
Questions-Answers
Answer
1. Delay time (td) : It is the time required for the response to reach 50 %
of the final value in first time.
2. Rise time (tr) : It is the time required for the response to rise from
10 % to 90 % of its final value for overdamped system and 0 to 100 % for
underdamped systems.
1 2
tan 1
tr =
n 1 2
3. Peak time (tp) : The peak time is the time required for the response to
reach the first peak of the time response or first peak overshoot.
tp =
n 1 2
4. Maximum overshoot (Mp) : It is the normalized difference between
the peak of the time response and steady output. The maximum percent
overshoot is defined as
c(t p ) c()
Maximum percent overshoot = 100
c()
% Mp = e 12
100
c(t)
c(t)max
Mp Steady state
1.0
1.5
2 3
t
td 2 2 2
tr n 1 – n 1 – n 1 –
tp
ts
Fig. 2.14.1.
Control System
www.aktutor.in 2–23 C (EN-Sem-5)
5. Settling time (ts) : The setting time is the time required for the response
to reach and stay within the specified range (2 % to 5 %) of its final value.
4
ts =
n
6. Steady state error (ess) : It is the difference between actual output and
desired output as time ‘t’ tends to infinity.
ess = lim r (t) c(t)
t
Answer
A. The rise time, tr :
1. The time response of a second order control system is given by
ent
c(t) = 1 sin[(n 1 2 ) t ]
1 2
2. Time tr is the time at which response reaches 100 % i.e., c(t) = 1.
entr
1 = 1 sin [(n 1 2 )tr ]
1 2
entr
sin [(n 1 2 )tr ] = 0
1 2
sin [(n 1 2 )tr ] = 0 [sin n = 0, n = 1]
[(n 1 2 )tr ] =
3. tr =
n 1 2
1
1 2
where, = tan
B. Maximum overshoot, Mp :
1. Mp = c(t)max – 1
c(t)max 1
% Mp = 100
1
ent
2. c(t) = 1 sin [(n 1 2 )t ]
2
1
Time Response Analysis
www.aktutor.in 2–24 C (EN-Sem-5)
dc(t) ent
= n 1 2 cos [n 1 2 )t ]
dt 1 2
( n ) ent
sin [(n 1 2 )t ]
1 2
dc(t)
3. For maximum, Put =0
dt
ewnt
{n 1 2 cos [(n 1 2 )t ]
1 2
n sin [(n 1 2 )t ]} = 0
2
n 1 2 cos [(n 1 2 )t ] = n sin [(n 1 )t ]
1 2
tan [(n 1 2 )t ] =
2
4. Putting, = tan – 1 1
1 2 1 2
tan (n 1 2 )t tan 1 = ...(2.15.1)
5. The general solution of eq. (2.15.1) is
(n 1 2 )t = n ...(2.15.2)
where n = 0, 1, 2...
At n = 1 , tp =
n 1 2
6. Putting t = tp in c(t),
n t p
e
c(t)max = 1 sin [( n 1 2 )t p ]
1 2
n
n 1 2
e
= 1 sin ( n 1 2 )
2 2
1 n 1
1 2 1 2
e e
= 1 sin ( ) 1 ( sin )
1 2 1 2
1 2
As, = tan 1 , sin = 1 2
1 2
e
Hence, c(t)max = 1 1 2
1 2
Control System
www.aktutor.in 2–25 C (EN-Sem-5)
2
c(t)max = 1 e 1
7. Mp = c(t)max – 1
Mp = 1 e 1
1 2
1 2
Mp = e
12
% Mp = e 100
Que 2.16. Derive expression for resonant frequency and resonant
Answer
1. For second order system
C ( s) 2n
= 2
R( s) s 2 n s 2n
2. Putting s = j
C( j) n2
=
R( j) – 2 n j 2n
2
1
M= 2
1 – j 2
n n
Let =x
n
1 1
M=
1 x 2 j 2x (1 x ) 4 2 x 2
2 2
dM
3. To maximize magnitude, put = 0,
dx
dM d
= [(1 – x 2 )2 4 2 x 2 ]–1/ 2
dx dx
1 d
= – [(1 – x 2 )2 4 2 x 2 ]–3/ 2 [(1 – x 2 )2 4 2 x2 ]
2 dx
1 1
= – [2(1– x 2 ) (– 2 x) 8 2 x ]
2 [(1 – x 2 )2 4 2 x 2 ]3/ 2
=0
4x[x2 + 22 – 1] = 0
As, x 0
x2 = 1 – 22
Time Response Analysis
www.aktutor.in 2–26 C (EN-Sem-5)
x= 1 – 2 2
r = n 1 – 2 2
1
4. |Mr| =
(1 x ) 4 2 x 2
2 2
1
|Mr| =
(1 ( 1 2 )) 4 2 ( 1 2 2 )2
2 2
1
|Mr| =
4 4 4 2 (1 2 2 )
1
Mr =
2 1 2
Answer
Given : G(s) = , H(s) = 1, Mp = from 0.6 to 0.2
s(1 s)
1 1
To Prove : 43.
1 2
1. Mp reduces from 0.6 to 0.2 due to change in only.
Transfer function
C (s ) G ( s) / s (1 s) /
= = = = 2 1
R( s) 1 G (s) H ( s) s 2
s
1 s s
s(1 s)
2. Compare denominator by s2 + 2 ns + n2 = 0
1 1
2 n = =
2 n
Control System
www.aktutor.in 2–27 C (EN-Sem-5)
n =
1
=
2 2
3. Maximum overshoot,
1 2
Mp = e
1
2
loge 0.6 = loge e 1 1
1
loge 0.6 = –
1 12
1
– 0.5 = –
1 12
(1 – 12) 0.25 = 12 2
(1 12 ) 9.87
= 39.48
12 0.25
39.48 12 + 12 – 1 = 0
40.48 12 = 1
1
1 = 0.157 =
2 1
1 = 10.142
1 – 1 = 9.142 ...(2.17.1)
4. Similarly for 2 at 2 = 0.2, we get
1
2 = 0.4538 = ...(2.17.2)
2 2
2 = 1.214
2 – 1= 0.214
5. Dividing eq. (2.17.1) and (2.17.2), we get
1 1
= 42.72
2 1
43
( s 2)
function G(s) = . Determine rise time, peak time and settling
s ( s 1)
Answer
(s 2)
Given : G(s) = , H (s) = 1
s(s 1)
To Find : tr, tp and ts.
1. The characteristic equation for unity feedback control system is given
by
1 + G(s) H(s) = 0
(s 2)
1+ =0
s(s 1)
2
s +s+s+2= 0
s2 + 2s + 2 = 0
2. Comparing by standard second order characteristic equation
s2 + 2ns + 2n = 0
n = 2 rad/sec
2n = 2
1
= 0.707
2
3. Rise time, tr =
d
1 2
= tan 1 = 45° = 0.7853 (in radian)
d = n 1 2 = 2 1 (0.707)2 = 1
tr = 2.3562 sec
4. Peak time, tp = 3.1415 sec
d
4 4
5. Setting time, ts (for 2%) = =
n (0.707) ( 2)
ts = 4 sec.
Que 2.19. Derive the expression for second order system response
when subjected to unit impulse input for damping ratio () < 1.
An unity feedback system is characterized by an open loop transfer
function.
K
G(s) =
s(s 10)
Determine the gain ‘K’ so that the system will have a damping
ratio of 0.5. For this value of ‘K’, determine the settling time, peak
overshoot and time to peak overshoot for a unit step input.
AKTU 2016-17, Marks 10
Control System
www.aktutor.in 2–29 C (EN-Sem-5)
Answer
Derivation for second order system response :
1. For unit impulse function, R(s) = 1 and the output of a second order
system is given by
n 2
C(s) = ...(2.19.1)
s 2n s n 2
2
n2 n 1 2
C(s) = ...(2.19.2)
n 1 2 [(s n )2 ( n 1 2 )2 ]
3. Taking inverse Laplace transform on both sides of eq. (2.19.2),
n
c(t) = e nt sin [( n 1 2 )t ] ...(2.19.3)
1 2
4. For < 1,
n
c(t) = e nt sin [( n 1 2 )t ] ...(2.19.4)
1 2
c(t) <1
0
t
Fig. 2.19.1.
5. The time response for < 1 is decaying exponential oscillations and the
output at times goes negative also.
Numerical :
K
Given : G(s) = , = 0.5, H(s) =1
s( s 10)
To Find : K, ts, tp and Mp.
C ( s) G (s )
1. =
R( s) 1 G (s) H ( s)
K
= ...(2.19.5)
s 2 10 s K
Time Response Analysis
www.aktutor.in 2–30 C (EN-Sem-5)
PART-4
Proportional, Derivative, Integral and PID Compensations, Design
Considerations for Higher Order System and Performance Indices.
Questions-Answers
Answer
Proportional control :
1. In proportional control the actuating signal for the control action in a
control system is proportional to the error signal.
2. The error signal being the difference between the reference input signal
and the feedback signal obtained from the output.
3. For the system considered as shown in Fig. 2.20.1, the actuating signal
is proportional to the error signal, therefore, the system is called
proportional control system.
4. Consider a second order system where controller input is error itself
and proportional constant is K = 1 as shown in Fig. 2.20.1.
H(s)
Fig. 2.20.1.
2n
G(s)H(s) =
s( s 2 n )
C ( s) 2n
= 2
R( s) s 2 n s 2n
For this system damping ratio is and natural frequency n.
5. And for steady state error
Kp = lim G( s) H ( s)
s0
n
and Kv = lim s G( s) H ( s)
s 0 2
Effects :
1. Steady state error is reduced.
2. Disturbance signal rejection occurs.
3. Relative stability is improved.
OR
What is the effect of PD (Proportional Derivative) Controller on
steady state error due to a unit ramp input in second order system ?
Prove mathematically. AKTU 2013-14, Marks 05
OR
Write a short note on PD controller and synchros.
AKTU 2015-16, Marks 10
Answer
1. A controller in the forward path, which changes the controller output
corresponding to proportional plus derivative of error signal is called PD
controller.
de(t)
2. Output of controller = K e( t) Td
dt
Taking Laplace transform
= K E(s) + sTd E(s) = E(s) [K + sTd]
+ n2
R(s) K C(s)
+ s (s + 2 n)
–
sTd
Fig. 2.21.1.
3. Assume, K = 1, we can write,
(1 sTd ) 2n
G(s) =
s( s 2 n )
C ( s) (1 sTd ) 2n
= 2
R( s) s s [2 n 2n Td ] 2n
4. Comparing denominator with standard form,
2n = 2 n 2n Td
n Td
=
2
n Td
5. Because of this controller, damping ratio increased by factor .
2
Kp = lim G( s) H ( s)
s0
n
Kv = lim s G(s) H (s)
s 0 2
without PD controller,
Control System
www.aktutor.in 2–33 C (EN-Sem-5)
2n
G(s) H(s) = [K =1]
s(s 2 n )
n
Kv = lim sG(s) H ( s)
s0 2
As there is no change in coefficients, error also will remain same.
6. Effects :
i. It increases damping ratio.
ii. ‘TYPE’ of the system remains unchanged.
iii. It reduces peak overshoot.
iv. It reduces settling time.
v. Steady state error remains unchanged.
7. Demerits :
1. It amplifies the noise signal.
2. It may cause saturation effect in the actuator.
3. It does not improve ess because of generation of a zero in the transfer
function.
Synchros : Refer Q. 1.8, Page 1–11C, Unit-1.
Answer
i. Derivative controller : Refer Q. 2.21, Page 2–31C, Unit-2.
Application of PD controller :
1. In an industrial plant, a PD controller is applied in applications where
overshoot cannot be tolerated. Such as batch pH neutralization.
2. In some positioning applications.
ii. Proportional integral controller (PI) :
1. A controller in the forward path, which changes the controller output
corresponding to the proportional plus integral of the error signal is
called PI controller.
Ki K
Taking Laplace = K E( s) E( s) E( s) K i
s s
+ n2
R(s) K C(s)
+ s (s + 2 n)
–
Ki/s
Fig. 2.22.1.
3. Assume, K = 1, we can write,
Ki 2
1 s n
( K s) n2
G(s) = 2 i
s (s 2 n ) s (s 2 n )
i.e., system becomes TYPE 2 in nature.
C ( s) ( K i s) n2
= 3
R( s) s 2 n s2 s 2n K i n2
i.e, it becomes third order.
4. Now as order increases by one, system relatively becomes less stable as
Ki must be designed in such a way that system will remain in stable
condition. Second order system is always stable.
while Kp = lim G( s) H ( s) , ess = 0
s0
5. Hence as type is increased by one, error becomes zero for ramp type of
input i.e., steady state of system gets improved and system becomes
more accurate in nature.
6. Applications of PI controller :
i. Pressure control
ii. Temperature control
7. Effects :
i. It increases order of the system.
ii. It increases TYPE of the system.
iii. Design of Ki must be proper to maintain stability of system.
vi. Steady state error reduces tremendously for same type of inputs.
Answer
PID controller (Proportional, integral and derivative
controller) :
1. It is a close loop system which has feedback control system and it
compares the process variable (feedback variable) with set point and
generates an error signal and according to that it adjusts the output of
system.
Control System
www.aktutor.in 2–35 C (EN-Sem-5)
Ki/s
Fig. 2.23.1.
Kts
Fig. 2.24.1.
Answer
20
s( s 2) 20
1. Inner loop =
20sK t s ( s 2 20 K t )
1
s(s 2)
20
C ( s) s( s 2 20 K t ) 20
2. = 2
R( s) 1 20 s s (2 20 K t ) 20
s(s 2 20 K t )
Time Response Analysis
www.aktutor.in 2–36 C (EN-Sem-5)
3. 2n = 20, n = 20
2m = 2 + 20 Kt
2 20 K t
= = 0.6 (given)
2 20
Kt = 0.1683
Que 2.25. Define about performance index in brief.
Answer
Performance Index : The goodness of the system performance can be
based on the performance index and helps in designing a control system.
a. Integral Square Error (ISE) :
I= e2 (t)dt
0
It provides a good compromise between the rise time to limit the effect
of a large initial error, reduction of peak overshoot and settling time.
When = 0.5, I is minimum. It is difficult to calculate for higher order
system.
b. Integral of Absolute Value of Error :
I= e(t) dt
0
There are fewer penalties for large error and large penalty for small
error. It requires use of computer and has poor sensitivity. A system
based on this criteria has enough damping and transient performance.
c. Integral of Time Multiplied Square Criteria :
I= te2 (t) dt
0
Large initial error is weighted lightly, while errors occurring late are
weighted heavily.
d. Integral of Square Time Multiplied Square Criteria :
I= t 2 e2 (t) dt
0
t2
input is r(t) = 2 + 4t + .
2
Ans. Refer Q. 2.9, Unit-2.
2(s + 1)
Control System 3–1 C (EN-Sem-5)
www.aktutor.in
3 Stability and
Algebraic Criteria
• Root Contour
• Construction of Root Loci
• Effect of Transportation Lag and Root Locus of Non Minimal
Phase System and Effect of Pole-Zero Cancellation
Questions-Answers
Answer
Stability :
1. A system is stable if its response (i.e., the transfer function) approaches
zero as time approaches infinity.
2. In other words, a system is stable if every bounded input produces a
bounded output.
Types of stability :
1. Absolutely stable : If a system output is stable for all variations of its
parameters, then the system is called absolutely stable system. It gives
the information about whether system is stable or unstable.
2. Relative stability : The system is said to be relatively more stable or
unstable on the basis of settling time. System is said to be relatively
more stable if settling time for the system is less than that of the other
system. Relative stability gives degree of stability or how close it to instability.
Control System 3–3 C (EN-Sem-5)
www.aktutor.in
Necessary conditions for absolutely stable :
i. All the coefficients of characteristics equation must have same sign.
ii. There should be no missing term.
iii. All poles of transfer function should be in left half of s-plane.
iv. The degree of denominator polynomial of transfer function is greater or
equal to that of numerator polynomial.
Answer
The Routh-Hurwitz stability criterion is a mathematical test that is a
necessary and sufficient condition for the stability of system.
Routh’s criterion :
1. The necessary and sufficient condition for the system to be stable is
that all the terms in the first column of Routh’s array (i.e. a0, a1, b1, c1,
d1 ……..an) must have same sign.
2. There should not be any sign change in the first column of Routh’s
array.
3. If there is any sign change then
i. System is unstable.
ii. Number of sign changes is equal to the number of poles lying in the
right half of the s-plane.
Method of forming a Routh array :
sn a0 a2 a4 a6 ......
sn–1 a1 a3 a5 a7
sn–2 b1 b2 b3
sn–3 c1 c2 c3
s0 an
i. Coefficients for first two rows are written directly from characteristics
equation,
a1 a2 a0 a3
b1 =
a1
Stability & Algebraic Criteria 3–4 C (EN-Sem-5)
www.aktutor.in
a1 a4 a0 a5
b2 =
a1
a1 a6 a0 a7
b3 =
a1
ii. From 2nd and 3rd row, 4th row can be obtained as
b1 a3 a1 b2
c1 =
b1
b1 a5 a1 b3
c2 =
b1
iii. This process is to be continued till the coefficient for s0 is obtained which
will be an. From this array stability of a system can be predicated.
Special cases :
1. When the 1st term in a row is zero, but all other terms is non-zero then
substitute a small positive number for zero & proceed to evaluate the
rest of the elements. When the 1st column term is zero, it means that
there is an imaginary root.
2. All zero row : In this case, write auxiliary equation from preceding
row, differentiate this equation & substitute all zero row by the coefficient
obtained by differentiating the auxiliary equation. This case occurs when
the roots are in pairs. The system is said to be limitedly stable.
Application : Routh’s criterion can be applied to determine range of
certain parameters of a system to ensure stability. It is usually to find
the range of the open loop gain K for closed loop stability.
Limitations :
1. It is valid only if the characteristics equation is algebraic.
2. If any coefficient of the characteristics equation is complex or contain
power of e then this criterion cannot be applied.
3. It gives information about how many roots are lying in the RHS of the
s-plane but values of the roots are not available. Also it cannot distinguish
between real and complex roots.
Answer
Answer
K
Given : G(s) H(s) =
s(s 1 ) (1 2s)
To Find : Range of K.
1. Using characteristics equation
1 + G(s) H(s) = 0
K
1+ =0
s(s 1)(1 2s)
s(s + 1) (1 + 2s) + K = 0
s(2s2 + 3s + 1) + K = 0
2s3 + 3s2 + s + K = 0
2. Using Routh array :
s3 2 1
s2 3 K
3 – 2K
s1 0
3
0
s K
3. For system to be stable :
3 – 2K
> 0 and K > 0
3
3
0<K<
2
Stability & Algebraic Criteria 3–6 C (EN-Sem-5)
www.aktutor.in
Que 3.5. Explain the working principle of stepper motor with
neat diagram.
The characteris tics equation for feedback control is ,
s3 + 5s2 + 12s + K = 0. Find the range of K for all the roots to lie to the
left of s = 1. AKTU 2016-17, Marks 10
Answer
Stepper motor : Refer Q. 1.9, Page 1–13C, Unit-1.
Numerical :
Given : Characteristics equation is s3 + 5s2 + 12s + K = 0
To Find : Range of K.
1. Putting s = s – 1
(s – 1)3 + 5 (s – 1)2 + 12 (s – 1) + K = 0
s3 – 1 + 3s – 3s2 + 5s2 + 5 – 10s + 12s – 12 + K = 0
s3 + 2s2 + 5s + (K – 8) = 0
2. Routh array :
s3 1 5
s2 2 ( K 8)
10 ( K 8)
s1 0
2
0
s ( K 8)
3. To lie all roots to the left side of s = –1, there should not be any sign
change in first column of Routh array.
10 ( K 8)
>0
2
10 – (K – 8) > 0
18 – K > 0
K < 18 ...(3.5.1)
Also, K–8> 0
K>8 ...(3.5.2)
From eq. (3.5.1) and (3.5.2), range of K is
8 < K < 18
Que 3.6. For a system having characteristic equation 2s4 + 4s2 +
1 = 0, find the following :
i. The number of roots in the left half of s-plane.
ii. The number of roots in the right half of s-plane.
iii. The number of roots on the imaginary axis. Use the Routh
Hurwitz Criterion. AKTU 2013-14, Marks 10
Control System 3–7 C (EN-Sem-5)
www.aktutor.in
Answer
Answer
s5 1 2 11
s 4 1 2 10
s3 0 1
s2
s1
s0
2. Putting ( very small positive value) in place of 0 in the first column
Routh array
s5 1 2 11
s4 1 2 10
s3 ( ve) 1 0
2 1
s2 ( ve) 10
2 1
10
s1 ( ve) 0
2 1
s0 10
3. Since there is two sign change in the first column of Routh array,
therefore system is unstable.
PART-2
Root Locus Techniques : Root Contour, Construction of Root Loci,
Effect of Transportation Lag and Root Locus of Non Minimal
Phase System and Effect of Pole-Zero Cancellation.
Questions-Answers
Que 3.8. Explain step by step procedure for plotting root locus.
OR
Write short note on :
i. Centroid
ii. Breakaway points
iii. Steady state error. AKTU 2017-18, Marks 10
Answer
Root locus :
1. Root locus is the plot of loci of the root of the complementary equation
when one or more parameter of open-loop transfer function is varied,
mostly the only one variable available is the gain K. The negative K has
no significance. Hence vary K from 0 to , the plot obtained is called root
locus.
2. It gives the complete dynamic response of system. It provides the
measure of sensitivity of roots of variation in the system.
Rules for construction of root loci :
1. The root locus is symmetrical about the real axis.
2. Each branch of the root locus originates form an open-loop pole where
K = 0 and terminates on an open loop zeros or on infinity where K = .
3. The number of branches of the root locus terminating on infinity is
equal to P – Z.
4. A point on the real axis lies on the locus, if number of open loop poles
and zeros on the real axis to the right of this point is odd.
5. Angles of asymptotes : The (P – Z) root locus branches that approaches
to infinity are along straight line is called asymptote. Asymptotes making
angles with the real axis are given by
(2q 1)180
q = ; where q 0,1,2,....., P Z 1
PZ
6. Centroid : The point of intersection of the asymptotes with the real
axis called centroid is at s = , where
Real part of poles Real part of zeros
=
Number of poles Number of zeros
7. Break points : The break points (breakaway and breakin points) of the
dK
root locus are determined from the roots of the equation = 0.
ds
Stability & Algebraic Criteria 3–10 C (EN-Sem-5)
www.aktutor.in
The r branches of the root locus which meet at a point break away at an
180
angle of ± .
r
8. Angle of departure : The angle of departure from an open loop pole is
given by
d = 180° + ;
where is the net angle contribution at the zero of all other open-loop
poles and zeros.
9. Angle of arrival : The angle of arrival at open loop zero is given by
a = 180° – ,
where is the net angle contribution of all other open loop poles and
zeros.
10. The point of intersection of the root locus branches with the imaginary
axis and the critical value of K can be determined by use of the Routh’s
criterion.
11. Gain : The open-loop gain K in pole-zero from at any point s0 on the root
locus is given by
Product of phasor lengths from s0 to open-loop poles
K=
Product of phasor lengths from s0 to open-loop zeros
Root contour :
1. The root locus technique for the study of stability of closed loop system
from open loop transfer function can be extended by varying parameter
other than K in system from 0 to .
2. Consider a system shown in the Fig. 3.8.1.
R(s) + K C(s)
s (s + )
–
Fig. 3.8.1.
K
3. The open loop transfer function is G(s)H(s) =
s (s )
4. The locus of roots of characteristic equation obtained by varying
parameter of the system other than K from 0 to is called root contour
of the system.
5. The parameters like , K are to be varied simultaneously, while sketching
the root contours.
Steady state error : Refer Q. 2.6, Page 2–11C, Unit-2.
Control System 3–11 C (EN-Sem-5)
www.aktutor.in
Que 3.9. Sketch the root locus for the open loop transfer function
of a unity feedback control system given below and determine value
K
of K at = 0.5, G(s) H(s) = AKTU 2014-15, Marks 10
s(s 1) (s 3)
OR
Construct root loci for open loop transfer function :
K
G(s) H(s) =
s(s 1 ) (s 3)
AKTU 2015-16, Marks 10
Answer
K
Given : G(s) H(s) = , H(s) = 1, = 0.5
s(s 1) (s 3)
To Sketch : Root locus and value of K at = 0.5.
A.
1. The open-loop poles are at s = 0, s = – 1, s = – 3.
2. There is no open loop zeros.
3. Number of poles, P = 3
Number of zeros, Z = 0
P – Z = 3 – 0 = 3 i.e., three branches of root locus end at infinity.
4. Angle of asymptotes :
(2 q 1)180
q = where q = 0, 1, 2, ...(P – Z –1)
PZ
(2 0 1) 180
0 = = 60°
(3 0)
(2 1 1) 180
1 = = 180°
(3 0)
(2 2 1) 180
2 = = 300°
(3 0)
5. Centroid of asymptotes :
Real parts of poles Real parts of zeros
=
PZ
(0 1 3) (0)
= = – 1.33.
30
6. Breakaway points : Between open-loop poles s = 0 and s = – 1, there
exists a breakaway point.
The characteristic equation is
s (s + 1)(s + 3) + K = 0
Stability & Algebraic Criteria 3–12 C (EN-Sem-5)
www.aktutor.in
K = – (s3 + 4s2 + 3s)
dK
= – (3s2 + 8s + 3) = 0
ds
3s2 + 8s + 3 = 0
s = – 1.33 ± 0.88 = – 0.42 and – 2.21
As the breakaway point has to lie between s = 0 and s = – 1, the valid
breakaway point is s = – 0.42.
7. Intersection with j axis :
Characteristic equation, 1 + G(s) H(s) = 0
s (s + 1) (s + 3) + K = 0
s3 + 4s2 + 3s + K = 0
Routh array :
s3 1 3
s2 4 K
s1 (12 K ) / 4
s0 K
The value of K at imaginary axis :
(12 K )
= 0 K = 12
4
Auxiliary equation,
4s2 + 12 = 0
s = ± j 1.73
Im
K=
(j1.732, K = 12)
– 0.35 + j0.6
180° 60°
K= K=0 K=0
–3 –1 0 Re
–1.33
300°
–0.45
(– j1.732, K = 12)
K=
Fig. 3.9.1. Root locus for G(s) = K/s(s + 1)(s + 3).
Control System 3–13 C (EN-Sem-5)
www.aktutor.in
B.
1. Since = 0.5
and cos =
cos = 0.5
= 60°
2. From the origin a line at angle = 60° is drawn as shown in Fig. 3.9.1,
which intersects the root locus plot at s = (– 0.35 + j0.6). As the point
s = (– 0.35 + j0.6) lies on the root locus the following equation is satisfied
K
|G(s) H(s)| = 1 or =1
s( s 1)( s 3)
3. Put s = (– 0.35 + j0.6), we get
K
=1
( 0.35 j 0.6)[( 0.35 j0.6) 1][( 0.35 j 0.6) 3]
K
= 1 or K = 1.66
1.66
Que 3.10. Sketch the root locus for a system having :
K (s + 1)
G(s) = and H(s) = 2
(s + 1) (s + 4s + 5)
and comment on the result. AKTU 2013-14, Marks 10
Answer
1. The characteristics equation of the system is
1 + G(s)H(s) = 0
K s 1
1 =0
s 1 s2 4 s 5
K
1 2 = 0
s 4 s 5
s2 + 4s+ K + 5 = 0
2. Two roots of the above equation are
s 1, s 2 = 2 j K 1
For various value of K the roots are
K 0 3 8
s1 – 2 + j1 – 2 + j2 – 2 + j3 – 2 + j
s2 – 2 – j1 – 2 – j2 – 2 – j3 – 2 – j
K=8 j3
K=3 j2
K=0 j1
Re
–2
K=0 – j1
K=3 – j2
K=8 – j3
– j4
K=
Fig. 3.10.1.
Conclusion :
1. For all values of K from 0 to the roots are complex with negative real
plot indicating that the time response is underdamped.
2. Whatever may be positive values of K there is no chance for the real
part of the roots to become positive, hence the system is inherently
stable.
Que 3.11. Sketch the root locus for the closed loop control system
K
with G(s) = . AKTU 2016-17, Marks 10
s(s 1)( s 2 4 s 5)
Answer
K
Given : G(s) =
s( s 1)(s 2 4 s 5)
To Draw : Root locus.
1. Poles i.e., s = 0, s = – 1, s = – 2 + j and s = – 2 – j
2. There is no open loop zero.
3. Number of poles, P = 4
Number of zeros, Z = 0
P – Z = 4 – 0 = 4 i.e., four branches of root locus terminates at infinity.
4. Angle of asymptotes :
Control System 3–15 C (EN-Sem-5)
www.aktutor.in
(2 q 1)
q = 180
PZ
where, q = 0, 1, 2, ...... upto (P – Z – 1)
201
0 = 180 = 45°
4 0
2 11
1 = 180 = 135°
40
221
2 = 180 = 225°
4 0
231
3 = 180 = 315°
40
5. Centroid of asymptotes :
Real parts of poles – Real parts of zeros
=
PZ
(0 1 2 j 2 j ) (0) 5
= = = – 1.25
40 4
6. Breakaway points : The characteristic equation is 1 + G(s) H(s) = 0
s(s + 1) (s2 + 4s + 5) + K = 0
s4 + 5s3 + 9s2 + 5s + K = 0
dK / ds = 4s3 + (5 × 3s2) + (9 × 2s) + 5 × 1 = 0
4s3 + 15s2 + 18s + 5 = 0
(s + 0.4) (4s2 + 13.4s + 12.5) = 0
(s + 0.4) {(s + 1.675 – j 0.565) (s + 1.675 + j 0.565)} = 0
Therefore three breakaway points are obtained
s = – 0.4 on real axis and
s = – 1.7 + j 0.6 and s = – 1.7 – j 0.6
7. Intersection points with imaginary axis :
The characteristic equation is s4 + 5s3 + 9s2 + 5s + K = 0
The Routh array :
s4 1 9 K
s3 5 5 0
2
s 8 K
5K
1
s 5 0
8
s0 K
Value of K at imaginary axis :
5K
5 = 0; K= 8
8
Solving auxiliary equation formed from the s2 terms in Routh array,
therefore
8s2 + K = 0
Stability & Algebraic Criteria 3–16 C (EN-Sem-5)
www.aktutor.in
8s2 + 8 = 0
s2 + 1 = 0
s=± j
8. Angle of departure from complex pole :
d = 180° – (P – Z)
P1 = 180° – tan–1(1/2) = 153.43°
P2 = 180° – tan–1(1/1) = 135°
P3 = 90°
(–2 + j) = 180° – (P1 + P2 + P3)
(–2 + j) = 180°– (153.43° + 135° + 90°) = – 198.43°
and – j) = + 198.43°
(– 2
Im
–2 + j
j1
13
5° 153.43°
–2 –1 Re
90°
– j1
p = 90 °
3
Fig. 3.11.1.
9. Root locus :
K= K=
Im
–2+j
j1
– 198.43°
–1
Re
+ 198.43° –2
– j1
–2–j
K ( s + 1)
system with OLTF G(s) = . AKTU 2017-18, Marks 10
s 2 (s + 9)
Control System 3–17 C (EN-Sem-5)
www.aktutor.in
Answer
K ( s 1)
Given : G(s) = , H(s) = 1
s 2 ( s 9)
To Draw : The root locus.
1. The open loop poles i.e., s = 0, s = 0, s = – 9
2. Open loop zero are at s = – 1
3. Number of poles, P = 3
Number of zeros, Z = 1
P – Z = 3 – 1 = 2 i.e., two branches of root locus terminates at infinity.
4. Angle of asymptotes :
(2 q 1)
q = 180
PZ
where, q = 0, 1, 2, ...... upto (P – Z – 1)
(2 0 1)
0 = 180 = 90°
3–1
2 1 1
1 = 180 = 270°
3–1
5. Centroid of asymptotes :
Real parts of poles – Real parts of zeros
=
PZ
0 9 1 8
= =–4
2 2
6. Breakaway points :
1 + G(s)H(s) = 0
K ( s 1)
1 2 1 = 0
s ( s 9)
s2(s + 9) + K(s + 1) = 0
s2 (s 9)
K=
(s 1)
dK (s 1)(3 s 2 18 s) s2 (s 9)(1) 2s( s 3)2
= = =0
ds (s 1)2 ( s 1)2
– 2s(s + 3)2 = 0
s = 0, s = – 3
7. Intersection with imaginary axis :
s3 + 9s2 + Ks + K = 0
Routh array :
s3 1 K
s2 9 K
s1 8 K / 9
s0 K
Stability & Algebraic Criteria 3–18 C (EN-Sem-5)
www.aktutor.in
Value of K at imaginary axis,
8K
=0
9
K=0
Now, 9s2 + K = 0
9s2 = 0
s=0
K=
Im
K=0
–9 –4 –1 Re
K=
Fig. 3.12.1. Root locus.
Que 3.13. Sketch the RL for a unity feedback system with OLTF
K (s2 + 2s + 10)
G(s) = . AKTU 2017-18, Marks 10
s2 + 4 s + 5
Answer
K ( s2 2s 10)
Given : G(s) = , H(s) = 1
s2 4 s 5
To Draw : Root locus.
1. The open loops poles are at s = – 2 – j, s = – 2 + j
2. The open loop zeros are at s = – 1 – j3, s = – 1 + j3
3. Number of poles = 2
Number of zeros = 2
P – Z = 2 – 2 = 0 i.e., none of branches of root locus terminates at
infinity.
4. Angle of departure and arrival can’t be calculated because the centroid,
= .
Control System 3–19 C (EN-Sem-5)
www.aktutor.in
Im
j3
j1
Re
–2
– j1
– j3
Fig. 3.13.1.
K ( s 4)
transfer function G(s) =
( s + 2)(s 1)
Answer
K (s 4)
Given : G(s) =
(s 2)( s 1)
To Draw : Root locus.
1. The open loop poles are at s = – 2, s = 1
2. The open loop zeros at s = 4
3. Number of poles = 2
Number of zeros = 1
P – Z = 1 i.e., one branch of root locus terminates at infinity.
4. Centroid of asymptotes :
Real parts of poles – Real parts of zeros 1 2 4
= = =–5
PZ 1
5. Breakaway point :
1 + G(s)H(s) = 0
1 K ( s 4)
=0
( s 2)( s 1)
2
s + s – 2 + K(s – 4) = 0
2 s s2
K=
(s 4)
Stability & Algebraic Criteria 3–20 C (EN-Sem-5)
www.aktutor.in
dK (s 4)(1 2s) (2 s s 2 )(1)
=
ds (s 4)2
s2 8 s 2
= =0
(s 4)2
s = 8.24 (Breakin point)
and s = – 0.24 (Breakaway point)
6. Angle of asymptotes :
(2q 1) 1
q = 180 × 180° = 180°
PZ 21
Im
8.24
–2 1 4 Re
– 0.24
Answer
Pole-zero cancellation :
1. Let us assume that some of the poles of G are at undesirable locations
and that we want to place them more appropriately.
2. We can insert a controller in series before the system in order to cancel
out the undesirable poles and keep only those which are suitable.
3. The undesirable poles of G are not removed by this technique, but that
their effect is cancelled out by the coincident zeros of the controller.
4. This is known as pole-zero cancellation. A difficulty with pole zero
cancellation is that it may not be exact.
Effects of pole-zero cancellation are :
1. If the specification on allowable steady state errors cannot be met, a low
frequency pole can be cancelled and replaced with a lower frequency
pole, yielding a large forward loop.
2. If poles with small damping ratios are present in the plant transfer
function, they may be cancelled and replaced with poles which have
larger damping rates.
Control System 3–21 C (EN-Sem-5)
www.aktutor.in
Que 3.16. Discuss the effects of adding poles to root locus.
Answer
Effects of addition of open loop poles :
1. Root locus shifts towards imaginary axis.
2. System stability relatively decreases.
3. System becomes more oscillatory in nature.
4. Range of operating value of K for stability of the system decreases.
Example :
1. Consider, G(s)H(s) = K / s(s + 2)
Corresponding root locus is shown in the Fig. 3.16.1.
K=
Im
Breakaway point
K=0 K=0
Re
–2 –1 0
K=
Fig. 3.16.1.
2. Now if pole at s = – 4 added to G(s)H(s) root locus becomes as shown in
the Fig. 3.16.2.
Im
K=
Answer
Effects of addition of zeros :
1. Root locus shifts to left, away from imaginary axis.
Stability & Algebraic Criteria 3–22 C (EN-Sem-5)
www.aktutor.in
2. Relative stability of the system increases.
3. System becomes less oscillatory.
4. Range of operating values of K for system stability increases.
Example :
K
1. Root loci for G(s) H(s) = is shown in the Fig. 3.17.1.
s(s 2)
K=
Im
Breakaway point
K=0 K=0
Re
–2 –1 0
K=
Fig. 3.17.1.
Re
–4 –2 0
Fig. 3.17.2.
3. It can be seen that root locus shift towards left. So as roots move towards
left half of s-plane, hence relative stability increases.
Que 3.18. Define transportation lag. Explain its effect with suitable
example.
Answer
1. The transportation lag is the delay between the time an input signal is
applied to a system and the time the system reacts to that input signal.
Control System 3–23 C (EN-Sem-5)
www.aktutor.in
2. Transportation lags are common in industrial applications. They are
often called dead time.
3. Consider a system having open-loop transfer function
Ke– sT
G(s)H(s) =
s(s 2)
where T is the transportation delay in seconds and is given as 1 s.
4. Let us draw the root loci for K varying in the range 0 < K < .
5. If the transportation delay is small, then we can assume
e–sT = 1 – sT
6. Rewriting G(s)H(s) as
K (1 – s)
G(s)H(s) = ; T=1
s( s 2)
– K ( s – 1)
G(s)H(s) =
s( s 2)
7. The characteristic equation becomes
1 – G(s)H(s) = 0
The root loci are drawn as shown in Fig. 3.18.1.
Im
j 2(K = 2)
– j 2(K = 2)
Control System
www.aktutor.in 4–1 C (EN-Sem-5)
4 Frequency
Response Analysis
PART-1
Frequency Response Analysis from Transfer Function Model,
Correlation between Time and Frequency Response.
Questions-Answers
Answer
A. Frequency response :
1. The magnitude and phase relationship between the sinusoidal input
and the steady state output of a system is termed as the frequency
response.
2. Linear system with sinusoidal input
r(t) = A sin t
3. Steady state output of
c(t) = B sin (t + )
B. Various methods of frequency response :
1. Polar plots 2. Inverse polar plots
3. Nyquist stability criterion 4. Bode plot
C. Advantages of frequency domain approach :
1. Without the knowledge of the transfer function, the frequency response
of stable open loop system can be obtained experimentally.
Control System
www.aktutor.in 4–3 C (EN-Sem-5)
2. These methods are easy to use for design of control and for finding
absolute as well as relative stability of the system.
3. Even if the system has moderate non-linearity it can be depicted by an
approximate transfer function.
4. There is a close relation between frequency response of a system and its
step response.
D. Disadvantages of frequency response :
1. For systems with very large time constants, the frequency test is
cumbersome to perform as the time required for the output to reach the
steady-state for each frequency of the test signal is excessively long.
2. The frequency response test obviously cannot be performed on
non-interruptible systems.
E. Limitations of frequency response methods :
1. These methods can be applied only to linear systems.
2. For an existing system, obtaining frequency response is possible only if
the time constants are up to few minutes.
3. Obtaining frequency response practically is fairly time consuming.
Answer
1. Resonant peak (Mr) : Resonance peak is the peak value of magnitude
of closed loop frequency response.
1
For second order, Mr =
2 1 2
2. Resonance frequency : The frequency at which resonant peak occurs
is called resonant frequency.
1
3 dB
0.707
r c
0
B.W
Fig. 4.2.1.
Frequency Response Analysis
www.aktutor.in 4–4 C (EN-Sem-5)
4. Cut-off rate : The cut-off rate is the frequency rate at which the
magnitude ratio decreased beyond the cut-off frequency c.
Answer
Closed loop frequency response :
1. Consider the transfer function for closed loop system,
C ( s) G ( s)
=
R( s) 1 G( s) H (s)
2. For unity feedback, H(s) = 1
C ( s) G (s )
= ...(4.3.1)
R( s) 1 G ( s)
Put s = j
C( j) G ( j)
= ...(4.3.2)
R( j) 1 G( j)
3. The polar plot of eq. (4.3.2) is shown in Fig. 4.3.1.
Im
– 1 + j0 0
Re
A r
(r – )
B
Fig. 4.3.1.
4. From Fig. 4.3.1, OB = G(j)
OA = – 1
AB = OB – OA = G(j) – (– 1)
AB = 1 + G(j)
5. From eq. (4.3.2)
C ( j) OB
= M() =
R( j)
AB
Control System
www.aktutor.in 4–5 C (EN-Sem-5)
C( j) OB r
=
= (r – )
R( j) AB
C( j)
= M() e j ()
R( j)
where M(j) is the magnitude and () = r – .
6. Frequency response consists of two parts :
i. Magnitude,
ii. Phase angle.
iii. Both can be plotted against different values of .
7. Hence frequency response of closed loop system is plot of magnitude
and phase angle.
Correlation :
1. For 2nd order system, the transfer function is,
C ( s) 2n
= 2
R( s) s 2 n s 2n
where, = Damping factor
n = Natural frequency of oscillations
+ n2
R(s) C(s)
2
– s + 2 ns
Fig. 4.3.2.
2. Closed loop frequency response is,
C( j) 2n
= T(j) =
R( j) ( j) 2 n ( j) n2
2
2n 1
= = 2
2 n ( j) n2
2
1 2 j
n n
1
= ...(4.3.3)
(1 u2 ) j 2u
where u = /n, normalized driving frequency.
1
|T(j)| = M = ...(4.3.4)
(1 u ) (2u)2
2 2
2u
and T(j) = = tan 1 ...(4.3.5)
1 u2
3. The steady state output is
1 2u
c(t) = sin t tan 1
2 2
(1 u ) (2u)
2 1 u2
Frequency Response Analysis
www.aktutor.in 4–6 C (EN-Sem-5)
ur = 1 2 2
r = n 1 2 2 ...(4.3.6)
1
Mr = ...(4.3.7)
2 1 2
5. The phase angle of T(j) at resonant frequency is
1 2 2
r = tan 1
M( = 1 < 0.707)
1.0
Magnitude
M
Mr
M( = 0.707)
0
u r = r / n
Normalized frequency, u
Fig. 4.3.3. Frequency response magnitude characteristics.
Control System
www.aktutor.in 4–7 C (EN-Sem-5)
0°
( = 1)
– 90°
Phase
angle
– 180°
0 1.0
Normalized frequency u
Fig. 4.3.4. Frequency response phase characteristic.
Bandwidth :
1
1. The frequency at which M has value of and is called cut-off frequency
2
c. The signal frequencies above cut-off are attenuated.
1
2. The range of frequencies for which M is known as bandwidth b.
2
The low pass filters has bandwidth equal to cut-off. b indicates the
noise filtering characteristics of the system.
M
Mr
1
0.707
Bandwidth
r c
Fig. 4.3.5.
3. Normalized bandwidth, ub = b/n
1 1
M=
2 2 2
(1 u ) (2ub )
b
2
1
ub = [1 2 2 2 4 2 4 4 ]2
1
b = [1 2 2 2 4 2 4 4 ]2
4. The damped frequency of oscillations d and peak overshoot Mp of the
step response for 0 1 are
Frequency Response Analysis
www.aktutor.in 4–8 C (EN-Sem-5)
d = n 1 2
Mp = exp
1 2
1
5. For > , Mr does not exists so the correlation breaks down
2
r 1 2 2
=
d 1 2
Que 4.4. The steady state output of the system for a sinusoidal
input of unit magnitude and variable frequency is given as
1 27u
c(t) = sin t tan 1
2 2 2 2
(1 u ) 4 z u 1 u2
Determine :
i. Resonant frequency ii. Resonant peak
iii. Bandwidth iv. Phase angle.
AKTU 2016-17, Marks 10
Answer
1 27u
Given : c(t) = sin t tan 1
2 2 2 2
(1 u ) 4 z u 1 u2
To Find : i. Resonant frequency
ii. Resonant peak
iii. Bandwidth.
1
1. M=
(1 u ) 4 z2u2
2 2
27u
2. = tan 1
1 u2
3. The frequency where M has a peak value is known as resonant
frequency. At this frequency, the slope of the magnitude curve is
r
zero. Let r be the resonant frequency and ur =
n
dM 1 [ 4(1 ur2 )ur 8 z2ur ]
=– =0
du u ur 2 [(1 4 2r )2 4 z2ur2 ]3/ 2
– 4(1 – ur2) = 8z2
1 – ur2 = 2z2
ur = 1 2z2
r = n 1 2 z2
Control System
www.aktutor.in 4–9 C (EN-Sem-5)
1 1 1
= = =
4 z4 4 z2 (1 2 z2 ) 2 z z2 1 2 z2 2 z 1 z2
1
5. The range of frequencies over which M is equal to or greater than :
2
b
Putting ub =
n
1 1
M=
(1 ub2 )2 4 z2ub2 2
(1 – ub2)2 + 4z2ub2 = 2
1 + ub4 – 2ub2 + 4z2ub2 – 2 = 0
ub = [1 2 z2 2 4 z2 4 z4 ]1/ 2
b = n [1 2 z2 1 4 z2 4 z4 ]1/ 2
PART-2
Construction of Polar and Inverse Polar Plots,
Nyquist Stability Criterion.
Questions-Answers
Que 4.5. Explain Polar plot. Discuss phase margin and gain
margin on polar plot.
OR
Define the following :
i. Resonance frequency ii. Bandwidth
iii. Cut-off rate iv. Phase margin
v. Gain margin
Also discuss the advantages of frequency domain analysis.
AKTU 2013-14, Marks 10
Answer
Resonance frequency : Refer Q. 4.2, Page 4–3C, Unit-4.
Bandwidth : Refer Q. 4.2, Page 4–3C, Unit-4.
Cut-off rate : Refer Q. 4.2, Page 4–3C, Unit-4.
A. Polar plot :
1. Polar plot of a transfer function G(j) is a plot of the magnitude of G(j)
versus the phase angle of G(j) on polar coordinates as is varied from
zero to infinity. Therefore it is the locus of vectors |G(j)|G(j) as
is varied from zero to infinity.
2. In frequency response we have
M = |G(j) H(j)| = Magnitude
= G(j) H(j) = Phase
Values of M and are obtained by varying the input frequency from 0
to .
Table 4.5.1.
M = |G(j) H(j)| = G(j) H(j)
0 M0 0
1 M1 1
2 M2 2
.. .. ..
.. .. ..
... ... ...
M
Control System
www.aktutor.in 4–11 C (EN-Sem-5)
3. If the polar plot passes through the point (–1 + j0) the system is
marginally stable and if the critical point is outside the polar plot, the
system is stable.
Im Im
Positive
gain Negative
margin phase
1/Kg margin
m
–1 m Re –1 Re
Positive
phase
margin
3
i. Type 0 system, G(s) =
( s + 1)( s + 2)
12
ii. Type 1 system, G(s) =
s( s + 1)( s + 2)
Answer
3 12
Given : i. G(s) = ii. G(s) =
( s 1)( s 2) s(s 1)(s 2)
To Sketch : Polar plot.
Control System
www.aktutor.in 4–13 C (EN-Sem-5)
i.
3
1. G(s) =
( s 1)( s 2)
3
2. Putting s = j, G(s) =
(1 j)(2 j)
3
3. M= and = – tan– 1 – tan– 1
1 4 2 2 2
For various value of , the value of M and are :
Table 4.6.1.
M
0 1.5 0°
5 0.109 – 146.89
10 0.029 – 178.28
0 – 180°
= =0
– 180° 0°
1.5
Fig. 4.6.1.
ii.
12
1. G(s) =
s(s 1)(s 2)
12
2. Putting s = j, G(j) =
j(1 j)(2 j)
12
3. M= and = – 90° – tan– 1 – tan–1
1 2 4 2 2
Table 4.6.2.
M
0 – 90°
5 0.087 – 236.8°
10 0.012 – 268.2°
0 – 270°
Frequency Response Analysis
www.aktutor.in 4–14 C (EN-Sem-5)
– 270°
Q =
– 180°
=0 – 90°
Fig. 4.6.2.
The co-ordinates of Q can be obtained by rationalizing G(j).
10
G(s) =
( s + 10)
Answer
10
Given : G(s) =
(s 10)
To Sketch : Inverse polar plot.
10
1. G(s) =
(s 10)
10 1
2. Putting s = j, G(j) =
( j 10) 1 j(0.1)
1
3. = 1 + j(0.1) = M
G( j)
where, M = –1
1 (0.1)2 and = tan (0.1)
Table 4.7.1.
M
0 1 0°
1 1.005 5.71°
10 1.414 45°
100 10.05 84.28°
. . .
. . .
. . .
. . .
+ 90°
Control System
www.aktutor.in 4–15 C (EN-Sem-5)
0° 1
=0
Fig. 4.7.1. Inverse polar plot.
10 e s
i. G(s) =
s +1
32
ii. G(s) =
( s + 4)( s2 + 4 s + 8)
and find its points of intersection with the real and imaginary
axes. AKTU 2017-18, Marks 10
Answer
10 e s
i. Given : G(s) =
s1
To Sketch : Polar plot.
1. Putting s = j
10 e j
G(j) =
j 1
2. Magnitude :
|10 e j | 10
|G(j)| =
2 1 2 1
3. Phase angle :
G(j) = – – tan – 1
4. Here the magnitude decreases from unity monotonically and the phase
angle also decreases monotonically and indefinitely, the polar plot of the
given transfer function is spiral as shown in Fig. 4.8.1.
Frequency Response Analysis
www.aktutor.in 4–16 C (EN-Sem-5)
Im
O 1
Re
180°
– 90°
Fig. 4.8.1.
32
ii. Given : G(s) =
(s 4)(s 2 4 s 8)
To Sketch : Polar plot.
1. Putting s = j
32
G(j) =
( j 4)[ j 2 4 j 8]
32
=
( j 4)[(8 2 ) 4 j ]
32(4 j)[(8 2 ) 4 j]
=
(16 2 )[(8 2 )2 162 ]
32[4(8 2 ) 42 j{(8 w2 ) 16}]
=
(16 2 )[(8 2 ) 2 162 ]
2. Saparating into real and imaginary parts
32[4(8 2 ) 4 2 ] j32[ (8 2 ) 16]
G(j) =
(16 )[(8 ) 16 ] (16 2 )[(8 2 )2 162 ]
2 2 2 2
4
G(j) = tan 1 tan 1
4 8 2
6. The value of G(j) at = 0, and =
32
|G(j0)| = =1
48
|G(j)| = 0
G(j0) = 0°
G(j) = – 270,
Im
=+2
=+ = 4.899
Re
=–
– 4.899
=0
=–2
Fig. 4.8.2.
Answer
1. Nyquist stability criterion calculates closed loop stability with the help of
open loop transfer function without any complex calculations.
2. Complex calculations are done in polar plot.
3. If open loop transfer function is G(s)H(s) then closed loop poles
[1 + G(s)H(s) = 0].
Im j
Re
– 1 + j0
Fig. 4.9.1.
Encirclement :
1. If – 1 + j0 point is encircled and in anti-clockwise direction.
N= P–Z
N = –ve (for clockwise)
Frequency Response Analysis
www.aktutor.in 4–18 C (EN-Sem-5)
Fig. 4.9.2.
We check Z,
If Z > 0 the given system is unstable.
If Z = 0 the given system is stable.
Phase margin and gain margin in polar plot : Refer Q. 4.5,
Page 4–10C, Unit-4.
Que 4.10. For the G(s) = 1/s(s – 2), H(s) = 1. Sketch the Nyquist plot
Answer
1
Given : G(s) = , H(s) = 1
s(s 2)
To Sketch : Nyquist plot.
1
1. G(s) =
s(s 2)
1
G(s) =
s(2 s)
2. Put s = j
1
G(j) =
j(2 j)
1
|G(j)| =
4 2
3. G(j) = 180° – 90° + tan–1
2
= 90° + tan–1
2
Control System
www.aktutor.in 4–19 C (EN-Sem-5)
Table 4.10.1.
M
0 90°
1 0.4472 116.56°
10 9.8 × 10–3 168.69°
100 9.99 × 10 – 5 178.85°
0 180°
Amplitude decreases Angle increases
4. This system is unstable because contour has infinite radius, so it encircle
(– 1, 0) point.
N = P – Z; N = – 1
P = 1 (number of open loop pole in right half of s-plane)
Z= 2
0
= Im
=
(– 1, 0) = – Re
=–0
Fig. 4.10.1.
Que 4.11. Draw the Nyquist plot for the open loop transfer function
given below and comment on closed loop stability
1.5 ( s + 4)
G(s)H(s) = . AKTU 2014-15, Marks 10
s( s – 2)
Answer
1.5( s 4)
Given : G(s)H(s) =
s( s 2)
To Draw : Nyquist plot.
1.5 ( j 4) – 1.5 ( j 4)
1. Put s = j, G(j)H(j) = =
j ( j – 2) j (2 – j)
1.5 16 2
2. M = |G(j) H(j)| =
4 2
www.aktutor.in
Frequency Response Analysis 4–20 C (EN-Sem-5)
3. = G(j) H(j) = 90 tan –1 tan –1
4 2
Table 4.11.1.
M
0 90°
0.1 29.97 94.29°
0.5 5.866 111.16°
10 0.1584 236.89°
100 0.015 266.57
0 270°
+0
s-Plane
–
–1 0.74 + Re
–0
Fig. 4.11.1.
Que 4.12. Sketch the Nyquist plot for the system with open loop
transfer function
60
G(s) H(s) =
(s 1 ) (s 2) (s 5)
and comment on stability. AKTU 2015-16, Marks 10
Control System
www.aktutor.in 4–21 C (EN-Sem-5)
Answer
60
Given : G(s) H(s) =
(s 1 ) (s 2) (s 5)
To Sketch : Nyquist plot.
60
1. Put s = j, G(j) H (j) =
( j 1 ) ( j 2) ( j 5)
60
2. Magnitude, M = |G (j) H (j)| =
( 1)( 4)( 2 25)
2 2
3. Phase = G(j)H(j)
1 1 1
= tan () tan tan
2 5
Table 4.12.1.
M
0 6 0°
1.12 3.402 – 90.11°
4.12 0.4769 – 180°
10 0.0523 133.58
0 90° (– 270°)
4. N=P–Z
P=0 (No pole on RH plane)
N=0 (No encirclement of – 1 + j0 point)
Z=0 (Number of roots of closed loop characteristic
equation having positive real parts).
So the given system is stable system.
Im
–0
Re
– 1 + j0 –
=+0
0.4769
Fig. 4.12.1.
Frequency Response Analysis
www.aktutor.in 4–22 C (EN-Sem-5)
Answer
(s – z1 )
Given : G(s) =
s(s p1 )
To Check : Stability of system.
j z1 2 z12
1. G(j) H(j) =
j( j p1 ) 2 p12
1 ( p1 z1 )
= 90 tan 2
( p1 z1 )
2. Now,
lim G ( j) H ( j ) = + 90°
0
1
G( j p1 z1 ) H ( j p1 z1 ) = 0
p1 p1 z1
3. Thus the locus comes down in the first quadrant, crosses the positive
real axis into the fourth quadrant, and approaches the origin from an
angle of – 90°.
4. Path def maps into the origin, and ija maps into on semicircle at
infinity. The resulting plot is shown in Fig. 4.13.1.
Im GH
a
GH(j )
= ± p1z1
j d , e, f Re GH
i
Fig. 4.13.1.
Control System
www.aktutor.in 4–23 C (EN-Sem-5)
Ke sT
G(s) H(s) = , determine the maximum value of the gain ‘K’ for
s( s 1)
stability. AKTU 2016-17, Marks 10
Answer
Ke sT
Given : G(s) H(s) =
s(s 1)
To Find : Maximum value of K.
1. Putting s = j
Ke jT
G(j) H(j) =
j( j 1)
K
2. |G(j) H(j)|=
1 2
180
G(j) H(j) = – tan–1 – 90° – T ×
= – 57.3T – 90° – tan–1
3. Intersection of Nyquist plot with negative axis of G(s) H(s) plane is
determined by using following relation,
G(j) H(j) = –180° (2k + 1)
where, k = 0, 1, 2 ....
4. It the first instant, the Nyquist plot intersects the negative real axis of
G(s) H(s) plane for k = 0, and the frequency at intersection point is 2.
5. Therefore, G(j2) H(j2) = – 180°
– 57.32T – 90° – tan–12 = – 180°
– 57.32T – tan–12 = – 90°
6. T = 0.5
– 28.652 – tan–1 2 = 90°
7. Using trial and error, 2 = 1.3075 rad/sec
K K
G(j2) H(j2) = = = 0.4646 K
2 1 22 1.3075 1 (1.3075)2
8. For stability the point (–1 + j0) be placed outside the Nyquist plot
0.4646 K < 1
K < 2.152
www.aktutor.in
Frequency Response Analysis 4–24 C (EN-Sem-5)
PART-3
Determination of Gain and Phase Margin from
Bode and Nyquist Plots, Nichol Charts.
Questions-Answers
Que 4.15. Write procedure to draw Bode plot from given open loop
transfer function.
Answer
A. Bode plot : In this plot we have two plots namely Magnitude plot;
which is drawn in between magnitude of transfer function in dB and
log10 ( is frequency) and phase plot; which is drawn in between
phase angle of transfer function in degree and log10 generally.
B. Procedure to draw Bode plot :
Step 1 :
Replace s = j from open loop transfer function G(s) to have G(j).
Step 2 :
Change all the factor of G(j) in standard form and make a table with
remarks for both magnitude plot and phase plot.
Some standard results and remarks is summarized here :
Control System
www.aktutor.in 4–25 C (EN-Sem-5)
Table 4.15.1.
S. Factors of Corner Remarks for Remarks for
No. G(j ) Frequencies Magnitude Plot Phase Plot
Step 3 :
From the table, draw the magnitude plot for each of the factor by
calculating resultant slope and starting with lowest corner frequency
and ending with the highest one.
Step 4 :
Add all the expression of phase angle to have different values of
(resultant phase).
Step 5 :
Draw with log10 to have phase plot.
Que 4.16. Explain gain margin and phase margin and how to
determine them using Bode plot.
Answer
A. Gain Margin (GM) :
1. It is the margin in gain which is allowed till the system reaches on
the verge of instability.
2. Mathematically, gain margin is reciprocal of magnitude of G(j)H(j)
at phase crossover frequency.
GM = 1 /|G ( j) H ( j)
pc
G( j) H ( j) can be directly known on phase angle plot.
gc
Que 4.17. Draw Bode plot (log magnitude plot) for the transfer
function.
20 s
G(s) =
s2 20 s (100)2
AKTU 2016-17, Marks 10
Answer
K = 20 / (100)2 = 0.002
1 2 3 4 5 6 7 8 9 101 2 3 4 5 6 7 8 9 102 2 3 4 5 6 7 8 9 103
Gain
Control System
in – 80°
40 dB
dB
– 20 dB/dec – 100°
Bode plot :
20 dB/dec
– 120°
– 140°
10 dB
– 160°
0 dB – 180°
– 200°
g = 1 rad/sec g = 100 rad/sec
– 10 dB
– 220°
– 240°
– 20 dB – 260°
www.aktutor.in
– 280°
Que 4.18. Draw the bode plot for the transfer function given as,
apply correction to the magnitude plot for the quadratic term and
comment on stability.
5
G(s) H(s) =
s (1 0.1s) ( s2 0.4 s 1)
2
Answer
5
Given : G(s) H(s) =
s 2 (1 0.1s)(s 2 0.4 s 1)
To Sketch : Bode plot.
1. G(s)H(s) is in time constant form.
2. Factors :
i. K = 5, 20 log K = 13.98 dB, straight line parallel to Log axis.
1
ii. , 2 poles at origin, straight line of slope – 40 dB / dec passing through
s2
intersection point of = 1 and 0 dB.
1
iii. Quadratic pole, , 2n = 1 i.e., c1 = n = 1
1 0.4s s2
Straight line of slope – 40 dB / dec for 1.
iv. Comparing middle term with 2 n = 0.4, = 0.2
Correction = – 20 Log 2 = + 7.95 dB at = n = 1.
Thus magnitude plot will shift upwards by + 7.95 dB at = n = 1.
1 1
iv. Simple pole, , T2 = 0.1, C2 = = 10, straight line of slope
1 0.1s T2
– 20 dB/dec for 10.
Table 4.18.1.
Range of 0<<1 1 < 10 10 <
Resultant slope – 40 – 40 – 40 = – 80 – 80 – 20 = – 100
in dB/dec
3. Phase angle :
5
G(j)H(j) =
( j)2 (1 0.1 j) [(1 – 2 ) 0.4 j]
Control System
www.aktutor.in 4–29 C (EN-Sem-5)
Bode plot :
1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
100
80
dB
0
–4
60
Due to correction
for quadratic factor
40
K = 5, 14 dB
20
–20
– 80 dB/dec
in deg
–90°
– 100 dB/dec
–180°
P.M = – 135°
–270°
–360°
–450°
gc = 1.9
0.1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
Frequency Response Analysis
www.aktutor.in 4–30 C (EN-Sem-5)
Table 4.18.2.
1 0.4
– tan –1 2
– tan–1 0.1
( j)2 1
0.1 – 180° – 2.31° – 0.57° – 182.88°
1 – 180° – 90° – 5.71° – 275.71°
5 – 180° + 4.76° – 180° = –175.23° – 26.56° – 381.79°
– 180° – 180° – 90° – 450°
4. Comment : From the Bode plot pc 0 and phase angle plot does not
cross – 180°. Thus system is unstable with G.M. = – dB and P.M. = –
135° at gc = 1.9 rad/s.
Que 4.19. A unity feedback control systems has :
40
G(s) =
s(s + 2) (s + 5)
Draw the Bode Plot. Find Gain Margin. AKTU 2013-14, Marks 10
Answer
40
Given : G(s) =
s(s 2) (s 5)
To Draw : Bode plot.
1. Magnitude plot :
40 40 4
G(s) = = =
s(s 2) (s 5) s s s s
s 2 1 5 1 s 1 1
2 5 2 5
First corner frequency (pole) (1) = 2 rad/s
Second corner frequency (pole) (2) = 5 rad/s
2. Phase plot :
4
G(s) =
s s
s 1 1
2 5
Put s = j
4
G(j) =
j(1 0.5 j) (1 0.2 j)
= – 90° – tan – 1 0.5 – tan – 1 0.2
Control System
www.aktutor.in 4–31 C (EN-Sem-5)
Magnitude plot :
Table 4.19.1.
S. No. Factor Corner Asymptotic log
frequency magnitude
characteristics
Straight line of constant
1. 1/s None slope (– 20 dB/dec) passing
through at = 1
Straight line of constant
2. 1/(1 + 0.5s) 1 = 2 slope (– 20 dB/dec)
originating from 1 = 2
Straight line of constant
3. 1/(1 + 0.2s) 2 = 5 slope (– 20 dB/dec)
originating from 2 = 5
Straight line of constant
4. 4 None slope of 0 dB/dec starting
from 20 log 4 = 12 dB point
Phase plot :
Table 4.19.2
S. No. (rad/sec) (degrees)
1. 0 – 90°
2. 0.2 – 98°
3. 1 – 127.874°
4. 2 – 156.8°
5. 5 – 203.198°
6. 10 – 232.125°
7. 15 – 243.97°
Frequency Response Analysis
www.aktutor.in 4–32 C (EN-Sem-5)
Bode plot :
1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
100
80
– 60 dB/dec
60
– 40 dB/dec
40
– 20 dB/dec
20 c1 = 2 rad/sec
12 dB
0
– 40 dB/dec
– 20 dB/dec – 60 dB/dec
c2 = 5 rad/sec
in deg
–80°
–100°
–120°
–140°
–160°
–180°
–200°
–220°
–240°
0.1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
Result :
i. Gain crossover frequency = 2.8 rad/s
ii. Phase cross over frequency = 3.3 rad/s
iii. Gain margin = 3 dB
iv. Phase margin = 8°.
Control System
www.aktutor.in 4–33 C (EN-Sem-5)
Que 4.20. For a unity feedback system, the open loop transfer
function is
2(s 0.25)
G(s) H(s) =
s2 (s 1) (s 0.5)
Draw bode plot and determine gain margin, phase margin.
AKTU 2015-16, Marks 10
Answer
2(s 0.25)
Given : G(s) H(s) =
s 2 ( s 1) ( s 0.5)
To Draw : Bode plot.
2 0.25 s
1
0.5 0.25 (4 s 1)
1. G(s) H(s) = = 2
2 s s (s 1) (2s 1)
s ( s 1) 1
0.5
2. K = 1, 20 Log K = 0 dB, No effect on Bode Plot.
1
3. , 2 Poles at origin, straight line of slope – 40 dB/dec
s2
1
4. (1 + 4s), simples zero, T1 = 4, C1 = = 0.25,
T1
straight line of slope + 20 dB/dec for 0.25
1 1
5. , simple zero, T2 = 2, C2 = = 0.5,
1 2s T2
straight line of slope – 20 dB/sec for 0.5
1 1
6. , simple zero, T3 = 1, C3 = = 1,
1 s T3
straight line of slope – 20 dB/sec for 1
Table 4.20.1.
Range of 0 < < 0.25 0.25 < 0.5 0.5 < 1 1 <
Resultant slope – 40 – 40 + 20 – 20 – 20 – 40 – 20
in dB/dec = – 20 = – 40 = – 60
(1 4 j)
7. Phase angle table : G(j)H(j) =
( j)2 (1 2 j) (1 j)
Frequency Response Analysis
www.aktutor.in 4–34 C (EN-Sem-5)
Bode plot :
1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
100
80
– 60 dB/dec
60
– 40 dB/dec
40
– 20 dB/dec
20
G.M. = – 21 dB
– 20 dB/dec – 60 dB/dec
in deg
–150°
–180°
P.M. = – 36°
–210°
–240°
–270°
0.1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
Control System
www.aktutor.in 4–35 C (EN-Sem-5)
Table 4.20.2.
1
+ tan–14 + tan–12 –tan –1
( j)2
Answer
A. Numerical : Refer Q. 4.13, Page 4–22C, Unit-4.
B. M-circle :
1. To obtain the constant magnitude loci (M-circle), let us first note that
G(j) is a complex quantity and can be written as follows :
G(j) = X + jY
where X and Y are real quantities
2. Then M is given by
X jY
M=
1 X jY
and M2 is
X2 Y2
M2 =
(1 X )2 Y 2
3. Hence, X2 (1 – M2) – 2M2X – M2 + (1 – M2) Y2 = 0 ...(4.21.1)
If M = 1, then from above eq. (4.21.1) we obtain X = –1/2. This is the
equation of a straight line parallel to the Y axis and passing through
point (–1/2, 0).
4. If M 1, equation can be written as
Frequency Response Analysis
www.aktutor.in 4–36 C (EN-Sem-5)
2M 2 M2
X2 2
X 2 Y2 = 0
M 1 M 1
5. If the term M2 /(M2 – 1)2 is added to both sides of this last equation we
obtain
2
2 M2 2 M2
X 2 Y = ...(4.21.2)
M 1 ( M 2 1)2
M 2
6. Eq. (4.21.2) is the equation of circle with centre X = , Y = 0 and
( M 2 1)
with radius |M/(M2 – 1)|.
7. The constant M loci on the G(s) plane are a family of circle. The centre
and radius of the circle for a given value of M can be easily calculated.
C. Significance : M-circle can be used to obtain closed loop frequency
response from open loop frequency response.
M = 1.6 Y
M=1
M = 0.5
M=3
Que 4.22. Discuss constant phase angle loci for frequency response
system.
OR
Discuss N-circle construction.
Answer
1. N-circles having constant phase.
G(j)= x + jy where x and y are real number.
C ( s) G( s) x jy
2. Also = ...(4.22.1)
R( s) 1 G( s) H ( s) 1 x jy
3. If is the phase shift of the output with respect to input,
y y
we have = tan 1 tan 1
x 1 x
Control System
www.aktutor.in 4–37 C (EN-Sem-5)
y y
1 y 1 y x 1 x
tan = tan tan tan
x 1 x y y
1
x 1 x
y
tan =
x2 x y2
4. If tan = N, we have
1
x2 + x + y2 – y=0 ...(4.22.2)
N
This can be modified to
2 2 2
1 1 1 1
x y = ...(4.22.3)
2 2 N 4 2N
1 1
5. This is again an equation of circle having centre , and radius
2 2N
1 1
.
4 4N 2
As x = 0, y = 0 and
x = –1, y = 0 satisfy this eq. (4.22.3) irrespective value of N. Each circle
passes through the origin and (–1, 0) point.
2
30°
60° 1
90°
120° X
–2 –120° 1 2
–90°
–60° –1
–30°
–2
Q. 5. Sketch the Nyquist plot for the system with open loop
transfer function
60
G(s) H(s) =
(s 1 ) (s 2) (s 5)
and comment on stability.
Ans. Refer Q. 4.12, Unit-4.
Q. 9. Draw the bode plot for the transfer function given as, apply
correction to the magnitude plot for the quadratic term
and comment on stability.
5
G(s) H(s) =
s (1 0.1s) ( s2 0.4 s 1)
2
Control System
www.aktutor.in 5–1 C (EN-Sem-5)
5 Introduction to Design
PART-1
The Design Problems and Preliminary Consideration of Lead,
Lag and Lead-Lag Compensation Networks, Design of Closed
Loop System using Compensation Techniques in Time
and Frequency Domain.
Eo (s) 1 sT
=
Ei ( s) 1 sT
• Phase lead-lag compensation :
i. Lead-lag compensator is combination of lead network and lag
network.
ii. Transfer function,
1 1
s T1
s
T2
Gc(s) = Kc
1
s T1
s
T2
Questions-Answers
Answer
1. Fig. 5.1.1 shows a phase-lead network where in the phase of output
voltage leads the phase of input voltage for sinusoidal input.
C
R1
e1 R2 eo
Re
1 1
– –
T T
Bode plot :
1. The two corner frequencies are
1
= , lower corner frequency
T
1
= , upper corner frequency
T
2. The maximum phase lead m occurs at mid-frequency m between upper
and lower corner frequencies
1 1 1
log10 m = log log10
2 10 T T
1
m =
T
+20 dB/decade
0 dB
m
45°
0° Log
1 1 1
m =
T T T
Fig. 5.1.3. Bode plot of lead compensator.
3. The phase angle Eo(j)/Ei(j) can be calculated as
E ( j)
o = tan–1 (T) – tan–1 (T)
Ei ( j)
1
4. At = m = ,
T
The phase angle is
1 1 1
m = tan T tan 1 T
T T
1
= tan 1 1
tan
1
tan m =
1
1
Control System
www.aktutor.in 5–5 C (EN-Sem-5)
1
tan m =
2
1
and sin m =
1
Answer
Step 1 : Consider equation
1
s
Gc(s) = K c T
1
s
T
Let Kc = K = DC gain
Ts 1
Gc(s) = K
Ts 1
This control system is shown in Fig. 5.2.1.
+
Gc(s) G(s)
+
Fig. 5.2.1.
10
Que 5.3. For the open loop transfer function, G(s) H(s) =
s(1 0.2s)
design a suitable compensator such that the system will have a
phase margin of at least 45º. AKTU 2014-15, Marks 10
Answer
10
Given : G(s)H(s) =
s(1 0.2 s)
To Design : Compensator.
1. Given transfer function is of type-1 system,
10
G(s)H(s) =
s
s 1
5
2. Starting point = K1/n [Here n = 1]
= 10 (cut 0 dB axis)
3. Starting slope = – 20 dB/decade
4. P.M = 180° + , at = g = 35.54°
Here, p = 0 rad/s
and G.M = dB
m = 45° – 35.54° + tolerance
m = 45°– 35.54° + 5° = 14.46°
1 sin m 1 sin14.46
= =
1 sin m 1 sin14.46
= 0.6
1
m =
T
Table 5.3.1.
1 –101.30°
3 –120.96°
10 –153.43°
100 –177.13°
–180°
1
5. 20 log10 |Gc(j)| = – 20 log10
= – 20 log10 () – 1/2 = 10 log10= – 2.218 dB
1
6. Lower corner frequency = m
T
Introduction to Design
www.aktutor.in 5–8 C (EN-Sem-5)
Answer
K
Given : GP(s) = , Phase margin = 45°
s(s 10) (s 1000)
Kv = 1000 sec – 1
To Design : Phase lead compensator.
K
1. Kv = lim sG p (s) = lim s
s0 s 0 s( s 10) (s 1000)
K
Kv =
10 1000
K = 10 × 1000 Kv = 107
107 103
2. Gp(s) = =
s s s s
10000 s 1 1 s 1 1
10 1000 10 1000
3. Starting point at 0 dB axis at = K1/n
Here, n = 1 (Type of the system)
Control System
www.aktutor.in 5–9 C (EN-Sem-5)
= 1000 rad/sec
4. First corner frequency, c1 = 10 rad/sec
Second corner frequency, c2 = 1000 rad/sec
Initial slope = – 20 dB/sec
1
Gp(j) = – 90 – tan – 1 tan
10 1000
5. We have got p = g = 100 rad/sec.
So, gain margin = 0 dB and phase margin = 180° – 180° = 0°
m = 45° – 0° + tolerance = 45° – 0° + 5° = 50°
6. The value of parameter of the phase lead network is given by
1 sin m
= = 0.1324
1 sin m
In decibel,
1
At = m, magnitude in dB = – 20 log10 = 10 log10 = – 8.78 dB
Table 5.4.1.
∠ Gp(j)
1 – 95.76°
– 135.57°
40 – 168.25°
90 – 178.80°
Gain curve
Bode plot :
– 80°
– 100°
Introduction to Design
– 120°
– 40 dB/dec Angle
– 20 dB/dec – 140°
in
degree – 160°
p = g = 100 rad/sec
– 180°
– 200°
Gain in dB
– 220°
Phase curve
Fig. 5.4.1.
– 240°
– 260°
www.aktutor.in
– 280°
– 60 dB/dec
Answer
A. Lead compensator : Refer Q. 5.1, Page 5–3C, Unit-5.
B. Lag compensator :
1.
R1
R2
ei eo
C
0º
– 45º m
– 90º
log
1 m 1
T T
Fig. 5.5.2. Bode plot of lag compensator.
Introduction to Design
www.aktutor.in 5–12 C (EN-Sem-5)
1
4. The two corner frequencies are = , upper corner frequency for
T
1 1 1
zero at s = – ,= , lower corner frequency for a pole at s = –
T T T
5. The maximum phase-lag, m occurs at mid frequency m between upper
and lower corner frequencies.
1 1 1
log10 m = log log 10
2 T T
1
m =
T
6. The phase angle Eo(j)/Ei(j) calculated as
Eo ( j)
= tan–1 (T) – tan–1 (T)
Ei ( j)
1
At = m = , the phase angle is m :
T
1
tan m =
2
1
sin m =
1
Pole-zero configuration is shown in Fig. 5.5.3.
Im
Re
1 1
– –
T T
s-plane
Fig. 5.5.3. Pole zero configuration.
C. Phase margin and gain margin : Refer Q. 4.16, Page 4–25C, Unit-4.
Answer
Basically, lag compensator is a low pass filter and its main function is to
provide attenuation in the high frequency range to give sufficient phase
margin.
Control System
www.aktutor.in 5–13 C (EN-Sem-5)
Answer
1.
C1
R1
R2
vi(t) vo(t)
i(t) C2
1/sC1
Z2
R1
R2
Vi(s) Vo(s)
1
I(s)
sC2
1 (sR2C2 1)
and Z2 = R2 ...(5.7.2)
sC2 sC2
2. Applying Kirchhoff’s voltage law to the shown Fig. 5.7.1 network
Vi(s) = (Z1 + Z2) I(s) ...(5.7.3)
V0(s) = Z2I(s) ...(5.7.4)
Hence the transfer function of the system
Vo ( s)
Gc(s) = ...(5.7.5)
Vi ( s)
3. Putting value of V0(s) and Vi(s) in eq. (5.7.5), we have
Z2 I (s)
Gc(s) =
(Z1 Z2 ) I (s)
Z2
Gc(s) = ...(5.7.6)
Z1 Z2
4. Putting value of Z1 and Z2 in eq. (5.7.6), we have
(sR2C2 1)
sC2
Gc(s) =
R1 1
R2
sR1C1 1 sC2
(sR2C2 1)( sR1C1 1)
Gc(s) = s 2 R R C C sR C sR C sR C 1
1 2 1 2 2 2 1 1 1 2
1 1
s R C s R C
2 2 1 1
=
1 1 1 1
s2 s
R C
1 1 R C
2 2 R C
2 1 R1 1 R2C2
C
1 1
s T s T
1 2
Gc(s) =
1
s s
T1 T2
where, T1 = R1C1, T2 = R2C2
1 1 1 1
=
T1 T2 R1C1 R2C2 R2C1
and T1T2 = R1R2C1C2
and = 1
1
5. Poles are at s= ,
T1 T2
1 1
Zeros are at s = ,
T1 T2
Introduction to Design
www.aktutor.in 5–16 C (EN-Sem-5)
j
1 1 1
– – – –
T1 T2 T1 T2
1
C2 = , for a simple zero
T1
1
C3 = – , for a simple zero
T1
1
C4 = – , for a simple pole
T2
+20 dB/dec
–20 dB/dec
0dB
+45°
0°
– 45°
Log
– 1 1 1
– – –
T1 T1 T2 T2
Effects :
1. Lag-lead compensator increases the low frequency gain which improves
the steady state.
Control System
www.aktutor.in 5–17 C (EN-Sem-5)
PART-2
The Concept of State and Space, State-Space
Model of Physical System, Conversion of State-Space to
Transfer Function Model and Vice-Versa.
Questions-Answers
Que 5.8. Define state and state variable ? What are the
Answer
i. State : The state of a dynamic system is the smallest set of variables
such that the knowledge of these variables at t = t0 with the knowledge
of the input for t t0 completely determines the behaviour of the system
for any time t t0.
ii. State variables : The variables involved in determining the state of
dynamics system are called state variables.
iii. State vector : If we need n variable to completely describe the behaviour
of a given system, then these n state variables may be considered as n
component of a vector x. Such a vector is called state vector.
iv. State space : The n-dimensional space whose coordinate axes consists
of the x1 axis, x2 axis ..., xn axis is called state space. Any state can be
represented by a point in the state space.
v. State equation : A state space representation is a mathematical model
of the physical system as a set of output, input and state variables
related by a differential equation is known as state equation.
Advantages :
1. The method takes into account the effect of all initial conditions.
2. It can be applied to non-linear as well as time varying conditions.
3. It can be conveniently applied to multiple input multiple output systems.
4. The system can be designed for the optimal conditions precisely by using
this modern method.
5. Any type of the input can be considered for designing the system.
6. As the method involves matrix algebra, can be conveniently adopted for
the digital computers.
Answer
1. In a state variable system, state variables are represented by x1(t),
x2(t) ...
2. Input variables by u1(t), u2(t), .....
3. Output variables by y1(t), y2(t) ...
where,
x1 (t) u1 ( t) y1 (t)
x2 (t) u2 (t) y2 (t)
x(t) = x3 (t) ; u(t) = u3 (t) ; y(t) = y3 (t)
x (t ) u (t ) y (t)
n n n
Control System
www.aktutor.in 5–19 C (EN-Sem-5)
x n–1 = xn
8. Now from eq. (5.9.1)
x n = – anx1 – an–1 x2 .... – a1xn + u
9. We can write this set of equations in matrix form as,
x = Ax + Bu
Output : y = Cx
x1
x
2
where, x= ;
xn1
x
n
0 1 0 0 0
0
0 0 1 0
A= ;B=
0 0 0 1 0
a a 1
n n 1 an 2 a1
th
10. For n order single-input-single-output system :
x = Ax + Bu
y= Cx
A= System matrix of order (n × n)
B= Input coupling matrix of order (n × 1)
C= Output coupling matrix of order (1 × n)
x= State vector of order (n × 1)
u= Scalar input of order (1 × 1)
y= Scalar output of order (1 × 1)
Introduction to Design
www.aktutor.in 5–20 C (EN-Sem-5)
Answer
Basic ways to represent state models are :
a. State space representation using differential equation :
1. The state model is represented using these state variables called as
phase variables.
2. The phase variable is defined as those particular state variables which
are obtained from one of the system variables and its derivatives.
3. Transfer function can be obtained by differential equation
d n y(t) d n1 y( t) dy(t)
n
a1 .... an 1 an y(t) = b u(t)
dt dt n1 dt
y ( n ) a1 y ( n1) ..... an1 y (1) an = b u(t)
Put x1 = y
x2 = y x1
x3 =
y x2
xn = y ( n 1) x n 1
x1 = y
Control System
www.aktutor.in 5–21 C (EN-Sem-5)
x1 = x2 = y
x 2 = x3 = y
: : :
x n = y ( n) bu an ....a2 y ( n 2) a1 y ( n 1)
x n = bu an x3 an1 x2 .... a2 xn 1 a1 xn
0 1 0 0 x1 0
x 1 0 x
x 0 1 0 2 0
2 = 0 u
0 0 0 1 xn1
x
n a an 1 an 2 a1 xn b
n
x = Ax + Bu
x1
x
y = 1 0 0 0 2 [0]u
xn
u(t) x· n xn xn–1 x3 x2 x1 = y
b +
–
a1
+
+
+ + an–1
an
b
c
Fig. 5.10.2.
Introduction to Design
www.aktutor.in 5–22 C (EN-Sem-5)
Y (s) (s b1 ) (s b2 )
= K
U ( s) (s a1 ) (s a2 )
Each group is then decomposed by direct decomposition and the state
diagram of the transfer functions is cascaded.
U(s) + + + + C(s)
1/s b1 1/s b2
+ + +
+
a1 a2
(I) (II)
Fig. 5.10.3.
Y ( s) b0 s n b1 s n 1 .... bn
T(s) =
U ( s) a0 s n a1 s n1 .... an
c1 c cn
Then factorized T(s) = 2 ....
s a1 s a2 s an
Simulate all terms using direct decomposition and add in parallel to have
complete representation.
u x· 1 x1 + y
+ 1/s c1
+ +
+
1 +
x· 2 x2
+ 1/s c2
+
2
x· n xn
+ 1/s cn
+
n
Fig. 5.10.4.
Control System
www.aktutor.in 5–23 C (EN-Sem-5)
dx d2 x
output are y1 = 4 + 3u1, y2 = + 4u2 + u3.
dt dt 2
AKTU 2017-18, Marks 10
Answer
1. Select the state variables as
x1 = x
x1 = x = x2
x2 = x
3
x3 = – 3x3 – 4x2 – 4x1 + u1(t) + 3u2(t) + 4u3(t)
x1 0 1 0 x1 0 0 0 u1
x2 = 0 0 1 x2 0 0 0 u2
4 4 3 x3 1 3 4 u3
x3
2. Output : y1 = 4x2 + 3u1
y2 = x3 + 4u2 + u3
x1 u1
y1 0 4 0 3 0 0
y = x2 u2
2 0 0 1 0 4 1
x3 u3
Que 5.12. What is Transfer Function ? Also derive the expression
for transfer function of a state model.
Answer
A. Transfer Function : Transfer function is defined as the ratio of Laplace
transform of output of the system to the Laplace transform of input,
under the assumptions that all initial conditions are zero.
x = Ax + Bu
and output, y = Cx
2. Now, taking Laplace transform with zero initial conditions
sX(s) = AX(s) + BU(s)
X(s) = [sI – A]–1 BU(s)
and Y(s) = CX(s)
Y(s) = C [sI – A]–1 BU(s)
3. For a single-input-single-output system, Y and U are scalars.
4. Now transfer matrix can be given as
Y (s )
Transfer matrix = C[ sI A ]1 B
U ( s)
C adj([ sI A]) B
5. Transfer function = C[sI – A]–1 B =
det [sI A]
6. Denominator part i.e. |sI – A| is called the characteristic equation.
|sI – A|= 0
th
7. n degree characteristic equation |sI – A| = 0 has n roots or eigen
values.
Que 5.13. Write the state variable formulation of the parallel RLC
network.
i A
iR iC iL
R C
L
i(t) = I 0 sin t VC
Fig. 5.13.1.
Answer
1. KCL at node A,
iR + iC + iL = I0 sin t
VC (t) 1 CdVC (t)
VC (t) dt = I0 sin t
R L dt
2. Differentiating w.r.t. t and dividing by C,
Control System
www.aktutor.in 5–25 C (EN-Sem-5)
Answer
i. State variables : Refer Q. 5.8, Page 5–17C, Unit-5.
ii. State : Refer Q. 5.8, Page 5–17C, Unit-5.
iii. State vector : Refer Q. 5.8, Page 5–17C, Unit-5.
iv. State space : Refer Q. 5.8, Page 5–17C, Unit-5.
v. State equation : Refer Q. 5.8, Page 5–17C, Unit-5.
State transition matrix : The matrix (t) = exp (At) is an n × n
matrix and it helps in transition from initial state X(0) to any other
state x(t) for t > 0, hence (t) is called state transition matrix.
Derivation :
1. Consider homogeneous equation, with u(t) = 0
x = Ax
Introduction to Design
www.aktutor.in 5–26 C (EN-Sem-5)
Laplace transform,
sX(s) – x(0) = A X(s)
(sI – A) X(s) = x(0)
I = Identity matrix
s = Scalar Laplace operator.
2. Multiplying both sides by [sI – A]–1
X(s) = [sI – A]–1 x(0)
X(s) = (s) x(0)
where,
(s) = [sI – A]–1
1
1 A I A A2
I 2 3 ...
(s) =
s s s s s
3. Taking inverse Laplace transform of X(s),
x(t) = L–1 [X(s)] = L–1 [sI – A]–1 x(0)
1 I A A2
= L 2 3 ... x(0)
s s s
A2 t2
= I At ... x(0)
2
x(t) = eAt x(0) = (t) x(0)
4. State transition matrix
(t) = L–1 [sI – A]–1 = L–1 [(s)]
5. So, state transition matrix (STM) can be given as,
At 1 2 2 1 3 3
(t) = e I At A t A t ...
2 3
Properties of STM :
1. (0) = eA0 = I
Proof : (0) = eA × 0 = I
2. –1(t) = (– t)
1 1
Proof : –1(t) = = e–At = (–t)
(t) e At
3. (t2 – t1) (t1 – t0) = (t2 – t0)
Proof :
(t2 – t1 )(t1 – t0 ) = e A(t2 – t1 ) e A(t1 – t0 ) = e A(t2 – t1 + t1 – t0 ) = e A(t2 – t0)
= (t2 – t0)
4. [(t)] k = (kt)
Proof : (t)k = (t).(t) ... k times = eAt eAt ... k times = eAkt = (kt)
5. (t1 + t2) = (t1) (t2) = (t2) (t1)
Proof :
(t1 + t2) = eA(t1 + t2 ) = eAt1 eAt2 = eAt2 eAt1 = (t1) (t2) = (t2) (t1)
Control System
www.aktutor.in 5–27 C (EN-Sem-5)
Que 5.15.
i. Derive the transfer function from state model.
ii. Obtain the complete solution of non-homogeneous state
equation using time domain method.
iii. Discus s the significance of lag network. Als o draw its
s-plane representation and Bode plot.
AKTU 2016-17, Marks 10
OR
What are homogeneous and non-homogeneous systems ? Derive
the solution of the two systems in terms of the state variables.
Answer
i. Derivation of the transfer function from state model : Refer
Q. 5.12, Page 5–23C, Unit-5.
ii.
A. Homogeneous system : If in a state model of a system, the matrix A is
a constant matrix and input control forces zero, then the equation takes
the form, x (t) = A X(f), such an equation is called homogeneous equation
and the system is called homogeneous system.
Solution of homogeneous state equation :
1. x = Ax {u(t) = 0 for homogeneous equation}
2. Taking Laplace transform
sX(s) – x (0) = A X(s)
Hence, [sI – A] X(s) = x (0)
where I is identity matrix and s is the scalar Laplace operator.
3. Premultiplying both side by [sI – A]–1
X(s) = [sI – A]–1 x(0) = (s) x(0)
4. Taking the inverse Laplace transform of X(s)
1 A A2
x(t) = L–1 2 3 ... x(0)
s s s
x(t) = (t) x(0)
(t) = eAt is (n × n) matrix and is called State Transition Matrix.
(t) is unique solution of
d(t)
= A (t) , (0) = I
dt
B. Non-homogeneous system : If in a state model of a system, if A is a
constant matrix and input control forces are applied to the system, then
the system is called non-homogeneous system.
Solution of non-homogeneous equation :
1. Consider state equation
x (t) = Ax(t) + Bu(t)
Introduction to Design
www.aktutor.in 5–28 C (EN-Sem-5)
U(s) = L[u(t)]
X(s) = L[x(t)]
2. Taking Laplace transform,
sX(s) – x(0) = AX(s) + BU(s)
[sI – A] X(s) = x(0) + BU(s)
3. Premultiplying by [sI – A]–1
X(s) = [sI – A]–1 x(0) + [sI – A]–1 BU(s)
4. Taking inverse Laplace transform
x(t) = L–1{X(s)}
= L–1[sI – A]–1 x(0) + L–1[sI – A]–1 BU(s)
t
Answer
Y (s) 1
Given : = 3 2
U ( s) s 2 s 3s 1
To Find : State representation.
Control System
www.aktutor.in 5–29 C (EN-Sem-5)
Y (s) 1
Given : =
U ( s) (s 1) (s 4)
To Find : State representation.
1. Y(s) [(s2 + 5s + 4)] = U(s)
Taking inverse Laplace
d 2 y 5 dy
+ 4y(t) = u(t)
dt 2 dt
d2 y
2. = y2 (t) = y3(t)
dt 2
dy
(t ) = y2(t) = y1 ( t)
dt
y(t) = y1 (t)
3. y2 (t) + 5y (t) + 4y (t) = u(t)
2 1
y2 (t) = – 5y (t) – 4y (t) + u(t)
2 1
y1 (t) = 0 1 y1 ( t) 0
– 4 – 5 y ( t) 1 u(t)
y2 (t) 2
Introduction to Design
www.aktutor.in 5–30 C (EN-Sem-5)
C ( s) 10(s + 4)
=
R( s) s(s + 1) (s + 3)
Construct the phase variable state model.
AKTU 2013-14, Marks 10
Answer
C (s ) 10 (s 4)
Given : =
R( s) s(s 1) (s 3)
To Find : Phase variable state model.
C ( s) 10 (s 4)
1. =
R( s ) s(s2 4s 3)
C ( s) 10 (s 4)
=
R( s ) s 3 4 s 2 3s
C ( s) C( s) X (s)
2. =
R(s) X ( s) R(s)
X ( s) 1
= 3
R( s) s 4 s 2 3s
d 3 x (t ) d 2 x(t) dx(t)
3. 3
4 3 = r(t)
dt dt 2 dt
Let x(t) = x1 (t)
dx(t)
= x2(t) = x1 (t)
dt
d 2 x (t )
2
= x3 (t) = x2 (t)
dt
3
x3 (t) = d x(t) = r(t) – 3x2(t) – 4x3(t)
dt3
x 1 (t ) 0 1 0 x1 (t) 0
4. = 0 0 1 x ( t ) 0 r( t )
x 2 (t ) 2
0 3 4 x3 (t) 1
x3 (t)
C( s)
5. = 10 s + 40
X ( s)
C(s) = 10 s X(s) + 40 X(s)
dx(t)
c(t) = 10 + 40x (t)
dt
Control System
www.aktutor.in 5–31 C (EN-Sem-5)
Answer
Y ( s)
1. = (s2 + 3s + 2)
X ( s)
Taking inverse Laplace both sides
d 2 x(t) dx(t)
3 2 x(t) = y(t)
dt 2 dt
X ( s) 1
2. = 3
U (s ) s 9s 2 26 s 24
Taking inverse Laplace both sides
d 3 x (t ) d 2 x( t) dx(t)
3
9 26 24 x(t) = u(t)
dt dt 2 dt
3. Choose state variables
x(t) = x1 (t)
dx(t)
= x 1(t) = x2(t)
dt
d 2 x (t )
= x 2(t) = x3(t)
dt2
So, x 3(t) = 4(t) – 9x3(t) – 26x2(t) – 24x1(t)
y(t) = 2x1(t) + 3x2(t) + x3(t)
4. In state space representation form
x 1 (t ) 0 1 0 x1 (t) 0
x
2 ( t ) = 0 0 1 x2 (t) 0 u(t)
24 26 9 x3 (t) 1
x3 (t)
Introduction to Design
www.aktutor.in 5–32 C (EN-Sem-5)
x1 (t)
y(t) = [2 3 1] x2 (t)
x3 (t)
3
. . .
u(t) x 3(t) x 3(t) = x 2(t) x2 = x1 x1
2
y(t)
9
26
24
PART-3
Similarity Transformation of the control System,
Concept of Controllability and Observability and their Testing.
Questions-Answers
Answer
1. The choice of state is not unique for a given system. Suppose that there
exists a set of state variables
x = [x1 x2 x3 .... xn] T
so that a linear or similarity transformation
x = Mz; |M| 0, M–1 exists
i.e., z = M –1 x ...(5.19.1)
transforms to another set of state variables
Control System
www.aktutor.in 5–33 C (EN-Sem-5)
i.e.,
Bu
z = Az
and output y = C M z C z
i.e., y = C z
where A = M–1 A M
B = M–1 B
C = C M
where M is a non-singular transformation matrix.
3. Hence by similarity transformation, the transformed system can be
represented by the vector-matrix differential equation as
Bu
z = Az
and output y = C z
where A = M–1 A M
B = M–1 B
C = C M
Que 5.20. State and explain controllability and observability in
view of Kalman and Gilbert test.
3 1 1
The state equation for a system is, x = x u.
2 0 1
Answer
1. Controllability : A system is said to be controllable if the state of the
system can be transferred to another desired state over a given time
period by using input.
2. Observability : A system is said to be observable if all the state of the
system can be determined on the basis of knowledge of the output of the
system at given time.
Kalman’s test for controllability :
1. Consider nth order multiple input linear time invariant system
represented by its state equation as,
Introduction to Design
www.aktutor.in 5–34 C (EN-Sem-5)
x = Ax + Bu
2. The necessary and sufficient condition for the system to be completely
state controllable is that the rank of the composite matrix Qc is n.
3. The composite matrix Qc is given by,
Qc = [B : AB : A 2B : ... A n–1B]
Gilbert’s test for controllability :
1. For the Gilbert’s test it is necessary that the matrix A must be in
canonical form. Hence the given state model is required to be
transformed to the canonical form first, to apply the Gilbert’s test.
2. Consider single input linear time invariant system represented by,
z = Ax+Bu
where A is not in the canonical form. Then it can be transformed to the
canonical form by the transformation,
x = Mz
where M = Model matrix
3. The transformed state model,
z = A z B u
where
A = M – 1 AM
B = M – 1 B
4. In such a case the necessary and sufficient condition for the complete
state controllability is that the vector matrix B should not have any
zero elements. If it has zero elements then the corresponding state
variables are not controllable.
5. If the eigen values are repeated then matrix A cannot be transformed
to Jordan canonical form. If A has eigen values 1, 1, 2, 2, 3, 4, ... n
then the transformation results Jordan canonical form shown in matrix
below
1 1 0 0 0 0 ..........0
0 1 0 0 0 0...........0
J=
0 0 2 1 0 0...........0
Jordan 0 0 0 2 1 0...........0
block
0 0 0 0 2 0...........0
0 0 0 0 0 3..........0
........
........
0 0........................................................... n
Control System
www.aktutor.in 5–35 C (EN-Sem-5)
where C = CM
For a single input single output system,
z1
z
C
y = Cz C 12 .... C 1 n 2
11
zn
= C 11 z1 C 12 z1 .... C 1n zn
3. For the system to be observable, each term corresponding to each
state must be observed in the output. Hence none of the coefficient of
C must be zero.
Numerical :
3 1 1
Given : x =
2 0 x 1 u
To Find : Controllability.
3 1 1
1. A= ,B= 1 , n = 2
2 0
Introduction to Design
www.aktutor.in 5–36 C (EN-Sem-5)
2. Qc = B AB
3 1 1 2
3. AB = 1 = 2
2 0
1 2
4. Qc =
1 2
1 2
5. |Qc| = =2–2=0
1 2
Rank of Qc is r n
Hence, system is not completely controllable.
Y ( s) 2
= . Find the state and output equation in
U ( s) s3 6 s 2 11s 6
matrix form and also test the controllability and observability of
the given system. AKTU 2015-16, Marks 15
Answer
Y (s) 2
Given : = 3
U ( s) s 6s 2 11s 6
To Find : State matrix; Controllability, Qc; Observability, Qo.
1. Y(s) [(s3 + 6s2 + 11s + 6)] = 2U(s)
2. x1 = x2, x2 = x3, x3 = – 6x1 – 11x2 – 6x3 + u
3. y = 2x 1
x1 0 1 0 x1 0
= x = 0 0 1 x2 0 u
2
6 11 6 x3 1
x3
2 x1
y = 0 x2
0 x3
4. Controllability :
0 0 1 0
B = 0 , A =
0 0 1
1 6 11 6
Control System
www.aktutor.in 5–37 C (EN-Sem-5)
0 1 0 0 0
AB = 0 0 1
=
0
1
6 11 6 1 6
Qc = [B : AB : A2B]
0 1 0 0 1 0
A2 = 0 0 1
0 0 1
6 11 6 6 11 6
0 1 0
= 6 11 6
36 60 25
0 0 1 0 1
A2B = 6 11 6 0 = 6
36 60 25 1 25
0 0 1
1 6 0 6 0 1
Qc = 0 1 6
=
0 –0 1
6 25 1 25 1 6
1 6 25
Qc = 0 – 0 + (– 1)
Qc = –1
|Qc| 0
It is of rank 3 and controllable.
5. Observability :
0 1 0 0 0 6
A= 0 0 1
T
, A
1 0 11
6 11 6 0 1 6
2 2
C = 0 , CT = 0
0 0
Q0 = [CT : ATCT : (AT)2 CT]
0 0 6 2 0
ATCT = 1 0 22 0 2
0 1 6 0 0
0 0 6 0 0 6
(AT)2 = 1 0 11 1 0 11
0 1 6 0 1 6
0 6 36
= 0 11 60
1 6 25
Introduction to Design
www.aktutor.in 5–38 C (EN-Sem-5)
0 6 36 2 0
(AT)2CT = 0 11 60 0 0
1 6 25 0 2
2 0 0
Q o = 0 2 0
0 0 2
2 0 0 0 2 2
= 2 –0 0 = 2[4 – 0] – 0 + 0
0 2 0 2 0 0
|Qo |= 8
Therefore, its rank is 3 and the system is observable.
Answer
0 0 0 40
Given : A = 1 0 3 , B = 10 , C = [0 0 1]
0 0 4 0
To Test : Controllability and observability.
1. Controllability test :
0 0
i. AB = 40 , A 2B =
30
0 0
ii. Controllability test matrix is given by,
40 0 0
Qc = [B : AB : A2B] = 10 40 0
0 0 0
|Qc| = 0
iii. Thus the rank of Qc is 2. Hence the system is not controllable.
2. Obervability test :
0 1 0 0
i. AT = 0 0
0 , C T = 0
0 3 4 1
Control System
www.aktutor.in 5–39 C (EN-Sem-5)
0 0
ii. ATCT = 0 , (AT)2 CT = 0
4 16
iii. The observability test matrix is given by
Qo = [CT : ATCT : (AT)2 CT]
0 0 0
Qo = 0
0 0
1 4 16
iv. Its rank is 1. Hence the system is not observable.
Answer
0 1 0 0 1
Given : A = 0 0 1 , B = 0 0 and C = [1 0 0]
0 2 3 1 1
To Test : Controllability and Observability.
A. Controllability test :
0 1 0 0 1 0 0
1. AB = 0 0 1 0 0 1 1
0 2 3 1 1 3 3
0 1 0 0 0 1 1
A2B = 0 0 1 1 1 3 3
0 2 3 3 3 11 11
2. Controllability test matrix is given by
Qc = [B : AB : A2B]
0 1 0 0 1 1
Qc = 0 0 1 1 3 3
1 1 3 3 11 11
Introduction to Design
www.aktutor.in 5–40 C (EN-Sem-5)
0 1 0
3. Consider 0 0 1 = 1 0 ; |Qc| 0
1 1 3
Hence rank of Qc is equal to its order i.e., 3. Therefore the system is
controllable.
B. Observability test :
1
1. C = [1 0 0] ; C = 0
T
0
0 1 0 0 0 0
A = 0 0 1 ; AT = 1 0 2
0 2 3 0 1 3
0 0 0 1 0
1 0 2 0 1
AT CT =
0 1 3 0 0
0 0 0 0 0
(AT)2 CT = (AT) (ATCT) = 1 0 2 1 0
0 1 3 0 1
2. Observability test matrix is given by
1 0 0
T T T T 2 T 0 1 0
Qo = [C : A C : (A ) C ] =
0 0 1
|Qo| = 1 0
Hence its rank is 3 equal to its order i.e., 3, therefore system is completely
observable.
Control System SQ–1 C (EN-Sem-5)
www.aktutor.in
1 Control System
Concepts
(2 Marks Questions)
1.1. Explain open loop and closed loop system with physical
examples. AKTU 2015-16, Marks 02
OR
Discuss open loop and closed loop system giving suitable
example. AKTU 2016-17, Marks 02
OR
Draw the block diagram and explain the open loop control
system and closed loop control system.
AKTU 2017-18, Marks 02
OR
What are the major types of control systems ? Explain
them in detail with examples. AKTU 2017-18, Marks 02
Ans. There are two types of control system :
1. Open loop control system :
Fig. 1.1.1
In open loop control systems the control action is independent of
the desired output. In this system the output is not compared with
the reference input.
Example : Washing machine, Immersion rod, Time operated
traffic control, DC shunt motor.
2. Closed loop control system :
Error signal
Reference Controlled
Amplifier Controlled process
input + – output
+
Feedback
Fig. 1.1.2.
2 Marks Questions SQ–2 C (EN-Sem-5)
www.aktutor.in
In a closed loop control system the output has an effect on control
action through a feedback as shown in Fig. 1.1.2.
Example : Automatic steering control system, Driving system of
an automobile, Home heating system, Ship stabilization system.
Control From
Vc Amplifier
winding controller
Reference winding
Vr
1.5. What is the basis for framing the rules of block diagram
reduction technique ? AKTU 2017-18, Marks 02
Ans.
1. Series connection :
R(s) Y(s) C(s) R(s) C(s)
G1(s) G2(s) G1(s)G2(s)
2. Parallel connection :
+ C(s)
R(s) G1(s)
R(s)
±
G1(s) ± G 2(s) C(s)
G2(s)
e1(t) C1 R2 e2(t)
Fig. 1.13.1.
1
R1 C s
1
Vi(s) = R2 I (s)
R1 1
C1 s
R1
+
1/C1s R2
Vi(s) Vo(s)
I(s)
–
Fig. 1.13.2.
2. And, V0(s) = R2 I(s)
R2 Vi ( s) R2 Vi ( s)
V0(s) =
1 R1
R1 C s R C s 1 R2
1
R2 1 1
R 1
1 C1 s
V ( s) R2 ( R1 C1 s 1)
Transfer function, 0 = .
Vi (s) R2 ( R1 C1 s 1) R1
2 Marks Questions SQ–6 C (EN-Sem-5)
www.aktutor.in
2.2. Define rise time and delay time for second order control
system. AKTU 2015-16, Marks 02
Ans.
1. Rise Time (tr) : Rise time is defined as the time required for the
response to rise from 10 % to 90 % (for overdamped) or from 0 %
to 100 % (for underdamped) of its final value.
1 2
tan 1
Rise time (tr) =
d
2. Delay Time (td) : Delay time is the time required for the response
to reach the half of the final value in the first attempt.
The delay time in terms of and n is given by :
1 0.7
td =
n
It is denoted by ess.
0.125 1 4
=
4a
1
e ss =
4a
Fig. 2.15.1.
we get, n2 = 16
n = 4
3. Damping factor is given as = 0.5
2 n = (0.8 + 16a)
2 × 0.5 × 4 = 0.8 + 16a
4 = 0.8 + 16a
a = 0.2
Control System SQ–11 C (EN-Sem-5)
www.aktutor.in
3 Stability and
Algebraic Criteria
(2 Marks Questions)
3.11. What will be the nature of impulse response when the roots
of characteristic equation are lying on imaginary axis ?
Ans. If the roots of characteristic equation lie on imaginary axis the
nature of impulse response is oscillatory.
3.13. How will you find the gain K at a point on root locus ?
Ans. The gain K at a point s = sa on root locus is given by,
Product of length of vector from open loop poles to the point sa
K=
Product of length of vector from open loop zeros to the point sa
2 Marks Questions SQ–14 C (EN-Sem-5)
www.aktutor.in
4 Frequency
Response Analysis
(2 Marks Questions)
0
Im
=0
Re
=
Fig. 4.13.1.
0 = =0
Re
= – tan –1 T
a
K
2 + a2
G(j ) – plane
Fig. 4.14.1.
4. Diameter = |G(j)|=0 – |G(j)|=
K K
= –0
a a
So, Radius = K / 2a
G ( j) 0 – G( j)
5. Center = = K/2a.
2
2 Marks Questions SQ–18 C (EN-Sem-5)
www.aktutor.in
5 Introduction to Design
(2 Marks Questions)
1 5 0 5 0 1
Ans. |Qc| = 0 0 1
5 11 1 11 1 5
|Qc| = [0 – 0] + [0 – 1] = – 1
Thus the rank of Qc is 3 and order is also 3. Hence, the system is
completely controllable.
Control System SP–1 C (EN-Sem-5)
www.aktutor.in
B.Tech.
(SEM. V) ODD SEMESTER THEORY
EXAMINATION, 2015-16
CONTROL SYSTEM
SECTION–A
Note : Attempt all sections. All sections carry equal marks. Write answer
of each part in short. (2 × 10 = 20)
1. a. Explain open loop and closed loop system with physical
examples.
1
e. Draw the polar plot of open loop transfer function .
s2
h. Define rise time and delay time for second order control
system.
SECTION–B
7. Sketch the Nyquist plot for the system with open loop
transfer function
60
G(s) H(s) =
(s 1 ) (s 2) (s 5)
and comment on stability.
SECTION–C
Solved Paper (2015-16) SP–4 C (EN-Sem-5)
www.aktutor.in
SOLUTION OF PAPER (2015-16)
SECTION–A
Note : Attempt all sections. All sections carry equal marks. Write answer
of each part in short. (2 × 10 = 20)
1. a. Explain open loop and closed loop system with physical
examples.
Ans. There are two types of control system :
1. Open loop control system :
Input Output
Controller Controlled process
Fig. 1.
In open loop control systems the control action is independent of
the desired output. In this system the output is not compared with
the reference input.
Example : Washing machine, Immersion rod, Time operated
traffic control, DC shunt motor.
2. Closed loop control system :
Error signal
Reference Controlled
Amplifier Controlled process
input + – output
+
Feedback
Fig. 2.
In a closed loop control system the output has an effect on control
action through a feedback as shown in Fig. 2.
Example : Automatic steering control system, Driving system of
an automobile, Home heating system, Ship stabilization system.
1
e. Draw the polar plot of open loop transfer function .
s2
Ans.
Given : G(s) H(s) = 1 / s2
To Find : Polar plot.
1. Putting s = j
1 –1
G(j) H(j) = = 2
( j)2
1
2. Magnitude, M = |G(j) H(j)|=
2
and phase angle, = – 180°
M
0
2 1/4
4 1/16
0
Im
=0
Re
=
Fig. 3.
we get, n2 = 16
n = 4
3. Damping factor is given as = 0.5
2 n = (0.8 + 16a)
2 × 0.5 × 4 = 0.8 + 16a
4 = 0.8 + 16a
a = 0.2
h. Define rise time and delay time for second order control
system.
Ans.
1. Rise Time (tr) : Rise time is defined as the time required for the
response to rise from 10 % to 90 % (for overdamped) or from 0 %
to 100 % (for underdamped) of its final value.
1 2
tan 1
Rise time (tr) =
d
2. Delay Time (td) : Delay time is the time required for the response
to reach the half of the final value in the first attempt.
The delay time in terms of and n is given by :
1 0.7
td =
n
SECTION–B
Fig. 5.
Ans. Step 1 : Shift S2 before G1, we get
H3 / G1
S1 – + S3
R(s)
G1 G2 G3 G4
– + – T2 C(s)
S2
H1 H2
Fig. 6.
Solved Paper (2015-16) SP–8 C (EN-Sem-5)
www.aktutor.in
Step 2 : G1 and G2 are in cascade.
H3 / G1
S1 – S3
R(s) + C(s)
G1 G2 G3 G4
– + – T2
S2
H1 H2
Fig. 7.
Step 3 : Interchanging summing points
H3 / G1
R(s) – + + S3 C(s)
+ G1 G2 G3 G4
– S1 – T2
S2
H1 H2
Fig. 8.
Step 4 : Solve the inner loop.
H3 / G1
R(s) – G1 G2 + S3 C(s)
+ G3 G4
1+G 1G 2H1 – T2
S2
H2
Fig. 9.
Step 5 : Shift T2 after G4.
H3/G1 G4
– S3
+ G1 G2 + C
G3 G4
R(s) 1+G1 G2 H1 – T2
S2 H2
Fig. 10.
Step 6 : Solve the inner loop.
H3/G1 G4
–
G1 G2 G3 G4
C(s)
R(s) + 1+G1 G2 H 1 1+G3 G4 H2
Fig. 11.
Step 7 : Forward blocks are in cascade and final result will be
G1G2G3G4
C ( s) (1 G1G2 H1 ) (1 G3G4 H2 )
=
R( s) G2G3 H3
1
(1 G1G2 H1 ) (1 G3G4 H2 )
Control System SP–9 C (EN-Sem-5)
www.aktutor.in
C ( s) G1G2G3G4
=
R( s) (1 G1G2 H1 ) (1 G3G4 H2 ) G2G3 H3
R(s) + 2n
s(s + 2 n) C(s)
–
Fig. 12.
The closed loop transfer function is
C ( s) 2n / s( s 2 n )
=
R( s) 1 2n / s( s 2 n )
where, = Damping factor or Damping ratio
C ( s) 2n
= 2
R( s) ( s 2 n s 2n )
Then output,
2n
C (s) = R s 2 ...(1)
s 2 n s 2n
2. For unit step input
r (t) = 1
1
R (s) =
s
1 2n
Then C (s) = ...(2)
s s2 2n s 2n
3. In eq. (2) putting [s2 + 2ns + 2n ] = [( s n )2 n2 (1 2 )] and
breaking it into partial fraction
1 s 2 n
C (s) = ...(3)
s [( s n )2 2n (1 2 )]
Put d = n 1 2
1 s 2 n
C (s) = ...(4)
s (s n )2 d2
4. Rewrite eq. (4)
1 s n n d
C(s) = ...(5)
s (s n )2 d2 d (s n )2 2d
5. Taking inverse Laplace transform of eq. (5),
Solved Paper (2015-16) SP–10 C (EN-Sem-5)
www.aktutor.in
t n nt
c(t) = 1 e n cos d t e sin d t
d
Put d = n 1 2
ent
c(t) = 1 [( 1 2 )cos d t sin d t ] ...(6)
1 2
6. Put sin = 1 2 ,
cos =
e n t
c(t) = 1 (sin cos d t cos sin d t)
1 2
e n t
c(t) = 1 sin( d t ) ...(7)
1 2
where d = n 1 2
1 2
and = tan 1
7. Eq. (7) is rewritten as
1 2
e n t
c(t) = 1 sin ( n 1 2 )t tan 1 ...(8)
1 2
c(t)
<1
t
Fig. 13. Underdamped oscillations.
Control System SP–11 C (EN-Sem-5)
www.aktutor.in
4. Using Routh’s stability criterion, determine the range of
K for open loop transfer function
K
G(s) H(s) =
s(s 1 ) (1 2 s)
Ans.
K
Given : G(s) H(s) =
s(s 1 ) (1 2s)
To Find : Range of K.
1. Using characteristics equation
1 + G(s) H(s) = 0
K
1+ = 0
s(s 1)(1 2s)
s(s + 1) (1 + 2s) + K = 0
s(2s2 + 3s + 1) + K = 0
2s3 + 3s2 + s + K = 0
2. Using Routh array :
s3 2 1
2
s 3 K
1 3 – 2K
s 0
3
s0 K
3. For system to be stable :
3 – 2K
> 0 and K > 0
3
3
0<K<
2
K
Given : G(s) H(s) = , H(s) = 1, = 0.5
s(s 1) (s 3)
To Sketch : Root locus and value of K at = 0.5.
1. The open-loop poles are at s = 0, s = – 1, s = – 3.
2. There is no open loop zeros.
3. Number of poles, P = 3
Number of zeros, Z = 0
Solved Paper (2015-16) SP–12 C (EN-Sem-5)
www.aktutor.in
P – Z = 3 – 0 = 3 i.e., three branches of root locus end at infinity.
4. Angle of asymptotes :
(2 q 1)180
q = where q = 0, 1, 2, ...(P – Z –1)
PZ
(2 0 1) 180
0 = = 60°
(3 0)
(2 1 1) 180
1 = = 180°
(3 0)
(2 2 1) 180
2 = = 300°
(3 0)
5. Centroid of asymptotes :
Real parts of poles Real parts of zeros
=
PZ
(0 1 3) (0)
= = – 1.33.
30
6. Breakaway points : Between open-loop poles s = 0 and s = – 1,
there exists a breakaway point.
The characteristic equation is
s (s + 1)(s + 3) + K = 0
K = – (s3 + 4s2 + 3s)
dK
= – (3s2 + 8s + 3) = 0
ds
3s2 + 8s + 3 = 0
s = – 1.33 ± 0.88 = – 0.42 and – 2.21
As the breakaway point has to lie between s = 0 and s = – 1, the valid
breakaway point is s = – 0.42.
7. Intersection with j axis :
Characteristic equation, 1 + G(s) H(s) = 0
s (s + 1) (s + 3) + K = 0
s3 + 4s2 + 3s + K = 0
Routh array :
s3 1 3
s2 4 K
1
s (12 K ) / 4
s0 K
The value of K at imaginary axis :
(12 K )
= 0 K = 12
4
Auxiliary equation,
4s2 + 12 = 0
s = ± j 1.73
Control System SP–13 C (EN-Sem-5)
www.aktutor.in
Im
K=
(j1.732, K = 12)
– 0.35 + j0.6
180° 60°
K= K=0 K=0
–3 –1 0 Re
–1.33
300°
–0.45
(– j1.732, K = 12)
K=
Fig. 14. Root locus for G(s) = K/s(s + 1)(s + 3).
1
M= 2
1 – j 2
n n
Let =x
n
1 1
M=
1 x 2 j 2x (1 x 2 )2 4 2 x 2
dM
3. To maximize magnitude, put = 0,
dx
Solved Paper (2015-16) SP–14 C (EN-Sem-5)
www.aktutor.in
dM d
= [(1 – x 2 )2 4 2 x 2 ]–1/ 2
dx dx
1 d
= – [(1 – x 2 )2 4 2 x 2 ]–3/ 2 [(1 – x 2 )2 4 2 x2 ]
2 dx
1 1
= – [2(1– x 2 ) (– 2 x) 8 2 x ]
2 [(1 – x 2 )2 4 2 x 2 ]3/ 2
=0
4x[x2 + 22 – 1] = 0
As, x 0
x2 = 1 – 22
2
x = 1 – 2
Resonant frequency,
r = n 1 – 2 2
1
4. |Mr| =
(1 x ) 4 2 x 2
2 2
1
|Mr| =
(1 ( 1 2 )) 4 2 ( 1 2 2 )2
2 2
1
|Mr| =
4 4 4 2 (1 2 2 )
Resonant peak,
1
Mr =
2 1 2
7. Sketch the Nyquist plot for the system with open loop
transfer function
60
G(s) H(s) =
(s 1 ) (s 2) (s 5)
and comment on stability.
Ans.
60
Given : G(s) H(s) =
(s 1 ) (s 2) (s 5)
To Sketch : Nyquist plot.
60
1. Put s = j, G(j) H (j) =
( j 1 ) ( j 2) ( j 5)
60
2. Magnitude, M = |G (j) H (j)| =
( 1)( 4)( 2 25)
2 2
3. Phase = G(j)H(j)
Control System SP–15 C (EN-Sem-5)
www.aktutor.in
1 1 1
= tan () tan tan
2 5
Table 1.
M
0 6 0°
1.12 3.402 – 90.11°
4.12 0.4769 – 180°
10 0.0523 133.58
0 90° (– 270°)
4. N= P–Z
P= 0 (No pole on RH plane)
N= 0 (No encirclement of – 1 + j0 point)
Z= 0 (Number of roots of closed loop characteristic
equation having positive real parts).
So the given system is stable system.
Im
–0
Re
– 1 + j0 –
=+0
0.4769
Fig. 15.
Fig. 16.
3. Assume, K = 1, we can write,
(1 sTd ) 2n
G(s) =
s( s 2 n )
C ( s) (1 sTd ) 2n
= 2
R( s) s s [2 n 2n Td ] 2n
4. Comparing denominator with standard form,
2n = 2 n 2n Td
n Td
=
2
n Td
5. Because of this controller, damping ratio increased by factor .
2
Kp = lim G( s) H ( s)
s0
n
Kv = lim s G(s) H (s)
s 0 2
without PD controller,
2n
G(s) H(s) = [K = 1]
s(s 2 n )
n
Kv = lim sG(s) H ( s)
s0 2
As there is no change in coefficients, error also will remain same.
B. Synchros :
1. It is a rotary transducer that converts angular displacement into
an AC voltage or an AC voltage into an angular displacement.
2. A synchros system consists of
i. A control transmitter (CX) and
ii. A control transformer (CT).
3. The control transmitter consists of a stator and a rotor. The rotor
is a dumb-bell shaped magnetic structure. The supply is given to
the rotor by means of slip rings, which are actually mounted on
the stator housing.
4. The secondaries are in the skewed slot all along the periphery of
the stator and are 120° apart because of their mechanical
displacement.
Control System SP–17 C (EN-Sem-5)
www.aktutor.in
5. The induced secondary voltage will depend upon the angle of the
rotor shaft. For reference the zero degree position of the shaft is
defined when the rotor is in alignment with the coil S2.
6. In this position, the voltage in coil S2 is maximum, and similarly
the maximum voltage in coils S1 and S3 will result at 120° and 240°
positions respectively. The voltage in S2 is a function of and so is
the voltage in S1 and S3. Thus,
E0s2 = A cos
E0s1 = A cos ( – 120°)
E0s3 = A cos ( – 240°)
S2
Stator
Slip rings
R1
AC supply
R2 S1
S3
Y (s) 1
Given : = 3
U ( s) s 2s2 3s 1
To Find : State representation.
1. Y(s) [s3 + 2s2 + 3s + 1] = U(s)
Taking inverse Laplace transform
d 3 y 2 d 2 y 3dy
y( t) = u(t)
dt3 dt 2 dt
2. Choosing state variable
y(t) = y1 (t)
dy( t)
= y1 ( t) = y2(t)
dt
d 2 y(t)
= y2 (t) = y3(t)
dt
d 3 y(t)
= y3 (t) = y4(t)
dt
3. y4 + 2y3 + 3y2 + y1(t) = u(t)
y3 = – 2y3 – 3y2 – y1(t) + u(t)
0 1 0 y1 0
y1 ( t)
y (t) = 0 0
1 y2 0 u(t)
2 –1 –3 –2 y3 1
y (t)
3
Y (s) 1
Given : =
U ( s) (s 1) (s 4)
To Find : State representation.
1. Y(s) [(s2 + 5s + 4)] = U(s)
Taking inverse Laplace
d 2 y 5 dy
+ 4y(t) = u(t)
dt 2 dt
d2 y
2. = y2 (t) = y3(t)
dt 2
dy
(t) = y2(t) = y1 ( t)
dt
y(t) = y1 (t)
3. y2 (t) + 5y (t) + 4y (t) = u(t)
2 1
y2 (t) = – 5y (t) – 4y (t) + u(t)
2 1
Control System SP–19 C (EN-Sem-5)
www.aktutor.in
y1 (t) = 0 1 y1 ( t) 0
– 4 – 5 y ( t) 1 u(t)
2
y2 (t)
SECTION–C
2(s 0.25)
Given : G(s) H(s) =
s 2 ( s 1) ( s 0.5)
To Draw : Bode plot.
2 0.25 s
1
0.5 0.25 (4 s 1)
1. G(s) H(s) = = 2
2 s s (s 1) (2s 1)
s ( s 1) 1
0.5
2. K = 1, 20 Log K = 0 dB, No effect on Bode Plot.
1
3. , 2 Poles at origin, straight line of slope – 40 dB/dec
s2
1
4. (1 + 4s), simples zero, T1 = 4, C1 = = 0.25,
T1
straight line of slope + 20 dB/dec for 0.25
1 1
5. , simple zero, T2 = 2, C2 = = 0.5,
1 2s T2
straight line of slope – 20 dB/sec for 0.5
1 1
6. , simple zero, T3 = 1, C3 = = 1,
1 s T3
straight line of slope – 20 dB/sec for 1
Table 2.
Range of 0 < < 0.25 0.25 < 0.5 0.5 < 1 1 <
Resultant slope – 40 – 40 + 20 – 20 – 20 – 40 – 20
in dB/dec = – 20 = – 40 = – 60
Solved Paper (2015-16) SP–20 C (EN-Sem-5)
www.aktutor.in
(1 4 j)
7. Phase angle table : G(jω)H(jω) =
( j)2 (1 2 j) (1 j)
Bode plot :
1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
100
80
– 60 dB/dec
60
– 40 dB/dec
40
– 20 dB/dec
20
G.M. = – 21 dB
– 20 dB/dec – 60 dB/dec
in deg
–150°
–180°
P.M. = – 36°
–210°
–240°
–270°
0.1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
Control System SP–21 C (EN-Sem-5)
www.aktutor.in
Table 3.
1
+ tan– 14 + tan– 12 – tan– 1
( j)2
Y (s) 2
Given : = 3
U ( s) s 6s 2 11s 6
To Find : State matrix; Controllability, Qc; Observability, Qo.
1. Y(s) [(s3 + 6s2 + 11s + 6)] = 2U(s)
2. x1 = x2, x2 = x3, x3 = – 6x1 – 11x2 – 6x3 + u
3. y = 2x 1
x1 0 1 0 x1 0
= x = 0 0 1 x2 0 u
2
6 11 6 x3 1
x3
2 x1
y = 0 x2
0 x3
4. Controllability :
0 0 1 0
B = 0 , A =
0 0 1
1 6 11 6
0 1 0 0 0
AB = 0 0 1
=
0
1
6 11 6 1 6
Solved Paper (2015-16) SP–22 C (EN-Sem-5)
www.aktutor.in
Qc = [B : AB : A2B]
0 1 0 0 1 0
A2 = 0 0 1
0 0 1
6 11 6 6 11 6
0 1 0
= 6 11 6
36 60 25
0 0 1 0 1
2
AB= 6 11 6 0
= 6
36 60 25 1 25
0 0 1
1 6 0 6 0 1
Qc = 0 1 6
=
0 –0 1
6 25 1 25 1 6
1 6 25
Qc = 0 – 0 + (– 1)
Qc = –1
|Qc| 0
It is of rank 3 and controllable.
5. Observability :
0 1 0 0 0 6
A= 0 0 1
T
, A
1 0 11
6 11 6 0 1 6
2 2
C = 0 , CT = 0
0 0
Q0 = [CT : ATCT : (AT)2 CT]
0 0 6 2 0
A C = 1 0 22 0 2
T T
R1
e1 R2 eo
Re
1 1
– –
T T
Fig. 19. Pole zero configuration.
Bode plot :
1. The two corner frequencies are
1
= , lower corner frequency
T
1
= , upper corner frequency
T
2. The maximum phase lead m occurs at mid-frequency m between
upper and lower corner frequencies
1 1 1
log10 m = log log10
2 10 T T
1
m =
T
+20 dB/decade
0 dB
m
45°
0° Log
1 1 1
m =
T T T
Fig. 20. Bode plot of lead compensator.
3. The phase angle Eo(j)/Ei(j) can be calculated as
E ( j)
o = tan–1 (T) – tan–1 (T)
Ei ( j)
1
4. At = m = ,
T
The phase angle is
1 1 1
m = tan T tan 1 T
T T
Control System SP–25 C (EN-Sem-5)
www.aktutor.in
1
= tan 1 1
tan
1
tan m =
1
1
1
tan m =
2
1
and sin m =
1
B. Lag compensator :
1.
R1
R2
ei eo
C
R1 R2
where > 1, =
R2
and T = R2 C
2. The transfer function given by eq. (1) can be expressed in sinusoidal
form as
Eo ( j) 1 jT
= ...(2)
Ei ( j) 1 jT
3. Bode plot for transfer function of eq. (2) is shown in Fig. 22.
1
4. The two corner frequencies are = , upper corner frequency
T
1 1
for zero at s = – ,= , lower corner frequency for a pole at
T T
1
s=–
T
Solved Paper (2015-16) SP–26 C (EN-Sem-5)
www.aktutor.in
– 20 dB/decade
0 dB
0º
– 45º m
– 90º
log
1 m 1
T T
Fig. 22. Bode plot of lag compensator.
5. The maximum phase-lag, m occurs at mid frequency m between
upper and lower corner frequencies.
1 1 1
log10 m = log log 10
T
2 T
1
m =
T
6. The phase angle Eo(j)/Ei(j) calculated as
Eo ( j)
= tan–1 (T) – tan–1 (T)
Ei ( j)
1
At = m = , the phase angle is m :
T
1
tan m =
2
1
sin m =
1
Pole-zero configuration is shown in Fig. 23.
Im
Re
1 1
– –
T T
s-plane
Fig. 23. Pole zero configuration.
C. Phase margin and gain margin :
a. Gain Margin (GM) :
1. It is the margin in gain which is allowed till the system reaches
on the verge of instability.
2. Mathematically, gain margin is reciprocal of magnitude of
G(j)H(j) at phase crossover frequency.
Control System SP–27 C (EN-Sem-5)
www.aktutor.in
GM = 1 /|G ( j) H ( j)
pc
Control System SP–1 C (EN-Sem-5)
www.aktutor.in
B.Tech.
(SEM. V) ODD SEMESTER THEORY
EXAMINATION, 2016-17
CONTROL SYSTEM
Time : 3 Hours Max Marks : 100
SECTION-A
g. Explain in brief :
i. Gain margin
ii. Phase margin.
SECTION-C
–H 1 – H2 –H3
H4
Fig. 2.
Fig. 3.
10 20 30 40
Fig. 4.
b. Sketch the root locus for the closed loop control system
K
with G(s) =
s(s 1)( s2 4 s 5)
3 1 1
The state equation for a system is, x = x u
2 0 1
Control System SP–5 C (EN-Sem-5)
www.aktutor.in
SOLUTION OF PAPER (2016-17)
SECTION-A
Fig. 1.
In open loop control systems the control action is independent of
the desired output. In this system the output is not compared with
the reference input.
Example : Washing machine, Immersion rod, Time operated
traffic control, DC shunt motor.
2. Closed loop control system :
Error signal
Reference Controlled
Amplifier Controlled process
input + – output
+
Feedback
Fig. 2.
In a closed loop control system the output has an effect on control
action through a feedback as shown in Fig. 2.
Example : Automatic steering control system, Driving system of
an automobile, Home heating system, Ship stabilization system.
Control From
Vc Amplifier
winding controller
Reference winding
Vr
g. Explain in brief :
i. Gain margin
ii. Phase margin.
Solved Paper (2016-17) SP–8 C (EN-Sem-5)
www.aktutor.in
Ans.
i. Gain margin : It is the reciprocal of magnitude |G(j)| at the
frequency at which the phase angle is – 180°.
Gain Margin (GM),
1
Kg =
|G( j c )|
where, c = Phase cross-over frequency.
ii. Phase margin : The phase margin is that amount of the additional
phase lag at the gain crossover frequency required to bring the
system to the verge of instability.
Phase margin is equal to 180° plus the angle of G(j) at the gain
crossover point.
m = 180° +
SECTION-B
Step 2 :
R(s) + G4 C(s)
+ G1 + G2 G3 +
– – – G2
H2
H1
Step 3 :
L1
R(s) G4
+ + G1 –
G2 G3 + C(s)
– – G2
H2
G
H1 G3 + 4
G2
R(s) G2G3 + G4
+ +– G1 C(s)
– 1 + (G2G3 + G4)H2
H 1G2
G3G2 + G4
n2 n 1 2
C(s) = ...(2)
n 1 2 [(s n )2 ( n 1 2 )2 ]
3. Taking inverse Laplace transform on both sides of eq. (2),
n
c(t) = e nt sin [( n 1 2 )t ] ...(3)
1 2
4. For < 1,
n
c(t) = e nt sin [( n 1 2 )t ] ...(4)
1 2
5. The time response for < 1 is decaying exponential oscillations
and the output at times goes negative also.
Control System SP–11 C (EN-Sem-5)
www.aktutor.in
c(t) <1
0
t
Fig. 5.
Numerical :
K
Given : G(s) = , = 0.5, H(s) =1
s( s 10)
To Find : K, ts, tp and Mp.
C ( s) G (s )
1. =
R( s) 1 G (s) H ( s)
K
= 2 ...(5)
s 10 s K
2. For second order system,
C ( s) 2n
= 2 ...(6)
R( s) s 2 n s 2n
3. Comparing eq. (5) and (6), we get
n = K rad/sec
2 n = 10
2 × 0.5 × K = 10
K = 10
K = 100
0.5
2 2
1 1 0.5
4. Peak overshoot, Mp = e = e = 0.1630
4 4
5. Settling time, ts = = = 0.8 sec
n 0.5 10
6. Peak time, tp = = 0.36 sec.
n 1 2 10 1 (0.5)2
Ke sT
Given : G(s) H(s) =
s(s 1)
To Find : Maximum value of K.
1. Putting s = j
Ke jT
G(j) H(j) =
j( j 1)
K
2. |G(j) H(j)|=
1 2
180
G(j) H(j) = – tan–1 – 90° – T ×
= – 57.3T – 90° – tan–1
3. Intersection of Nyquist plot with negative axis of G(s) H(s) plane is
determined by using following relation,
G(j) H(j) = –180° (2k + 1)
where, k = 0, 1, 2 ....
4. It the first instant, the Nyquist plot intersects the negative real
axis of G(s) H(s) plane for k = 0, and the frequency at intersection
point is 2.
5. Therefore, G(j2) H(j2) = – 180°
– 57.32T – 90° – tan–12 = – 180°
– 57.32T – tan–12 = – 90°
6. T = 0.5
– 28.652 – tan–12 = 90°
7. Using trial and error, 2 = 1.3075 rad/sec
K K
G(j2) H(j2) = = = 0.4646 K
2 1 22 1.3075 1 (1.3075)2
8. For stability the point (–1 + j0) be placed outside the Nyquist plot
0.4646 K < 1
K < 2.152
– 1 + j0 0
Re
A r
(r – )
Fig. 6.
4. From Fig. 6, OB = G(j)
OA = – 1
AB = OB – OA = G(j) – (– 1)
AB = 1 + G(j)
5. From eq. (2)
C ( j) OB
= M() =
R( j)
AB
C( j) OB r
=
= (r – )
R( j) AB
C( j)
= M() e j ()
R( j)
where M(j) is the magnitude and () = r – .
6. Frequency response consists of two parts :
i. Magnitude,
ii. Phase angle.
iii. Both can be plotted against different values of .
Solved Paper (2016-17) SP–14 C (EN-Sem-5)
www.aktutor.in
7. Hence frequency response of closed loop system is plot of
magnitude and phase angle.
Correlation :
1. For 2nd order system, the transfer function is,
C ( s) 2n
= 2
R( s) s 2 n s 2n
where, = Damping factor
n = Natural frequency of oscillations
+ n2
R(s) C(s)
2
– s + 2 ns
Fig. 7.
2. Closed loop frequency response is,
C( j) 2n
= T(j) =
R( j) ( j) 2 n ( j) n2
2
2n 1
= = 2
2 2 n ( j ) 2n
1 2 j
n n
1
= ...(3)
(1 u2 ) j 2u
where u = /n, normalized driving frequency.
1
|T(j)| = M = ...(4)
(1 u ) (2u)2
2 2
2u
and T(j) = = tan 1 ...(5)
1 u2
3. The steady state output is
1 2u
c(t) = sin t tan 1
(1 u ) (2u) 2 2 12 u2
From eq. (4) and (5) when
u = 0, M = 1 and = 0
1
u = 1, M = and = –
2 2
u = , M = 0 and =–
4. The frequency where M has a peak value is called the resonant
frequency. At this frequency the slope of magnitude curve is zero.
If r = Resonant frequency.
ur = r/n is normalized resonant frequency.
dM 1 [ 4(1 ur2 )ur 8 2ur ]
= 0
du u ur 2 [(1 ur2 )2 (2ur )2 ]3/ 2
Control System SP–15 C (EN-Sem-5)
www.aktutor.in
– 4(1 – ur2) ur + 82ur = 0
– 4ur(1 – ur2 – 22) = 0
ur = 1 2 2
r = n 1 2 2 ...(6)
1
Mr = ...(7)
2 1 2
5. The phase angle of T(j) at resonant frequency is
1 2 2
r = tan 1
1.0
Magnitude
M
Mr
M( = 0.707)
0
ur = r / n
Normalized frequency, u
Fig. 8. Frequency response magnitude characteristics.
0°
( = 1)
– 90°
Phase
angle
– 180°
0 1.0
Normalized frequency u
Fig. 9. Frequency response phase characteristic.
Solved Paper (2016-17) SP–16 C (EN-Sem-5)
www.aktutor.in
Bandwidth :
1
1. The frequency at which M has value of and is called cut-off
2
frequency c. The signal frequencies above cut-off are attenuated.
1
2. The range of frequencies for which M is known as bandwidth
2
b. The low pass filters has bandwidth equal to cut-off. b indicates
the noise filtering characteristics of the system.
M
Mr
1
0.707
Bandwidth
r c
Fig. 10.
3. Normalized bandwidth, ub = b/n
1 1
M=
(1 ub2 )2 (2ub )2 2
1
ub = [1 2 2 2 4 2 4 4 ]2
1
b = [1 2 2 2 4 2 4 4 ]2
4. The damped frequency of oscillations d and peak overshoot Mp of
the step response for 0 1 are
d = n 1 2
Mp = exp
2
1
1
5. For > , Mr does not exists so the correlation breaks down
2
r 1 2 2
=
d 1 2
R2
ei eo
C
0º
– 45º m
– 90º
log
1 m 1
T T
Fig. 12. Bode plot of lag compensator.
SECTION-C
G10
G9 G7
X2 G8
X3 X4 X5 X6
X1 X7
G1 G2 G3 G4 G5 G6
–H 1 – H2 –H3
H4
Fig. 13.
Ans.
1. Forward path and gains are :
P1 = G1 G2 G3 G4 G5 G6, 1 = 1
P2 = G1 G9 G4 G5 G6, 2 = 1 – G7
P3 = G1 G9 G10 G6, 3 = 1 – G7 – G8 + G7 G8
P4 = G1 G2 G3 G10 G6, 4 = 1 – G8
2. Loops and gains are :
L1 = – G2 H1
L2 = – G4 H2
L3 = – G5 H3
L4 = – G3 G4 G5 H4
L5 = G 7
L6 = G 8
L7 = H1 G9 G10 H4
L8 = H1 G9 G4 G5 H4
Solved Paper (2016-17) SP–20 C (EN-Sem-5)
www.aktutor.in
L9 = – H4 G3 G10
L10 = H2 G10 H3
= 1 – (L1 + L2 + L3 + L4 --------- + L10) + (L1 L2
+ L1 L3 + L1 L6 + L1 L10 + L7 L6 + L5 L6 + L5 L2
+ L5 L3 + L5 L10 + L9 L5 + L9 L3)
4
X7
Pk k
k 1
3. Using Mason’s gain formula, T(s) = =
X1
G1G2G3G4G5G6 G1G9G4 G5G6 (1 G7 ) G1G9G10G6 (1 G7 G8 G7G8 )
=
1 (G2 H1 G4 H2 G5 H3 G3G4 G5 H4 G7 G8 H1G9G10 H4
H1G9 G5 H4G4 H4 G3G10 H2G10 H3 ) (G2 H1G4 H2 G2 H1G5 H3
G2 H1G8 G2 H1 H2G10 H3 H1G9G4 G5 H4 G8 G7G8 G7G4 H2
G7G5 H3 G7 H2G10 H3 H4G3G10G7 H4 G3G10G5 H3 )
Fig. 14.
Ans. Sensitivity : Sensitivity is the change in variable due to variation
in parameters of control system.
e.g.
T
100
% Change in T
T
S = T
G
% Change in G G
100
G
Numerical :
K
Given : G(s) = , H(s) = 1
s(s 5)
K ( s2 5s K ) 1 K
=
K ( s 2 5s K ) (s 2 5 s K ) 2
s 2 5s K K s(s 5)
= ( s 2 5s K ) 2 2
= 2 2
( s 5 s K ) ( s 5s K )
4. Putting s = j, K = 1
j( j 5)
=
2
( j)2 5 j 1
( 25 j 25)
=
( j 25 24)2
252 252
| SKT | = = 1.02
252 24 2
10 20 30 40
Fig. 15.
Ans. Types of test signal :
1. Unit step : Signals which start at time t = 0 and have magnitude of
unity are called unit step signals.
They are represented by a unit step function u(t).
Solved Paper (2016-17) SP–22 C (EN-Sem-5)
www.aktutor.in
u(t)
t
0
Fig. 16. Unit step.
They are defined mathematically as :
1 ; t 0
u(t) =
0 ; t 0
2. Unit ramp : Signals which start from zero and are linear in nature
with a constant slope m are called unit ramp signals.
They are represented by a unit ramp function r(t).
They are defined mathematically as :
mt ; t 0
r(t) =
0 ; t0
r(t)
r(t) = mt
m = tan
t
0
Fig. 17.
3. Unit impulse : Signals which act for very small time but have
large amplitude are called unit impulse functions.
They are represented by (t).
They are defined mathematically as,
0; t0
(t) =
1 ; t0
(t) dt =1
x(t)
(t)
–t t
Fig. 18.
Control System SP–23 C (EN-Sem-5)
www.aktutor.in
4. Unit Parabolic Signal : The continuous-time unit parabolic
function p(t), also called acceleration signal starts at t = 0, and is
defined as :
t2
; for t 0
p(t) = 2
0 ; for t 0
t2
or p(t) = u(t)
2
p(t)
0 t
Fig. 19.
Relation :
1. Relation between impulse and step signal :
d
(t) = u(t)
dt
2. Relation between step and ramp signal :
d
u(t) = r (t )
dt
3. Relation between ramp and parabolic signal :
d t2
r(t) = .
dt 2
Numerical :
400( s 1)
Given : G(s) = , H(s) = 1
( s 2)( s 8)
To Find : Steady state error.
1. Laplace of periodic waveform using formula :
R1 (s)
R(s) =
1 e TS
2. First we have to find Laplace of r(t),
r(t)
t
10
Fig. 20.
By integration method
Solved Paper (2016-17) SP–24 C (EN-Sem-5)
www.aktutor.in
10
R1(s) = 2 e st dt
0
2 10
= e st
( s) 0
2
R1(s) = (1 e 10 s )
s
R1 (s)
4. Here, T = 20, So R(s) =
1 e TS
2
(1 e10 s )
R(s) = s
(1 e20 s )
2
(1 e 10 s )
R(s) = s
(1 e10 s )(1 e10 s )
2
R(s) =
s(1 e10 s )
5. Using formula
sR( s)
ess = lim
s0 1 G( s) H (s)
2
s
s (1 e10 s )
= lim
s0 400( s 1)
1
( s 2)( s 8)
2
= lim (1 e10 s )
s 0 ( s 2)( s 8) 400( s 1)
( s 2)( s 8)
2(s 2)( s 8)
= lim
s 0 [( s 2)( s 8) 400(s 1)](1 e10 s )
2(0 2)(0 8)
=
[(0 2)(0 8) 400(0 1)](1 e0 )
2 28 32
=
2 8 400 416
ess = 0.0769
SB
SD D B
C
Stator
Vd SC
b. Sketch the root locus for the closed loop control system
K
with G(s) =
s(s 1)( s2 4 s 5)
Ans.
K
Given : G(s) =
s( s 1)(s 2 4 s 5)
To Draw : Root locus.
1. Poles i.e., s = 0, s = – 1, s = – 2 + j and s = – 2 – j
2. There is no open loop zero.
3. Number of poles, P= 4
Number of zeros, Z= 0
P – Z = 4 – 0 = 4 i.e., four branches of root locus terminates at
infinity.
4. Angle of asymptotes :
(2 q 1)
q = 180
PZ
where, q = 0, 1, 2, ...... upto (P – Z – 1)
201
0 = 180 = 45°
4 0
2 11
1 = 180 = 135°
40
221
2 = 180 = 225°
4 0
231
3 = 180 = 315°
40
5. Centroid of asymptotes :
Real parts of poles – Real parts of zeros
=
PZ
(0 1 2 j 2 j ) (0) 5
= = = – 1.25
40 4
Solved Paper (2016-17) SP–28 C (EN-Sem-5)
www.aktutor.in
6. Breakaway points : The characteristic equation is 1 + G(s) H(s)
= 0
s(s + 1) (s2 + 4s + 5) + K = 0
s4 + 5s3 + 9s2 + 5s + K = 0
dK / ds = 4s3 + (5 × 3s2) + (9 × 2s) + 5 × 1 = 0
4s3 + 15s2 + 18s + 5 = 0
(s + 0.4) (4s2 + 13.4s + 12.5) = 0
(s + 0.4) {(s + 1.675 – j 0.565) (s + 1.675 + j 0.565)} = 0
Therefore three breakaway points are obtained
s = – 0.4 on real axis and
s = – 1.7 + j 0.6 and s = – 1.7 – j 0.6
7. Intersection points with imaginary axis :
The characteristic equation is s4 + 5s3 + 9s2 + 5s + K = 0
The Routh array :
s4 1 9 K
s3 5 5 0
s2 8 K
5K
s1 5 0
8
s0 K
Value of K at imaginary axis :
5K
5 = 0; K= 8
8
Solving auxiliary equation formed from the s2 terms in Routh array,
therefore
8s2 + K = 0
8s2 + 8 = 0
s2 + 1 = 0
s=± j
8. Angle of departure from complex pole :
d = 180° – (P – Z)
P1 = 180° – tan–1(1/2) = 153.43°
P2 = 180° – tan–1(1/1) = 135°
P3 = 90°
(–2 + j) = 180° – (P1 + P2 + P3)
(–2 + j) = 180°– (153.43° + 135° + 90°) = – 198.43°
Control System SP–29 C (EN-Sem-5)
www.aktutor.in
and – j) = + 198.43°
(– 2
Im
–2 + j
j1
13
5° 153.43°
–2 –1 Re
90°
– j1
p = 90 °
3
Fig. 22.
9. Root locus :
K= K=
Im
–2+j
j1
– 198.43°
–1
Re
+ 198.43° –2
– j1
–2–j
K= K=
1 27u
Given : c(t) = sin t tan 1
2 2 2 2
(1 u ) 4 z u 1 u2
To Find : i. Resonant frequency
ii. Resonant peak
iii. Bandwidth.
iv. Phase angle
1
1. M=
(1 u ) 4 z2u2
2 2
27u
2. = tan 1
1 u2
3. The frequency where M has a peak value is known as resonant
frequency. At this frequency, the slope of the magnitude curve is
r
zero. Let r be the resonant frequency and ur =
n
r = n 1 2 z2
4. Maximum value of magnitude is known as resonant peak is given
by
1 1
Mr = 2 2 2 2
=
(1 u ) 4 z u
r r (1 1 2 z ) 4 z2 (1 2 z2 )
2 2
1 1 1
= = =
4 2 2 2 2
4 z 4 z (1 2 z ) 2z z 1 2z 2 z 1 z2
5. The range of frequencies over which M is equal to or greater than
1 b
: Putting ub =
2 n
1 1
M=
(1 ub2 )2 4 z2ub2 2
(1 – ub2)2 + 4z2ub2 = 2
Control System SP–31 C (EN-Sem-5)
www.aktutor.in
1 + ub4 – 2ub2 + 4z2ub2 – 2 = 0
ub = [1 2 z2 2 4 z2 4 z4 ]1/ 2
b = n [1 2 z2 1 4 z2 4 z4 ]1/ 2
b. Draw Bode plot (log magnitude plot) for the transfer
function.
20 s
G(s) =
s2 20 s (100)2
OR
Using Nyquist stability criterion, investigate the stability
of a unity feedback system with open loop transfer function.
( s z1 )
G(s) = , z1, p1 > 0
s( s p1 )
Also discuss the significance of M-circle.
Ans. Bode plot :
20 s
1. G(s) = 2
20 s s
(100)2 1 2
(100) 100
K = 20 / (100)2 = 0.002
1 2 3 4 5 6 7 8 9 101 2 3 4 5 6 7 8 9 102 2 3 4 5 6 7 8 9 103
Gain
in – 80°
40 dB
dB
– 20 dB/dec – 100°
Bode plot :
20 dB/dec
Solved Paper (2016-17)
– 120°
– 140°
10 dB
– 160°
0 dB – 180°
– 200°
g = 1 rad/sec g = 100 rad/sec
– 10 dB
– 220°
Fig. 24.
– 240°
– 20 dB – 260°
www.aktutor.in
– 280°
j z1 2 z12
1. G(j) H(j) =
j( j p1 ) 2 p12
1 ( p1 z1 )
= 90 tan 2
( p1 z1 )
2. Now,
lim G ( j) H ( j) = + 90°
0
1
G( j p1 z1 ) H ( j p1 z1 ) = 0
p1 p1 z1
3. Thus the locus comes down in the first quadrant, crosses the positive
real axis into the fourth quadrant, and approaches the origin from
an angle of – 90°.
4. Path def maps into the origin, and ija maps into on semicircle at
infinity. The resulting plot is shown in Fig. 25.
Im GH
a
GH(j )
= ± p1z1
j d , e, f Re GH
i
Fig. 25.
Significance of M-circle :
M-circle can be used to obtain closed loop frequency response from
open loop frequency response.
Solved Paper (2016-17) SP–34 C (EN-Sem-5)
www.aktutor.in
M = 1.6 Y
M=1
M = 0.5
M=3
3 1 1
The state equation for a system is, x = x u
2 0 1
Ans. Kalman’s test for controllability :
1. Consider nth order multiple input linear time invariant system
represented by its state equation as,
x = Ax + Bu
2. The necessary and sufficient condition for the system to be
completely state controllable is that the rank of the composite
matrix Qc is n.
3. The composite matrix Qc is given by,
Qc = [B : AB : A 2B : ... A n–1B]
Gilbert’s test for controllability :
1. For the Gilbert’s test it is necessary that the matrix A must be in
canonical form. Hence the given state model is required to be
transformed to the canonical form first, to apply the Gilbert’s test.
2. Consider single input linear time invariant system represented
by,
z = Ax+Bu
where A is not in the canonical form. Then it can be transformed
to the canonical form by the transformation,
x = Mz
where M = Model matrix
3. The transformed state model,
z = A z B u
where A = M – 1 AM
B = M – 1 B
Control System SP–35 C (EN-Sem-5)
www.aktutor.in
4. In such a case the necessary and sufficient condition for the
complete state controllability is that the vector matrix B should
not have any zero elements. If it has zero elements then the
corresponding state variables are not controllable.
5. If the eigen values are repeated then matrix A cannot be
transformed to Jordan canonical form. If A has eigen values 1, 1,
2, 2, 3, 4, ... n then the transformation results Jordan canonical
form shown in matrix below
1 1 0 0 0 0 ..........0
J= 0 1 0 0 0 0...........0
0 0 2 1 0 0...........0
Jordan
0 0 0 2 1 0...........0
block
0 0 0 0 2 0...........0
0 0 0 0 0 3..........0
........
........
x=Ax+Bu
and y=Cx
2. The system is completely observable if and only if the rank of the
composite matrix Qo is n.
Composite matrix Qo is given by,
Qo = [CT : AT CT :.......... (AT)n–TCT]
where CT = Transpose of matrix C
AT = Transpose of matrix A
Gilbert’s test for observability :
1. For Gilbert’s test, the state model must be expressed in the
canonical form. Consider the state model of linear time invariant
system as,
Solved Paper (2016-17) SP–36 C (EN-Sem-5)
www.aktutor.in
x=Ax+Bu
and y=Cx
2. Use the transformation x = Mz where M is the model matrix.
y = CM z = C z
where C = CM
For a single input single output system,
z1
z
C
y = Cz C 12 .... C 1 n 2
11
zn
= C 11 z1 C 12 z1 .... C 1n zn
3 1 1
Given : x = x u
2 0 1
To Find : Controllability.
3 1 1
1. A= , B = 1 , n = 2
2 0
2. Qc = B AB
3 1 1 2
3. AB = 1 = 2
2 0
1 2
4. Qc =
1 2
1 2
5. |Qc| = =2–2=0
1 2
Rank of Qc is r n
Hence, system is not completely controllable.
Control System SP–37 C (EN-Sem-5)
www.aktutor.in
b. Design a phase lead compensator for a negative unity
feedback system with plant transfer function.
K
Gp (s) = to satisfy the conditions :
s(s 10)( s 1000)
phase margin in atleast 45°, static error constant = 1000 s–1
Ans.
K
Given : GP(s) = , Phase margin = 45°
s(s 10) (s 1000)
Kv = 1000 sec – 1
To Design : Phase lead compensator.
K
1. Kv = lim sG p (s) = lim s
s0 s 0 s( s 10) (s 1000)
K
Kv =
10 1000
K = 10 × 1000 Kv = 107
107 103
2. Gp(s) = =
s s s s
10000 s 1 1 s 1 1
10 1000 10 1000
1
Gp(j) = – 90 – tan – 1 tan
10 1000
1
= m = ( 0.1324) × 165 = 60 rad/sec
T
8. The upper corner frequency
1 1
= × 60 = 453.1722 rad/sec
T 0.1324
1 jT
9. Lead compensation network, Gc(s) = Gc(j) =
1 jT
1
1 j
60 (1 j 0.0166)
= =
1 (1 j 0.022066)
1 j
453.1722
10. So, the open loop transfer function of the compensated system is
1000 (1 j 0.0166)
G(s) = Gp(s) Gc(s) =
s(s 10) (s 1000) (1 j 0.022066)
1000 (1 0.0166 s)
=
s (s 10)( s 1000)(1 0.22066 s)
1 2 3 4 5 6 7 8 9 101 2 3 4 5 6 7 8 9 102 2 3 4 5 6 7 8 9 103
Control System
Gain curve
– 80°
– 100°
Bode plot :
– 120°
– 40 dB/dec Angle
– 20 dB/dec – 140°
in
degree – 160°
p = g = 100 rad/sec
– 180°
– 200°
Gain in dB
– 220°
Fig. 27.
Phase curve
– 240°
– 260°
www.aktutor.in
– 280°
– 60 dB/dec
Note : Attempt all sections. If any missing data is required, then choose
suitably.
SECTION-A
e1(t) C1 R2 e2(t)
Fig. 1.
f. Draw the block diagram and explain the open loop control
system and closed loop control system.
SECTION-B
2(s + 1)
1
Fig. 2.
SECTION-C
Control System SP–5 C (EN-Sem-5)
www.aktutor.in
SOLUTION OF PAPER (2017-18)
Note : Attempt all sections. If any missing data is required, then choose
suitably.
SECTION-A
Fig. 1.
In open loop control systems the control action is independent of
the desired output. In this system the output is not compared with
the reference input.
Example : Washing machine, Immersion rod, Time operated
traffic control, DC shunt motor.
2. Closed loop control system :
Error signal
Reference Controlled
Amplifier Controlled process
input + – output
+
Feedback
Fig. 2.
In a closed loop control system the output has an effect on control
action through a feedback as shown in Fig. 2.
Example : Automatic steering control system, Driving system of
an automobile, Home heating system, Ship stabilization system.
e1(t) C1 R2 e2(t)
Fig. 3.
Ans.
1. Apply KVL in the loop,
1
– Vi(s) + R1 || I(s) + R2 I(s) = 0
C1 s
1
R1 C s
1
Vi(s) = R2 I (s)
R1 1
C1 s
R1
+
1/C1s R2
Vi(s) Vo(s)
I(s)
–
Fig. 4.
2. And, V0(s) = R2 I(s)
Control System SP–7 C (EN-Sem-5)
www.aktutor.in
R2 Vi ( s) R2 Vi ( s)
V0(s) =
1 R1
R1 C s R C s 1 R2
1
R2 1 1
R 1
1
C1 s
V ( s) R2 ( R1 C1 s 1)
Transfer function, 0 = .
Vi (s) R2 ( R1 C1 s 1) R1
0 = =0
Re
= – tan–1 T
a
K
2 + a2
G(j ) – plane
Fig. 5.
4. Diameter = |G(j)|=0 – |G(j)|=
K K
= –0
a a
So, Radius = K / 2a
G ( j) 0 – G( j)
5. Center = = K/2a.
2
f. Draw the block diagram and explain the open loop control
system and closed loop control system.
Solved Paper (2017-18) SP–8 C (EN-Sem-5)
www.aktutor.in
Ans. There are two types of control system :
1. Open loop control system :
Fig. 6.
In open loop control systems the control action is independent of
the desired output. In this system the output is not compared with
the reference input.
Example : Washing machine, Immersion rod, Time operated
traffic control, DC shunt motor.
2. Closed loop control system :
Error signal
Reference Controlled
Amplifier Controlled process
input + – output
+
Feedback
Fig. 7.
In a closed loop control system the output has an effect on control
action through a feedback as shown in Fig. 7.
Example : Automatic steering control system, Driving system of
an automobile, Home heating system, Ship stabilization system.
2. Parallel connection :
+ C(s)
R(s) G1(s)
± R(s) C(s)
G1(s) ± G2(s)
G2(s)
SECTION-B
2(s + 1)
1
Fig. 8.
Ans. Step 1 :
+ 20
R(s) 20 2 C(s)
– s 2 +47s + 50
Step 2 :
+ 800
R(s) C(s)
– s 2 +47s + 50
800
So, G(s) = and H(s) = 1
s 2 47 s 50
800
G(s) H(s) =
s0 (s 2 47 s 50)
Type of system is zero.
i. r(t) = 6
6 1
So, R(s) = A = 6, which is coefficient of
s s
800
Now, Kp = lim G( s) H ( s) =
s0 50
A 6 6 50
Now, ess = = 0.353
1 K p 1 800 850
50
ii. r(t) = 8 t
8 1
So, R(s) = A = 8, which is coefficient of 2
s2 s
Now, Kv = lim sG( s) H ( s)
s0
s 800
= lim =0
s0 s 2 47 s 50
Control System SP–11 C (EN-Sem-5)
www.aktutor.in
A 8
Now, ess = =
Kv 0
iii. r(t) = 10 + 4t + 15t2
t2
= 10 + 4t + 30
2
Now, So, A1 = 10, A2 = 4, A3 = 30
K = lim G (s) H (s)
p s0
800
= = 0.941
850
Kv = lim
s0
sG( s) H (s)
=0
Ka = lim s2 G(s) H (s)
s 0
=0
A1 A A
Now, ess = 2 3
1 K p Kv Ka
10 4 30
= =
1 0.941 0 0
0 0 0 40
Given : A = 1 0 3 , B = 10 , C = [0 0 1]
0 0 4 0
To Test : Controllability and observability.
1. Controllability test :
0 0
i. AB = 40 , A 2B = 30
0 0
ii. Controllability test matrix is given by,
Solved Paper (2017-18) SP–12 C (EN-Sem-5)
www.aktutor.in
40 0 0
Qc = [B : AB : A B] = 10 40 0
2
0 0 0
|Qc| = 0
iii. Thus the rank of Qc is 2. Hence the system is not controllable.
2. Obervability test :
0 1 0 0
T
i. A = 0
0 0 , CT = 0
0 3 4 1
0 0
ii. ATCT = 0 , (AT)2 CT = 0
4 16
iii. The observability test matrix is given by
Qo = [CT : ATCT : (AT)2 CT]
0 0 0
Qo = 0 0 0
1 4 16
iv. Its rank is 1. Hence the system is not observable.
8.24
–2 1 4 Re
– 0.24
35(s 4)
= lim 1=
s0 s( s 2)(s 5)
Kv = lim
s0
s G(s) H ( s)
Control System SP–15 C (EN-Sem-5)
www.aktutor.in
s 35 ( s 4)
= lim 1
s0 s (s 2) (s 5)
35 4
= = 14
10
Ka = lim s 2 G( s) H ( s)
s0
35( s 4)
= lim s2 1 = 0
s( s 2)(s 5)
s0
K=0
–9 –4 –1 Re
K=
Fig. 10. Root locus.
Control System SP–19 C (EN-Sem-5)
www.aktutor.in
6. Attempt any one part of the following : (10 × 1 = 10)
a. Sketch the RL (root locus) for a unity feedback system
K (s 2 + 2s + 10)
with OLTF G(s) =
s2 + 4 s + 5
Ans.
K ( s2 2s 10)
Given : G(s) = , H(s) = 1
s2 4 s 5
To Draw : Root locus.
1. The open loops poles are at s = – 2 – j, s = – 2 + j
2. The open loop zeros are at s = – 1 – j3, s = – 1 + j3
3. Number of poles = 2
Number of zeros = 2
P – Z = 2 – 2 = 0 i.e., none of branches of root locus terminates at
infinity.
4. Angle of departure and arrival can’t be calculated because the
centroid, = .
Im
j3
j1
Re
–2
– j1
– j3
Fig. 11.
K ( s 2)
Given : G(s) = , r(t) = t2u(t)
s(s3 7 s 2 42 s)
To Find : Ka and ess.
t2
1. r(t) = t2u(t) = 2 u( t)
2
A = 2. [Coeff. of t2/2]
K ( s 2)
2. Now, G(s)H(s) =
s(s3 7 s 2 42 s)
K ( s 2)
=
s 2 ( s2 7 s 42)
3. Now, Ka = lim s 2 G( s) H(s)
s0
s2 K (s 2)
= lim
s 0 s ( s2 7 s 42)
2
2K K
=
42 21
A 2 42
4. ess =
Ka K K
21
10 e s
i. Given : G(s) =
s1
To Sketch : Polar plot.
1. Putting s = j
10 e j
G(j) =
j 1
2. Magnitude :
|10 e j | 10
|G(j)| =
2 1 2 1
Control System SP–21 C (EN-Sem-5)
www.aktutor.in
3. Phase angle :
G(j) = – – tan – 1
4. Here the magnitude decreases from unity monotonically and the
phase angle also decreases monotonically and indefinitely, the polar
plot of the given transfer function is spiral as shown in Fig. 12.
Im
O 1
Re
180°
– 90°
Fig. 12.
32
ii. Given : G(s) =
(s 4)(s 2 4 s 8)
To Sketch : Polar plot.
1. Putting s = j
32
G(j) =
( j 4)[ j 2 4 j 8]
32
=
( j 4)[(8 2 ) 4 j ]
32(4 j)[(8 2 ) 4 j]
=
(16 2 )[(8 2 )2 162 ]
32[4(8 2 ) 42 j{(8 w2 ) 16}]
=
(16 2 )[(8 2 ) 2 162 ]
2. Saparating into real and imaginary parts
32[4(8 2 ) 4 2 ] j32[ (8 2 ) 16]
G(j) =
(16 )[(8 ) 16 ] (16 2 )[(8 2 )2 162 ]
2 2 2 2
4
G(j) = tan 1 tan 1
4 8 2
6. The value of G(j) at = 0, and =
32
|G(j0)| = =1
48
|G(j)| = 0
G(j0) = 0°
G(j) = – 270,
Im
=+2
=+ = 4.899
Re
=–
– 4.899
=0
=–2
Fig. 13.
Control System SP–1 C (EN-Sem-5)
www.aktutor.in
B.Tech.
(SEM. V) ODD SEMESTER THEORY
EXAMINATION, 2018-19
CONTROL SYSTEM
Time : 3 Hours Max Marks : 70
Note : Attempt all sections. If any missing data is required, then choose
suitably.
SECTION-A
1. Attempt all parts of the following : (2 × 7 = 14)
a. What is Mason’s gain formula ?
SECTION-B
c. Draw the Nyquist plot for the unity feedback system whose
open loop transfer function is
K
G(s) H(s) =
s2 (1 + sT )
Solved Paper (2018-19) SP–2 C (EN-Sem-5)
www.aktutor.in
d. Determine the transfer function of the circuit given below :
G2
+ + + + +
R G1 G3 G4 C
– – –
H1 H2
H3
Fig. 1.
SECTION-C
Solved Paper (2018-19) SP–4 C (EN-Sem-5)
www.aktutor.in
SOLUTION OF PAPER (2018-19)
Note : Attempt all sections. If any missing data is required, then choose
suitably.
SECTION-A
It is denoted by ess.
SECTION-B
1
P2 = 180° – tan–1 = 135°
1
P3 = 90°
(–4 + j3) = 180° – (126.87° + 135° + 90°)
= – 171.87
Root locus :
K= d = –171.87° j
K=0
4
(– 3 + j4)
K=
3 = 90°
K=
K=0
4
d = 171.87°
K=
Fig. 1.
Reference winding
Vr
Fig. 2.
2. Field controlled DC motor (DC Servomotor) transfer
function :
1. A schematic diagram of a field controlled DC motor (DC servo
motor) shown in Fig. 3.
Control System SP–9 C (EN-Sem-5)
www.aktutor.in
Rf Ra
ia = Constant
ef if Lf J ea
T
B
Fig. 3.
2. Here, Rf = Field winding resistance
Lf = Field winding inductance
If = Field winding current
Ra = Armature resistance
ia = Armature current
= Angular displacement.
3. The torque T developed by the motor is proportional to product of
the air-gap flux and armature current ia so we get
T = K1 i a ...(1)
where K1 is constant.
4. But the air gap flux and the field current if are proportional for the
usual operating range of the motor and ia is assumed to be constant,
we can rewrite the above equation as
T = K2 i f ...(2)
where K2 is a constant.
5. The equations for this system are
dif
Lf R f if = ef ...(3)
dt
Jd 2 Bd
and = T = K2 i f ...(4)
dt 2 dt
6. By taking the Laplace transform on both sides of eq. (3) and (4)
where all initial conditions are zero, we get
(Lfs + Rf) if(s) = Ef(s) ...(5)
2
(Js + Bs) (s) = K2if(s) ...(6)
7. From the above equations, the transfer function of this system is
obtained as
(s) K2
= ...(7)
E f (s ) s( Lf s Rf )( Js B)
8. Block diagram is shown in Fig. 4.
Solved Paper (2018-19) SP–10 C (EN-Sem-5)
www.aktutor.in
1 If(s) K2
Ef(s) (s)
L fs + R f s(Js + B)
Fig. 4.
c. Draw the Nyquist plot for the unity feedback system whose
open loop transfer function is
K
G(s) H(s) =
s2 (1 + sT )
Ans.
1. Putting s = j
K
G(j) H(j) =
( j)2 (1 jT )
K
2. |G(j) H(j)| =
2 1 2T 2
3. G(j) H(j) = – 180° – tan–1 T
4. At = 0, G(0)H(0) = – 180°
At = , G()H() = 0 – 270° = 0 90°
Imj
=0 =
( –1,0) 0 Re
Fig. 5.
Fig. 6.
Control System SP–11 C (EN-Sem-5)
www.aktutor.in
Ans. Step 1 :
G 2/G 1
+
G1 + + a C
G3 G4
R 1 + G 1H 1 – – +
H2
H3
H 3/G 3
H2
+
G1 + G 3 G4
R C
1 + G1H1 – 1 + G 3H 2 +
H3/G4
+ +
G1 + G 3G 4
C
1 + G 1H 1 – 1 + G 3H 2
H 3/G 4
G1 G3 G4 + G2 + G2 G3 H2
R C
1 + G 1H 1 + G 3 H2 + G 4 H 3 + G1 G3 H1 H2 + G1 G4 H 1 H 3
+ n2
R(s) C(s)
2
– s + 2 ns
Fig. 7.
2. Closed loop frequency response is,
C( j) 2n
= T(j) =
R( j) ( j) 2 n ( j) n2
2
2n 1
= = 2
2 2 n ( j ) 2n
1 2 j
n n
1
= ...(1)
(1 u2 ) j 2u
where u = /n, normalized driving frequency.
1
|T(j)| = M = ...(2)
(1 u2 )2 (2u)2
2u
and T(j) = = tan 1 ...(3)
1 u2
3. The steady state output is
1 2u
c(t) = sin t tan 1
(1 u ) (2u) 2 2 1 u2
2
ur = 1 2 2
r = n 1 2 2 ...(4)
1
Mr = ...(5)
2 1 2
5. The phase angle of T(j) at resonant frequency is
1 2 2
r = tan 1
1.0
Magnitude
M
Mr
M( = 0.707)
0
ur = r / n
Normalized frequency, u
Fig. 8. Frequency response magnitude characteristics.
0°
( = 1)
– 90°
Phase
angle
– 180°
0 1.0
Normalized frequency u
dFe ( s) d s 2 2s
5. C1 = lim lim
s0 ds s 0 ds (s 1)2
( s 1)2 (2s 2) ( s2 2s) (2) (s 1)(1) 2
= lim = lim =2
s0 ( s 1) 4 s 0 (s 1)3
6. C2 = lim
s0
d 2 Fe (s)
ds 2
lim
s 0 ds
d dFe ( s)
ds
lim
d 2
s 0 ds ( s 1)3
3 2
( s 1) (0) 2(3) (s 1) (1) 6
= lim = lim =–6
s0 (s 1)6 s0 (s 1)4
d 3 Fe (s) d d 2 Fe ( s) d 6
7. C3 = lim lim = lim
s 0 ds
4
s0 ds 3 s 0 ds ds 2 ( s 1)
( s 1)4 (0) 6 4(s 1)3 (1) 24
= lim = lim 5 = 24
s0
( s 1)8 s 0 (s 1)
X =
O Re
=0
= 0.35
rad/sec
Fig. 10.
7. Polar plot makes an intercept on the negative real axis at point X.
The intercept OX is calculate as given below :
1 4 j
G(j) H(j) =
( j)2 (1 j)(1 2 j)
Separating into real and imaginary parts
1 10 2 j( 83 )
G(j) H(j) = 2 2 2 2
2 2 2 2
[(1 2 ) 9 ] [(1 2 ) 9 ]
Equating imaginary part to zero
( – 83) = 0
83 = 1
= 0.35 rad/sec
i.e., At = 0.35 rad/sec,
plot crosses negative real axis.
H(s)
Fig. 11.
2n
G(s)H(s) =
s( s 2 n )
C ( s) 2n
= 2
R( s) s 2 n s 2n
For this system damping ratio is and natural frequency n.
5. And for steady state error
Kp = lim G( s) H ( s)
s0
n
and Kv = lim s G(s) H (s)
s 0 2
Advantages :
1. Steady state error is reduced.
2. Disturbance signal rejection occurs.
3. Relative stability is improved.
2. Proportional integral controller (PI) :
1. A controller in the forward path, which changes the controller
output corresponding to the proportional plus integral of the error
signal is called PI controller.
2. Output of controller = K e(t) + Ki e(t) dt
Ki K
Taking Laplace = K E( s) E( s) E( s) K i
s s
+ n2
R(s) K C(s)
+ s (s + 2 n)
–
Ki/s
Fig. 12.
3. Assume, K = 1, we can write,
Ki 2
1 s n ( K s) n2
G(s) = 2 i
s (s 2 n ) s (s 2 n )
i.e., system becomes TYPE 2 in nature.
C ( s) ( K i s) n2
= 3
R( s) s 2 n s2 s 2n K i n2
Control System SP–17 C (EN-Sem-5)
www.aktutor.in
i.e., it becomes third order.
4. Now as order increases by one, system relatively becomes less
stable as Ki must be designed in such a way that system will remain
in stable condition. Second order system is always stable.
while Kp = lim G( s) H ( s) , ess = 0
s0
Ki/s
Fig. 13.
7. The output of PID controller is time domain is,
de(t)
Controller output = Ke(t) K i e(t) dt Td
dt
8. Taking the Laplace transform, controller output in s-domain is,
K
Controller output = E(s) K i s Td
s
Solved Paper (2018-19) SP–18 C (EN-Sem-5)
www.aktutor.in
Advantages :
1. Flexible and reliable.
2. Improves stability.
Fig. 14.
The closed loop transfer function is
C ( s) 2n / s( s 2 n )
=
R( s) 1 2n / s( s 2 n )
where, = Damping factor or Damping ratio
C ( s) 2n
= 2
R( s) ( s 2 n s 2n )
2n
Then output, C (s) = R s ...(1)
s 2 n s n2
2
Put d = n 1 2
1 s 2 n
C (s) = ...(4)
s (s n )2 d2
4. Rewrite eq. (4)
1 s n n d
C(s) = ...(5)
s (s n )2 d2 d (s n )2 2d
5. Taking inverse Laplace transform of eq. (5),
t n nt
c(t) = 1 e n cos d t e sin d t
d
Control System SP–19 C (EN-Sem-5)
www.aktutor.in
Put d = n 1 2
ent
c(t) = 1 [( 1 2 )cos d t sin d t ] ...(6)
1 2
6. Put sin = 1 2 ,
cos =
e n t
c(t) = 1 (sin cos d t cos sin d t)
1 2
e n t
c(t) = 1 sin( d t ) ...(7)
1 2
where d = n 1 2
1 2
and = tan 1
7. Eq. (7) is rewritten as
1 2
e n t
c(t) = 1 sin ( n 1 2 )t tan 1 ...(8)
1 2
8. The term n is called natural frequency of oscillations. Term
c(t)
<1
t
Fig. 15. Underdamped oscillations.
5. Attempt any one part of the following : (7 × 1 = 7)
a. Construct the state model of a system characterized by
the differential equation. Give the block diagram
representation of the state model.
d3 y d2 y dy
3
+ 6 2 + 11 + 6y = u
dt dt dt
Solved Paper (2018-19) SP–20 C (EN-Sem-5)
www.aktutor.in
Ans.
d3 y d2 y dy
1. 3
+ 6 2 + 11 + 6y = u
dt dt dt
2. Taking Laplace on both sides : Y(s)(s3 + 6s2 + 11s + 6) = U(s)
Y(s)/U(s) = 1/s3 + 6s2 + 11s + 6
. . .
x3 x2 = x3 x1 = x2
+ + + x1
u y
+ + + 1
–6
– 11
–6
Fig. 16.
3. Choosing state variable : x1 = x2, x 2 = x3
x 3 = – 6x1 – 11x2 – 6x3 + u and y = x1
x1 0 1 0 x1 0
4. State model, x 2 = 0 0 1 x2 0 u
x 3 6 11 6 x 1
3
x1
y = [1 0 0] x2
x3
0 1 0 0 1
Given : A = 0 0 1 , B = 0 0 and C = [1 0 0]
0 2 3 1 1
To Test : Controllability and Observability.
A. Controllability test :
0 1 0 0 1 0 0
1. AB = 0 0 1 0 0 1 1
0 2 3 1 1 3 3
0 1 0 0 0 1 1
A2B = 0 0 1 1 1 3 3
0 2 3 3 3 11 11
Control System SP–21 C (EN-Sem-5)
www.aktutor.in
2. Controllability test matrix is given by
Qc = [B : AB : A2B]
0 1 0 0 1 1
Qc = 0 0 1 1 3 3
1 1 3 3 11 11
0 1 0
3. Consider 0 0 1 = 1 0 ; |Qc| 0
1 1 3
Hence rank of Qc is equal to its order i.e., 3. Therefore the system
is controllable.
B. Observability test :
1
1. C = [1 0 0] ; C = 0
T
0
0 1 0 0 0 0
A = 0 0 1 ; AT = 1 0 2
0 2 3 0 1 3
0 0 0 1 0
1 0 2 0 1
AT CT =
0 1 3 0 0
0 0 0 0 0
T 2 T T T T
(A ) C = (A ) (A C ) = 1 0 2 1 0
0 1 3 0 1
2. Observability test matrix is given by
1 0 0
Qo = [CT : AT CT : (AT)2 CT] = 0 1 0
0 0 1
|Qo| = 1 0
Hence its rank is 3 equal to its order i.e., 3, therefore system is
completely observable.
6. Attempt any one part of the following : (7 × 1 = 7)
a. For the given transfer function G(s) H(s) = 2/s(1 + 0.5s) (1 + 0.05s)
Determine phase crossover frequency and gain margin.
Ans.
2
Given : G(s) =
s(1 0.5s) (1 0.05 s)
To Draw : Bode plot.
Solved Paper (2018-19) SP–22 C (EN-Sem-5)
www.aktutor.in
1. Magnitude plot :
2 2
G(s) = =
s(1 0.5s) (1 0.05 s) s s
s 1 1
2 20
First corner frequency (pole) (1) = 2 rad/s
Second corner frequency (pole) (2) = 20 rad/s
2. Phase plot :
2
G(s) =
s s
s 1 1
2 20
Put s = j
2
G(j) =
j(1 0.5 j) (1 0.05 j)
= – 90° – tan – 1 0.5 – tan – 1 0.05
Magnitude plot :
Table 1.
S. No. Factor Corner Asymptotic log
frequency magnitude
characteristics
Straight line of constant
1. 1/s None slope (– 20 dB/dec) passing
through at = 1
Straight line of constant
2. 1/(1 + 0.5s) 1 = 2 slope (– 20 dB/dec)
originating from 1 = 2
Straight line of constant
3. 1/(1 + 0.05s) 2 = 20 slope (– 20 dB/dec)
originating from 2 = 20
Straight line of constant
4. 2 None slope of 0 dB/dec starting
from 20 log 2 = 16 dB point
Phase plot :
Table 2.
S. No. (rad/sec) (degrees)
1. 0 – 90°
2. 0.2 – 96.28°
3. 1 – 119.42°
4. 2 – 140.71°
5. 5 – 172.23°
6. 10 – 195.25°
7. 15 – 209.27°
Control System SP–23 C (EN-Sem-5)
www.aktutor.in
Bode plot :
1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
100
80
– 60 dB/dec
60
– 40 dB/dec
40
20
GM =22 dB
– 20 dB/dec
in deg
–80°
–100°
–120°
–140°
–160°
–180°
pc = 6.35 rad/sec
–200°
–220°
–240°
0.1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
Fig. 17.
Results from bode plot :
Phase crossover frequency (p) = 6.35 rad/sec
Gain margin = 22 dB
Solved Paper (2018-19) SP–24 C (EN-Sem-5)
www.aktutor.in
b. The forward path transfer function of unity feedback
control system is G(s) = 100/s(s + 6.45). Find the resonance
peak Mr, resonant frequency r and bandwidth of the closed
loop system.
Ans.
100
1. G(s) = and H(s) = 1
s(s 6.45)
100
= ...(1)
s2 6.45s 100
n2
3. Compare eq. (1) with
s 2 n s n 2
2
1 1
5. Resonant peak, Mr = = 1.65
2
2 1 2 0.32 1 0.322
2n
G(s) =
s(s 2n ) (1 Ts)
2. We can obtain a qualitative indication on the bandwidth properties
by referring to Fig. 19, which shows the plots of |M(j)| versus
for n = 1, = 0.707, and various values of T. Since the system is
now of the third order, it can be unstable for a certain set of system
parameters.
Control System SP–25 C (EN-Sem-5)
www.aktutor.in
|M(j )|
1
T=5 G(s) =
s(s + 1.414)(1 + Ts)
T=1
T = 0.5
1.0
T=0
(rad/sec)
Fig. 18. Magnification curves for a third-order system with a
forward-path transfer function G(s).
3. The effect of adding a pole to the forward-path transfer function is
to make the closed-loop system less stable, while decreasing the
bandwidth.
4. The unit-step responses of Fig. 19 show that for larger values of T,
T = 1 and T = 5, the following relations are observed :
i. The rise time increases with the decrease of the bandwidth,
ii. The larger values of Mr also correspond to a larger maximum
overshoot in the unit-step responses.
y(t)
T=5
T=1
T = 0.5
1
G(s) =
s(s + 1.414)(1 + Ts)
Time (sec)
Fig. 19. Unit-step responses of a third-order system with a
forward-path transfer function G(s).
R2
ei eo
C
0º
– 45º m
– 90º
log
1 m 1
T T
Fig. 21. Bode plot of lag compensator.
1
4. The two corner frequencies are = , upper corner frequency
T
1 1
for zero at s = – ,= , lower corner frequency for a pole at
T T
1
s=–
T
5. The maximum phase-lag, m occurs at mid frequency m between
upper and lower corner frequencies.
1 1 1
log10 m = log log 10
2 T T
1
m =
T
6. The phase angle Eo(j)/Ei(j) calculated as
Eo ( j)
= tan–1 (T) – tan–1 (T)
Ei ( j)
Control System SP–27 C (EN-Sem-5)
www.aktutor.in
1
At = m = , the phase angle is m :
T
1
tan m =
2
1
sin m =
1
Pole-zero configuration is shown in Fig. 22.
Im
Re
1 1
– –
T T
s-plane
Fig. 22 Pole zero configuration.
Control System SP–1 C (EN-Sem-5)
www.aktutor.in
B.Tech.
(SEM. V) ODD SEMESTER THEORY
EXAMINATION, 2019-20
CONTROL SYSTEM
Time : 3 Hours Max Marks : 70
Note : Attempt all sections. If any missing data is required, then choose
suitably.
SECTION-A
1. Attempt all parts of the following : (2 × 7 = 14)
a. The impulse response of unity feedback control system is
C(t) = te – t + 2be – t, find transfer function.
b. Give example of type zero, type one and type two systems.
H1
H2
Fig. 1.
Solved Paper (2019-20) SP–2 C (EN-Sem-5)
www.aktutor.in
b. The close loop transfer function is given by
T(s) = k (s + z)/ s2 + 4s + 8, where k, z is adjustable.
i. If r(t) = t finds k and z so that steady error is zero.
ii. for the value of k, z obtain in part (i) find e() for input
r(t) = t2/2
SECTION-C
3. Attempt any one part of the following : (7 × 1 = 7)
a. Explain the effect of sensitivity on feedback system,
Determine sensitivity of negative close loop system where
G(s) = 20/ s(s + 4), H(s) = 0.5 wrt forward path transfer
function, feedback path transfer function.
Solved Paper (2019-20) SP–4 C (EN-Sem-5)
www.aktutor.in
SOLUTION OF PAPER (2019-20)
Note : Attempt all sections. If any missing data is required, then choose
suitably.
SECTION-A
1. Attempt all parts of the following : (2 × 7 = 14)
a. The impulse response of unity feedback control system is
C(t) = te – t + 2be – t, find transfer function.
Ans. Numerical :
Given : r(t) = (t) , c(t) = te– t + 2be– t
To Find : Transfer function.
1. Laplace transform gives
R(s) = 1
c(t) = t e– t + 2b e– t
2. Laplace transform of c(t) gives
1 2b
C(s) =
(s 1)2 (s 1)
1 2b (s 1) (2bs 2b 1)
C(s) = ...(1)
(s 1)2 (s 1)2
3. For unity feedback control system
C ( s) G (s )
=
R( s) 1 G (s )
G (s )
C(s) = R(s) ...(2)
1 G (s )
4. Putting value of C(s) and R(s) in eq. (2)
(2 bs 2b 1)
(1 + G(s)) = G(s)
( s 1)2
(2 bs 2b 1) (2 bs 2b + 1)
= G( s) 1
( s 1)2 (s 1) 2
s2 1 2 s 2bs 2b 1
= G( s)
( s 1)2
2bs 2b 1
5. Transfer function =
s2 2s 2bs 2b
Control System SP–5 C (EN-Sem-5)
www.aktutor.in
b. Give example of type zero, type one and type two systems.
Ans. Type zero system :
k(s z1 ) (s z2 )
G(s) H(s) =
s0 (s p1 ) (s p2 )
Type one system :
k(s z1 ) (s z2 )
G(s) H(s) =
s1 (s p1 ) (s p2 )
Type two system :
k(s z1 ) (s z2 )
G(s) H(s) =
s2 (s p1 ) (s p2 )
0
1
Real 0°
– 90°
Fig. 1.
H1
H2
Fig. 2.
Ans.
i. Masson’s gain formula :
1. The overall gain can be determined by Mason’s gain formula given
below
k
Pk k
T=
k1
H2
Fig. 3.
G3
G2
R(s) + G4 C(s)
G1 G2
1 + G4H 1
– –
H2
Fig. 4.
G3
G2
R(s) + G4 C(s)
G1G2 1
1 + G4H 1
– –
H2
Fig. 5.
After simplication
G1 + G 2
R(s) G1G2 G G4 C(s)
1+ 3 G2
1 + G1G2H2 G2 1 + G4H 1
–
Fig. 6.
Fig. 7.
R(s) G1G4(G2 + G 3) C(s)
(1 + G1G2H2) (1 + G 4 H 1)
–
Fig. 8.
Control System SP–9 C (EN-Sem-5)
www.aktutor.in
G1G4(G2 + G3)
R(s) (1 + G1G2H2)(1 + G 4H1) C(s)
G1G4(G2 + G3)
1+
(1 + G1G2H2)(1 + G4H1)
Fig. 9.
C (s ) G1G4 (G2 G3 )
=
R( s) 1 G1G2 H2 G4 H1 G1G2G4 H1 H2 G1 G4 (G2 G3 )
0 = lim sE(s)
s0
or R(s) = C(s)
k( s z)
R(s) =
s2 4s 8
s2 + 4s + 8 = k(s + z)
3. Compare the coefficients of s in the both sides.
s1 coefficient 4 = k
s0 coefficient 8 = kz
or, z=2
ii.
1. Given, r(t) = t2/2
1
2. Taking Laplace transform, R(s) = –
s3
3. We know ess = lim sE(s)
s0
k( s z)
E(s) = R(s) – C(s) = R(s) – R(s)
(s 2 4 s 8)
4. Putting the value of k and z then we get,
4( s 2)
E(s) = R(s) – R( s)
(s 2 4 s 8)
Solved Paper (2019-20) SP–10 C (EN-Sem-5)
www.aktutor.in
(s 2 4 s 8 4 s 8)
E(s) = R(s)
s2 4 s 8
s2 1 1
ess = lim s 2
3 = = 0.125
s0 s 4s 8 s 8
s
4 6 6 k
3 11 1 0
s
2 60 k 0
s
11
60
–11k 0
1 11
s 60
11
0
s k
2. For system to be stable there should not be sign change in the first
column.
k > 0 for s0
Control System SP–11 C (EN-Sem-5)
www.aktutor.in
and 60 – 121 k > 0
60 > 121 k
0.49 > k, k < 0.49
3. So the range of k is 0 < k < 0.49
kmar = 0.49
4. To find frequency, find out roots of auxiliary equation at marginal
value of k
60 2
A(s) = s k =0
11
60 s2 + 11 k = 0
60 s2 + 11 × 0.49 = 0
60 s2 + 5.39 = 0
60 s2 = – 5.39
5.39
s2 =
60
s2 = 0.089
s = 0.089 = ± j0.298
5. Comparing with s = ± j
so = Frequency of oscillations = 0.298 rad/sec.
A =
= 0 –2 –1, 0 0 RE
C
E
Fig. 10.
6. For stability
N = P–Z
where N = Number of encirclement
P = Number of poles enclosed by the contour
Z = Number of zeros enclosed by the contour
N=2
P= 0
2= 0–Z
Z= 2
So the system is unstable because its value lies in the right half
s-plane.
–7
–14
–8
Fig. 11.
SECTION-C
– 2.25 j6 – 2.25 j 6
= =
– 2.25 j6 10 7.75 j 6
(– 2.25 j6) (– 7.75 j 6)
=
( j6 7.75)( j6 – 7.75)
17.43 – j 46.5 – j13.5 – 36
=
– 36 – 60.06
– 18.57 – j60
=
– 96.06
V ia E T J
if
Lf
Constant
Fig. 12.
2. Let Ra = Armature resistance
La = Armature inductance
ia = Armature current
E = Induced emf in armature
V = Applied armature voltage
= Angular displacement of the motor shaft
T = Torque developed by motor
J = Equivalent moment of inertia of motor shaft
and load referred to the motor
B = Equivalent viscous friction coefficient.
3. Apply KVL in armature circuit
Solved Paper (2019-20) SP–16 C (EN-Sem-5)
www.aktutor.in
La dia
V= Ra ia E ...(1)
dt
4. In the armature-controlled DC motor, the field current is held
constant. For a constant field current, the flux becomes constant,
and the torque becomes directly proportional to the armature
current so that
T ia
T = Ki a ...(2)
5. When armature is rotating, an emf is induced
K b d
E= ...(3)
dt
6. The armature current produces the torque which is applied to the
inertial mass and friction hence the force balance equation is
Jd 2 Bd
= T = KIa ...(4)
dt 2 dt
7. Taking the Laplace transform on both sides of eq. (4) and (1)
(sLa + Ra) Ia(s) + E(s) = V(s) [Initial condition is zero]
(Js2 + Bs) (s) = T(s) = KIa(s) ...(5)
8. The transfer function of this system is obtained as
( s) K
= ...(6)
V (s ) s[ La Js 2 ( BLa JRa ) s Ra B KK b ]
9. Block diagram is shown in Fig. 13.
+ 1 Ia(s) K (s)
V(s)
– sLa + Ra s(Js + B)
E(s)
sK b
Fig. 13.
ii. Field controlled DC motor :
1. A schematic diagram of a field controlled DC motor (DC servo
motor) shown in Fig. 14.
Rf Ra
ia = Constant
ef if Lf J ea
T
B
Fig. 14.
Control System SP–17 C (EN-Sem-5)
www.aktutor.in
2. Here, Rf = Field winding resistance
Lf = Field winding inductance
If = Field winding current
Ra = Armature resistance
ia = Armature current
= Angular displacement.
3. The torque T developed by the motor is proportional to product of
the air-gap flux and armature current ia so we get
T = K1 i a ...(1)
where K1 is constant.
4. But the air gap flux and the field current if are proportional for
the usual operating range of the motor and ia is assumed to be
constant, we can rewrite the above equation as
T = K2 i f ...(2)
where K2 is a constant.
5. The equations for this system are
dif
Lf R f if = e f ...(3)
dt
Jd 2 Bd
and = T = K2i f ...(4)
dt 2 dt
6. By taking the Laplace transform on both sides of eq. (3) and (4)
where all initial conditions are zero, we get
(Lfs + Rf) if(s) = Ef(s) ...(5)
(Js2 + Bs) (s) = K2if(s) ...(6)
7. From the above equations, the transfer function of this system is
obtained as
(s) K2
= ...(7)
E f (s ) s( Lf s Rf )( Js B)
8. Block diagram is shown in Fig. 15.
1 If(s) K2
Ef(s) (s)
L fs + R f s(Js + B)
Fig. 15.
where, Kp = lim
s0
G( s) H ( s)
Kp = Position error constant.
1
For step input, ess =
1 Kp
Case 1 : For type ‘0’
Kp = Constant
ess = Constant
Case 2 : For type ‘1’
Kp =
1
ess = = 0
1
Case 3 : For type ‘2’
Kp =
1
ess = =0
1
10 (1 4 s)
B. Numerical : G(s) =
s2 (1 s)
H(s) = 1
10 (1 4 s)
1. kp = lim G( s) H (s) lim =
s0 s2 (1 s)
s0
10 (1 4 s)
2. kv = lim G (s)· H ( s) lim s 2 =
s0 s0 s (1 s)
10 (1 4 s)
3. ka = lim s 2 G (s)· H ( s) lim = 10
s0 s0 (1 s)
Control System SP–19 C (EN-Sem-5)
www.aktutor.in
t2
r(t) = 1 t
2
4. Steady state error is given by,
R(s)
ess = lim s
s0 1 G(s) H (s)
R1 R R
ess = 2 3
1 k p kv ka
1 1 1/2
= = 0.05
1 10
0 s-plane
>1 – 1
>1 =–
= 0 – 1 –
=1
0 > > – 1
=0
Fig. 16.
ii. A. Effect of adding poles to the forward path transfer
function :
1. Adding a pole at s = – 1/T to the forward-path transfer function of
lead to
2n
G(s) =
s(s 2n ) (1 Ts)
2. We can obtain a qualitative indication on the bandwidth properties
by referring to Fig. 17, which shows the plots of |M(j)| versus
for n = 1, = 0.707, and various values of T. Since the system is
now of the third order, it can be unstable for a certain set of system
parameters.
Solved Paper (2019-20) SP–20 C (EN-Sem-5)
www.aktutor.in
|M(j )|
1
T=5 G(s) =
s(s + 1.414)(1 + Ts)
T=1
T = 0.5
1.0
T=0
(rad/sec)
Fig. 17. Magnification curves for a third-order system with a
forward-path transfer function G(s).
3. The effect of adding a pole to the forward-path transfer function is
to make the closed-loop system less stable, while decreasing the
bandwidth.
4. The unit-step responses of Fig. 18 show that for larger values of T,
T = 1 and T = 5, the following relations are observed :
i. The rise time increases with the decrease of the bandwidth,
ii. The larger values of Mr also correspond to a larger maximum
overshoot in the unit-step responses.
y(t)
T=5
T=1
T = 0.5
1
G(s) =
s(s + 1.414)(1 + Ts)
Time (sec)
Fig. 18. Unit-step responses of a third-order system with a
forward-path transfer function G(s).
B. Effect of adding a zero to the forward path transfer
function :
1. Consider unity feedback system having forward path transfer
function as,
2n
G (s) = ...(1)
s(s 2 n )
2. Let us add a zero at s = – 1/T to the transfer function so that
eq. (1) becomes
2n (1 Ts)
G(s) =
s(s 2 n )
Control System SP–21 C (EN-Sem-5)
www.aktutor.in
3. The closed-loop transfer function for unity feedback system is
C ( s) 2n (1 Ts)
= 2
R( s) s (2 n T 2n ) s n2
1 2
B.W. = b b 4 4n
2
where b = 42 2n + 4 3n + T – 22n – 4n
for = 0.707 and n = 1
B.W
1
T
0
Fig. 19.
4. The general effect of adding a zero to the forward path transfer
function is to increase the bandwidth of the closed-loop system.
|M(j )|
T=0 For n = 1
T = 0.2
= 0.2
T = 0.5
T=1
1 T=2 T=5
0.707
T
b
1
b
2
Fig. 20.
5. So bandwidth increases as T-increases
c(t)
T=5 T = 0.2
1
T=1
t
tr tr tr
3 12
Fig. 21.
tr3 < tr2 < tr1
So as T increases rise time decreases.
Solved Paper (2019-20) SP–22 C (EN-Sem-5)
www.aktutor.in
5. Attempt any one part of the following : (7 × 1 = 7)
a. A unity feedback system has an open loop transfer
function.
Draw the root locus for the system.
G(s) H(s) = k (s + 2)/(s + 3) (s2 + 2s + 2)
k ( s 2)
Ans. G(s) H(s) =
(s 3) ( s 2 2s 2)
Step-1 :
Poles, P1 = – 3, P2 = – 1 + j, P3 = – 1 – j
Zeros, Z = – 2
No of poles = 3, No of zeros = 1
Step-2 : No of asymptotes = P – Z = 3 – 1 = 2
Angle of asymptotes
(2 k 1) 180 (2 0 1) 180
1 = = 90°
PZ 2
(2 1 1) 180
2 = = 270°
2
Step-3 : Centroid
R.P of poles R.P of zeros
=
PZ
( 3 1 1) ( 2) 5 2
= = – 1.5
2 2
Step-4 : Breakaway point :
Characteristics equation,
1 + G(s) H (s) = 0
k ( s 2)
1
(s 3) (s 2 2s 2)
(s + 3) (s2 + 2s + 2) + k ( s + 2) = 0
s3 + 5s2 + 8s + 6 + k (s + 2) = 0
s3 5s 2 8 s 6
k=
s2
dk
=0
ds
Control System SP–23 C (EN-Sem-5)
www.aktutor.in
d s3 5s 2 8s 6
=0
ds s2
(s + 2) [–3 s2 – 10s – 8] – [– s3 – 5s2 – 8s – 6] × 1 = 0
– 2s3 – 11s2 – 20s – 10 = 0
2s3 + 11s2 + 20s + 10 = 0
s = – 0.80, – 2.34 + 0.84j – 2.34 – 0.84]
That breakaway point is not valid
Step-5 : Intersection point : Characteristics equation,
1 + G(s) H (s) = 0
k (s 2)
1 =0
(s 3) (s 2 2s 2)
s3 + 5s2 + 8s + 6 + ks + 2k = 0
Apply Routh’s criteria,
s
3 1 k+8
s
2 5 6+2k
s
1 34 + 3k 0
5
s
0 6 + 2k 0
34 3k
=0
5
34 + 3k = 0
3k = – 34
k = – 34/3 = – 11.33
Auxiliary equation
5s2 + 6 + 2k = 0
5s2 + 6 + 2 × – 11.33 = 0
s2 = 3.33 = s = ± 1.82
That intersection point is not valid.
Step-6 : Angle of departure :
d = 180° – [P – Z]
P1 = 90°, P2 = 26.56°
Z = 45°
d = 180° – [90° + 26.56° – 45°] = 108.44°
Solved Paper (2019-20) SP–24 C (EN-Sem-5)
www.aktutor.in
Imj
k=
1
–1 + j
j
P2 – 108.44
Z 90°
–2 –1.5 –1 Re
–3 k=0
RL 270°
P1
–j
–1 – j
108.44
k=
2
Fig. 22.
b. What is the necessary condition for stability ? Explain
limitation of Routh’s stability method, construct Routh
array and determine the stability of the system whose
characteristics equation, s5 + 2s4 + s3 + 2s2 + s + 4 = 0
Ans.
i. Necessary condition :
i. All the coefficients of characteristics equation must have same
sign.
ii. There should be no missing term.
iii. All poles of transfer function should be in left half of s-plane.
iv. The degree of denominator polynomial of transfer function is greater
or equal to that of numerator polynomial.
ii. Limitation :
1. It is valid only if the characteristics equation is algebraic.
2. If any coefficient of the characteristics equation is complex or contain
power of e then this criterion cannot be applied.
3. It gives information about how many roots are lying in the RHS of
the s-plane but values of the roots are not available. Also it cannot
distinguish between real and complex roots.
Control System SP–25 C (EN-Sem-5)
www.aktutor.in
iii. Numerical :
1. Routh array is :
s5 1 1 1
s4 2 2 4
s3 0 1
s2
s1
s0
2. Putting ( very small positive value) in place of 0 in the first column
Routh array
s5 1 1 1
s4 2 2 4
s3 1 0
2 2
s2 ( ve) 4
2 2
4
s1 ( ve) 0
2 1
s0 4
3. Since there is sign change in the first column of Routh array,
therefore system is unstable.
6. Attempt any one part of the following : (7 × 1 = 7)
a. Draw Bode Plot of unity feedback control system having
10
OLTF G(s) = .
s(1 0.2s) (1 0.2s)
Determine GM, PM, gain cross over frequency, phase cross
over frequency and discuss stability of closed loop system.
Ans.
10
Given : G(s) =
s(1 0.2 s) (1 0.2s)
To Draw : Bode plot.
1. Magnitude plot :
10 10
G(s) =
s(1 0.2s) (1 0.2 s) s s
51 1
5 5
Corner frequency (1or 2) = 5 rad/s
Solved Paper (2019-20) SP–26 C (EN-Sem-5)
www.aktutor.in
2. Phase plot :
10
G(s) =
s s
s 1 1
5 5
Put, s = j
10
G(j) =
( j) (1 0.2 j) (1 0.2 j)
Phase plot :
Table 2.
80
– 60 dB/dec
60
47.95 – 40 dB/dec
40
20
– 20 dB/dec
–20
GM =46 dB
–40
–60
–80°
–100°
in deg
–180°
–225°
=5 rad/sec
pc
–270°
–360°
0.1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 102
Fig. 23.
Results from bode plot :
A. Phase crossover frequency (pc) = 5 rad/sec
B. Gain margin = 46 dB
C. Phase margin = 47°
Solved Paper (2019-20) SP–28 C (EN-Sem-5)
www.aktutor.in
D. Phase crossover frequency = 0.2 rad/sec
E. Stability : The gain margin as well as phase margin are both
positive therefore, the system is stable.
– 1 + j0 0
Re
A r
(r – )
Fig. 25.
4. From Fig. 25, OB = G(j)
OA = – 1
AB = OB – OA = G(j) – (– 1)
AB = 1 + G(j)
5. From eq. (2)
C ( j) OB
= M() =
R( j)
AB
Control System SP–29 C (EN-Sem-5)
www.aktutor.in
C( j) OB r
=
= (r – )
R( j) AB
C( j)
= M() e j ()
R( j)
where M(j) is the magnitude and () = r – .
6. Frequency response consists of two parts :
i. Magnitude,
ii. Phase angle.
iii. Both can be plotted against different values of .
7. Hence frequency response of closed loop system is plot of
magnitude and phase angle.
Correlation :
1. For 2nd order system , the transfer function is,
2
C ( s) n
= 2
R( s) s 2 n s 2n
where, = Damping factor
n = Natural frequency of oscillations
+ n2
R(s) C(s)
2
– s + 2 ns
Fig. 26.
2. Closed loop frequency response is,
C( j) 2n
= T(j) =
R( j) ( j)2 2 n ( j) n2
2n 1
= = 2
2 n ( j ) 2n
2
1 2 j
n n
1
= ...(3)
(1 u2 ) j 2u
where u = /n, normalized driving frequency.
1
|T(j)| = M = ...(4)
(1 u2 )2 (2u)2
2u
and T(j) = = tan 1 ...(5)
1 u2
3. The steady state output is
1 2u
c(t) = sin t tan 1
(1 u2 )2 (2u)2 1 u2
Solved Paper (2019-20) SP–30 C (EN-Sem-5)
www.aktutor.in
From eq. (4) and (5) when
u = 0, M = 1 and = 0
1
u = 1, M = and = –
2 2
u = , M = 0 and =–
4. The frequency where M has a peak value is called the resonant
frequency. At this frequency the slope of magnitude curve is zero.
If r = Resonant frequency.
ur = r/n is normalized resonant frequency.
2 2
dM 1 4(1 ur )ur 8 ur
= 0
du u ur 2 (1 u2 ) 2 (2u )2 3/ 2
r r
– 4(1 – ur2) ur + 82ur = 0
– 4ur(1 – ur2 – 22) = 0
ur = 1 2 2
r = n 1 2 2 ...(6)
1
Mr = ...(7)
2 1 2
5. The phase angle of T(j) at resonant frequency is
1 2 2
r = tan 1
1.0
Magnitude
M
Mr
M( = 0.707)
0
ur = r / n
Normalized frequency, u
Fig. 27. Frequency response magnitude characteristics.
Control System SP–31 C (EN-Sem-5)
www.aktutor.in
0°
( = 1)
– 90°
Phase
angle
– 180°
0 1.0
Normalized frequency u
Fig. 28. Frequency response phase characteristic.
10
0
|D(j )| dB
– 10 20 log
– 20
– 20 dB/decade
0°
– 30°
D(j )
– 60°
– 90°
0.01 0.1
1 1
(rad/sec)
Fig. 29.
0 1 1 3 3
AB = 1 1 1 2 2
1 1 0 1 1
0 1 1 0 1 1 0 0 1
A2 = 1 1 1 1 1 1 2 3 0
1 1 0 1 1 0 1 2 0
0 0 1 3 1
A 2B = 2 3 0 2 0
1 2 0 1 1
QC = [B : AB : A2B]
3 3 1
= 2 2 0
1 1 1
QC = 0
Hence rank of QC is equal to its order i.e., 3 Therefore the system
is controllable.
ii. observability :
0 1 1
A T = 1 1 1
1 1 0
1
CT = 0
0
Solved Paper (2019-20) SP–34 C (EN-Sem-5)
www.aktutor.in
0 1 1 1 0
ATCT = 1 1 1 0 1
1 1 1 0 1
0 1 1 0 1 1 0 2 1
(AT)2 = 1 1 1 1 1 1 0 3 2
1 1 1 1 1 1 1 0 0
0 2 1 1 0
(AT)2CT = 0 3 2 0 0
1 0 0 0 1
Q0 = [CT : ATCT : (AT)2 CT]
1 0 0
= 0 1 0
0 1 1
=1
Hence rank of Q0 is equal to its order i.e., 3, so the system is
observable.