Tutorial 8 - Stability and The S' Plane: © D.J.Dunn 1
Tutorial 8 - Stability and The S' Plane: © D.J.Dunn 1
Tutorial 8 - Stability and The S' Plane: © D.J.Dunn 1
Dunn 1
TUTORIAL 8 STABILITY AND THE s PLANE
This tutorial is of interest to any student studying control systems and in particular the EC
module D227 Control System Engineering.
On completion of this tutorial, you should be able to do the following.
Define Poles and Zeros
Explain the Characteristic Equation of a Transfer Function.
Explain and interpret Root Locus Diagrams.
Use The Rules for Graphical Construction of Root Loci.
Construct Root Loci.
Explain and apply the Routh-Hurwitz criteria of stability.
If you are not familiar with instrumentation used in control engineering, you should
complete the tutorials on Instrumentation Systems.
In order to complete the theoretical part of this tutorial, you must be familiar with basic
mechanical and electrical science.
You must also be familiar with the use of transfer functions and the Laplace Transform
(see maths tutorials).
1. POLES AND ZEROS
Most higher order transfer functions may be expressed as a polynomial on the top and bottom line such
that
)..... p )(s p (s
).... z )(s z (s
(s)
G(s)
2 1
2 1
i
o
+ +
+ +
= =
The values of s that make the denominator zero are called poles and the values of s that make the
numerator zero are called zeros. The zeros are hence -z
1
, -z
2
and the poles are -p
1
, -p
2
If s =p in any of the factors, then the denominator is zero and G(s) =. If this happens, or if any value of
z makes the numerator infinity, the system is unstable. The analysis however, is conducted not on the
closed loop but on the characteristic equation defined next.
2. CHARACTERISTIC EQUATION
The characteristic equation is the bottom line (denominator) of the closed loop transfer function
when equated to zero.
The diagram shows a system with an open loop transfer function G(s) and feedback transfer function
H(s).
When H(s) = 1 we have UNITY FEED BACK.
Figure 1
1 H G
G
G
cl
+
= or
1/G H
1
+
The characteristic equation is GH +1 =0 or H +1/G =0
If GH has poles and zeros such that
)..... p )(s p (s
).... z )(s z (s
GH
2 1
2 1
+ +
+ +
= then instability occurs when
1
)..... p )(s p (s
).... z )(s z (s
0 1 GH
2 1
2 1
+
+ +
+ +
= = +
It is important not to confuse the poles and zeros of the characteristic equation with those of the
closed loop transfer function.
We learn a lot about the stability and response of a system by examining the poles and the roots of the
characteristic equation.
D.J .Dunn 2
3. THE AFFECT OF POLES AND ZEROS ON DYNAMIC RESPONSES
REAL POLES
Consider the simple open loop transfer function
) (s
1
(s)
G(s)
i
o
+
= =
There is only one pole and this is s =-
It is normal to use an impulse input to study how a system responds, especially regarding stability.
If the input to the system is an impulse, the output in the time domain is simply the reverse Laplace
transform giving
o
(t) =e
-t
Plotting the output for different values of p show that :-
(i) The larger the value of , the quicker it decays.
(ii) The smaller the value of , the longer it takes to decay.
(iii) When the values of become negative the output grows with time.
(iv) The more negative the value of , the quicker it grows with time.
D.J .Dunn 3
This is demonstrated on the following plot
of output against time for various values of
.
- is the pole of the transfer function so we
learn that the position of the pole on the s
plane greatly influences the system.
Negative poles (positive ) are stable and
positive poles are unstable.
The same evidence is produced for a step
input later. Figure 2
REAL ZEROS
Lets examine the closed loop transfer function
( )( ) 5 s 1 s
A s
(s)
G(s)
i
o
+ +
+
= =
There is a zero at s =-A and two poles at s =-1 and s =-5
Again examine the time response for an impulse
input. This is simply the case of inverse Laplace
transform.
4
e
A
4
e
5
4
e
4
e
A (t)
5t 5t t t
o
+ = This is plotted in
the next diagram.
Clearly if A =1 then the zero will cancel one pole. If
A =5 it will cancel the other. If A is between 1 and 5
the result is mid way between the other two. We
conclude that the zero moves the response toward one
pole or another.
Figure 3
The affect of zeros on the root locus plot is more difficult to understand and is covered later.
COMPLEX POLES
Complex poles in the characteristic equation produce oscillatory responses to impulse and step inputs.
Consider the following open loop transfer function.
3) 2)(s 1)(s (s
k
G(s)
+ + +
=
With unity feed back we get the following.
1 G(s)
G(s)
G(s)
cl
+
=
The characteristic equation of the closed loop transfer function is 0 1
3) 2)(s 1)(s (s
k
= +
+ + +
Hence 0 k 3) 2)(s 1)(s (s = + + + +
If any of the brackets are zero then the system is unstable. If k =0, this will clearly happen when s =-1 or
-2 or -3. For all other values of k we must solve the equation
s
3
+6s
2
+11s +6 +k =0
The condition we are examining is s
3
+6s
2
+11s +6 +k =0. If this is plotted for various values of k we
get the result below. It can be seen that when k =0 there are three real roots at -1, -2 and -3. As the value
of k increases the roots at -1 and -2 converge to a single value (around k =0.4) and after that they become
complex roots. The other root is always real.
Figure 4
If we expressed G(s) as
) p - )(s p )(s p - (s
k
3 2 1
p
1
, p
2
and p
3
are the POLES and these are the roots of
the characteristic equation when k =0. It follows that if
3) 2)(s 1)(s (s
k
G(s)
+ + +
= the poles are at -1, -2
and -3.
Consider the case k =1. The characteristic equation has a zero values at one real point (root) around
s =- 3.4. By trial and error with a calculator it is possible to arrive at the solution.
We are finding the roots of s
3
+6s
2
+11s +7 =0 Guessing s =-3.5 we get
s =- 3.5 -42.875 +73.5 38.5 +7 = - 0.875
First correction s =- 3.4 -39.3 +69.36 37.4 +7 = - 0.34
next correction s =- 3.3 - 35.937 +65.34 - 36.3 +7 = 0.103
s =-3.35 - 37.593 +67.335 36.85 +7 = - 0.108
s =-3.33 - 36.926 +66.533 36.63 +7 = - 0.023 0
So the real root is -3.33. Now we need to determine the remaining complex roots.
D.J .Dunn 4
We now know that (s +3.33)(remaining quadratic) =s
3
+6s
2
+11s +7
The quadratic will have the form (as
2
+bs +c) so
(s +3.33)(as
2
+bs +c) =s
3
+6s
2
+11s +7
Multiply out and as
3
+s
2
(b +3.33a)
+s (c +3.33b) +3.33c =s
3
+6s
2
+11s +7
Comparing coefficients we see that a s
3
=s
3
so a =1
(b +3.33a) =6 hence b =2.67
(c +3.33b) =11 hence c =2.1
Checking 3.33c =7 is correct.
The equation is (s +3.33)(s
2
+2.67s +2.1 ) =0
One root is at -3.33 the other roots are found by solving (s
2
+2.67s +2.1) =0
Using the quadratic equation we find
0.564j 1.335
2
1.27 2.67
1 x 2
2.1 x 1 x 4 2.67 2.67
s
2
=
=
=
Note that the two complex roots are conjugate numbers, in other words they are mirror images about the
real axis. This is always true for complex poles; they only exist in conjugate pairs.
To summarise, the characteristic equation is 0 1
3) 2)(s 1)(s (s
k
= +
+ + +
In which case s
3
+6s
2
+11s +6 +k =0
For the case k =1 we have found that the values of s that produces zero values are:
s
1
=- 3.33 s
2
=-1.335 +0.564 j s
3
= -1.335 - 0.564 j and these are the roots.
Note that it is quite possible to have complex poles (when k =0) as well as complex roots.
A root locus plot would require you to evaluate the roots for many values of k and without a suitable
calculation aid this would be very difficult. A graphical method is outlined later. The next section shows
the complete plot for the example just examined.
D.J .Dunn 5
4. ROOT LOCUS DIAGRAMS The s PLANE
Real and complex roots can be plotted on an Argand diagram for different values of k in the previous
section. Note that a locus can be plotted for other parameters such as damping ratio (see tutorial 5). If the
points are joined we obtain a locus known as the ROOT LOCUS DIAGRAM. This is laborious unless a
computer package is used but as this is not available under exam conditions it may be as well to learn
how to sketch them. The result for the case in hand is shown below and because there are three roots, we
have 3 loci. The plane is usually referred to as the s plane.
Figure 5
The three loci are plotted for k values between zero and -. One locus runs from -3 to - along the real
axis. The other two loci approach each other from -2 and -1 and depart from the real axis at a point
between them to form mirror images.
We will look at rules for the construction of root loci later in the tutorial.
D.J .Dunn 6
5 INTERPRETATION OF ROOT LOCUS DIAGRAMS
5.1 INTERCEPT WITH THE IMAGINARY AXIS
Referring to the root locus plot in section 4, the value of k where the locus cuts the imaginary axis is the
critical value where the system changes from stable to unstable. It is the point corresponding to -1 on the
Nyquist diagram. In the example, this occurs at k =60 and the NYQUIST plot of the system confirms it.
The Nyquist plot for
3) 2)(s 1)(s (s
k
G
+ + +
= is shown with k =1 and k =60. When k =60 the plot
passes through the -1 point.
Figure 6
5.2 POSITION OF THE POLES
The reason for sketching these loci is to give an insight into the stability of a system. If a pole appears to
the right of the imaginary axis, the output will grow indefinitely and the system is unstable. If all the poles
are to the left of the imaginary axis, the system is stable. The poles closest to the origin on the left side of
the axis are the MOST SIGNIFICANT as they indicate the output will take a long time to decay (die
away) and dominate the system. The poles furthest from the origin on the positive (right) side indicate
greatest instability. The reasoning behind this was covered earlier in this tutorial but here is a more
detailed look at it.
Remember that the pole p is the value of s that makes the bracket zero. Hence if we had a bracket (s +2)
the pole is at s =-2 so we might write the bracket as (s p).
Consider that a closed loop transfer function is of the form:
( )( )( ) ( )
n 3 2 1 i
o
p s .... p s p s p s
K
G(s)
= =
For a unit step input
i
(s) =1/s the output
o
(s) is
( )( )( )
3 2 1
o
p s p s p s s
K
(s)
=
Change into partial fractions
3
3
2
2
1
1 0
o
p s
A
p s
A
p s
A
s
A
(s)
+ =
Inverse Laplace Transform
t p
3
t p
2
t p
1 0 o
3 2 1
e A e A e A A (s) + + + =
Repeating the process for a unit ramp we get
t p
4
t p
3
t p
2 1 0
o
3 2 1
e B e B e B tB B (t) + + + + =
D.J .Dunn 7
The initial terms A
0
B
0
and tB
1
or others resulting from the solution, are the steady state components
(constant or proportional to time).
With real poles the exponential terms such as Ae
pt
are transient terms. At t =0, Ae
pt
=A. If p is negative,
it will decay with time. If p is positive it will grow with time. The rate at which it decays or grows
depends on the value of p. This is illustrated on the diagram below.
Figure 7
It follows that if p lays on the negative real axis of the s plane, the transient will die away and the further
it is from the origin, the quicker it will take. If p lays on the positive real axis, the transient will grow and
the further it is from the origin, the faster it will grow.
When p is a complex number, the forgoing still applies but to the real part of the number. Without proof,
it can be shown that the smaller the imaginary part of the number, the larger the damping and vice versa.
The affect of the position of the poles on the s plane may be summarised as shown on the diagram.
Figure 8
D.J .Dunn 8
5.3 AFFECT OF THE ZERO ON THE ROOT LOCUS
We have studied the root locus of the closed loop system with an open loop transfer function of
3) 2)(s 1)(s (s
k
G(s)
+ + +
=
Now consider the open loop transfer function
( )
3) 2)(s 1)(s (s
1.5 s k
G(s)
+ + +
+
=
We now have a zero at s =- 1.5. What affect does this have on the root locus?
The closed loop transfer function is
( )
( ) 3) 2)(s 1)(s (s 1.5 s k
1
1.5 s k
3) 2)(s 1)(s (s
1
1
G
1
1
1
or
G 1
G
loop) G(closed
+ + + + +
=
+
+ + +
+
=
+
+
=
The characteristic equation of the closed loop transfer function is now:
( ) 0 1.5 s k 3) 2)(s 1)(s (s = + + + + +
Any value of k that satisfies this equation will produce instability.
We are solving the roots of ( ) 0 1.5 s k 3) 2)(s 1)(s (s = + + + + +
The result is shown below. The first thing we notice is that no value of k will produce instability as the
root loci never cross into the positive real regions. The first root starts at -3 and ends at -1.5 (the zero).
One complex root starts at -2 and moves to -3 before becoming complex. The other starts at -1 and moves
to 2.5 before becoming complex.
Figure 9
The real locus runs from -3 to 1.5. The complex loci have asymptotes at 90
o
and 270
o
and intercept the
real axis at -2.25. The affect of the zero is stabilizing the system in this case by reducing the affect of the
complex roots. Increased gain makes the system more oscillatory.
D.J .Dunn 9
6 GENERAL RULES FOR CONSTRUCTING ROOT LOCUS DIAGRAMS.
In an examination situation, it is unlikely that a student will have access to a calculator capable of solving
complex roots to give the answers listed above. Rules for general construction are given in the EC exam
and are as follows. The red comments are the authors. For some reason they use w instead of s in the
published rules.
1 Arrange the closed-loop characteristic equation into the form 1 +KF(w) =0
where F(w) =M(w)/N(w),
M(w) =(w - z
1
) (w - z
2
) ... (w z
m
)
and N(w) =(w
-
P
1
)( w -P
2
) ... (w -p
n
),
so that z
1
, z
2
, ...,z
m
, are the m zeros of F(w) and P
1
,
p
2
, . . . p
n
are the n poles of F(w).
2 For positive K every point on the root loci satisfies the magnitude and angle criteria:
|F(w) | =1/K
F(w) =(1 +2h) x 180, where h =0, 1, 2.....
3 (a) The number of separate loci is equal to n. There is one locus for each root.
(b) For m =n the n loci start at poles when K =0 and terminate at zeros
when K
(c) For m <n, and r =n - m, m of the loci start at poles when K =0 and terminate at zeros when
K, and r of the loci start at poles when K =0 and approach asymptotes as K
(d) Pairs of loci are mirrored in the real axis of the w-plane.
4 The asymptotes intersect the real axis at a single point whose coordinate is
r
z p
w
n
i
m
i
i i
o
=
= = 1 1
When they are complex use the real part in this calculation.
5 The angles of the asymptotes to the positive real axis are
( )
r
180 x 2q 1
q
+
=
where q =0, 1, 2; ..., (r- 1).
For 1 pole the angle is 180
o
.
For 2 poles the angles are 90
o
and 270
o
.
For 3 poles the angles are 60
o
and 180
o
and 300
o
For 4 poles the angles are 45
o
, 135
o
, 225
o
and 315
o
and so on.
6 The loci, or segments of loci, coincide with the real axis at points for which the total number of
poles and zeros on the axis to the right is odd.
7 The loci branch to and from the real axis at points given by the real solutions to
0 =
dw
w dN
w M
dw
w dM
N
w
) (
) (
) (
that satisfy rule 6.
8 The locus departs from a complex pole, p
i
, at an angle
Di
=180 +F'(p
i
), where
F'(w) =(w -p
i
)F(w).
9 The locus arrives at a complex zero, z
i
, at an angle
Ai
=180 - F"(z
i
)
where F" (w) =F(w)/( w -z
i
).
D.J .Dunn 10
Explanation of rule 8
CASE 1
Consider the case of 4 poles and no zeros. Two poles are complex.
Choose a point s just to the right of the complex pole as shown. Draw a line from each pole to this point.
If it is close to a pole as shown, the angle
3
is
zero. Measure or calculate the other angles.
Add the angles together and then calculate the
next largest angle as a multiple of 180
o
.
The breakaway angle of the locus from the pole
is -
p
1
=-10 p
2
=-1
p
3
= -4 +4j and p
4
=p
3
= -4 - 4j
1
=tan
-1
4/6 =33.7
o
2
=180
o
- tan
-1
4/3 =126.9
o
3
=0
o
4
=90
o
=250.6
o
The next largest multiple of 180
o
is 360
o
so the
break away angle is 360 - 250.6 =109.4
o
Figure 10
CASE 2 WITH A ZERO
p
1
=0 p
2
=-10
p
3
= -4 +4j and p
4
=p
3
= -4 - 4j
z
1
=-1
1
=180
o
- tan
-1
4/4 =45
o
2
=tan
-1
4/6 =33.7
o
3
=0
o
4
=90
o
=250.6
o
=180 - tan
-1
4/3 =126.9
o
- =126.9 - 250.6 =-123.7
o
D.J .Dunn 11
Figure 11
3 pole system given p
1
=0,
1
=180
o
- tan
-1
4/4 =180
o
- 45
o
=135
o
o
he next largest multiple of 180
o
is 360
o
so the break away angle is
Figure 12
The next multiple of 180
o
is 0
o
so the break away
angle is 0 (-123.7) =-123.7
o
CASE 3 WITH 3 POLES
Now consider the case of a
p
2
=-4 +4j and p
3
=-4 - 4j
2
=0
o
3
=90
o
=225
T
360 - 225 =135
o
D.J .Dunn 12
WORKED EXAMPLE No.1
A system has an open loop transfer function of
3) 2)(s 1)(s (s
k
G(s)
+ + +
= Sketch the root locus of the
characteristic equation.
SOLUTION
This is the same problem covered earlier with a full plot.
1 G
G
(s) G
cl
+
=
The characteristic equation is G +1 =0 so the correct form is 0 1
3) 2)(s 1)(s (s
k
= +
+ + +
One locus is along the real axis from - to -3. The others are loci are asymptotes to lines at 60
o
and
300
o
respectively that intersect the real axis at (p
1
+p
2
+p
2
)/3 =(-3 2 -1)/3 =-2
The diagrams show the simplified and actual plots.
Figure 13
WORKED EXAMPLE No.2
D.J .Dunn 13
A system has an open loop transfer function of
( )
5) 3)(s 1)(s (s
2 s k
G(s)
+ + +
+
= Sketch the root locus of
the characteristic equation.
SOLUTION
1 G
G
(s) G
cl
+
= The characteristic equation is
( )
0 1
5) 3)(s 1)(s (s
2 s k
+
+ + +
=
+
m =number of zeros =1 n =number of poles =3. r =n - m =2
Intercept ={-1 -3 - 5 - (-2)} r =-7/2 =-3.5
Angles =180 x (1+2q)/r = 180 x (1+2q)/2 =90(1+2q)
For q =0 =90
o
For q =1 =270
o
For q =2 =450
o
The poles are at -1, -3 and -5 the zero is at -2. This is enough information to sketch the root locus
Figure 14
WORKED EXAMPLE No.3
16) 4s 5)(s (0.5s
k
G(s)
2
+ + +
= A system has a forward transfer function of and a feed-back transfer
function
1) (s
H(s)
+
= . Sk
1
etch the root locus of the characteristic equation and discuss the
importance of the pole positions.
SOLUTION
1 GH
(s) G
cl
+
= The characteristic equation is
G
0 1 GH = +
( )
D.J .Dunn 14
( )
( )
( )
form. correct in the now is This 0 1
3.46j) 2 3.46j)(s 2 10)(s (s 1 s
5k
16) 4s 10)(s (s 1 s
16 4s
value. s unity obtain to 2 by bottom and op t
= +
+ + +
+ + + +
+
multiply 0 1
16) 4s 5)(s (0.5s 1 s
k
2
= +
+ + + +
quadratic t Factorise 0 1
) 10)(s (s 1 s
2
= +
+ + +
he
2k
quadratic the Factorise 0 1
5k
2
= +
There are 4 poles and no zeros m =0 n =4 r =4
The intercept is {-1 - 10 - 2 - 2} 4 =-15/4 =-3.75
Two loci start on the real axis at -1 and -10 and break away to two asymptotes. The other two loci
start at the complex poles -2 3.46j
The asymptotes for 4 loci are 45
o
, 135
o
, 225
o
and 315
o
The complex poles break away at an angle found as follows.
Figure 15
1
=180 tan
-1
(3.46/1) =106
o
2
=tan
-1
(3.46/8) =23.4
o
3
=0
4
=90
o
=219.4
o
Break away angle =360 - 219.4 =140.6
o
Now we can attempt to sketch the loci. Without further work we cannot determine the break away
point or the precise point where the loci meld with the asymptotes and the student should be aware
that creating accurate plots requires a lot of practise and experience.
Figure 16
The dominant pole is p
1
at -1 when the values of K are small. As the value of K increases the real
roots have less dominance but the complex pair of poles p
3
and p
4
becomes more dominant and at
some critical point where they cross the imaginary axis, the system becomes unstable. More work is
needed to find this point.
D.J .Dunn 15
D.J .Dunn 16
SELF ASSESSMENT EXERCISE No.1
1. Sketch the approximate root locus diagram for a closed loop system with unity feed back when the
forward loop transfer function is:
3) 2)(s (s s
k
G
+ +
=
Determine the poles and the approximate point where the locus cuts the imaginary axis. Comment on
the stability.
(0, -2, -3 and j 3 )
2. Plot the root locus with respect to k for a system with a forward transfer function
1 5 3s
k
G
2
+ +
=
s
and unity feed back. (Answer on next page)
3. Plot the root locus with respect to k for a sy forward transfer function stem with a
1) 5s s(3s
k
G
2
+ +
= and unity feed back. Comment on the difference to the result for Q2
(Answer on next page)
4) 2s 4)(s (s
k
G(s)
2
+ + +
= 4. A system with unity feed back has a forward transfer function of
Determine the largest value of k for which the system is stable.
5. A system with unity feed back has an open loop transfer function
( )
80 8s s
s 1 80
G
2
d
ol
+ +
+
=
Plot the root locus of the characteristic equation with being the variable.
D.J .Dunn 17
Answer Q2
The two loci run from poles at -0.232 ans -1.434 and converge at -0.833 running off at 90
o
. The
system is always stable.
Figure 17
Answer Q 3
There is one extra pole at 0,0 giving an intercept at 0.55 and angles of 60 and 300 . The extra s term
o o
produces the possibility of instability at high gains.
Figure 18 The Complete Plot
Figure 19 The Graphical Construction
D.J .Dunn 18
Solution Q4
Poles are t a -4, (-1 +3) and (-1 +3)
Intercept is (-4 -1 -1)/3 =-2
to the real axis Angles for 3 poles are 60
o
The roots are given by
R1( ) k ( ) 8 k
1
3
2
R2( ) k
.
1
2
( ) 8 k
1
3
2
. . .
1
2
i 3( ) 8 k
1
3
R3( ) k
.
1
2
( ) 8 k
1
3
2
. . .
1
2
i 3( ) 8 k
1
3
k =50 is the limit for stability.
D.J .Dunn 19
6. ROUTH-HURWITZ CRITERION FOR ST BILITY A
if the system is unstable. We know that stability depends on
criterion is based on this and offered here with out
arranged as a polynomial with descending orders of s. It
a
n
s
n
+a
n-1
s
n-2
+a
n-2
s
n-2
+ a
3
s
2
+a
2
s
2
+a
1
s
1
+s
o
The coefficients should then be set up as a marix as follows.
a
n
a
n-2
a
n-4
a
n-6
......
(the first, third, fifth and so on)
a
n-1
a
n-3
a
n-5
a
n-7
(the second, fourth, sixth and so on)
(a new row created from the others)
d
1
d
2
.(a new row created from the others)
We calculate
b
1
= a
n-2
(a
n
/a
n-1
) a
n-3
b
2
= a
n-4
(a
n
/a
n-1
) a
n-5
b
3
= a
n-6
(a
n
/a
n-1
) a
n-7
This is a way of visualising it.
This is a mathematical way of determining
where the roots lay on the s plane. The Routh-Hurwitz
proof. The characteristic equation should be
should have a form of:
t
b
1
b
2
b
3
(a new row created from the others)
c
1
c
2
c
3
Next we calculate
c
1
= a
n-3
(a
n-1
/ b
1
) b
2
c
2
= a
n-5
(a
n-1
/ b
1
) b
3
This is a way of visualising it.
This process is repeated until all the coefficients are found.
d
1
=b
2
- (b
1
/c
1
)c
2
hen the matrix is set up, the criterion is that the system is stable if the coefficients in the first
column are all positive. One or more negative values in the first column indicate instability.
SIMPLEST CASE
W
If the highest power of the characteristic equation is 3 the criteria may be simplified as follows.
Characteristic equation =(as
3
+bs
2
+cs +d)
If the next coefficient is negative the system is unstable and this is given by R =c - ad/b
D.J .Dunn 20
WORKED EXAMPLE No.4
A system has a characteristic equation 2s +4 s +2s s +2s 2 =0
6 5 4 3
Determine if it a stable system.
SOLUTION
Note that the coefficie could be written in full as nt for s
2
is zero since it doesnt exist. The equation
2s
6
+4 s
5
+2s
4
s
3
+0s
2
+2s 2 =0 The highest power is 6
a =2 a
4
=2 a =0 a
0
=-2
6 2
a
5
=4 a
3
=-1 a
1
=2
b
1
=a
4
5
) a /4)(-1) = 2.5 (a
6
/a
3
=2 - (2
b
2
=a
2
5
) a (2/4)(2) = -1 (a
6
/a
1
=0 -
b
3
=a
5
) a
-1
=-2 0 =-2
0
(a
6
/a
c =a (a
5
/ b
1
) b
2
=-1 - (4/2.5)(-1) =0.6
1 3
c
2
=a
1
(a
5 1
/ b ) b
3
= 2 (4/2.5)(-2) =5.2
c
3
=a
3
(a
5
/ b
1
) b
2
=-1 - (4/2.5)(-1) =0.6
d
1
=b
2
- (b
1
/c
1
)c
2
=-1 (2.5/0.6)(5.2) =-22.67
The complete matrix is like this.
2 2 0 -2
4 -1 2
2.5 -1 -2
0.6 5.2 0.6
-22.67 -2
The first column has a negative value so the system is unstable.
WORKED EXAMPLE No.5
The characteristic equation of a closed loop system is 3s
3
+s
2
+0.2s +1. Determine if the system is
stable.
SOLUTION
The highest power is 3 so use the simplified test.
a =3 b =1 c = 0.2 d =1
R =c - ad/b =0.2 - 3 x 1/1 =-2.8
ce R is negative the system is unstable. Sin
WORKED EXAMPLE No.6
A closed loop system with velocity feed back has a transfer function of
200K )s 4k 50(1 12s s
200K
G(s)
2 3
+ + + +
=
Determine the value of which makes the system stable when K =8.
SOLUTION
The characteristic equation is 12s s
2 3
+ 200K 4ka)s 50(1 + + +
a =1 b =12 c =50(1 +4K) d =200K
R =c - ad/b = 50(1 +4K) - 1 x 200K/12 Put K =8 R =50(1 +32) - 1600/12
At the limit of stability R =0 so 12 x 50(1 +32) =1600
1 +32 =1600/600 =2.6667
32 =1.6667 =1.667/32 =0.05208
WORKED EXAMPLE No.7
Find the value of k where the root locus cuts the imaginary axis for worked example No.3. where
1) (s
1
H(s)
+
=
16) 4s 5)(s (0.5s
G(s)
2
+ + +
= and
k
SOLUTION
( )
0 1
16) 4s 10)(s (s 1 s
5k
2
= +
+ + + +
The characteristic equation is
D.J .Dunn 21
( ) 10)( (s 1 s 5k + + +
( ) 0 5k 160 216s 30s 15s s
0 16) 4s s
2 3 4
2
= + + + + +
= + +
The highest power is 4
a =1 a =30 a =(160+5k
4 2 0
)
a
3
=15 a
1
=216
b = 30 - (1/15)(216) = 15.6 b = (160+
1 2
5k)
c
1
=216 - (15/15.6)(160+5k)
D.J .Dunn 22
The complete matrix is like this.
1 30 160+5k
15 216
15.6 160+5k
216 (15/15.6) (160+5k)
The only term in the first column that can be negative is the last so limit of stability is reached when
216 (15/15.6) (160+5k) =0
224.64 =(160+5k)
5k =64.64
k =13
The root locus crosses the imaginary axis when k =13
SELF ASSESSMENT EXERCISE No.2
1. n, verify that k =60 where the root locus cuts the imaginary axis Using the Routh Hurwitz criterio
for the characteristic equation s
3
+6s
2
+11s +6 +k =0 (the example used earlier)
2. Determine if the closed loop systems described by the following transfer functions are stable.
a. G(s) =1/(2s
3
+3s
2
+0.8s +1) (R =0.113)
b. G(s) =1/(s
3
+12s
2
+2s +1000) (R =-81.3)
3. A system with velocity feed back has a transfer function of
5K )s 2K 200(1 10s s
30K
G(s)
2 3
+ + + +
=
Determine the value of which makes the system stable when K =40. (-0.011)
4. A system with velocity feed back has a transfer function of
20K )s K 500(1 10s 2s
3
10K
G(s)
2
+ + + +
=
Determine the value of which m akes the system stable when K=4. (-0.242)