Unformatted text preview: RootLocus Design
The rootlocus can be used to determine the value of the loop gain K , which results in a
satisfactory closedloop behavior. This is called the proportional compensator or
proportional controller and provides gradual response to deviations from the set point.
There are practical limits as to how large the gain can be made. In fact, very high gains
lead to instabilities. If the rootlocus plot is such that the desired performance cannot be
achieved by the adjustment of the gain, then it is necessary to reshape the rootloci by
adding the additional controller Gc ( s ) to the openloop transfer function. Gc ( s ) must be
chosen so that the rootlocus will pass through the proper region of the s plane. In many
cases, the speed of response and/or the damping of the uncompensated system must be
increased in order to satisfy the specifications. This requires moving the dominant
branches of the root locus to the left.
The proportional controller has no sense of time, and its action is determined by the
present value of the error. An appropriate controller must make corrections based on the
past and future values. This can be accomplished by combining proportional with integral
action PI or proportional with derivative action PD . One of the most common controllers
available commercially is the PID controller. Different processes are suited to different
combinations of proportional, integral, and derivative control. The control engineer's task
is to adjust the three gain factors to arrive at an acceptable degree of error reduction
simultaneously with acceptable dynamic response. The compensator transfer function is KI
+ KDs
s
For PD or PI controllers, the appropriate gain is set to zero.
Gc ( s ) = K P + (1) Other compensators, are lead, lag, and leadlag compensators. A firstorder compensator
having a single zero and pole in its transfer function is s + Z0
s + P0
The pole and zero are located in the left half splane as shown in Figure 1.
s1
s1
Gc ( s ) = × − p0 θp 0 θz θz 0 0 θp ×
− p0 0 − z0
− z0
(a) Phaselead
(b) Phaselag
Figure 1 Compensator phase angle contribution 7A.1 (2) For a given s1 = σ 1 + jω1 , the transfer function angle given by θ c = (θ z0 − θ p0 ) is positive
if z0 < p0 as shown in Figure 1 (a), and the compensator is known as the phaselead
controller. On the other hand if z0 > p0 as shown in Figure 1 (b), the compensator angle
θ c = (θ z0 − θ p0 ) is negative, and the compensator is known as the phaselag controller
In general, the openloop transfer function is given by
K ( s + z1 )( s + z2 ) ( s + zm )
KG ( s ) H ( s ) =
( s + p1 )( s + p2 ) ( s + pn )
where m is the number of finite zeros and n is the number of finite poles of the loop
transfer function. If n > m , there are (n − m) zeros at infinity. The characteristic equation
of the closedloop transfer function is
1 + KG ( s ) H ( s ) = 0
Therefore
( s + p1 )( s + p2 ) ( s + pn )
= −K
( s + z1 )( s + z2 ) ( s + zm )
From the above expression, it follows that for a point in the s plane to be on the rootlocus, when 0 < K < ∞ , it must satisfy the following two conditions.
K=  s + p1  s + p2 
 s + z1  s + z2   s + pn 
 s + zm  or
(3) product of vector lengths from finite poles
K=
product of vector lengths from finite zeros
and ∑ of zeros of G (s) H (s) − ∑ angle of poles of G (s) H (s) = r (180), r = ±1, ±3, or
m n ∑θ zi − ∑θ pi = 180r ,
i =1 r = ±1, ±3, (4) i =1 The magnitude and angle criteria given by (3) and (4) are used in the graphical rootlocus
design.
In addition to the MATLAB control system toolbox rlocus(num, den) for root locus plot,
MATALB control system toolbox contain the following functions which are useful for
interactively finding the gain at certain pole locations and intersect with constant ω n
circles. These are:
sgrid generates a grid over an existing continuous splane root locus or polezero map.
Lines of constant damping ratio ζ and natural frequency ω n are drawn. sgrid('new') clears
the current axes first and sets hold on. 7A.2 sgrid(Z, Wn) plots constant damping and frequency lines for the damping ratios in the
vector Z and the natural frequencies in the vector Wn.
[K, poles] = rlocfind(num, den) puts up a crosshair cursor in the graphics window
which is used to select a pole location on an existing root locus. The root locus gain
associated with this point is returned in K and all the system poles for this gain are
returned in poles.
rltool or sistool opens the SISO Design Tool. This Graphical User Interface allows you
to design singleinput/singleoutput (SISO) compensators by interacting with the root
locus, Bode, and Nichols plots of the openloop system.
1. Gain Factor Compensation or PController Design The proportional controller is a pure gain controller. The design is accomplished by
choosing a value K 0 , which results in a satisfactory transient response. The specification
may be either the step response damping ratio or the step response time constant or the
steadystate error. The procedure for finding K 0 is as follows:
• Construct an accurate rootlocus plot
• For a given ζ draw a line from origin at angle θ = cos −1 ζ measured from
negative real axis.
• The desired closedloop pole s1 is at the intersection of this line and the rootlocus.
• Estimate the vector lengths from s1 to poles and zeros and apply the magnitude
criterion as given by (3) to find K 0 .
Example 1 The openloop transfer function of a control system is given by
K
KGH ( s ) =
s ( s + 1)( s + 4)
(a) Obtain the gain K 0 of a proportional controller such that the damping ratio of the
closedloop poles will be equal 0.6. Obtain rootlocus, step response and the timedomain
specifications for the compensated system.
The rootlocus plot is shown in Figure 2. For ζ = 0.6 , θ = cos −1 0.6 = 53.13
The line drawn at this angle intersects the rootlocus at approximately, s1
The vector lengths from s1 to the poles are marked on the diagram 7A.3 −0.41 + j 0.56 . Figure 2 PController Design . From (1), we have K = (0.7)(0.8)(3.65) = 2.04
This gain will result in the velocity error constant of K v = 2.04
= 0.51 . Thus, the steady4 1
1
=
= 1.96 .
K v 0.51
The compensated closedloop transfer function is state error due to a ramp input is ess = 2.05
C (s)
=3
2
R( s ) s + 5s + 4s + 2.05
(b) Use the MATLAB Control System Toolbox functions rlocus and sgrid(zeta, wn) to
obtain the rootlocus and the gain K 0 for ζ = 0.6 . Also use the ltiview function to obtain
the system step response and the timedomain specifications.
The following commands
num=1;
den=[1 5 4 0];
rlocus(num, den);
hold on
sgrid(0.6, 1) % plots constant line zeta=0.6 & constant line wn=1 7A.4 result in Figure 3
Zoom in at the area of intersection, click at the intersection, hold and move the mouse at
intersection and adjust for Damping: 0.6. The gain is found to be 2.04. In addition, the
percentage overshoot and natural frequency are obtained, i.e., PO = 9.48% and
ω n = 0.697 .
To obtain the step response and timedomain specifications, we use the following
commands.
numc=2.04;
denc=[1 5 4 2.04];
T=tf(numc, denc)
ltiview('step', T) The result is shown in Figure 4. Rightclick on the LTI Viewer, use Chracteristics to
mark peak response, peak time, settling time, and rise time. From File Menu use Print to
Figure to obtain a Figure plot. 7A.5 Figure4
2. PD Compensator Design Here both the error and its derivative are used for control
Gc ( s ) = K P + K D s
(5)
or K
Gc ( s ) = K D s + P = K D ( s + z0 )
KD where
K
Z0 = P
(6)
KD
From above, it can be seen that the PD controller is equivalent to the addition of a simple
zero at Z 0 = K P / K D to the openloop transfer function, which improves the transient
response. From a different point of view, the PD controller may also be used to improve
the steadystate error, because it anticipates large errors and attempts corrective action
before they occur.
The procedure for the graphical rootlocus PD compensator design is as follows:
• Construct an accurate rootlocus plot
• From the design specifications; the desired damping ratio and time constant of the
dominant closedloop poles, obtain the desired location of the dominant closedloop poles. 7A.6 ζω n = 1 τ s1
• and θ = cos −1 ζ b = ζω n tan θ and b s1 = −ζω n + jb −ζω n ωn
θ
0 • Mark the poles and zeros of the openloop plant transfer function. Find the
location of the compensator zero Z 0 such that the angle criterion as given by (4)
is satisfied.
θ z 0 + (θ z1 + θ z 2 + ) − (θ p1 + θ p 2 + ) = −180 • Estimate the vector lengths from s1 to all poles and zeros and apply the magnitude
criterion as given by (3) to find K D . Find K P from (6) Example 2
Consider the control system shown in Figure 5. R( s )
− Gc ( s ) C ( s)
1
s ( s + 2)( s + 5)
Figure 5 (a) Assume the compensator is a simple proportional controller K , obtain all pertinent
pints for root locus and draw the rootlocus. Determine the location of the dominant poles
to have critically damped response, and find the time constant corresponding to this
location. Also determine the value of K and the corresponding time constant for
dominant poles damping ratio of 0.707. Obtain the compensated system step response.
(b) Gc ( s ) is a PD compensator. Design the compensator for the following timedomain
specifications.
•
• Dominant poles damping ratio ζ = 0.707
Dominant poles time constant τ = 0.5 second (a) First we construct the root locus
• The rootloci on the real axis are to the left of an odd number of finite poles and
zeros.
• n − m = 3 , i.e., there are three zeros at infinity.
• Three asymptotes with angles θ = 180 , and ±60 .
• The asymptotes intersect on the real axis at
finite poles of GH ( s ) − ∑ finite zeros of GH ( s ) −(2 + 5)
σa = ∑
=
− 2.33
n−m
3
• Breakaway point on the real axis is given by
7A.7 dK d 3
= ( s + 7 s 2 + 10s ) = 0
⇒
3s 2 + 14 s + 10 = 0
ds ds
The roots of this equation are s = −3.7863 , and s = −0.8804 . But s = −3.7863 is not
part of the rootlocus for K > 0 , therefore the breakaway point is at s = −0.8804 . The
Routh array gives the location of the jω axis crossing.
s3 1 10 2 7 K s ⇒ for stability 0 < K < 70 and s = ± j 3.16
s 70 − K 0
0
s0 K
The rootlocus is shown in Figure 6.
1 Figure 6
For the dominant poles to have critically damped response, the dominant poles are at the
breakaway position A, i.e., s1 = s2 = −0.8804 . The time constant and the gain K are
1
τ=
= 1.136 second
0.8804
K = (0.8804)(1.1296)(4.1196) = 4.06
For dominant poles damping ratio of 0.707, s1 is at position B. The time constant and the
gain K are
1
τ=
= 1.24 second
0.8074
K = (1.14)(1.44)(4.3) = 7.06 7A.8 (b) The PD controller design
1
1
ζω n = =
= 2 , and θ = coc −1 (0.707) = 45
τ 0.5
Therefore
s1 = −2 + j 2
The desired location of s1 requires the rootlocus to be shifted towards the left half splane, which requires the addition of zero by the PD controller as shown in Figure 7.
s1
j2
13 8
2 4.16
×
5 θz 33.69 ×
− z o x 2
o
0 135 ×0 Figure 7
The position of z0 is found by applying the angle criterion given by (4) θ z 0 − (135 + 90 + 33.69) = −180 ⇒ θ z 0 = 78.69 K
2
⇒ x = 0.4, and z0 = 2.4 = P
x
KD
The compensated openloop transfer function is
K ( s + 2.4)
Gc ( s )GH ( s ) = D
s ( s + 2)( s + 5)
The vector lengths from s1 are marked on the diagram as shown. Applying the magnitude
criterion, we have
( 8 )(2)( 13 )
KD =
= 10
4.16
KP KP
=
= 2.4
⇒ K P = 24
K D 10
Therefore, the controller transfer function is
Gc ( s ) = 24 + 10s
We use the following commands to obtain the closedloop transfer function and the step
response.
tan 78.69 = Gp = tf([0 0 1],[1 9 14])
Gc = tf([10 24],[0 1])
GpGc = series(Gp, Gc)
T = feedback(GpGc, 1)
ltiview('step', T) % Plant transfer function
% PD compensator
% Openloop transfer function
% closedloop transfer function
% obtains the step response 7A.9 The result is shown in Figure 8. Figure 8 Step response for the system of Example 2
3. PI Compensator Design The integral of the error as well as the error itself is used for control, and the compensator
transfer function is
K
Gc ( s ) = K P + I
(7)
s
or
K
K P (s + I )
K ( s + z0 )
K s + KI
KP
Gc ( s ) = P
=
=P
s
s
s
where
K
Z0 = I
KP
Integral control bases its corrective action on the cumulative error integrated over time.
The controller increases the type of system by 1 and is used to eliminate the steadystate
errors.
Example 3 For the control system shown in Figure 9 design a PI compensator for the following
specifications: 7A.10 R( s) 1
( s + 3)( s + 7) Gc ( s ) − C (s) Figure 9
• Zero steadystate error due to a step input
• A pair of dominant closedloop poles with a time constant of 0.25 seconds and a
damping ratio of 0.8.
Obtain the compensated system step response. ζω n = 1 τ = 1
= 4 , and θ = tan −1 (0.8) = 36.87
0.25 Therefore
s1 = −4 + j 4 * tan 36.87
⇒
s1 = −4 + j 3
The poles of the openloop transfer function and the controller pole at origin are marked
in Figure 10.
s1
j3 18
11.25 5 108.435 × 10 45 ×o 4 7 x 3 θz 0 − zo 143.13 ×0 Figure 10
The position of controller zero for the desired location of s1 is obtained by applying the
angle criterion given by (4) θ z 0 − (143.13 + 108.435 + 45) = −180
tan(180 − 116.56) = 3
x ⇒ ⇒ θ z 0 = 116.565 x = 1.5, and z0 = 4 − 1.5 = 2.5 Therefore
KI
= 2.5
KP
The compensated openloop transfer function is
K P ( s + 2.5)
Gc ( s )GH ( s ) =
s( s + 3)( s + 7) 7A.11 The vector lengths from s1 are marked on the diagram as shown. Applying the magnitude
criterion, we have
(5)( 10 )( 18 )
KP =
= 20
11.25
KI KI
=
= 2.5 ⇒ K I = 50
K P 20
Therefore, the controller transfer function is
50
Gc ( s ) = 20 +
s
The PI controller increases the system type from zero to 1. That is, we have a type 1
system and the steadystate error due to a step input is zero. We use the following
commands to obtain the closedloop transfer function and the step response.
Gp = tf([0 0 1],[1 10 21])
Gc = tf([20 50],[1 0])
GpGc = series(Gp, Gc)
T = feedback(GpGc, 1)
ltiview('step', T) % Plant transfer function
% PI compensator
% Openloop transfer function
% closedloop transfer function
% obtains the step response The result is shown in Figure 11. Figure 11 Step response for the system of Example 2
4. PID Compensator The PID controller is used to improve the dynamic response as well as to reduce or
eliminate the steadystate error. With a proportional controller increasing the controller
gain will reduce the rise time and the steadystate error. However, in systems of third
7A.12 order or higher, large gain will make the system unstable. Derivative action contributes
phaselead and will improve the transient response, reducing the overshoot and settling
time. The integral action increases the system type by 1 and eliminates the steadystate
error, but it may make the transient response worse. When you are designing a PID
controller, first set K P to a large value to produce a fast response without loosing
stability. Then add derivative gain K D and adjust its value to meet the transient response
specifications. If required introduce the Integral gain K I to eliminate the steadystate
error. Repeat the design and finetune the gains to obtain the desired response.
5. PhaseLead Design In the phaselead controller z0 < p0 , thus the controller contributes a positive angle to
the rootlocus angle criterion and tends to shift the rootlocus of the plant toward the left
in the splane. Since z0 < p0 , the compensator is a highpass filter. The phaselead
compensator has the same purpose as the PD compensator. It is utilized to improve the
transient response, to raise bandwidth and to increase the speed of response. A lead
compensator approximates derivative control and reduces the highfrequency noise
present in the PD compensator. The procedure or the graphical rootlocus design is as
follows:
•
•
•
•
• From the timedomain specifications obtain the desired location of the closedloop dominant poles.
Select the controller zero. Place the zero to the left of the smallest plant’s pole (or
on the pole for polezero cancellation)
Locate the compensator pole so that the angle criterion (3) is satisfied.
Determine the compensator gain K c such that the magnitude criterion (4) is
satisfied.
If the overall response rise time, overshoot and settling time is not satisfactory,
place the controller zero at a different location and repeat the design Moving the controller zero to the left away from the origin in the splane results in a
faster response with increase in overshoot. Moving the controller zero to the right
towards the origin will result in a slow response and reduces or eliminate the overshoot.
The compensator angle θ z 0 − θ p 0 must be positive Therefore, there is a limit on how far
to the left along the real axis the compensator zero may be moved and still be able to
satisfy the angle criterion.
Example 4 For the control system of Example 2 design a phase lead compensator to meet the
following timedomain specifications:
• Dominant poles damping ratio ζ = 0.707 7A.13 • ζω n = Dominant poles time constant τ = 0.5 second
1 τ = 1
= 2 , and θ = coc −1 (0.707) = 45
0.5 Therefore
s1 = −2 + j 2
The desired location of s1 requires the rootlocus to be shifted towards the left half splane, which requires the addition of phase lead controller as shown in Figure 12.
s1 j2 40 8 2 13
×
 po
8 θp 0 x × 33.69 4.0625
θz ×o
 zo 5 0 2 1.75 Figure 12 135 ×0 Let the controller zero be located at z0 = 1.75 . The position of the controller p0 is found
by applying the angle criterion given by (4)
θ 0 = 180 − tan −1 (2 / 0.25) = 97.125
97.125 − (135 + 90 + 33.69 + θ p 0 ) = −180 ⇒ θ p 0 = 18.435 2
⇒ x = 6, and z0 = 2 + 6 = 8
x
The compensated openloop transfer function is
K c ( s + 1.75)
Gc ( s )GH ( s ) =
s ( s + 2)( s + 5)( s + 8)
The vector lengths from s1 are marked on the diagram as shown. Applying the magnitude
criterion, we have
( 8 )(2)( 13 )( 40 )
Kc =
= 64
4.0625
Therefore, the controller transfer function is
64( s + 1.75)
Gc ( s ) =
( s + 8)
The controller dc gain is
(64)(1.75)
a0 = Gc (0) =
= 14
8
We use the following commands to obtain the closedloop transfer function and the step
response.
Gp = tf([0 0 0 1],[1 7 10 0])
% Plant transfer function
tan 18.435 = 7A.14 Gc = tf(64*[1 1.75],[1 8])
GpGc = series(Gp, Gc)
T = feedback(GpGc, 1)
ltiview('step', T) % PI compensator
% Openloop transfer function
% closedloop transfer function
% obtains the step response The result is shown in Figure 13. Figure 13 Step response for the system of Example 4. 6. Phaselag compensator approximate design
The lag compensator is an approximate integral control. The phaselag compensator is
used when the system transient response is satisfactory but requires a reduction in the
steadystate error. Since p0 < z0 , the compensator is a lowpass filter. It adds a negative
angle to the angle criterion and tends to shift the rootlocus to the right in the splane.
In the phaselag control, the controller poles and zeros are placed very close together, and
the combination is located relatively close to the origin of the splane. Thus, the rootloci
in the compensated system are shifted only slightly from their original locations. The
compensator contributes a magnitude of
 Gc ( s ) = K c  s1 + z0 
 s1 + p0  Kc The gain to satisfy the desired damping ratio is given by
7A.15 1
K0
For the compensated system, the magnitude criterion requires that
K 0  GH ( s1 ) = 1 ⇒  GH ( s1 )= K  GH ( s1 )  Gc ( s1 ) = 1 ⇒ K 1
Kc = 1
K0 or
Kc = K0
Gain to satisfy the desired damping ratio
=
K Gain to satisfy the desired steadystate error (7) For a given desired location of a closedloop pole s1 , the design can be accomplished by
trial and error. The procedure for approximate phaselag design is as follows:
• Obtain the rootlocus and determine the gain K 0 to satisfy the desired damping
ratio. • Determine the gain K to satisfy the desired steadstate error. • Evaluate the controller gain
K
Gain to satisfy the desired damping ratio
Kc =
=0
Gain to satisfy the desired steadystate error K • Select the controller zero z0 close to origin. • Based on the compensator DC gain of unity, K c z0
= 1 , find the controller pole
p0 p0 = K c z0
Example 5
Consider the control system shown in Figure 14. R( s)
− Gc ( s ) C (s)
130
( s + 10)( s + 30) Figure 14
(a) Assume the compensator is a simple proportional controller K , obtain all pertinent
pints for root locus and draw the rootlocus. Determine the gain K 0 for the step response
damping ratio of 0.8. Obtain the steadystate error and the system step response. 7A.16 • The rootloci on the real axis are to the left of an odd number of finite poles and
zeros. • n − m = 2 , i.e., there are two zeros at infinity. • Two asymptotes with angles θ = ±90 . • The asymptotes intersect on the real axis at
finite poles of GH ( s ) − ∑ finite zeros of GH ( s ) −(10 + 30)
=
− 20
σa = ∑
n−m
2
Breakaway point on the real axis is given by
dK d 2
= ( s + 40s + 300) = 0
⇒
2s + 40 = 0
ds ds
Therefore the breakaway point is at s = −20 . • The rootlocus is shown in Figure 15 Figure 15
For ζ = 0.8 ⇒ θ = cos −1 (0.8) = 36.87
The intersection of the line drawn from origin at this angle with root locus gives the
desired complex pole s1 = −20 + j15 . Applying the magnitude criterion (3), the gain K 0
is found
130 K 0 = 325 325
⇒ K 0 = 2.5
The position error constant is
7A.17 (130)(2.5)
= 1.08333
(10)(30)
The steadystate error is
1
1
ess =
=
= 0.48
1 + K p 1 + 1.08333
The step response is shown in Figure 16.
Kp = Figure 16 The step response for Example 5 (a).
(b) It is required to have approximately the same dominant closedloop pole locations and
the same damping ratio ( ζ = 0.8 ) as in part (a). Design a phaselag compensator such
that the steadystate error due to a unit step input ess will be equal to 0.0845. Obtain the
step response, and the timedomain specifications for the compensated system. The gain K , which results in ess = 0.0845 is given by
1
130 K
ess = 0.0845 =
⇒
K p = 10.8343 =
1+ K p
(10)(30)
Thus the gain to realize the steadystate error specification is K = 25
Using the approximate method, the controller gain is given by
Kc = Gain to satisfy the desired damping ratio
2.5
=
= 0.1
Gain to satisfy the desired steadystate error 25 Next choose a small value for the compensator zero, e.g., z = 1.5 7A.18 Based on the controller dc gain of unity K c z0 / p0 = 1 , the controller pole is found
p0 = K c z0 = (0.1)(1.5) = 0.15 . Thus the controller transfer function is
0.1( s + 1.5)
Gc ( s ) =
( s + 0.15)
and the compensated openloop transfer function is
Gc ( s ) KG ( s ) = (0.1)( s + 1.5)(130)(25)
325s + 478.5
=3
( s + 0.15)( s + 10)( s + 30) s + 40.15s 2 + 306 s + 45 The compensated closedloop transfer function is
C (s)
325s + 478.5
=3
R( s ) s + 40.15s 2 + 631s + 532.5
The compensated characteristic equation roots are −19.63 ± j14.5 , and –0.894. The
compensated step response is shown in Figure 17. The complex poles are shifted slightly
to the left from the specified value of −20 ± j15 . Figure 17 The compensated step response for Example 5 (b).
Note that that the complex poles are located approximately in the same location as in part
(a). The steadystate error is greatly reduced, but because of the addition of the root at –
0.894, the step response rise time and settling time are increased. If a faster response is
desired, select the controller zero further to the left away from the origin. This would
move the complex pole s1 to the right further away from the specified value. 7A.19 7. Phaselead Compensator Analytical Design
The DC gain of the compensator Gc ( s ) =
a0 = Gc (0) = K c ( s + z0 )
is
( s + p0 ) K c z0
p0 (8) In the analytical design the controller dc gain a0 is specified, usually in accordance to the
steadystate error specification. Then, for a given location of the closedloop pole
s1 = s1  ∠β ,
z0 , and p0 are obtained such that the equation
1 + Gc ( s1 )GH ( s1 ) = 0
is satisfied. It can be shown that the above parameters are found from the following
equations
z0 = a0
,
a1 a1 = sin β + a0 M sin( β −ψ )
s1 M sinψ p0 = 1
,
b1 and Kc = a0 p0
z0 (9) where (10) sin( β + ψ ) + a0 M sin β
b1 = −
s1 sinψ
where M and ψ are the magnitude and phase angle of the openloop plant transfer
function evaluated at s1 , i.e.,
GH ( s1 ) = M ∠ψ (11) For the case that ψ is either 0 or 180 , (10) is given by
a1 s1 cos β ± b1  s1 
1
+ a0 = 0
cos β ±
M
M (12) where the plus sign applies for ψ = 0 and the minus sign applies for ψ = 180 . For this
case the zero of the compensator must also be assigned. 7A.20 8. PID Compensator Analytical Design
For a desired location of the closedloop pole s1 , as given by (3), the following equations
are obtained to satisfy
− sin( β + ψ ) 2 K I cos β
KP =
−
M sin β
 s1 
(13)
KI
sinψ
KD =
+
 s1  M sin β  s1 2
For PD or PI controllers, the appropriate gain is set to zero. The above equations can be
used only for the complex pole s1 . For the case that s1 is real, the zero of the PD
controller ( z0 = K P / KD) and the zero of the PI controller ( z0 = K I / K P ) are specified
and the corresponding gains to satisfy angle and magnitude criteria are obtained
accordingly. For the PID design, the value of K I to achieve a desired steady state error is
specified. Again, (13) is applied only for the complex pole s1 . 9. GUI program for rootlocus compensator design (rldesigngui)
Based on the above equations, a Graphical User Interface program has been developed
for the design of a firstorder controller in the forward path of a closedloop control
system for proportional, phaselag, phaselead, PD, PI, and PID controllers. The GUI
program named “rldesigngui”, which has the following options, can invoke these
programs: Pushbutton P Controller – This option is used for the design of gain factor
compensation. K 0 is obtained for the specified damping ratio ζ .
Pushbutton Phase Lag Controller – This option is used for the design of a phaselag
K
controller using the approximated method, K c = 0 . Gc ( s ) is designed for a desired
K
damping ratio ζ and the gain K required for the steadystate error specification. The user
must estimate the compensator zero. z0 is selected far away from s1 and close to origin.
Pushbutton Phase Lead Controller – This option is used for the design of a phaselead
controller for a desired location of the dominant complex closed loop poles. The DC gain
Kz
of the controller Gc (0) must be specified. Gc (0) = c 0 is found from the steadystate
p0
error requirement.
Pushbutton PD Controller – This option is used for the design of a PD controller for a
desired location of the dominant complex closed loop poles. 7A.21 Pushbutton PI Controller – This option is used for the design of a PI controller for a
desired location of the dominant complex closed loop poles.
Pushbutton PID Controller – This option is used for the design of a PID controller for a
desired location of the dominant complex closed loop poles. The integral gain K I must
be specified.
For each case the open loop and the closedloop compensated system transfer functions
are displayed. Also, the variables Gc (controller transfer function), Tfo (compensated
openloop transfer function), and TFc (compensated closedloop transfer function) are
sent to the workspace. For each design the pushbutton System Responses can be used to
obtain the timedomain and frequencydomain responses of the compensated system 7A.22 Example 6
Use the rldesigngui to design a phaselead controller for the system of Example 2 and
the design specifications outlined in Example 4.
The openloop transfer function of Example 2 is
1
GH ( s ) =
s ( s + 2)( s + 5)
The specification of ζ = 0.707 , and τ = 0.5 for dominant closedloop poles as specified
in Example 4 resulted in the closedloop pole location s1 = −2 + j 2 . The analytical
phaselead controller design requires the specification of the controller dc gain. This is
often obtained by specifying the steadystate error. We are going to use the controller dc
gain obtained in the Example 4, i.e., a0 = 14 .
In MATLAB set the Current Directory to the folder where rldesigngui and the related
files are located. At the MATLAB prompt type
>> rldesigngui
The following graphical window is displayed 7A.23 Enter the plant transfer function numerator and denominator coefficients. Select the
Phase Lead Controller pushbutton. This opens the phase Lead Controller Design; enter
the desired closedloop pole and the controller dc gain. Pressing the Find Gc ( s ) button, the controller transfer function, the compensated open
loop and closedloop transfer function, and the roots of the compensated characteristic
equation are obtained as shown in the Figure. Pressing the System Responses pushbutton
will activate the ltiviewer, which enables you to obtain all system response, and their
characteristics.
The phaselead controller is
64( s + 1.75)
Gc ( s ) =
( s + 8)
The compensated openloop transfer function is
64 s + 112
s + 15s 3 + 66s 2 + 80s
and the closedloop transfer function is
C (s)
64s + 112
=4
3
R( s ) s + 15s + 66 s 2 + 144s + 112
Gc ( s )GH ( s ) = 4 7A.24 The compensated step response is as shown. You can use the rldesigngui to design the controllers for the remaining examples. 7A.25 ...
View
Full
Document
This note was uploaded on 01/23/2012 for the course EE 4580 taught by Professor Gu during the Fall '10 term at LSU.
 Fall '10
 GU

Click to edit the document details