Introduction To Control System Design
Introduction To Control System Design
DESIGN
Performance specifications
3
Transient response
When,
6
Unit ramp response of the first order systems
When,
9
Relation between the characteristic equation roots of the
second order system and the damping ratio and the nature
frequency.
Overdamped
Critically damped
10
Underdamped
Undamped
11
Divergent
12
13
Delay time
Rise time
Peak time
Maximum
Overshoot
Settling time
14
Zn
2
G(s)
C (s) R( s) 2 R( s)
1 G(s) s 2]Z n s Zn
2
e ]Zn t 1 ] 2
c(t ) 1 sin(Z d t tan 1 )
1 ] 2 ]
15
Rise time:
Peak time:
16
c(t p ) c(f)
Mp e (V / Z ) S
d
e (] / 1] 2 ) S
c (f )
Settling time: 4
ts (2%criterion)
]Z n
3
ts (5%criterion)
]Z n
17
Example:
18
Design consideration related to transient
response (High Order System)
19
20
a q aj r
b ( s ] k Zk ) ck Zk 1 ] k2
C (s) k
s j 1 s pj k 1 s 2 2] k Zk s Zk2
q r
a a je bk e] k Zk t cos Zk 1 ] k2 t
p jt
c(t )
j 1 k 1
r
ck e ] k Zk t
sin Zk 1 ] k2 t
k 1
21
3. If a zero is located very near the
dominant complex poles, the zero will
affect the transient response of the system.
22
23
24
Typical steady-state error due to a ramp-function input.
25
26
27
The steady state error in unity feedback system
is defined as E(s).
sR ( s )
ess (t ) lim e(t ) lim sE ( s ) lim
t of s o0 s o0
1 G ( s)
28
29
sR ( s ) s 1 1
ess (t ) lim lim
s o0
1 G ( s ) so0 1 G ( s ) s 1 G (0)
Kp lim G ( s ) G (0)
s o0
1
ess (t )
1 Kp
30
Static velocity error constant.
31
sR ( s) s 1 1
ess (t ) lim lim lim
s o0 1 G ( s ) s o0 1 G ( s ) s3 s o0 s 2G ( s )
Ka lim s 2G ( s)
s o0
1
ess (t )
Ka
32
33
K (Ta s 1)(Tb s 1)(Tm s 1)
G(s)
s N (T1s 1)(T2 s 1)(Tp s 1)
34
35
36
2. It is generally desirable to increase
the error constants, while
maintaining transient response
within an acceptable range.
37
38
39
2. The value of the loop gain effects the
locations of the closed loop poles (roots
of the characteristic equation).
3. The Characteristic equation: 1 KG ( s ) H ( s ) 0
C (s) KG ( s )
R( s) 1 KG ( s ) H ( s )
40
41
KG ( s ) H ( s ) 1
42
Procedures
43
44
45
Example:
K ( s 1.5)
G ( s)
s( s 1)( s 10)
46
Design methods
Compensators
Lead compensators
Lag compensators
PID compensators
47
Compensation configuration
Series compensation
Parallel compensation
State feedback
48
49
Series compensator
2
Effects of the addition of poles
The addition of a pole to the open loop
transfer function has the effect of pulling the
root locus to the right, tending lower the
systems relative stability and to slow down
the settling of the response.
PID Controllers
KI
Gc ( s ) Kp KDs
s
PD Control
PI Control
Gc ( s ) K p KI / s
5
Example: Consider the following system:
Motor
Vi
1
Z
G(s)
+
- s 2 5s 6
Figure 1
8
Step response for P and PD control:
10
11
Step response for P, I and PI control:
Blue: original
Green: I control
Red: PI control
12
13
numc=[1];
denc=[1 0];
sysc=tf(numc,denc);
sys2=series(sysc,sysg);
figure(2);
rlocus(sys2);axis([-4 0 -6 6]);
grid;
hold on
plot([0 -6],[0 -6*tan(acos(0.5))],'--',[0 -6],[0 6*tan(acos(0.5))],'--');
hold off
[ki,ploes]=rlocfind(sys2)
14
numc=[1 1];
denc=[1 0];
sysc=tf(numc,denc);
sys3=series(sysc,sysg);
figure(3);
rlocus(sys3);axis([-4 0 -6 6]);
grid;
hold on
plot([0 -6],[0 -6*tan(acos(0.5))],'--',[0 -6],[0 6*tan(acos(0.5))],'--');
hold off
[kpi,ploes]=rlocfind(sys3)
15
figure(4)
sys1_o=kp*sys1;sys1_cl=feedback(sys1_o,[1]);
sys2_o=ki*sys2;sys2_cl=feedback(sys2_o,[1]);
sys3_o=kpi*sys3;sys3_cl=feedback(sys3_o,[1]);
[y1,t]=step(sys1_cl,t);
[y2,t]=step(sys2_cl,t);
[y3,t]=step(sys3_cl,t);
plot(t,y1,'-',t,y2,'--',t,y3,':');grid;
xlabel('time [sec]'),ylabel('y(t)')
title('Gc(s):proportional (solid), integral (dashed) & PI (dotted)')
16
Lead compensation
17
18
19
20
Pole-zero configurations
21
22
23
2. By drawing the root locus plot of the
uncompensated system ascertain
whether or not the gain adjustment alone
to yield the desired closed loop poles. If
not, calculate the angle deficiency I ,
which must be contributed by the lead
compensator, if the new root locus is to
pass through the desired locations for the
dominant closed loop poles.
24
25
4
s2 +2s
Step Transfer Fcn Scope
s2
Answer:
4
s4
26
Lag compensation
6
The procedures of designing a lag
compensator are:
1
s
Gc ( s ) K c T 1 E
1
s
ET
9
Example Consider the system shown below.
10
11
12
13
14
15
16
17
Parallel compensation
10 1
s+2 s
Step Transfer Fcn Transfer Fcn1 Scope
Gain
(K=0.18)
18
Example: Design the value of K such that the
damping ratio of the dominant closed loop pole is
0.4.
20 1
s2 +5s+4 s+1
Step Transfer Fcn Transfer Fcn1 Scope
Gain
(k=0.449)
19
Bode Diagram
Relative stability
J 180 0 I
Gain Margin
The phase crossover frequency is defined to be the
frequency at which the phase angle of the open loop
transfer function equals 180 degree.
The gain margin is the reciprocal of the magnitude
|G(jw)| at the phase crossover frequency.
1
Kg
G ( jZ1 )
Z1 is the phase crossover frequency.
In terms of decibels,
1
Zr Z n 1 2] 2 Mr
2] 1 ] 2
The magnitude of the resonant peak gives an
indication of the relative stability of the
system.
A large resonant peak magnitude indicates a
small damping ratio.
Phase margin and the damping ratio are directly related.
J
]
100
C ( jZ ) G ( jZ )
Me jD
R ( jZ ) 1 G ( jZ )
K
G(s) K 1 :
s ( s 1)(0.5s 1)
CONTROL SYSTEM DESIGN BY FREQUENCY
RESPONSE
1
s
Ts 1 T s z1
Gc ( s ) K cD Kc Kc
DTs 1 s
1 s p1
DT
j ZT 1
K cD 0 D 1
jZDT 1
1D
2 1D
sin I m
1 D 1 D
2
Bode plot of a
lead
compensator
Kc 1, D 0.1
1 1 1 1
log Z m (log log ) , Z m
2 T DT DT
j ZT 1 1
jZDT 1 Z Z 1 D
DT
m
The lead compensator is basically a high
pass filter.
1
s
T Ts 1
Gc ( s ) Kc K cD
s
1 DTs 1
DT
K K cD
Lag compensation
1
s
Ts 1 T s z1
Gc ( s ) KcE Kc Kc
ETs 1 s
1 s p1
ET
jZT 1
Kc E E !1
jZET 1
1
s
T Ts 1
Gc ( s ) Kc Kc E
s
1 ETs 1
DT
K Kc E
Ts 1 Ts 1
Gc ( s )G ( s ) K G( s) G1 ( s )
ETs 1 ETs 1
G1 ( s ) KG ( s )
Determine gain K
to satisfy the requirement on the given
static error constant.
Bode plots of compensated and uncompensated system
Unit step response
Examples:
x1 vo
x
x2 vo
x 0 1 x 0
xx1 1 R 1 1 vi
x LC
L 2 LC
x
2
x1 vo
x2 i
1
x 0 0
xx1 C x1 1 v
1 R x2 i
x L
2
L L
Note:
The state variables need not be physically
measurable or observable quantities.
x1 ia
x2 Z
Ra Kb 1
x L
La x1 La
0
vt
xx1 a
x B x 1 TL
2 Km 2 0
J J J
b2 s 2 b1s b0
T ( s)
s a2 s 2 a1s a0
3
x
x1 x2
x
x2 x3
x
xx1 0 1 0 x1 0
x 0 0 1 x2 0 r
x2
x3 a0 a1 a2 x3 1
x1
y >b0 b1 b2 @ x2
x3
Deriving the state space model from the
system block diagrams
5( s 1) 1 6
R( s) o o o o Y (s)
s5 s2 s3
x
xx1 3 6 0 x1 0
x 0 2 4 x 5 r
x2 2
x3 0 0 5 x3 5
x1
y >1 0 0@ x2
x3
T ( s ) C( sI A) 1 B D
Q( s)
sI A
adj ( sI A )
( sI A ) 1
det( sI A )
x 0 1 x1 1 1 u1
xx1
x
24
4 x 2 0 1 u 2
2
y1 1 0 x1 0 0 u1
y
2
0
1 x 2 0
0 u 2
A >0 1 ; 24 4 @
B >1 1 ; 0 1@
C >1 0 ; 0 1@
D >0 0 ; 0 0@
[ num , den ] ss 2 tf ( A , B , C , D ,1 )
num
0 1 4 U1
0 0 24
dem
1 4 5
Y1 ( s ) s 4
U 1(s) s 2
4 s 25
Y2 (s) 24
U 1(s) s 2
4 s 25
t
x (t ) e At x(0) e A ( t W ) Bu (W )dW
0
Laplace transform approach to solve the time
invariant state equation
k k
1 2 2 1 f A t
e At I At A t Akt k
2! k! k 0 k!
If the system is completely state controllable, then the
following matrix has a rank of n.
>B AB x x x A n1B @
Observability
C
CA
.
.
. n 1
CA
x
x Ax Bu
y Cx
Solving pole placement problem with Matlab
Meeting steady state response