Solution 5.8.3.2 For the system of Figure 1 let G H +R C Figure 1: Standard Closed Loop Con guration 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