ME 325 ControlSystems Lecture RootLocusMethod
ME 325 ControlSystems Lecture RootLocusMethod
1
Subject to unit step disturbance D(s) =
s
1
eD ( ) = −
1
+ lim s →0 G1 ( s )
lim G ( s )
s →0 2
1
= −
lim s ( s + 2)
+ lim s →0 K
s →0 0.3162
1
=−
K Equivalent block diagram
with D(s) as input
The steady-state error produced by the step disturbance is And − E(s) as output
inversely proportional to the DC gain of G1(s).
3/12/2024 ME 325 Control Systems
Root Locus Technique
K
T (s) =
s + 18s + 77 s + K
3 2
The root locus can be used to describe qualitatively the performance of a system as various
parameters are changed. For example, the effect of varying gain upon percent overshoot,
settling time, and peak time can be vividly displayed. The qualitative description can then be
verified with quantitative analysis.
NG ( s ) NH ( s) KN G ( s ) DH ( s )
G(s) = H (s) = T (s) =
DG ( s ) DH ( s ) DG ( s ) DH ( s ) + KN G ( s ) N H ( s )
where N and D are factored polynomials and signify numerator and denominator terms,
respectively.
KN G ( s ) DH ( s )
T (s) =
DG ( s ) DH ( s ) + KN G ( s ) N H ( s )
The poles T(s) are not immediately known without factoring the
denominator, and they are a function of K.
s−plane s−plane
M M
j j
+a
s = + j s+a
j j
s−plane s−plane
M M
j j2
−a 5
−7
s+a ( s + 7 ) s →5 + j 2
K = 1 G(s) H (s)
1 + KG ( s ) H ( s ) = 0
KG ( s ) H ( s ) = ( 2n + 1)180
0
F (s) = i =1
=
( s + p )
m
denominator's complex factors
j
j =1
j
and a pole length, ( s + pi ,) is the magnitude of the vector
j =1 drawn from the pole of F(s) at − pj to the point s.
Block diagram of a security camera with auto tracking can used to follow moving
objects automatically
where K = K1K 2
K ( s + 3)( s + 4 )
Example KG ( s ) H ( s ) =
( s + 1)( s + 2 )
K ( s + 3)( s + 4 )
Closed loop TF T (s) =
(1 + K ) s 2 + ( 3 + 7 K ) s + ( 2 + 12K )
j3
L1 L2 L3 L4
1 2 3 4
−4 −3 −2 −1
Vector representation of G(s) at −2+j3
K=
1
=
1
=
pole length K=
L3 L4
= 0.33
G(s) H (s) M zero length L1L2
G ( s ) H ( s ) = zeros − poles P1
P4 P3 P2
NG ( s ) NH ( s) KN G ( s ) DH ( s )
G(s) = H (s) = Closed loop TF: T (s) =
DG ( s ) DH ( s ) DG ( s ) DH ( s ) + KN G ( s ) N H ( s )
Clearly, this open—loop transfer function has three poles, 0, −1, −2. It has no
finite zeros. For large s, we can see that
K
KG ( s ) H ( s )
s3
a =
finite poles − finite zeros
#finite poles − #finite zeros
a =
( 2k + 1)1800 k = 0, 1, 2,
#finite poles − #finite zeros
In this example,
poles = {0, −1, − 2, − 4},
zeros = {−3},
( 2k + 1)1800
a =
4 −1
= 600 ,1800 ,3000
KG ( s ) H ( s ) = =
( )(
s + 1 s + 2 ) ( s 2 + 3s + 2)
So on the real-axis segments we have
K ( ) = −
( + 1)( + 2 )
=−
( 2 + 3 + 2 )
( − 3)( − 5) ( 2 − 8 + 15)
Taking the derivative,
dK ( ) 11 2 − 26 − 61
=−
d ( )
2
2 − 8 + 15
dK ( )
and setting =0 we find
d
1 = −1.45, 2 = 3.82
Alternatively, poles = {−1, −2},
zeros = {+3, +5} so we must solve
1 1 1 1
+ = + 11 2 − 26 − 61
( − 3) ( − 5) ( + 1) ( + 2 )
Kj + 3K
KG ( j ) H ( j ) = 4
− j 7 3 − 14 2 + j8
KG ( j ) H ( j ) = −1
− 4 + j 7 3 + 14 2 − j ( 8 + K ) − 3K = 0
2 =
(8 + K )
7 − ( 8 + K ) = 0
3
7
2
8+ K 8+ K
− + 14 − 3K = 0 K + 65K − 720 = 0
2
7 7
Angles of Departure
−1 + 2 + 3 − 4 − 5 + 6 = ( 2k + 1)1800
1 = 2 + 3 − 4 − 5 + 6 − ( 2k + 1)1800
Angles of Arrival
−1 + 2 + 3 − 4 − 5 + 6 = ( 2k + 1)1800
2 = 1 − 3 + 4 + 5 − 6 + ( 2k + 1)1800
a =
( ( sum of the poles ) − ( sum of the zeros ) )
(n − m)
( ( −1 − 2 − 4 ) − ( 0 ) ) = −
1 7
=
3 3
−7ω2 + 8 + K = 0
−ω3 + 14ω = 0
giving ω = 0, √14
K = −8, 90
Since the root locus is defined only for K > 0 we conclude that the system will become
unstable for K > 90, and the locus will cross the imaginary axis at s = j 14 rad/s
( ( −1 − 2 − 4 ) − ( −6 ) ) = −0.5
1
The centroid will be a =
2
These two cases are sketched below
AC D E
K= = 1.71
B
In summary, we search a given line for the point yielding a summation of angles (zero
angles-pole angles) equal to an odd multiple of 180°. The gain at that point is then found
by multiplying the pole lengths drawn to that point and dividing by the product of the zero
lengths drawn to that point.
d = n , d = n 1 − 2 , tan =
1− 2
We can rewrite the step response as
1
1− e− d t cos (d t − ) j
1− 2
d
=
n s − plane
d −n = − d
1− 2 =
n
0
d 1− 2
tan = =
d
cos =
− jn 1 − 2 = − jd
3/12/2024 ME 325 Control Systems
Transient response and pole location
Peak time
Tp =
n 1 − 2
Ts = −
(
ln 0.02 1 − 2 ) 4
n n
( approximation valid for 0 0.9 )
1− 2
tan =
3/12/2024 ME 325 Control Systems
%OS2
j
Tp is inversely proportional to
Tp2
the imaginary part of the pole.
%OS1
s − plane Since horizontal lines on the
s−plane are lines of constant
Tp1
imaginary value, they are also
0 lines of constant peak time
Ts2 Ts1
Settling time is inversely proportional to the real part of the pole. Since vertical lines
on the s − plane are lines of constant real value, they are also lines of constant
settling time.
Since = cos , radial lines are lines of constant . Since percent overshoot is only
a function of , radial lines are thus lines of constant percent overshoot, %OS.
As
If the given root locus does not allow the
• Rise time Tr (slower)
desired transient characteristics to be
• Settling time Ts (slower)
achieved, then we must modify the root locus
• Peak time Tp (slower)
by adding poles/zeros (compensator design)
• Overshoot %OS (smaller)
Tp =
n 1 − 2
− ln ( %OS 100 )
=
2 + ln 2 ( %OS 100 )
4
Ts =
n
Gc ( s ) =
( s + zc )
( Gc ( s ) ) = −ve
( s + pc )
zc pc
Consider 1 1
zc = , pc =
zc zc
1
s + (1 + j zc )
Gc ( s ) = z c
=
(1 + zc s ) Gc ( j ) =
1 (1 + zc s ) (1 + j zc )
s +
z
c
K 1 1
G(s) = k p = lim s →0 G ( s ) = 8.23 e() = = = 0.108
( s + 1)( s + 2 )( s + 10 ) 1 + K p 1 + 8.23
Compensate the system whose root locus is shown below, to improve the steady-state
error by a factor of 10 if the system is operating with a damping ratio of 0.174.
1 1 − e()
e() = = 0.0108 K p = = 91.59
1+ K p e()
zc
K pN = K pO K pO
pc
zc K pN 91.59
= = = 11.13
pc K pO 8.23
4 4
= = = 3.613
Ts 1.107
Just as the active ideal integral compensator can be approximated with a passive lag
network, an active ideal derivative compensator can be approximated with a passive lead
compensator. When passive networks are used, a single zero cannot be produced; rather, a
compensator zero and a pole result. However, if the pole is farther from the imaginary axis
than the zero, the angular contribution of the compensator is still positive and thus
approximates an equivalent single zero. In other words, the angular contribution of the
compensator pole subtracts from the angular contribution of the zero but does not preclude
the use of the compensator to improve transient response, since the net angular contribution
is positive, just as for a single PD controller zero.
The advantages of a passive lead network over an active PD controller are that (1) no
additional power supplies are required and (2) noise due to differentiation is reduced.
The differences are in the values of static error constants, the gain required to reach the
design point on the compensated root locus, the difficulty in justifying a second-order
approximation when the design is complete, and the ensuing transient response.