Cse UNIT 4 Theory
Cse UNIT 4 Theory
UNIT-4
Poles and Zeros – Stability of Systems – Conditions for Stability – Routh –
Hurwitz Stability Criterion – Root Locus Analysis – Properties – Root Locus Plots
Part –A (2 Marks)
break away or break in points, form a equation for K from the characteristic equation and
ot
differentiate the equation of K with respect to s. Then find the roots of the equation dK/dS = 0.
hN
The roots of dK/dS = 0 are break away or break in points provided for this value of root the gain
ec
In equation (i) all the coefficient are positive and so the necessary condition for stability
is satisfied. In equation (ii), (iii) and (iv) some of the coefficient are negative and some of the
es
coefficient are missing. Hence the necessary condition for stability is not satisfied
ot
hN
ec
Part –B
M
CONCEPTS OF STABILITY:
Stability in a system implies that small changes in the system input, in initial conditions
(or) in system parameters do not result in large changes in system output. The linear time
invariant system is stable, if the following two notions of system stability are satisfied.
▷ When the system is excited by the bounded input, the output is bounded.
▷ In the absence of the input, the output tends towards zero irrespective of initial
conditions, [This stability concept is known as asymptotic stability]
▷ If the system is subjected to an unbounded input and produces an unbounded response,
the system is unstable.
But if it is subjected to a bounded input and produce an unbounded response the
system is unstable. The output of the unstable system may increase to certain extent and then
the system may breakdown (or) become non-linear, after the output exceeds certain magnitude
so that a linear Mathematical model no longer applies.
Absolute stability:
The linear system is characterized as absolutely stable system with respect to a
parameter of the system if it is stable for all value of this parameter.
Conditionally stable system: (qualitative Stability)
A system is said to be conditionally stable with respect to a parameter, if a
system is stable for only certain bounded ranges of values of this parameter
Relative stability: (Quantitative stability)
Relative stability is the Quantitative measure of how fast the transients die out in the
system. It may be measured by Relative settling times of each root or pair of roots.
As a root [or a pair of roots] moves far from the imaginary axis as shown in figure, the
relative stability of the system improves
equation ( if ao ≥ 0) indicates that the system is either unstable (or) atmost limitedly stable.
ot
hN
▷ The positiveness of the coefficients of the characteristic equation ensures the negativeness
of real roots but does not ensure the negativeness of the real parts of complex roots for III &
higher order system .Therefore it cannot be a sufficient condition for stability of III & higher
order systems.
▷ If the characteristic equation of the system is of degree higher than second, the possibility
of its instability cannot be excluded even when all the coefficients of characteristic equation
are positive. The First step in analysing the stability of the system is to examine its
characteristic equation. If some of the coefficient are zero (or) negative it can be concluded
that the system is not stable.
▷ On the other hand if all the coefficients of the characteristic equations are positive (or
negative), the possibility of stability of system exists and one should proceed further to
examine the sufficient condition of the stability.
2. Explain about the Routh stability criterion in details
Routh Array:
S n a0 a2 a4 a6
S n-1 a1 a3 a5
S n-2 b1 b2 b3
S n-3 c1 c2
S n-4 d1 d2
. . .
. . .
. . .
s 2 e1 an
s 1 f1
S o an
b1 = (a1 a2 – a0 a3) \a1;
b2 =(a1 a4-a0a5)\a1;
c1= (b1 a3 – a1 b2)\b1;
c2=(c1 b5 – a1 b3)\b1;
d1= (c1 b2 – b1 c2)\c1;
d2= (c1 b3 -b1 c3)\c1
This process is continued till we get zero as the last coefficients, in a particular row. In
the process of generating the routh Array, the missing terms are regarded as zero, also all the
elements of any row can be divided by a positive constants during the process to simplify the
es
work.
ot
hN
For a system to be stable, it is necessary and sufficient that all elements of I Colum of
st
M
Routh array of its characteristic equation be positive if ao>o. If this condition is not met, the
system is unstable and number of sign changes of the terms of the I column of routh array
corresponds to the number of roots of the characteristic equation in the right half of the s-plane.
characteristic equation
S +8s +18s +16s+5 =0.
4 3 2
Necessary condition:
All the coefficients given in the polynomial are positive (or of same sign) and none of the
terms are missing. Hence the system may be stable
Sufficient condition:
For the system to be stable it has to satisfy the necessary condition as well as the
sufficient condition. The sufficient condition can be checked using Routh’s stability condition.
Routh array:
S4 1 18 5
S3 8 16 0
S2 16 5 0
S1 1.7 0
S0 5
On examining the elements of I column of Routh array is observed that all the elements
have positive sign and there is no sign change. Hence all the roots are lying on left half of its
plane.
▷ Stable System
▷ All four roots lie on left half of s-plane.
Necessary condition:
All the coefficient given in the polynomial are positive and none of the term are missing.
Hence the system may be stable.
Sufficient condition:
For the system to be stable it has to satisfy the necessary condition as well as the
sufficient condition. The sufficient condition can be checked using Routh stability condition.
Routh array:
S6 1 8 20 16
S5 1 6 8 0 (Dividing throught by 2)
S4 1 6 8
S3 0 0 0
S2 3 8
S1 0.333 0
es
ot
S0 8
hN
A (s) = 1s +6s +8
4 2
M
dA /dS = 4s +12s]
3
After forming a Routh array if we examine I column all elements seems to be positive & it
satisfies the Routh stability criterion. Hence the system is stable
All the roots of the given system are located in the left half of the s-plane. Find out the root
that passes entire row all zero
The roots of the Auxillary polynomial are
s =x
2
X +6x+8=A(s)
2
as conjugate pairs).
ot
Number of Branches There are n branches of the locus, one for each pole of the loop
hN
gain.
ec
Starting and Ending Points The locus starts (when K=0) at poles of the loop gain, and ends
M
(when K→∞ ) at the zeros. Note: there are q zeros of the loop gain
as |s|→∞ .
Locus on Real Axis The locus exists on real axis to the left of an odd number of poles
and zeros.
Asymptotes as |s|→∞ If q>0 there are asymptotes of the root locus that intersect the real
axis at
sum of poles-sum of zeros
σ =
A n-m
0
180 (2q+1)
, and radiate out with angles φ = , where q=0,1,2...
A n-m
Break-Away and -In Points on There are break-away or -in points of the locus on the axis
Real Axis wherever dK/ds =0,if r numbers of branches of root locus meet at a
point,then they break away at an angle of ±1800/r.
Angle of Departure from Complex Angle of departure from pole pj is
Pole 0
φp =±180 2q+1 +φ; q=0,1,2,....
Where Φ is the net angle contribution at the pole by all other open
loop poles and zeros
Angle of Arrival at Complex Zero Angle of arrival at zero zj is
0
φz =±180 2q+1 +φ; q=0,1,2,....
Where Φ is the net angle contribution at the zero by all other open
loop poles and zeros
Locus Crosses Imaginary Axis Use Routh-Horwitz to determine where the locus crosses the
imaginary axis.Alternatively they can be evaluated by letting s=jω in
the characteristic equation and equating the real and imaginary part
to zero,to solve for ω and K. The values of ω are the intersection
points on imaginary axis and K is the value of gain at the
intersection points.
Determine Location of Poles, Rewrite characteristic equation as D(s)+KN(s)=0. Put value of K into
Given Gain "K" equation, and find roots of characteristic equation.
Determine Value of "K", Given n
Pole Locations sa + p
i
i=1
Rewrite characteristic equation as K= m
=
sa + z
i
i=1
, replace "s" by the desired pole location and solve for K. Note: if "s"
is not exactly on locus, K may be complex, but the imaginary part
should be small. Take the real part of K for your answer.
Determination of the open loop gain for a specified damping of the dominant root. A
hN
damping line making an angle θ= cos– 1 Σ with negative real axis is drawn for the specified
ec
damping in the region where the rough root locus edge intersects the delta line, the trial error
M
procedure is adopted along the delta line to determine the point of intersection accurately. This
point is then the desired root at which the open loop gain is computed. Roots along other
branches can then be obtained for this open loop gain using the procedure outline above.
The breakaway point of the root locus are determined from the roots of the equation dk/ds = 0
M
R Branches of root locus which meet at a point, break away at an angle +or - 180 /r
1 + k /[ S(s+1) (s+2)] =0
S(s+1) (s+2)+k=0
- (s +S) (s+2) = K
2
K= -s +s -2s -2s
3 2 2
K = - (s +3s +2s)
3 2
S1 = - 0 . 423
S2 = - 1 . 577
In b\w -1 and -2 there is no root locus exists.So no break away ,breakin poin exists.
Step 5:To find the angle of departure and angle of arrival:
No complex poles in loop gain, so no angles of departure.
No complex zeros in loop gain, so no angles of arrival.
Routh array:
S3 + 3s +
2 2s + k = 0
S3 1 2
S2 3 k
S1 (6-k)/3 0
0<K <6
Sub. K = 6 in Routh array
IN S 1 - row, all element becomes zero
Auxiliary polynomial is,
A = 3s2 + 6 = 0
S = ±j √2
The crossing point on root locus is
S1 = + j 1. 414
S2 = - j 1.414
es
ot
hN
ec
M
1
7. Consider the transfer function 2 . Sketch the Root locus.
s( s +5s+6)
We can rewrite the open loop transfer function as G(s)H(s)=N(s)/D(s) where N(s) is the
numerator polynomial, and D(s) is the denominator polynomial.
N(s)= 1, and D(s)= s3 + 5 s2 + 6 s.
0
180 (2q+1)
Angle of asymptotes =φ = where q=0,1,2...
A n-m
es
ot
hN
ec
M
Here n= 3 and m= 0.
ØA=± 180 (2q+1) / (n-m) q = 0,1,2, (n-m-1)
q= (3-0-1)
q=2
When q=0,
ØA0 = ±180(2*0+1)/ 3-0
ØA0=±60*
When q=1,
ØA1 = ±180 (2*1+1) /3-0
ØA1=±180*
When q=2,
ØA2= ±180(2*2+1)/ 3-0
ØA2= ±30*
sum of poles-sum of zeros
Centroid =σ = =0-3-2-0/3 =-1.87
A n-m
R Branches of root locus which meet at a point, break away at an angle +or - 180 /r
ot
hN
1 + k /[ S(s+1) (s+2)] =0
D(s)+KN(s) = s3 + 5 s2 + 6 s+ K( 1 ) = 0
ec
K= -s3-5s2 -6s
M
K = - (s3+5s2 +6s)
dk/ds = -(3s2 + 10s + 6) => dk/ds =0
From these 2 roots, there exists 2 real roots at s = -2.5, -0.78. These are highlighted on the
Since K is positive and real for s=-0.78,this point is actual breakaway point.
when s=-2.5,the value of K is given by,
K= -[(-2.5) +5(-2.5) +6(-2.5)]=
3 2
S 1
3 6
S2 5 k
S1 (30-k)/5 0
es
K>0
ec
M
0 < K < 30
Sub. K = 30 in Routh array
IN S - row, all element becomes zero
1
No. of branches = 4
There are 4 root locus braches they originate at each open loop pole with k= 0
ec
M
S1 =-2.28
S2 =- 0.7301 + j0 . 348
S3 =-0.731 – j0 . 348
To Check the real part:
The magnitude condition states that s= sa will be a point on the root locus if for that value of
root locus magnitude of |G (s) .H (s) | = 1
The angle criterion states that s = sa will be a point on the root locus for that value of s
k = - (-4.331)
k = 4.331
k is positive and real
Sub s = - 0.731 + 0 . 348i in equ.1
k = - [(0.731+0.348i) + 5 (-0.731+0.348i)
4 3
es
+ 8[ (-0.731+0.348i) + 6 (-0.731+0.348i)]
2
ot
K = 1.793 – j 0.18.
hN
Where Ø is the contribution at pole of all other open loop poles and zeros similarly the angle
of arrival at an open loop zero is given by
Øz = ±180 (2q+1) - Ø ; q=0,1,2…………….
0
Ø=-251 6 0 /
Øp =±180 - 251 6
0 0 /
Øp =±71 6 0 /
8. The intersection of root locus branches with the imaginary axis can be determined by the use
of routh criterion .
S(s+3) (s +2s+2) + k = 0
2
8(s+3) (s +2s+2) + k = 0
2
S +5s +8s +6s
4 3 2 +k = 0
Routh array;
S4 1 8 k
S3 5 6
S2 6.8 k
S1 (40.8-5K)/ 6.8
S0 k
1)(40.8-5k)/6.8 > 0
40.8-5k = 0
5k = 40.8
K = 8.16
2)K > 0
O < k < 8.16
Sub k = 8.16 in s row, all elements becomes zero.
1
S1 = -1.4
ot
S2 = -2.85
hN
S3 = -0.4+j0.7
ec
S4=-0.4-j0.7
M