Control Principles
Control Principles
Y(s)=T(s)R(s)
Ec ( s ) = (1 T ( s )) R( s )
Step input
r (t ) = A
A
R( s) =
s
Ramp input
r (t ) = At
A
R( s) = 2
s
Parabolic input
r (t ) = At 2
2A
R( s) = 3
s
R( s) = ( s ) = 1
(t ) dt = 1
(t a ) g (t )dt = g (a )
Why?
1
, t ; ( 0, t 0)
f (t ) =
2
2
0, otherwise, f (t ) (t )
(t )dt
= lim f (t )dt
0
-/2
1/
/2
/2 1
= lim / 2 dt = 1
1
, t a ; ( 0, t a )
f (t a ) =
2
2
0,
otherwise, f (t a ) (t a )
f
a
-/2
1/
/2
t
6
1
, t a ; ( 0, t a )
f ( t a ) =
2
2
0,
otherwise , f (t a ) (t a )
-/2
(t
1/
/2
= lim
a + / 2
0 a / 2
g (t ) dt = g (a ) lim
( t
a + / 2
0 a / 2
dt = g (a)
a ) g (t ) dt = g ( a )
7
g (t ) r ( ) dt = 1{G ( s ) R ( s )}.
r ( ) ( )
y (t ) =
g (t ) ( ) dt.
= t g (t ) ( )dt = g (t ) {if = 0}
For R(s)=(s)=1
Y ( s) = G ( s) R( s) = G ( s)
n!
R( s) =
s
n +1
2
/ s( s + 2n )
G( s)
n
R( s)
Y ( s) =
R( s) =
2
1 + G( s)
1 + n / s( s + 2n )
n2
= 2
R( s)
2
s + 2n s + n
n2
1
= 2
s + 2n s + n2 s
Transfer function:
Characteristic equation: s 2 + 2n s + n2 = 0
Poles:
Im(s )
s1, 2 = n jn 1 2
n 1 2
Re(s )
cos 1
time constant
oscillation
frequency
10
n2
Transfer function:
T ( s) = 2
s + 2n s + n2
Laplace transform of step response (assuming zero states at t=0):
n2
1
Y ( s) = T ( s) =
s
s ( s 2 + 2n s + n2 )
12
y (t ) = 1
e nt sin( n t + ),
Figure 5.5 (a) Transient response of a second-order system for a step input.
(b) The transient response of a second-order system for a step input as a
13
function of and n t .
Natural Frequencies
For an under-damped system, the step response will be:
y (t ) = 1
1 nt
sin( n t + ),
e
1 2 ,
0 < <1
= cos 1 cos =
sin = 1 cos = 1 =
2
n 1 2
(time constant)
(oscillation/vibration frequency)
dy (t ) n nt
=
e
sin n t
dt
14
y (t ) =
n nt
sin n t , =
e
1 2 ,
dy / dt = n e nt sin n t
15
16
Second-Order System
y (t ) = 1
1 nt
sin( n t + ), (0 < < 1)
e
Percent Overshoot =
M p fv
fv
100
(usually corresponding
response within 2%)
17
Second-Order System
y (t ) = 1
1 nt
sin( n t + ), (0 < < 1)
e
Percent Overshoot =
M p fv
fv
100
(usually corresponding
response within 2%)
18
Setting Time
y (t ) = 1
e nt sin( n t + ),
corresponding to 2% of input
y (t ) 1
nTs
1 nt
e
=
sin(n t + ) 0.02
sin(n t + )
0,sin 0, 0
sin x
sin(n t + ) 0
1
=
lim
=
x
x
0
. { = 1 /( )}
n
n
19
e nt sin( n t + ), =
n nt
=
e
sin n t 0 t =T
p
1 2 ,
dy (t ) n n t
=
e
sin(n t + ) n e n t cos (n t + )
dt
n nt
[ sin(n t + ) cos (n t + )]
e
= cos 1 cos =
sin = 1 cos 2 = 1 2 =
=
dy (t ) n n t
=
e
[ sin(n t + )cos sin cos (n t + )]
dt
n nt
e
[ sin(n t + )]
= n e n t sin n t
20
dy (t )
dt
sin n T p = 0 n T p = n , n = 0,1,
for
=
n 1 =
Tp
n 1 2
1 nT p
Mp =
y (T p ) =
1 e
sin(n T p + ),
= 1 e
nT p
sin ( + )
M pt = 1 + e
nT p
= 1+ e
1 2
= 1+ e
nT p
sin
= cos 1 cos =
,
sin = 1 cos 2 = 1 2 =
M pt fV
P.O. =
100%
fV
P .O . = 100e
1 2
21
Second-Order System
Swiftness of response: rise time and peak time
Tp =
n 1
Tpn =
1 2
M pt = 1 + e
/ 1 2
4
=
=
T
4
. Note e - 4 0.02
With = 0.02 f v then s
n
Second-Order System
T pn =
M pt = 1 + e
/ 1 2
23
Second-Order System
0.2
24
Second-Order System
25
Second-Order System
Normalized rise time
26
=
1 2
Percent Overshoot =
M p fv
fv
100
27
n2
1
Y ( s) = T ( s) =
s
s ( s 2 + 2n s + n2 )
1
1 ent sin(n t + cos 1 ),
y (t ) =
Ts = 4 =
( for = 2% )
PO = 100e /
Tp =
Tr1
1 2
1 2
=
Settling time
Percentage
overshoot
Peak time
n 1 2
2.16 + 0.6
Im
Re
Ts = 4 =
29
n2
( s 2 + 2n s + n2 )(s + 1)
3rd system
(n = 1)
Experiments show as far as percent
overshoot and settling time are
concerned, it can be approximated
by second-order when
1 / 10 n
1/
n
30
n2
(5.15)
( s 2 + 2n s + n2 )(s + 1)
R ( s ) = 1 / s = 0.45 n = 1
| n |= 0.45
P .O . = 100e /
1 2
P .O . = 20%
, Ts =
0.45
Ts = 8.9 sec .
31
n2
(5.15)
( s 2 + 2n s + n2 )(s + 1)
R ( s ) = 1 / s = 0.45 n = 1
P .O . = 100e /
| n |= 0.45
1 2
P .O . = 20%
, Ts =
0.45
Ts = 8.9 sec .
32
Ts =
( for = 2% )
1 2
4
<4
n
n > 1
?
?
33
( for = 2% ) P .O . = 100e /
1 2
4
<4
n
n > 1
= cos 1 = cos 1 0.707 = 450
34
Second-Order System
PO = 100e
T pn =
overshoot<5%
35
= 1 / 2 and n = 1 / = 2
4 and
1
n =
2
2 =1
= n2
= 2 and p = 2n = 2
Proposed design
36
Y ( s ) Pi ( s ) i ( s )
=
T (s) =
,
R( s)
( s)
1
R( s ) =
s
N
Bk s + C k
1 M Ai
Y (s) = +
+ 2
,
2
2
s i =1s + i k =1 s + 2 k s + ( k + n )
M
y (t ) = 1 + Ai e
i =1
i t
+ Dk e k t sin(k t + k )
k =1
37
38
39
E ( s ) = R( s ) Y ( s ) = (1 G ( s )) R( s )
is the dc gain.
40
For
(1 + GH G )
R( s)
1 + G (s) H (s)
H (s) = 1
E (s) =
(1 + GH G )
1
R( s) =
R( s)
1 + G (s) H (s)
1 + G (s)
41
42
( for H ( s ) = 1)
We can calculate
ess for different
types of standard
signals
System type
Type zero
Type 1
Type 2
43
A) Step Input
R( s ) =
A
s
s( A / s )
A
=
s0 1 + Gc ( s )G ( s ) 1 + lim Gc ( s )G ( s )
ess = lim
s 0
Gc ( s )G ( s ) =
K ( s + zi )
i =1
Q
s ( s + pk )
N
For N=0
ess =
k =1
A
=
1 + Gc (0)G (0)
A
M
i =1
k =1
1 + K zi / p k
A
1+ K p
i =1
k =1
A
ess =
1+ Kp
44
R( s ) =
A
s
sR( s )
s0 1 + G ( s )G ( s )
c
s( A / s )
A
=
ess = lim
s0 1 + Gc ( s )G ( s ) 1 + lim Gc ( s )G ( s )
s 0
Gc ( s )G ( s ) =
K ( s + zi )
i =1
Q
s N ( s + pk )
s 0
i =1
k =1
= lim
1 + K zi / ( s p k )
N
k =1
As N
s 0 N
i =1
k =1
=0
s + K zi / ( p k )
45
R( s) =
A
s2
s( A / s 2 )
A
A
= lim
= lim
ess = lim
s0 1 + Gc ( s )G ( s )
s0 s + sGc ( s )G ( s )
s0 sGc ( s )G ( s )
ess = lim
s 0
A
M
sK ( s + zi ) / [ s ( s + pk )]
N
i =1
= lim
s 0
k =1
A
M
K ( s + zi ) / [ s
i =1
N 1
( s + pk )]
k =1
s 0
A
M
i =1
k =1
sK zi / ( s 0 pk )
46
s 0
A
M
sK ( s + zi ) / [ s ( s + pk )]
N
i =1
= lim
s 0
k =1
A
M
K ( s + zi ) / [ s
N 1
i =1
( s + pk )]
k =1
s 0
A
M
i =1
k =1
sK ( s + zi ) / [ s ( s + pk )]
i =1
k =1
K zi / pk
i =1
k =1
K v = lim sGc ( s )G ( s ) = K zi / pk
s 0
For N2 system,
ess = lim
s 0
A
M
sK ( s + zi ) / [ s ( s + pk )]
i =1
k =1
= lim
s 0
A
Kv
velocity constant
As N 1
i =1
k =1
K zi / pk
=0
47
s( A / s 3 )
A
ess = lim
= lim 2
s0 1 + G ( s )G ( s ) s0 s G ( s )G ( s )
c
c
A
(r (t )=At 2 / 2)
3
s
s( A / s 3 )
ess = lim
= lim
s 0
1+
s 0
K ( s + zi )
i =1
Q
s N 2 A ( s + pk )
k =1
k =1
i =1
s N ( s + pk ) + K ( s + zi )
s N ( s + pk )
k =1
For type 0 and type 1 systems, the steady state error is infinite.
For type 2 systems the steady state error is given by
ess =
A
M
i =1
k =1
K zi / pk
A
Ka
acceleration error
constant
M
i =1
k =1
K a = lim s Gc ( s )G ( s ) = K zi / pk
2
s 0
48
Input
step
ess =
A
1+ Kp
ess = 0
ess = 0
ramp
parabola
A
Kv
A
Ka
Conclusions:
1.Adding integrators (increasing system type) eliminates steady state error.
2. If steady state error is finite and not zero, then increasing the system dc
gain (increasing controller proportional gain), reduces steady state error.
49
s 0
K a = lim s 2Gc ( s )G ( s )
s 0
Then for
Step input
Ramp input
Acceleration input
ess =
A
1+ Kp
A
ess =
Kv
A
ess =
Ka
50
Gc ( s ) = K1 + K 2 / s
K p = lim Gc ( s )G ( s )
s 0
Gc ( s )G ( s ) =
K2 K
K1K
+
s + 1 s(s + 1)
K v = lim sGc ( s )G ( s )
s 0
ess =
A
A
=
1 + K p 1 + KK1
ess =
A
A
=
K v KK 2
51
A
A
ess =
=
K v KK 2
52
E ( s ) = R ( s ) Y ( s ) = [1 T ( s )]R ( s )
1
For a step input ess = lim s[1 T ( s )] = [1 T (0)]
s 0
s
G( s)
K ( s + 4)
T ( s) =
=
1 + G ( s ) H ( s ) ( s + 2)( s + 4) + 2 K
4K
T ( 0) =
8 + 2K
This is obtained when K=4.
Achieving e = 0 requires T(0)=1.
ss
53
T ( s) =
where
Y ( s)
Gc ( s )G ( s )
Z (s )
=
=
R( s ) 1 + H ( s )Gc ( s )G ( s ) 1 + Z (s )
Z (s ) =
Gc (s )G (s )
1 + Gc (s )G (s )(H (s ) 1)
The loop transfer function of the equivalent unity feedback system is Z(s).
The error constants for non-unity feedback system are given by:
K p = lim Z ( s ) K v = lim sZ ( s ) K a = lim s 2 Z ( s )
s 0
s 0
Noted that: H ( s ) = 1 Z ( s ) = Gc ( s )G ( s )
s 0
K p = lim Z ( s ) = lim Gc ( s )G ( s )
s 0
s 0
54
When the fourth floor button is pressed on the first floor, the
elevator rises to the fourth floor with a speed and floor level
accuracy designed for passenger comfort.
55
56
57
59
Performance Indices-Example
A suitable performance index is the integral of the square
of the error, ISE, which is defined as
T
ISE = e 2 (t ) dt
0
The upper limit T is a finite time, usually convenient to choose T as the settling
time Ts .
A step input and response for ISE is shown here:
60
ISE = e 2 (t )dt
IAE = | e(t ) | dt
ITAE = t | e(t ) | dt
0
Integrated time-multiplied
absolute error
ITSE = te 2 (t )dt
0
Integrated time-multiplied
squared error
61
Type=?
T ( s) =
1
s 2 + 2s + 1
ITAE = t | e(t ) | dt
0
IAE = | e(t ) | dt
0
64
Normalized Responses
65
Normalized Responses
66
Normalized Responses
The best
among
ISE,IAE
and ITAE
67
Design Goal:
select K,K1,
so that:
1) Overshoot less than 10% for step input, 2) minimizing SS error for a ramp
input, and 3) reducing the effect of a step disturbance.
69
G(s)
R(s)
d ( s ) =0
R(s)=0
Td(s) +
-
Y ( s)
R( s ) T
G(s)
= T ( s) =
KG ( s )
L( s )
=
1 + KG ( s ) 1 + L( s )
Y ( s)
G( s)
G( s)
= TTd ( s) =
=
= T ( s) / K
Td ( s) R ( s )=0
1 + KG( s) 1 + L( s )
1
G (s )
R (s )
Td (s ).
1 + L (s )
1 + L (s )
70
KG ( s )
R( s)
1 + KG ( s )
1
1
1
1
=K 2
.
s( s + K1 ) + K s
s + K1s + K s
0
+ 2n s + n2 =
K1 = 2 n K = n2
/ 1 2
P .O . = 100e
= 0.6 O.P % = 9.5%
Goal: 1.
71
we have
K v = lim sG ( s )
s 0
B
B
B
= lim
=
.
ess =
s 0
Kv
sKG ( s ) ( K / K1 )
72
G( s)
G( s)
Td ( s ) =
Td ( s )
1 + L( s )
1 + KG ( s )
G (s ) 1
=
eTd () lim=
sE (s ) lim s
s 0
s 0 1 + KG (s ) s
1
1
=
lim
=
s 0 s (s + K ) + K
K
1
73
Y ( s) =
K
1
.
2
s + K1s + K s
B
B
B
= lim
=
.
s
0
Kv
sKG ( s ) ( K / K1 )
G ( 0)
1
=
1 + KG (0) K
s 2 + 2 n s + n2 = s 2 + 2(0.6)n s + K = s 2 + k1s + K .
K1 = 2(0.6) n n = K
K
K
K
=
=
.
K1 1.2 K
1 .2
K
K1
K
K1
K = 25
25
=
= 4.17 ,
1.2
K =100
ess =
K1 = 1.2n = 1.2 K
B
B
B
= lim
=
.
s
Kv
sKG ( s ) ( K / K1 )
B
B
ess |k = 25 =
=
= 0.24 B .
K / K1 4.17
B
B
100
=
=
= 0.12 B .
e
|
=
= 8.33. ss k =100 K / K
8.33
1
1.2
75
76
77
G1 =
5000
5
=
1 = 1 / 1000 = 0.001
s + 1000 1s + 1
G1 5
Dominate poles
78
Y ( s) =
=
5K a
R( s ),
s( s + 20) + 5K a
5K a
s + 20 s + 5K a
2
R ( s ),
n2
s + 2 n s
2
+ n2
R ( s ),
n2 = 5K a , and 2n = 20.
or
79
5K a
R( s ) ,
s 2 + 20 s + 5K a
Kv
The best compromise (Ka=40) still does not meet all the specifications. (Ts=400ms)
80
What to do?
Ts = 4 =
( for = 2% )
PO = 100e /
Tp =
n 1 2
Tr1
2.16 + 0.6
(0 <
1 2
=
< 1)
Settling time
Percentage
overshoot
Peak time
Rise time (10%
- 90%)
81
5K a
R( s ) ,
2
s + 20 s + 5K a
Kv
Using PD controller, or using velocity feedback to increase 2n
82
Summary
In this chapter, we have concerned with the definition and
usefulness of quantitative measures of the performance of
feedback systems.
Transient response: (for a unit step input signal:) in terms:
overshoot, peak time, and settling time
Steady state errors for different standard test signals.
Optimal control indices.
83
Y(s)=T(s)R(s)
Ec ( s ) = (1 T ( s )) R( s )
Step input
r (t ) = A
A
R( s) =
s
Ramp input
r (t ) = At
A
R( s) = 2
s
Parabolic input
r (t ) = At 2 / 2
A
R( s) = 3
s
84
Second-Order System
y (t ) = 1
1 nt
sin( n t + ), (0 < < 1)
e
Percent Overshoot =
M p fv
fv
100
(usually corresponding
response within 2%)
85
Second-Order System
y (t ) = 1
1 nt
sin( n t + ), (0 < < 1)
e
Percent Overshoot =
M p fv
fv
100
(usually corresponding
response within 2%)
86
( for H ( s ) = 1)
The general form of the open loop transfer function of a system is given
by
M
Here N is known as the
K ( s + zi )
system type, and corresponds
i =1
Gc ( s )G ( s ) =
Q
N
to the number of integrators in
s ( s + pk )
k =1
the system.
N
We can calculate
ess for different
types of standard
signals
System type
Type zero
Type 1
Type 2
87
K v = lim sGc ( s )G ( s )
s 0
K a = lim s 2Gc ( s )G ( s )
s 0
Then for
ess =
Step input
Ramp input
Acceleration input
A
1+ Kp
A
ess =
Kv
ess =
A
Ka
88
Input
step
ess =
A
1+ Kp
ess = 0
ess = 0
ramp
parabola
A
Kv
A
Ka
Conclusions:
1.Adding integrators (increasing system type) eliminates steady state error.
2. If steady state error is finite and not zero, then increasing the system dc
gain (increasing controller proportional gain), reduces steady state error.
89
T ( s) =
Y ( s)
Gc ( s )G ( s )
=
R( s ) 1 + H ( s )Gc ( s )G ( s )
T ( s) =
Y ( s)
Z ( s)
=
R( s) 1 + Z ( s)
where
Z ( s) =
Gc ( s )G ( s )
1 + Gc ( s )G ( s )( H ( s ) - 1)
The loop transfer function of the equivalent unity feedback system is Z(s).
The error constants for non-unity feedback system are given by:
K p = lim Z ( s ) K v = lim sZ ( s ) K a = lim s 2 Z ( s )
s 0
Noted that:
s 0
H ( s ) = 1 Z ( s ) = Gc ( s )G ( s )
s 0
K p = lim Z ( s ) = lim Gc ( s )G ( s )
s 0
s 0
90
ISE = e 2 (t )dt
IAE = | e(t ) | dt
ITAE = t | e(t ) | dt
0
Integrated time-multiplied
absolute error
ITSE = te 2 (t )dt
0
Integrated time-multiplied
squared error
91