Solution 5.8.3.2
For the system of Figure 1 let
G
H
+R
C
Figure 1: Standard Closed Loop Conguration
GH =
K(s+4;j4)(s+4+j4)(s+10)
s(s+2;j2)(s+2+j2)
The rst step in drawing the root locus is to plot the poles and zeros of
GH.Thepoles of GH are not the closed loop pole locations, but they can
be used to nd the closed loop poles. The closed loop zeros can be found
immediately: they are the zeros of G and the poles of H.Thepoles and
zeros of GH serves as landmarks that help in nding the poles of the closed
loop system.
The portion of the root locus on the real axis is the shaded regions shown
in Figure 2. This region is determined byinvoking the rule that states that
root locus on the real axis is found to the left of an odd countofpoles and
zeros of GH.
The root locus has three poles and three nite zeros. All of the limbs
of the root locus will end at the nite zeros, and there are no asymptotes.
The most probable root locus is shown in Figure 3. Break-in and break-out
points seem unlikely since all the poles and zeros are close together.
Im(s)
Re(s)
X
X
X
-10
Figure 2: Root Locus on Real Axis
1
Im(s)
Re(s)
-10
X
X
X
Figure 3: Probable root locus
The angle of departure can be found using Figure 4.
1
+
2
+
3
;
1
;
2
;
3
= ;180
:
As the circle shrinks in radius all the angles except
1
can be computed.
That is,
1
=
1
+
2
+
3
+180
;
2
;
3
= tan
;1
(2=8)+ tan
;1
(6=2)+ (;tan
;1
(2=2)+ 180
[180;tan
;1
(2=2];90
= 14:03
+71:57
;45
+180
;135
;90
= ;4:4
The MATLAB dialogue
EDU>gh = zpk([-4+j*4 -4-j*4 -10],[0 -2+j*2 -2-j*2],100)
Zero/pole/gain:
100 (s+10) (s^2 + 8s + 32)
--------------------------
s(s^2+4s+8)
EDU>rlocus(gh)
2
Im(s)
Re(s)
-10
1
θ
1
θ
2
α
1
θ
3
?1
α
3
α
2
Figure 4: Calculation of angle of Departure
EDU>print -deps rl5832f.eps
EDU>
draws the root locus shown in Figure 5. The MATLAB root locus program
"rlocus' is not very good so the angle of departure appears to be positive
even though we knowitisnegative.
The MATLAB program
z1 = 4-j*4
z2 = 4+j*4
z3 = 10
p1 = 0
p2 = 2-j*2
p3 = 2+j*2
gh = zpk([-z1 -z2 -z3],[p1 -p2 -p3],1)
K=linspace(0,0.5,1600);;
[R,K] = rlocus(gh,K);;
figure(1)
plot(R,'k.')
print -deps 5832rl.eps
Draws the root locus in Figure 6. which shows the angle of departure some-
what more clearly.
3
-10 -8 -6 -4 -2 0
-5
-4
-3
-2
-1
0
1
2
3
4
5
Real Axis
Imag Axis
Figure 5: MATLAB generated root locus
4
-2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5
1.8
1.85
1.9
1.95
2
2.05
2.1
2.15
2.2
Figure 6: Angle of Departure
5