File: FRELEG.FC of Tape: Various/ETH/fc1
(.FC,.FO Focal format converted to 'WRITE' listing)
C-FRELEG.FC APR-77 01.10 I (LX),1.2;G LX 01.20 L N FOCCOM.ZE 01.30 L P FOUPIC 01.35 E 01.40 S LX=2.1 02.10 A !"LEGENDRE FIT: LOWEST ANGLE ",LA," HIGHEST ANGLE ",HA 02.20 A "ANGLE INCREMENT",IA 02.25 A !" LOWEST FREQU ",Z," HIGHEST FREQU ",Y;S OF=(Z+Y)/2 02.27 S SO=0 02.30 S NP=(HA-LA)/IA;C # OF INTERVALS 02.35 S XP=NP*3;I (FITR(NP/2)-NP/2)2.99 02.40 S IP=FITR(256/NP);S HP=IP*NP 02.45 F A=LA,IA,HA;T !"ANGLE= ",%3,A," FREQU";A Z;S Y=FCOM((A-LA)*IP/IA,Z-OF);D 13.5 02.50 S P=0;S P(1)=0;S P(2)=0;S P(3)=0;S P(4)=0 02.55 F A=LA,IA,HA;D 3 02.60 S P=P/XP;S P(1)=P(1)*3/XP;S P(2)=P(2)*5/XP;S P(3)=P(3)*7/XP;S P(4)=P(4)*9/XP 02.65 T !"COEFFS ",%6.04,P,P(1),P(2),P(3),P(4) 02.70 S SL=0;S U=0;F I=0,HP;D 4 02.73 I (-SO)2.75;A !"ONLY SECOND ORDER=1 ",SO;I (SO),2.75,2.7 02.75 F I=0,HP;D 23 02.80 F A=LA,IA,HA;T !"ANGLE= ",%3,A," DLNA/DPHI= ",%6.03,FCOM((A-LA)*IP/IA)*180/PI 02.90 S LX=2.1;Q 02.99 T !"BITTE EIN WINKEL MEHR ODER WENIGER!!!!";G 2.1 03.20 S Z=FCOM((A-LA)*IP/IA) 03.30 S J=(2*A-LA-HA)/(HA-LA) 03.31 C NEXT TWO LINES SIMPSON RULE 03.32 I (-FABS((A-LA)*(A-HA)))3.33,3.35 03.33 S Z=Z*2;I ((A-LA)/2/IA-FITR((A-LA)/2/IA)),3.35;S Z=Z*2 03.35 S P=P+Z 03.40 S P(1)=P(1)+Z*J 03.45 S P(2)=P(2)+Z*(3*J^2-1)/2 03.50 S P(3)=P(3)+Z*(5*J^3-3*J)/2 03.55 S P(4)=P(4)+Z*(35*J^4-30*J^2+3)/8 04.10 S A=I*IA/IP+LA;S J=(2*A-LA-HA)/(HA-LA) 04.17 S X=0;S Y=0;I (SO),4.2,4.4 04.20 S X=(35*J^4-30*J^2+3)/8 04.30 S Y=(5*J^3-3*J)/2 04.40 S Z=P+P(1)*J+P(2)*(3*J^2-1)/2+P(3)*Y+P(4)*X 04.45 I (-SL*(U-Z))4.47,4.47;T " MINIMUM AT PHI= ",%6.04,I*IA/IP+LA 04.47 S SL=U-Z;S U=Z 04.50 I (FABS(FITR(I/IP)-I/IP)),4.6;S Y=FCOM(I,Z);D 13.5;R 04.60 T !"ANGLE= ",%3,A," DELTA NU= ",%6.04,FCOM(I)-Z 05.01 C FOUR 05.04 S DL=DK*.512;I (DL),5.07,5.05 05.05 T "FREQ. IN TESLA / FMAX=",%6.02,500/DK 05.06 A !,"BEG"SI,"END"SF;S SI=SI*DL;S SF=SF*DL;G 5.1 05.07 A ?SI,SF? 05.10 S X=(SF-SI)/16;I (X-FITR(X)-.5)5.14;S X=X+.5 05.14 S SS=FITR(X)/16;S S2=SI+SS*128;S SF=SI+SS*256 05.16 I (-DL)5.17;T %3,"'"SS*32"'";G 5.2 05.17 T %4.01,"'"SS*32/DL,"'" 05.20 T !;D 13.9;D 12 05.30 S Z=FOUR(SI,SS);S ZM=FX(1,760) 09.01 C TAU + NUMAX 09.10 S X=256+(NU-SI)/SS 09.20 I (FCOM(X+1)-FCOM(X))9.25,9.25;S X=X+1;G 9.2 09.25 I (FCOM(X-1)-FCOM(X))9.3,9.3;S X=X-1;G 9.25 09.30 S XM=X;S AH=.5*FCOM(XM);S NU=(XM-256)*SS+SI 09.40 S X=X-1;I (FCOM(X)-AH)9.45,9.45,9.4 09.45 S X1=X;S A1=FCOM(X);S A2=FCOM(X+1);S X=XM+1 09.50 I (FCOM(X)-AH)9.55;S X=X+1;G 9.5 09.55 S X3=X-1;S A3=FCOM(X3);S A4=FCOM(X) 09.60 S X1=X1+(AH-A1)/(A2-A1);S X3=X3+(A3-AH)/(A3-A4) 09.70 S TA=(X3-X1)*SS/2 09.71 S TA=TA*N/256 09.72 I (TA-.95)9.74;S TA=TA^4-.75;S TA=FSQT(FSQT(TA));G 9.8 09.73 I (-TA)9.74,9.74;S TA=.1;G 9.8 09.74 I (TA-.885)9.9;S TA=FSQT((TA-.885)/.28) 09.80 S TA=TA*256/N;R 09.90 T "!TAU!";S TA=0.01 10.01 C PEAKEVAL 10.04 A "FREQ"NF;S NU=NF*DL;I (-FABS(NU))10.05;A ?I?;S NU=NF(I)*DL 10.05 I (NU-256)10.06;T !!,"FREQ. TOO BIG!!",!;Q 10.06 S SI=NU-8;S X=1;D 5.14;D 5.2;D 5.3 10.10 D 9;D 20;T %5.03, "NF",NU/DL," PH",PH," AF",AF," TF",TA*2*PI/DL,! 10.20 S NY=NU;F I=152,155;D I/10 10.30 A !"STORE AS NF(I)? ",?I?;I (-FABS(I))10.4;R 10.40 S NF(I)=NU/DL;S PH(I)=PH;S AF(I)=AF;S TF(I)=TA*2*PI/DL 12.01 C FOURSCALE 12.04 S S1=SI 12.05 I (DL),12.1;S S1=SI/DL;S S2=S2/DL;S SF=SF/DL 12.06 I (SF-1000)12.1;S Z=10^(FITR(.4343*FLOG(SF)-2));S S1=S1/Z;S S2=S2/Z;S SF=SF/Z;T "FREQS*",%4,Z,! 12.10 S X=194;S A=S1/100;S Y=FITR(A+.00001) 12.12 I (-Y)12.14;S Z=FDIS(X+2,-96);G 12.16 12.14 S L=12.2;G 12.52 12.16 D 12.8;I (-Y)12.18;S Z=FDIS(X+1,-96);G 12.8 12.18 D 12.9;G 12.8 12.20 S X=201;S A=S2/100;S Y=FITR(A);S L=12.3;G 12.5 12.30 S X=211;S A=SF/1000;S Y=FITR(A);I (-Y)12.32;D 12.54;G 12.34 12.32 D 12.9;D 12.8;D 12.9;S L=12.4;G 12.56 12.34 D 12.8;S L=12.4;G 12.5 12.40 R 12.50 I (-Y)12.52,12.54 12.52 D 12.9;G 12.56 12.54 S Z=FDIS(X,-96);S X=X+1 12.56 D 12.8;D 12.9;D 12.8;D 12.9;G L 12.80 S A=(A-Y)*10+1E-4;S Y=FITR(A) 12.90 S Z=FDIS(X,Y-80);S X=X+1 12.92 G 12.2 13.01 C DIS SWITCH 13.40 S Z=FDIS(0,0) 13.50 S Z=FX(0,500) 13.60 S Z=FDIS(1,760) 13.70 S Z=FDIS(A,FABS(FDIS(A)));S Z=FDIS(B,-FABS(FDIS(B))) 13.80 S A=4;S B=5;D 13.7 13.90 S A=5;S B=4;D 13.7 15.01 C ALIASED FREQUENCIES 15.05 A ?DK?,! 15.10 A ?NU?;S NY=NU 15.20 T " FREQ=",%5.02 15.30 T NY/DL 15.40 S NY=512-NU;D 15.3 15.50 S NY=512+NU;D 15.3 15.90 T !;G 15.1 16.01 C PLOT 16.02 C WITH AXES AND TEXT: PA=0, ELSE PA=2 16.05 S PA=2;A ?PA? 16.10 S PF=4;I (FDIS(5))16.2;S PF=5 16.20 F I=2,19;S Z=FDIS(I,-FABS(FDIS(I))) 16.70 S X=PF+PA;D 16.9;S Z=FDIS(21,1) 16.80 D 16.9;S X=PF;D 16.9;S $=1024;Q 16.90 S Z=FDIS(X,-FDIS(X)) 18.01 C SQ.ROOT OF POWERSPECTRUM(FCOM(1)) 18.10 F I=256,511;S Z=FCOM(I,FSQT(FCOM(I)));S Z=FX(1,760) 19.10 A !,X;T FCOM(X);G 19.1 20.01 C FOUGNAU 20.10 S N1=NU*PI/256;S A=0;S B=0 20.20 F I=0,N-1;S A=A+FCOM(I)*FCOS(I*N1) 20.25 I (NU*(NU-256))20.3,20.4,20.3 20.30 F I=0,N-1;S B=B+FCOM(I)*FSIN(I*N1) 20.40 S AF=2*FSQT(A^2+B^2)/N 20.50 I (A)20.7,20.6;S PH=0;G 20.8 20.60 S PH=90*FSGN(B);R 20.70 S PH=180*FSGN(B) 20.80 S PH=PH+180*FATN(B/A)/PI 20.95 S AM=AF*TA*PI*N/(256*(1-FEXP(-TA*N*PI/256))) 21.01 C SUB PEAK 21.10 D 13.8 21.15 A ?NF?;I (-FABS(NF))21.2;A ?I? 21.16 S NF=NF(I);S PH=PH(I);S AF=AF(I);S TF=TF(I);G 21.22 21.20 A ?PH AF TF MP? 21.22 I (-FABS(MP))21.25;S MP=1;T " MP=1!!" 21.25 A !"ADD",Z;I (-FABS(Z-0YES))21.3;S AF=-AF 21.30 S NO=NF*2*PI*DK*1E-3;S PO=PH*PI/180;S TO=TF*DK*1E-3;S A1=HI*DK*1E-3 21.31 S A=AF*TO*N/(1-FEXP(-N*TO));S H6=MP*NO*TM*.0048/A1 21.32 F J=0,N-1;S H1=1+J*A1;D 21.4;D 13.5 21.33 R 21.35 C IF MODULATION CHANGE FSQT(H1) 21.40 S Z=FCOM(J,FCOM(J)-A*FCOS(J*NO-PO)*FEXP(-J*TO)/((1-FEXP(-H1*H6))*FSQT(H1))) 23.10 S A=I*IA/IP+LA;S J=(2*A-LA-HA)/(HA-LA) 23.17 S X=0;S Y=0;I (SO),23.2,23.4 23.20 S X=(140*J^3-90*J^2)/8 23.30 S Y=(15*J^2-3)/2 23.40 S Z=(P(1)+P(2)*3*J+P(3)*Y+P(4)*X)*2/(HA-LA) 23.50 S Z=FCOM(I,Z/(FCOM(I)+OF));D 13.5