De Silva Clarence W Sensors and Actuators Control Systems Instrumentation
De Silva Clarence W Sensors and Actuators Control Systems Instrumentation
De Silva Clarence W Sensors and Actuators Control Systems Instrumentation
Design of
MODERN
CONTROL
SYSTEMS
Edited by
D J Bell PA Cook N Munro
l
l
CONTROL
SYSTEMS
Other volumes in this series
Volume 1 Multivariable control theory J. M. Layton
Volume 2 Elevator traffic analysis, design and control G. C. Barney and
S. M. Dos Santos
Volume 3 Transducers in digital systems G. A. Woolvet
Volume 4 Supervisory remote control systems R. E. Young
Volume 5 Structure of interconnected systems H. Nicholson
Volume 6 Power system control M. J. H. Sterling
Volume 7 Feedback and multivariable systems D. H. Owens
Volume 8 A history of control engineering, 1800-1930 S. Bennett
Volume 9 Modern approaches to control system design N. Munro (Editor)
Volume 10 Control of time delay systems J. E. Marshall
Volume 11 Biological systems, modelling and control D. A. Linkens
Volume 12 Modelling of dynamical systems—1 H. Nicholson (Editor)
Volume 13 Modelling of dynamical systems—2 H. Nicholson (Editor)
Volume 14 Optimal relay and saturating control system synthesis E. P. Ryan
Volume 15 Self-tuning and adaptive control: theory and application C. J. Harris
and S. A. Billings (Editors)
Volume 16 Systems modelling and optimisation P. Nash
Volume 17 Control in hazardous environments R. E. Young
Volume 18 Applied control theory J. R. Leigh
Volume 19 Stepping motors: a guide to modern theory and practice P. P. Acarnley
Volume 20 Design of modern control systems D. J. Bell, P. A. Cook
and N. Munro (Editors)
Volume 21 Computer control of industrial processes S. Bennett and
D. A. Linkens (Editors)
Volume 22 Digital signal processing N. B. Jones (Editor)
Volume 23 Robotic technology A. Pugh (Editor)
Volume 24 Real-time computer control S. Bennett and D. A. Linkens (Editors)
Volume 25 Nonlinear system design S. A. Billings, J. O. Gray and
D. H. Owens (Editors)
Volume 26 Measurement and instrumentation for control M. G. Mylroi and
G. Calvert (Editors)
Volume 27 Process dynamics estimation and control A. Johnson
Volume 28 Robots and automated manufacture J. Billingsley (Editor)
Volume 29 Industrial digital control systems K. Warwick and D. Rees (Editors)
Volume 30 Electromagnetic suspension—dynamics and control P. K. Sinha
Volume 31 Modelling and control of fermentation processes J. R. Leigh (Editor)
Volume 32 Multivariable control for industrial applications J. O'Reilly (Editor)
Volume 33 Temperature measurement and control J. R. Leigh
Volume 34 Singular perturbation methodology in control systems D. S. Naidu
Volume 35 Implementation of self-tuning controllers K. Warwick (Editor)
Volume 36 Robot control: theory and applications K. Warwick and A. Pugh
(Editors)
CONTROL
SYSTEMS
Reprinted 1988
While the author and the publishers believe that the information and guidance
given in this work are correct, all parties must rely upon their own skill and judgment
when making use of them. Neither the author nor the publishers assume any liability
to anyone for any loss or damage caused by any error or omission in
the work, whether such error or omission is the result of negligence or any other
cause. Any and all such liability is disclaimed.
Bell, D. J.
Design of modern control systems.
(IEE control engineering series; v. 20)
1. Control theory
I. Title II. Cook, P. A. III. Munro, N.
ISBN 0 906048 74 5
Preface
1 State-space theory 1
1.1 Introduction 1
1.2 Controlability and observability 2
1.3 Standard forms
1.4 Calculation of transfer function
matrix 9
1.5 Realizations 11
1.6 Inverse Systems 14
References 15
Problems 17
D.J. Bell
P.A. Cook
N. Munro
Manchester, 1982.
Chapter 1
State-space theory
Dr. P. A. COOK
Synopsis
1.1 Introduction
x = Ax + Bu (1.1)
y = Cx + Du (1.2)
x(0) = xQ (1.3)
y(t) =
- * r.
C exp(At)x n + I C exp {A (t—r) }Bu (T )dx + Du(t)
(1.4)
R = W = [w1,...,wn]
so that
w"'1AW = diag (A1,...,An) (1.6)
x(t) = 0
exp(-AT)Bu(i)dT
0
exp(-AT)B
rank $ = n (1.10)
where
$c = [B / AB / ...,A n " 1 B] (1.11)
C exp (At)
rank $ = n (1.12)
o
where
$ = [C ,A C ,...,(A ) n ~ C ] (1.13)
X X X X
0 A 0 X B
CO c
0 0 B =
A
uu X 0
0 0 0 A J 0
uo
[0 x ]
rank B = m (1.14)
H = b ] (1.15)
m
where B = [b,
l ,...,b m ] and the integers (y,i , . . .mf ii_.) are
known as the controllability indices of the system. Because
of (1.14), they are all strictly positive and, since H has
n columns,
rank [B,AB,...,AV~1B] = n
A -* H~ 1 AH, B + H" 1 B
"•• A lm
H~ 1 AH
A
ml mm J
A =
3i
v \ 0
0. . .0 '1 a 0 ... 0 a .
l
(1.16)
= v = v
c ' y r+l c ~
n = m(v -1) + r
T T —1
Q T A(Q T )
T
Q B
rank C = p
max
as follows: since
R Z
n-1
-1 =
D — ID 4- a T
R + a
n-2 " n-1 n-lx
! (1.18)
R = AR + a
0 1 lI
0 = AR Q + a Q I
tr R o = tr
(1.19)
tr R Q = tr (ARX) + n a x
O = tr (ARQ) + n a Q
we get
(n-1) an—l, =, tr Rn—z
o
i
a tr R
l = 0
0 0 -a r 0. , ...
or
0 0 1 -a -a r -a_
from which det (sI-A) and adj (sI-A) may readily be ob-
tained. The transformation can be performed with O(n )
multiplications and divisions, but it is also subject to
numerical difficulties and there does not seem to be any
convenient way of checking the results.
1.5 Realizations
n
i
d
s io
n ± -l
3
xo
B
l
A = , B = , c = [<:,
0 ^ A 0
where
0 1 0 0
X
I V
B
i =
-d.
iO
and d
i0 J
G(z) = D + -i
K.z
K
2 ••' K
n .. . K
K K =
. . Ko
n 2n-l -1
_. T _ ~
(1.21)
K = 0 $ , K =
rank K = n
MKN = (1.22)
0
10
by using, for example, Gauss reduction . From (1.21) and
(1.22), it follows that, for some nonsingular matrix L,
N = x ]
LB x CL" 1 x
MK = KN =
and
LAL" 1 x
MKN =
F(s)G(s) = s" q I
G(s)F(s) = s" g I
sI-A B
P(s) =
-C D
l
= [0 l m ]P (s) (1.25)
14
and, by suitable elementary transformations , P(s) can be
reduced to the form
References
Problems
2 1 1 0 1 0
1 0 0 -1 0 1
A = B =
1 0 -1 -1 0 1
0 1 1 -1 1 0
0 -1 -1 0
C =
1 0 0 1
Synopsis
2.1 Introduction
Nyquist-Bode approach); or
(ii) study closed-loop frequency as a function of gain (the
Evans Root-Locus approach).
£ : x = Ax + Bu
(2.1)
y = Cx + Du
linear m inputs
dynamical and
system m outputs
signal
7TFT
signal
vector vector
injection return
Fig 2.1
Fig 2.2
Fig 2.3
S(g) = A
i—
G(5)
Fig. 2.4
I
g
F(s) = I m + g""1G(s)
det[sl n - S(g)] = 0
det[gl m - G(s)] = 0
G(s) —
Fig 2.5 "j
H(s) —
E'. portion of
s-plane Voot locus
g-plane
D1
D1, portion of
nyquist locus
E, traced from
g = oo to g =
F i g . 2.6
2.3 Zeros
In physical terms zeros are associated with the vanishing
of vector gain; that is with the existence of non-zero in-
put exponential signal vectors which result in zero output .
The discussion of zeros in the literature is confusing be-
cause one can talk about zeros in connection with various
objects: £(A,B,C,D), or G(s), or g(s) for example.
The object £(A,B,C,D) has associated with it a larger set
of zeros than the object G(s) because in passing from the
representation £(A,B,C,D) to the representation G(s) one
discards phenomena (uncontrollability and unobservability)
which can be discussed in terms of zeros (the so-called de-
coupling zeros 1 ) . In turn the object G(s) has a larger
set of zeros than the object g(s) because vanishing of all
characteristic gains does not necessarily imply the vanish-
ing of the vector gain for G(s). For the purposes of this
paper we will take the system representations being discussed
to be completely controllable and completely observable, and
the zeros to be what are usually called the transmission
zeros . A useful physical interpretation of a zero then
comes from the following result.
u(t) = 1 (t)e Zt u z
and an output
zl - A -B
= 0 (2.6)
- C -D
Zero Pencil; From (2.6) we have that, for the case when
D = 0,
= Bu (2.7)
Cx z = 0
= Mv z (2.8)
ation matrix fN
B+
(z NM-NAM)v = 0
uz = B + (zI-A)Mv z
v *
A = S 2 A S1
det[S*S 1 - S*AS1] = 0
and one may say that: the finite zeros of £(A,B,C) are
the spectrum of the restriction of A to Ker C in domain
and Ker B in range.
D = Mo]
rp rp rp
where [U M ] and [V N ] are the eigenvector and dual
o o
eigenvector matrices for D and J is the Jordan block
associated with the non-zero eigenvalues of D. Further,
suppose that there is a complete set of eigenvectors spanning
the null space of D. Define A ,B ,C by
J
^ o o' o
= A -
31
BQ = BM
C N
o = oC
In case D has full rank one has that there are n finite
zeros of £(A,B,C,D) which are given by a(A-BD" 1 C), the
spectrum of the matrix (A-BD~ C ) .
S = — —E-
cp - a
28
on the complex frequency it has been shown by Edmunds
that it is possible to find a correspondingly transformed
system £(A,B,C,D) such that
G(p) = G(s)
l
A = (cA + dl) (aA + bl)
B = (cA + dI)""1B
D = D - cC(cA + dl)"1B
a g + b dq + b
Q = 2 Q = — —
^ cg + d ^ c q - a
S'(q) = S(g)
A1 = A - cB(cD + dI)~1C
C1 = ( c D + dI)"1C
G G9
G(s) = -^ + -=- +
S c,
G = CB G = CAk
l k ~lB k
= 2'3'
V
and take G. = G, .
v
Let d± be the nullity (rank defect) of G±. Then the
high-gain asymptotic behaviour of the various branches of
the generalized root-locus diagram can be determined as foll-
ows.
V
where
± 1 z a[G ]\{0} and g -> 0
v
The collection of symbols co[G.]\{0}
is to be read as
1
v
"belongs
2nd Order to the non-zeros 22 spectrum
branches: of^-tl G^±".
s o = = 1,2 , . . . (d 9 -d, )
v i,2 g z i
where A. 0 E Q [ G ] \ { 0 } and g -• 0
i/2 2
34
s
i,k
V
where A. , ea[G,]\{0} and g + 0
1 ,K K
A .
i
kth Order branches: g. , - -
1,K
v
where A. , ea[G,]\{0} as before, and
1 / .K. K.
= Im B © Ker C
(Sl n _ m - NAK)
and
ai(CB)
g i /.I " JOJ
1 m - d m - d.
i
2 d d 2(d1 - d 2 )
i - 2
3 d d 3(d2 - d 3 )
2 " 3
V d vd
v-l v-l
Total v-1
Number of m m + I d.
Asymptotes i=l
v-1
Number of finite zeros = n - m - £ d.
x Ax + Bu
y Cx
u = - - Ky
S K (g) = A - | BKC
[SR(g) - s i l]x i = 0
i = 1,2,...,n
y*[SK(g) - s.,1] = 0
A-s ± I x.
-1
-C -gK w.
i = 1,2,
A-s ± I
= 0
-C -gK -1
-y.BKC x.
arg(ds ± ) = arg - - 0 < g <_ °° (2.10)
i = 1,2, ... ,p
or * -l
n.K w
arg(ds ± ) = arg — ~ 0 £ g < «> (2.11)
y x
i i
i = 1,2,. ..,p
The angles of departure from poles are found using (2.10)
with g = «, and the angles of arrival at zeros are found
using (2.11) with g = 0. Note that for g > 0 we will
have p = n and for g = 0 we will have 0 <_ p <_ n - m,
where p is the number of finite zeros and m is the number
of inputs and outputs of the vector feedback loop.
(i) All the finite zeros of the optimal feedback loop lie
in the left-half complex plane,
(ii) All the generalized Nyquist diagrams for the optimal
feedback loop have infinite gain margin at least 60°
phase margin,
(iii) When p in the performance index (2.12) is °° the n
branches of the optimal root locus (i.e. the locus of
38
2.9 Conclusions
References
Synopsis
3.1 Introduction
plant
7 8
of Doyle and Stein , Postlethwaite, Edmunds and MacFarlane
introduced a multivariable phase margin which in conjunction
with Doyle and Stein's gain margin, allows a larger class of
perturbations to be identified for which the feedback system
remains stable. These results are discussed in the latter
half of this lecture.
A system which has a large degree of stability in terms of
these multivariable gain and phase margins will have charac-
teristic gain loci which although not necessarily insensitive
to small perturbations will nevertheless be insensitive
enough for the number of encirclements of the critical point
to remain unchanged. In the first half of this lecture a
direct study is made of the sensitivity of the characteristic
gain loci to small perturbations. Sensitivity indices are
defined in terms of the left-hand and right-hand eigenvectors
of the return-ratio matrix corresponding to the point at
which the system is perturbed. It follows that the charac-
teristic gain loci are at their least sensitive when the
corresponding return-ratio matrix is normal.
Theorem 1 (Gerschgorin)
s
i = Zi 2Si i = 1,2,...,m
where T
^. is y. transposed, and y.,x. are (respectively)
left-hand and right-hand eigenvectors of T, normalized so
that the sum of the squares of the moduli of their elements
is equal to unity.
3.2.2 Analysis
W""1TW = diag(A ± )
taken to be Z^/Sj_ •
With this definition for W it follows that
e 11/s 1 e 12/s 1 e lm /s l
e 21 /s 2 e 22 /s 2 ... e 2m / s 2
W -1 (T+eE)W = diag(X i ) + e
eml /sm em2 /sm emnr m
T =
1
1
[l+(t ll"t2 t
22
<
t
l l "t22 0
=
V
1 > 1
€E(ju
K(jou) (joa)
-I
(jui)
3.3.1 Normality
K (s) I • AG ( s ) G(s)
R Q (s) A
x
R ± (s) = G(s) R Q (s)G(s)
R = UH D (3.1)
R = HLU (3.2)
and
c
l
and
C
2 (w
(a)) 4 tan
(i) A G ( S ) is stable,
(ii) {Q^(w)+e•(w) : i,j = 1,...,m} have a
spread of less than TT , for all w,
(iii) [c 1 (u))-1]C2 (w) < 1/ for all u>,
(iv) £;L (a)) + 6 1 (w) ~^ m (w) > -TT, for all w, and
e +G +l a) < 7Ti for
(v) m^) m^^ ^m^ '
1.0
0.01
0.01
3.6 Conclusion
REFERENCES
PROBLEMS
-47s+2 56s
G(s)
" (s+1)(s+2)
-42s 50s+2
0
K(s) =
0 k.
where nominally k, = k~ = 1.
0.10 10.0
0.31
0.10
0.03
-50-
-100-
-150
0.10 1.0 10.0
Synopsis
4-i\ Introduction
I main feed
mixing quench
chamber stream
x = Ax + Bu
y = Cx + Du
distu rbance^
C0
TI
'out Tm
TF feed L—- bed I r l mixing Cm bed II
T C
heater o X1 out chamber Q X"
input
III
T, T2 T.
n II 1
T 5 T 6 T7
," 'a,
control
inputs
Fig. 4.2 Reactor system components and input-output
variables. Notation: T = temperature,
C = concentration, Q = flow rate, X = state
vector.
Feed heater
-1 T + -1 T _
T o T F
Bed 1 + B1
+ D1
3
I
out
I
out
Mixing chamber
Q
= M I
out
I
"out
Bed 11
II
69
T
5
T
6
11
T
7
T
out
C
out '
-2-
Fig. 4.3 Myquist array for the second bed of the reactor.
Outputs ( T o u t , c o u t ) '' Inputs (Q, T Q ) .
71
-1.16 1.04
(4.1)
-3.35 2.34
0.896 -0.571
(4.2)
0.382 0.785
0.4r
-0.6L
-20
K 3 (s) = 4 s+O* 1
75
1.5
1.0
0.5
-0.5
-1.0
bed II
compensation •out
K3(s)K2(s)K1(s)Kh TQ ^
2-bed -out
reactor
single-loop
PI controller
s+0.2
K 4 (s) = 0.07
s+0.2
K 3 (s) = 4
O.Ar
072
0.47-
0.06
0.30-
016
/ uncontrolled T ou(
i.o-
/ controlled T ou j
time
^Controlled Cout
\
-1.0- V uncontrolled C,
6. Or
-2.0
-2.00
Bed 1
2.5936 0.0179
1.5082 0.0645 0.1760 0.0000
-0.1333 0.0946 -0.1088 0.0000
1
B = -0.5473 0.0902 D = 0.0450 0.0000
0.2840 0.1007 0.0027 0.0000
0.6519 0.1258 0.0000 0.7789
0.6540 0.1356
Mixing Chamber
Bed 11
References
Synopsis
5.1 Introduction
diag {f } (5.1)
e(s) u(s)
r(s)- j———- Ms) G(s) y(s)
F(s)
Fig. 5.1
= diag i = 1, . • . ,m (5.3)
Fig. 5.2
1 2
s+1 i+3
Q.(s) = 1 (5.6)
1
s+1 s+1
20
10 unstable
asymptotically
stable
10 20
Fig. 5.3
(s+1)(s+3)
M(s) = (5.7)
0
s-1 5s+l
(s+1)2 (s+1) 2
Q(s) = (5.8)
-1 s-1
(s+1)' (s+1)'
M(s) = (5.9)
s+2
s+1
10
asymptotically
stable
-A -2 8 10 f1
Fig. 5.4
5.3 Stability
F(s)
Fig. 5.5
T G (s) O G (s)
P G (s) = (5.10)
-V G (s)
89
T
K U
K 0 0 0 0 0 o • 0 0
V W 0 0 0 0 0 ' 0 -e 0
K K
0 0 T
G U
G
0 0 0 0 ' 0 iG 0
0 0 -vG W
G 0 -Im 0 0 f
0 -u 0
0 0 0 0 T
L U
L 0 o • 0 = 0 (5.11)
0 0 0 0 -vL W
L
0
~Im '
0 -y 0
0 0 0 0 0 0 T U 0 0
F F '
0 0 0 0 0 ~VF W
F f
-Jm
-z 0
0 0 0 -v -z
T
K
(5.12)
-VT. w,H J
The closed-loop system poles are the zeros of |T H (s)|, and
Rosenbrock has shown that
= |Q(s)|/|I+QF| (5.16)
Equation (5.5),
H.(s) = (I+Q(s)F)"1Q(s)
ni
|q ± i (s)| > I |q i ± (s)| (5.18b)
(5.19)
m
d i (s) = I |q-(s) (5.20)
I
\
L
-@-
\
\
)
\ /
V if
\ /
\ d
\
Fig. 5.6
NQ = I N± (5.21)
m A m A
NR - N = - pQ + pc (5.24)
Hence, Nn - Nu = p^ (5.25)
y n o
replaces the relationship used with direct polar plots, given
by equation (5.15).
the q..(s)
do not represent anything directly measurable
11
3
on the system. However, using a theorem due to Ostrowski
h^(s) = h ^ ( s ) - f± (5.26)
d. (s)
<i>. (s) = max ^ (5.27)
3 |f i +q i i (s)|
3 3D
-i-
Fig. 5.7
96
R = R + aR
2 2 l (5.29)
Fig. 5.8
(5.30)
(5.31)
m
(5.32)
I
k=l
m
(5.33)
ji
m
k=l
(5.34)
m m 2 fm 2
l I k .1 .a l- J ic..31
li=l -*1 1
1
k=l i=l -* i=l -
(5.35)
—Y> — =
I
k^3
2
fm
'i-l
^ji a ik a
)
• ^ k
j
- 2Xk..
(5.36)
= 0 for SI =
98
Now writing
m
(ct (5.37)
J.
k=l il
ik 6 *k
and (5.38)
m m
A.k. - xk^ = 0 (5.39)
since I |q 1k <j«)
k=l
= X (5.40)
F i g . 5.9
+ +
+ +
F i g . 5.10
100
|qj:JO)| = 1 (5.41)
5.5 Conclusions
visual insight. All the results for inverse plots apply with
suitable changes to direct Nyquist plots (see Reference 1,
pp. 174-179). We can also obtain multivariable generalisa-
tions of the circle theorem for systems with nonlinear, time-
dependent, sector-bounded gains ' . Not only do the Gersh-
gorin bands give a stability criterion; they also set bounds
on the transfer function hi:L(s) seen in the ith loop as the
gains in the other loops are varied.
References
APPENDIX
(5.44)
q . (a,s) = aq ± . (s) , j/i
m
0 = N - I N± (5.46)
±
which is (5.21).
104
1
s+1
G(s) = 1
s-1
|TG| = (s-Dts+i)'
1 1
i+T s+1
G(s) =
l 1
s+2 I+Z
(i) State one reason why the inverse systems G~ (s) and
H (s) are used in Rosenbrock's inverse Nyquist array
design method,
(ii) Sketch the inverse Nyquist array for G~ (s), and
comment on the diagonal dominance of the uncompensated
system,
(iii) Determine a wholly real forward path compensator K
such that the closed-loop system H(s), with unity
feedback, is decoupled.
105
Synopsis
6.1 Introduction
adj(I+Q(s)).Q(s) (6.1)
det(I+Q(s))
form:
g
ll (s) g
l2 (s) g
13 (s) n s (s) g 1 4 (s)
g
21 (s) g
22 (s) g
23 (s) g24(s) (6,5)
g
31 (s) g
32 (s) g
33 (s) '34'
Variables ng, yf and yf are the controllable inputs to
the plant, whilst y. represents a disturbance input. The
response of the plant following a disturbance in y is
representative of its performance following a change in
station power output demand.
Ptref J
5—- Ms)
Vfv
*v$—- k 2 (s) G(s) l
WD
A pv ref
> Att Att
&—-
t
k 3 (s)
80s
ki (S, = (6.6)
k 2 (s) (6.7)
= .114 (6.8)
alone has little effect on the boiler feed flow. The over-
all system stability characteristics cannot readily be ass-
essed from an observation of this array; the off-diagonal
elements are far too large for the diagonal elements to pro-
vide a useful guide to stability and the characteristic loci
were calculated explicitly.
20
*- *
M22 H
-2 23
-15
30
q32 ^33
A3
k =0.114
1000
-100
-500
plant
compensated
921 923
IF^oo
-1000
9
31 932
-400
was that the plant input which most influenced a given output
signal should be used for its control. In the open loop fre-
quency response array of Figure 6.6, the largest element of
row 2 is 9o3' implying that drum level control is influenced
most strongly by the boiler feedpump turbine (BFPT) actuator
adjustment. Thus the most effective way of controlling drum
level is via the feedpump turbine, rather than the feedwater
regulator valves.
i.e.
det(I+Q(s)) - (l+q1]L) (l+q 22 ) (l+q 33 ) for e small (6.10)
k,(s)
l
WD
G(s)
ally high.
f 20 30 40
-20
-30
10 15
-10L
-2
-2
200
0.6
0.3
The control over the pressure drop across the feed valves
is not critical^ the main requirement being that reasonable
feed valve openings are maintained. Control merely needs to
be such that the pressure drop across all of the feed regula-
ting valves is above a specified minimum.
Such analysis has shown that the designed scheme can with-
stand all combinations of nonfunctioning controller elements
without loss of stability. Consequently, any combination of
automatic and manual control of loops will provide stable
operation and complete freedom of choice over the order of
loop closure is permitted when commissioning.
6.7 Conclusions
References
Acknowledgment
Optimal control
Dr. D. J. BELL
Synopsis
3F
(7.1)
= 1
'2 n)
3F_
F - x -*-F- dt + dx. + dG = 0
i
3x. 3x.
(7.2)
3_F_ 3F
and F - x. (i = 1,2,...,n) (7.3)
3x. ax.
must be continuous across a corner of the extremal (i.e. the
trajectory satisfying the necessary conditions).
F - x. ^-- = constant J t e (t r
ax
2
+ [(F. - x. F v )dt: + 2F v dx. dt]
U- I X . X, 1
1 1
{F 6x
k.kv i6xk 2F
x.xv 6x
i6xk
1 K IK
i k
6x
6x
Yx i
3x X
i
Writing
6x
i6xk 6x
i6xk F
x.xv 6x
i6xk
84) .
IK
6x 6x
-•
d_
= 0
dt
x.x v
1 JC
subject to the constraints
*6 •
—•-- 6 x . = 0
3X. >•
129
= f±(x,u,t)
(Uj-dj)(bj-Uj) - a* = 0
f (f f f
* = = l 2 n)T
(i = 1,2,...,n)
d =
H = H"' ^T- constant (j = l,2,...,m)
9x
3x.J i 3z.J
1 3 (7.6)
The transversality conditions are
F — x. —— dt + - 3 — dx. + - — dz.) + dG
X D
O 1
3x, 3x.
l
dz.3 It
;
o
.
= 0 (7.7)
d*q = 0
^T- = 0 (j = 1,2,...,m)
3zj
te (7 8)
~r- = w: - ° ^''i' -
dz. D
to
" A. = |-- (L + X±f±) (7.9)
= - ix- U - 1,2,...,n)
|H dt - A i dx. dG = 0 dip = 0
o
A
ivul; " 3x ± (t 1 )
j
On u s i n g e q u a t i o n s ( 7 . 1 0 ) - (7.11) , t h i s i n e q u a l i t y r e d u c e s t o
L(x,Z,t) + X.f.(x,Z,t) - L(x,z,t)
- Xifi(x,z,t) >_ 0
(7.12)
(7.13)
H
.u^ u.uh638h)dt
subject to = 0
3f
i i
+
\ Fu~
Definition
i2 = -j (n Gxx n)
tf +
f'4 • z
\- 6 ux
H n)dt
(7.14)
Jo
subject to the equations of variation
= f
- x" + (f f
x u" V '
The second variation becomes
T ( G +
2 = 1 ^ xx « «T G
xx f
u * +
\ ? f
u G xx f
u *>\
tf
+ Hui v fu )?) + (-9- C R C+C Q a+ 4 a X P a) dt
t J 0
r\2 \
—.K- H > 0 t e (t ,t~) (7.16)
2 U f
dt J * °
This inequality is a generalized Legendre-Clebsch necessary
condition in that it is similar in form to the classical
Legendre-Clebsch necessary condition
136
i~ H > 0 t e (t ,tf)
9u u - or
If (7.16) is met with equality for all t e ^o^f^ then
the functional I2 is totally singular and another transfor-
mation must be made. In this way the generalized necessary
conditions are obtained:
H t€
k £P u =° J W
p an odd integer
+ 3 T H u x n)dt
1
t n +AT
-I. s\
T
(H
UX
+f
U
Q)n dt, n (t )
X
= 0 (n continuous)
i2 = AT
0(AT 2 ))
Similarly, we let
Clearly,
S(x(N),Nh) = <f)(x(N)) (7.20)
i? 7 (IS. f i (x ' u)
3S
u must be chosen to minimize the sum f. -r— resulting in a
1 oX:
control law u(x,t). The partial differential equation
(7.21) may then be solved for S(x,t) subject to the bound-
ary condition
S(x(t 1 ),t 1 ) = 0
x = Ax + Bu
(xTQx + u T Ru)dt
o
(Q and R constant, symmetric matrices with R positive
definite) the solution of the Hamilton-Jacobi equation can
be obtained in closed form. The Hamiltonian is
H = x T Qx + u T Ru + VS T (Ax+Bu)
Assuming u is unbounded,
1-3 = 2Ru + B T VS = 0
oU
whence u° = - -• R " 1 B T VS(x,t) (7.25)
X
(VS) AX - ^(VS) BR B 1 VS + X X Qx + ~ = 0
T
It can be shown that S(x,t) = x P(t)x (P > 0) is a sol-
ution to this equation provided that the matrix P satisfies
the Matrix Biooati equation
P + Q - PBR~ 1 B T P + PA + A T P = 0 (7.26)
u°(x,t) = -R""1BTPx
= -K T (t)x say.
Q - PBR~ 1 B T P + PA + A T P = 0
References
Problems
xx = cos x 3 , x 1 (0) = 0, x x (t f ) = /2
x3 = u , x 3 (0) = 0, x 3 (t f ) free,|u|<l
Synopsis
8.1 Introduction
where
6(z,p,u)) 4 [I6(z,p,u))] - a[R6 (z,p,aj) ] 2 + 6
and
9(z,p,u>) 4 det[T(z,p,jw)]
or, equivalently:
T h
An example of such a function is x •> (x Px) where P is
positive definite. For all x in Rn let the set B(x) be
defined by:
then:
(a) x ( t , z , x Q ) e cB for all x Q e B , all t e [0,oo)f
x e Y
i i-l Y i-2 ••' ^ O B ( X O }
ty ( z ) 4 max { <>
f (z , a ) | a £ A}
and
ty ( z ) 4 max { 0 , i> ( z ) }
F = {z e R n | xp ( z ) <_ 0}
= {ze R 11|> ( z ) , = 0}
and
r
A (z) A_ {a A (z) | a is a local minimizer of
$ (z, •) on A}
(z)
F
A A {z|<j)(z,a) £ 0 for all a e A^}
12 13
8.4.4 Cut-map algorithms '
References
Synopsis
9.1 Introduction
where
n
4 = J Iq-jJ i = 1,2,...,n
where
J |q±j| i = 1,2,...,n
9.2.1 Decomposition
Q (jw) = G (joa) K
We first note that both |q.^| and <S? are positive first
order homogeneous functions of K in the sense that for any
A ± 0
(9.3)
162
= G(u r )K
163
dc.
P, : min max{J?[P,Q]}
Pi1
An equivalent problem to Pj^ is that of choosing a set of
n elements of DC which form a 'perfect matching1 of D C
and additionally minimise the maximum element in the match-
ing. A set of n elements forms a perfect matching of D C e
/? n x n if there is exactly one element of the set in every
row and column of D . A simple heuristic for finding such
a set (which usually suffices for 'by eye 1 solution of prac-
tical problems) is as follows:
Algorithm^ _(to_find P^
_ , _c n nxn
Data: D eR
Step 0: Set P 1 = D° ; l =1
Step 1: 'Colour' the smallest 'uncoloured' element of P
Step 2: As each additional element is coloured, check to
see if the coloured set contains a perfect matching.
If not, return to Step 1.
Step 3: Colour the element in D c corresponding to the last
£
element coloured in P . If £ = n, stop. Other-
0—1 ( Q—1 W (o - "\ )
wise form a new matrix p ^ - L e / ? ^ - L ^ v ) c ± ; b y
removing the row and column of P£ containing the
last coloured element. Set £ = £+1.
Return to Step 2.
y(s) = Gu(s)
Z = SG
min max{J?(Z)}
[s] i *
166
otherwise
P^ '• min(n}
s
subject to: £ m..s. <_ n s. ; i = l,2,...,n
31 3 1
s.>0,
= 1
Data: e > 0
Step 0: Set i = 1; n, = max J m.
i jl D
z£ >_ 0
£
Step 2: d = n1/2 ; if d < e then stop.
167
Step 3: If zz = 0 then n £ + 1 = n£ - d
[K]
z = Qk
l|k|| = 1
• j I 2 4 = /a?j+ B?,
p
j j = l,2,...,n and
c n
6"
: mln
k
where w., i = l,...,r represent selection frequencies
168
P-3 : min n
3
[k]
«2 ~ T)\ l Z il I 1 ° (9.4a)
Z = QK (9.4b)
||k|| = 1 (9.4c)
(Conceptual)
Data: e > 0
Step 0: Set £ = 1, n, = S^ik), k i = 1, i = l,...,n
Step 1: If inequalities (9.4) are consistent, set d = r^/2
n = n - d
= n +
otherwise set ^o+i c ^
Step 2: If d < e stop, otherwise set £ = £+1.
Return to Step 1.
6
L V = (0.0968, 527.4, 16.06)
(7 92 6
L = ' ' ' 85 ' 0-194)
(0
= '9'
where U = S
2P2S1P1 and L =
with
0 0 2.51990E-2
U = 0. 5707 0 0
- 0 0. 330706 0
The forward gain factor for stable and good responses is:
u y
U Kc L G
J
w
Fig. 9.1
i\r\ _ . . .
F i g . 9.2
F i g . 9.3
173
9.7 Discussion
KC = K + K I /S
References
PROBLEMS
1 + 3 , 2 + 1 ,3 + 5 ,4 + 2 , 5 + 4
1 2
3 4
n s |g i .|
176
n
n si
^ 1 I T~ m i-i' j = 1
i=l s j 1 D
S
2
Chapter 10
Pole assignment
Professor N. MUNRO
Synopsis
10.1 Introduction
pole assignment
state ^ output
Feedback feedback
I
canonical non-canonical
r dyadic
full rank
I
spectral mapping
constant
I
dynamic
Fig 10.1 Pole Assignment Algorithms
Ax + Bu (10.1)
y = Cx (10.2)
where A,B,C are real matrices with dimensions n*n, nxm, £*n,
respectively, or by its corresponding transfer-function matrix
representation
179
when and how can we apply feedback to this system such that
the resulting closed-loop system has a desired set of arbit-
rary poles, or eigenvalues (subject to complex pairing). First,
we consider some of the well established results, starting
with the state-feedback case.
u = r - K x (10.4)
x = (A-BK)x + Br (10.5)
K = fmt (10.6)
m1 = I 6. v* (10.7)
1 x
q
n
(10.8)
P n U.-X.)
1 3
j=l
p = <
i -i'B i > (10.9)
AQ(S) =
(10.11)
5
then we can define a difference polynomial (s) as
Ms) Ad(s) - A Q ( S )
0
1
a,
X = (10.15)
a a
n-2 n-l
a a
n-l n-2
and bf = Bf (10.18)
N 2 (s)
(10.19)
N n (s)
A,(s) = |sI-A+BK|
|sI-A+bfm
(10.21)
n
I m.N.(s) = 6(s) (10.24)
x ±
i=l
m
l
1 1
1
=
1
(10.25)
1 1
6
nn • n •
i.e. Xm = 6 (10.26)
= X-1 (10.27)
for any 6_ .
11
A = (10.30)
I A ml A
mm J
where the diagonal blocks have the companion form
0 1 0
0 0 1
(10.31)
0 0 1
X X X
with dimensions v • *v • ,
where the v. are the controll
ability indices of the system; i.e.
0
(10.32)
X X
B = (10.33)
where
(10.34)
0 ... 0 1 0 ... 0
column i
in _ n
n det[si-Ai:L + (BK)i:L] = n (s- Yi ) (10.35)
K = K T (10.36)
K - (10.37)
K = FM (10.38)
186
be considered as
(2) Since r (1) (s) only has rank one, we can define
s=^
as
rU)(yi)f(2) = O (10.43)
f(2) andd(1)(s)
m + I - 1 > n (10.47)
189
n- (10.48)
max{m,£}
N 1 (s)
£(s) = A o (s)
(10.51)
Vs)
190
where _ _ ,
A (s) = s + a.s + ... + a (10.DZ)
^t(s) =
A \s) CM
l(s) '•• * V S ) J (10.54)
M
i(s) = 6
iosr + G
il sr "" 1 +
""• + e
ir (10-56
14
Then, Chen and Hsu have shown that the resulting closed-
loop system characteristic polynomial A, (s), where
x
r Er = £r (10.61)
9X
X =
rr ^r -^-r (10.62)
g
l 17
where X is a g -inverse of the matrix Xr defined by
X
rXr r = X
r (10.63)
91
£r = X r 1 lr (10.65)
(10.67)
r nl (10 69)
) i ± -
We can now see that for the case of state-feedback, since
then i = n, equation (10.69) yields r = 0; i.e. F is a
constant matrix.
r = p(Xr,6_r) (10.70)
F = I f(i) m ( i ) (10.71)
i=l
]i = min{m,U (10.72)
193
F * Fc F Q (10.73)
We can ensure that for the case where m < £ the matrix F
has full-rank by its construction, but F will have full
rank if and only if FQ and FQ both have full rank. It is
tedious but straightforward to show that if F is gener-
ated by the method shown in the following algorithm, then the
corresponding matrix F simultaneously has full rank. The
case where m > £ is obtained by duality.
Algorithm
Let n0
be the number of modes influenced by b = Bf
(2) ~
where 1 < n o < n - 1 . T h e n , f is determined to move qo
2'n2-^ poles to desired locations by solving the equa-
tion (10.74) :
where
r(1>(s)f<2> = o
= [adj (10.76)
10.3.3 Example
Consider a system with transfer-function matrix
s2(s+l)
G(s) = ±j (10.77)
s 0 sz(s+l)
(10.78)
195
So, the system seen through input 1 has two poles at the
origin and one zero at -1. Now, since n 1 = 2 and £ = 1,
we get q = 1. Thus, we determine f (1) to move one pole
to a new location at s = - 2 , say. This is achieved by
fi 1} = [4 0]; i.e.
(10.79)
s (s+l) (s+l)
(10.81)
s2(s+2)2 (s+l)(s+2)'
(s-2)
r<2) (s) = (s+2) (s + l) (10.83)
(s+2)
i.e.
(20s+37)/2 (-97/8)
F 2 (s) = (10.85)
-2(20s+37) (97/2)
, , (28s+45)/2 -97/2
F(s) = -—- I (10.86)
1
-2 (20s+37) 97/2 j
References
Appendix 1
The parameters in the vector matrix equation
Xr£r = «r
I -
•<f <=>
^• • • • * CD
I I
I • • "l 2
ca o
-f J
^ o
o ^
"o
O ~ -3 O ef o o o o
o o ;
tl-
200
5 2-j -6-2j
5 2+j -6+2j
0 - 4 2
0 1 0 0
0 0 1/3
T = 1 4/3 5/3
0 1 -2/3
1 4/3
and Q =
0 1
s +6s+4
G(s) =
s 3 +6s 2 +7s+2 s+2
P.5 Given
s2+s+l
G(s) = ±j
Nonlinear systems
Dr. P. A. COOK
Synopsis
u(t) = U sin(a)t)
a
0
y(t) = y- + I {a cos (nut) + b sin (nut)}
n n
n=l
27T
y(t) cos(na)t)d(a)t)
1 f 27T
bn = -TT y(t) sin(nwt)d(a)t)
JQ
y(t + 1) = - y(t)
just like the input, and hence a and b vanish for even
n, leaving only the odd harmonics (n = odd integer).
y( t) =
^ - y{t)
Example
An ideal relay.
y = sgn (u)
206
i.e.
1, u > 0 ;
y = 0, u = 0 ;
-1, u < 0 .
ira - f 7T
0
c o s (ne)de
-I
ZTT
7T
cos(n6)de = 0
-1 -I
TT 2TT
Tib s i n (ne)de sin(n6)de
0 IT
2 (1 - COS nir)
n
so that
— , n odd ;
0 , n even.
f 27r 2
A = (y(t) - y o (t) } d(wt)
and minimise this with respect to a1 and b^ (regarded as
207
= 2
f
~ lYnW ~ y(t)> cos(o)t)d(o)t)
i Jo
2{TT a, - y(t) cos (wt) d (wt) }
1
Jo
~~ = 2 I {y o (t) - y(t)}sin(ait)d(a)t)
2TV
= 2{TT b - y(t) sin(ut)d(a)t) }
j
0
Y(
N(U,uO = ^a)) exp{j|(U)} 1 ai
so that
y Q (t) = Im{UN(U,uO
u(t) = Im (Uexp(jwt)}
A i_ __L_
s ' 2 ' s+1 '
s
s+1
s2+l
u(t) = U sin(ut) + V
a
0
y o (t) = a1cos(a)t) + b1sin(o>t) + y-
Example
A quadratic nonlinearity,
2
y = u -
= ~- {1-cos(2wt)} + 2UVsin(wt) + V 2
hl = 2UV
aQ = U 2 + 2V 2
b
l
N(U) = ~- = 2V
f(u)
nonlinearity •• 9( S H
Fig. 11.1
u(t) = U sin(a)t)
U = -g(ja))N(U,u))U
-1
g(juO
N(U,a>)
_,., x _ - 1
^VJUJ; " N T U T
In this form, it is convenient for graphical solution, by
plotting -1/N(U) on the Nyquist diagram of g(s), and
looking for intersections.
Fig. 11.2
212
y = h(x,u,t)
x = f(x)
x = f(x)
V(x) = I §2-
3x
fX (X)
i i
V(0) = 0,
Example
X = X
l 2
X =
2 ~X1 " YX
2
V(x) = x^ + 2x2
X X
so that V(x) = 4(x^x x + 2 2^
V(x)
Ax
217
V(x) = x T Lx
LA + A T L = -Q
L = f exp(ATt)Q exp(At)dt
m
V = -xQx
v it:
positive; denoting this by minf follows that the re-
gion
V(x) < V
min
is a domain of attraction.
l+gg(s)
l+ag(s)
219
For simplicity, suppose that & > a > 0 and that g(s)
is the transfer function of an asymptotically stable system.
Then, the circle criterion requires that the Nyquist plot of
g(jw) should not intersect or encircle the disc on
[-a" ,-3 1 as diameter.
(l+qs)g(s) + e""1
model
input
— controller plant
adaptive law
Fig. 11.3
References
f (u) =
and g.(s) =
s(s+l)(s+4)
x = X
l 2
Synopsis^
12.1 Introduction
As a result of the enormous impact of microprocessors,
electronic engineers with sometimes only a cursory back-
ground in control theory are being involved in direct dig-
ital control (ddc) system design. There is now a real need
for an easily understood and simply implemented design tech-
nique for single-input DDC systems which is the objective
of this tutorial paper. The proposed design procedure is
shown as a flow chart in Figure 12.1, and it contains appro-
priate references to sections of text that treat the more
important issues in detail. It is hoped that this diagram
will provide a systematic approach to DDC system design, as
well as an appreciation for the organisation of the text.
entry
design filters;
references 13, 14 I calculate pulsed frequency
note common- I response of plant, a d . and d a .
mode levels I equation 7 — section 3
4
Knowles and Edwards show that the average power of the con
trol error component due to imperfect side-band attenuation
is bounded for a unity-feedback, single-input ddc system by:
wg = 2TT/T (12.3)
(12.4)
*
where H (jw) defines the overall pulsed frequency response
of the stabilised closed-loop system, 0 (w) is the input
x
*
signal power spectrum, and P (jw) denotes the pulsed trans-
fer function of the plant and zero-order hold.
12.3 Frequency Domain Compensation Method
* *
X K(z)P(z) 4 XK (jf)P (jf) (12.5)
z=exp( j2irfT)
p s
P*(jf) = k
1
( ) (12.7)
S=j27Tf
z-a
X K(z) = x (r—) with 0 < a 9 < a, < 1 (12.9)
1.5
; a =
KD(z) = <V*l><i=T> K ^ <12-U>
With these values, the gain and phase of K(s) for frequen-
cies above u>_ are a and zero respectively. Hence,
hi
inclusion of the network improves the Gain Margin by the
factor 20 log l Q (I/a), without altering the original sat-
isfactory open-loop zero-frequency gain. Both the above
design techniques for phase-retard networks are based on
sacrificing phase in a relatively unimportant portion of
the Nyquist diagram for gain in a separate region that mark-
edly influences control accuracy or stability. As the
frequency range over which the network produces an appreci-
able phase-lag widens, this simple design procedure becomes
evidently complicated by the progressive overlap of these
regions of gain and phase changes.
6 = oj - IDQ (12.19)
Where:
1 - exP(-j2o>nT)
U (12.21)
- r Q exp(-j2a)0T)
B = ± { U r
N0 " 0)/T} (12.22)
z
N(z) = K ~aZ+1 7 (12.23)
N
z - rbz+r
coefficient rounding on {a, rb, r } is seen to modify both
the centre-frequency and the bandwidth of a notch filter
design. Ideally with no rounding the coefficients in equa-
tion (12.23) are specified by:
(r2) = (r2)
exp(j6)-l * " e X ? [ -3
exp(j6+je)-r - r exp[-j
(12.27)
N (12.30)
|N
((l-r)/6T)
where
K = K (12.31)
N N 1 - r exp (-J2(JO
= 2 (12.32)
B
NT
(12.33)
N
3B
N
6B
< + Sr (12.34)
NO 3r
0 0
_ 2 "(wordlength - 1) (12.36)
(12.37)
NO
aQ = 0.724716 (12.38)
P(s) = — — ~ (12.45)
1 2 3 A 5 6 7
time s
1 2 3 4 5 6 7
time s
3.0r
2.0
1.0
f T = 0.16 (12.47)
B
fT = 0.03 (12.48)
*K(Z) = * T K T I T (12 49
->
which with a gain constant (x) of 0.4 realises a gain-
margin of approximately lOdB. Figure 12.5 shows that the
step and ramp function responses are now both adequately
damped, and the additional control error due to multiplica-
tive rounding errors is derived from equation (12.42) as:
|e Q | 1 37.5q (12.50)
<b
REFERENCES
1. IEE Coloquium: Design of Discrete Controllers,
London, December 1977.
Synopsis
13.1 Introduction
This chapter shall attempt to look at some aspects of con-
troller synthesis for the multivariable servomechanism pro-
blem when the plant to be controlled is subject to uncertain-
ty - in this case, a controller is to be found so that sat-
isfactory regulation/tracking occurs in spite of the fact
that the parameters of the plant may be allowed to vary (by
arbitrary large amounts), subject only to the condition that
the resultant controlled perturbed system remains stable.
This type of controller design will be called robust con-
troller design in this chapter. It is to be noted that
other definitions of "robust controller design" may be found
in the literature.
The topic of robust control is an extremely active area of
research at present, and so it is impossible to include all
results/approaches obtained to date. Some recent literature
dealing with various aspects of robust control design by us-
ing a number of different approaches may be found, e.g. in
Sain 1 7 ; Spang, Gerhart 19 . No attempt shall be made in this
short chapter to survey all such approaches. Instead, the
chapter shall concentrate only on some of the structural
247
X = Ax + Bu + Ew
y = Cx + Du + Fa) (13.1)
+
m = Hh D
mu V , e = y
" y ref
n
n A
l
l " lnl ' n
i£R (13.2)
and that the reference inputs yref arise from the system:
n.
n2
a (13.3)
y = Ga
ref
n 2 (0)eR ^
and all bounded controller initial conditions.
u = + K y + K y +
V l m 2 ref' * = V A ^ + A 2 y ref
X = AQ0X + A l x ( t - A ) + Bu + EGO
C
mx(t~A> (13.4)
y = cQx ^xCt-A) + Du + Fu)
Definition
I 0B
n
2
rank
I ®C I ®D
n n
2
= rank
I ®C ®G)C
~29
n 2
2
Definition
FC +CT
F
mCl+CmT 1
(F
mCl+CmT)Al = rank
rank
n,-! <FmC mT)A 1
(F C +C T)A n^-1
m l m l
(FC T)A
Definition
Theorem I 5 ' 6
Corollary 1
Theorem 2 1 0 ' 1 3
Theorem 4
Theorem 5 7 ' 3
(ii) m > r
Theorem 6 7 ' 3
hold:
(i)
(Cm,A,B) is stabilizable and detectable
(ii) m _> r
Theorem 7
distinct.
256
<*
i = 1,2,...,v (13.5)
y = C
i ix D
iui F
iw
Definition
Let
C . c1 >
m
B 4 (B1#B2#...fBv), C
m A c2
m
cV
m '
and let
257
0 0 .. 0
0 .. 0
r
l
0 0 .. 0
0 X 0
r2 '•
0 0 .. 0
0 0 .. I
o
Theorem 9
There exists a solution to the robust decentralized servo-
mechanism problem for (13.5) if and only if the following
conditions all hold:
Definition
Let 0
P+l
Theorem 10 14
Assume that the plant (13.1) is open loop asymptotically
stable; then there exists a solution to the robust decen-
tralized servomechanism problem for (13.5) with sequential
reliability if and only if the following conditions all hold:
A
l ' X 2 / # * * / X p' X p+l"
(ii) there exists a solution to the robust decentralized
servomechanism problem for (13.1) with control agents
1, 2 for the class of disturbance/reference signal
X X
poles 1' 2 ' * * *' A p' X p+l*
Definition
xt = A x t + Bu
, dx t
de , e e [-A ,0)
Ax.
A Q x t (0)+A ± x t (-A) , e =0
0 [-A ,0)
CD
u)
Bu
Bu 0
CD
II
Theorem 11
(ii)
rank = n, V A e sp(A) with Re(A)>0
(iii) B
rank -A, A = n+r, i = 1,2,...,p
D
m
(iv) y is measurable.
Definition
c 4
pJ
where the coefficients ±, i = l,2,...,p are given by the
coefficients of the polynomial (X-X i ) # i.e.
A P + 6 A*
n
The following compensator, called a servo-aompensator is
of fundamental importance in the design of robust controllers
for (13.1):
Definition (Servo-Compensator)
* (13.7)
= C
B e
where
C 4 block diag(C,C,...,C), B A. block diag(y/ ,..,y)
r r
u = KQx + K (13.8)
0
(13.9)
EC C B D
C
m
u, u = (KQ,K)
* B D
B C
(13.10)
= (13.11)
V K
lym
12
Theorem 12
Assumptions
(i) The plant (actual) can be described by a linear time-
invariant model and is open loop asymptotically stable.
(ii) The disturbance/tracking poles have the property that
Re(X ± ) = 0 , i = 1,2,...,p.
(iii) The outputs y can be measured.
Definition
The tracking steady-state gain parameters of (13.1) are
defined as follows:
Remarks
2. The tracking steady-state gain parameters of a system
can be easily identified by sequentially applying m linearly
264
Existence Result
o
Theorem 13
There exists a solution to the robust servomechanism pro-
blem for (13.1) if and only if
rank T^ = r i=l,2,...,p
Controller Structure
Controller 0
u = K(yret^-y) + {K,
l + Ku
A (sI-A n )
u A,}y
l •* m, y m Ar" (ym-D mu)
Controller 1
u = -( " " . . - ! - -
Controller 2
u =
1 2 2
Here, K/s, K +K s/s is a servo-compensator for (13.1) for
the case of constant, ramp signals respectively, and
K , + K Q ( S I - A Q ) "* A.^ is a stabilizing compensator for the sys-
tem.
Let e A, y ~-y ; then the following steady-state errors
are obtained for the controlled system when reference sig-
nals y yty1 t 1 , i = 0,1,2
r; = where y is a constant, and
_ i _
disturbance signals u = wt , i = 0,1,2 where w is a
constant, are applied to the plant controlled by controllers
0,1,2 (see Tables 13.1, 13.2).
266
=
y v ^ref=? t2
Controller 0 G y 00 00
t
0
Controller 1 0
V 00
Controller 2 0 0
where:
G.r A {I+C[A+B(K -K n A~ X A )Cni]'"1BK}"1 Tracking Position
0 l u u I m Error Gain
Tracking Velocity
Error Gain
. A " 11 A 1 )C ]~
2{C[A+B(K 1 -K n A" ~ 11BBK
K 11}}"
" 11 Tracking Accelera-
2
'2 tion Error G a m
where:
C A C + D(K1-KQAQ A1)C.
m
TABLE 13.2 Generalized Steady-State Errors for
Various Disturbance Signals
-4.2 1
Controller 0 00
; Controller 1 0 G j 0) oo
d
l
Controller 2 0 0 G
d ^
where
F 4 F + D(K 1 -K 0 A- 1 A 1 )F m
267
Problem Statement
u = K Q y m + K i n + K£ , n = AQn + A ^ (13.12)
Definition
Let z A D£ where £ is the output of the servo-compensa-
tor (13.7) where
Performance Index
Let
r °°
I nil *P
J £ E (AzTOAz+yAu RAu)dx , Q > 0 , R > 0 ,
J 0 (13.13)
-2
10
opt
0.01
10
' ^0.02 0.04
/ -0.01 t.s
0
0.02 0.04
t.s a
-0.005
0 "~" 0.02
t.s
0.04
b
-f 0.02 0.04
t.s
0005 • 0.005 -
(c)
*ref ) , E(D = (1 1 1 1)
13.6 Conclusions
References
Synopsis
14.1 Introduction
Example 1
2
| | (x,t) = a i - | (x,t) + b(x)u(t)
0(0,t) = 0U,t) = T
0(x,O) = 0Q(X)
0(x,t) = T(x,t) + T
e Q (x) = T Q (x) + T
then
T(x,0) = T Q (x)
h(t) = f c(x)6(x,t)dx
J
0
If h(t) = y(t) + T
- [[ c(x)dx,
i
275
then
y(t) = c(x) T(x,t) dx (14.2)
Jo
Equations (14.1), (14.2) represent a single input-single out-
put distributed parameter system.
The results to be described are more easily formulated by
abstracting the system. To do this let us ask what proper-
ties one would expect of the solution of a dynamical system
on a Hilbert space, H. If it is assumed that the dynamics
which govern the evolution from the initial state z to
z(t) are linear and autonomous, then for each time t, we
can define a linear operator T., such that
TQ = I (14.3)
z(t) = T t zQ (14.4)
T z = T
z(t+s) = t+s o tz(s) = T
tTszo
T = T
t+s t Ts (14
*5)
with solution
ft
y(t) = T.
<c, T b> ul
t _ p b> u(p)dp
Az = —
ax
3z
D(A) = (zeH : *-% eH, z (o) = z (£) = 0)
(sI-A)z = Z]L
d2z
sz(x) - cA-~(x) = z. (x)
Z 1
z(0) = zU) = 0
The solution is
sinh(s/a)^(x-a)z (a)da
c(x) I
sinh(s/a)^(£.a)b(a)da
fx 3k 1
sinh(s/a) Mx-a)b(a)da dx
J
0 JJ
It is evident that this is not a rational function of s
so that new theories for irrational transfer functions need
to be developed for distributed parameter systems. An al-
ternative approach is to show that the solution is equiva-
lent to the solution of an infinite number of ordinary
differential equations. Let us suppose that the spectrum
of A consists entirely of distinct real points A^ ,
i = 1,2,.... where
e D(A*)
* (14.9)
z(t) = I a n (t)* n
then
a bu
z(t) = I a (t)4 = A I (t)<f> + < fc )
n=l n n=l
n=l
I Vn(t)*n + bu(t)
If we assume
b c =
= n=l
^Vn ' J Vn
n=l
°° 2 °° 2
where £ b < «> , £ c <»
n=l n=l
then by taking inner products with $ we obtain the un-
coupled system of ordinary differential equations
Jl
For Example 1 we have 4 (x) = (>
j (x) = /ijl sin
0 0 0 n n
X = n Z ir Z a/£
14.2 Controllability
Let us consider a question of controllability for Example
1. If the initial state is zero, we have
ve
-f "f
(t-s)
n
an(t) = b u(s)ds
n
J0
Jn
Suppose we wish to steer the system to a temperature distri-
bution T.(x) at time t., where
279
then t. - .
1 x (t
n l~s)
e
b n u(s)ds
9 ^1 2X (t.
(t.-s)
b e ds
u 2 (s)ds
nl J j
for n = 1,2,...
or
2|tJ2 r .
\ | < f u i(s)ds for A < 0
n
*" In
where b z = <b,z>.
1
T g b b* T g z> ds >• Y| | Z | I 2
b* T* z = 0 on (0,t1) (14.11)
implies z = 0
For Example 1
b* T* z = I b exp(-an2Tr2t/£2)a
r
n=l
where z = } a <j>
n=l n n
So Example 1 is approximately controllable on (0,t ) if
b ¥ 0 for all n = 1,2,..., where
n
281
14.3 Observability
implies z = 0
= G
= R -1
z(t) =
z(t) = V -1 T
t-s Q(s) c dp
Q(t) =
J
j=l k=l "
Q(t)
(14.17)
V R
P G
References
[l 2 \tl 2
z
J(u) = J T (x,t_)dx
1 + j0 u (t)dt
J0 1 j
y(t) = f 1 -c(x)w(x,t)dx
j
0
Decentralised control
Professor MADAN G. SINGH
synopsis
15.1 Introduction
N
x(t) = Ax(t) + I B.u.(t) + £(t) (15.1)
i=l 1 1
N
yjt) = cx(t) + I D.u. (t) + jKt) (15.2)
i=l
H =
H = I
u±(t) =
A lot of work has been done in the last few years to seek
out special cases where linear solutions are possible. Ko
and Chu found that the so-called "partially nested informa-
tion structures" lead to decision rules which are linear in
information.
y = K K + D u
i i i
k e E.
z± = 1 (15.3)
where E± = { k t D i k ¥ 0}
A lot of the work that has been done on design methods for
stochastic decentralised control ' ' has relied on seeking
fixed structure controllers which are usually constrained to
be linear. Parallel to this, a lot of work has been done on
the deterministic decentralised control problem. Next we
will examine the basis of this work,
V
x(t) = Ax(t) + I
(15.4)
= C
i ~(
B i y i (t)
(15.5)
u ± (t) = Q i Z i (t) + K ^ U ) + v ± (t)
i = 1,2,... v
m.
wa
Va/ —La •• g J * «& ml~
where z. (t) e R " and v. (t) e R is an external input.
Definition
Comment
Definition
A(C,A,B,K) = f[ a (A+BKC)
K e K
Remark
The reason for the name fixed modes is now clear, since
these modes do not depend on the values assumed by K c K
293
iHX,C,A,B,K) =
Theorem 15.1
(C,A,B,K) a c"
12
15.4.1 The algorithm of Geromel and Bernussou
=
x, (o) = xio i 1/2,... N
where
the form
I
A r.
x^ c R "*", u^ e R . This system can be rewritten in
N N
I n^ r = I
we seek a linear static decentralised feedback controller of
295
the type
Hi = - k±x± ; i = 1,2,... N (15.9)
1
k. e R }
min (x Qx + u Ru)dt
J
Keft(s) 0
s.t. x = Ax + Bu
u = - Kx
T
where it is assumed that Q = CC ^ 0 and R > 0 .
f00 m m
J = x(o)T <T(t)(Q + KTRK)(j>(t)dt x(o)
0
We can write a new cost function by considering x(o) to be
a random variable which is uniformly distributed over an
n-dimensional sphere of unit radius, i.e.
min J(K) = tr P
KtR(S)
H i £ ! = 2(KK-BTP)L
(A-BKTP+P(A-BK)+Q+KTRK = 0 (15.10)
(A-BK)L+L(A-BK)T+I = 0
Lemma 15.1
Remarks
a = va otherwise
J
N ,
lfi±l R. }dt
s.t. (15.11)
N
min J = i I }dt
i-1
R.
s.t. A B (15.12)
i i Hi
where
-i
1
Hi = Ki *^! (15.13)
x± = (A± - B.G il )x i + C ± z ± -
300
or
x± = A.x. + C.z. - E.G. 2 z. (15.16)
A = A
where i i ~ B i G il
x
_\ is the suboptimal state of subsystem i resulting from
—i
using z^. instead of _z. and
= A + C
~i i-i i-i " B i G i2-i (15.17)
z. = z. - z. (15.18)
* * —1 T
z. = z. — S. C.K.x.
= (A±
x± = A.x. + ( C ^ B ^ ^ ) ^
ii = A zii
and we could represent these equations by Figure 15,2 •
C:Z:
-B|
Fig. 15.2
References
Chapter 1: Solutions
S.I.I H = 1 -2 2
0 -2 0
0 -2 -1
1 -1 -2
1 -2 2 2
1 0 0 1
2 1 0 2
0 1 -1 0
r T -T
-1 0
0 -1
0 0 "
1 0
0 0
0 1
305
R2 = A + 41 2 1 1 0 "
' a
2 = - ~tr(AR 2 ) 9,
1 4 0 -1
1 0 3 -1
0 1 1 3 _
AR1+1OI 1 1 0 - 1 = - £tr(AR0) = 3
-1 5 - 3 -2
0 - 3 3 0
- 1 2 0 1
4s2+8s+3 - (2s3+7s2+7s+2)"
G(s) =
s 4 +4s 3 +9s 2 +10s+3 2s 3 +5s 2 +3s 4s +8s+3
306
Chapter 3: Solutions
WG(s)W -1
2
s+2
t 8
where W =
6 7
1 -8
-1
and W
-6 •7
and
7
1
/85 -2
6 /ITT
1 7 6
/ITT -8 /85- -7
s
0.0102
/85 /fl3
s
—=^ 0.0102
/85
Chapter 5: Solutions
S. 5.1 We have
|TGI|H|/|G|
(s+D+l s+1 0
0 s-1
(s-1)(s+1)
0 s-1 ' 0 1
P o (8) =
-1 0 ' 0 0
0 -1 ' 0 0
gi;L(s) =
g2 ' h2 * Also p Q = 1, so
ZN__. - IN = 1 = Pr
hi
S.5.2
(i) For a conventional feedback system with unity feedback.
iir ° 1 -1
G(s) =
1 1
Therefore,
1
1
1
f[-1 M
x)
s+1
0
0
s+2
(s+1) (s+2)
1
2
-(s+1) (s+2)
0.5 1.0
-0.5
y
22
1.0
INA of given system
An examination of this inverse Nyquist array reveals that
for row diagonal dominance - row 1 is not dominant
row 2 is dominant
and for column diagonal dominance - col 1 is not dominant
col 2 is not dominant
although in the latter case, we are on the margin between
dominance and non-dominance.
In either case, dominance can be determined graphically by
superimposing circles of radius
-1
K =
1
which results in
s+1 0
Q(s)
0 s+2
1
1
Thus, K = K"1
2 -1 1
i jUJ JUJ
-2
Chapter 7: Solutions
S.7.1 H = ;^x 2 + X2(-x
Canonical equations
h - - \
{
2 =
For u = 1,
x = X
l 2 ' X2 =
~xl+ 1
dx X
l 2 2 2 2
i.e. gj- = — f whence (xfl) + x2 = R
i.e. Circles with centre (1, 0 ) , radius R.
x
Also, i~1 = ~ R cos(t+Y) and x 2 . = R sin(t+y)
X^ +k22 = R2
x = X X =
l 2 ' 2 "xl ~ 1
2 2 2
Phase trajectories are (x^l) + x 2 = R : Circles centre
(-1, 0)
(p odd).
s n X + U
H = X, cos x 3 + A 2 i 3 *3
X, = - H X i = 0, X2 = -H^ = 0
X 3 (t f )
If X3 > 0 then u = -1
x3 < 0 u = +1
x
Suppose u = +1. Then x- = t, o = l ~ c o s t, x, = sin t.
H(tJ = - ?2- , G
f 3t f
Then a = b = --.
/2
315
On singular subarc
x3 = \ , x2 = i +£ (t - i - ?-)
x ( t
i - ^
x = anc X 1
l ^f ^ ^ ^ 2 ^f ^ ~ satisfied when
u = C. J < t < 1 + J .
u 3 u 3 1 3 2 3
H = - i ( - + i ) = -1 < 0
u /2 /2 /2
316
f =
then
bf = Bf =
and , A
0 3-2
1 0 1
1 2 - 2
= <v
Pi -i' bf > = -4 - j3
_ -.5(3 + j'
" (3 - J4)
= -0.14 + jO.48
m (
"1"4 0#4 3
'6)
A O (S) IsI-Al = s 3 + 2s 2 + s + 2
{-2, - 2 , -2}
6(s) = 4s 2 + 11s + 6
=
S.10.3 For (A-BK) to have eigenvalues (YJ_) {~4/ - 4 , -4}
we proceed as follows. The matrix BK has the parametric
form
0 0 0
BK = k k
K k
l 12 *
k.2l
, k
22 k
23
-80/9+k lx s-26/3+k 12
i.e.
224/9 50/3 10/9
K =
0 0 13/3
K Q"1 K T
319
1 4/3 0 0 1/3
K 1 4/3 5/3
0 1 0 1 -2/3
0 13/3 -26/9
A d (s) s 3 + 7s 2 + 14s + 8
we have (for r = 0)
^ = ( 1 7 6 )
Here, 1 0
6 1
4 2
Now, a g^inverse of XQ is
320
1 0 0
-6 1 0
therefore, p = X x 6
o —o
A d (s) = s 3 + 6s 2 + 11s + 6
which yields 6 = ( 0 4 4)
r 4 3 2
Then, since A (s).s = s + 6s + 7s + 2s
o
we obtain (4 28 48 24)
1 0 0 0
6 1 1 0
and
4 6 2 1
0 4 0 2
A g^inverse of is
1 0 0 0
0 0 0 0
3 -1 h -i
-9 2 i
-7 2 -1 1
Y
• l • 4 '
9
10 0
e
ll 2
X
El = «i
1
|
i
1
2
CD
to
O
4
CD
to
r-
i.e. m t (s) = [2 2(s+2)]
0 -1 f 3
l
1 -3 f = 2
2
0 -2 e
l
0
322
which yields e1
0 0
i.e. f (D f (D
—c —o -6 -7
s2+s-6 1
(s+1)(s+2)(s-3)
s +s+6
fi 2) r (1) (-D - o t
i.e.
[1 1]
Thus, N(2)(s) .
= [2s(s+l) (s+1)]
2 0 f
l 8
0 1 f
2
i.e.
—c 16
(2)f(2)
4 4
=
Thus, F2 = f
z —c —u 16 16
323
Now , 4 4
10 9
- sin(wt)
IT
N(U) = iL
.*. the graph of -1/N(U) lies along the negative real
axis, tending to -« as U -* °°.
The Nyquist plot of g(s) is obtained from
U = 1/5TT
and p e r i o d
2TT
325
V =
Lyapunov function is
326
W = 0 , G = I, R = l
p
mn ( t l> = X m = n
= 0 m f n
= 0
So
P n n (t) = exp(-2-rr2an2/£) (t± - t)
For n = 1
or 2
2
at
— = Ae _
But p l x (tx) = 1, so
327
-t)
f* ,
= -p (t) /2/£ sin TTX/A z(x,t)dx
V=l, c = b , P = I
9t
•V
f
y(t) = I /2/fc sin TTX/£ T(x,s)dxds + v(t)
T(xf0) = TQ(x)
ft
T ( x£ ,/ t ) =
= I /2/£ e P n ^ , - 3 ) s i n TTX/£ dp ( S )
in
where
f ]
p (t) = y(t) - J2/1 sin TTX/£ z(x,s)dxds
f
328
2-
s 2 w(x,s) = i~ (x,s) + b(x)u(s) (14.18)
y(s) = c(x)w(x,s)dx
Then
s 2 w(x f s) + b(x)5(s)
(a+3 - I e s p b(p)dp) = 0
(a ++
(es(a
J
[ e" s p b(p)dp) + e" s 0) = 0
0
so that
l-e zs
f
J
0
(e sp + e- sp )b(p)d P
J
1-e " 0
+ e" s x (e - [ e s p b(p)dp) dx
•'x
In order to find the semigroup we have to solve the
329
y = <c,z>
c(x) = (c(x),0)
Az =
and dx
D(A) {z e H J ( O , 1 ) , Z 2 e H 2 (0,l) A H*(0,l)}
w(x,t) = T(t)X(x)
d 2 X(x)
then (t) =
T(t)" d t 2 X(x)
dx 2
So that
T(t) = A sin At + B cos At
mrA = /2
n
330
So
w(x,t) = cos n7Tt > sin
n=l n
§?(*,t>
group is given by
T z = T (<w o , V cos
t o t n 1
sin nirtH
(-mr t+<.w 1 f4>R>
c o s mrt)<|)
Realization, 11,13
Reference input, 248,255
Regulation, 247
Regulator, 2 63
Reliability, 247
Robustness, 247
Root locus, 21,25,32,34,36,196
Schwartz inequality, 279
Second variation, 128,133-135
Semigroup, 275
Sensitivity index, 49,60
Servomechanism, 246,248
Singular
control, 133
value, 56,147,150
value decomposition, 55
Spectrum, 20
Design of
MODERN CONTROL SYSTEMS
About the book The Editors
The IEE Control Engineering Series of David Bell obtained his B.Sc. in
which this book is a member contains a Mathematics (1955) from the University of
volume (number 9) devoted to the modern Manchester and his M.Sc, Ph.D. degrees
approaches to control system design. The (1962,1970) in Mathematics from the
present text may be considered a University of London. From 1955 to 1961
companion volume which presents new he worked as a research mathematician for
the aircraft firm A. V. Roe and Co. Ltd.,
developments in analysis and design followed by a period as lecturer in
techniques for control systems. Included Mathematics at the Bristol College of
are the exciting new results for feedback Science and Technology (1961-1964), and
systems using complex variable methods the University of Bath (1964-1967). His
as developed by Professor A. G. J. present joint appointment between the
Macfarlane and his coworkers at the Department of Mathematics and the
University of Cambridge, the important Control Systems Centre at UMIST was
concept of robustness in controller design taken up as a lecturer in 1967 and senior
reviewed by Professor E. J. Davison of the lecturer in 1973.
University of Toronto, and the increasingly
important topic of decentralized control for Peter Cook obtained his B.A. in
large scale systems as seen by professor Mathematics (1962) from the University of
M. G. Singh of the Control Systems Cambridge and his Ph.D. (1966) in
Centre, UMIST. These and many other Theoretical Physics from the University of
Manchester. He was then appointed
contributions illustrate the great activity and
lecturer in Applied Mathematics (1968) at
rapid progress which has taken place in the University of Kent. From 1971 to 1974
the subject over the past few years. Only he was a Senior Research Fellow in the
by bringing these contributions together Department of Mathematics at the
under one cover can the practising Loughborough University of Technology.
engineer in industry and indeed the Since 1974 he has been a lecturer in
engineer in university or polytechnic keep Control Engineering at the Control Systems
fully informed on the “state of the art” on a Centre, UMIST.
number of different fronts. Application of
the theoretical developments and practical Neil Munro obtained his B.Sc. in
aspects of the subject are not forgotten; Electrical Engineering (1962) from
analysis and design of a nuclear boiler and Strathclyde University and his M.Sc. (1967)
some direct digital control system design and Ph.D. (1969) in Control Engineering
procedures are but two topics discussed in from the University of Manchester. He was
the present book. Several of the chapters then appointed as a lecturer (1969), senior
lecturer (1973) and Professor of Applied
are followed by problems on the subject
Control Engineering (1979) in the Control
matter and worked solutions to most of Systems Centre, UMIST. He has also spent
these problems are given at the end of the some 10 years in industry as an electronics
book. This aspect will find favour with technician, development engineer and
many readers since such contributions are design engineer working for Barr and
often a great help in the understanding of Stroud Ltd., Ferranti, Parkinson-Cowan
the subject matter. Measurement and ICL.