Robust Control System
Robust Control System
1
Why Robust?
• When we design a control system, our ultimate goal is to control a
particular system in a real environment.
D(s) Disturbance
+
R(s) Prefilter + Controller Plant Y(s)
GP(s) GC(s) + G(s)
- Output
+
Sensor + N(s)
1 Noise
3
Why Feedback Control Systems?
• Decrease in the sensitivity of the system to variation in
the parameters of the process G(s).
4
Sensitivity of Control Systems to Parameter Variations
Y (s) T ( s ) / T ( s ) T / T
T (s) ; S (sensitivity)
R( s) G ( s ) / G ( s ) G / G
G(s) T G 1 G
T (s) ; STG . .
1 GH ( s) G T 1 GH 2 G / 1 GH
1
STG Less S for larger GH
1 G ( s) H ( s)
5
The sensitivity of the feedback system to changes in the
feedback element H(s) is
2
T T H G H GH
SH . .
H T 1 GH G / 1 GH 1 GH
Often we need to determine S αT , where α is a parameter
within the transfer function of a block G. Use the chain rule
T T G T / T
S S G S (System Sensitivity)
/
ri ri
S (Root sensitivity)
/
6
Robust Control Systems and System Sensitivity
A control system is robust when: it has low sensitivities, (2) it is stable
over the range of parameter variations, and (3) the performance
continues to meet the specifications in the presence of a set of
changes in the system parameters.
dT / T
System sensitivity is : ST ( is the parameter;T is the tranfer function)
d /
dri
Root Sensitivity is : Sri zeros of T ( s) are independent of the parameter
d /
n 1
S Sri .
T
i 1 s ri s 1
The root is r1 1
1
ST ; ST Sri
s 1
+
R(s) 1/s+ Y(s)
-
7
Let us examine the sensitivity of the following second-order system
K
T ( s) 2
s sk
We know from Eq. (4.12) that
T 1 s s 1
SK 2
1 GH ( s ) s s K
+
R(s) K/s(s+1) Y(s)
-
8
Example 12.1: Sensitivity of a Controlled System
10
Disturbance Signals in a Feedback Control System
11
The Steady-State Error of a Unity Feedback Control System (5.7)
• One of the advantages of the feedback system is the reduction of the steady-
state error of the system.
• The steady-state error of the closed loop system is usually several orders of
magnitude smaller than the error of the open-loop system.
• The system actuating signal, which is a measure of the system error, is denoted
as Ea(s).
Ea(s) Y(s)
R(s) G(s)
H(s)
G ( s) 1
E ( s) R( s) Y ( s) R( s) R ( s ) When H ( s) 1
1 GH ( s ) 1 G ( s )
12
Compensator
• A feedback control system that provides an optimum performance
without any necessary adjustments is rare. Usually it is important to
compromise among the many conflicting and demanding
specifications and to adjust the system parameters to provide
suitable and acceptable performance when it is not possible to
obtain all the desired specifications.
13
Root Locus Method
• The root locus is a powerful tool for designing and analyzing
feedback control systems.
• It is possible to use root locus methods for design when two or three
parameters vary. This provides us with the opportunity to design
feedback systems with two or three adjustable parameters. For
example the PID controller has three adjustable parameters.
• The root locus is the path of the roots of the characteristic equation
traced out in the s-plane as a system parameter is changed.
• Read Table 7.2 to understand steps of the root locus procedure.
• The design by the root locus method is based on reshaping the root
locus of the system by adding poles and zeros to the system open
loop transfer function and forcing the root loci to pass through
desired closed-loop poles in the s-plane.
14
The root Locus Procedure
1
K s 1
Step 1 : The characteristic equation 1 GH ( s) 1 2
1
s s 1 0
4
2 K s 2
Step 2 : The transfer function GH ( s) is written in terms of poles and zeros :1 0
s s 4
The multiplicative gain parameter is 2 K . To determine the locus of roots for the gain 0 K (Step3)
we locate the poles and zeros on the real axis.
Step 4 : The angle criterion is satisfied on the real axis between the points 0 and - 2, because the angle p1
at the origin is 180o , and the angle from the zero and pole p 2 at s - 4 is zero degrees.
The locus begins at the poles and ends at the zeros.
Step 5 : Find the number of separate loci (equal to the number of poles).
Step 6 : The root loci must be symmetrical with respect to the horizontal real axis.
Step 7 : The loci proceed to the zeros at infinity along asymptotes centered at A and with angle A .
Step 8 : Determine the point at which the locus crosses the imaginary axis.
Step 9 : Determine the breakway point on the real axis.
Step 10 : Determine the angle of departure of the locus from a pole and the angle of arrival at a zero.
15
Example
1
G (s)
z1=-3+j1 s 2 s 3
R(s) Y(s)
Controller Plant
+ GC(s) G(s)
-
j2
-z1
j1
G ( s)GC ( s ) K 3 s z1 s zˆ1
T (S ) -2 -1
1 G ( s )GC ( s) s r2 s r1 s rˆ1 -z1
16
Analysis of Robustness
System goals : Small tracking error [e(t ) r (t ) y (t )] for an input r (t )
and keep the output y (t ) small for a disturbance d (t ).
Sensor noise n(t ) must be small to r (t ) so r n
+
Prefilter + Controller Plant
Y(s)
GP(s) GC(s) + G(s)
- Output
+
Sensor + N(s)
1 Noise
17
The Design of Robust Control Systems
• The design of robust control systems is based on two tasks:
determining the structure of the controller and adjusting the
controller’s parameters to give an optimal system performance. This
design process is done with complete knowledge of the plant. The
structure of the controller is chosen such that the system’s response
can meet certain performance criteria.
18
PID Controllers
PID stands for Proportional, Integral, Derivative. One form of controller
widely used in industrial process is called a three term, or PID controller.
This controller has a transfer function:
A proportional controller (Kp) will have the effect of reducing the rise time
and will reduce, but never eliminate, the steady state error. An integral
control (KI) will have the effect of eliminating the steady-state error, but it
may make the transient response worse. A derivative control (KD) will
have the effect of increasing the stability of the system, reducing the
overshoot, and improving the transient response.
KI
GC (t ) K p KDs
s
The controller provides a proportional term, an integratio n term, and a derivative term
de(t )
u (t ) K p e(t ) K I e(t )dt K D
dt
19
Proportional-Integral-Derivative (PID) Controller
kp
+
e(t) ki/s u(t)
+
+
kis
e(t )
u (t ) K p e(t ) K I K D se(t )
s
e(t ) r (t ) y (t ) is the error between the reference signal
and the system output; K P , K I , and K D are the proportional,
integral, and derivative feedback gains, respectively.
KI
U ( s) ( K p K D s) E ( s)
s
U (s) K D s 2 K P s K I
GPID ( s )
E (s) s 20
Time- and s-domain block diagram of closed loop system
Y ( s) Gsys ( s )GPID ( s)
G (s)
R( s ) 1 Gsys ( s )G PID ( s )
21
PID and Operational Amplifiers
A large number of transfer functions may be implemented using
operational amplifiers and passive elements in the input and feedback
paths. Operational amplifiers are widely used in control systems to
implement PID-type control algorithms needed.
22
Inverting amplifier
Vo (t ) R2
VS (t ) R1 Figur
e 8.5
23
Op-amp Integrator
Z 2 ( s)
Z1 ( s )
Figure
Vout ( s ) Z 2 ( s) 1 8.30
G(s)
Vs ( s ) Z1 ( s ) RS C F s
24
Op-amp Differentiator
The operational differentiator performs the differentiation of the input signal.
The current through the input capacitor is CS dvs(t)/dt. That is the output voltage
is proportional to the derivative of the input voltage with respect to time, and
Vo(t) = _RFCS dvs(t)/dt
Z 2 (s)
Figure 8.35Z
V ( s) Z ( s) 1 (s)
G(s) o 2 _ RF C S s
VS ( s ) Z1 ( s )
25
Linear PID Controller
Z2(s)
C1 C2
R2
R1
Z1(s)
vs(t) vo(t)
R1C1 R2 C 2 1
R2 C1s 2 s
G ( s)
Vo ( s ) R C s 1 R2C2 s 1
1 1
R1C 2 R1C 2
VS ( s ) R1C 2 s s
KDs2 KPs KI R C R2 C 2 1
G PID ( s ) ; KP 1 1 ; KI ; K D R2 C1 26
s R1C 2 R1C 2
Tips for Designing a PID Controller
When you are designing a PID controller for a given system, follow the
following steps in order to obtain a desired response.
27
The popularity of PID controllers may be attributed partly to their robust
performance in a wide range of operation conditions and partly to their
functional simplicity, which allows engineers to operate them in a simple
manner.
K2 K 3 s 2 K1s K 2
GC (t ) K1 K3s
s s
K 3 s 2 as b
3
K s z1 s z 2
s s
Where a K1 / K 3 ; and b K 2 / K 3 . Accordingly the PID introduces
a transfer function with one pole at the origin and two zeros that can be
located anywhere in the left - hand s - plane
28
Root Locus
Root locus begins at the poles and ends at the zeros.
1
G(s)
s 2 s 5
Assume we use a PID controller with complex zeros, we can plot the root locus. As K3 of
the controller increases, the complex roots approaches the zero.The closed loop transfer function is
G ( s )GC ( s )GP ( s ) K 3 s z1 s zˆ1 K G ( s)
T ( s) GP ( s ) 3 P
1 G ( s )GC ( s) ( s r2 ) s r1 ( s r1 ) ( s r2 )
Because the zeros and the complex roots are approximately equal.Setting G P ( s) 1, we have
K3 K3
T ( s) ; If K 3 is large, the system will have a fast response and zero steady state error.
s r2 s K 3
j4
K3 increasing
r1
z1 j2
r2
-2
z1
r1 29
Design of Robust PID-Controlled Systems
The selection of the three coefficients of PID controllers is basically a search
problem in a three-dimensional space. Points in the search space correspond to
different selections of a PID controller’s three parameters. By choosing different
points of the parameter space, we can produce different step responses for a
step input.
The first design method uses the (integral of time multiplied by absolute error
(ITAE) performance index in Section 5.9 and the optimum coefficients of Table
5.6 for a step input or Table 5.7 for a ramp input. Hence we select the three PID
coefficients to minimize the ITAE performance index, which produces an
excellent transient response to a step (see Figure 5.30c). The design
procedure consists of the following three steps.
T
ITAE t e(t ) dt
0
30
The Three Design Steps of Robust PID-Controlled System
Y ( s) bo
T (s) n
R( s ) s bn 1s n 1 ... b1s bo
31
Input Signals; Overshoot; Rise Time; Settling Time
• Step: r(t) = A R(s) = A/s
• Ramp: r(t) = At R(s) = A/s2
• The performance of a system is measured usually in terms of step response.
The swiftness of the response is measured by the rise time, Tr, and the peak
time, Tp.
• The settling time, Ts, is defined as the time required for the system to settle
within a certain percentage of the input amplitude.
• For a second-order system with a closed-loop damping constant, we seek to
determine the time, Ts, for which the response remains within 2% of the final
value. This occurs approximately when
4
e nTs 0.02; nTs 4; Ts 4 (ωn : undamped natural frequency; ξ : damping ratio)
n
1 2
Tp (Peak Time); M p 1 e / (Peak Response)
2
n 1
/ 1 2
Percentage Overshoot (PO) 100e
32
Example 12.8: Robust Control of Temperature Using PID Controller employing
ITAE performance for a step input and a settling time of less than 0.5 seconds.
D(s)
+
R(s) + E(s) U(s) +
GP(s) GC(s) G(s) Y(s)
-
1 1 1 A 1 2
G (s) ; ess ; K p lim G ( s ) 1; ess 50%; 1
s 1 2 s 2 2s 1 s 2 2n s n2 1 K p s 0 2 2
If GC ( s ) 1, the steady - state error is 50%, and the settling time (2% criterion) is 3.2 seconds for a step input.
We desire to obtain an optimim ITAE performance for a step input for a settling time of less than 0.5 seconds.
K 3s 2 K1s K 2
Using a PID controller: GC ( s ) ; The closed - loop transfer function with GP ( s ) 1 is
s
Y ( s) GC G ( s ) K 3s 2 K1s K 2
T1 ( s )
R( s ) 1 GC G ( s ) s 3 2 K 3 s 2 1 K1 s K 2
The optimum coefficients of the characteristic equation for ITAE (Table 5.6) : s 3 1.7ωn s 2 2.15ωn2 s ωn3
33
Cont. 12.8
We need to select n in order to meet the settling time requirement. Ts 4 / n .
is unknown but near 0.8, we set n 10. Equate the denominator of the equation
to the desired equation, we obtain the three coefficients as K1 214,
K 3 15.5, and K 2 1000.
15.5s 2 214s 1000 15.5 s 6.9 j 4.1 s 6.9 j 4.1
T (s)
s 3 17.5s 2 215s 1000 s 3 17.5s 2 215s 1000
We select a prefilter GP ( s ) so that to achieve the desired ITAE response
Gc ( s )GG P ( s ) 1000
T (s)
1 GGc ( s ) s 3 17.5s 2 215s 1000
64.5
Therefore we require GP ( s )
2
s 13.8s 64.5
in order to eliminate the
zeros in the previous equation and bring the overall numerator to 1000.
34
Results for Example 12.8
35
E12.1: Using the ITAE performance method for step input, determine the required GC(t).
Assume n = 20 for Table 5.6. Determine the step response with and without a prefilter
GP(s)
D(s)
+
R(s) + E(s) U(s) +
GP(s) GC(s) G(s) Y(s)
-
1 K
G (s) ; Use a PI controller given by GC K1 2
s 1 s
GC G ( s )
T (s) Find it s 2 ( K1 1) s K 2
1 GC G ( s )
The ITAE characteristic equation is : s 2 1.4n s n2
When n 20 then we have K1 27 and K 2 400
Y (s) 27 s 400
Without a prefilter, the closed - loop system is
R( s ) s 2 28s 400
Y ( s ) GC ( s )GGp ( s ) 400
With a prefilter, the closed - loop gain is
R( s) 1 GC G ( s ) s 2 28s 400
14.8
Where GP ( s ) ; Draw the step response without and with the prefilter 36
s 14.8
E12.3 A closed-loop unity feedback system has
9 T
G (s) ; p 3; Find S P and plot T( jω) and S( jω) on a Bode plot
s s p
9
The closed - loop transfer function is T ( s )
s 2 ps 9
s 2 ps
The sensitivity function is S ( s )
s 2 ps 9
The sensitivity of T to changes in p is determined by
dT p ps
S PT
dp T s 2 ps 9
Then plot the relationship
37
15,900
A system has a plant G ( s ) and a negative unity
E12.5 s s
s 1 1
100 200
feedback with PD compensato r GC ( s ) K1 K 2 s. Design GC ( s ) so that
the overshoot to a step is less than 20% and the settling time is less than 60 ms.
K
15900 K 2 s 1 2 10 4
K2
The open - loop transfer function is GGc
s s 100 s 200
K
K s 1
K2
where K 3.18 108 K 2 . Select K1 / K 2 100
s s 100 s 200
K
GGc ( s)
s s 200
K
The closed - loop transfer function is T ( s )
s 2 200 s k
Let 0.5 for P.O. 20%. 2n 200; n 200 and K n2 40000
4 4
The settling time is TS 40 ms; The controller Gc ( s ) 0.00012 s 100
n 100
38
DP12.10
1
G(s) (The space robot transfer function)
s ( s 10)
GGc ( s ) K
(a) Consider Gc ( s ) K ; T ( s )
1 GGc ( s ) s 2 10s K
For PO 4.5%; 0.702; K 50.73; Gc ( s ) 50.73
(b) Consider PD controller : Gc ( s ) K1 K 2 s
K1 K 2 s
T ( s) ; Use ITAE method : K1 100; K 2 4
2
s (10 K 2 ) s K1
100
Gc ( s ) 4 s 100; G p ( s )
4 s 100
K 2 K1s K 2
(c) Consider the PI controller : Gc ( s ) K1
s s
K1s K 2
T ( s) ; Use ITAE : n 5.7; K1 70.2; K 2 186.6
3 2
s 10s K1s K 2
186.6
Gc ( s ) 70.2 186.6 / s; G p ( s ) 39
70.2s 186.6
K1s 2 K 2 s K 3
(d) Consider PID controller: Gc ( s )
s
K1s 2 K 2 s K 3
T (s) ; Use ITAE with n 10
3 2 2
s 10s K1s K 2 s K 3
7.5s 2 215s 1000
K1 7.5; K 2 215; K 3 1000; Gc ( s )
s
1000
Gp ( s )
7.5s 2 215s 1000
Find a summary of the performance for the four cases : K, PD, PI, and PID
Performance means overshoot,settling time, and peak time
l 1 2
PO 100e (Eq 5.15); t p (Eq. 5.14)
n 1 2
40