File: NOISY.FC of Tape: Various/ETH/fc1

(.FC,.FO Focal format converted to 'WRITE' listing) 

---

C-OS/8 FOCAL, 1972

01.01 C NOISY
01.10 A ?Q?,?NU?,?TA?
01.20 A ?N?;D 10
01.30 D 2;D 3
01.40 T #,:30,%4.03,"VM=",VM,:46,"VV=",VV
01.50 R

02.01 C V MESSEN
02.10 F I=0,N-1;D 2.2
02.20 S Z=FADC(Q);S V=Z/10^T2;S Z=FCOM(I,V)

03.01 C STATISTIK
03.10 S VA=0;F I=0,N-1;S VA=VA+FCOM(I)
03.20 S VM=VA/N
03.30 S VB=0;F I=0,N-1;S VB=VB+(FCOM(I)-VM)^2
03.40 S VV=FSQT(VB/(N-1))

04.01 C CORRELATION ZW. V(N) UND V(N-1) >>> DIS
04.10 S Z=FDIS(1,9);S Z=FDIS(10,1);S Z=FDIS(11,128);
04.20 S Z=FDIS(12,128);S Z=FDIS(13,6);F I=14,17;S Z=FDIS(I,10)
04.22 S Z=FDIS(18,10)
04.30 S A=50/VV;S I=0;D 4.5;D 5.1;D 5.3
04.40 F I=1,N;D 4.5;D 5
04.50 S X=A*(FCOM(I)-VM)

05.01 C POINTPLOT SR
05.10 S X=4*(X-256)
05.20 S Z=FDIS(18+2*I,X)
05.30 S Z=FDIS(19+2*I,X)

06.01 C REPETIEREN
06.10 A ?M?;F J=1,M;D 6.2
06.20 D 2;D 3;D 1.4;T !

07.01 C DIS-CORREL.-PROG.
07.10 F I=-510,-19;S Z=FDIS(-I,0)
07.20 D 2;D 3;D 4

08.01 C AUTOCORRELATION
08.10 F I=0,N-1;S Z=FCOM(I,FCOM(I)-VM)
08.20 A ?M?;S R=0;D 8.3;S C=C1;D 8.4;F R=1,M;D 8.3;D 8.4
08.25 R
08.30 S C1=0;F J=0,N-R-1;S C1=C1+FCOM(J)*FCOM(J+R);S Z=FCOM(256+R,C1)
08.40 T %3,100*C1/C,"/",%4.03

10.01 C GATE TIME
10.10 A ?T?;S Z=FADC(T);S T2=4-T

11.01 C VV AUF ANDERE ART
11.10 S V1=0;F I=0,N-1;S V1=V1+(FCOM(I))^2
11.20 S V2=N*VM^2;S V3=V1-V2;S VS=FSQT(V3/(N-1))

---