Computer Control of Process 2UPs
Computer Control of Process 2UPs
1 Figure 1.1 shows the basic elements ofa typical closed - loop control
system with sampled data. The sampler simply represents a device or
operate that outputs a pulse train. No information is transmitted between '.
two consecutive pulses. Figure 1.2 illustrates typical input and output -of
a sampler. A continuous input signal e(t) is sampled 'by the sampler, and
the output is a sequence of pulses. In the illustrated case, the sampler is
assumed to have a uniform sampling rate. The magnitudes of the pulses
at the sampling instants represent the values of the input signal ett) at
INTRODUCTION the corresponding instants,
In general, the sampling scheme has many variations, some of these
are the periodic, cyclic-rate, multirate random, and pulse width-modulated
1.1 . INTRODUCTION samplings. The most eom~on type of sampling found in practical systems
is the single-rate periodic sampling, such as that shown in Figure 1.2.
In recent years significant progress has been made in 'the analysis
and design of discrete data and digital control _systems. There systems c(t)
r(t)
have gained popularity and importance in industry due to in part to the
advances in digital computers for controls and,more recently, in
miroprocessors (MP) and digital signal processor.
Discrete - data and digital control system differ from the
.Fig. 1.1: Closed - loop sampled-date. control system
conventional. continuous - data ~n analog ,systems in that the signals in
one ormore parts of these systems are in the form of either-pulse trains e(t) e~ (t)
We shall show that, from an analytical stand ..p oint, both types of Future Generation Control Systems
systems ,are treated by the same analytical tools, The study of emerging applications 'shows that Artificial Intelligence
Digital __- - _ (AI) will affect the design and. application of control systems, as profoundly
code Digital..; To-
input Analogconverter Output as the impact of micro processors in the last two decades. It is clear that
future generation control systems will have a significant AI component,
the list of, applications of computer-base control will continue to' expand.
Analog - To-
Digital converter
The ability to 'redesign' the contoller by changing software (ratherthan Digital- To-!
Output.
hardware' is a.n important feature of digital control against analog control). , Analog converter
However'; with the advent of inexpensive digital computers with In signal processing, sampling is the reduction of a continuous-
virtually limitless computing power, the techniques of modern control theory. time signal to' discrete time signal. Sample is a value or set of
.may now be put to 'practice. For example, in multivariable control systems value, at a point in timaand/or space.
with more than one input' and one output modern .techniques for optimizing
system performance or reducting interaction between feedback loops can
now be implemented.
2.2 ComputerControl of Processes )
The q outputs of the systems are related to the state variables and
gk[Xl (t), X2
.:W i
where x(t) is an n x 1 column matrix and is called the state vector; that is.
2.2 STATE EQUATIONS AND STATE'TRANSIT EQUATIONS OF CONTINUOUS
DATA SYSTEMS
Consider that an nth-order continuous data system with 'p'; inputs ... (2.5)
and 'q' outputs, as shown in Figure 2.1 is characterized by the following
set of n first-order' differential equations (state equations)
The input vector U(t) is a P x 1 column matrix and is denoted as
dXl(t)= f .
[Xl (t), X2 (t), Xn (t), U 1 (t), U 2 (t), ', Up (t), t] (2.1)
i
dt
i = 1,2, ... n
U(t) -- ... (2.6)
where Xl (t), X2(t), ... , Xn (t) are state variables, U 1 (t), U 2 (t) ... Up (t) are
input variables, and h denotes the i th functional relationship. III general,
n can be linear or. a .nonlinear function.
2.4 ComputerControl of Processes ) ,
fTheState-variable Technique : 2.3
To determine state transition matrixott), we take_laplace transform
The output vector C(t) is defined as
on both sides -of Equation (2,.jO ), we get
S xes) - x(D) = A xes)
Solving, Xes) from last equation, we get
OCt) ... (2.7)
xes). = (SI - Ar" x(o) ... (2.13)
Cq(t) - where we assume that matrix (81 - A) is. non singular. Taking the
inverse laplace transform.
Which isa q x 1 column matrix or vector.
On both sides of the last equation yields
If the system is linear and time-in variant the dynamic equations xtt) = L-1 [eSl - A)=-l] xeD) t ~ 0 ..(2.1)
are written as
Comparing Equations (2.12 ) .and (2-.10 ), the state transit matrix is
, dx(t) written
A x(t) + B U (t) ... (2.8)
dt ... (2.14)
and An alternative method of solving the homogeneous state equation is
C(t) = D xCt) + E U (t) ... (2.9) to assume a solution, as in the classical method of solving linear ordinary
differential equations, and. then substitude it in the state equation to show
where A is an n x n matrix, B is an n x p matrix, D is a q~ xn matrix,
that it is indeal a correctsolution.
-E is a q x P matrix, all with constant coefficients.
Let the solution of equation (2.10) be
2a3 -lH,ESTATE TRANSITIDNMATRIH SOLUTIONS OF HOMOGENEOUS STATE x(t) = eAt xeD) (2.15)
EDUATIDNS - for t :2: 0, where eAt represents the power series of the matrix At, '
The state transition matrix, <I>(t) is defined as . n x n matrix that 1 2 2
-satisfies the homogenous state equation, = I+At+ + ... + ... ... (2.16)
2iAt
dx(t) = A x (t) Taking the derivative of last equation with respect to t, we get,
... (2.10)
- dt
d(e At ) -
Thus 4>Ct) statisfies the equation ---
dt
=' Ae At
<!>(tl - t2) <I> (t2 - ta) <!>(t 1 - ta) for any t l , tz, ta' 2.'6 STATEEDUATIONSO,f 'OIS.CRETE-DAJASYST,EMWITHSAMPLE,AN:D -
HOLD
2.5 SO,LUTION OF NON HOMOGENEOUS STATE EOUATION
A. typical multichannel discrete data system with sample-and-hold
The state transition equation (S'I'H) devices or zero-order holds (Zoh) is shown in Figure 2.2.
The. n~n homogeneous state equation
e1(t) - - _ X
u1{t) c1 (t)
dx(t) . , Linear system
- - = A x(t) + BU(t) ... (2.18)
u2(t)
dt e2(t) ,- _-X ,~(t),
Using the definition of state transition matrix in .equation (2.14 ) The outputs of the ZOH are described by
and the convolution integral, the last equation is written as UiCt):::: Ui (kT) = ei (kT) kT ~ t <k +i ... (2.23)
t
where k = 0,1,2,and 1,2, ... p
x(t) Ht)Jl:(O) + fHt-'t)BUCr)d't t~O .. (2.21)
Let the dynamics of the linear system be 'represented by the state
o
equation in equation (2.-8). The state transition equation is given by
This; solution of state equation' is called state transition equation. equation (2.21) for all t and to. Since the inputs are constantbetween any
By setting the initial time at t" = to and using x (to) as initial state we , two consecutive sampling periods, the input vector U(T) in equation (2.21)
can show by using equation (2.20) and some of properties of q,(t) that the has, the properly
state transition equation is written generally as IJ( 't)" == TJ(kT) for kT ~ t < (k + 1) T ... (2.24)
Thus, in equation-(2.21) Utt) can be placed outside the integnal, and
t .
x(t) = <p(t-to)x(to) + f~(t- 't)BU('t)d't ...t .~O ...( 2.21) equation (2.21) becomes
to ~
2.8 computer Control of Processe:!}
.(The state..\ta.riable Technique . 2.7
T
aCt - kT) f<P(t-1:)Bd1: ... (2.27) The discrete state equation in equation (2.29) along wiht the output
kT equation in equation (2.33) are defined to be the dynamic equations of the
system in Figure (2.2).
, Thus equation (2',26) is written as
xtt) = <I> (t- kT) x (kT) + e (t -kT) U (k'I') .~. (2.28) In certain situations, the sampling operation is not done wihtrespect
th
to time, and the Iatter is not the independent variable. Under these
Notethat although the i component of U(kT) , Uifk'I') is held constant -
circumstances the discrete dynamic equations describe the discrete
only for thetimeinterval from t =k to t= (k+ 1)T, the solution in.equation
'events' and are expressed as
(2.28) is valid for the entire interval including t = (k + 1) T, since x(t) is a
continuous' function of t. For numerical iteration, it is more conveilientto x (k--L) = <I> (1) x (k) + e (1) U (k) . ... (2.34)
describe x(t) only at sampling instants. We let t = (k + 1) T, and equation C(k) = D x(k) + EU(k) ... (2.35)
(2.28} becomes .
which are obtained from equations (2.29) and (2 ..33), by setting T = 1.
x[(k _+ 1) T] = <f>(T) x(kT) + 8(1') U(kT) ... (~.29) Another way of looking at the notation is that" in equations (2.34) and
where from equation (2.16) (2.35) the sampling period T has been normalized to be one.
th(T) =. eAT = I +AT + l.-A2T2 + l.-A3T3 + ... ... (2.30) 2.1 STATE EQUATIONS OF DIGITAL SYSTEMS WITH ALL-DIGITAL ELEMENTS
'I' . 2i 3i
From equation (2.27) When a digital system is composed of all-digital, the system may be
described by the following discrete dynamic equations
(k+l)T
e (T) = f <p[(k+l)T-1:]Bd1:
KT
... (2.31)
x (k + 1) =' A x (k) + BU (k)
C(k):: D x (k) + EU (k)
... (2.36)
... (2.37)
By setting m =- ~T -'t, we get dm = dr, (K + 1) T ..;.. 't =T- ID, the last
equation is written as
(TheState-Variable Technique 2..9 ( 2.1'0 ) Computer Control" of pr8ess~
The difference between (2.29) and (2.36) is that, in the for:mer, the Since from the properties of state transition matrix
coefficient .matrix <!>(k) is the state transition. matrix of the A matrix of <p (tl - t2) <I> (t2 - t"s) <p (tl- ts)
the continuous-data p~ocess, and 8Ck)js. a function of $(t) and B; thus 'they
:. Which implies <p(T) <p(T) - <p(2T)_ ... (2~40)
must all confirm with the .properties of the st~te transition matrix'. In
equation (2.36), however,' the matrices A.anda can be arbitrary and depend Equation (2.39) can bewritten as
entrely on the characteriatics of the digital system.
x(2T) = $(2T)x(O) + ~(T)8(T) U(O) + e(T) U (T) ... (2.41)
The solutions of the discrete state equations in equation (2.29), (2.34) x(NT) =<f>([(N-l) T] x [(N-l)Tl + e [(N -1)T] U [(N - 1)' T] ... (2.42)
or (2.36) are the state transition equations of the discrete - data system, Substituting the first N - 1 equations for x[(N-l)Tl; x[(N-2)TJ ... x(T),
Since these equations of all same form, their solutions will be very similar. into last equation sucessively, we' have
However, equations (2.29) and (2.34) contain statetransit.ion.matrix N-l
<peT) of the A matrix of equation (2.8), which is .nonsingular. Thus, ~he x(NT) = cj> (NT) x(0) + L<t>[(N-i-l)T]e(T)U(iT) ... (2.43)
.transit,ion of state for the systems described by equations (2.29) and (2.34) L=O
can take place in the forward or backward direction.
which is the solution of equation (2.29) given the initial state x(O)
For the discrete state equation' in' equation (2.36), generally, the and the input U(iT), for i = 0, 1, ... N-l.
matrix A has no restriction. However if the state transition is to progress
Inderiving the equation (2.43) we have use the property
in the negati~e direction. A must be non singular.
We'.shall derive the state transition equation of the state equation <p(T) <p(T) = <PCT) = <p(NT) .... (2.44)
in equation (2.29) and then extend these to equation (2.34) and (2.36). which follows directly from the property of the state transition
matrix <per), equation (2.44) implies that if <P(t) is a state transition matrix
.2.8.1 The Recursive Method satisfies all the properties of state transition matrix, the multiplying <P(T)
by itself N times is' equivalent to replacing all the T in <p(T) by NT.
The most straight, forward method of solving equation (2.29) is by
, recursion. The following 'equations, are written by substituting The time frame of equation (2.43). Equation (2.44) implies that if(T)
k = 0, 1, 2, ... in equation (2.29) sucessively, is a state transition matrix that satisfies all the ,properties of state
transition Matrix, then multiplying <p(T) by itself N times is equivalent to
k =0 replacing all the T in$(T) by NT .
x(T) ~(T) x(O) +, SCT) U (0) ... (2.38) The time frame of 'equation (2.43) c~n be shifted forward by any
positive integer M, So that solution of x[(N+M)] T is given in terms of the
k = 1
initial state x(MT) and inputs U(iT) for i=M,
x(2T) ~(T) x(T) + SCT) U(T)
M + 1 ..... N +M - 1. 'I'his is accomplished by solving x(O) from
= ~(T) [<!>(T)x(O) + H(T) UCO)]' + SeT) U(T) equation (2.43).
<p(T) <p(T) , xeD) + <pT S(T) UCo) + S(T) U(T) ... (2.39) N-l
xeD) cj>{-NT)x(NT)- 2:<1>[( -i-l)TJ8(T}U{iT) ... (2.45)
T _1'1 '
( 2.12 ) computer Controlof Processes )
Where we have utilized the properties of, 2.8.2 The Z transform Method
[4>CNTJ]'-l The Z transform method can be applied to solve the linear time-in
<PC-NT
varios discrete state equations in equation (2'.29)" (2.34) and (2.36). This is
$(-NT)$. ~IN "- i - DT] = <I> [(-i-t)TJ ... (2.46)' carried' out in the following for Equation (2.29).
Replacing N with N+ M in equation (2.45), we have The state equation in equation (2.29) are represented as
. N+M-l
x[(K+l)T] = q, (T) x(KT)+8(T)U(KT) ... (2.29)
x [(N+M)~J = <I>[(N +M)T]x(O)+ L <t>[(N+M-ic-l)T]8(t)U(iT)
i=O TakingZ transform on both sides of the last equation, vee have,
.., (2.46) z.~x(Z) = z x(O) = cp(T) X(z). + e (T) Ll(z) h_ (2.53)
Substituting equation (2.45) into the last equation and simplifying
where x(Z) is defined as the Z-transform of the state vector x(KT)
'. the .state transition equation for any initial state x(MT) is written as
and is given as
N-l '
00
x [(N+M)T]'=cI>[(NT)T]x(MT)+L<t>[(N -i-l)T.J.S{t)U[(M + i) TJ x(Z) Lx(KT)Z-K ... (2.54)
i=O'
K=O
... (2.46)
and the same applies to U(z). Note'that taking the Z-transformof a matrix
The solution of Equations (2.34) and (2.36) are easily obtained and is equivalent to taking Zvtransformofevery element of the matrix. Solving
are presented as follows.' forx(Z) from equation (2'.53),. we have
State equation
-xfz) = '[Zr-$(T)J-1Zx(O)+[ZI-$(T)J-1e(T)U(z) ... (2.55)
x(k+l) = <1>(1) x(k) ,+ 8(1) UCk) , ...(2.34)
The inverse Z - transform of the last equation is
State Transition Equation
1
N-l x(KT) = Z-1[(ZI-HT)f Z]x(O)+Z-l {(ZI-HT)tEl(T)U(z)}
x(N) <I> (N)'x(0) + L q> [(N- i .; 1)] e(1) U (i) ... (2.49).
......(2.56)
i=O
WOe shall show that inverse Z -transform of [Zl - <f>(T)]-l Z is the
State equatdon
discrete state transition matrix q, (KT).
x(K+l) = A x(K) + B U(K) ... (2.50)
The Z _. transform of q,(KT) is defined in the usual manner as
State Transitio~ Equation
00
Thus,
IExample 2.1 I
... (2.59)
In this example we shall illustrate the analysis of an open-loop
Taking invers Z transform on both sides. discrete-data system by the state variable method presented' above. The
block diagram of the sytem under consideration is shown in 'Figure 2.3 .
... (2.60) The dynamic equation that describe the linear process are.
x(KT) = cp
'x-i
(kT).X(O) + LcP[(K-i-l)T]e(T)U(iT) ... (2.62)
A = [~2' _~] B = [~] ... (2.67) .
~ -t -e -2t
ze e-t- e-2t] , ,[ 2e-NT _e-2NT e -NT -e -2NT]
~(t) = L- l [SI.-Al-l = ,[ -2e-t+ 2,
e-3t ,-t '2 -2t
-e+e .. (2.70) where ~(NT) = -2e-N'r + 2e- 3NT _e- NT + 2e- 2NT '
(2.74)
-
-
T
f
-,
0 -
'2e:-(T-t) - e-2(T-t)
'2 -(T-t)
e
2 -3(T-t)
+ e
-(T-t)
:-: ~(T-<) ::e -2(T-<) 1 dt
.. [ X2(NT) r
L
_2e~NT + 2e-2NT 2e-2NT
, e-(T-t)_ e-2(T,- t) ]
. , ,dt (1- e-T)e-(N-K-.l)T -0-.5(1-. e-2T)e -2(N-K-l)Tj
=
R
o -e
-(T-t) 2 -2(T-t)
'+ e
... (2.71)
+~
N-l[
LJ
K=O -
(1 -e-T)'e -(N-K-l)T (1 -e-2T) .e -2(N-K-l)T U(KT)
T 2T
_ [0,.5 _e- ,+ 0.'5 e- ]
_... (2.72)
-
e -T -e-2T ,2g9 METHODS OF COMPUTING STATEIRANSITION MATRIX
Now substituting equation (2.70) with t = T and equation (2.72) into It is worthwhile to emphasize the difference between the formulation
equation (2.29), the discrete state equation of the system is written; of the state transition matrix under two different conditions. For the
i.e., x[(K+1)T] = <I>(T) x K(T) + SCT) U(KT) sampled data system.
x(t) = A x (t) + B U(t) "Q(2.75)
'where
U(t) = U(KT) KT 5: t :::; (K+l) T ~ .. (2.76)
2T
o .5 - e- T +0.5 e- ] U(kT)
The discretize state equations are express . as
+ .... (2.73)
[ e- T _e-2T x[(K+l) T] = <P(T) X (KT) + aCT) U(KT) ... (2.77)
where <peT) is a state transition matrix of A with t replaced by T, that is
Using equation (2.43), the solution of equation (2.73) is, .
Equation (2.43) implies
N-l
The state transition matrix <p(T) 'is given by,
, x(NT) = ~(NT)x(O)+L~[(N -K-1)T]8(T)U(KT)
k=O
... (2.78)
[The State-Variable Technique (2.18 )
Since q.(t) is the state transition matrix of A, $(NT) is determined 2.9.-2 . TheZ-tra:n'sform' method
from $(t) by substituting NT for t. For the sampled -data system of equation (2.77), the state transition
The state equations of digital control system are expressed as matrix of <j>{T) is expressed in terms of the Z - transform in equation (2.60) .
as,
x (K+1) = Ax(K) + B U(K) (2.79)
~(NT) =. Z-l {(ZI - <j>(T)-1 Z} '... (2.85)
The state transition matrix is given as
For the digital control system of equation (2.79), the state transition
matrix of A is expressed as,
... (2.80)
~(N) = Z-1 [eZI - A)-l z] ... (2.86)
In this case <peN) is obtained by multiplying A by itself N times.
In. the last two equation the evalution of STM involves the matrix
We shall first present. two methods of computing ~(NT) giveser'I') in inverse and then inverse Z-transform. For second-order systems, these
the sample data case, or <j>(NT) given A in. the digital. case. analytical steps can be carried out with ease.
However for higher order systems, the amount of work involved can
2.9.1 The Cayley-Hamilton T~~oremMethod be quite tedious if carried out by hands.
Given Matrix $(T) or A, as the case maybe equation (2.78) and (2.80) The task of finding {ZI-A)-l Z in equation (2.86),- given A, can be
can be computed using the Cayley-Hamilton theorem. simplified by the following method. Let
The theorem states that every square matrix must satisfy its own
F =' (ZI - Ar' Z ... (2.87)
characteristic equation. For example, consider that the characteristic
equation. of A (n x n) is written as, Pre multiplying both sides of the last equation by (ZI _. A) and
. rearranging we get
... (2.81)
ZF = ZI +AF "~ (2.88)
then,
Premultiplying both sides of last equation by ZI+A and rearranging
A D +a n-1.An-1 +a n -:-2 A n - Z + +a1 A + a.-o I ... (2.82) we g e t Z2F = A2 F + ZA.+ Z2 I ..,.(2.89)
By repeatedly multiplying by ZI + A,
... (2.83)
we' have
Similarly, if the stateequation is given in the formof equatidll (2.77),
A 3, F + Z A 2 + Z2 A + 'Z3 I ... (2.90)
the state transition matrixotn'I') must satisfy the characteristic equation
2
of $(T) -(nxn) that is, A 4F + ZA3 + Z2 A +'Z3 A + Z4I ... (2.91)
. r
n 'n
IZI-A!
Izj Iai Ai-j
j=i i=j = a3 Z
31
+ (a2 I+a3 A)Z2 + (al 1+ a2 A+ a3 A 2)z
F = =
IZI-AI
... (2.95)
jZI-A!
Substituting the coefficient of a3, 82, a1 and A into the last equation, '
we have
The numerator of the last equation is also known as the matrix.
[adj (ZI - A)] Z
Once F is evaluated from equation (2.95), <P(N) is given by,
<p(N) = Z-l (F) ... (2.96) F ... (2.96)
(Z-1)(Z-2) (Z-3)
For the case.in equation (2.85), simply replace A by <p(T) in equation (2.95).
Performing partial -:- fraction expansion Equation (2.96) gives
( 2.22 ) Computer Centrolof Processes]
(TheState-Variable Technique . [ 2.21 )
The Laplace Transform Method
r2(ZZ~ 1)[ ~: :]
7 8
7 11 J z [27 16
---'-6
Z-2
Repeating equation (2.14)
f(t) = L-1 I(SI-A)-l] ... (2.99)
-17 1 . -12 -28
:]
F = ... (2.97) The matrix inverse of SI-A can be conducted using the same procedure
l +2(Z~3)-6
[38
-18
equations we get
j=l i=j
(SI-A)-l = ... (2.100)
~ ]-[~ }~O.693K
7 8
! [ -6
18
7 16 : where n is, the dimension of A, and
2 -6 -17 1 -12 -28 lSI - A I = sn+ a n- 1 sn-l + + ajS + ao ... (2.101)
ep(K) = 9 -1
... (2.98)
Direct Power Series Expansion Method
+~
[ 38
-6 27 3 ]e-
1 1K
.
. The power-series representation of 4>(T) is
-18 -39 9
A2T 2
We can check this by verifying that q>(O) = I, the identity matrix. $(t) = eAT. = I + AT + 2! + ... (2.102)
Applying the Cayley - Hamilton theorem to the last equation, we have, .' IExample 2.4
To illustrate the eigenvalue method described above, consider the
... (2.105)
. A . A= [3 2] .
matrix z gIven In as
2 .3
Where n is the dimension of A and cq are constants that depend on
the coefficients of the characteristic equation of A and can be computed The eigen values of A are 81 = 5 and 8 2 = 1 which are distinct from
recursively-in general N andn are not related, so N 'can be greater than n. equation (2.108).
Ck A k ... {2~104)
[ 0.5 0.5 J
k=l
! A-SjI
... (2.114)
-3
Fig. 2.4
. z-
'x2(kT) 4 Storage ~ x1 (kT} - x2(kT) 0----.----0 x1(z) Figure 2.4 A state diagram of the digital' system described by
(equation 2.122).
(e)
The output equation is,
Figure 2.3 Basic Elements of a digital state diagram and. the
corresponding digital computer operations.
... (2.123) .
The state' transition equation of the system, which are the solutions
of the state equation in Equation (2.212) are obtained from the state
(The State-Variable Technique Computer Controlof Processes )
( 2.27 )
diagram using the gain form when ,Xl(Z) and' X2(Z) are regarded as output 2.11.1' Direct decomposition
nodes and R(z), Xl(O), X2(O) are input nodes, we have, Consider that the transfer function of a digital controller or process
is,
The state diagram of Figure (2.5) can also be used for analytical ?11.2 Cascade decomposition
purposes., If the numerator and the denominator of a transfer function D(Z)
By defining the variables at the nodes of all the time-delay units are expressed as products of first- order factors, each of these factors
state variables, the dyna~ic equations and the state transition equations can be realized by simple digital program on state diagram. The entire
can all be determined from the state diagram by applying gain formula. D(Z) is then represented by the cascade connection of the individual
programs.
Consider that the transfer functionof a digital process is written as.
*0_ (2.133)
C(z)
R(z)n---~rq.---""""--4;C}--""""'---I----""~--------I"'-~
Where K, Zi, i = 1, 2, ..... m and Pj, j = 1, 2, .... n are real constant.'
.D(Z) to be physically realizable n ~' m
Writting n(Z) as a product of the. constant Kand.n first order transfer"
funtion yields.
DeZ) = KD l (Z) D2 (Z).~ .... Dn(Z) ... (2.134)
where,
. Figure 2.5state' diagram representation of the transfer function of
Z+Zj 1+ Zi Z~~
equation (2.127) with n = m = 3'by direct decomposition. Di(Z) = L = J,2,.....m ... (2.135)
(Z+Pi ) l+Pi Z-l
The discrete state equations are obtained from Figure 2.5'by applying
the gain formula with Xl (Ks-L), X2 (Ks-L) Xg (Ks-L) as outputs. After deleting 1 Z-l
the time - delay branches with gain Z-l, the inputs are r'(K), xi(K), x2(K) DiZ) = - - j =m+l n (2.136)
Z+P j l+P Z-l
j
and Xg (K) ~
The state diagram representation of equation (2.135) is ,constructed
The state equations' are,
in Figure 2.6 (a) and that of equation (2.136) .is in Figure 2.6 (b), all by
Xl (Ks-L) = X2 (K) direct decomposition.
'X2 (K+l) = Xg (K) The overall state diagram for D(Z) is obtained by connecting the m
diagrams of Figure 2.6 (a) and n - m diagram of Figure 2.6 (b) with gain
X:l (K+1) = -a2 Xl (K) - al X2 (K)- (K)+~r(K)
ao Xs ... (2.182) branches of K, all in series.
ag' ag ag' ag .
From ,these state. equations, we can see clearly that direct The advantages with the cascade decomposition is that gain K and
decomposition will always yield a state variable model in the phase. - the constant Zi and P j are all preserved in the state diagram. In design
variable canonical form. This is apparently one advantage of using direct situation these parameters may be varied individually, and the effects
can be. observed independently.
decomposition.
(TheState-Variable Techniq ue [ 2.31] (2.32 ) computer Control of Processes] _
.CoD,1llon~nt.ofD(Z)~re all the form of figure 2.6 (b) and. the over all
state diagram is 'obtained by connecting thesein parallel. Fig. 2.6': A sampled - data control system
(The State-Variable Technique [ 2.33,J ( 2.34 ]
... (2.144)
described by the following dynamic equations: 'i=O
E 0 0 0 0 U(K)
2.15.2 Theorem on Observability
DB E 0 0 0 U(K.+,1)
The linear digital system described by equation (2.141) and (2.142)
+ DAB DB E 0 0 U(K+2)
is completly observable if and only if the following n x PN .matrix is .of ... (2.159)
rank n.
D A N - 2B D A N- 3 R' D A N - 4B DB E U(K +N -1)
... (2.156)
Where n is 'the dimension of x(K) and p is the dimension of 9(K). To determine x(K) from the last equation, given C(K + j) and U(ks-j)
The matrxL is known as observability matrix. for j = 0, 1, 2... :...N- 1, the pN x n
matrix.
Proof:
Following e9-uation '(5.64), the state transition equation of equation '
(2.141) is written as,
... (2.160)
j-l
Aj x(K) +'LAj-i-l BU(K + i]: ... (2.157)
i=O
for J == 1, 2 N- 1 must have n independent rows. (or assuming that pN ~n). The condition
Substituting Equation (2. 157) into output equation, equation (2.142), in equation (2.160) is equivalent to the matrix L in equation (2.15) must
we get have a rank n.
. J-l For a discrete - data system that is described by the state equation
C(K+j) = DAjx(K)+ LDAj-i-1BU(K+i)+EU(K+j) ... (2.158) in equation (2.147) and the out equations of equation (2.142) the condition
i::::O of complete .observabifity is that the matrix.
must be of rank n.
( 2.42 ] Cort!pu:ter,Control'.ofP,rocesses )
( The State-Variable Technique ( 2.41 )
In this chapter we shallgive definition of' stability and, investgate Definition 2.6: Bounded - Input - Bounded - StateSt'ability
the methods of testing the 'stability of linear discrete. - data and digital Consider a linear time - Invariant discrete - data system, that is
control systems. described by the following dynamic equations.
The state controlability matrix is x (K+1) =A x,(K) + B U(K) - ... (2.169)
... (2.166) C (K) = D x (K) + E U(K) ... (2.170)
The system is said to be bounded Input - bounded - state (BIBS) stable
if for any bounded input; the, state x (K) is also bounded.
[~x~ (K)J
Since L is.singlular, the system is unobservable.
Ilx (K)II = - ... (2.171)
We observe that all the coefficients are of the same sign, and none of
coefficients is zero, which are the necessary and sufficient conditions for
all the roots of the second - order equation in equation (2.179) to be in the
left .half r - plane. Thus, all the roots of equation (2.178) are insidethe
unit circle.
I Example,2.7 I
Consider that the characteristic equation of a third - order discrete
- data control system is given as
@ F(Z) = Z3- 1.25 Z2- 1.375 Z - 0.25 = 0 ... (2.180)
Z-Plane y .. Plane, Which is known to. have roots at Z = 0.25, 0.5 and -2.0
Figure 2.11 Mapping of the unit circle I Z "=1 onto the imaginary axis . Transforming F(Z) into the r - domain using Equa~ion.(2~180), we get
of the r - plane through, the r - transformation Z = (r + 1) '/, r-1 '-1.875 r + 3.87 r
3 2
+ 4.875 r +1.125 = 0 ... (2.181)
I Example 2.6 I Since the coefficient of the last equation do not all have same sign,
we know that the system has at least one zero outside the unit circle in
Consider that the matrix A of the state equations in equation (2.169) is the- Z - plane. Nevertheless applying the Routh - Hurwitz Criterion to
equation (2.181), we have the following Routh tabulation.
'[-0.5 0]
A = ,0 0.5 ... (2.176) r3 - 1.875 4.875
Sign change
- The characteristic equation of the system is
r2 3.875 .' 1.125
IZI - A I = (Z - 0.5) (Z + 0.5) =0 ... (2.177) r 1
5.419 o
F(Z)= Z2-0.25=O ... (2.178) 1.125
Transforming F (Z) into' a polynomial in r using Since there is one sign change in first column of the Routh tabulation,
equation" (2.181) has one root' in the right-half of r - plane or 'equation
r+1 (2.180) has one root outside the unit circle in the Z - plane, which in this
Z = --'
r-1
we get
case is known to be at Z' = -2.
, Singular Cases
(r+lr -0.25
r-1
= 0 For the first singular case in ther -plane, the same remedy is applied
by replacing the zero element in the first row by ail arbitarily small
(r + 1)2'_ 0.25 (r - 1)2 = 0 number E.
0.75z-2-2.5r + 0.75 = 0 ... (2.179) When a row of zeros occur in the Routh tabulation of the equation'
.looking from the transformation relationship between Sand ,Zand r we
conclude that
.
the, following conditions may exist.
'
(TheState-Variable Technique ( 2.48 ) ComputerControl of Process!!]
1. Pairs of real roots in the Z - plane that are the inverse of each Taking the derivative of ~(r) with respect to r we get,
other; eg Zl = 2 and Z2 = 1/Z l = 0.5.
dA(r)
2. Pairs of roots on the unit circle 1Z,I= 1, including Z = 1 and --'-
dr
= I.Sr ... (2.185)
Z = -1 simultaneously.
The second condition corresponds to roots on the imaginary in the r plane. dA(r)
.
. . . The coefficient of the r 1 row are now filled with the coefficient .~
The' remedy to the singular case is simila~ tothat for the S -.domain.
and the remaining Routh tabulation is,
1. Form the auxilary equation A (r) = 0 using the coefficient of the r1 1.8 0
row in the Routh tabulation just above the row of Zeros.
rO - 0.1
.. ,2. Performd A (r)/dr = 0
Since one sign change occurs in the first column of Routh tabulation
3. Replace the row of zeros in the tabulation with the coefficien~s of the characteristic equation has one root in the right half of the r - plane,
d A(r) I dr. or Ffz) has one root outside the unit circle in the Z - plane. Solving forthe
4. Continue with the Routh tabulation in theusual manner. roots of the auxiliary equation A(r) = 0, we have,
r = 0.333 and - 0.333 ... (2.186)
IExample 2.8 ~.
Which are two of the roots of the characteristic equation. These two
Consider the equation roots correspond to, the root at Z =- 0.. 2 and - 0.5 respectively.
F(Z) = Z3 + 3.3 Z2 + 3Z + 0.8=.0 . ... (2.182)
Which has roots at Z =- 0.5, - 0.8 and - 2.0. The roots at ., 0.5 and
IExample 2.9 I
- 2.0 form an inverse symmetry about Z = 1, and should create a singular Consider adiscrete -data system that has the following equations.
case for the Routh test. F(Z) = Z3 + Z2 + Z + 1 = 0 .~- (2.187)
I.
(The State-Variable technique ( 2.51 ) ,(2.52 ) Computer Control of Process:!!J
One advantage of the Jury stability test is that, for systems 'of any Since the equation. is of the second order, under Jury's test the
order the eonditions on F(1), F(-l) and between ao and iriin conditions in Equation (2.192) are necessary and sufficient for the system
equation (2.191) from necessary conditions of-stability that are very simple to be stable.
to check without carrying out Jury tabulation. For second order .systems, Thus F(l) > 0 F(1) = 1- 0.25= 0.75> 0
the conditions in Equation (2.192)-are necessary and sufficient.
F(-I 0 F(-I)= 1- 0.25 = 0.75 > 0
For equations that pass the necessary conditions and would require
carrying out Jury's tabulation, the inequality tests of the calculated laol> (a2) laol = 0.25 a2 = 1
coefficients in Equation (2.191) could be quite tedius if the equation Thus, laol < a2 The system is stable expected.
contains variable parameters.
Singular Cases IExample 2.11 I
When some or all of the elements of a row in the Jury tabulation are Consider the characteristic equation, in Equation (2.180)
zero', the tabulation ends prematurely. We refer to this situation as a F(z) = Z3 -1.25Z2 -1.375Z -''0".25= 0 ... (2.195)
singular case. A singular case can be remedised by expanding or
contracting the unit circle I z I = 1 infinitesimally, which is equivalent to ag = 1 a2 = -1.25 al;: -1.375 8.0 = -0.25
moving t'hezeros of Ffz) of the unit- circle. The transformation for this Checking the first three conditions in Equation (2~ 191)
purpose is
F(1) > 0
Z = (1 + s) Z . ... (2.193)
.Ft L) 1 - 1.25 - 1.375 - 0.25 = ~1.875 ... (2.196)
where e is a very small real number. When is a positive number in
Equation (2.193), the radius of the unit circle is reduce to (1 + E), and F(-l)< 0
when. e isnegative, the radius of the unit circle is reduced to 1 + e. The F(-1)= -1- 1.25 + 1.375 - 0.25 = -1.125 ... (2.197)
difference. between the number of zeros found inside (or outside) the unit
circle when the circle is expanded or contracted by E is the number of laol < a3
zeros on the circle. I ao I 0.25 < 1 ... (2~198)
The transformation of Equation (2.193) is actually quite'simple to Since F(l) is negative, not all the roots of Equation (2.195) are not
apply, since inside the unit circle and system is instable. There is no need to carry out
(1 + s)" zn == (1 -i- ns) zn ... (2.194) the Jury tabulation.
For positive or negative E. This means that the coefficient of the zn
terms i~ multiplied by (1 + fie). IExample 2. 12 1
Consider the following characteristic equation of a discrete-data
1 Example 2.10 I system.
Consider the characteristic equation given in Equation (2.178) Z3 + 3.3 Z2 + 4Z+ 0.8=0 ... (2.199)
F(z) = Z2 ~ 0.. 29 = (2 - 0.5) (2 + 0.5) = 0 The equation has roots at Z = u.2463, Z = -1.5268 + j 0.9574 and
With reference to Equation (2.189),'the coefficient of the equations are' Z =-1.5268 - jO.9574. Thus the system is unstable.
. RZ = 1, at = 0, . ao = -0.25
( 2~54 ) ComputerControlof Processes )
(The State-Variable Technique
The Liapunov Function
Applying the necessa~y condit.ions of the Jury test in Equation'
(2.191), we have Any function V = V[x(k)] of definite sign (positive definite "or
negative definite is called a Liapunov function. If V( 0) = O. and x(k) is the
F(l) > 0 F(l) = 9.1 ... (2.200) solution of the state equation.
F(-l)< 0 F(-l)= -0.9 ....(2.201) x(k +1) = f(x(k .... (2.206)
laol < a3 .1 ao I =0.8 < 1 ... (2'.202) We define a difference opeartion 8 so that 8V [x(k)] is the solution of
Thus, we must carry out the Jury tabulation as follows to determine the state equation.
the stability: of the system x(k + 1) = f[x(k)] ... (2.207)
Row ZO Zl Z2 Z3
We define difference operation ~ so that ~v [x(k)] is defined. as
1 0.8 4.0 3.3 1.0
LiV [x(k)] =. V[x(k +1)] - V[x(k)]
2 1.0 3.3 4.0 0.8
Stability Theorem of Liapunov
3 .bo bI b2 , Consider a discrete-data system described by
h2 = [ao a1=aoa2-a1a3=
] . -1.36 ... (2.205)
such that
a3 a2' .
' . 1. V[x(k)]= O} = V(O) == 0
Since I bo I is not greater than Ib 2J Equation (2.199) has at least one 2. V[x(k)] = Vex) > 0 for. x "* 0 ... (2.210)
root out side the unit circle.
3. vex) approaches infinity as II x II 4 00
Vex) = X1 2 (k) + X22 (k ) ... (2.213)- L\V(x) = - XI 2 (k) - X22(k ) ... (2.20)
Which is positive for all values of XI(k) and A2(k) not equal to zero. So that it is negative definite for all Xl (k) and x2(k) *- 0
Thenthe function ~V(x)is evaluated, using Equation (2.208) FormingAv'(x) according to Equation (2.208)'
L\Vex) V[x(k + 1)] - V[x(k)] We have
= V[x(k + 1)] ~ V[xCk)] LlVex) = V[x(k + 1)] - Vlxtk)] ... (2.221)
= x~2 (k + 1) + x2 2 (k + 1) -'X12 (k ) - X2 2 (k) .. .{2~214) = 1.25al X1 2
(k) - 0.5 a2 xICk) x2(k) - 0.75 ag ,x2
2
(k ) ... (2.222)
Bubstitutefor x1(k +1) and x2(k + 1) in equation (2.214) Comparing Equations (2~222) and (2.220)
= 0.25 X1
2 (k ),+ 0.25 2 2 k) - X22(k )
X2 (k ) - XI C ag = 1.333 ... (2.2'23)
= -0.75 X1 2 (k) - O. 75 X22(k)~ ... (2.215) - Thus from Equation (2.219) .
Since ~V(x) is negative for all x(k) *- 0, the. system is asymptotically Vex) = -O~8' X1 2 (k) + 1.33 X2
2
(k ) (2.223)
stable. ','
Since Vex) is indefinite there is again no conclusion of the stability
I Example 2. 14 I condition.
Consider the digital system described by the state equations Liapunov Stability Theorem for Linear Digital Systems
xI(k + 1) = -1.5 xI(k) ... (2.16) The stability and instability theorem of Liapunov are valid for both
linear and non linear systems. The last two examples iflustrate that a
X2 (k + 1) = -0.5 x2(k) ... (2.17)
successful execution of the Iiapunov tests depends upon 'the correct guess
It is simple to show that the eigen values of the A matrix are at -1.5 or assignment of Vex) or LiVex), whichis always a' difficult task. However,
and -0.5 and that system is unstable. However without prior knowledge for linear digital systems a simple procedure is available.
of the stability conditions of the system, the stability theorem of Liapunov
is applied.
(TheState-Variable Technique' Computer Control of Processes J
(2.57 )
Consider that a linear time - invariant digital system is described IExample 2.15 f
by the difference equation
Consider the discrete-data system discribed by
x(k + 1) = Ax(k) ... (2.224) xICk + 1) = -0.5 xI(k) ... (2.230)
where x(k) is n x 1 and A is an n x n matrix. The equilibrium state x2(k+ 1) = -0.5 x2(k) ... (2.231)
Xc = 0 is asymptotically stable, if and only if, given any positive definite
The coefficient matrix is .
real symmetric matrix Q, there exists a positive definite real symmetric
matrix P such that
... (2.225) .
A= [-0.50 -0.50]' .~. (2.232)
[~~]
... (2.226)
Q =. ... (2.333)
is a liapunov function for the system,and further
and let P be of the form
dV(X) = -xT(k) Qx(k) (2.227)
.and the characteristics equation of the closed - loop system is Fig. 2.13: Control configuratiori' for a state regulator
ISI-(A-Bk)] = 0 ... (2.340) In this structure, these is no command input (r = 0). Control systems
When evaluated, this yields' an nth - order polynomial in s containing in which the output mustfollow the command signals (called-servo system)'
the n gains k}, k 2 , ... k n The control . . law design then consists' of picking. will be considered later.
the gains. So that the roots of Equation (2.340) are in desirable location. Selection of desirable locations for the closed-loop poles requires
In the 'next section we find that a mildly restrictive condition some iteration by the designer. For now, we will assume that the desired
(namely, the system (2.337) must "Qe completely controlable),all the eigen locations are known, say,
values of (A - B]{) can be arbitrarily-located in the complex plane, by S = AI, A2' ... An
choosing k suitably (with the restriction that complexeigen values, occur
Then the desired characteristics equation is
in complex . . conjugate pairs).
... (2.341)
If all the eigen values of (A - Bk) are placed in the left-half plane,
the closed- loop system is, ofcourse, asymtotically stable; x.(t) 'will decay <The required elements of k are obtained by matching coefficients in
zero irrespective of the value ofx(O) the initial perturbation in the state'. Equations (2.340) and (2.341) then forcing the system characteristic
The system. state is thus maintained at zero value in spite of distrubance equation to be identical with the desired equation.
that act upon the system. The canonical form used in control-law design is the controllable
Systems with this property are called 'regulator systems'. The origin canonical form. Consider, a system represented by the transfer function
of state space is the 'equilibrium state' of the system.
Y (8) ~I Sn-l + f32sn~2 + ... + f3 n
Control configuration, for, a state regulator is shown in Figure 2.13. U(S) Sn+XlSri-l+ .. '+Un
(The State-Variable Technique [ 2.63,) ( 2.64 ) Computer Control of Processes )
A companion-form realization of this. transfer function is given below 2.18.3 Necessary and Sufficient conditions for Arbitrary Pole-
x = Ax. + Bu Placement
y = ex ..(2-.342) Consider the linear time-invariant system (2..3 37) with, state-feedback
control law (2.338); the resulting closed-loop system is given by
Equation (2.339)~ In the following, we shall prove that necessary and
~1
0 1 0 0
sufficient condition for arbitrary placement of closed-loop eigen values
0 0 1 0
in the complex plane (with the. restriction that complex eigen values occur
.where A b
in complex-conjugate pairs) is that the system (2.337) is completely
,0 1
~J
0 0 controllable. We shall first prove the sufficient condition, i.e., if the system
:-a.n -U n-1 -cx,n-2 -0.1 (2.337) is completely controllable, all the eign values of (A - Bk) in Equation
(2.339) can be arbitrarily placed.
C = [~n Bn - 1 ... B2 , B l ] In proving the sufficient condition on arbitray pole-placement, it is
The matrix A in Equation (2.342) has a very special structure. The convenient to transform the state equation (2.337) in to. be controllable
coefficients of the denominator of the transfer function ,preceded by minus canonical form (2.342).
signs form a string along the bottom row of the ~atrix.The rest of the Let us assume that such a transformation exists and is given by,
matrix is zero except for the superdiagonal terms, which are all unity. It
can easily be proved that the pair (A, B) is completely controllable forall x = Px
values of aiS. For this reason, the companion form realization given by
equation (2.342) is referred to as the controllable canonical form.
PI 1 P12 Pin 1
:~j
Pi1 P22
One of the advantages of the controllable canonical formis that the
controller gains can be" obtained from it, just by inspection. The closed- Pn l Pn2
loop system matrix.
0 1 0 0
Pi = [~i1. Pi2 ... Pin] := 1~ 2, .... n.
(2.337)
0 0 1 0
A-Bk = Under the transformation (2.343), system (2.337) is transformed to
the. following controllable canocical model:
0 0 0 1
-an -k1 . ~a1, -kn
-cx,n-l -k2 -a n - 2 - k a x = Ax+BU ... (2.344)
~I
0 1 0
s" + (a~ +kn ) sn-1 + ... +(an~2-+ kg) '8 2 + (cx,n-l + k 2) s + Un + k~ - 0. 0 0 -I
and the controller gains can be found by comparing the coefficients where, K= PAP-l =
-~J
of this characteristic equation with equation (2~341). 0 0 0
-an -Un-l -U n - 2
(The State-Variable Technique, [ 2.65 ) ( 2.66 ) Computer Control ofProces$es,J
B= PB =
o 0
0
r 'PlAB
PIB
=
lPIA~_2B
'0 PB
1 0
1 PlAn-lB
I~I - AI = s" + n l
U1S - + ... + an-l s + an = 'I SI - A I
Or
. (characteristic polynomial is invariant under equivaIerice
Pl~B n 2B' An-lbB]
transformation), l. AB ... A - [ 0 ... 01]
=0
The first equation in the. set (2.343) is given by
This gives
Xl = PII Xl -f:'P12 X2 + ....+ Pin Xn = :PIX
Taking the derivative on' both sides of this equation, we get where
=. PI X .=
,~
xl PlAx +.PIBu is the controllability matr-ix, which is non s ingular because of the
But xl (= x2) is a function of X only as per the canonical model (2.344) , assumption of controllability of the system (2.337).
r PI u = -kx kx .. ~ (2.346)
x= Px ,= l .
PIA X where k-- k p~l= [k1 k2 .. k n ]
PI~n-1
...
~
It is chosen, an The eigen values of'{A - Bk) can be arbitrarily placed in the complex
plane (with the restriction that complex eigen ~alues' occur in complex-
conjugate pairs) by choosingk suitably if, and only if, the system (2.350)
Transforming the feedback controller (2.346) to the original is completely controllable. The important result on pole placement was
proved in the previous sectiou. The following design steps' for pole
coordinates, w'e obtain
placement emerge from the proof, '
.... (2.349)
This proves that -if (2.337) is controllable, the closed loop poles can
be orbitrarily assigned (sufficient condition). .
(TheState-Variable Technique ( 2.69] ( 2.70 ) Computer Control of Processes)
,Step 2:
Determine the transformation .matrix P that transforms the system 1
= [0 0 .. , 0 1]uc- [(a l -u1)A
n-1 +(a2 -u n-2 + ... +(a -un)]
2)A n
(2.350) into controllable canonical form.
(2.357)
The characteristic polynomial of matrix A is equation (2.353)
p I. SI - A I = s
n
+, alS rr-L + ans n-2 + .... + an-l S + Ctn
and determine the vaues of ai, a2, ... , an-I, an. Equations (2.359 (a describe the Ackermann's formula for the
Step 4: determination of the state-feedback gain matrix.
= (A-mc)x(t) ... (2.366) A companion-form realization of this transfer function is given below
It can easily be proved that the pair (A, C) is completely observable Since
for all values ofcq's. For this reason, the conipanion form realization given
det W
by Equation (~.a68) is referred to as observable canonical form.
One obtains
One of the advantages of hte observable canonical form is that the.
observer gainsmcan be obtained from it, just by inspection. The observer- , det lSI - (AT - CT mT)] = det lSI - (A - mC)]
error matrix is (i.e.,) the eigen values of (AT - C T mT) are same as the eigen values of
(A- me)
0 0 0 -an - f i l
Table 2.1 Duality between and control and Estimation
1 0 0 -an-l -m2
A AT
0 0 1 ~al ~mn
B CT
which has the characteristic' equation
k mT
sn + (al + mn)sn-l + ... + (an- 2+ ms)s2 + (Un-l + m2) s + an + mj = 0 and
the observer gains can be found by comparing t.he coefficients of this By comparing the characteristic equation of the closed-loop system
equation' with equation (2.367 (b. (2.352) and that of the auxilIary system (2.370), we obtain the duality
A procedure for observer design, therefore, consists of transforming relations given in Table 2.1 between the control and estimation problem. '
the given state variable. model to observable canonical form, solving for The Ackermann's control-design formula given by Equations (2.359)
the observer gains, and transforming the gains back to the original state. becomes the observer-design formula if the substitutions of Table 2.1 are,
made.
V!e can, howevervdirectly use' the equations of the control-law design
for computing the observer gain matrix ID, if we examine the resemblance A necessary and sufficient condition for determie of the observer
between the estimation and control problems. In fact, the two problems gain matrix ill for the desired eigen values of (A - me), is that the auxillary
are mathematically equivalent. This property is called duality. The design system (2.369) be completley controllable. The controllability condition
of a full.. order observer requires the determination of the gain matrix 'm' T T 1 T
for this system is that the rank of [C A TC ... (A Tr- C ] is n. This
such that (A _. me) has desired eigen values Ai, i = 1, 2, ... n. This is
mathematicaly equivalent to .designing a full jstate.. feedback controller is the condition for complete observability of the original system defined
for the 'transposed auxillary system'. by Equation (2.362). This means that a ne~essary and sufficient condition
for estimation of the state of the system defined by Equations (2.362) is
~(t) = AT l;(t) +CTll(t) ...(2.369a)
that the system be completely observable.
With feedback
Again by duality, we can say that for the can of single-output system,
l1(t) = -inT ~(t) ... (2~369b) the gain matrix m, which places the observer poles at desired locations,
So that the closed-loop. auxiliary system. is unique.. In the multi-output case, the same pole configuration can be
achieved by various feedback gain matrix.
. ~(t)= (AT _CTmT)~(t) ... (2.370)
has eigen values A.i; i = 1, 2, ... n
[ 2.78 ) Computer Control of Processes )
(Th,e State-Variable 1echnique ( 2.77 )
where the two rightmost terms are known and can be considered as
2.18.7 Reduced .. Order StateObserver an input into the Xe dynamics. Since Xl = y, measured dynamics are given
The observer developed' in the previous. subsection reconstruct the entire by the scalar equation
state vector. However, the measurements usually available are some of the states
... (2.373)
of the plant. For example, for the satellite-attitude control problem considered
in the previous subsection, the measurement is orientation of the satellite, which If we collect the known terms of Equation (2.373) on one side, we get
is Xl(t). The measurement of a state, in general will be more accurate than any
y' -- -allY - hI U = aie x, ... (2.374),
estimate of the 'state based on the measurement. Hence, it is not logicalin most -v--'
Known Measurement,
cases to estimate states that we are measurement.' One possible exception is
Note that Equations (2.372) (2.374) have the same relationship to
the case in which a measurement is very noisy. The state observer for this case
the state X e that the original Equation (2.371) had to the entire state x,
may furnish some beneficial noise fillter,
Following this line of reasoning, we can establish the following'
Since we will not usually want to estimate any state that we are measuring, substitutions in the original observer design equations, to obtain an
we prefer' to design an observer that estimates only those states that are not (reduced order) observer of X e :
measured. This type of observer is called a "reduced-order state observer".
x -(- x,
We develop design equations for such an observer in this subsection. A +- A ee
We consider only the case of one measurement. It is assumed that the
bu~ aai Y + beu
state variables are always chosen such that the state measured is ~l(t);
we can, do this without loss of generality. The output equation is given by Y <f- Y- allY - b l u
yCt) = Xl(t) = Cxtt}. C ~ aie ... (2.375)
where C = [1 0 0 ... 0] Making these substitutions into the equation for full-order observer
(Equation 2.364) we obtain the equation of the reduced -order observer.
To derive the reduced - order observer, we partition the state vector
into two parts; one part is Xl which is directly measured and the other
part is Xe, representing the state variables that need to be estimated. i e= Aeex e + "aeiY + bell + m[y-allY-blU-aiexe)
-------v-----' ~--'
... (2.376)
input . measurement
[X~le J = [al~
a
aie J[X~J + [blJu
A ee x,
el be .
...(2.371a) ... (2.378)
observer are same as for the full-order observer - namely observability of Where x is the nx 1 state vector, u is the scalar input, y is the scalar
the pair (A, C). output; F g and c are, respectively, n x n, n x 1 and 1 x nreal constant
matrices and k = 0, 1, 2_ ~ ....
Let us now look at theimplementational aspects of the reduced-order
observer given by Equation (2.376). This equationcan be rewritten as We will carry out the design of digital control system for-the plant
(2.383) in two steps. One step assumes that we have all the elements of
xe = (A ee -male)xe~(ael -mall)y+(b e -mb1) u + mY .. (2.381)
the state vector at our disposal for feedback purposes. The next step is to
The fact that the reduced-order observer requires the derivative of design a state observer which estimates the entire state vector, when
yet) as aninput, appears topresenta practical difficulty. It is known that provided with the measurement of the 'system indicated by the output
diffenentiation amplifier noise. So if y is noisy the use of y is unacceptable. equation in (2.383).
To get around this difficulty, we define the new state as, The final step will consist of combining the control law and the
... (2.382a) observer, where the control law calculations are based on the estimated
state variables rather than the actual state.
Then, in terms of this new state, the impelementation of the reduced
order observer is given by
2.18.9 State Regulator Design
... (~.382b) - Consider the nth - order, single-input, linear time- invariant system
and if no longer appears directly. A block diagram representation of the (2.383) with state feedback control law .
reduced order observer is shown in Figure 2.16. u(k) = -kx(k) ... (.2~384)
y------.--,..---------.. .
The resulting closed-loop system is
x(k + 1) = (F- gk) x (k) ... (2~385)
If all the eigen values of (F - gk} are placed inside the unit circle in
the complex plane, the state x(k) will decay to the equilibrium state x = 0
u
irrespective of the value of x(O) the initial perturbation. in the state.
n-1 parallel
integrators A necessary and sufficient condition f~arbitray placement of closed-
loop eigen values (with the restriction thatcomplex eigen values occur in
conjugate pairs), is that the system (2.383) is completely controllablet
Dc = [g Fg ... Fn~lg ]
where F = eAT = [~ ~]
... (2.388)
g
( 2.84) Computer Control of Processes J
(The,State..v ariable Technique "f 2.83 )
2.18~11 Prediction Observer
The characteristic equation of the open - loop system is
An estimation, scheme employing a full-code observer is shown in
JZI-FI
r Z - 1 _T l
Z- 1 J I=(Z _1)2' = 0
,
Figure, and the equation for it is
! 0 x(k+l) = Fx(k) +gu(k) + m(y(k)-Cx(k))'
'L ... (2.394)
with, the control law where m is an n x 1 real constant gain matrix. We will call this a
u(k) = -kx(k) = -[k1 k 2]x(k) predication observer because the estimate x(k+l) is one sampling period
ahead of the measuremet y(k).
the clos sd loop system becomes
x(k + 1; = .F(- g k) x(k) u(K) y(K)
... (2.893b)
Fig. 2.17: Prediction Observer
Matching coefficients in Equations (2.393 (a) and 2.334 Eb) we obtain
A difference equation describing the behaviour of the error is
k 1 = 10; k 2 = 3.5 obtained by subtracting Equation (2.394) from Equation (2.383).
The control law designed in the last subsection assumed, that all where x = x - X
states were available for feedback since, typically, not all states are
The characteristic 'equation of the error is given by
measured; the purpose of. this subsection is to show how to determine
algorithms which will reconstruct all the states, given measurements or IZI - (F - mel] == 0 ... (2.396a)
a portion of them. If the' state is x, then the estimate is x and the idea is Assuming that the desired- characteristic equation is
to let u':-= --k X; replacing the true states by their estimates in the control ... (2.396b)
law.
( 2.86 ) Computer Control' of Processes]
(The State-Variable Technique ( 2.85 )
The required elements ofm are obtained by matching coefficients in The improved x(k+l) is
i(k~+l). 'I'he state observer based on this
Equation (2.396 a) and Equation (2.396 b). A necessary and ~ufficient formulation is calledthe current observer. The current observer equation
condition for the arbitrary assignment of eigen values of (F - me), is that are given by
the system (2.383) is completely observable. x(k+l) = Fi(k) + gu{k] ... (2.398a)
The problem of designing a full - order, observer is mathematically
equivalent to designing a full state - feedback controller for the 'transposed x(k+1) = 'x(k+l) + m[y(k + 1)-cx(k + l)J ... (2.398b)
auxiliary 'system'.
In practice, the current observer cannot be implemented exactly
t;(k +,1) = FTt;Ck) + C T llCk) ... (2.397a) because it is impossible to sample, perform calculations, and output with
With feedback absolutely no time elapse. However, theerro.cs introduced due to
computational delays will be negligible if the computation time is quite
ll(k) = _mT l;(k) ... (2.397b)
short-compared to the sample period.
So that the closed-loop auxiliary system
The error Equation for the current observer is similar to the error
t;(k + 1) '=' (F T - C T m T) s(k) ... (2.397c) eqution for the prediction observer given in (2.395). The current-estimate
error equation is obtained by subtractly Equation (2.399) from (2~383).
has eigen values Ai, i = 1, 2 ... n.
This duality provides may be used to design full-order state ob-servers .x(k+l) = x(k+l) - x(k+1)
by Ackermann's formula (2.387) (2.390).
= Fx(k) + gu(k) - Fi(k) --, gu(k) - mc[x(k+l)-x(k)]
,2.18.12 Current Observer
FX(k) - mcFx(k)
The prediction observer given by Equation (2.394) arrives at the state,
estimatex(k) after receiving measurements up through yCk-1). Hence (F-mcF)x(k) ... (2.399)
the control u(k) = -lci.(k) doesnot utilize the information on the' current
Therefore, the gain matrix m is obtained exactly as before, except
output y(k). For higher- order system controlled with a slow computer, that C is replaced by CF.
or any time the sample rates are fast compared to the computation' time,
this delay between .making a .measurement and using it.in control law
2.18.13 Reduced - Order Observer
may be a blessing. In many system, however, the computation time
required to evalute Equation (2.394) is quite short-compared to the sample The observers discussed so far, are .designed to reconstruct the entire
period- and the control based on prediction observer may not be as state vector, given measurement of some of the states. To pursue an
accurate as it could be. observer fOl only the unmeasured states, we partition the state vector
into two parts. One part is Xl which is directly measured, and the other
An alternative formulation of the state observer is to use y(k) to
part is Xe , representing the state, variables that need to be estimated. If
obtain the state estimate x(k) . This can be alone by separating the
we partrtion the system matrices accordingly, the complex description of
estimation process into two steps, In the first step we determine x(k+l),
the system (50)i8 given by
an approximation of x(k+ 1) based on x(k) and ufk),using the model of the
plant. In the second step, we usey(k+1) to improve x(k+l). ,
rXl (k+l)ll
... (2.400a)
Lx, (k+l)-J
(The State-Variable Technique (2.88 ) Computer Control of Processes )
( 2.87 )
IExample 2.16 I
y(k) = [1 0] [Xl (k)] ... (2.400b)
xe{k) As an example of complete design, we will add a state observer to
the satellite-attitude control, The system equations of motion are (refer
The portion describing the dynamics of unmeasured states is to Equation (2.392).
The measured dynamics are given by the scalar equation y(z) _ a(Z) = C(ZI _ F)-l g+ d ... (2.408)
u(z) -
y(k + 1) - f 11y (k ) - gju'(k) = FIe x, (k) ... (2.4~2)
Applying this result to the model (2.407)~ Wa obtain
Equations (2.401) and (2.402) have the same relationship to the state
x, thatthe original equation (2.38---) had to the entire state x. Following . u(z) _ D(Z) = k(ZI - F + gk + me)"! m ... (2.408)
this reasoning, we arrive at the desired observer by making the following -y{z) -
, . substitutions into the observer equation.
x(k + 1) = Fx(k) + gu(k)
Fee xe(k)+ FeIy(k)+ geu(k) + m(y(k + 1)-fI1y (k)- glU(k)- fIe x(k)) The partitioned matrices are
\- J \ -I
. input . . measurement
... (2.404)
SubtractingEquation (2.404) from (2.401) yields the error equation
From Equation (2.406), we find the characteristic equation interms
... (2.405)
of'm:
IZ - (1 - mT) I= 0
The characteristic equation is given' by For the observer to be about four times faster than the control; we
... (2.406) place the observer pole at
We design the dynamics of this observer by selecting ill so that Z = 0.5 (= 0.835)4
Equation (2.406) matches a desired reduced-order characteristic equation. ' .. I-roT = 0.5
'rhe design may be carried out directly or by using duality principle.
( 2.90 ) Computer Control of Processes)
(The State-Variable Technique [ 2.89 )
TWO MARKS QUESTIONS AND ANSWERS
For T = 0.1 sec,'
1. What is state space analysis?
m=',5.
The state space analysis is a modern approach ;or powerful
The observer equation 'is.
technique for the analysis and design of control systems, and also
easier for analysis using digital computer.
2. Write the advantages 'or state space analysis?
.. This analysis can be carried with initial conditions and can be
= O.5X 2 (k ) + 5(y(k+l)-y(k)}+O.075u(k) carried on multiple input and multiple output systems.
.Substituting for u(k) from the control law This approach is easily amenable to solution through digital
computer.
u(k) = ~10y(k)-3.5X2(k) ... (2.409a)
The modern control are specially apparent when we design
We 'obtain controllers for systems with more than one control input or
semed output.
,X2 (k + l ) = O.2375X 2 (k) + 5y(k+l)-5.75y(k) ... (2.409b)
3. Define state.
To relate the observer based state-feedback design to a classical The state of a dynamic system is a minimal set of variables such
design, one needs to compute the Z-transform of Equations(2.409a) and "that the knowledge of the inputs for t 2:: to, comp.letely determines
(2.409b) obtaining the behaviour of the system for t > to.
9. How the state dtagram of the system. is represented? 14. What are the different methods available for computing A k?
Adder It is states that the system is static if all the poles of the Z - transfer
function of the system lies inside the unit circle in Z -. plane.
Unit delay
16. Define characteristic equation of sampled data system,
11. Define state transition matrix.
The characteristic equation is the denominator polynomial or the
The solution of the 'homogeneous state equation.
Z - .transfer function of a sampled data control system.
x(k + 1) = Ax (k) ... (1)
17~ Write the types of stability.analysis of sampled data control
as x(k) =~ (k) x (0) ... (2) systems,
where, 4> (k) is a unique n x n matrix satisfying the condition. . Jury's stability test.
<I> {k -+- 1) =A <I> (k) $ (0) = I ... (3) Bilinear transformation.
k ... (4)
<1>-. (k)= A Root locus techniques.
From equation ,(2) the solution of equation is simply a 18. What are the necessary and sufficient conditions' to be
transformation of the initial state. Therefore, the unique matrix satisfied for the stability of sampled data control system
<I> (k) is' called the state transistion matrix. It is also called the using Bilinear transformation and Jury's stability test"
fundamental matrix..
Jury's stability test
12. Write the pr-oper-tf es of "the state transistion matrix of
Necessary condition
discrete time system?
Let F(Z) be the characteristic equation. of the system F(l) >0,
<I> (0) I and (_l)n F (-1) > 0
$-1 (k) <I> (- k)
or
ep (k, lto) = <I> (k - k o) A(k...,. kO), where k > k o
F(-I 0 Even power
13. What is state transition matrix of discrete time systemf
F(-1)< 0 Odd power :
The.matrixq> (k) (or A k ) is called state transition matrix of discrete
time system.
It is given by, A k = Z-l {(ZI - A)-l Z} It is used to find the state
of the system, at any discrete time instant k.
(The.State.;VariableTechnique [2.93 ). Computer Control of Processes )
27. Write the concepts. of, obser-vabtlfty? 320, State pole - placement design.
The concept of observability is useful in solving the problem of In pole. - placement. method, the closed loop poles may be placed
reconstructing unmeasurable state variables. at any desired locations by means of state feedback through an
appropriate sta~e feedback gain matrix 'k'. The necessary and
It is the complement of controllability. sufficient condition to satisfied by the system for .arbitrary pole
It is necessary to estimate the unmeasurable state variables in placement in that system be completely state controllable.
order to contruct the feedback control signals. 33. Write the expression .for state feedback gain matrix?
28. Definition ofobservability?
K = [k 1 k 2 knl
Complete observability: The digital system describe in equation
x (k + 1) =Ax (k) + Bu (k), y(k)' ~ ex (k) +Du (k) is said to he
completely observable if for .any initial time (stage) k = 0, any
state 'x (0) can be determined from knowledge of the output C(k)
and input u (k), for 0 ~ k :::; N, where N is some finite time.
29. Write the features of observability?
It can be used to find whether the state variable is observable
or measurable.
It is very useful in solving the problem of reconstructing
unmeasurable state variables from measurable ones in the
minimum possible length of time.
In state feedback control the estimation of unmeasurable state
variables is estimated in order to construct the control signals.
30. What are the methods used to find observabifityf
Gilbert's' test.
Kalman's test.
31. State the Kalman's test for observabffity? Consider a system "
with stata model.
x (k + 1) = Ax (k). + Bu (k) and
~
3
h
INTRODUCTION TO Volume V
Concentration C
~ FlowF
SYSTEM IDENTIFICATION Concentration C
A model is an image of a dynamic system. A dynamic system is shown Suppose we, want to design a controller which acts on the flows Fl(t)
in .Fig.' 3.1. The dynamic system is driven by input variables tr(t) and and F 2(t) using the measurements of F(t) and c(t).,The purpose of the
disturbances vet). The user can control u(t) to obtain desired output yet). controller is to ensure that F(t) or C(t) remain conta.nt even if the
For a dynamic system the control action at time t will influence the output concentrations Cl(t) and C 2( t ) 'vary considerably.
at time instants s > t. For. such a controller design we need some form of mathematical
model which' describes how the input, output and disturbing variables
Disturbance
v(t) are related.
Many industrial pr,?cess, for example for production of Iron, Sugar,
Glass, Paper must be controlled in order to run safely and efficiently..
Input Output, To design controllers, some type of model of the process is needed.
u(t) y(t)
The models can be of following various types.
. i
CHAPTER 4.2 Computer Control of Processes )
2. Frequency Analysis,
The input is sinusoidal. For a linear system the output will also be
sinusoidal at steady state. The change in amplitude and' phase will give
the frequency response (Example: bode plot). From the frequency response
the model will be estimated.
parameter vector. In this, method estimated parameters will have some yet) Lh(K),u(t-'K) -+- vet) ... (4.1)
, physical insight of the process. K=O
, Thefour different nonparametric methods for system identification are Where utt) Input '
KI-----~--~~-
Equation (4.1) is used as a model in spectral analysis. The model can
be estimated for.arbitrary inputs by dividing the cross-spectrum (between
output and Input) to the input .spectrum.
Where,Y(s) Laplace Transform of the ~utput Signal Yet) Where, yes) = Laplace Transform of the Output Signal yet)
Ll(s) = Laplace Transform 'of the, Input signal urt) Uts) = Laplace Transform of the Input signal uit)
G(s) Transfer function of FODT system G(s) Transfer function of second order. system
T Time constant K = Steady state gain
r' = Dead time =
0)0 Undamped natural' frequency
K = Steady state gain o= Damping ratio
Apply a, unit step as a input to FODT system as shown in Figure 4.'3 Apply a unit step asa input to second o~der system as shown in
and obtain a step response which is shown in Figure 4.4. Figure 4.5 and obtain a step response which is shown in Figure 4.6.
FOOT system
'2
K Output
Unit step input Kmo Output
Unit step input ""'"tS
---e 2 2
T8'+ 1 u(t} .8 + 2omoS + 0)0 Y(t) .
u(t) Y(t) ,
Fig. 4.3: Unit step Inpu.t to FOOT system Fig. 4.5: Un.it step Il1put to second order system
Flg. 4.6: Step response of a second order system to a' unit s~ep input'
The gain K is given by final value (after'convergence) as shown in Apply a following sinusoidal Input utt) to a system (described in (4.9)
-Figure 4.6. "I'he maxima and Minima of the step response occur- at times. as shown in Figure 4.7.
ri(t) a sin (rot) ... (4.10)
... (4.4)
Where a -amplitude of the sinusoidal inpututt)
co = Frequency of the sinusoidal input utt) in rad/sec
n = 1, 2, ... ~
Sinusoidal input .Sinusoidal output
and that ... (4.5)
u(t) Y(t)
where the overshoot M is given by
Fig. 4.-7: SinusoidallnputtoG{s)
~S1t ] If the system G(s) is asymtotically stable, then the output y (t) is
M = exp [ ~1- 32 ... (4.6)
'also a sinusoidal signal.
yet) b sin (cot + <p) ... (4.11)
Note: First maxima occurs at n = 1, First minima' occurs -at n = 2,
second maxima occurs at n= 3 and so on.. Where, b = amplitude of sinudsoidal Output yet)
From the step response as shown in Figure 4.6, the first peak time <P Phase diff~rence between Input u(t} and output
tl and first peak overshootM can be determined from (4.5). yet)' (Shown in Figure 4.8)
-logM Y(t)
s = 1
... (4.7)
/" . . , '.r: Y(t)
[1[2 + (log M)2 J2 I 'f .
\
\
Then roo can be- determined from (4.4) at tl \
\
7t I
(since n =. 1) ... (4.8) J
t 1 ~1-o2 I
I
u(t) I
I
From (4.7) and (4.8) the parameters S and roo can be determined. \ I
\ I
' ......... /
t
= b sin (rot + ~) ... (4.16)
fh( '1:) a sin (ro(t-'l:))d'l:
o From above equation (4.11) and (4.12) are proved.
t (eiro(t-T) _ e-iro(t-T))
By measuring the amplitudes a and b as well as the phase difference ~,
= fh('I:)a 2i dt one can draw a bode plot (or Nyquist 'or equivalent plot) for, different co
o values. From the bode plot, one can easily estimate the transfer function
;i t
fh( '1:)( eiro(t":'-r) - e-iro(t-T))d'l: .
model Gts) of a system.
= ;i t t
e irot Jh(1:)e(-iroT)d1: - e-irot Jh( '1:) e-(-irot)d'l: Where . Yet) = Output, Signal
ult) Input Signal
o
'---y-------I
0
"-------y-----
G(iro) G(-iro) ,h(K) = Weighting sequence
vet) = Disturbance term
( 4.10] Computer Control of Processes )
( Non ParametricMethods of System Identification 4.9
't = 0,1,2, ...
V(t)
Where
r--~-------------- --~ N = Number of experimental data
I I
i I 00 I
nput :.-.... L h(K)u (t-K) + + ~.""'11_ ...Output r ('t) =
Yu Estimated cross covariance function between Y(t)
u{t) I K=O I Y(t)
I i
and utt) from input and Output data
~--------------------~ ru ('t) = Estimated covariance function of urt) from Input
System
~~O)""J'
/\
ryu (0) /\
1 N-max(t,o)
= N . L.
Y(tH)U(t)
t:=l-min(~,o)
r =0, 1, 2, ... ... (4.19)
For white noise,
rue 't) = 0 for t =J, 2, ...
r u (r) = Constant value for t = 9 (i.e., non zero value)
1 N-t
= N LU(tH)U(t) r u (0)= Constant value
t=1.
(Non, Parametric Methodsof System Identification ( 4.11 ) ( 4.12 J ComputerControl; of Processes)
.Si~ce input is awhite noise,the equation (4.20) can be written as <!>Yu(CO) ... (4.24)
Where
h(K) = ryu (K) K=O, 1, ... ... (4.22)
ru (0) Cross spectral density between Input u'(t) and
Output yet)
ryu (0)
(0)
. ru
'rh(O) ryn.(l)
h(1) <Pu (00) Spectral density of input u (t)
ru'(O)
,lh(~) ... (4.25)
From the above. equation, the weighting function {h(K)} (i.e, h(O), hf l), ... ,
L h (K)e
co
-iKro
h( 00) can be easily estimated.
K=O
4.5 SPECTRAL ANALYS'IS Now the transfer function H(e-iro) can be estimated from (4.24) as
'Theform of model used in spectral' analysis is
co
<Pu (co) = -
21t-r=_N
L ru{'t)e-
N /\ .
r rm
_
... (4.28)
the following relation for the spectral densities can bederi~ed'from cross
Using (4.19), equation (4.27) can be written as
covariance function described in (4.18).
, N N-max(t,o)-
$yU (00) = _1_.
21tN
L L N t=l-IDln
' . (-r,O )
Y(t+ -r)u(t)e-itro
,
... (4.29)
t="-
( Non ParametricMethods of SystemIdentification,
( 4.14 ) Computer Control of Proc~
[ 4.13 )
Similarly, .
Next make the substitution s= t + 'to Figure 4.10 illustrates how to
derive the limits for the new summation index.
1 N N . . "
Since --~" u(s)U(t)e-lSIDe-ltco
21tN LJL...J
s=.i t=l
We get,
... (4.31)
-N
Fig. 4.10: Change of summation variables
(summetien is over the 'shaded area)
... (4.30)
... (4.32)
Where,
N
YN(ro) = Ly(s)e- isro = 'Discrete fourier transformofylt) This transfer function is called empirical transfer function estimate.
s=1
N .
UN(ro) = L u(s)e- isro
= Discrete fourier transform of net)
5=1
TWO MARKS QUESTIONS AND ANSWERS 4. What is correlation analysis method for system identification?
1. What is .non parametr-ic method for system identification? The model used in correlation analysis is
. property that the resulting models are curves or functions,which yet) Lh(K)u(t-K) + v(t)
K=O
are not necessarily parametrised by a finite dimensional parameter
vector. The four different non parametric methods for system Where, Y(t) = Output Signal
identification are net) = Input Signal
1. Transient Analysis h(K) = Weighting sequence
2. Frequency 'Analysis vet) = Disturbance term
3. Correlational Analysis The input u(t) is a white noise signal, a normalised.crossconvariance
function between 'output and input will provide an estimate of the
4. Spectral Analysis
weighting function h(K).
2. What is Transient analysis merhodfor system identification?
5. What is spectral analysis' method for system Identdficatdon?
, In this method the input is taken as a step or impulse' and the
The model used in spectral analysis is
recorded output (example transient response ora FODT system to a
unit step input) helps to Identify a model. C()
'. o
<>
Ir
I' "
Jransient
response to
.~'J a.unit step}nput
h(t)
vet)
Weighting sequence
Disturbance term
J .' t.
The model can be estimated for arbitarary inputs by dividing the
j~ ~~I~ td ~I
0
cross spectrum (between output and input) to the input spectrum.
.' 'The above transient response ~(~Qorded datal-helps to find out FOl::>T
system parameter~ s~c~ asK" tJand td. e, (,
The method of least squares is about estimating parameters by Yet) == h o u(t) + hI u (t - 1) + .. ~~+Jlrti~l u(t - M + 1)
minimizing the squared error between observed data .a~d their expected values. The input. signal ut t ) ... u(t ~ '.M + .1) are recor~ed during thee
The linear regression is' the simple type of, parametric model. experiment. Hence. regression variables
This model structure can be written as
~ (t) '=, (urt) u(t - ,1) ... u (t'- M + 1) is a )\1 - vector of
Yet) = ~T (t) e ... (5.1)
.known quahtitie~. and, para~eter vector.
Where,
'0 ~ 'eho hI' .... '- 'hni_1)T is a lVI-Ve~tor of unknown
Yet) = Measurable Quantity " parameters to be estin.1at-ed.
4>(t) = n - Vector of known quantities
[~Y(t-l), -yet - 2), ... ,.-Y(t - n a) utt - 1) ... -u(t- nb)]T
The .problem: is to. find an tistiJl}.ate e of the parameter vector e as
shown in (5.1) from experimental measurements Yi I), ~(1), Y(2), $(2), ...
e n . Vector . of unknonwn parameters Y (N) ' ~(N). Here 'N' represents number of .experimenta~ data and 'n'
The following two examples shows how a model can be represented represents number' of known quantities in <p(t) or nu.mberof unknown
using linear regression" model" form; parameters in O.
(Parametric Methods of System'Identification 5.4 Computer Control of Processes"]
5.3
and stack these in a vector E defined as
Given experimental measurements, a system of linear equations is
obtained as
Y(l) = <l>T (1) o
Y(2) = <f>T(2) 8
Note:
an (N x 1) vector ." (5.5a)
The other form of loss functions are
... (5.8a)
One way to find e from (5.4) would of course be to choose the number where IHI denotes the Euclidean vector norm.
of measurements, N to be equal to n. Then <I> become's a square matrix. If
this matrix is non singular the linear system of equations (5..4) could easily Thee will be estimated from experimental measurementsYt I), p(l);,
be solved for O. In practice, however noise, disturbances and model misfit Y(2), ~(2) ... , Y (N), ~(n) by minimizing the loss function V(S) in OJ.7) and
are good reasons for using a number of experimental data 'N' greater than (5.6). The solution to this .optimized problem is "
en'. With the possible to get an improved estimate. When N> n,exact
solution for linear system of equations (5.4) in general not exist.. ... (5.9)
To find an estimate e, introduce the equation of errors as For this solution, the minimum value of Vee) is "
Note: Disturbances
1. The matrix $Tep is positive definite.
2. The form (5.9) of the least squares' estimate can be rewritten in
the equivalent form
... (5.11)
5~3 RECURSIVE IDENTIFICATION METHOD Fig. 5.1: A general scheme for adaptive control
In recursive (also called on-line) identification methods, the Most adaptive systems, for example adaptive control systems as '
parameterestimates are computed' recursively in time. This means that shown in Figure 5.1 are based (explicity or implicity) on recursive
,,'.'" _:';"'. C'-'.;~<'",' ' , ' 1\
if' thereis an e-stimate"6(t-l) based on data upto time (t - 1), then '9(t). is identification.
, 1 \ '
computed by sOll1e':(simp!~ modification' of_,9(~-1). Then a current estimated model of the process is available at all
times. This time varying model, is used to determine the parameters' of
The counter parts to one-line methods are the so-called off-line or the (also time-varying) regulator (also called controller).
batch 'methods, in which' all the recorded data are 'u~ed simultaneously to
find the parameter estimates. In this way the regulator will be dependent on the previous' behaviour
of the process (through the information flow : Process ~ model -7'
Recursive identification methods have the following .general features:
regulator).
They are central part of adaptive systems (used, for example,. for If an appripricate principle is used to design the regulator, then,
control or signal processing) wherethe action is based on ,the most
the regulator should adopt to the changing characteristics of the
recent model.
.process.
Their requirement on primary memory is quite less compare to offline
The various identification methods are
identification methods' which, require large memory.tostore entire data,
a) .Recursivc least squares method
They c~n be easily 'modified into real-time algorithms" aimed at
tracking time-varying parameters. b) Real time identification method
, They can be first step in a fault detection algorithm, which is used c) Recursive instrumental variable method
to find out whether the system has changed significantly. ' d) Recursive prediction error method
( Parametric Methods of System.Identification 5.7 5.8 ) Computer Control of Processes )
_Then the least squares parameter estimate is given by S(t) = P(t)[~HS)Y(S) + Ht)Y(t)] ... (5.19)
The argument t has been used to stress the dependence ofe on time.
The (5.14) can be computed in a recursive fashion. By substituting (5.16)
Note that in (5.22) there is now a scalar 'division (scalar inversion) The LS estimate of parameter vector q is
instead of matrix inversion. From (5.21b) and (5~22),
= P(t-l)~(t)
K(t)
[1 + <j)T(t)P(t-1)~(t)J
(Parametric Methods of System Identification' (5.11 ) Computer Control of Processes )
4. Where the limitations of LS method? The various recursive identification methods are
The LS method gives consistent .parameter estimates only under 1. Recursive least squares method
following restriye conditions,
2. Real time identification method
If E<I>(t) <l>T(t) is non singular
3. Recursive instrumental variable. method
If E <!>(t) v(t) =0 4.' Recursive prediction error method
If .input is persistently exciting of order nb, ,7. What are the advantages of recursive identification method? .
, .
&.1 ,REVIEW OF Z-TRANSFORMATION If we introduce a new variable Z = eST in equation (6. 7) we get
co
"We found laplace tr ansforrns to be helpful in' the analysis of
conventional - control systems. In subsequent chapters we will see that f*(s)lz~eST = Lf(nT)e- nST .. (6.8)
n=O . '
theanalysis of sampled-data control problems is conveniently handled in
terms ofZ-transforms. This result is defined to 'the Zvtransform of f(t) and is denoted as
F(z). Thus
In this chapter we shall develop some of the important concepts of
Z-transforms. In this .chapter we shall develop some of the important 00
concepts of Z-transforms. To begin, recall that the laplace transform of 'a F(Z) = z{f(t)}= f~(s)lz=eST = Lf(nT)z-n ... (6.9)
function f(t) is defined as . n=O
a:>
0Cl
By definition
f*(t) = f(t) L8(t-nT) ... (6.2) .
L U(nT)z-n
C1J
o
Z{U(t)} = F(Z) = ... (6.11)
00 n=O
f*(t) = Lf(nT)8(t.... nT) 00
t<Ol
... (6.14)
f(t)
{Si:at t z OJ
... 0 (6.20)
.. ~ (6.16)
l'
IZ-l! < (6.17)
~.
J n=O
f( e JDaT _e-janT)z-n .0. (6.20)
F(Z} = 1- e- Z-l 'at eaT
Ramp Function 1(
2j 1- e j aT z'
1
l-e-~T z' J
f(t) ... (6.18)
jaT
[ 1- e Z-l .,.. 1+ e jaT Z-l ]
=:
(e
1 + Z-2 - Z-l j aT +e -jaT)'
F(Z) = Z{KT} :::: 'IKnTz-,n
n=O .
z' (e jaT _. ~jaT)
_ KT {Z-l +2 Z-2 + ... }
2j
KT Z-l {l + 2Z- 1 + 3Z- 2 + 4Z- 3 +...} ..' (6.18)
1+ Z-2_ 2Z- 1 (jaT
e ~e -jaT)
KTZ-1
=
(l-Z~l t Z-l sin aT
= 1 + Z-2 - 2~-lcos aT
Z-l Biua~
Z{sin at} .00 (6.21)
1 - 2ZcosaT + Z-i
F(Z) =KTZ ... (6.19)
(Z_1)2 ZsinaT
Zlsin at} ... (6.2,2)
Z'-2 - 2ZcosaT + 1
(Digital Controller Design 6.. 5 6.6 ) ComputerControlof Processes ]
F(Z) _ Z(Z2 +Z+1) Now, recognize that the presence of areal number other than one in
... (6.24) the denomination of the Z-transform expension willgive rise to an
Z - (Z-O.5)(Z-1)(Z-(O.5+ jO.74))(Z-(O.5- jO.74))
exponential function, in its inverse.
F(Z) _ _A_+~+ . C + D Therefore, we must express the number 0.5 in exponential form.
Z - Z-O.5 Z-l .Z-(O~5+.jO.74) Z-(O.5-jO.74) Similarly, complexnumber give rise to exponential and sine/cosine terms.
Thus, by converting from rectangular to polar coordinates.
'By comparison ao = f(O), al = f(T) and an = f(nT) which is the value of We .know that Z{e-~Ts}= Z-N
f(t) at the nth sampling instant. ,. Gp(Z) = Z-N Z{G(s) e-9S} ... (6.35)
Z-l e- amT
Zm{e-at } = 1 - Z-l e-a
T ... (6.41)
{ I .} {e-as}
Z'-=z-
m s+a s+a
... (6.39)
y (Z) =
X(Z)
G(z) ... (6.43)
IU(nT-8)e-a (nT- a)
6.5..1 Complex Series Representation of the Sampler
n=O
Recallthat the sampling operation gives us the values of the output
=e-a (T - 9) Z-l + e-a (2T - 9) Z-2 + ... ... (6.39) at each sampling instant and can be represented as
LC
00
4- f
, and the Fourier Coefficients are given by
OT(t) = ejnwsl ... (6.52)
n=-oo
... (6.47)
and X*(t)= 4f n=-'oo
X(t)~jnWsl ... (6.53)
From this Figure By the same reason in which led to Equation (6.54), this can be written
as
yes) = x*(s) G(s) ... (6.55)
Y*(s) =, X* (s) G*(s) , ~ .. (6.63)
By analog to Equation (6.54)
Taking Z transform .of this Equation gives
~ I
n=-oo
X*(s+jnws)g(s+jnws) ... (6.59)
2. For X (t) = S(t) ( i.e., the impulse function)
Y*(s)= ~ i
n=-<X)
G(s+jnws)X*(s) . ,
n=O
n=-co
.
... (6.62) function, it may be split into partial fractions then steps 2 and 5 can be
applied to the individual partial fractions.
Computer Control of Processes) -
(Digital Controller Design [ 6.15 )
Xes)
1 1
Y(s) fig. 6,,'1: Block diagram of a closed loop system
s+1 s+2
1- e Ts
G(Z) =
00
Ly(nT)z-n L(
co
e-nT'_e-2nT) Let Gho(s)=
s
... (6.70)
n=O . n=O
From Figure
1 1 C(s) = G(s) Drs) E*(s) ... (6.71)
GCZ) = 1 _ e-aT Z-l 1- e-2aT Z-l
Or
After ,simplification we get C*(s)= G*(s) D*(8) E*(s) .~. (6.72)
( ~Z-1)Z(s~ J
TZ-1
=
(l-Z-l)
computer Control of processe;],.
[Digital ControllerDesign ( 6.19 J
The inverse of which is e-at
Substitution ofR(Z) in the closed-loop pulse transfer function given
where
C(Z) = (Z-l) [Z+(K c T-l}J = ZC 1( z-)
1 [
(l-I~c T)J
1 l
a ln
T
Now we expand C 1(Z) in partial fraction., Thus,
, Kc T A B ., C(t)
(Z-l) [Z+(Kc T-l)J= (Z-l) + Z-(l-K c T)
u(t) _ e1n(1 - KcT)tlT
Fig. 6.4:. Closed loop sampled data system 6. Dead time compensations
We examined the detrimental effects of dead time on the response'
2. ResetProblem of first order system. The presence of dead time necessitates lowering of
The algorithm must guarantee offset-free performance in the controller gain to. maintain stability.
presence of modeling errors .. Modeling errors arise when the nominal W e observed that w hen the apparen t dead tim e ed exceeded the
model upon which the algorithm is based is different from the Teal process. dominant time constant of the system t the peak offsets following a step
3. Controller Tuning change in load could, approach those of the uncontrolled situation even
with best PID controller tuning. Under these conditions, the settling time
It is desirable that the algorithm. contain a small number of
approaches g8d.
independent (adjustable) tuning parameter. The larger the number, the
more difficult will be the task of keeping the controller tuned in the field. Since many chemical engineering process exhibit apparent dead time
characteristics and since dead time is detrimental to control, there is
4. Robustness considerable incentive to develop control algorithm that can compensate
A digital control algorithm is based on a nominal plant model that 'is for such time delays.
developed from first principle or 'obtained from experimental data.
The concept of dead-time compensation involves the use of a
Modeling error willfnvariably be present. Further more, new modeling
mathematical model of the process to in effect remove the dead time from
, errors will creep in due to equipment' fouling changing production
the feedback signal so that the performance of the system is improved.
strategies, etc:
The algor-ithm must maintain stable proceas operation in the 7. Working with RHP zeros
presence cf such modeling errors. Robustners issues having to do with When the roots of the numerator polynomial of the laplace' domain
the ability 'of the control system to maintain.stability in the presence of'a transfer function contain one that lies to the right of the imaginary axis
plant model mismatch, have received consideratble _attention in recent in the s-plane, the system is said to have a RHP(right half plane) zero.
years. Such a system containing all RHP zero. begins in a direction that is
,:';. ConstraintHandling opposite to the direction in which it eventually settle out. Inverse response
occurs in boiler level systems and in some distillation base level systems.
Industrial process must often be operated in the neighborhood of
operating costraints. There are two types of constraints. There are two
types of constraints:
(6.24] ComputerControl.ofProcesses)
{Digital Controller Design ( 6.23 )
To develope the algorithm for set point chages recall from study
It should be clear that such systems will pose'difficulties for the that the closed loop pulse transfer function of the system of Figure 6.4 for
PIn type controller in that the controller will take wrong action intially set point changes is
with digital control. A facility to accommodate inverse response can be
built into the algorithm and this should be considered wherever applicable.' C(Z) D(Z)G(Z)
... (6.79)
8. -Manipulated variable -movements R(Z) 1 + D(Z)G{Z)
Excessive movement of the manipulated variable should be avoided where, G(Z)' = zro., (s)Gp(s)}
to reduce actuator wear.
The solution of Equation (1) for n(Z) is
9. Shaping of closed leep dynamics
In the absence of modeling errors the algorithm should yield a closed C(Z)
loop response having desired dynamics. 1 R(Z)
D(Z) = G(Z) 1 _ C(Z) ... (6~80)
10.' Reset windup R(Z)
If the algorithm contains- integral mode the phenomen.in of reset
windup may occur. Anti-reset windup features should be provided' The design, procedure is to specify the desired response characteristic
wherever 'appropriate. C(Z) , ..... - .
R(Z) , for example the controlled variable shall ,reach the new .set point
11. Bumpless Transfer
in one sampling period. Then Equation (6.80) may be solved for D(Z),
The facility for bumpless transfer from manual to automatic should
provided the process transfer function Gp(s) is known.
be provided.
T 2T 3T 4T Time
Q)
IExample 6.4 ,
C/)
c: Design a dead beat control algorithm for a process whose transfer
8.
CJ)
(\)
0:::,
function is given by Gp(s)' = 1 . , with T = 1 sec.
0.4.8 + 1
1.0
Solution:
1
Gp(s)=
0.48 +1
C(Z)
n(Z)
1 R(Z)
= - - x _--:.--:.
G(Z) 1- C(Z)
R(Z)
1 Z-l
... (1)
G(Z) 1 -- Z-l
T 2T 3T 4T 5T Time We know that,
(b)
~.6:
Fig. Response specifications for dead beat design
with process dead time (a) i d =T (b) td =21 =
(l_e-
zJ---x--.
TS 1}
l s 0.48+1
C(Z)
= Z-2
R(Z) =Z 1 _Z_e_ -TS_ }"
{ s (0.48 + 1) } { s (0.48 + 1)
E(Z)
or . R(Z) = Z-3 etc.
z{ 1
8(0.48 +1)
}-z-lzf, 1
l8 (0.48 +1)
}
Yielding following controller design equations.
I
,1 Z-2,
n(Z) = ---- .t, (6.91) l-Z-1)Z 1. I
J ... (2)
G(Z) 1- Z~2 ( { (O.4s + 1)
'S
1 .Z-3 2.25 A B
or ... (6.92)
. D(Z) = G(Z) 1 -' Z-3 8(8+2.25). = s + 8+2.25
( Digital.ControUerDesi,gn ', (' 6.30) ComputerControl of Processes )
Solution:
Now (1) =>
10
1 z' where t = Is Gp'(S)=
0.58 + 1.
D(Z) = G(Z) 1 _ Z-l
e(z)
= (1_e-Z.5T Z- 1) x _ ~ 1 R(Z)
~(1_e-2.5T) 1- Z-l
D(Z) = G(Z) =
1 + C(Z)
R(Z)
(1- e- Z-l )
2 .5
1 Z-2
(1-Z-1)(1- e-Z.5 ) = --x
G(Z) 1- Z-2
... (1)
TS 1 Z-2
orz: = z{ l-se O.:SO+l} (1) =>
neZ):::: g(Z) 1- Z-2
Z-1{1_e-2T Z - 1)
10
10
~=
20 10(I-Z-1 )(1 ~ e-2T )
8(s+2) s(s+2)
Substitute T =1
20 A
-+-
B (1- e--2 Z-1)
. Z-1
s(s+2) - 8 8+2
10(1-Z-2 )(1- e-2 )
20 ~ A(8 + 2)+ Bs
Z-1 (1- 0.1353 Z-1 )
Put s= 0 Put =-2 10 (1--Z-2) (0.8646)
20= A2
M(Z) Z-l (1 -0.1353Z-1 )
I A = 10 I I B = -10 I . E(Z) 8.646(1- Z-2)
(2) => {l_.Z-l)z{lO _~} M(Z) 8.646 - M(Z) 8.646 Z-2 = E(Z) - E(Z) 0.01353 ~2
- .. s s+~
M(Z) = M(Z) Z-2 + 0.1156 E(Z) Z-1_ 0.01353 E(Z) Z-2
= (1 _Z-1) Z{lO U(t) - 10 e-2t } Mn = Mn - 2 + 0.1156 E1:<-1 - 0.01353 E n - 2
10Z.-1(1- e-2T ) F(Z) represents first order lag and N is the number of sampling
periods in the process dead time 0d. To derive the expression for F(Z)
(1- e-2T Z-l) consider the differential equation for a first - order process.
fDigit~f ControllerDesign ( 6.34 J'-' Computer Control of Processes )
[6.33 )
.1 EXample 6.6
. 'tf
.
dy +
dt.
yet) = X(t). . (6.85) e-1. 4 s
Design a Dahlin's controller for Gp(s) = T =1 sec.
3.34s+1
Note that in Equation (6:85) the steadygain have gain of one. The
Solution:'
numerical solution of Equation (6.85) is
Dahlin's Control algorithm,
Yn (l~af)Xn+.afYn-l ... (6.89)
where ar = e- Thr (O < af< 1)' D(Z) =
1
--x----'--------
1- aZ-(N+l)
As we will see, Dahlin's algorithm is essentially the same as the -NTs -es 1
recently developed internal model control for first order plus dead time . :. GCZ) = (1-Z-1IZ e e
types of proceses. It turns out that elf can be used to enhance robustness
\ J { 3.348+1 f
of the loop in the 'presence of modeling errors. We know that, Z{e- NTs } = Z-N
A high value of elr close to 1 gives a high degree of robustness, but
( -6s }
the response become quite sluggish. Smaller values of ar close to zero G(Z) = (l_Z-1)Z-l zJ__e _
improve dynamic response, but the system becomes more sensitive to .\ . . l8(3.34s+ 1) "
modeling error.
[Digital Controller Design ( 6.36 ) Computer Control of Processes)
( 6.35]
(l~Z-l)Z-lZ {_ _
m s(3.348+1)
l_} D(Z)
1
G"(Z) 1 ~ aZ-
1
(1 -
= - - x----C.---=---------'-
-
a)Z-(N+l)
(1 - a)Z-(N+l)
We know.that,
m s(1+0.3)
f -0.741z-1 l[0.6321Z-
2
]
[ O.165Z- 2,+O.094Z-1 'I -O.3678Z-1 - O.6321Z- 2 xz
0.6321Z-1 -'o.468Z-2
0.094+0.1305Z-1 - O.12 Z- 2 - 0.1043Z-2
Where a = 0.3, 8 = 0.4
m= 1 - ~ =. 1 - 0.4 = 0 n(z) = M(Z)
T l' E(Z)
Substitute the above values into above equation
G(Z) = (l_'Z-l)Z-lZ-l
, , {I e
-(O.3xO.6 xl)
_
}
M(Z)
E(Z)
0.6321Z-1 - 0.468Z-2
0.094 + O,.1305Z-.1 - 0.12Z-2 - 0.1043Z-3
. ,..' . 1- Z":'1 1- e-(O.3xl) Z-1
00094 M(Z) + 0.1305 Z-l M(Z)
(
1 _ Z-1)Z-2{_ 1 _ - . 0.835
; 1-Z-1 1-0.741Z-1 J
1 - 0.12 Z-2 M(Z) - 0.1403 Z-3 0.031 Z-1 E(Z) - 0.468 Z-2 E(Z)
TakingInverse Z transform
1 1
'Mn - 1 = Controller output at the (n_l)th sampling instant
G(Z) = Z"'2[1-O.741Z- - 0.835 + O.835Z- ]
1 -O.741Z-1
, .
Mn - 3 = Controller output at the (n_3)th sampling instant
Z~2 [ 0
1 65 + 0.094'Z-1 1
1 - O.741Z-1 -J
En
E n- 1
= Error at the nth sampling instant
Error at the (n_l)th sampling instant
2 3
G(Z) = O.165Z- + 0.094Z-
1- 0.741Z-1
(Digital Controller Design [ 6.37 ) 'tomputerC9ntrol of Processes J
~ample 6.7 1
We know that,
. . e-O8s
Z{e-:-NTs } = Z-N
Design a Dahlin's controller algorithmn for Gp(s) = -0.6s+1
- - ; T = 0.4 'sec
. Solution: .
acz) = (1- Z-l )Z-2
.
z{ 8(0.6s+1)
1 .}.
Dahlin's Control Algorithmn
= (1-Z-1)Z-2 Z{ 1.67.. }
1 (1 - a) Z-(N+l) . s(8+1.67)
n(Z) = --x--~--..:...---
. G(Z) 1 - aZ""! ~ (1 - a)Z-(N+l) By using partial fraction method,'
where,
z{ 1.67 }
8(s+1.67)
_ A +_B_
- s s+1.67
G(Z) =
Z{l-e..s-TS 0.6s+1
e-O.. 8~.}
x
1.67. = A(s + 1.67) + Bfs)
=> A= 1; B = -1
-O.8s }
_ 1-Z-1 Z e
- ( ).{ s(0.6s+1)
:. GCZ) = (l~Z-l
.
)Z-2 Z{-!_ _
S 8+1.67
l_}
In this example Gp(s) contains dead time Od.
ed= NT
(1 - Z-1)Z-2{ Z-l
Z - Z}
In this. example: Z_e-l.67xT.
. (1-Z-1)Z-2 [~_ . Z ]
Z-l Z-O.51
1 - Z~lZ -:NTs
e. ]
. = .Z2[1_ Z-0.51
Z-l ]
(- ) [ s(0.6s+1)
= Z-2[Z-0.51- Z+l]
Z - 0.51.
(Digital Controller Design (6~40 ] ComputerControlef Processes J
Let D(Z) is the ratio of two polynomials in z-', as demonstrates
Z-2[ 0.49 J C(Z) f30 + f3IZ- 1 + (32 z-2 + + r\. Z-k
'Z - 0.51
D(Z) = -()
E Z
= . 12m
1+ al Z- + a2 Z- + + am Z-
... (6.95)
O.49,Z-2
The roots, of the polynomial in the' denominator are the poles of the
Z - 0.51
controller and for the stable controller, these poles should be inside or
1, (1 - a )Z-(N+l) on the unit circle inthe complex plane. Let us examine what the effect is
DCZ) = G(z) x 1 _ aZ-1 _ (1 _ a)Z.,.CN+l) of a pole's location on the controller outputCrt).
Consider a controller with the following "discrete transfer function \
where, a = e~Tlt ''t' is not given.so.r ~T e"':' l = 0.3678
C(Z) 1
D Z = Z-O.51 [ " (1 ~O.3678)Z-3 , ] DCZ) = E(Z) 1-PZ-l
... (6.96)
( ) . 0.49Z-2 1-0.3678Z-1-(1-0.3678)Z-3
An er~9r equal to' a .unit impulse enters the controller at t = o.
O.6321Z-2 -O.322Z-3 Then from the above equation, C n - P Cn- l = En
=
0.49Z-2 -O.18Z-3 -O.31Z-5
where En::;: Lfor ri e G and Ej e D forn=1,2,3 ...
Let us now invert this algorithm into thetime domain. To accomplish
LetP = -0.9, then Cn + 0.9 C n-1 = En
this we note that
Controller output, Cn
n(z) = M(Z)
E(Z) n En ForP = -0.9 For P = ~10.9
Imaginary
p=~().9
+j
Unit circle
-1 +1
(e) (d)
IExample 6.8 I
Fig~ 6.7: EHect of pele.Ieeefien on liehaviour of controller output
Eliminating the ringing of a digital controller. Suppose that for a"
For positive poles, the controller outp.ut exhibits no ringing (as shown given process" Dahlin's method leads to a eorrteoller with the transfer
in Figure 6.7 (c) and (d). Thus from above discussion: . function.
'Any negative pole of the controller's discrete transfer function
will cause ringing. .... (6.97)
The closer a ringing pole is to the' unit circle, the higher the ringing
. of the controller output will .be. The term (1 + 0.6 Z-l) in the denomimator implies the presence of a
pole at -0.6, which will cause controller ringing.
. Hence, we should design a digital controller, so that all its' poles are'
poaitive, preferably with medium - size absolute values as shown in Dahlin has suggested that the ringing pole can be eliminated by
, .Figure 6.8. multiplying D(Z) with (1 + 0.6 Z-l), thus canceling th~ ringing term from
the denominator. But inorder to keep the static gain K the same.we should
Unfortunately, the deadbeat and Dahlin's algorithms usually contain
poles that causesevere ringing of the controller output. This may be the' divide DCZ) by the term lim(1+0.6Z-1 ) .
Z~l
= 1.6.
case of any' response, - .specification design technique 'and is not unique to
Thus the ringing - free controller has the transfer function'
the foregoing two algorithms.
,,(Digital Controller Design ( 6.43 ) [6.44 ) Computer Control of Proces'ses.] .-
tI 0 " dt
The digital controller transfer function is,
To convert this expression to its digital equivalent, we use the
following finite difference approximations
t 'n DCZ) = ~g}= Kc[1+~:(1_~_1)+ :~ (1-Z-1)] ... (6.106)
Equation (6.102) can also be written aes a Z- transform ,expression. The velocity form of the ,controller' can be found by shifting (6.102) ,to
obtain Pn-l and subtracting it from (6.102).
Let pIn = Pn -P be defined as'a deviation variable.. Recall that the
Z-transform translation theorem yields
~Pn = x, [(en - en-l)+ AttI en + t~tn (en - 2e n. _1 + e n-2)]
'
... (6.109),
Z(en ) = E(Z)
,Z(en-l) = Z,...l E(Z) ... (6.103) Taking the Z-transform of (6.109) given
AP(-r)= Kef (1- Z-l )E (Z)+ ~: E(Z)+ ~ {1,..;2Z-1 +Z-2)E(Z)] ... (6.111) .
(DigitalController Design .1 645 ) ( 6.46 ] Computer Control ofProcesses )
'Note that Equation .(6.11) can be obtained by multiplying both 'sides For the position form algorithm, there are several modifications that
of (6.105) by (1 - Z-l). Another variation of the digital PID controller is can be made .to reduce reset windup.
based, on the more accurate trapezoidal approximation for the integral a) Place an upper limit on the. value of the summation. When the
(in Eqution 6.99) . - controller saturates, suspend the summation unt.il the controller
output moves away from the limit.
!e(t)dt "' ~(ek +2ek- }t 1
b) Back - calculate the value of en that just causes the controller to
saturate. If the saturation occurs, and this actual as error term
After ,this expression is substituted into(6.99), the velocity form of en-l in the next controller calculation. Field testing has' indicated
the control law becomes, that approach (b) is supperior to (a), although it is somewhat more
complicated.
APn"' ~ x, .[(en - en-l) + ~t.
'tI'
(en + en-I) + 'tn (en -
2 8t,
2 en-l + .en- 2 )]. Note that in the velocity' form of the algorithm, Equation (6.109) no
summation appears avoiding the windup. problem. The controller remains-
or as a Z-tranform (Equati~1?>_Q .ll_3J'
saturated until the error (en) decreases to a point where the control action
bn and the signal, to the final control element P n are computed in terms C - em is 'the error Em, which arises because of modeling errors or load
of percentage full range K, willbe dimensionless. This will also be true if upsets.
the computer input .and output have the same. units (Example rnA or V).
(DigitalController Design ( 6.49 )
[6.50. J Computer Control of Processes )
L
L
rr L
G'm(s) = G'p(s)
c
and 8 m = ad; equation (6.116) reduces to.
1 - e- st
Where, Transfer function of zero-order hold - - -
s
Kp = Process gain
- o K e-9 dS
ad = dead time Transfer function of process, -..:.P_-
rs-i-I
t = Process time constant
Substituting the transfer functions in Equation (6.119) we get
The smith predictor scheme for this computer control application is
shown in Figure 6.14 .. The implementation of this scheme would require
that we develops equations for BM , eM and UK. ' ... (6.120)
(Digital Controller Design " [ 6.53.) ( 6.54. ) Computer.Control of Processes]
If we denote the integral .number of sampling periods in dead time Cross multiplying and inverting gives'
, as N, then the following equality holds
... (6.124)
ad = (N + (3) T ... (6.121)
From Figure (6.14) the error is computed from the equation
Where, f3 = a fraction
between 0 and 1. With this expression .for ed,
equation (6.120) can be written as ... (6.125)
[:_,_1_~
Thus, dead time has been compensated. The digital controller must .
C,M(Z)=
U(Z)
K p
Z-N(l_Z-l)Z-l
1- Z-l
e-(l-fl)T ]
!
... (6.123)
bea conventional control algorithmCi.e., P, PI or PID). If the controller is
based on Z transforms, it incorporates the effect of dead time, and the
. .' 1-etZ~
. . smith predictor algorithm is not applicable.
Cross multiplying and inverting gives an equation for the model
output containing true delay. 6.5.10 Internal Model Control
Thus, Internal model control (IMC) is a comprehensive model based des~gn
method developed byMorari and coworkers. The I~,1C method is based on
CM,K = A ZCM,(K-1)+A zk p (;3-l)Uk-(N-Z) +kp (1-1:)U
k-(N+li an assumed process model and leads to analytical expressions for the
controller settings. Both IMC method and direct synthesis method
.. ~ (6~124)
produce identical' controllers if the design parameters are specified in a
Where consistent manner. However, the IMC approach has the advantage that it
allows .model uncert aining and. trade offsbetween performance and
A2 = e-TIt robustness to be considered irra niore systematic fashion.
A3 = e- J3Th
The IMC method is based on the following simplified block diagram
Now, we need to get the model output containing no dead time. For as shown in figure 6.15.
this' case, .
A process model G and the controller output P are used. to calculate
1 . .:. e- sT K p ] the model response, Y. The model response is subtracted from. the actual
Z [ s. r s j I response. Y and the difference, Y - Y is used as the input signal to the
IMC controller, G*c. I~ general, due to modeling errors (G;tG) and
kp(1-Z-
.'
1 ) [_1_ --
l-Z-1
' 1. J
l~e-T/'t. Z-l
... (6.123) unknown disturbances (D:;t:O) that are not accounted for in the model.
(Digital Controller Design (6.55 ) ( '6.56 J Computer Control of Processes )
Where G+ contains any time delays and right half plane zeros.
In addition, G is required to have a steady state gain equal to one inorder
(a) Classical feedback control to ensure that the two factors in (6.131) are unique'.
o
Step 2:
The controller is specified as
I------,r------+ y
... (6.132)
y
1
f ... (6.133)
(b) Internal model control In analog with the direct synthesis method,
Fig. 6.15: Feedback centrel strategies t c is the desired closed loop time constant
From Figure (6.15) it can be shown that the two block diagrams-are r is positive integer (usually r = 1)
dentiealifcontrollers Gc andG*c' satisfy the relation
The IMC controller in (6.132) is based on the invertible part of the
G* process model; G_ rather than the entire proces model, G. If G had been
Gc = ' c ... (6.128) used, the controller could contain a prediction term e 9s (if G+ caintains a '
1- G:G
time delay 8), or an unstable pole (if' ,0-+ contains a right half plane zero).
Thus any IMC 'controller Gc* is equivalent to a standard feedback Thus, by employing the factorization 'of (6.131) and using a filter of
controller G c and vice' versa; From Figure (6.15 (b)) , the closed loop the form of (6.133), the resulting controller o', is guaranteed to be
relation for TMC is physically realizable and s,table. -
In general, the noninvertible part 'of the model, 0-+, places limitations
... (6..129)
on the performance that can be achieved by any control system. Because
the standard IMC design method is based 'on pole-zero cancellation, the
For a special case -of perfect model, G-= G ' IMe approach' must be modified for process that are open loop unstable.
Thus the closed loop transfer function for setpoint changes is:
G(s) = ... (6.140)
y -
-'= G f ... (6.135) .
Ysp +
The IMC' and direct synthesis (DS) methods can produce equivalent Factor, this. model as
controllers and identical closed loop responses, .even when modeling .. .. .. _. ' e
errors are present. This equivalence occurs if the desired transfer function G = G+ G - where G+ = 1 - -s
2 .. 0 (6.141)
Gc = ... (6.136)
(1 + ~S)(~S+l) ... (6.143)
,
l-(lJ Ysp d
K(t cs+l)
IExample'.' I
Use IMC design method to design two controller for the FOPTD ... (6.144)
(First order plus time delay) model given by ,
_ . Ke-es
G(s) = - , - .. ~ (6.137) Rearranged into the PID controller of
rs-tI
Ke= .l[_
K
. 't c
't
+8
), ... (6.147)
Considering only the feedforward path (and ignoring the feedback
loop), the appropriate equation for perfect control (with R = 0, C = 0 and
L *- 0) is
LGL (Z) + LGt (Z) GF (Z)HGr Gp(Z) = 0 .: (6.148)
A comparison of (6.14.6) and (6.147) indicates that the type of
controller that is designed depends on the time delay approximation. The Solving for Gf to obtain the FF controller gives
IMC controller in {(6.~47) is identical to DB controller for the" FOPTD
model. ... (6.149)
Gt<Z) = LGt (Z)HGv Gp (Z)
6.5.11 ~igital Feed Forward Control 'I'his result is net as attractive as for continuous control, because
L . 'LGL{Z) GL{Z)
the load variable cannot be factored area cancelled (i.e~,). LGt{Z) of. G (Z)
Computer t
_. HG L LG L (Z)
in gen.eral. However, the ratio HG .may be used in place of LG (Z)
t . .t
where L is a step input or a sequence of steps (piecewise constant input).
I
I
This assumption removes the requirement that Equation (6.149) be
I
I revaluated for every different input L(Z). The resulting expression will
. I -
RJX~+
be reasonably accurate except when the load change is not well
:.~~- I
I c approximated by the piecewise constant function.
.1L _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ JI
B
Gt<Z) = -HGL (Z) ... (6.150)
Kt Z-T HGv G p (Z)
Fig. 6.16: Block diagram for digital feedforward controller
[Digital Controller Design [ 6.61 ) ( 6.62 ) ComputerControlof Processes)
The feedforward controllerwill not be realized unlessBj, > (St + 8p ) , 'fWD MARKS QUESTIONS AND ANSWERS
where Sp is' the process time delay. An unrealizable controller is indicated 1. Define Z -transform.
bya positive power of z premultiplying Equation (6.150).
The Z - transform plays an important role 'in the communication
For effective feedforward control, the manipulated variable must engineering .and control engineering. It is used in the analysis and
be able to act on the process output before' the disturbance affects the representation of discrete. time linera shift invariant system. The Z
response through the load transfer function, hence the time delay - transform. provide, a method for the analysis of discrete time
restriction. systems in the frequency domain.
Another way to handle the approximation of Gr is to tune Gf in the 2. What ,is region of 'convergence (ROC)?
field for typical disturbances using adjustble parameters, such as 'done The Z - transform of a discrete sequence is an infinite power series,
for continuous control. In this way, the mathematical intractability of hence the Z - transform exists only for these values of Z forwhich
Equation (6.149) can be avoided. A lead-lag digitalmodelwith time delay the series convergen. If theF (Z) is 'a Z.- t.ransform of f(k) then the
is given. by region of convergence of F (Z) is the set of all values of Z, for which
F(Z) attains a finite values.
We can easily implement a pure time delay in Gr. This is an advantage 4. State final value theorem.
for digital control because a time. delay term can' only be approximated If f(k) is causal and .stable signal and F(Z) exists with Z = 1 included
using lead-lag components for continuous feed forward control. in the ROC then the final value theorem is given by
CONTROL Proce-:s
-----y1
Fig. 7.1: 5150 Controi system 1-----.----....
In industry, many process may have more number of output variables Process
which can be controlled and may have more number of input variables
which can be manipulated. Such a system is called Multi-Input,
Multi-Output (MIMO) process.
Since each controller uses only one measured controled varibale and
adjustes only one manipulated variable, the actions of the controllers are
relatively easy to monitor. A third advantage is that hardware and Xsp
software are readily available for multiloop controllers (using multiple SISO
controllers).
A MIMO process is said to have interaction when process input
,<manipulated) variables affect more than one process output
(controlled) 'variable. Fig. 7.. 5: Blending process
In most industrial MIMO process, each manipulated variable (input
signal) may affect severalcontrolled variables (output signals) causing Loop 1
interaction between the input/output loops. For this reason, control of 1- - - ~ - - _. - - - - - - - - - - - - - - -
I
output flow rate and output composition. Let the outputsbe represented
by Y2 (total flow) 'and Yl (composition), and the inputs by U1 (flow rate of
streamFg) and U2 (flow rate of stream FB).
How should outputs and in~ts be paired? If _we assumethat the
output-input ,pairings are Yl"U1' and Y2-U2, the instrumentation diagram is (b) Y1-U2 and
.
Y2 -u 1 pairing
, ,
shown in Figure 7.6(a)., If we assume that the output-input pairings are Fig. 7 . 6: Blending system alternatively closed loop pairings
Yl"u2and Y2-U h the instrumentation diagram is shown in Figure7.6(b).
(Multi-L{)op Reguiatory Control' 7.5 Computer Control of Processes)
Assume that a control engineer has decided to blend a concentrated :r 3.2 ::r
4.6 0.35 8
4 0.2 2
3.8 0.15-+--...,...--....,....--------1 0
4 6 8 '10 0 2 4 6. 8 10
o 2 ,4 6 8 10 0 2 4 6 8 10
Time,min --+ Time.min ---Ito-
Time,min ---+ Time,min ---+
3.6 ,8 2.5
2
3.4 ~6
:::s ~ 1.5
::r 3.2 4
1
, 3 I--------------~ 2 0.5
o -+-_~----....,---.l-J'_t_-.u. __-J...I_I O-+---,-I-L-.a.,--..1...,-01...........1,---1
2.8 L.-_oI...-_oI...-_'---_L..-..~
o 2 4' 6 8 10 2 4 6 '8 10 o 2 4 6 8 .10 o 2 4. 6 8 10
Time,min - + Time,min --+
Time,min----ll> Time,min ----+
Fig. 7.. 9: Response to diluted ethanorc~ncentrafion setpelnt change
Fig. 7.8: Resoponse to total output flow rate setpoint change: with both loops closed. Lo.op2 pairs Y2: Uland leep 1.pairs Y1- u2
Loop2 has lotal output flow rate as the measured oulpul Y~
and pure water flow _rate as the manipulated input u 2 Loop 1 is open. If you select variable pairing as Yl - U2 and Y2 - ,lib would this system
be stable?
We have seen that the individual control loops have quite good
performance. Since We desire to control both the diluted" ethanol The response .to step increase in 80%- concentrated ethanol flow rate
concentration (Yl) and total output flow' rate (Y2), we need to have both Ul under two conditions:
loops closed simultaneously. The response to a setpoint change in diluted 1. With the flow controller (loop2) open (solid)
ehanol concentration with both loops closed is shown in Figure 7.9.
2,. With the flow controller closedIdashed) are shown in Figure 7.10__
Notice that the closed loop system is unstable. What characteristics
From the response we can understand that the closed loop system is
of this ,control strategy caused the instabilityto occur when both loops stable. -
,are closedr. W.9uld this behaviour have occured ifvatiable pairing of
Yl-Ul and Y2-U2 have been used?
From the responsesshown in Figure 7.9 and Figure 7.10, we can
understand that selection variable pairing YI-U2 and Y2.;.Ul will make the
blending system to be stable whereas selection of variable pairing YI-Ul
and Y2-U2 will make the blending. system to become unstable.
[Multi:.toopRe,gulatory'Control 7.9 ( 7.10 ) ComputerControl of Processes, )
[~xample 7.2
Diluted ethanol concentration Y1 Total output flow rate Y2
Interaction of control Ioops in Stirred T.ank Reactor
0.33
4.6 In the CSTR of Figure 7.11, the temperature is controller by the
/ ' - - Flow loop open
0.32 I - - ~ Flow loop closed flow of coolant in the jacket (Loop 2) while the efficient concentration is ..
I
4.4
I controlled by the inlet feed flowrate (Loop l), Assume that initially both
~ 0.31
~ 4.2
1~---------f effluent concentration and temperature are at their desired values.
0.3 r,
4 ,-------------- FtGAi'
Feed
,0.29
3.8 t---~-....,--__-__..-____J
o 2 4 6 8 10 o 2 4' 6 8 10
Time,min ---+ Time,min~
(/_l..---....
::r
~AI",,----,"---'
3.2 "
Coolant
,
\
3.0 --------_.-._---- Effluent
..
2.8 1-----,.----.-----...--_-~ 0.5 I
o 2 4 6 8 10 o 2 4 6 8 10 Fig. 7.11: Control loops in CSTR
Tirne,min ---...: Time.min ~
1. Consider a change in inlet concentration (load change) or the
desired effluent conceritration (setpoint change). Loop I will
Fig. 7.10: Response of diluted ethanol" concentration to step change:
compensated for these changes by manupulating the feed flow
Loop2 open (solid curve) and loop 2 closed (dashed curve).
rate. However, this change in the feed flow rate also disturbsthe
Hence choice of pairing is crucial since, a bad choice" will make the reactor temperature away from the desired value. Then Loop '2
closed loop system to become unstable, This problem could arise due to attempts to compensate for the change in temperature by varying
interaction' between the different loops. ,the coolant flow rate, which inturn affects the effluent concentration.
Eventually there is a need for a measure that can both give some 2. On the otherhand, an attempt to compensate for changes in feed
advise when solving th~ pairing problem and the also quantifies the le~el temperature (load change) or the desired setpoint of "reactor
of interaction occuring in the .system. , temperature (setpoint change), it 'also causes the effluent
concentration to vary. The loopl attempts to compensate for the
change in effluent concentration by varying the feed flow rate,
which in turn disturbs the reactor temperature.
It is clear form the above that, loopl interacts with loop 2 and also
loop 2 interacts with loopI (i.e., in both directions).
(Multi-Loop Regulatory Control .[ 7.11') ( 7.12 ) Computer Control of Processes )
gIl Gel
1---_ _...1...-. . . . . . Y2 (5) Y1 = Yisp ... (7.5)
1+ gIl Gel
Process
g21 Gel
Yz = 1 + gil Gel
Yls p ... (7.6)
"Loop 2
(b) Closed loop From the above equation, it is clear that any change in the set point
y 1 sp will affect not only controlled outpu y 1 but also the uncontrolled
.fig. 7.1.2: Block diagram of 2 x 2 process output Y2.
rMulti~LoopRegulatory Control ('7.13 ) ( 7.14 ) Computer Control of Processes )
Loop 1
Y1 sp
U1 Direct effect
:: - - -+- - - - ,,-+----+-+ Y1
.... ..... -' "
.... .".
..... ./
/./'"
"',"
./ "
..... -1
I
I
I Indirecteffect :
--------------+-------------
(b) Both loops Closed
The dashed lines in Fi~re 7.13 (a) indicate schematically the routes
through which Ylsp affects the two outputs. Simillar conclusions are drawn if
we consider loop 1 open and loop 2 closed. The situation become more complex
when both loops are. closed.
..~ (7.7)
2. Both Loops Closed
This case is shown-in Figure 7.13 (b), Initially the process is at steady
state with both outputs at their desired values. Consider a change in the Yis) = g 21 (s) Gel [Ylsp - Y1 (s)] + g 22 (s) G e2 [ Y2sp (S) - Y2 (s)] '
setpoint Y1 sp only, and keep the setpoint of loop2 the same (i.e Y2sp = 0).
Y2 (s) +[1 + g22 (s)G c2 ] +g 21 (s) Gel Y1 (s)
g21 (8) Gel Ylsp + gZ2 (S)Ge2 Y2sp (8) ...(7.8)
Solve equation (7.7) and (7.8) with respect to' controlled outputs y 1 and
Y2 and take the following closed loop input-output relationships:
(Multi-Loop RegulatoryControl . ( 7.15)
( 7.16 ) Computer ControlofProc~sses J
1~5 STEADY STATE EFFECTIVE GAIN-llf TWO INPQI AN,D TW,D OUTPUT
,y1 =' Pi l (s) Y1 sp + Pl 2 (s) Y2 sp (7.9)
!
PROCESSES
. Y2 = P21 (S)YI sp +P22 (S)Y2sp ... (7.10)
Let us consider the relationship between Y1 - u l and Y2 - u 2 under a
Where;
number of conditions for a 2 x 2 (two input - two output) process as shown in
= gil Gel + Gel Gc2 (gil g2~ - gl2 g2l) Figure 7.12 (a).' The' open loop input-output relationships are
'Q(s)
... (7.13)
g22 G c2. + Gel Gc2 (gIl g22 - g12 g21) ... (7.14)
Q(s)
"'. (7.11) The corresponding feedback control system, if pairing are' lil - Yland
liz - Y2 is shown in Figure 7.12 (b).
Remarks:
, 1. WIlen Hi.2 = H 2l = 0 (Refer diagram 7.1) in equation (7.7) and (7.8), The open-loop steady state relationship for ,a 2 ?C 2 process is
there is no interaction between the two loops. The closed loop outputs
are given by, the following equations.
H' (7.15)
... (7.12)
a) Steady state effective gain between u 1 and Yl with loop2 closed
Tuning, two control loops without interaction is' very simple, and we that is Y2 = 0 (Refer Figure 7.14).
can tune each loop separately.
2. Equation (7.9) and (7.10) describe the response of outputs Y1 and Y2
when both loops are closed (i.e they have accounted for the interaction
between two loops). ,Tuning two control loops with interaction is very u2 = ...-(7.16)
difficult and we cannot able to tune each loop seperatly.
3. Interaction between control loops is a. significant factor which affects and we can then find the following relationship from equation (7.15)
the "goodness" of a control system. By a suitable pairings 'of
manipulatd variables (u, and u 2) and controlled outputs (Y1 and Y2)'
'we can minimize interaction between two control loops 'and then we
can easily tune the controllers. If the interaction between the two
control loops is more, we have to design special control system which
eliminates the interaction (i.e decoupling the loops).
(MultiLoop Regulatory Control
.@s] Computer Control-of Processes ]
[ 7.17,)
From equation (7.15),
I
Let us ass~e perfect control of YI (loop is closed) that i~ Y1 '=' O.
J
u l + ~2 u 2 ~ 0
I
I Y1 = Kn
I
I
I
I u1 = , ... (7.19)
I
I
f From equation (7.15), we can write
I
I
I y 2 ,= ~l
+ ~2 u 2 UI , ... (7.20)
I-
I Substitute equation (7.19) in equation (7.20)
J
= x, (~ KK12 U2 J+ K 22 U2
I
~~-------------~---------~----3 Y2
11
911 eff
1------.,.. Y2 ~ ~
............. oIIIIII!II ..... _,
"'""'"'"----.
~ _ -.-,;. ~ - _ - - -.. -'- - - -
I
-.--..
Fig. 7.1-5: Block diagram of loop 1 and loop 2 open with u2 =0 ,Fig. 7.16: Block diagram of Y2 -u 2 relcltionship:'with loop 1 closed
( 7.20 ) ,Computer Control of Processes )
[Multi-Loop Regulatory Control ( 7.19)
From equation (7.15), we can write
d) Steady state gain between u 2 and Y2 with loop! open
If the first input u l is constant (~~opJis',open), then output Y2 depend
Y1 = K n u 1 + K12 u 2
only on the second input u 2 as shown inFigure 7.17. K n (0) + K 12 U
2
-. K 12 U2 ... (7.23)
This is the steady state gain relation between u 2 and Y1 with loop 1
is open. -
f) Relative gain between u 2 and Y1 with loop 2 closed ,
u =0 - - . - - - - - - - . . - - -.....
I
I
Fig. 7.17: Block diagram of loop1 and loop2 open with u 1 =0 I
I
!
From equation (7.15)' I
I
i
Y2 = ~1 u 1 + ~2 u 2 - i
!
K 2l (0) + K 22 u 2 I
I
I
~2 u 2 (7~22) i-
r
This is the steady state gain relation between u 2 and Y2 with loopl I
I
open.
I
e) Relative gain between u 2 and Y1,with loop 1 open 912eff
If the first inputu, is constant (loop 1 is open) then output Y1 depend Fig. 7.19: Block diagram of Yl- u2 relationship with I~op 2 closed
only on th~ second input u 2 a~ shown in Figure 7.18.
From equation (7~15), we can write
y2 = K21 U l + K22 U2 =0
Y1 =(
~
KIt K 22 + K 12
K 2I
J
U2
Let us assume perfect control of Y1 (loop 1 is closed) that is Y1 = 0
~-------------------------~----~
I I
: I
1-
: ~
Y1 =- K 12 eff u 2 ..:(7.24)
This is the steady state effective gain relation between u 2 and- y 1
with loop 2 closed.
g) Relative gain between u l and Y2 with loop 2 open -----~--~----~-~----------~-~~~
. I
I
If the second input u 2 is constant f loop 2 is open) then output Y2 921 eft
depend only on the ,first input u l as shown in Figure 7.20.
Fig. 7.21: Blockdiagram of Y2- u t relationship with loop 1 elesed
u 2--
Y2 = ~1 u l + ~2 u 2
~lUI + ~2 (0)
Y2 = ~l u l
This is. the steady state gain relation betweenu1and Y2 with the
loop. 1 open.
(. 7.24 ) Computer Control ofProcess~s] ,
(Multi-loop Regulatory Control. ( 7.23 )
b) Relative gain between input 1 and output 1 for a 2 x.2 process
The relative gain between input land output 1 is
Gain between HI and Y1 with U2 constant (ie loop 2 'is open)
Y2 = .....~l eff .u l Gain between lil and YI with Y2 constant (ie loop 2 is closed)
This is the steady state effective gain relation between u, and Y2
with loop 1 closed. - ( ~ JU2 = constant
7.& RELATIVE GAIN ARRAY-AND SELECTIONOF VARIABLE PAIR lJVITH (
OJI JY2 = constant"
aUI .
MINIMUM INTERACTION
(Refer equation 7.18)
For a 2 x 2process (2 controlled outputs and 2 manipulated inputs) = (Refer equation 7.17)
there are 2ldifferent input-output variable pairs are possible. The four
different input-output variable pairs are u 1 - Yl ' u 2 - Y2 , u 2 - y1and U 2 - Y2: %
Which .one is best? one way to answe.r this question is to consider the
interactions between input - output variable pair of ali 4 combinations
and select one where the interactions are minimaL
1
The Relative Gain Array (RGA) is a heuristic, technique to predict K I I K 22 - K I 2 K 21
possible interactions between input-output variable pairswhen multiple K I I K 22
8180 loops (each loop is constructed using one input and one output
"variable) are used in MI1VIO process. It was first proposed by Ed -Bristol
... (7.28)
and today is a popular tool to determine optimuminput-output pairs of a K I I K 22 - K l 2 K 21
MIMO process. c) Relative gain between input 2 and output 1 for a 2 x'2 process"
a) Definition of the Relative Gain The relative gain between input 2 and output 1 is
The relative gain (Ai) between input j and output i is defined in the Gain betv/ee,n U2 and YI with-u'l constant (ie loop 1 is open)
following fashion: ' Gain between U2 and y 1 with y 2 constant (ie loop 2 is 'closed)
(8ylJU1= constant
Gain between inputj and output i with aU"other loops open
A..
I;
Gain between input j and output! with all 'other loops closed OU2
(~)UK =constantK;tj
A..
IJ
... (7.27) K 12 (Refer equation 7.23)
= (:: )YK = constant'K;t i =
. -[ KllK22~2~12K 21}Refer equation 7.24).
-K12 K 2I
... (7.29)
(Multi.loop Regulatory Control ( 7.26) Computer Control of Processes )
d) Relative gai:il between input, 1 and output 2 for a 2 x 2 process f) 'Relative Gain Array (RGA)
The relative gain between input. 1 and output 2 is Arrange the four relative .gains "'12' A21 an~ "'u' "'22 (for a 2 x 2 process)
into a matrix form, which is known as RGA.
Gain between Ul and y 2,' with u2 co;nstant (ie loop 2 is open)
Gain between Ul and Y2 withy, constant (ie loop 1 is closed)
( Oy2
BUI
J
U2 = constant
'
(~ J Ul = constant
All + 1...12 = 1
0.75 + A,12 =1
A12
0.25 +
+ A22
"'22
=1
=1
1...21 + =1
A22
"'21 +0.75 =1
0.25 0.75 Y2
[ (Refer equation7.21)
KI I
Note:
... (7.31) For a 3 x 3 process, only four of the nine elements would need to be
calculated.
[Multi.Lo()p Regulatory Control [ 7.27 ) ( 7.28 ) Computer Control of Processes)
1 O]Yl Since all the elements in' RGA are equal" the amount of interaction
[0 1 Y2. between the two loops is the same in both configurations of Figure 7.22
and Figure 7.23. We can use two possible different combination of variable
The 1's in the diagonal elements indicate that we can form two non- , pairing.
interacting control loops by coupling u l with YI and u 2 with Y2 as shown in Combination 1 : Yl - u l and Y2 - u 2
figure 7.22.
Combination 2 : Y1 - u 2 and Y2 - ul
The degree of interaction remains the same in both combination.
2x2
I O<An < 0.51 (say Au = 0.25). Then the RGA is
Process Ul u2.
I\, =
[0.25 O.75]Yl
0.75 0.25 Y2
The two larger numbers (i.e 0.75) indicate the recommended variable
Fig. 7.22: 2 X 2 proces with Y1 - u 1 and Y2 - u 2 paring pairing with the smaller amount of interacti?n. Thus we can couple u l
with Y2 and u 2 with Y1 as shown in Figure 7.23.
I
Au = 0 Then the relative gain array is
I O.5<A.n < 1 I(say An = 0.8). Then the RGA is
Ul U2
A = [o.?
0.2]Yl
0.2 O.S Y2 ,
The 1's 'in the off-diagonal elements .indieate that we can form two non- In this case also the. two larger numbers '(ie.O.S) indicate the
interactingloops by coupling u l with Y2 and u 2 with YI as shown in figure 7.23. recommended variable pairing with the smaller amount of interaction.
Thus we can couple ill with Y1 and u 2 with Y2 as shown in figure 7.22. This
is opposite to previous case.
Y1 SP
2x2
AU >1 I
Then A22 = Au > 1 and A12 = A21 =1 - An < 0 (say An = 5).
Process Then the RGA is,
Y2SP
(Multi-Loop Regulatory Control (' 7.29 ) ( 7.30 ) Computer Control of Processes )
In this case we can couple u 1 with Y1and u 2 with Y2as shown in Hence we can pair Yl with .0 1 and Y2 with u 2i n this case because their
figure 7.22. The larger the values of the relative gains above unity, the corresponding gains (0.95, 0.95) are positive and closer to unity.
"larger the "holding back" effect will be. Thus we need larger values for Note:
the controller gains. .
U1 U2
Note:
However, if A 0 .05 ' O.95]Yl
Never forms variable pairing by coupling inputs to outputs with [ 0.95 0.05 Y2
negative relative gains.
In this case, if you couple y1with U 2 and Y2 withu., the corresponding Then we would pair Y1 with u 2 and Y2 with u 1
gains, A. 12 and A21 are negative. Then the interaction will take the controlled
outputs in the opposite direction from that desired by the controlled effort IExample 7.51
and control will be lost altogether. - Consider a blending process with steady state gain matrix as shown
We can summarize all the above observations with the following rule below.
for selecting' the input-output variable pairs.
A...
IJ
Possible pairing
A..
IJ
= 0 Avoid 'pairing yi,with u j How would you choose input - output pairings, for this process.
pr~cess gain matrix isgiven as
A...
IJ
=- I Pair s. and u j
The steady state'
A... < 0 Avoid pairing Yi withuj " K = fK.'11 K12 J = [0.025 -0.0751
IJ
LK 21 K 22 . 1 1 J
A...::;
IJ
0.5- Avoid pairing Yi with u j
The RGA can be calculated from K matrix as
A...
IJ
>-1 Pair v. with u j
IExample 7.4 I
Determine optimum- input-output pair for a 2x2 process whoseRGA
K 11"K22 - K 12 K 21 K 11 K 22 - K 12 K 21
is,
. -K 21 K 12 K 11K22
Ul U2
A = [0.95 ,O.05]Y1'
U1 u2
0.05 0.95 Y2
0.25 0.75]Y1
We can select an optimum input-output variable pairs in such a way [ 0.75 0.25, Y2
that the relative "gains are positive and as close as possible to' unity. \ .-
1 1 0.2 8 +1
Y1(8) = -.-ul8+.
8+1
- ---
0.1 s-r I 0.8 0.58+1
Ul(8)
... (A)
1
Substitute uI(s)= - and find the final steady state value of Y1(s) .as,
s .
, -0.2. ( ) . 0.8 ()
Yz (s).= ---Ul s +--U2 s ... (B)
0.5s+1 s+1. Y1ss(S) ;::: lim(s Yl (s))
s~o
'Step 8:
yes) = -.1_(!J+O
1 8 +1 8 Determine All
By using final value theorem, we can find new steady state value of
~
Y1 as,
( 12=constant _1__ R
;::: limi_1_.~ = 1.25 - O.v
040 (s + 1) (t) 1 .(fJyl J
.au 1 Y2 =constant
Step 4:
Determine A12 > A21 andA22 from A. n
Step'2: Since
Keep yz = constant (ie., yz ;::: 0).. Introduce .a unit step change in ul
'All + A12 1 Au + ~2l 1
(ie., U1 (s) =.; J . The.equation (B) can be written as 0.8 + Al 2 1 0.8 + A21 1
.. "'12 = 0.2 .. A,21 0.2
-0.2 ()0.8 () Similarly
o = -O-.5-5';"'-
'+-lUl 8 +-8-+-1 U2 s
A.12 + 11..22 =1
() 0.2. s + 1 ( ) 0.2 + A.22 =1
Uz 8 = 0.8 0~55+1 Ul ,8
.. A,22 = O~8
( 7.34 ) computer Control of Processes )
[ 7.33 )
The steady state mass b~lance equations of the mixer process are,
step 5:
F =F I + F2 ... (A)
Ul U2
Fx= F 1 Xl + F 2 X 2 ... (B)
The RGA is A= [0.8 0.2]Yl Given data:
0.2 0.8 Y2
F = 200mole/hr
Step 6:
x 0.6
We can draw a conclusion fromRGA as
We can pairy, with u l and Y2 with u 2 to form two loops withminimum Xl = 0.8
interaction. ' x 2 ,= 0.2
Note: with given data and equations (A) & (B)
If you couple differently ie YI with Uz and Y2 with ul' then'the 200 = FI + F2 ... (C)
interaction of the loops would have been four times ~arger (ie., ~:~ = 4)' (200)(0.6) = F (0.08) + F (0.2)
1 2
F = 201 arid
x 0..6012
(:~ 12 =
(201- 200) 1
.. eonstan t
=
(134.4 - 133.4) "1= 1
2. ChangeP'l by one unit (i.e. Fl = 133.4 + 1= ~34.4) while holding Second, we donot wantto pair with a relative gain of 0 because that
, x = 60 % constant. Solve equations (C) and (D) and find means that the particular input does not have an effect on theparticular
F = 201.67 and output, when all of the other loops are open.
Look at Row 3 in relative gain array, this corresponds to output 3,
F 2 = 67.27
we would not pair Ys with Us because of the 0 term. We also would not
(~FJ .
(201.67 - 200) 1.67
-1-= 1.67
pair Yg with u1because of -3 term. This means that Yg must be paired
with. u 2 " Let us indicate this- choice in our RGA.
8 Fl x = constanf = . (134~4 - 133.1) =
Therfore relative gain between F and F 1 is
A= . [0.0.4.6 O
..4]F
0.6 x r{l} 1
A = 3 4
l-3
I
From relative gain array,. the optimum input - output pairing is
{4}
F -F1 and x - F 2
(Multi-Loo.p Regulatory Control [ 7.37,J ( 7.38 ] Computer Control ofProcesses )
NOV-Twe have eliminated outputs Land 3 (row -1 and row 3 as well as Note: RGA for Non-square systems
input's 1 and 2 (column 1 and column 2) we have no choice but to pair Y2 The relative gain array is a square matrix, which implies that the
with Us and we make _that choice below: number of manipulated va~ables,is equal to number of controlled outputs.
Now, suppose that we have a process with two outputs (Y1 and Y2) and
three manipulated inputs (up u 2 and us) How to pair three inputs (up u 2
and us) with two outputs (Y1 and Y2) there are three possible pairs of
manipulated variables (up u 2 ) (u 2 , us) and (us' u.). Therefore, we can from
three .different relative gain arrays:
U1 u2
How would you choose input-output 'pairings for this process? and we need to examine all of the above RGAs arid we can select set of
two variable pairs with minimum interaction.
Looking at columri B, we see only one relative gain that is not
negative, ~herefore Y'e must pair Y2 withu.. IExample 7.10 I
Looking at row 4, we notice that there is only one relativegain that Non - Square Process
, is not negative,therefore, we must pair Y4 with u r
Consider process with 2 measurements and 3 inputs as shown below:
Looking at column 2, there are 2 non negativerelative gains, but we
have already used output ~, so we must pair Y3 with u 2
ry~] =
, Y 0.5 0.07 0.04 U1]
This leaves Y1 and u 4 which fortunately, has a favourable relative
[0.004 -0.003 -0.001] [ ::
gain (0.86).
'Note: How to pair input-output variable pairs for this process.
{ If the Y1 - u 4 relative gain had not been favourable, we would have , Step 1:
been forced to drop output 1 and input 4 and simply have three control
loops for this pr.ocess. Forma three possible pairs of manipulated variables (Up-.u2 ) , (u 2 , u g)
and (U
3,U1).
( Multi-~oop Regulatory Control'
( 7.40)
Step 2:
1.1 MULTI LOOP'PIO CONT.ROLLER
In multiloop PID controller scheme, each controlJoop is controlled
Form a three possible 2' x 2 process as shown below:
by a single PID controller. After selection of input-output variable pairings
[Y2Yl] = [0.5 .0.07] [Ul] with minimum interactions (using RGA) one can design PIDcontrollers
... (A) for individual loops. This multiloopPID controller scheme was the .first
0.004 ....-Q.003 U2
approach used for MIMO control in process' industries.
Y1
[y
,].=. [0.07 ~"O.04] [U 2] ... (B)
Advantages of Multiloop PID controller scheme are,
2 ,-0.003 -0.001 _u3 i) Use of simple Pl Di.control algorithm.
Ii) Ease of understanding by plant operator.
[Yl] . [0.'5 '0.04 _] [U1] ... (C) iii) The actions of the controllers are relatively easy to monitor.
Y2 =0.004 -0.001 u 3
The matrix notation for multiloop PID controller scheme is explained
Step 8:' as shown below.
Determine 'three RGAs such as RGA 1, RGA 2 and RGA 3 for the For a single loop PID control scheme as shown in Figure 7.25.
three 2 x 2 process (A), (B) and-(C) respectively.
D
Ul U2
u2. Ug
=0 ... (7.33)
0.24 0.76 Y2
For a multiloop PID control scheme for a nxn process, the closed
Step 4: loop transfer function Y = (I + Gp Gc)-l G p Gc Ysp (7.34)
For a 2 '5 ~. precess.vas shown in Figure 7.1~. Steps involved in tuning multiloop PID controller scheme using BLT
method are~.
P'],
Y Sp = .,[. Y1S.
Y2sp . Y 2
Y = [Yl.]" G=
p
Step 1:
By ignoring process interactions, each PID controllers are first.
where designed using Ziegler - Nichols tuning method (as shown in Table below).
Step 5: The RGA indicates how the inputs should be paired with outputs to
form loops with the smaller amount of interactions. In some process, even
Determine multi loop .closed loop log modulus
t,hough the interactions between the loops are small, the two control loops
still affect each other's operation very seriously and the overall control
LIDaX =
C
max{2010gl~I}'
W 'l+W system is characterised as unacceptable.
When the designed process is given with two mildly/strongly
Layben suggested that Lrgax_ = max.
w
L, = (2n)
.
dB
interacting loops, he or she introduces in the control system special new
where n is the dimension of the multiloop process.' elements called decouplers. The purpose of d-ecouple is to cancel the
interaction effects between the two loops and fhus rreduce to two
Step 6:
noninteracting. control loops.
If L~ax =2n, stop the iteration and take Kci andT'Ii as calculated in
Let us now study how we, can 'design the decoupler for a 2 ,>.c 2 process
step 3as detuned PI controller parameters. If L~ax* (2n) dB, select a
with two strongly interacting loops as shown in Figure 7.26. In this figure'
newv~.lue of F and return to step 2 until L c max = (2n) dB.
the decoupling matrix is restricted to the form,
7.9 DECOUPLER
D(s) = [d ll (s()=) 1 d 1(2 ()s) .] ... (7.36)
We have seen the problems with multi loop interaction. These occur d 21 S d 22 S = 1
because a manipulated input affects more than one controlled output. One
Since the purpose of the -decoupler D(s),.ls to remove the interaction
approach toha~dling this problem is known as decoupling. The idea is to
in the process as shown in Fig. 7.26. The, interaction between loop 1 and
develop "Synthetic" manipulated inputs that, affect only one process output
loop 2 is zero means g12 = g21 = 0 in the 2 x 2 process. So we can write
.eac1l,.This approach is shown in Fig~ 7.26.
Loop 1
"iIIIIIIIIII __
u*1 I
~_"""'''''' __ ''''
- - - - - - - - - - -- - - - -.
I
... (7.37)
I
1 _
00. (7.39)
Decoupler Process
12.8e-s -18.ge-35
XD(S)] = 16.7s + 1 218+1 [R(S)]
[ XB(s) 6.6e-78 -19.4e-8
0 0
S(s)
10.9s + 1 14.4s+1
where
Y1 = Xn Distillate composition
Y2 = XB Bottom composition
Fig. 7.27: Equivalent block diagram of Fig 7.~6 with complete decoupiing
(i.e., two non interacting loops) u1 =R Reflux flow rate
Remarks:
u2 =S Steam flow rate in Reboiler
Design a deeoupler for this process
i) Two interacting loops are perfectly decoupled only when the
o process transfer function is perfectly known ie gll" g12' g21 and g22 0
Step 1:
are known correctly. Identify steady state gain' matrix for this process
i i) For a non linear and non stationary process, adaptive decouplers
are required. '
x, XB
~K21
/\ = [2.01 -l.OI]R
d 21 -1.01' aoi S
K.22
where KIt' K 12 ," K 21 and ~2 are steady state gains of the process .
..
(Multi.Loop Regulatory Control [ 7.47 )
Step 8: .IExomple:i.'12,
Identify the .input-output pairing using RGA, fromRfkA, the best Consider a following 2 x 2 process
.input - output variable pair is ~ ~ Rand 'XB - s.
Step 4: 1-18 .9 e-3s -12.8 e-s
Steady state decoupler is = l.-19.4e-
21s +J
s
16..78 +.1. [So (s)]
6.6e-7s R(s)
-K12 . 14.48 + 1 10.98 + 1
D 12 =
K 11
Design a decoupler for this process
-(-18.9)
= = 1.477
12.8 Step 1:
-K21 Identify steady state, gain matrix for this process
D 21 =
K22
x, XB
= -(6.6) K = [-is.9 12.8] S
= O~34
-19.4 -19~4 -6.6 R
Step 5:
. Step 2:
The dynamic decoupler is
Identify the RGA matrix
-(-18.9 e~3S) XD XB
21s+1 A _= [~1.01 2.01 ],8
12.8e-s 2.01 -1.01- R
16.7s+1
Step 8:
= (1.477) 16..7-s+ 1 e-2s Identifying input - output, pairing using RGA .
12.8s+1
From RGA, the best input - output variable pair isX, - R and ~ - S
-6.6 e-7s Step 4:
10.9s+1 Steady state decoupler is
, = , -19.4'e-s
.14.48+1
du = -K12 = -0. 677.
K 11
= (0.34) 14.48 + le-6s
10.9s+1
(Multi.L~~p Regulatory Control ( 7.49 )
( 7.50 ) Computer Control of Processes )
Step 5:.
TWO MARKS OUESTIDNS AND ANSWERS
The dynamic decoupler is
1. What are the advantages of RGA?
_ . -g12 (s) i) 'Provides two types .ofuseful information
- -gll (s) a) Measure of process interaction
b) Recommendation about best pairing of manipulated and
-12.8e-s controlled variables.
l-6.7s+1 ii) Requires knowledge of steady state gains but not dynamics of a
=
-18.ge-38 process.
21s+1. 2. What are the problems arising from interactions?
i) Closed loop system may become destabilized.
= (0.677) 21s + 1 e2s ii) Controller tuning becomes more difficult.
16.78+1
3. How do you pair input u with output y for a process?
.. Select a control . loop by pairing the controlled outputs Yi with
-g2l (8)
g22(S) manipulate variables mj in such a way that the relativ~.gains Iij are
positive, and as close as 'possilble to unity.
-19.4 e-8 4. What are properties' of. RGA?
. (-)14.4s+1 The properties of RGA are,
= 6.6e-78 i) Summation of elements in any row = 1
10.9s + 1 i i) Summation of elements in any column =1
iii) RGA is independent of scaling used for input' and output variables.
10.9 s+ 1 68
= (2.94) ----e 5.. Comment about interaction exisiting a system in the followipg .
14.48+1
RGAval~es. .
Thedecoupler elements Dj, and D 21 are physically unrealizable
and 'e6s respectively. This is due to current decoupler
because the term s e 2s ii) A. = [0.5 0.5J
0.5 0.5
output to depend on the' future' inputs.
iii) At = [o.s 0.2] Iv)
0.2 0.8 .
A. ~ [01 11
OJ
i) No interaction (pairing should be Yl - u l and Y2 - u 2) .
ii) Strong interactio existing. in the system (pairing can be either
(Yl - u l an Y2 - u 2 ) or (Y1 - u 2 and Y2 - u l )
iii) Mild interaction exists in the system (pairing should be Yl~ - u l
and Y2 - u 2 ) .
iv) No interaction (pairing should be Y2 - ~1 and Yl -u2 )
8.2 computer Control of Processes]
CHAPTER
Distrubances
MULTIVARIABLE ..
REGULATORY CONTROL
Centralized
1 .. .- ..... oW 11IIo""
control . . . . . _ ... t
8.1 INTRODUCTION
The chemical plant usually consist of Multi-Inputs and Multi-Output Set point
(MIMO) variables,
Fig. 8.2: Centralised control scheme
Distrubances
, \
Decentralized f+
control,
Manipulated 1---"-'" Controlled
variables Process variables
2. Easy to understand by' Does not optimize the use of the error diagonal decoupler process
plant personal. manipulated variables. controller
3. Standard control designs There are many. possible control ... (8.1)
Loop 1
have been developed for configuration.
common' unit operation.
4. Maintaining. a failed loop Does not handle input constraint
doesnot affect the. other efficiently.
loops.
The closed loopTiF is" 'I'hestability of the closed Ioop system is determined by, the poles of
Y(s) = Gp(s) U(s) ... (8.3) the characteristics equaiton. '
Frosn the' above equations ", 'Usually process with equal number of i nputs and outputs are
controlled by methods as shown in Figure 8.4 and Figure 8.5. Such a system
Y(s) = (1+ des) Gc(s~i ,G(s)'G~(s) ~(s) ... (8.5)
is called square system process with unequal number of inputs and outputs
Where, (called norisquar systems) often arise in the process industry. Such
nonsquare systems may have either more outputs than 'inputs or more
Y,' 1(8)] inputs than outputs. Examples of nonsquare systems are mixing tank
yes) [
,Y2 (s) process: 2 x 3 system, shell standard control problem: 5 x 7 system, crude
distillation system: 4 x 5 system etc.
The two simple methods to control non square system are
G(s)
1. Davison's method
2. Tanttu and Lieslehto method-
Davinson's methods of designing centralised PID controller
Davision has proposed a centralised multivariable PID controller
Res)
'[r (s)]'
1 tuning method for square systems. Here proportional Ingegral and
r2 (s) derivative gain matrixare given by
Kn [G(s= 0)]-1
where [G(s = 0)] is called the rough tuning matrix' and 8 and E are
the fine tuning parameters, which generally range from 0 to 1. This method
can also be extended to non-square systems. As inverse not exist for non
square system, Moore-pe.nrose pseudo inverse is used. For matr-ix A,
Moore-penrose pseudo inverse is
A+= A H (A x AH)-l
where A H is the Hermitian matrix of A.
Full matrix Process So for a non-square system, PIn controller gains are
controller Kc = 8[G (s= 0)]+
Loop 2 K1 = c;[G (8 = 0)]+
Fig. 8.5: 2 X 2,Centralised controller coti1prising ,'ullma'trix controller Kn = [G(s) = 0]+ ... (8.8)
(MultivariableRegulatory Control 8.7 8.8 computer Control of Processes]
8.3 MODEL PREDICTIVE CONTROL (MPCJ However, the fundamental framework of MPC algorrthma is in
common for any kinds of MPC schemes. The basic elements of MPC are
MPC .is a multi'variable. controller but also a wide class of optimal
illustrated in Figure 8.. 6 and can be defined as follows: ' e
:
y
Model prediction
1
p
~:... Prediction Set point
current horizon
step
.................. max
I
.
u
: y(k+i)
... "' .IiI_Willllllll.fIII_
~
I
N-l
Well understood tuning parameters .
~Uk-i+l + SN Uk--N+1
"C
Yk+1 LSi + dk +1 ... (8.16)
- Prediction horizon i=l .
Development time much shorter than for competing advanced Y~+l = Sl d U k + L S i dUk:"'i+1 + SN Uk-N+l + ak+~.
i=2 ." ~
control methods.
Easier to maintain: changing model or specs does not require So, for the jth step into the future, we find
complete redesign, sometimes can be done on the fly. "C "
Yk+j = Yk+j +dk +j
8.4 DYNMIC MATRIX CONTROLLER IOMC) N-j N-l .
DMC is amultivariable controller which .was developed by' Shell Oil L S i L\Uk_i+j + L s, L\Uk_i+j + s~ Uk-N+j+ dk +i
Company in the 1960s and 1.970s. It is based on a step response model, ~ ~=j+1 . 'Co~rc;;;term
effect of future effect of pasteontrol moves
which has the form control moves .
'Yk = 8 1 ~Uk-1 + 82 AU k- 2 +.... + SN-1 LiUk-N+l+ 8N Uk-N ... (8.12) ... (8.17)
and we can separate the effects of past and future control moves
which is written in the form
"C . ' . ( effect of current and
N-l
Yk+j = 8 1 L\uk+j-l + 8 2 L\u k +j - 2 + ... + SJo L\uk ~.
Yk = LSi~Uk-i + SNuk-N ... (8.13) 1future moves
i=l
where. Ykis the model prediction at time step k, and Uk-N is the
+'SNUk-N+j + Sj+l tiuk_1 + Sj+2. LlU k _2.
. + +S
...
AU
N-l Ll k-N+j+1
I
.
J effect of past moves
manipulated input N steps in the past.
+dk +j {correction term
Note that the model-predicted output is unlikely to .be equal to the
actual ..measured'output at .time step k. The difference between the ... (8.18)
~ measured output (Yk) and the model prediction is called the additive
The most common assumption is that the correction term is constant
disturbance. in the future (this is the "constant additive disturbance assumption"):
... (8.14)
... (8.19)
The "corrected prediction" is then equal to the actual measured
output at step k, Also, realize that there are no control 'moves beyound the control
horizon of M steps, so '
"C
Yk =Yk + d k .. (8.15) L\Uk+p-1 =. 0
Similarly, the corrected predicted output .at the first time step' in
the future. can be found from,
"C
Yk+l = "
Yk+1 + d" k+l
(Multivariable Regulatory Control [ 8.13 J computer Control of Processes' )
In matrix-vector form, a prediction horizon ofP steps and a control In equation. (8~21) the correct-predicted output response is naturally
horizon of M.steps, yields composed ora "forced response" (contributions of the current andfuture
control moves) and a "free response" (the output changes that are;
predicted if there are no future control moves). The difference between
the setpoint trajectory, r, and the future predictions is
.. c 81 0 0 o o
Yk+l
.. c
Yk+2'
82 81 O o 0' r _ y C= ... (8.22)
"---v--'
corrected
= S predicted
.. c J Sj-l Sj-2 error.E?
Yk+j
Uk-N+P d k +P
~ '---v---' ... (8.25)
Pxl past inputs, up Pxl predicted
disturbances.d
... (8.20) .
and
which we write using matrix-vector notation
Notice that the current and future control .move vector (AUf) is
proportional to the unforced error vector (E). That is, a controller gain
matrix,K, multiplies the unforced error vector (the future errors that
would occur if there were ~o control move changes implemented).
Because only the current control move is actually _implemented, we
use the first row of the K matrix, and
AUk = K1E ... (8.31)
where K 1 represents the first row of the K matrix,
where K =(STf Sf + W)-l STf
Perhaps it is worth summarizing the steps involved in implementing
DMC on a process.
1. Develop a discrete ste:presponse model with length NC8.13), based
on a. sample time, ~t.
2. Specify the prediction CP) and control (M) horizons, N ~ P ~ M.
I'