File: TDKT.FC of Tape: Various/ETH/fc2
(.FC,.FO Focal format converted to 'WRITE' listing)
C-OS/8 FOCAL, 1972 01.01 C DINGLE FIT: TDKT.FC 01.02 C FIT VON TD UND MU 01.03 C MIT K-TERM (KT) 01.05 E 01.10 T "FREQ:";S I=0 01.20 S I=I+1;A NU(I);I (-NU(I))1.2;S NF=I-1;G 1.3 01.30 A !?HI,DK? 01.40 T !" TEMP NU NUEX AMP TAU DAMP DTAU KT" 01.50 S J=J+1;A !,X(J);I (-X(J))1.6;S N=J-1;G 1.65 01.60 F I=1,NF;T !,%3,:7, NU(I);S L=NF*J+I;D 10 01.62 G 1.5 01.65 F I=1,NF;D 15 01.70 F K=1,NF;D 12;T !,"AMPFIT:";D 1.75; 01.72 Q 01.75 D 5;D 7;D 9;T "TAUFIT:";D 6;D 9;T "MU*PER",MU*512*DK/NE(K),! 02.01 C SUMS 02.10 S A=0;S B=0;S C=0;S D=0;S E=0 02.20 S DE=DY(1)^2;F I=1,N;S P(I)=DE/DY(I)^2 02.40 F I=1,N;S A=A+P(I)*X(I);S B=B+P(I)*Y(I);S C=C+P(I)*X(I)*X(I) 02.50 F I=1,N;S D=D+P(I)*X(I)*Y(I);S E=E+P(I)*Y(I)*Y(I) 02.55 S NP=0;F I=1,N;S NP=NP+P(I) 02.60 S G=NP*C-A*A 02.70 S AA=(B*C-A*D)/G;S BB=(NP*D-A*B)/G 02.80 S MU=BB/KK;S TD=AA/BB 03.01 C FEHLER UND PRINT 03.10 F I=1,N;S V(I)=(Y(I)-X(I)*BB-AA) 03.20 S Q=0;S QA=0;F I=1,N;S Q=Q+V(I)*V(I)*P(I);S QA=QA+(V(I)/Y(I))^2 03.30 S M=FSQT(Q/(N-1));S MA=M*FSQT(C/G);S MB=M*FSQT(NP/G) 03.40 S D1=MA/AA;S D2=MB/BB;S DT=TD*FSQT(D1*D1+D2*D2) 04.10 S H=NP*E-B*B;S RR=BB*BB*G/H;S RR=FSQT(RR);T %5.04 !,"RR =", RR 05.01 C AMPFIT 05.05 F I=1,N;S M=NF*I+K;D 5.1 05.07 G 5.2 05.10 S Y(I)=-FLOG(FCOM(M)/X(I));S DY(I)=FCOM(128+M)/FCOM(M) 05.20 S KK=1.47E5/HI;D 2;D 3 06.01 C TAUFIT 06.10 F I=1,N;S M=NF*I+K;S Y(I)=FCOM(64+M);S DY(I)=FCOM(192+M) 06.20 S KK=DK/.838;D 2;D 3 07.01 C TD AUS AMPFIT 07.05 S KK=DK/.838;S TM=0;F I=1,N;D 7.1 07.07 G 7.2 07.10 S M=NF*I+K;S TD(I)=FCOM(64+M)/(KK*MU)-X(I);S TM=TM+TD(I) 07.20 S TD=TM/N;S TQ=0 07.30 F I=1,N;S TQ=TQ+(TD-TD(I))^2 07.40 S DT=FSQT(TQ/(N-1)) 09.01 C PRINT 09.10 T " TD="%4.03,TD," (+/-",DT,")" 09.20 T " MU="%4.03,MU," (+/-",MB/KK,")" 09.30 T " RMSD=",%3.01,100*FSQT(QA/(N)),"%" 09.40 D 4;T ! 10.05 A :13,NE;S Z=FCOM(32+L,NE) 10.10 A :20,AM;S Z=FCOM(L,AM) 10.20 A :27,TA;S Z=FCOM(64+L,TA) 10.22 S Z1=FEXP(-TA*PI/(256*DK*1E-7*HI)) 10.24 S Z=FCOM(96+L,AM*FSQT(HI*1E-4)/(X(J)*Z1)) 10.30 A :33,DU;S Z=FCOM(128+L,DU) 10.40 A :39,DU;S Z=FCOM(192+L,DU) 10.50 T :44,%5.03,FCOM(96+L) 10.60 S Z=FX(0,500) 11.05 A ?TM,TA,AM,DK,HI,? 11.10 S AF=AM*(1-FEXP(-TA*PI))/(TA*PI) 11.20 A !?TA? 11.30 S Z1=FEXP(-TA*PI/(256*DK*1E-7*HI)) 11.35 S AM=AF*TA*PI/(1-FEXP(-TA*PI)) 11.40 S K1=AM/(TM*Z1) 11.50 T %6.05,K1;G 11.2 12.10 T !"FREQUENCY",%5.02,NE(K) 12.20 T " PERIOD",512*DK/NE(K) 12.30 T " K.TERM.",KT(K),"(+/-",KV(K),")",! 13.10 S TA=T1 13.20 S X=1-FEXP(-TA*PI);S Y=256*DK*1E-7*HI 13.30 S T1=-FLOG(AF*FSQT(HI*1E-4)*PI*TA/X(J)/X/KT(I))*Y/PI 13.40 I (1E-5-FABS(T1-TA))13.1;S TA=T1 15.05 S NE(I)=0 15.10 S KT(I)=0;S KV(I)=0;F J=1,N;S L=NF*J+I;D 15.95;D 15.9;D 15.98 15.15 S NE(I)=NE(I)/N 15.20 S KT(I)=KT(I)/N;S KV(I)=FSQT((KV(I)-N*KT(I)^2)/(N-1)) 15.30 F J=1,N;S L=NF*J+I;D 15.8;D 13;D 15.4;D 15.5 15.35 R 15.40 S AM=AF*TA*PI/(1-FEXP(-TA*PI)) 15.50 S Z=FCOM(64+L,TA);S Z=FCOM(L,AM) 15.80 S AM=FCOM(L);S T1=FCOM(64+L);S AF=AM*(1-FEXP(-T1*PI))/(T1*PI) 15.90 S KT(I)=KT(I)+FCOM(96+L) 15.95 S KV(I)=KV(I)+FCOM(96+L)^2 15.98 S NE(I)=NE(I)+FCOM(32+L)