586 MULTIVARIABLE PROCESSES ‘) log w FIGURE 16.9 SISO closcdloop plots. If we plotted just the denominator <strong>of</strong> the right-hand side <strong>of</strong> Eq. 16.42, the curve would be the reciprocal <strong>of</strong> the closedloop servo transfer function and would look something like that shown in Fig. 16.9b. The lower the dip in the curve, the lower the damping coefficient and the less robust the system. B. MATRIX MULTIVARIABLE SYSTEMS. For multivariable systems, the Doyle- Stein criterion <strong>for</strong> robustness is very similar to the reciprocal plot discussed above. The minimum singular value <strong>of</strong> the matrix given in Eq. (16.43) is plotted as a function <strong>of</strong> frequency o. This gives a measure <strong>of</strong> the robustness <strong>of</strong> a closedloop multivariable system. Doyle-Stein criterion: minimum singular value <strong>of</strong> u + &,$ (16.43) The Q matrix is defined as be<strong>for</strong>e Quo) = GM,,-, &,,, . The tuning and/or structure <strong>of</strong> the feedback controller matrix B~ic, is changed until the minimum dip in the curve is something reasonable. Doyle and Stein gave no definite recommendations, but a value <strong>of</strong> about - 12 dB seems to give good results. The procedure is summarized in the program given in Table 16.4 which calculates the minimum singular value <strong>of</strong> the B + &$J matrix <strong>for</strong> the Wood and Berry column with the empirical controller settings. Figure 16.10 gives plots <strong>of</strong> the singular values as a function <strong>of</strong> frequency. The lowest dip occurs at 0.23 radians per minute and is about - 10.4 dB.
ANALYSIS OF MULTIYARIABLE SYSTEMS 587 TABLE MA Doyle-Stein criterion <strong>for</strong> Wood and Berry column REAL KP(4,4),KC(4) DIMENSION RESET(4),TAU(4,4,4),D(4,4),WK(4),WKEIG(12),SVAL(4) COMPLEX GM(4,4),B(4,4)&(4,4),QIN(4,4),YID(4,4) COMPLEX A(4,4),ASTAR(4,4),H(4,4),CEIGEN(4),WA(24),VECT(4,4) C PROCESS TRANSFER FUNCTION DATA DATA KP(1,1),KP(1,2),KP(2,1),KP(2,2)/12.8,-18.9,6.6,-19.4/ DATA TAU(1,1,1),TAU(1,1,2),TAU(1,2,1),TAU(1,2,2)/4~0./ DATA TAU(2,1,1),TAU(2,1,2),TAU(2,2,1),TAU(2,2,2) + /16.7,21.,10.9,14.4/ DO 1 Id,2 DO 1 J=1,2 TAU(B.I,J‘)=O. 1 TAUi4,i,Jj=O. DATA D(l,l),D(1,2),D(2,1),D(2,2)/1.,3.,7.,3./ C EMPIRICAL CONTROLLER SETTINGS DATA KC(l),KC(2)/0.2,-0.04/ DATA RESET(l),RESET(2)/4.44,2.67/ C SET INITIAL FREQUENCY AND NUMBER OF POINTS PER DECADE W=O.Ol DW=(lO.)w(1./20.) N=2 IA=4 M=N IB=4 IJOB=O C MAIN LOOP FOR EACH FREQUENCY WRITE(6,lO) 10 FORMAT(’ FREQUENCY SINGULAR VALUES’) 100 CALL PROCTF(GM,W,N,KP,TAU,D) CALL FEEDBK(ti,W,k,Ik,RESET) CALL MMULT(GM,B,Q,N) CALL IDENT(QIN,N) CALL LEQ2C(Q,N,IA,QIN,M,IB,IJOB,WA,WK,IER) CALL IDENT(YID,N) CALL MADD(YID,QIN,A,N) CALL CONJT(A,ASTAR,N) CALL MMULT(ASTAR,A,H,N) IZ=4 JOBN=O CALL EIGCC(H,N,IA,JOBN,CEIGEN,VECT,IZ,WKEIG,IER) DO 6 I=l,N 6 SVAL(I)=SQRT(REAL(CEIGEN(I))) WRITE(6,7)W, (20.*ALOGlO(SVAL(I)),I=l,N) 7 FORMAT(lX,4F10.5) W=W*DW IF(W.LT.l.) GO TO 100 STOP END FREQUENCY SINGULAR VALUES .OlOOO -.00281 -.02777 .01122 -.00410 -.03454 .01259 -.00573 -.04311 .01413 -.00775 -.05402 .01585 -.01024 -.06790 .01778 -.01323 -.08560
- Page 1 and 2:
I WILLIAM WIlllAM 1. LUYBEN PROCESS
- Page 3 and 4:
The Series Bailey and OUii: Biochem
- Page 5 and 6:
PROCESS MODELING, SIMULATION, AND C
- Page 7 and 8:
ABOUT THE AUTHOR William L. Luyben
- Page 9 and 10:
CONTENTS Preface Xxi 1 1.1 1.2 1.3
- Page 11 and 12:
CONTENTS .. . xlu 5.7 Batch Reactor
- Page 13 and 14:
CONTENTS xv 10.3 10.4 Performance S
- Page 15 and 16:
comwrs xvii 151.3 Transpose of a Ma
- Page 17 and 18:
CONTENTS Xix 20.4 Sampled-Data Cont
- Page 19 and 20:
. Xxii PREFACE Another significant
- Page 21 and 22:
PROCESS MODELING, SIMULATION, AND C
- Page 23 and 24:
2 PROCESS MODELING, SIMULATION. AND
- Page 25 and 26:
4 PROCESS MODELING, SIMULATION, AND
- Page 27 and 28:
6 PROCESS MODELING. SIMULATION, AND
- Page 29 and 30:
8 PROCESS MODELING, SIMULATION, AND
- Page 31 and 32:
10 PROCESS MODELING, SIMULATION, AN
- Page 33 and 34:
12 PROCESS MODELING, SIMULATION, AN
- Page 35 and 36:
CHAPTER 2 FUNDAMENTALS 2.1 INTRODUC
- Page 37 and 38:
FUNDAMENTALS 17 big, complex system
- Page 39 and 40:
FUNDAMENTALS 19 z=o I z + dz z=L FI
- Page 41 and 42:
FUNDAMENTALS 21 The minus sign come
- Page 43 and 44:
FUNDAMENTALS 23 Molar flow of A lea
- Page 45 and 46:
FUNDAMENTALS 25 For liquids the PP
- Page 47 and 48:
FUNDAMENTALS 27 Flow of energy (ent
- Page 49 and 50:
FUNDAMENTALS 29 The units of this f
- Page 51 and 52:
FUNDAMENTALS 31 FIGURE 2.9 Laminar
- Page 53 and 54:
FUNDAMENTALS 33 Then Eq. (2.48) bec
- Page 55 and 56:
FUNDAMENTALS 35 Wenzel in Chemical
- Page 57 and 58:
FUNDAMENTALS .37 This exponential t
- Page 59 and 60:
FUNDAMENTALS 39 a direction 0 = 45
- Page 61 and 62:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 63 and 64:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 65 and 66:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 67 and 68:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 69 and 70:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 71 and 72:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 73 and 74:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 75 and 76:
EXAMPLES OF MA~EMATICAL MODELS OF C
- Page 77 and 78:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 79 and 80:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 81 and 82:
EXAMPLES OF MATHEMATICAL MODELS OE
- Page 83 and 84:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 85 and 86:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 87 and 88:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 89 and 90:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 91 and 92:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 93 and 94:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 95 and 96:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 97 and 98:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 99 and 100:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 101 and 102:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 103 and 104:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 105 and 106:
EXAMPLES OF MATHEMATICAL MODELS OF
- Page 107 and 108:
88 COMPUTER SIMULATION Our discussi
- Page 109 and 110:
90 COMPUTER SIMULATION lation packa
- Page 111 and 112:
92 COMPUTER SIMULATION For this rea
- Page 113 and 114:
94 COMPUTER SIMULATION TABLE 4.1 Ex
- Page 115 and 116:
96 COMPUTER SIMULATION Clearly, the
- Page 117 and 118:
98 COMPUTER SIMULATION TABLE 4.2 Ex
- Page 119:
I()() COMPUTER SIMULATION Settingfl
- Page 589 and 590: ANALYSIS OF MULTNARIABLE SYSTEMS 57
- Page 591 and 592: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 593 and 594: ANALYSIS OF MULTIYARIABLE SYSTEMS 5
- Page 595 and 596: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 597 and 598: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 599 and 600: FIGURE 16.5 Multivariable INA (- 1,
- Page 601 and 602: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 603: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 607 and 608: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 609 and 610: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 611 and 612: ANALYSIS OF MULTIVARIABLE SYSTEMS 5
- Page 613 and 614: DESIGN OF CONTROLLERS FOR MULTIVARI
- Page 615 and 616: DESIGN OF CONTROLLERS FOR MULTIVARI
- Page 617 and 618: DESIGN OF CONTROLLERS FOR MULTIVARI
- Page 619 and 620: DESIGN OF CONTROLLERS FOR MULTIVARI
- Page 621 and 622: DESIGN OF CONTROLLERS FOR MULTIVARI
- Page 623 and 624: DESIGN OF CONTROLLERS FOR MULTIVARI
- Page 625 and 626: DESIGN Ok CONTROLLERS FOR MULTIVARI
- Page 627 and 628: DESIGN OF CONTROLLERS FOR MULTIVARI
- Page 629 and 630: DESIGN OF CONTROLLERS FOR MULTIVARI