Toevallige afwijkingen bij het meten van eigenschappen van random signalen van Heck, J.G.A.M.

Transcription

1 Toevallige afwijkingen bij het meten van eigenschappen van random signalen van Heck, J.G.A.M. Gepubliceerd: 01/01/1982 Document Version Uitgevers PDF, ook bekend als Version of Record Please check the document version of this publication: A submitted manuscript is the author's version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website. The final author version and the galley proof are versions of the publication after peer review. The final published version features the final layout of the paper including the volume, issue and page numbers. Link to publication itation for published version (APA): van Heck, J. G. A. M. (1982). Toevallige afwijkingen bij het meten van eigenschappen van random signalen. (DT rapporten; Vol ). Eindhoven: Technische Hogeschool Eindhoven. General rights opyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal? Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Download date: 18. Nov. 2017

2 i

3

4

5

6

7 4 n

8 - I- b

9 A T r 8. ' P

10 Y o a % P

11

12

13 b W ".

14 k BI c.

15 x u...- e-%

16 9 a U

17 3 Q

18

19 P

20 P

21

22

23 I I I I I I l l I I I I I ,000 2,000 5, Number of Averages, nd I

24

25

26 0.20 y 0.10 g L- w E 0.05 m [r ,02 i c

27 u, Em-

28 t i 6

29

30 1 A B

31

32 IS ANALYZER ONNETED AND HP-IB SWITH ON ADRESSABLE? njij : IS BASI LOADED? R OK? BASI LOADED IS BASI PROGRAM LOADED? IB B

33 REALL ADDITIONAL OHERENE ;FUNTION FROM' WP5420A TAPE LOATION 90,95 OR 99% ONF INTERVAL? 9 E SOLUTION?

34 NO ALGORITHM IMPLEMENTED FOR PEAK AVERAGING

35

36 0 I lll i A,:.

37 P8-J

38 ONFIDENE INTERVALS FOR REASUREMENTS DIGITAL SIGNAL ANALY LOGIAL YESNO,NVA$A REAL MAXEPS,NSIG,MAXI,MIMI,DATIJ1,DATIJ2 INTEGER*4 ONFI4 INTEGER 1,12,SIZE,FUNOD,AVGS,PEAK,LEN,D1(9), *BUFFER(40),Q,OMPL DIMENSION D2 (5),HEAD1 (16),HEAD2(16>,DAT(512),EPSR(512), *OH (51 2) ALL INIT ********** REALL FUNTION FROR ANALYZER ALL AOM( '1T',3) ALL RLNTP('AUTO,TRANS,OHER,ROSS',22) ALL SRQ(Q) IF (Q.NE.69) GO TO 90 ALL AOH('RS',2) ALL WA IT ( ALL S501SA(HEADI,DAT,Q) IF (Q.NE.96) GO TO 90 ALL RHEAD(HEADl,SIZE,Dl (l),dl(2),d1(3),d1(4), *ORPL,D1(6),DI (7),PEAK,D1(8),AVGS,D2(1), *D2 (2),O2 (3),D2 (4),Dl (9),D2 (5) 1 FUNOD=DI (8) IF (FUMOD.LT.7.0R.FUNOD.GT~ll) GO TO 91 IF (PEAK.NE.0) GO TO 92 I2 = SIZE/2 WRITE (1,99995) I2 FORNAT('SIZE=',I6) WRITE (1,99996) FUNOD FORP!AT('FUNOD=',I6) WRITE( ) OMPL FO RWAT (' ORPL= ', 16) IF (FUMOD.NE.9.AND.FUNOD.NE.10) GO TO 10 ********** REALL ADDITIONAL OHERENE FUNTION ALL AOM('2T',3) ALL RLMTP('ADDITI0MAL OHER',l6) ALL SRQ(Q) IF (Q.ME.69) GO TO 90 ALL AOM('RS',2) ALL WAIT(2000) ALL S501SA(HEAD2,OH,Q) IF (Q.NE.96) GO TO 90 DO 9 I=1,256 IF (OH(I).LT.O.O) WRITE (1,99993) OH(I),I FORMAT ('NEGATIVE OHERENE ',F10.4,' ON DATAPOINT ',I41 IF (OH( I). LT.O.O) OH ( I) =ABS (OH (I) 1 ONTINUE ONTINUE ********** DETERMINE NUMER OF SIGMAS 20 ALL MSTR$A('90,95 OR 99% ONF INTERVAL',26, *BUFFER,80) LEN = 26 ALL QUEST(BUFFER,LEN)

39 IF (NVA$A (1,BU FFER,LEN,ONFI4) 1 GO TO 30 ALL TNOU('ANSWER 90,95 OR 99',18) GO TO 20 ONTINUE IF (INTS(ONFI4).EQ.90) frfsig=l.645 IF (INTS(ONF14).EQD95) NSIG=1.960 IF (INTS(ONF14).EQD99) NSIG=2.576 WRITE (1,99998) MSIG FORRAT ( I NSIG=', F1 O. 2) ********** BENDATgPIERSOL ALGORITHNS FOR RANDOM ERRORS IF (FUNOD.NE.7.AND.FUNOD.NE.8) GO TO 40 ********** FIRST Of? SEOND AUTO-POWER SPETRUM DO 35 1=1,256 EPSR(I)=I.O/SQRT(FLOAT(AVGS)) ONTINUE GO TO 70 O WTI MU E IF (FUNOD.NE.9) GO TO 50 ********** ROSS POUER SPETRUM DO 45 I=1,256 IF(OH(1).EQmO.O) OH(I)=O~OOOOI EPSR (2*1-1)= 1.O/ (SQRT ( FLOAT (AVGS)*OH (I) EPSR(Z*I) = EPSR(2*1-1) 45 ONTINUE GO TO ONTINUE IF (FUNOD.NE.10) GO TO 60 ********** TRANSFER FUNTION DO 55 1=1,256 IF (OH(1).EQ.O.O) OH(I)= IF (OH(I).GT.I.O) OH(I)=l.O EPSR (2* = SQRT ( 1-OH ( I) 1 / (SQRT (OH( I) *SQRT (FLOAT * (2*AVGS))) EPSR(Z*I) = EPSR(Z*I-I) 55 ONTINUE GO TO ONTINUE IF (FUNOD.NE.11) GO TO 91 ********** OHERENE FUNTION DO 65 1=1,256 IF (DAT(I).EQ.O.O) DAT(I)=O.OOQOI EPSR ( I) =SQRT (2.l) *(1.O-DAT ( I) 1 / * (SQRT( DAT ( I) 1 *SQRT ( FLOAT (AVGS) ONTINUE 70 ONTINUE ********** ALGORITHM ONTROL DI (1)=0 DO 75 I=1,12

40 IF(EPSR(I).GT.O.Z) DI (l)=di (1)+1 75 ONTINUE D2 (1 1 =FLOAT (DI (1 1) /FLOAT( 12>*100.0 IF (DI (l).ne.o) WRITE (1,99997) DI (1),D2(1) FORMAT('ALGOR1THM ONLY VALID FOR EPSR < 0.2',/, *14,'DATAPOINTS WITH EPSR > 0.2 (',F5.1,'%>') MAXI=DAT (I I* (I.O+NSIG*EPSR ( I 1) MINI=MAX I IF (YESNO('1 TRAE SOLUTION',Ió)) GO TO 80 DO 72 1=1,12,2 DAT(I)=DAT( I>*(l.O+NSIG*EPSR(I) ********** 2 TRAE SOLUTION ON HANNEL A AND B DAT(I+I)=DAT(I+I)*(I.O+NSIG*EPSR(I+1)) IF (OPIPL) 93,78,79 DAT1 J I=DAT (I) DAT IJ 2= DAT ( I+1) GO TO 71 DATIJI =SQRT(DAT (1)**2+DAT( I+1)**2> DATI J 2=DATI J 1 ONTINUE IF (DATIJI.LT.HINI)MINI=DATIJ1 IF (DATIJ2.LT.MINI)MINI=DATIJ2 IF (DATIJI.GT.MAXI)MAXI=DATIJI IF (DATIJ2.GT.MAXI)MAXI=DATIJ2 ONTINUE ALL AOM('IT',3) ALL UAIT(1000) ALL SS01 RA (HEAD1,DAT,Q) IF (Q.NE.112) GO TO 90 DO 74 1=1,12,2 DAT ( I) =DAT ( I) * (1,O-NSIG*EPSR ( I) 1 (1.O+NSIG*EPSR ( I) ) DAT(I+~)=DAT(I+I)*(~.O-NSIG*EPSR(I+I~)/(~.~+~SIG*EPSR~I+~~~ IF (OMPL) 93,76,77 DATIJ 1 =DAT ( I) DATIJ2=DAT(I+I) GO TO 73 DATI J 1 =SQ RT (DAT ( I) **2+ DAT ( I+1) **2) DATI J 2=DATI J1 ONTINUE IF (DATIJI.LT.MINI)MINI=DATIJ1 IF (DATIJ2.LT.NIMI)MINI=DATIJ2 IF (DATI J 1.GT. MAX I) FlAX I=DATI J 1 IF (DAT1 52.GT. MAX I)HAXI=DATI 52 ONTI NU E ALL WAIT(1000) ALL AOM('ZT',3) ALL WAIT(1000) ALL S501 RA (HEAD1,DAT,Q) IF (Q.NE.112) GO TO 90 GO TO 85 ONTINUE ********i* 1 TRAE SOLUTION ON HANNEL B DO 82 I=1,12,4 DAT ( I) =DAT ( 1) * (1.O+MSIG*EPSR ( I) 1 IF (OIVIPL) 93,83,84

41 c DAT (I+1)=DAT (I+1)*(1.O-NSIG*EPSR (I+I)) DAT (I+2) =DAT (I+2>*(1.O+NSIG*EPSR (I+2) 1 GO TO DAT(I+1)=DAT(I+1)*(1.O+NSIG*EPSR(I+I~) DAT(I+2)=DAT(I+2)*(1.O-NSIG*EPSR(I+2)) 86 DAT(I+3)=DAT(I+3)*(1.O-NSIG*EPSR(I+3)) DO 87 J=û,3,2 IF (OMPL) 93,88,89 8% DATIJI=DAT(I+J) DATIJ2=DAT(I+J+I) IF (DATI J 1.LT.MINI 1 HINI=DATI J 1 IF (DATIJ2.LT.MINI)HINI=DATIJ2 IF (DATI J1.GT.MAX I)MAXI=DATIJl IF (DATIJ2.GT.MAXI)MAXI=DATIJ2 GO TO 81 DATIJI=SQRT(DAT(I+J)**2+DAT(I+J+1)**2) I F (DATI J 1. GT. MAX I) MAX I=DATI J 1 IF (DATIJ 1.LT. NINI) MINI=DATIJ 1 ONTINUE ONTINUE ONTINUE ALL AOM('2T'J) ALL WAIT( 1000) ALL S501RA(HEADI,DAT,Q) IF (Q.NE.112) GO TO 90 ONTI HU E ********** Y EXPAND ON BOTH TRAES WRITE(l,99999> RIN1,MAXI FORHAT( '~INIMUR VALUE=',E10.3,/, El O. 31 ALL AOM('3FM1,3) ALL AOM('OT',3) ALL WAIT(500) ALL AOH('RE',2) ALL WAIT(5OO) IF (OMPL.EQ.l.O> ALL AOM('MG',2) ALL WA IT( 500) ALL XYEXP('Y',MINI,WAXI) GO TO 99 * MAX IHU PI VALUE=, ********** ERROR HESSÄGES 90 ONTINUE ALL TNOU('WR0NG STATUS SERVIE REQUEST',28) GO TO ONTINUE ALL TNOU('N0 ALGORITHM IMPLEMENTED FOR FUNTION',37) GO TO ONTINUE ALL TNOU('N0 ALGORITHM IMPLEWENTED FOR PEAK AVERAGING',42) 93 ONTINUE ALL TNOU('0MPL IN HEADER NEGATIVE',24) 99 ONTINUE ALL EXIT END

42 SUBROUTINE RLNTP (BUF,LEN) INTEGER LEN INTEGER BUF (40) INTEGER BUFFER (40) ALL MSTR$A('REALL ',7,BUFFER,80) ALL MSUB$A (BUF,80,1,LEN,BUFFER,80,8,LEN+7) ALL MSUB$A (' FROM HP5420A TAPE LOATION', *28,1,28,BUFFER,80,LEN+8,LEN+35) LEN = LEN +34 ALL QUEST(BUFFER,LEN) ALL MSUB$A('RA1,2,1,2,BUFFER,8O,LEN+1,LEN+2) ALL AOW (BUFFER,LEN+2) RETURN END SUBROUTINE NVR4 (R4,BU FFER,LEN) REAL R4 INTEGER LEN, BUFFER(40) INTEGER*2 EXP,L,IBUF(5) INTEGER*4 IIUIAMT LOGIAL NVB$A EXP = 1 + IFIX(ALOGlO(ABS(R4))) INANT = INTL(IABS(INT(1000.0*(R4/lO~O**EXP)+.5))~ IF (.NOT.NVB$A(4,If+4AMT,IBUF,3)) GO TO 90 L= 1 IF (R4.LT.O.O) ALL MSTR$A('-',I,BUFFER,80) IF (R4.LT.O.O) L=2 ALL MSUB$A('. ',I,I,BUFFER,80,L,L) ALL MSUB$A(IBUF,10,1,3,BUFFER,8O,L+l,L+3) ALL IUISUB$A (' E',I,I,I,BUFFER,80,L+4,L+4) IF (.NOT.NVB$A(I,INTL(EXP),IBUF,3)) GO TO 90 ALL MSUB$A(IBUF,10,1,3,BUFFER,80,L+5,L+7) LEN = L+7 G(3 TO ALL TNOU('NVB$A ERROR IN NVR4',21) 99 ONTINUE RETURN END SUBROUTINE XYEXP(XY,WINVAL,P!AXVAL) INTEGER XY REAL RI MVAL,PlAXVAL LOGIAL SUB$A INTEGER*2 LEM1,L1,L2,L3,LEN2,IBUF(40),BUFFER(40) ALL NVR4(NINVAL,IBUF,LENI) ALL RsTR$Ä ipbuf,8u,buffer88ui ALL MSUB$A(',',1,1,1,BUFFER,8O,LENl+l,LENl+l~ ALL NVR4 (IViAXVAL, IBU F,LEN2) ALL MSUB$A~IBUF,80,1,LEN2,BUFFER,8O,LENl+2,LENl+LEN2+1~ IF (XY.EQ.l.OR.XY.EQ.'X') ALL MSUBSA *( ' X I,2,1,2,BU FFER,80,LENI +LEN2+2,LENI +LEN2+3) IF (XY.EQ.2.OR.XY.EQ.'Y1I ALL MSUB$A *( ' Y ',2,1,2,BU FFER,80,LENI +LEN2+2,LENI +LEN2+3) L2=1 L3=LEN1 +LEN2+3 DO 10 L1=1,L3 IF (SUB$A (BUFFER,80,LI,LI, ',I,I,I 1) GO TO 20 ALL MSUB$A (BUFFER,80,LI,LI,IBUF,80,L2,L2) L2=L ONTINUE 10 ONTINUE

43 ALL AOM(IBUF,L2-1) ALL WAIT(500) ALL AOM('EX',2) RETURN END