toevallige afwijkingen bij het meten van … · toevallige afwijkingen bij het meten van...
TRANSCRIPT
Toevallige afwijkingen bij het meten van eigenschappenvan random signalenvan Heck, J.G.A.M.
Gepubliceerd: 01/01/1982
Document VersionUitgevers 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 differencesbetween the submitted version and the official published version of record. People interested in the research are advised to contact theauthor 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
General rightsCopyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright ownersand 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 policyIf you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediatelyand investigate your claim.
Download date: 29. Aug. 2018
i
4 n
b I- -
A T
r 8. '
P
Y
P % o a
b
W
91. --"
.
k BI
c c
.
x
.. u
.-
e-%
9 a
U
3
Q
P
P
0.011 I I I I I I l l I I I I I 1 1 1 1
20 50 100 200 500 1,000 2,000 5,000 10.000 Number of Averages, nd
1
I
0.20
y 0.10 L-
w g E 0.05
m [r
0.03
0.,02
i c
Em-
u,
t i
6
A B
1
IS A N A L Y Z E R CONNECTED AND HP-IB SWITCH ON ADRESSABLE?
n J i J :
IS B A S I C LOADED?
R
OK?
B A S I C LOADED
IS B A S I C PROGRAM LOADED?
IB
B
RECALL ADDITIONAL COHERENCE ;FUNCTION FROM' WP5420A TAPE LOCATION
90,95 OR 99% CONF INTERVAL?
9
E SOLUTION?
NO ALGORITHM IMPLEMENTED FOR PEAK AVERAGING
0
lll
I i
A
, : .
P8-J
C
C C C
CONFIDENCE INTERVALS FOR REASUREMENTS DIGITAL SIGNAL ANALY
LOGICAL YESNO,CNVA$A REAL MAXEPS,NSIG,MAXI,MIMI,DATIJ1,DATIJ2 INTEGER*4 CONFI4 INTEGER 1,12,SIZE,FUNCOD,AVGS,PEAK,LEN,D1(9),
*BUFFER(40),Q,COMPL DIMENSION D2 (5) ,HEAD1 (16),HEAD2(16> ,DAT(512) ,EPSR(512),
CALL INIT *COH (51 2)
89995
99996
99994
C C C
99993
********** RECALL FUNCTION FROR ANALYZER
CALL ACOM( '1TC',3) CALL RCLNTP('AUTO,TRANS,COHER,CROSS',22) CALL SRQ(Q) IF (Q.NE.69) GO TO 90 CALL ACOH('RS',2) CALL WA IT ( 2000 1 CALL S501SA(HEADI,DAT,Q) IF (Q.NE.96) GO TO 90 CALL RHEAD(HEADl,SIZE,Dl ( l ) ,D l(2) ,D1(3) ,D1(4) ,
*CORPL,D1(6) ,DI (7) ,PEAK,D1(8) ,AVGS,D2(1), *D2 (2) ,O2 (3) ,D2 (4) , D l (9) ,D2 (5) 1
9 10
C C C
20
FUNCOD=DI (8) IF (FUMCOD.LT.7.0R.FUNCOD.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) FUNCOD FORP!AT('FUNCOD=',I6) WRITE(1.99994) COMPL FO RWAT ( ' CORPL= ' , 16) IF (FUMCOD.NE.9.AND.FUNCOD.NE.10) GO TO 10
********** RECALL ADDITIONAL COHERENCE FUNCTION
CALL ACOM('2TC',3) CALL RCLMTP('ADDITI0MAL COHER',l6) CALL SRQ(Q) IF (Q.ME.69) GO TO 90 CALL ACOM('RS',2) CALL WAIT(2000) CALL S501SA(HEAD2,COH,Q) IF (Q.NE.96) GO TO 90 DO 9 I=1,256 IF (COH(I).LT.O.O) WRITE (1,99993) COH(I),I FORMAT ('NEGATIVE COHERENCE ',F10.4,' ON DATAPOINT ',I41 IF (COH ( I ) . LT .O .O) COH ( I) =ABS (COH (I) 1 CONTINUE CONTINUE
********** DETERMINE NUMER OF SIGMAS
CALL MSTR$A('90,95 OR 99% CONF INTERVAL',26,
LEN = 26 CALL QUEST(BUFFER,LEN)
*BUFFER,80)
- 4 2 -
30
99998 C C C
C C C
35
40
C C C
45
50
C C C
I F (CNVA$A (1 ,BU F F E R , L E N , C O N F I 4 ) 1 GO T O 30 C A L L TNOU('ANSWER 90,95 OR 9 9 ' , 1 8 ) GO TO 20 CONTINUE I F ( I N T S ( C O N F I 4 ) .EQ.90) frfSIG=l.645 I F ( I N T S ( C O N F 1 4 ) . E Q D 9 5 ) NSIG=1.960 I F ( I N T S ( C O N F 1 4 ) . E Q D 9 9 ) NSIG=2.576 WRITE ( 1 , 9 9 9 9 8 ) MSIG FORRAT ( I NSIG=' , F1 O. 2)
********** B E N D A T g P I E R S O L ALGORITHNS FOR RANDOM ERRORS
IF (FUNCOD.NE.7.AND.FUNCOD.NE.8) GO TO 40
********** F I R S T Of? SECOND AUTO-POWER SPECTRUM
DO 35 1=1,256 E P S R ( I ) = I .O/SQRT(FLOAT(AVGS)) CONTINUE GO TO 70 CO WTI MU E I F (FUNCOD.NE.9) GO TO 50
********** CROSS POUER SPECTRUM
DO 45 I=1,256 I F ( C O H ( 1 ) .EQmO.O) C O H ( I ) = O ~ O O O O I EPSR (2*1-1)= 1 .O/ (SQRT ( F L O A T (AVGS)*COH (I) E P S R ( Z * I ) = EPSR(2*1-1) CONTINUE GO T O 70 CONTINUE I F (FUNCOD.NE.10) GO T O 60
********** TRANSFER FUNCTION
DO 55 1=1,256 IF (COH(1) .EQ.O.O) COH(I)=0.00001 IF ( C O H ( I ) . G T . I . O ) C O H ( I ) = l . O EPSR (2* 1-1 1 = SQRT ( 1-COH ( I) 1 / (SQRT (COH ( I) *SQRT ( F L O A T
* (2*AVGS))) E P S R ( Z * I ) = E P S R ( Z * I - I )
55 CONTINUE GO TO 90
60 CONTINUE IF (FUNCOD.NE.11) GO T O 91
C C ********** COHERENCE FUNCTION C
DO 65 1=1,256
E P S R ( I) =SQRT (2 .Cl) * (1 .O-DAT ( I) 1 / IF (DAT(I) .EQ.O.O) DAT(I)=O.OOQOI
* (SQRT ( DAT ( I) 1 *SQRT ( F L O A T (AVGS) 1 1 65 CONTINUE 70 CONTINUE
C C ********** ALGORITHM CONTROL C
DI (1)=0 DO 75 I=1,12
IF(EPSR(I).GT.O.Z) DI (l)=DI (1)+1 D2 (1 1 =FLOAT (DI (1 1) /FLOAT( 12>*100.0 IF (DI (l).NE.O) WRITE (1,99997) DI (1),D2(1)
75 CONTINUE
99997 FORMAT('ALGOR1THM ONLY VALID FOR EPSR < 0.2' ,/ , *14,'DATAPOINTS WITH EPSR > 0.2 ( ' , F5 .1 , '%>')
C C C
78
79
71
72
76
77
73
74
80 C C C
MAXI=DAT ( I I* (I .O+NSIG*EPSR ( I 1) MINI=MAX I I F (YESNO('1 TRACE SOLUTION',Ió)) GO TO 80
********** 2 TRACE SOLUTION ON CHANNEL A AND B
D O 72 1=1,12,2 DAT(I)=DAT( I>*(l .O+NSIG*EPSR(I) 1 DAT(I+I)=DAT(I+I)*(I.O+NSIG*EPSR(I+1)) IF (COPIPL) 93,78,79 DAT1 J I=DAT (I) DAT I J 2= DAT ( I+1) GO TO 71 DATIJI =SQRT(DAT (1)**2+DAT( I+1)**2> DATI J 2=DATI J 1 CONTINUE I F (DATIJI.LT.HINI)MINI=DATIJ1 IF (DATIJ2.LT.MINI)MINI=DATIJ2 IF (DATIJI.GT.MAXI)MAXI=DATIJI IF (DATIJ2.GT.MAXI)MAXI=DATIJ2 CONTINUE CALL ACOM('ITC',3) CALL UAIT(1000) CALL 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) ) D A T ( I + ~ ) = D A T ( I + I ) * ( ~ . O - N S I G * E P S R ( I + I ~ ) / ( ~ . ~ + ~ S I G * E P S R ~ I + ~ ~ ~
IF (COMPL) 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 CONTINUE I F (DATIJI.LT.MINI)MINI=DATIJ1 I F (DATIJ2.LT.NIMI)MINI=DATIJ2 I F (DATI J 1 .GT. MAX I) FlAX I=DATI J 1 IF (DAT1 52 .GT. MAX I)HAXI=DATI 52 CONTI NU E CALL WAIT(1000) CALL ACOM('ZTC',3) CALL WAIT(1000) CALL S501 R A (HEAD1 ,DAT,Q) IF (Q.NE.112) GO TO 90 GO TO 85 CONTINUE
********i* 1 TRACE SOLUTION ON CHANNEL B
DO 82 I=1,12,4 DAT ( I) =DAT ( 1) * (1 .O+MSIG*EPSR ( I) 1 I F (COIVIPL) 93,83,84
- @ 4 -
83
84
86
8%
89
81 87 82
85 C C C
99999
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 86 DAT(I+1)=DAT(I+1)*(1.O+NSIG*EPSR(I+I~) DAT(I+2)=DAT(I+2)*(1.O-NSIG*EPSR(I+2)) DAT(I+3)=DAT(I+3)*(1.O-NSIG*EPSR(I+3)) DO 87 J=û,3,2 IF (COMPL) 93,88,89 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 I F (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 I F (DATIJ 1 .LT. NINI) MINI=DATIJ 1 CONTINUE CONTINUE CONTINUE CALL ACOM('2TC'J) CALL WAIT( 1000) CALL S501RA(HEADI,DAT,Q) IF (Q.NE.112) GO TO 90 CONTI HU E
********** Y EXPAND ON BOTH TRACES
WRITE(l,99999> RIN1,MAXI FORHAT( '~INIMUR VALUE=',E10.3,/,
* MAX IHU PI VALUE= , E l O . 3 1 CALL ACOM('3FM1,3) CALL ACOM('OTC',3) CALL WAIT(500) CALL ACOH('RE',2) CALL WAIT(5OO) IF (COMPL.EQ.l .O> CALL ACOM('MG',2) CALL WA IT( 500) CALL XYEXP('Y',MINI,WAXI) GO TO 99
C c ********** ERROR HESSÄGES C
90 CONTINUE CALL TNOU('WR0NG STATUS SERVICE REQUEST',28) GO TO 99
91 CONTINUE CALL TNOU('N0 ALGORITHM IMPLEMENTED FOR FUNCTION',37) GO TO 99
92 CONTINUE CALL TNOU('N0 ALGORITHM IMPLEWENTED FOR PEAK AVERAGING',42)
93 CONTINUE
99 CONTINUE CALL TNOU('C0MPL IN HEADER NEGATIVE',24)
CALL EXIT END
SUBROUTINE RCLNTP (BUF,LEN) INTEGER LEN INTEGER BUF (40) INTEGER BUFFER (40) C A L L MSTR$A('RECALL ',7,BUFFER,80) C A L L MSUB$A (BUF,80,1 ,LEN,BUFFER,80,8,LEN+7) C A L L MSUB$A (' FROM HP5420A TAPE LOCATION',
LEN = LEN +34 C A L L QUEST(BUFFER,LEN) C A L L MSUB$A('RA1,2,1,2,BUFFER,8O,LEN+1,LEN+2) C A L L ACOW (BUFFER,LEN+2) RETURN END SUBROUTINE CNVR4 (R4,BU FFER,LEN) REAL R4 INTEGER LEN, BUFFER(40) INTEGER*2 EXP,L,IBUF(5) INTEGER*4 IIUIAMT L O G I C A L CNVB$A EXP = 1 + IFIX(ALOGlO(ABS(R4))) INANT = INTL(IABS(INT(1000.0*(R4/lO~O**EXP)+.5))~ I F (.NOT.CNVB$A(4,If+4AMT,IBUF,3)) GO TO 90 L= 1 I F (R4.LT.O.O) C A L L MSTR$A('-',I,BUFFER,80) I F (R4.LT.O.O) L=2 C A L L MSUB$A('. ' , I ,I ,I ,BUFFER,80,L,L) C A L L MSUB$A(IBUF,10,1,3,BUFFER,8O,L+l,L+3) C A L L IUISUB$A ( ' E' ,I ,I ,I ,BUFFER,80,L+4,L+4) I F (.NOT.CNVB$A(I,INTL(EXP),IBUF,3)) GO TO 90 C A L L MSUB$A(IBUF,10,1,3,BUFFER,80,L+5,L+7) LEN = L+7 G(3 TO 99
90 C A L L TNOU('CNVB$A ERROR I N CNVR4',21) 99 CONTINUE
RETURN END SUBROUTINE XYEXP(XY,WINVAL,P!AXVAL) INTEGER X Y REAL R I MVAL,PlAXVAL L O G I C A L CSUB$A INTEGER*2 LEM1,L1,L2,L3,LEN2,IBUF(40),BUFFER(40) C A L L CNVR4(NINVAL,IBUF,LENI)
C A L L MSUB$A(',',1,1,1,BUFFER,8O,LENl+l,LENl+l~ C A L L CNVR4 (IViAXVAL, IBU F,LEN2) C A L L MSUB$A~IBUF,80,1,LEN2,BUFFER,8O,LENl+2,LENl+LEN2+1~ I F (XY.EQ.l.OR.XY.EQ.'X') C A L L MSUBSA
I F (XY.EQ.2.OR.XY.EQ.'Y1I C A L L MSUB$A
L2=1 L3=LEN1 +LEN2+3 DO 10 L1=1,L3 I F (CSUB$A (BUFFER,80,LI ,LI, ' , I ,I ,I 1 ) GO TO 20 C A L L MSUB$A (BUFFER,80,LI ,LI ,IBUF,80,L2,L2) L2=L2+1
*28,1,28,BUFFER,80,LEN+8,LEN+35)
C A L L R s T R $ Ä iPBuF,8U,BuFFER88UI
*( ' X C I ,2,1,2,BU FFER,80,LENI +LEN2+2,LENI +LEN2+3)
*( ' Y C ' ,2,1,2,BU FFER,80,LENI +LEN2+2,LENI +LEN2+3)
20 CONTINUE 10 CONTINUE
C A L L ACOM(IBUF,L2-1) C A L L WAIT(500) C A L L ACOM('EX' ,2) R E T U R N E N D