[Luyben] Process Mod.. - Student subdomain for University of Bath

DESIGN Ok CONTROLLERS FOR MULTIVARIABLE PROCESSES 607
TABLE 17.3
TLC for Ogunnaike and Ray column
REAL KP(4,4),KC(4j,KL(4)
DIMENSION RESET(4),TAU(4,4,4),D(4,4),TAUL(4,4),DL(4)
COMPLEX GM(4,4),B(4,4),Q(4,4),YID(4,4),CHAR(4,4),A(4,4),C(4,4)
COMPLEX GL(4)
C GM PROCESS TRANSFER FUNCTION DATA
DATA KP(l,l).KP(l,2),KP(1,3)/0.66,-0.61,.0.0049/
DATA KP(2,1),KP(2,2),KP(2,3)/1.11,-2.36,.0.012/
DATA KP(3,1),KP(3,2),KP(3,3)/-34.68,46.2,0.87/
DATA TAU(1,1,1),TAU(l,1,2),TAU(l,l,3)/3+0./
DATA TAU(l,2,l),TAU(l,2,2),TAU(1,2,3)/3+0./
DATA TAU(1,3,1),TAU(1,3,2),‘~AU(l,3,3)/2+0., 11.61/
DATA TAU(2,l,l),TALJ(2,1,2),TAU(2,1,3)/6.7,8.64,9.06/
DATA TAU(2,2,l),TAU(2,2,2),TAU(2,2,3)/3.25,5.0,7.09/
DATA TAU(2,3,1),TAU(2,3,2),‘TAU(2,3,3)/8.15,10.9,3.89/
DATA TAU(3,1,l),TAU(3,1,2),TAU(3,1,3)/3*0./
DATA TAU(3,2,1),TAU(3,2,2),TAU(3,2,3)/3+0./
DP.TA TAU(3,3,l),TAU(3,3,2),TAU(3,3,3)/2*0.,18.8/
DO 1 1=1,3
DO 1 J=1,3
1 TAU(4,I,J)=O.
DATA D(l,l),D(1,2),D(1,3)/2.6,3.5,1./
DATA D(2,1),D(2,2),D(2,3)/6.5,3.,1.2/
DATA D(3,1),D(3,2),0(3,3)/9.2,9.4,1./
C GL PROCESS TR.ANSFER. FUNCTION DATA
DATA KL(l),KL(2),KL(3)/0.14,0.53,-11.54/
DATA TAUL(l,l),TAUL(1,2).TAUL(1,3)/3*0./
DATA TAUL(2,1),TAUL(2,2),TAUL(2,3)/6.2,6.9,7.01/
DO 2 I=1,3
TAUL(3,1)=0.
2 TAUL(4,1)=0.
DATA DL(l),DL(2),DL(3)/12.,10.5,0.6/
C CONTROLLER DATA - EMPIRICAL SETTINGS
DATA KC(l),KC(2),KC(3)/1.2,-0.15,0.6/
DATA RESET(l),RESET(2),RESET(3)/5.,10.,4./
N=3
CALL TLC(N,KP,TAU,D,KL,TAUL,DL,KC,RESET)
STOP
END
C
SUBROUTINE TLC(N,KP,TAU,D,KL,TAUL,DL,KC,RESET)
REAL KP(4,4),KC(4),KL(4)
DIMENSION RESET(4),TAU(4,4,4),D(4,4),TAUL(4,4),DL(4)
COMPLEX GM(4,4),B(4,4),Q(4,4),HIN(4,4),H(4,4)
COMPLEX GL(4),WA(24),XTLC(4)
DIMENSION DB(4),WK(4)
W=.Ol
DW=(lO.)**(l./lO.)
NN=N
IA=4
M=N
IB=4
IJOB=O
WRITE(6,77)
7 7 FORMAT(’ W TLC1 TLC2 TCL3’)
100 CALL PROCTF(GM,W,N,KP,TAU,D)
CALL FEEDBK(B,W,N,KC,RESET)
CALL MMULT(GM,B,Q,N)