E I G H T

Root Locus Techniques

SOLUTIONS TO CASE STUDIES CHALLENGES
Antenna Control: Transient Design via Gain

a. From the Chapter 5 Case Study Challenge:

G(s) =

Since Ts = 8 seconds, we search along - , the real part of poles with this settling time, for 180o.

We find the point to be - 0.5+j6.9 with 76.39K = 7194.23, or K = 94.18. Second-order

approximation is OK since third pole is much more than 5 times further from the imaginary axis

than the dominant second-order pair.

b.
Program:
numg= 1;
deng=poly([0 -150 -1.32]);
'G(s)'
G=tf(numg,deng)
rlocus(G)
axis([-2,0,-10,10]);
title(['Root Locus'])
grid on
[K1,p]=rlocfind(G)
K=K1/76.39

Computer response:
ans =

G(s)

Transfer function:
1
-----------------------
s^3 + 151.3 s^2 + 198 s

Select a point in the graphics window

selected_point =
8-2 Chapter 8: Root Locus Techniques

-0.5034 + 6.3325i

K1 =

6.0690e+003

p=

1.0e+002 *

-1.5027
-0.0052 + 0.0633i
-0.0052 - 0.0633i

K=

79.4469

>>
ans =

G(s)

Transfer function:
1
-----------------------
s^3 + 151.3 s^2 + 198 s

Select a point in the graphics window

selected_point =

-0.5000 + 6.2269i

K1 =

5.8707e+003

p=

1.0e+002 *

-1.5026
-0.0053 + 0.0623i
-0.0053 - 0.0623i

K=

76.8521
 Solutions to Case Studies Challenges 8-3

UFSS Vehicle: Transient Design via Gain

a. Push -K1 to the right past the summing junction yielding G(s) = , where K = 0.25K1. Combine the
parallel feedback paths and obtain H(s) = (s+1). Hence, G(s)H(s) = . The root locus is shown below
in (b). Searching the 10% overshoot line ( = 0.591; ), we find the operating point to be
-1.07 ± j1.46 where K = 3.389, or K1 = 13.556.

b.
Program:
numg= [1 0.437];
deng=poly([0 -2 -1.29 -0.193]);
G=tf(numg,deng);
numh=[1 1];
denh=1;
H=tf(numh,denh);
GH=G*H;
rlocus(GH)
pos=(10);
z=-log(pos/100)/sqrt(pi^2+[log(pos/100)]^2);
sgrid(z,0)
title(['Root Locus with ' , num2str(pos), ' Percent Overshoot Line'])
[K,p]=rlocfind(GH);
pause
K1=K/0.25
T=feedback(K*G,H)
T=minreal(T)
step(T)
title(['Step Response for Design of ' , num2str(pos), ' Percent'])
8-4 Chapter 8: Root Locus Techniques

Computer response:
Select a point in the graphics window

selected_point =

-1.0704 + 1.4565i

K1 =

13.5093

Transfer function:
3.377 s + 1.476
---------------------------------------------
s^4 + 3.483 s^3 + 6.592 s^2 + 5.351 s + 1.476

Transfer function:
3.377 s + 1.476
---------------------------------------------
s^4 + 3.483 s^3 + 6.592 s^2 + 5.351 s + 1.476
 Solutions to Case Studies Challenges 8-5

ANSWERS TO REVIEW QUESTIONS

1. The plot of a system's closed-loop poles as a function of gain
2. (1) Finding the closed-loop transfer function, substituting a range of gains into the denominator, and
factoring the denominator for each value of gain. (2) Search on the s-plane for points that yield 180 degrees
when using the open-loop poles and zeros.
3. K = 1/5
4. No
5. At the zeros of G(s) and the poles of H(s)
6. (1) Apply Routh-Hurwitz to the closed-loop transfer function's denominator. (2) Search along the
imaginary axis for 180 degrees.
7. If any branch of the root locus is in the rhp, the system is unstable.
8.If the branch of the root locus is vertical, the settling time remains constant for that range of gain on the
vertical section.
9. If the root locus is circular with origin at the center
10. Determine if there are any break-in or breakaway points
11. (1) Poles must be at least five times further from the imaginary axis than the dominant second order
pair, (2) Zeros must be nearly canceled by higher order poles.
12. Number of branches, symmetry, starting and ending points
13. The zeros of the open loop system help determine the root locus. The root locus ends at the zeros. Thus,
the zeros are the closed-loop poles for high gain.

SOLUTIONS TO PROBLEMS
1.
a. No: Not symmetric; On real axis to left of an even number of poles and zeros
8-6 Chapter 8: Root Locus Techniques

b. No: On real axis to left of an even number of poles and zeros

c. No: On real axis to left of an even number of poles and zeros
d. Yes
e. No: Not symmetric; Not on real axis to left of odd number of poles and/or zeros
f. Yes
g. No: Not symmetric; real axis segment is not to the left of an odd number of poles
h. Yes

2.
j j

s-plane X
s-plane

X X O  X 

(a) (b)

j j

X s-plane X s-plane

O O  

X X

(c) (d)
 Solutions to Case Studies Challenges 8-7

j j

X s-plane X s-plane

O O  

X X

(c) (d)
j j

s-plane
s-plane

X X X X  O O X X 

(e) (f)

3.
π 3π
θ a= ,
a. There are two asymptotes with 2 2 and real axis intersection
2
0+ 0−6−(− )
3
σ a= =−2 .67
3−1 . The break-in and breakaway points are obtained by finding
2 3 2
σ (σ +6) σ +6 σ
K=− =−
2 2
σ+ σ+
3 3 . Obtaining

dK
=−
( )
2
σ + ( 3 σ 2 +12 σ )−( σ 3 +6 σ 2 )
3
=−
2 σ ( σ + 2)
2

( ) ( )
dσ 2 2 2 2
σ+ σ+
3 3 and solving for the roots of the

numerator we get: σ =0 , −2 , −2 . So the root locus looks as follows:

8-8 Chapter 8: Root Locus Techniques

( 23 )|
K s+
1+ s=−2=0
b. We can obtain K from s 2 ( s+6 ) , resulting in K=12 . Note that the open loop

T (s )=
( 32 )
12 s +

zero will appear as a closed loop zero, so the closed loop transfer function is (s+2)3

4.

Breakaway:  = -2.43 for K = 52.1

5.
The root locus is:
 Solutions to Case Studies Challenges 8-9

There are no break-in or breakaway points which can be verified by computing

N ( s ) D ( s )−N ( s ) D ( s )=0, or ( s2 +11 s+10 ) ( 2 s−2 )−( 2 s +11) ( s 2−2 s−24 ) =0, or
' '

2
13 s +68 s +244=0which has only complex conjugate roots not on the locus.

To find the range of K for closed loop stability we use the characteristic equation
K (s +1)(s +10)
1+ =0
(s +4)(s−6)
Which can be rewritten as ( 1+ K ) s 2 + ( 11 K−2 ) s + ( 10 K−24 )=0. The Routh array is

2
s 1+ K 10 K−24

s 11 K−2

1 10 K−24

2 12
The first row demands K >−1; the second K > and the third K > . The intersection
11 5
12
of the three requirements is K > =2.4 .
5
8-10 Chapter 8: Root Locus Techniques

6.

Convert the denominator to the following form: and thus identify

Plotting the root locus yields

Im Im
s-plane s-plane

Re Re
 Solutions to Case Studies Challenges 8-11

7.

Closed-loop poles will be in the left-half-plane when rhp pole reaches the origin,

or .

8.
8-12 Chapter 8: Root Locus Techniques

Closed-loop poles will be in the right-half-plane for (gain at the origin).

9 9
Therefore, stable (although marginally stable) for K< ; unstable for K > .
16 16

9.
System 1:

(a)

a. Breakaway:  = 1.41 for K = 0.03; Break-in:  = -1.41 for K = 33.97.

b. Imaginary axis crossing at j1.41 for K = 1. Thus stable for K > 1.

 Solutions to Case Studies Challenges 8-13

c. At break-in point, poles are multiple. Thus, K = 33.97.

d. Searching along 1350 line for 1800, K = 5 at 1.414 1350.

System 2:

(b)

a. Break-in:  = -1.41 for K = 28.14.

b. Imaginary axis crossing at j1.41 for K = 0.67. Thus stable for K > 0.67.

c. At break-in point, poles are multiple. Thus, K = 28.14.

d. Searching along 1350 line for 1800, K = 4 at 1.414 1350.

10.
a.

Root locus crosses the imaginary axis at the origin for K = 6. Thus the system is stable for K > 6.

b.
8-14 Chapter 8: Root Locus Techniques

Root locus crosses the imaginary axis at j0.65 for K = 0.79. Thus, the system is stable for K < 0.79.

11.

The system is closed loop unstable for all K >0 . The value of K at which there will be only
two roots in the RHP can be found by calculating the point at which the real RHP pole
crosses the origin. Namely

1+ 2
K (s+6)
(s +1)(s−2)(s+ 4) s=0 |
=0

4
It can readily be seen that there will be only two RHP roots for K < .
3
 Solutions to Case Studies Challenges 8-15

12.
System: G
Gain: 517
Root Locus for Prob. 8-16 Pole: 0.00268 + 6.3i
8 Damping: -0.000425
Overshoot (%): 100
Frequency (rad/s): 6.3

6
System: G
Gain: 74.2
Pole: -1.66 + 2.21i
4 Damping: 0.602
Overshoot (%): 9.38
Imaginary Axis (seconds )
-1

Frequency (rad/s): 2.76

-2

System: G
-4 Gain: 121
Pole: -1.37 - 3.14i
Damping: 0.401
Overshoot (%): 25.3
-6 Frequency (rad/s): 3.43

-8
-9 -8 -7 -6 -5 -4 -3 -2 -1 0 1
-1
Real Axis (seconds )

Root locus crosses the imaginary axis at j6.3 with a gain of 517. Real axis breakaway is at –
−13 1
=−4
2.0 at a gain of 36.0. Real axis intercept for the asymptotes is 3 3 . The angles of
the asymptotes are: , , . Some other points on the root locus are (see figure):
 = 0.401: -1.37 - j3.14, K = 121
 = 0.602: -1.66 + j2.21, K = 74.2

13.
a. Root locus crosses the imaginary axis at at K = 52.
8-16 Chapter 8: Root Locus Techniques

b. Since the gain is the product of pole lengths to -5, .

14.
a.

Root Locus

6 System: G
Gain: 69
Pole: -0.0387 + 2.89i
Damping: 0.0134
4 Overshoot (%): 95.9
System: G Frequency (rad/s): 2.89
Gain: 1.93
Pole: -0.492
2
Imaginary Axis (seconds )
-1

Damping: 1
Overshoot (%): 0
Frequency (rad/s): 0.492
0

-2

-4

-6

-8

-10
-8 -6 -4 -2 0 2 4
-1
Real Axis (seconds )

0−1−4−6−(−3) −8 (2 k +1)π π 5π
b. σ a= = ; θ= = ,π ,
4−1 3 3 3 3
c. The characteristic equation is s4 + 11s 3+ 34 s 2 + ( 24+ K ) s+3 K =0

The corresponding Routh array is:

4
s 1 34 3K
3
s 11 24+ K
2
350−K
s 11 3K
(−K +75)(K +112)
s 350−K

1 3K

The system will have poles on the jω axis when K=75 .

d. We calculate the breakaway points from N ( s ) D' ( s )−N ' ( s ) D ( s )=0 where N ( s )=(s +3);

D ( s )=s 4 +11s 3 +34 s2 +24 s; N ' ( s )=1 ; D' ( s )=4 s 3+ 33 s 2+ 68 s+ 24; giving
 Solutions to Case Studies Challenges 8-17

3 s +34 s + 133 s + 204 s +72=0 which has roots at −0.4915 , −2.93 ± j 1.099, −4.973.
4 3 2

The only viable breakpoint is -0.4915. We use the characteristic equation to solve for K=1.9256 .

15.

Assume that root locus is epsilon away from the asymptotes. Thus, a = ≈ -1; Angle = = , . Hence 

= 7. Checking assumption at –1 ± j100 yields -180o with K = 9997.02.

16.

Root Locus

1
Imaginary Axis (seconds )
-1

0.5

-0.5

-1

-1.5
-3 -2 -1 0 1 2 3 4
-1
Real Axis (seconds )

a. The breakaway points are found from N ( s ) D ' ( s )−N ' ( s ) D ( s )=0 or
( s2−5 s+ 6 ) ( 3 s 2+ 10 s+ 6 ) −( 2 s−5 ) ( s 3 +5 s 2+ 6 s )=0
or
4 3 2
s −10 s +13 s + 60 s+ 36=0

whose roots are 10.662, -2.5, 2.4, -0.56; the latter two being on the root locus.

b. The characteristic equation is

8-18 Chapter 8: Root Locus Techniques

K (s−2)(s−3)
1+ =0
s(s+2)( s+3)

or
3 2
s + ( 5+ K ) s + ( 6−5 K ) s+ 6 K =0

The Routh array is

3
s 1 6−5 K
2
s 3+ K 6K
−(K −1)(K + 6)
s 5+ K

1 6K

The intersection with the j-axis will occur when the third row is 0, namely K=1. At this point we

can find the intersections by solving 6 s2 +6=0; which gives jω=± j .

c.
From the Routh Hurwitz table., the last row demands K >0 , which in turn satisfies the
second row. Thus the third row will be positive if K <1. The range for closed loop
stability is 0< K < 1.

d. Searching 120, we find the point -0.366j0.637 for K=0.41

17.
 Solutions to Case Studies Challenges 8-19

a. Asymptotes: int = = - ; Angle = = , , ,

b. Breakaway: -1.38 for K = 1 and -3.62 for K = 1

c. Root locus crosses the imaginary axis at ±j2.24 for K = 126. Thus, stability for K < 126.
d. Search 0.7 damping ratio line (134.427 degrees) for 1800. Point is 1.4171 =
- 0.992 ± j1.012 for K = 10.32.
e. Without the zero, the angles to the point ±j5.5 add up to -265.074o. Therefore the contribution of

the zero must be 265.074 - 180 = 85.074o. Hence, tan 85.074o = , where - zc is the location of the
zero. Thus, zc = 0.474.
8-20 Chapter 8: Root Locus Techniques

f. After adding the zero, the root locus crosses the imaginary axis at ±j5.5 for K = 252.5. Thus, the

system is stable for K < 252.5.

g. The new root locus crosses the 0.7 damping ratio line at 2.7318134.427o for K = 11.075

compared to 1.4171134.427o for K = 10.32 for the old root locus. Thus, the new system's settling

time is shorter, but with the same percent overshoot.

18.
 Solutions to Case Studies Challenges 8-21

19.
a.

b. Root locus crosses 20% overshoot line at 1.8994 117.126o = - 0.866 ± j1.69 for K = 9.398.

c. Ts = = 4.62 seconds; Tp = = 1.859 seconds

d. Root locus crosses imaginary axis at ±j3.32 for K = 60. Therefore stability for K < 60.

e. Other poles with same gain as dominant poles: = -4.27

20.

a
8-22 Chapter 8: Root Locus Techniques

Root Locus
20

15

Imaginary Axis (seconds -1) 10

-5

-10

-15

-7 -6 -5 -4 -3 -2 -1 0 1 2
-1
Real Axis (seconds )
b.
(−5−4−3−1 ) −(2)
σ a= =−7.5
4−2
(2 k +1)π π 3 π
θ a= = ,
4−2 2 2
c. The characteristic equation is s4 + 13 s 3+( K +59)s 2+ (107−2 K ) s+(60+ 2 K )=0

The corresponding Routh array is:

4
s 1 K+59 60+2 K
3
s 13 107−2 K
2
660+15 K
s 13 60+2 K
−(K + 45.792)( K−44.025)
s 22+ 0.5 K

1 60+2 K

So the system is closed loop stable when 0< K < 44.0252.

d. The breakaway points are found from N ( s ) D ' ( s )−N ' ( s ) D ( s )=0 or
5 4 3 2
2 s +7 s −44 s −147 s +166 s +334=0
 Solutions to Case Studies Challenges 8-23

Which has the roots: 4.324, 1.6573, -3.4207, -4.5698 and -1.4908. Only the last two roots
are in the root locus.

e. For 20% overshoot, the Equation yields  = 0.456. Searching along this damping ratio line, we find
the 1800 point at –0.566 + j1.09 where K = 12.5.
f. –0.566 + j1.09
g. Second-order approximation not valid because of the existence of closed-loop zeros in the rhp.

h.
Program:
s=tf('s');
G=12.5*(s^2-2*s+2)/(s^4+13*s^3+59*s^2+107*s+60);
T=G/(1+G);
step(T)
Computer response:

Step Response
0.4

0.35

0.3

0.25

0.2
Amplitude

0.15

0.1

0.05

-0.05

-0.1
0 2 4 6 8 10 12
Time (seconds)

Simulation shows over 35% overshoot and nonminimum-phase behavior. Second-order

approximation not valid.

21.
a. Draw root locus and minimum damping ratio line.
8-24 Chapter 8: Root Locus Techniques

Minimum damping ratio

145.55
o

Minimum damping ratio is  = cos (180 - 145.55) = cos 34.45o = 0.825. Coordinates at tangent point

of  = 0.825 line with the root locus is approximately –1 + j0.686. The gain at this point is 0.32.

b. Percent overshoot for  = 0.825 is 1.019%.

c. Ts = = 4 seconds; Tp = = 4.57 seconds

d. Second-order approximation is not valid because of the two zeros and no pole-zero cancellation.

22.

Since the problem stated the settling time at large values of , assume that the root locus is
−25+α −4
approximately close to the vertical asymptotes. Hence, σ a= = . Since specified T s=1
2 Ts

sec, ∝=17 . The root locus is shown below

 Solutions to Case Studies Challenges 8-25

Root Locus
50

40

30
Imaginary Axis (seconds -1)

20

10

-10

-20

-30

-40

-50
-25 -20 -15 -10 -5 0 5
-1
Real Axis (seconds )
23.
The design point is – 0.6 ± j1.2. Excluding the pole at – , the sum of angles to the design point is –
156.08o. Thus, the contribution of the pole at – is 156.08 – 180 = – 23.92o.

The following geometry applies:

1. 2
tanθ= =tan 23 . 92=0 . 4436 .
Hence, α−0 .6 Thus  = 3.305. Adding this pole at –
3.305 yields 180o at – 0.6 ± j1.2 with K = 51.8.
8-26 Chapter 8: Root Locus Techniques

24.
a.

15

10

5
Imag Axis

-5

-10

-15
-10 -8 -6 -4 -2 0 2
Real Axis

b. Searching along the 10% overshoot line (angle = 126.239o), the point - 0.7989 + j1.0898 yields

180o for K = 81.74.

c. Higher-order poles are located at approximately –6.318 and –7.084. Since these poles are more than
5 times further from the imaginary axis than the dominant pole found in (b), the second-order
approximation is valid.
d. Searching along the imaginary axis yields 180o at j2.53, with K = 394.2.
Hence, for stability, 0 < K < 394.2.

25.
a. For a peak time of 1s, search along the horizontal line, Im = Tp to find the point of

intersection with the root locus. The intersection occurs at –2 ± j at a gain of 11.
 Solutions to Case Studies Challenges 8-27

10

Imag Axis
0

-2

-4

-6

-8

-10
-5 -4 -3 -2 -1 0 1 2
Real Axis


b.
Program:
numg=11*[1 4 5];
deng=conv([1 2 5],poly([-3 -4]));
G=tf(numg,deng);
T=feedback(G,1);
step(T)
8-28 Chapter 8: Root Locus Techniques

26.
a.

b. Searching the j axis for 180o, we locate the point j6.29 at a gain of 447.83.
c. Searching for maximum gain between -4 and -5 yields the breakaway point, -4.36. Searching for
minimum gain between -2 and -3 yields the break-in point, -2.56.

d.
j

j2


x j1

    
x x x O O 
-6 -5 -4 -3 -2 -1

90o
x -j1
 Solutions to Case Studies Challenges 8-29

To find the angle of departure from the poles at -1±j1: - - 2 - 3 + 4 + 5 - 6 - 900

= - tan-1(1/5) - tan-1(1/4) - tan-1(1/3) + tan-1(1/2) + tan-1(1/1) - 6 - 90o = 1800 . Thus, 6 = - 242.22o

e. Searching along the  = 0.3 line ( = 180 - cos-1() = 107.458o) for 180o we locate the point

3.96107.458o = -1.188±j3.777. The gain is 127.133.

27.

a.

b. To find the j crossing we will use Routh Hurwitz. The characteristic equation can be written
as

The Routh array is

3
s
2
s

For closed loop stability, the fifth row demands , the 3d row ; satisfying these

requirements the 4th row will be positive as long as . Thus the system will have roots in the
8-30 Chapter 8: Root Locus Techniques

j axis when . Substituting this value, we get a row of zeros so the poles are located at the

solution of ; namely .
c. The breakaway points are obtained from the solution of N ( s ) D ' ( s )−N ' ( s ) D ( s )=0 or
, which has roots , and
the breakaway point occurring at the latter. At this point .

e. Searching for the crossing of the =0.3 line ( )for we locate

the point , when .

28.
a.

15

10

5
Imag Axis

-5

-10

-15
-5 -4 -3 -2 -1 0 1 2 3 4
Real Axis

b. Searching the j axis for 180o, we locate the point j1.69 at a gain of 4.249.
c. Searching between -2 and -3 for maximum gain, the breakaway is found at -2.512.
d.
 Solutions to Case Studies Challenges 8-31

j

s-plane

 
X j2


 X 
X
-3 -2 -1 2

 
X -j2

To find the angle of arrival to the zero at 2 + j2:

q 1 +q 2−q3 −q 4 −q5 −q 6 =q1 + 90o −0 o −tan −1 ( 43 )−tan ( 24 )−tan ( 25 )=180

−1 −1 o

Solving for 1, the angle of arrival is 1 = –191.5 .

0

e. The closed-loop zeros are the poles of H(s), or –1 ± j2.

f. Searching the  = ) for 180o, we find the point

= -0.6537+j1.705. The gain, K = 0.8764.

g. Higher-order poles are at –2.846 ± j1.731. These are not 5 times further than the dominant poles.
Further, there are closed-loop zeros at –1 ± j2 that are not cancelled any higher-order poles. Thus, the
second-order approximation is not valid.
8-32 Chapter 8: Root Locus Techniques

29.
Root Locus for Problem 8-39
40

30

20

Imaginary Axis (seconds -1)

10

-10

-20

-30

-40
-40 -35 -30 -25 -20 -15 -10 -5 0
Real Axis (seconds -1)
 Solutions to Design Problems 8-33

Zoomed Root Locus for 8-39

5
System: untitled1
Gain: 5.94e+03
4 Pole: 0.00403 + 2.76i
Damping: -0.00146
Overshoot (%): 100
Frequency (rad/s): 2.76
3 System: untitled1
Gain: 1.09e+03
Pole: -0.978 + 0.979i
2 0.707 Damping: 0.707
Overshoot (%): 4.33
Frequency (rad/s): 1.38
Imaginary Axis (seconds-1)

0
System: untitled1
Gain: 690
-1 Pole: -1.13 - 0.0584i
Damping: 0.999
Overshoot (%): 0
-2 0.707 Frequency (rad/s): 1.13

-3

-4

-5
-2 -1.5 -1 -0.5 0 0.5
-1
Real Axis (seconds )

a. As could be seen from the zoomed graph, the root locus crosses the imaginary axis at j2.76 with K
= 5940. Therefore, the system is stable for 0 < K < 5940.

b. Search the 0.707 damping ratio line for 180o and find –0.978 + j0.978 with K = 1090.
c. Assume critical damping where root locus breaks away from the real axis. Locus breaks away at –
1.13 with K = 690.

30.
Program:
numg=1;
deng=poly([0 -3 -7 -9]);
numh=[1 30];
denh=[1 20 225];
G=tf(numg,deng)
Gzpk=zpk(G)
H=tf(numh,denh)
rlocus(G*H)
pause
K=0:10:1e4;
rlocus(G*H,K)
sgrid(0.707,0)
axis([-2,0.5,-5,5]);
pause
for i=1:1:3;
[K,P]=rlocfind(G*H)
8-34 Chapter 8: Root Locus Techniques

end
T=feedback(K*G,H)
step(T)

Computer response:
G=
1
------------------------------
s^4 + 19 s^3 + 111 s^2 + 189 s
Continuous-time transfer function.
Gzpk =
1
-------------------
s (s+9) (s+7) (s+3)
Continuous-time zero/pole/gain model.
H=
s + 30
----------------
s^2 + 20 s + 225
Continuous-time transfer function.
Root Locus
100

80

60

40
Imaginary Axis (seconds -1)

20

-20

-40

-60

-80

-100
-100 -80 -60 -40 -20 0 20 40 60 80
Real Axis (seconds -1)

Select a point in the graphics window

 Solutions to Design Problems 8-35

Root Locus
5

2 0.707
Imaginary Axis (seconds-1)

-1

-2 0.707

-3

-4

-5
-2 -1.5 -1 -0.5 0 0.5
Real Axis (seconds -1)

selected_point =
-0.9729 + 0.9783i
K=
1.0891e+03

P=
-9.9610 +11.1994i
-9.9610 -11.1994i
-8.5608 + 1.6058i
-8.5608 - 1.6058i
-0.9782 + 0.9799i
-0.9782 - 0.9799i

Select a point in the graphics window

selected_point =
0.0030 + 2.7484i
K=
5.9000e+03
P=
-9.7971 +11.2966i
-9.7971 -11.2966i
-9.7016 + 3.2277i
-9.7016 - 3.2277i
-0.0013 + 2.7518i
-0.0013 - 2.7518i
Select a point in the graphics window
selected_point =
-1.1295 + 0.0155i
K=
8-36 Chapter 8: Root Locus Techniques

688.7886
P=
-9.9753 +11.1923i
-9.9753 -11.1923i
-8.3945 + 1.2181i
-8.3945 - 1.2181i
-1.1302 + 0.0155i
-1.1302 - 0.0155i
T=
688.8 s^2 + 1.378e04 s + 1.55e05
---------------------------------------------------------------------
s^6 + 39 s^5 + 716 s^4 + 6684 s^3 + 28755 s^2 + 4.321e04 s + 2.066e04
Continuous-time transfer function.

The step response for the critically damped case is shown below:

Step Response
8

5
Amplitude

0
0 2 4 6 8 10 12
Time (seconds)

31.

a. Search j=j5 for 180and find with .

b.

c. A settling time of 0.8 sec yield a real part of -5. Thus if the zero is at the origin

, which yields complex poles with -5 as the real part. At the design point

, .
 Solutions to Design Problems 8-37

32.
a. Searching along n = -1 for 180o, find –1 + j2.04 with K = 170.13.

b. Assume critical damping when root locus breaks away form the real axis. Searching for maximum

gain, the breakaway point is at -1.78 with K = 16.946.

33.

The damping ratio corresponding to the required percent overshoot is:

−ln 0. 205
== √ π 2+ ln2 0 . 205 = 0.45;

The following MATLAB M-file was written to plot the root locus; zoom it to find the gain K at that

value of  , then plot the closed-loop unit-step response, c(t).

numg = [1 0.02];

deng = poly([0 0 -4 -10 -25]);

G = tf(numg, deng);

rlocus(G);

title('Full Root Locus')

pause

axis ([-5, 0, -1, 6]);

z=0.45;

sgrid(z,0)

title('Root Locus Zoomed-in for Dominant Poles at a 0.45 Damping

Line')

[K1,p]=rlocfind(G);

pause

K = K1;

T=feedback(K*G,1); %T is the closed-loop TF of the system

T=minreal(T);
8-38 Chapter 8: Root Locus Techniques

step(T);

grid

ylabel('Unit-Step Response, c(t)')

title('Unit-Step Response at K = K1')

The full root-locus obtained is shown below, followed by the zoomed-in root locus and the unit-step

response with all of the requested important characteristics marked on it.

For  = 0.45, K = 2,183.8 and the dominant poles were found to be at –1.2263 ± 2.4239i.
 Solutions to Design Problems 8-39

Root Locus Zoomed-in for Dominant Poles at a 0.45 Damping Line

6
0.45

4
Imaginary Axis (seconds-1)

-1
-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0
Real Axis (seconds -1)

Full Root Locus

80

60

40
Imaginary Axis (seconds-1)

20

-20

-40

-60

-80
-100 -80 -60 -40 -20 0 20 40 60 80
Real Axis (seconds -1)
8-40 Chapter 8: Root Locus Techniques

Unit-Step Response at K = K1
1.4 System: T
Peak amplitude: 1.21
Overshoot (%): 20.5
At time (seconds): 1.41
1.2
System: T
Final value: 1
1
Unit-Step Response, c(t)

System: T
Settling time (seconds): 3.04
System: T
0.8 Rise time (seconds): 0.592

0.6

0.4

0.2

0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
Time (seconds)

34.
K (s−1)
1+ =0
a. The characteristic equation is ( s+2 )(s+ 3) or s2 +(5+ K )s+(6−K )=0

The Routh array is

s2 1 6−K
s 5+ K

1 6−K

It follows that −5< K < 6

b. The locus for K >0 is

 Solutions to Design Problems 8-41

( σ +1)( σ +2) σ 2 + 5 σ +6
K= =
c. To find the break-in, breakaway points let σ −1 σ−1 . Then
dK ( σ−1 )( 2 σ +5 )−( σ +5 σ +6 ) σ 2 −2 σ−11
2
= =
calculate
dσ ( σ −1 )2 ( σ−1 )2

The roots of the numerator are σ =−2.4641, 4.4641

The locus for K <0 is

8-42 Chapter 8: Root Locus Techniques

d. The smallest settling time for the system will occur when both roots are as far away as possible
to the left of the jω axis. This will happen when the system has identical roots at -2.4641. To find
the value of K we use the characteristic equation
K (s−1 )
1− | =0
( s+2)(s +3 ) s=−2. 4641 . Solving for K we get K=−0 .0718

K ( s−1) K
K p =Lim =−
e. The proportional error constant s→0 ( s +1)( s+2 ) 6
1 6
e ss = = =0 . 9882
Then
1+K p 6−K
4
T s≈ =1 . 6 sec
f. The system in this case is critically damped, the settling time 2. 4641 .The step
response will approximately be:

35.
1300 K
L( s )=
a. The open loop transmission is s 2−739600 .
−860+860 π 3π
σ a= =0 θ a= ,
There are two asymptotes with 2−0 and angles 2 2
2 2
σ −860 dK 2σ
K=− =− =0
To find the breakaway points let 1300 . Then dσ 1300 , so the

breakaway points occur when σ =0 . The gain at this point is given by the solution of
1300 K
1+ |s=0 =0
s 2−8602 or K=568 . 9

It is obvious from the figure that the system is unstable for all values of K >0 .
 Solutions to Design Problems 8-43

1300 K (s +200 )
L(s )=
b.The open loop transmission is (s +1000)( s 2−739600 ) .
−1000−860+860−(−200 )
σ a= =−400
There are two asymptotes with 3−1 and angles
π 3π
θ a= ,
2 2
To find the breakaway points let
(σ +1000 )(σ 2 −8602 ) σ 3 +1000 σ 2 −739 .6×103 σ−739 . 6×10 6
K=− =−
1300(σ +200 ) 1300( σ +200) .

dK 2 σ 3 +1600 σ 2 + 400×103 σ +591 .68×106

=− =0
Then dσ 1300( σ +200) The numerator has two

complex conjugate solutions, and a real σ =−928 . The gain at this point is given by the solution
1300 K (s +200 )
1+ | =0
of ( s+1000 )(s 2−8602 ) s=−928 or K=9 .25

We use Routh-Hurwitz to find the range of K for which the system is closed loop stable. Let
1300 K (s +200 )
1+ =0
( s+1000 )(s 2−8602 ) or

s3 +1000 s 2 +(1300 K −739600)s+(260000 K −739. 6×106 )=0

The Routh array is:
8-44 Chapter 8: Root Locus Techniques

s3 1 1300 K−739600

s2 1000 260000 K−739 . 6×106

s 1040 K

1 260000 K−739 . 6×106

The dominant requirement given by the fourth row: K >2844 . 6

36.
V
s+
ϕ 0. 6
G( s )= ( s )=0 . 5V 2
a. After substituting numerical values δ s −12. 25

The system’s characteristic equation is:

V
s+
0. 6
1+0 . 5V 2 =0 2 2
s −12. 25 or s +0 .5 Vs+(0 .833 V −12 .25 )=0
 Solutions to Design Problems 8-45

The Routh array is

s2 1 0 . 833 V 2−12 . 25
s 0 . 5V

1 0 . 833 V 2−12 . 25
m
V >3 . 83
From which we get s
b.The characteristic equation cannot be written in the form
1+VG eq ( s)=0
−0 .5±√−3 . 08V 2 +49
s1 , 2=
c. Solving for the two roots we get 2

A simple script that will plot the root locus is

>> v=linspace(0,4,100000);

>> s1 = (-0.5+sqrt(-3.08*v.^2+49))/2;

>> s2 = (-0.5-sqrt(-3.08*v.^2+49))/2;

>> plot(real(s1),imag(s1),real(s2),imag(s2))
8-46 Chapter 8: Root Locus Techniques

0.3

0.2

0.1

-0.1

-0.2

-0.3

-0.4
-4 -3 -2 -1 0 1 2 3 4

37.
a. The system’s characteristic equation is found by calculating det (sI− A )=0 . This results in
s 4 −12. 3415 s 3 +(54 .5414 K−256 . 9538 )s2 +1250. 2 Ks−1995 . 2 K =0

Which can be manipulated into

54 .5414 s 2 +1250 . 2 s−1995 .2 54 .5414 ( s−1. 498 )(s +24 . 42 )
1+ K =0 1+K =0
4 3
s −12 .3415 s −256 . 9538 s 2
or s 2 (s+11)( s−23 .3472 )
π 3π
θ a= ,
b.There are two asymptotes with 2 2 and real axis intersection
0+0−11+23 .3472−1. 498+24 . 42
σ a= =17 .63
4−2 .

To find the break in-breakaway points let

σ 4 −12 . 3415 σ 3 −256 . 9538 σ 2
K=−
54 .5414 σ 2 +1250 . 2 σ−1995 . 2 , then
dK ( 54 . 5414 σ 2 +1250 . 2 σ−1995 . 2 )( 4 σ 3 −37 .0245 σ 2−513 . 91 σ )
=−
dσ ( 54 . 5414 σ 2 +1250 . 2 σ−1995 . 2 )2
( σ 4 −12 .3415 σ 3 −256 . 9538 σ 2 ) ( 109 . 0828 σ +1250 .2 )
−
( 54 . 5414 σ 2 +1250 . 2 σ−1995 . 2 )2
 Solutions to Design Problems 8-47

109 . 0828 σ 5 +3077 .5 σ 4 −38840 σ 3−247380 σ 2 +1025400 σ

=−
( 54 . 5414 σ 2 +1250 . 2 σ−1995 . 2 )2
The numerator of this expression has roots at σ=0,−36.13, 11.95, −7.1, 3.0664
The root locus is:

c. For any value of K >0 there are always closed loop poles in the RHP.

38.
u
a. Using Mason’s rule it can readily be found that (with d=0) the open loop transmission from r to
uo is
( K ν τs+1 ) 1 1
KT Km
u0 τs Ls Cs K T K m ( K ν τs+1 )
= =
ur 1 1 C τs [ L(C +C L ) s 2 + K m α Cs+1 ]
1+ K m α+ 2
+ L
Ls LCs C

So the system’s characteristic equation is:

K T K m ( K ν τs+1)
1+ β =0
τs [ L(C +C L ) s2 + K m α Cs+1 ]

or
8-48 Chapter 8: Root Locus Techniques

τL( C +C L ) s3 + K m αCτs2 + [ τ + βK T K m K ν τ ] s+ βK T K m =0

b.Substituting numerical values the characteristic equation becomes

15 .2 μ(11 μ+C L )s 3 +14 . 52 μs 2 +73 . 7375 m s+2 . 39125=0

The Routh array is

s3 15 .2 μ(11 μ+C L ) 73.7375m

s2 14 . 52 μ 2.39125
36 . 347 μ (11μ +C L )−1. 07067 μ
s
−14 .52 μ
1 2.39125

To have all the first column positive it is required that

C L >−11 μF and
C L <29 . 45 mF . It follows that
0<C L <29 . 45 mF due to physical restrictions in the capacitance.

The characteristic equation can also be expressed as

90909 s 3
1+C L =0
s 3 +86842 s 2 + 4 . 4101×108 s+1 .5212×1010 or
90909 s 3
1+C L =0
( s+ 35 ) (s +5379 )( s+ 81428) .

There are no asymptotes. The break-in and breakaway points are calculated by first obtaining
σ 3 + 86842σ 2 + 4 . 4101×108 σ +1 .5212×1010
C L =−
90909 σ 3 . Then differentiating,
dC L σ 3 ( 3 σ 2 + 173684 σ + 4 . 4101×10 8 )−( σ 3 +86842 σ 2 +4 . 4101×108 σ +1. 5212×10 10) 3 σ 2
=−
dσ 90909 σ 6
 Solutions to Design Problems 8-49

86842σ 4 +7 . 893×108 σ 3 + 4 . 564×1010 σ 2

=
90909 σ 6 The roots of the numerator are

σ =0, 0, −9030 .7 , −58 .2

The value of
CL at -9030.7 is obtained from
90909 s 3
1+C L 3 |s=−9030 . 7=0
s +86842 s 2 + 4 . 4101×108 s+1 .5212×1010 giving

C L =3. 5523×10−5

The root locus is:

8-50 Chapter 8: Root Locus Techniques

39.
a. Substituting values we have
K s−0.5 0.1538
Gd ( s ) =0.003(1+ ) , and G t ( s ) =
s s+ 0.5 s(s+ 0.833)

The systems characteristic equation is

1+Gd Gt =1+4.614 × 10−4 1+ ( K

) s−0.5
s s ( s +0.833 )( s+0.5 )
=0

Which can be rewritten as:

−4
4.614 ×10 (s−0.5)
1+ K 4 3 2 −4
=0
s +1.333 s +0.417 s −2.308 ×10 s

or
−4
4.614 ×10 (s−0.5)
1+ K =0
s ( s−0.0006 ) ( s +0.5029 ) (s+ 0.8306)

There are 3 asymptotes with a real axis intersection at

0+0.0006−0.5029−0.8306−0.5
σ a= =−0.61
4−1

The breakaway points are found by defining

N ( s )=s−0.5 and D ( s )=s 4 +1.333 s 3+ 0.417 s 2−2.308 ×10−4 s . Then N ' ( s )=1 ;

D ' ( s )=4 s3 +4 s 2+ 0.834 s−2.308× 10−4

Then N ( s ) D' ( s )−N ' ( s ) D ( s )=0 becomes
4 3 2 −4
3 s + 0.667 s −1.583 s −0.417 s+1.154 ×10 =0
That has roots at s=0.7417 ,−0.694 ,−0.2703 ,0.0003 . The only solution in the locus is -

0.2703.

The roots locus is:

 Solutions to Design Problems 8-51

b. The system is unstable for all values of K, although with relatively large time constants
for small K. It is very unlikely that the driver-train system is unstable.

40.

1
Rr ¿
Cr s Rr 5
Zr = = =
eq 1 Rr ¿ C r s+ 1 0 .04 s +1
Rr+
a. Cr s &
Z cable=R cable +Lcable s=0 . 06+0 . 00005 s

Hence:
1 1 0 . 04 s+1
= =
Z cable +Z r 5 ( 0. 06+0. 00005 s )( 0 . 04 s+1 ) +5
eq 0. 06+0 .00005 s+
0 . 04 s+1

0 . 04 s+1 20000 s+500000

=
2 E−6×s2 +0. 00245 s +5. 06 s2 +1225 s+2530000
=
8-52 Chapter 8: Root Locus Techniques

V s ( s) s 2 +1225 s+2530000
G( s)= =
I s ( s) ( s 2 +1225 s+2530000 ) 0 . 008 s+ ( 2 0000 s+5 00000 )
Thus: =

s 2 +1225s+253 0000
0.008 s 3 +9.8 s 2 +40240 s +500000
=
b.The transfer function of the forward loop, KG(s):

V s ( s) K ( s 2 +1225s+253 0000 )
KG(s)= =
E v ( s) 0 . 008 s 3 +9 .8 s 2 +40240 s +5 00000

200
H ( s )=
s+ 200
The transfer function of the feedback loop is given by: .

The system’s characteristic equation is:

1+ ( K ( s 2 +1225s+253 0000 ) 200

0 . 008 s +9 .8 s +40240 s +5 00000 s+200
3 2 )(
=0 )
Or equivalently

s 4 +1425 s3 + ( 5275000+25000 K ) s 2 + ( 1. 0685×109 +3. 0625×107 K ) s + ( 1 . 25×10 10+6 . 325×1010 K )

=0

The following MATLAB M-file was written to plot the root locus for the system and to find the

required-above operational parameters and functions:

numg = [1 1225 2.53E6];

deng = [0.008 9.8 40240 500000];

G = tf(numg, deng);

numh = 200;

denh = [1 200];

H = tf(numh, denh);

rlocus(G*H);

title('Full Root Locus for DC Bus Voltage Control System')

pause
 Solutions to Design Problems 8-53

axis ([-150, 0, -150, 150]);

z=0.707;

sgrid(z,0)

title('Root Locus Zoomed-in around Dominant Poles with a 0.707

Damping Line')

[K1,p]=rlocfind(G*H);

pause

K = K1;

T=feedback(K*G,H); %T is the closed-loop TF of the system

T=minreal(T);

step(750*T);

grid

title('Step Response of DC Bus Voltage')

The first figure shown below is the full root locus for that system.
8-54 Chapter 8: Root Locus Techniques

4
x 10 Full Root Locus for DC Bus Voltage Control System
1
T(s)

0.8

0.6

0.4

0.2
Imaginary Axis

-0.2

-0.4

-0.6

-0.8

-1
-1200 -1000 -800 -600 -400 -200 0 200 400 600
Real Axis

Zooming into the locus by setting the x-axis (real-axis) limits to -150 to 0 and the y-axis (imaginary-

axis) limits to -150 to 150, we get the following Root locus with the line corresponding to a damping

ratio, ζ = 0.707. That plot was used to find the gain, K, at which the system has complex-conjugate

closed-loop dominant poles with a damping ratio ζ = 0.707.

 Solutions to Design Problems 8-55

Root Locus Zoomed-in around Dominant Poles with a 0.707 Damping Line
150
0.707 T
Dominant Poles

100

50
Imaginary Axis

-50

-100

0.707
-150
-150 -100 -50 0
Real Axis

From the locus we found that:

i. The gain, K, at which the system would have complex-conjugate closed-loop dominant poles with a
damping ratio ζ = 0.707 is: K = 1.6832;
ii. The coordinates of the corresponding point selected on the root-locus are:

-1.085 E+002 +1.090 E+002i

iii. The corresponding values of all closed-loop poles are:

p= 1.0 E+003 *

-0.6040 + 2.1602i

-0.6040 - 2.1602i

-0.1085 + 0.1090i

-0.1085 - 0.1090i
iv. The output voltage vs(t) for a step input voltage vdc-ref (t) = 750 u(t) was plotted in (c) below.
8-56 Chapter 8: Root Locus Techniques

c. The output response is shown below with the required characteristics noted on the graph and listed
below.

Step Response of DC Bus Voltage

800
T
System: T
Settling Time (sec): 0.0338
700 System: T
Peak amplitude: 721
Overshoot (%): 7.41 System: T
At time (sec): 0.0206 Final Value: 671
600
System: T
Rise Time (sec): 0.00996
DC Bus Voltage, Volts

500

400

300

200

100

0
0 0.01 0.02 0.03 0.04 0.05 0.06
Time (sec)

i. The actual percent overshoot = 7.41% and the corresponding peak time, Tp = 0.0206 sec;
ii. The rise time, Tr = 0.00996 sec, and the settling time, Ts = 0.0338 sec;
iii. The final steady-state value is 671 volts.

41.
a. The open loop transfer function has an excess of 1 pole over zeros, so there is 1
asymptote at θ=180° .
The breakaway points are calculated from N ( s ) D ' ( s )−N ' ( s ) D ( s )=0 or
( s2 +0.342 s+0.023 )( 3 s2 +0.428 s+ 0.001449 )−( 2 s+0.342 ) ( s3 +0.214 s 2+ 0.001449 )=0
which can be simplified to
4 3 2 −5
s + 0.684 s +0.140739 s +0.009844 s+3.3327 × 10 =0
The roots for this polynomial are -0.1487j0.048, -0.0036, and -0.3831; the latter two
being the breakaway points since they are on the root locus. The values for gain at these
points are found solving:
 Solutions to Design Problems 8-57

1+
−3
0.0415 k c (s +0.092)(s +0.25)
s( s+ 0.007)(s +0.207) |s=−0.0036
=0

giving k c =2.754 ×10 . And

1+
0.0415 k c (s +0.092)(s +0.25)
s( s+ 0.007)(s +0.207) |
s=−3831
=0

giving k c =15.78.
The corresponding root locus is:

b. The characteristic equation can be rewritten as

3 2
s + ( 0.214 +0.0415 k c ) s + ( 0.001449+ 0.01419 k c ) s+ 0.0009545 k c =0
Substituting k c =5.35 we get
3 2
s +0.436 s +0.0774 s +0.0051=0
The roots of this polynomial are -0.1457j0.1185 and -0.1447. They are marked in the
root locus above. Note that all roots have the same real value so the system will be the
fastest possible; as per the overshoot the complex conjugate pair results in a damping
0.1457
factor ¿ cos θ= =0.76 , with negligible overshoot.
√0.11852 +0.1457 2

42.

Ω L (s ) 250 (s 2 +1 . 2 s+12500)
G P (s )= =K M D (s )=
M ( s) s3 +8 .1 s 2 +62003 s+31250
8-58 Chapter 8: Root Locus Techniques

The MATLAB M-file is:

num1=250;

num2=[1 1.2 12500];

num = conv(num1, num2);

den=[1 8.1 62003 31250];

G=tf(num, den);

rlocus(G)

pause

K=0:0.1:5;

rlocus(G,K)

sgrid(0.4,0.2)

axis([-50,0.5,100,250]);

pause

for i=1:1:3;

[K,P]=rlocfind(G)

end

T=260*feedback(K*G,1)

step(T)

grid;

The MATLAB full root-locus obtained is shown below. After a few attempts, a step

response with an overshoot of 14.6% was achieved. The figure, shown below the root locus,

illustrates the step-response obtained with all of the requested important characteristics

marked on it.
 Solutions to Design Problems 8-59

Root Locus
300

200

Imaginary Axis (seconds -1) 100

-100

-200

-300
-250 -200 -150 -100 -50 0 50
Real Axis (seconds -1)

System: T
Peak amplitude: 298
Overshoot (%): 14.8
At time (seconds): 0.0423
Step Response
300
System: T
Settling time (seconds): 0.102

250 System: T
Final value: 259

System: T
Rise time (seconds): 0.0272
Output Angular Speed, rad/sec

200

150

100

50

0
0 0.05 0.1 0.15
Time (seconds)

This output response was obtained at

8-60 Chapter 8: Root Locus Techniques

K=

3.1319

P=

1.0e+02 *

-7.2280 + 0.0000i

-0.3414 + 1.1144i

-0.3414 - 1.1144i

T=

2.036e05 s^2 + 2.443e05 s + 2.545e09

--------------------------------------------------

s^3 + 791.1 s^2 + 6.294e04 s + 9.819e06

Continuous-time transfer function.

SOLUTIONS TO DESIGN PROBLEMS

43.

G(s) =
 Solutions to Design Problems 8-61

a. Root locus crosses the imaginary axis at ±j6.755 with 61.73K equal to 134892.8. Thus for
oscillations, K = 2185.21.
b. From (a) the frequency of the oscillations is 6.755 rad/s.
c. The root locus crosses the 20% overshoot line at 6117.126o = - 2.736 + j5.34 with 61.73K =
23323.61. Thus, K = 377.83 and Ts = = = 1.462 seconds.

44.
π 3π
θ a= ,
a. There are two asymptotes with 2 2 and real axis intersection
0+ 0−10−(−1)
σ a= =−4 . 5
3−1 . To find the breakaway and break-in points, write
2
σ (σ +10)
K=
3 . 333×10 4 ( σ +1) . The derivative of this expression is
dK (σ + 1)(3 σ 2 + 20 σ )−(σ 3 +10 σ 2 ) σ (2 σ 2 +13 σ +20 )
=− =−
dσ 3 .3333×104 (σ +1 )2 3 .3333×104 (σ + 1)2 . The denominators

roots for this expression are σ =0, −2. 5,−4

The values of K at -2.5 and -4 are calculated from

3 . 333×10 4 (s+1)
1+K |s=−2. 5=0
s 2 ( s+10 ) −4
giving K=9 .375×10 and
3 . 333×10 4 (s+1)
1+K |s=−4 =0
s 2 ( s+10 ) giving K=9 .6×10
−4
8-62 Chapter 8: Root Locus Techniques

The root locus is:

b. The line corresponding to ξ=0 . 7 must be at an angle θ=cos

−1
0. 7=45 . 57∘ so the line

must lie along the points s=a+ jb=a− ja tan θ=a(1− j1 . 0202) . The angle condition for

the root locus must be satisfied. Namely

arg
( 3 . 333×10 4 (s+1)
s2 (s +10) ) (
|s=a( 1− j1. 0202 ) =arg
3 .333×10 4 (1+a− j1 .0202 )
a 2 (1− j1 . 0202)2 (10+a− j1 . 0202) )
1. 0202 1. 0202
¿−tan−1 +tan−1 −2(−45 .57 ∘)=180∘
1+a 10+a
A numerical search gives a=−1 .136 . So the root locus and the line intersect when

s=−1. 136(1− j1 . 0202)=−1. 136+ j1 .1589 . Then the value of K can be found from
3 . 333×10 4 (s+1)
1+K |s=−1 . 136+ j1 .1589 =0
s 2 ( s+10 ) giving K=6 .05×10
−4

.45.
a. Using MATLAB and the Symbolic Math Toolbox, the open-loop expression that yields a root
locus as a function of N2 is

0.2284x107N2 (s2 + 3.772e-05s + 66.27) (s2 + 49.99s + 8789)

Gdt(s) = ________________

s(s+45.12) (s2 + 4.893s + 8.777e04)

 Solutions to Design Problems 8-63

Program:
syms s N KLSS KHSS KG JR JG tel s
numGdt=3.92*N^2*KLSS*KHSS*KG*s;
denGdt=(N^2*KHSS*(JR*s^2+KLSS)*(JG*s^2*[tel*s+1]+KG*s)
+JR*s^2*KLSS*[(JG*s^2+KHSS)*(tel*s+1)+KG*s]);
Gdt=numGdt/denGdt;
'Gdt in General Terms'
pretty(Gdt)
'Values to Substitute'
KLSS=12.6e6
KHSS=301e3
KG=668
JR=190120
JG=3.8
tel=20e-3
8-64 Chapter 8: Root Locus Techniques

numGdt=3.92*N^2*KLSS*KHSS*KG*s;
numGdt=vpa(numGdt,4);
denGdt=(N^2*KHSS*(JR*s^2+KLSS)*(JG*s^2*[tel*s+1]+KG*s)
+JR*s^2*KLSS*[(JG*s^2+KHSS)*(tel*s+1)+KG*s]);
denGdt=vpa(denGdt,4);
'Gdt with Values Substituted'
Gdt=numGdt/denGdt;
pretty(Gdt)
Gdt=expand(Gdt);
Gdt=vpa(Gdt,4);
'Gdt Different Form 1'
pretty(Gdt);
denGdt=collect(denGdt,N^2);
'Gdt Different Form 2'
Gdt=collect(Gdt,N^2);
pretty(Gdt)
[numGdt,denGdt]=numden(Gdt);
numGdt=numGdt/0.4349e10;
denGdt=denGdt/0.4349e10;
denGdt=expand(denGdt);
denGdt=collect(denGdt,N^2);
Gdt=vpa(numGdt/denGdt,4);
'Gdt Different Form 3'
pretty(Gdt)
'Putting into Form for RL as a Function of N^2 using previous results'
numGH=[1 49.99 8855 3313 582400];
denGH=[41.87 2094 0.3684e7 0.1658e9 0];
denGH=denGH/denGH(1)
GH=tf(numGH,denGH)
GHzpk=zpk(GH)
'Zeros of GH'
rootsnumGH=roots(numGH)
'Poles of GH'
rootsdenGH=roots(denGH)
K=0:1:10000;
rlocus(GH,K)
sgrid(0.5,0)
pause
axis([-10,0,-20,20])
[K,P]=rlocfind(GH)

Computer response:
ans =

Gdt in General Terms

98 2 / 2 2 2
-- N KLSS KHSS KG s / (N KHSS (JR s + KLSS) (JG s (tel s + 1) + KG s)
25 /

2 2
+ JR s KLSS ((JG s + KHSS) (tel s + 1) + KG s))

ans =

Values to Substitute

KLSS =
 Solutions to Design Problems 8-65

12600000

KHSS =

301000

KG =

668

JR =

190120

JG =

3.8000

tel =

0.0200

ans =

Gdt with Values Substituted

16 2 /
.9931 10 N s / (301000.
/

2 2 8 2
N (190100. s + .1260 10 ) (3.800 s (.02000 s + 1.) + 668. s)

13 2 2
+ .2396 10 s ((3.800 s + 301000.) (.02000 s + 1.) + 668. s))

ans =

Gdt Different Form 1

16 2 / 10 2 5 12 2 4 14 2 3
.9931 10 N s / (.4349 10 N s + .2174 10 N s + .3851 10 N s
/

14 2 2 16 2 12 5 13 4
+ .1441 10 N s + .2533 10 N s + .1821 10 s + .9105 10 s

17 3 18 2
+ .1602 10 s + .7212 10 s )

ans =

Gdt Different Form 2

8-66 Chapter 8: Root Locus Techniques

16 2 / 10 5 12 4 14 3
.9931 10 N s / ((.4349 10 s + .2174 10 s + .3851 10 s
/

14 2 16 2 18 2 12 5
+ .1441 10 s + .2533 10 s) N + .7212 10 s + .1821 10 s

13 4 17 3
+ .9105 10 s + .1602 10 s )

ans =

Gdt Different Form 3

7 2 /
.2284 10 N s / (
/

5 4 3 2 2
(1.000 s + 49.99 s + 8855. s + 3313. s + 582400. s) N

9 2 5 4 7 3
+ .1658 10 s + 41.87 s + 2094. s + .3684 10 s )

ans =

Putting into Form for RL as a Function of N^2 using previous results

denGH =

1.0e+006 *

Columns 1 through 4

0.0000 0.0001 0.0880 3.9599

Column 5

Transfer function:
s^4 + 49.99 s^3 + 8855 s^2 + 3313 s + 582400
--------------------------------------------
s^4 + 50.01 s^3 + 8.799e004 s^2 + 3.96e006 s

Zero/pole/gain:
(s^2 + 66.27) (s^2 + 49.99s + 8789)
---------------------------------------
s (s+45.12) (s^2 + 4.893s + 8.777e004)

ans =

Zeros of GH
 Solutions to Design Problems 8-67

rootsnumGH =

-24.9950 +90.3548i
-24.9950 -90.3548i
-0.0000 + 8.1404i
-0.0000 - 8.1404i

ans =

Poles of GH

rootsdenGH =

1.0e+002 *

0
-0.0245 + 2.9624i
-0.0245 - 2.9624i
-0.4512

Select a point in the graphics window

selected_point =

-3.8230 + 6.5435i

K=

51.5672

P=

-21.1798 +97.6282i
-21.1798 -97.6282i
-3.8154 + 6.5338i
-3.8154 - 6.5338i

b. From the computer response, K = 0.2284x107N2 = 49.6. Therefore, N is approximately 5/1000.

46.
a.The characteristic equation is given by:
4
s
1+ K 2 2
=0
(s +s +1)
or
( 1+ K ) s 4 +2 s3 +3 s 2 +2 s+ 1=0

The Routh table is

4
s 1+ K 3 1
8-68 Chapter 8: Root Locus Techniques

3 2 2
s
s
2
2−K 1
s 2(1−K )
2−K
1 1

Clearly for closed loop stability K <1 .

b. There is no locus on the real axis, and no asymptotes. The root locus starts at the loci of the

complex poles and ends at the zeros in the origin:

Root Locus
1.5

0.5
Imaginary Axis

-0.5

-1

-1.5
-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3
Real Axis

47.
% Parameters
Jl=10;Bl=1;k=100;Jm=2;Bm=0.5;a=0.25;%a is the location of the zero
%numerator and denominator of the open loop transfer function
numo=k*[1 a];
deno=[Jl*Jm (Jl*Bm+Jm*Bl) (k*(Jl+Jm)+Bl*Bm) k*(Bl+Bm) 0];
syso=tf(numo,deno);

%Pole-Zero map for the open loop transfer function

pzmap(syso);
%Root Locus
rlocus(syso);axis([-1 0 -0.3 0.3]);
 Solutions to Design Problems 8-69

zgrid(0.707,[]); %grid for zeta=0.707 for approx. 5% overshoot

[KD,poles]=rlocfind(syso);
%Choose the appropriate location of the poles in the window and multiply the factor by the open
loop
%transfer function
syso=KD*syso;
%Close the loop
sysc=feedback(syso,1);
figure;
%Obtain the response
step(sysc)

Root Locus

0.707
0.25

0.2

0.15

0.1

0.05
Imaginary Axis

-0.05

-0.1

-0.15

-0.2

-0.25
0.707
-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0
Real Axis

Select a point in the graphics window

selected_point =

-0.1031 + 0.0978i
8-70 Chapter 8: Root Locus Techniques

Step Response
1.4
System: sysc
Peak amplitude: 1.05
1.2 Overshoot (%): 5
At time (sec): 25.8

0.8
Amplitude

0.6

0.4

0.2

0
0 10 20 30 40 50 60
Time (sec)

Note the 5 % overshoot specified.

 Solutions to Design Problems 8-71

48.

Front Panel

Block Diagram

Details of the Case Structure

8-72 Chapter 8: Root Locus Techniques

49.

The following MATLAB M-file was written to plot the root locus for the system and to find the

required-above operational parameters and functions:

numg = poly ([-0.071-6.25i -0.071+6.25i]);

deng = poly ([-0.047 -2 -0.262+5.1i -0.262-5.1i]);

G = tf(numg, deng);

rlocus(G);

pos=(16);

z=-log(pos/100)/sqrt(pi^2+(log(pos/100))^2);

sgrid(z,0)

title(['Root Locus with ', num2str(pos), ' Percent Overshoot Line

for Synchronous Machine with Te = 0.5 sec'])

[K1,p]=rlocfind(G);

pause

K=0.936*K1;

T=feedback(K*G,1); %T is the closed-loop TF of the system

T=minreal(T);

step(T);

grid
 Solutions to Design Problems 8-73

Root Locus w ith 16 Percent Overshoot Line for Synchronous Machine w ith Te = 0.5 sec
20

15

10

5 0.504
Imaginary Axis

-5 0.504

-10

-15

-20
-3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2
Real Axis

Root Locus with 16% Overshoot Line for Synchronous Machine

a.The gain at which the system becomes marginally stable is:

K = 7.1045

b.The closed-loop poles, p, and transfer function, T(s), corresponding to a 16% overshoot
are:

p= -0.2300 + 4.9446i

-0.2300 - 4.9446i

-1.0555 + 1.8634i

-1.0555 - 1.8634i

2.633 s^2 + 0.3739 s + 102.9

T (s) = --------------------------------------------------------

s^4 + 2.571 s^3 + 29.88 s^2 + 53.81 s + 105.3

8-74 Chapter 8: Root Locus Techniques

c. The coordinates of the point selected on the root-locus were: -1.0557 + 1.8634i and the
corresponding unit-step response, with  (t) in p.u, was found.
d.MATLAB was used to plot that unit-step response and to note on that curve the
required characteristics:

Step Response
1.4 System: T
Peak amplitude: 1.13
Overshoot (%): 15.3
1.2 At time (sec): 1.46
System: T
Settling Time (sec): 5.17
1

System: T
Final Value: 0.977
0.8 System: T
Amplitude

Rise Time (sec): 0.703

0.6

0.4

0.2

0
0 2 4 6 8 10 12
Time (sec)

50.
a. The open loop transfer function can be described as:
3 2 −4
0.5232 b p (s +1.3179 s +1.2354 × 10 )
G ( s )= 3 2 −4
s + 3.373 s +2.321 s+8.715 ×10
The system has open loop poles at -0.0004, -0.9625 and -2.4101; open loop zeros at 0, -
0.001, -1.3178
 Solutions to Design Problems 8-75

Since the number of open loop poles equals the number of open loop zeros there are no
asymptotes.
To calculate the breakaway points we form N ( s ) D ' ( s )−N ' ( s ) D ( s )=0 which in this case
is
( s3 +1.3179 s2 +1.2354 × 10−4 ) ( 3 s2 +6.746 s+2.321 ) −( 3 s2 +2.6358 s +1.2375 ×10−4 ) ( s3 +3.373 s 2 +2.3
Simplifying this results in
4 3 2 −7
2.0551 s + 4.6418 s + 3.0610 s + 0.0023 s+ 1.0786× 10 =0
The solutions for this equation are -1.129j4.618, -0.0007, -0.00005033 the latter two
points being on the root locus.
From

|
3 2 −4
0.5232b p (s +1.3179 s + 1.2354 ×10 )
1+ 3 2 −4
=0
s +3.373 s +2.321 s +8.715 ×10 s=−0.00005073
5
b
b p=5.0092983 ×10 and similarly for s=-0.0007, p =2570 .
The sketch for the roots locus is:

b. The system is overdamped when 0< b p <2570 and 5.0092983 ×10 < b p <∞ .
5

c. The system is critically damped when b p=2570 and b p=5.0092983 ×10 .

5

5
d. The system is underdamped when 2570< b p <5.0092983 ×10 .
e. For larger values of b p the system has dominant closed loop poles very close to the j
axes.
8-76 Chapter 8: Root Locus Techniques

51.
a.

The breakaway points are found from N ( s ) D ' ( s )−N ' ( s ) D ( s )=0 or
s ( 2 s−b )−( s2−2 bs +b2 ) =0
Which can be simplified to s2−b2=0 , from which s=±b , only the negative solution being
on the root locus. To find the value of gain at this point we use the characteristic equation

1+
ks
(s−b)2 |
s=−b
=0

From which we get that at the breakaway point k =b .

b. The characteristic equation can be written as s2 + ( k −2 b ) s +b2=0 from the quadratic

equation solution is can be readily observed that the roots will be in the left half plane
as long as k > 2b .

c. From the quadratic equation it can also be observed that the two roots will be complex
when k =2 b. At this point the characteristic equation becomes s2 +b 2=0 the roots
are s=± jb so the system will oscillate with a frequency of b rad/sec.

52.

Substituting the values given for the capacitance, inductances and resistances and reducing the

internal loops to single blocks, yields the following feedback loop:

 Solutions to Design Problems 8-77

The transfer function of this feedback loop is:

V α ( s) 1017637 .5 (s+2200)2
G P (s)= =
M α ( s) 16445 .0 (s+114. 3 ) (s+2200)+(s+2200 )((s+220 )(s 2 +114 .3 s+17636700 ) =

1017637 . 5 ( s+2200 ) 1017637 . 5 ( s+2200 )

= =
s3 +334 .3 s 2 +17678291 s+3 .88195×10 9 ( s+219 . 9) (s 2 +114 . 4 s+17653251 . 48)

The MATLAB M-file is:

numg = [1017637.5 2238802500];

deng1 = [1 219.9];

deng2 = [1 114.4 17653251.48];

deng = conv(deng1, deng2);

G = tf(numg, deng);

rlocus(G);

title('Full Root Locus for Grid Connected VSC')

pause

axis ([-300, 0, -50, 5000]);

z=0.012;

sgrid(z,0)
8-78 Chapter 8: Root Locus Techniques

title('Root Locus Zoomed-in for Dominant Poles at a 0.012 Damping

Line')

[K1,p]=rlocfind(G);

pause

K = K1;

T=feedback(K*G,1); %T is the closed-loop TF of the system

T=minreal(T);

step(208*T);

grid

ylabel('Output Voltage of the VSC, Va(t), in volts')

title('Output Response for a Grid Connected VSC')

The full root-locus obtained is shown below. The close-up of that locus (from -300 to 0 on the real

axis and from -50 to 5000 on the imaginary axis) follows. From that close-up, K and the coordinates

of the dominant poles corresponding to ζ = 0.012 were found to be:

Selected_point = -5.0590e+01 + 4.2172e+03i; K = 0.1278;

p= 1.0e+03 *

-0.0500 + 4.2169i

-0.0500 - 4.2169i

-0.2344 + 0.0000i

The corresponding closed-loop transfer function is:

1.3e05 s + 2.86e08

T= ---------------------------------------
 Solutions to Design Problems 8-79

s^3 + 334.3 s^2 + 1.781e07 s + 4.168e09

The figure, shown below the close-up, illustrates the step-response obtained at that value of the gain

(with all of the requested important characteristics marked on it) when a step input, r(t)= 208 u(t),

volts, was applied at t = 0.

4 Full Root Locus for Grid Connected VSC
x 10
2.5

1.5

1
Imaginary Axis (seconds -1)

0.5

-0.5

-1

-1.5

-2

-2.5
-3000 -2500 -2000 -1500 -1000 -500 0 500 1000 1500
Real Axis (seconds -1)
8-80 Chapter 8: Root Locus Techniques

Root Locus Zoomed-in for Dominant Poles at a 0.012 Damping Line

5000
0.012

4500

4000

3500
Imaginary Axis (seconds-1)

3000

2500

2000

1500

1000

500

0
-300 -250 -200 -150 -100 -50 0
Real Axis (seconds -1)

53.
a. The open loop transfer function is

K (520 s +10 .3844 )

KG (s ) P(s )=−
s 3 +2 . 6817 s 2 +0 . 11s +0 . 0126
−520 K ( s+0 . 02)
¿
( s+0. 02+ j0 . 0661)( s+0. 02− j 0. 0661 )(s+2. 6419 )
There are two asymptotes with a real axis intersection given by
−2 . 6419−2(−0 . 02)−(−0 .02 ) 2 πk
σ a= −1. 33 θ a=
3−1 and angles 2 . For k =0 , θ a=0 and

for k =1 ,
θ a=π .

To obtain the breakaway points let

1 σ 3 +2 . 6817 σ 2 +0. 11σ +0 .0126
K=− =
G( σ ) H (σ ) 520 σ +10 .3844
And calculate and solve
 Solutions to Design Problems 8-81

2 3 2
dK (520 σ +10. 3844 )(3 σ +5 .3634 σ +0 .11)−520 (σ +2. 6817 σ +0 . 11σ +0 . 0126 )
=
dσ (520 σ +10. 3844 )2
1040 σ 3 +1426 σ 2 + 55. 7 σ −5 . 41
¿ =0
(520 σ + 10. 3844 )2
Giving σ =−1.33, −0.0879, 0.0446 with only the latter in the root locus. The value of K at
σ =0 . 0446 is given by:
σ 3 +2. 6817 σ 2 +0.11σ +0. 0126
K= |σ =0.0446 =6 .82×10−4
520 σ +10.3844
It was already found in Problem 6.? That the system is closed loop stable for K <2 . 04×10−4 . The

root locus is:

b. Now the open loop transfer function is:

K (520 s +10. 3844 )
KG (s ) P(s )=
s 3 +2 . 6817 s 2 +0 . 11s +0 . 0126
520 K (s +0 . 02)
¿
( s+0. 02+ j0 . 0661)( s+0. 02− j 0. 0661 )(s+2. 6419 )
( 2 k +1 ) π
θ a=
There are two asymptotes with a real axis intersection as in part a. but with angles 2
π 3π π
θ a= θ a= =−
. For k =0 , 2 and for k =1 , 2 2.

The breakaway point calculation is similar to the one in part a. giving

σ =−1.33, −0.0879, 0.0446 with the first two points in the root locus.
The value of K at σ =−1 .33 is given by:
8-82 Chapter 8: Root Locus Techniques

σ 3 +2 . 6817σ 2 +0. 11σ +0 .0126

K=− |σ=−1.33=0 . 0033
520 σ +10 .3844
The value of K at σ =−0 . 0879 is given by:
3 2
σ +2 .6817σ +0.11σ +0 .0126
K=− |σ=−0.0879 =6 .5×10−4
520 σ +10 .3844
We use Routh-Hurwitz to show that the system is closed loop stable for all K >0 . The

characteristic equation is:

520 s +10. 3844
1+ K =0
s 3 +2 .6817 s 2 +0 .11 s+ 0 .0126

or

s3 +2. 6817 s2 +0 .11 s+0 . 0126+K ( 520 s+10 .3844 )=0

or

s3 +2. 6817 s2 +(0 .11+520 K )s +(0 . 0126+10 . 3844 K )=0

The Routh Array is:

s3 1 0 . 11+520 K
s2 2 .6817 0 . 0126+10 . 3844 K
0 .2824 +1384 . 1 K
s 2. 6817
1 0 . 0126+10 . 3844 K

It can easily be verified that all the entries in the first column are positive for all K >0 .

The root locus is:

 Solutions to Design Problems 8-83

54.

a. With the speed controller configured as a proportional controller [KISC= 0 and GSC(s) =
KPSC], the open-loop transfer function is:

Expanding the denominator of this transfer function, gives:

DG ( s)=6 s2 + 3 .613 s + 0 . 05724 .

Solving for the roots shows that there are two open-loop poles: – 0.5858 and – 0.0163. Thus, the

open-loop transfer function may be re-written as:

(1)

K P ¿ 0 .11
sc
K 1=
In this equation: 6 (2)

The following MATLAB M-file was written to plot the root locus for the system and to find the

value of the proportional gain, K1, at the breakaway or break-in points.

8-84 Chapter 8: Root Locus Techniques

numg = [1 0.6];

deng = poly ([-0.0163 -0.5858]);

G = tf(numg, deng);

rlocus(G);

pos=(0);

z=-log(pos/100)/sqrt(pi^2+(log(pos/100))^2);

sgrid(z,0)

title(['Root Locus with ', num2str(pos) , ' Percent Overshoot

Line'])

[K1,p]=rlocfind(G);

pause

T=feedback(K1*G,1); %T is the closed-loop TF of the system

T=minreal(T);

step(T);

axis ([0, 8, 0, 1]);

grid

The root locus shown below was obtained. Using MATLAB tools, the gain at the break-in point was

found to be larger and, hence, would yield a faster closed-loop unit-step response. The following

repeated real poles were found, which indicated that the step response is critically damped: p = -

0.6910, - 0.6910. These poles corresponded to: K1 = 0.78 (which corresponds to KPSC = 42.54). The

closed-loop transfer function, T(s), was found to be:

0.78 s + 0.468

T(s) = -----------------------

s2 + 1.382 s + 0.4775
 Solutions to Design Problems 8-85

Therefore, it was used to find the closed-loop transfer function of the system, to plot its unit-step

response, c(t), shown below, and to find the rise-time, Tr, and settling time, Ts.

Root Locus with 0 Percent Overshoot Line

0.1
0.992 0.982 0.966 0.93 0.86 0.6

0.08

0.997
0.06

0.04
System: G
0.999 Gain: 0.78
Pole: -0.691
0.02
Imaginary Axis

Damping: 1
Overshoot (%): 0
Frequency (rad/sec): 0.691
0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1
0

-0.02
0.999
-0.04

-0.06
0.997

-0.08

0.992 0.982 0.966 0.93 0.86 0.6

-0.1
-0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1
Real Axis
8-86 Chapter 8: Root Locus Techniques

System: T
Step Response Final Value: 0.98
1

System: T
0.9 Settling Time (sec): 4.69
System: T
0.8 Rise Time (sec): 2.69

0.7

0.6
Amplitude

0.5

0.4

0.3

0.2

0.1

0
0 1 2 3 4 5 6 7 8
Time (sec)

As could be seen from the graph, these times are:

Tr = 2.69 sec and Ts = 4.69 sec

b. When integral action was added (with KISC/KPSC = 0.4), the transfer function of the
KI K P ( s+0 . 4 )
sc sc
G SC ( s )= K P + =
speed controller became: sc s s and the open-loop transfer
function obtained was:

Where or
 Solutions to Design Problems 8-87

The following MATLAB M-file was written to plot the root locus for the system and to find the

gain, K1, which could result in a closed-loop unit-step response with 10% overshoot.

numg = poly ([-0.4 -0.6]);

deng = poly ([0 -0.0163 -0.5858]);

G = tf(numg, deng);

rlocus(G);

pos=(10);

z=-log(pos/100)/sqrt(pi^2+(log(pos/100))^2);

axis ([-1, 0, -0.5, 0.5]);

sgrid(z,0)

title(['Root Locus with ', num2str(pos) , ' Percent Overshoot

Line'])

[K1,p]=rlocfind(G);

pause

T=feedback(K1*G,1); %T is the closed-loop TF of the system

T=minreal(T);

step(T);

axis ([0, 20, 0, 1.5]);

grid

The root locus shown below was obtained. Using MATLAB tools, the gain at the point selected on

the locus (- 0.275 + j 0.376) was found to be K1 = 0.526 (which corresponds to KPSC = 28.7). The

corresponding closed-loop transfer function, T(s), is:

0 . 526 s2 + 0 . 526 s+ 0 . 1262

T (s)=
s3 + 1. 128 s 2 + 0. 5355 s + 0 . 1262
8-88 Chapter 8: Root Locus Techniques

T(s) has the closed-loop poles: p = – 0.580, – 0.275 ± j 0.376 and zeros at – 0.4 & – 0.6. Thus, the

complex conjugate poles are not dominant, and hence, the output response, c(t), obtained using

MATLAB, does not match that of a second-order underdamped system. Note also that the settling

time, Ts = 15 sec, , the rise time, Tr = 2 sec, the peak time, Tp = 5.03 sec, and the overshoot is 24.5%

(higher than the 10% corresponding to the dominant poles).

Root Locus with 10 Percent Overshoot Line

0.5
0.591

0.4
System: G
0.3 Gain: 0.527
Pole: -0.275 + 0.376i
Damping: 0.591
0.2 Overshoot (%): 9.98
Frequency (rad/sec): 0.466
0.1
Imaginary Axis

-0.1

-0.2

-0.3

-0.4

0.591
-0.5
-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0
Real Axis
 Solutions to Design Problems 8-89

Step Response
1.5
System: T
Peak amplitude: 1.24
Overshoot (%): 24.5
At time (sec): 5.03

System: T
Final Value: 1
1

System: T
Output, c(t), volts

Settling Time (sec): 15

System: T
Rise Time (sec): 2

0.5

0
0 2 4 6 8 10 12 14 16 18 20
Time (sec)

It should be mentioned that since we applied 1 volt-unit-step inputs (as compared to 4 volts in the

Hybrid vehicle progressive problem in Chapter 5) in both parts (a) and (b) above, we should not be

surprised that the final (steady-state) value of output voltage of the speed transducer was 1 volt,

which corresponds to a change in car speed of only 5 km/hr.

55.
a. After substitution of the delay approximation the resulting open loop transfer
function is:

There are 2 asymptotes, using equation 8.66 we obtain that the angles are and
. The intersection with the real axis (although irrelevant given the asymptote
8-90 Chapter 8: Root Locus Techniques

angles) is

b. The break-in and breakaway points can be obtained from

where ,

, and

. After substituting and simplifying one gets

. This polynomial has solutions
at . Only first root is in the root locus so it’s a
break-in point.
c. The characteristic equation is:

or

The corresponding Routh array is:

It is readily seen that the fourth row entry is positive as long as . The third row

entry is positive as long as . Therefore the range for closed loop stability is

.
d. The system will oscillate when . At this point the auxiliary equation
becomes:

Equating to zero the equation has a solution at . So the oscillation

frequency is 0.025 rad/sec.
 Solutions to Design Problems 8-91

The resulting root locus is shown below: