19.07.2014 Views

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

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

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

ANALYSIS OF MULTIYARIABLE SYSTEMS 587<br />

TABLE MA<br />

Doyle-Stein criterion <strong>for</strong> Wood and Berry column<br />

REAL KP(4,4),KC(4)<br />

DIMENSION RESET(4),TAU(4,4,4),D(4,4),WK(4),WKEIG(12),SVAL(4)<br />

COMPLEX GM(4,4),B(4,4)&(4,4),QIN(4,4),YID(4,4)<br />

COMPLEX A(4,4),ASTAR(4,4),H(4,4),CEIGEN(4),WA(24),VECT(4,4)<br />

C PROCESS TRANSFER FUNCTION DATA<br />

DATA KP(1,1),KP(1,2),KP(2,1),KP(2,2)/12.8,-18.9,6.6,-19.4/<br />

DATA TAU(1,1,1),TAU(1,1,2),TAU(1,2,1),TAU(1,2,2)/4~0./<br />

DATA TAU(2,1,1),TAU(2,1,2),TAU(2,2,1),TAU(2,2,2)<br />

+ /16.7,21.,10.9,14.4/<br />

DO 1 Id,2<br />

DO 1 J=1,2<br />

TAU(B.I,J‘)=O.<br />

1 TAUi4,i,Jj=O.<br />

DATA D(l,l),D(1,2),D(2,1),D(2,2)/1.,3.,7.,3./<br />

C EMPIRICAL CONTROLLER SETTINGS<br />

DATA KC(l),KC(2)/0.2,-0.04/<br />

DATA RESET(l),RESET(2)/4.44,2.67/<br />

C SET INITIAL FREQUENCY AND NUMBER OF POINTS PER DECADE<br />

W=O.Ol<br />

DW=(lO.)w(1./20.)<br />

N=2<br />

IA=4<br />

M=N<br />

IB=4<br />

IJOB=O<br />

C MAIN LOOP FOR EACH FREQUENCY<br />

WRITE(6,lO)<br />

10 FORMAT(’ FREQUENCY SINGULAR VALUES’)<br />

100 CALL PROCTF(GM,W,N,KP,TAU,D)<br />

CALL FEEDBK(ti,W,k,Ik,RESET)<br />

CALL MMULT(GM,B,Q,N)<br />

CALL IDENT(QIN,N)<br />

CALL LEQ2C(Q,N,IA,QIN,M,IB,IJOB,WA,WK,IER)<br />

CALL IDENT(YID,N)<br />

CALL MADD(YID,QIN,A,N)<br />

CALL CONJT(A,ASTAR,N)<br />

CALL MMULT(ASTAR,A,H,N)<br />

IZ=4<br />

JOBN=O<br />

CALL EIGCC(H,N,IA,JOBN,CEIGEN,VECT,IZ,WKEIG,IER)<br />

DO 6 I=l,N<br />

6 SVAL(I)=SQRT(REAL(CEIGEN(I)))<br />

WRITE(6,7)W, (20.*ALOGlO(SVAL(I)),I=l,N)<br />

7 FORMAT(lX,4F10.5)<br />

W=W*DW<br />

IF(W.LT.l.) GO TO 100<br />

STOP<br />

END<br />

FREQUENCY SINGULAR VALUES<br />

.OlOOO -.00281 -.02777<br />

.01122 -.00410 -.03454<br />

.01259 -.00573 -.04311<br />

.01413 -.00775 -.05402<br />

.01585 -.01024 -.06790<br />

.01778 -.01323 -.08560

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!