

c 



FAST FUUHI.EX TK1NSFUHM SUBROUTINE 

 SUbHUUTlNE FFT(FH»FI«K f ICO) 

 UlnENSION FR(l).FHl) 



5 





IMKO.EQ.O) GO TU 10 



oo e i»i,n 





b 



FI(I)b-FKI) 





10 



CONTINUE 



10 



l 



MK«0 



NNSN» 1 



DO 2 M»1,NN 



L'N 



l»L/J 



15 





1F(MH*L.GT,NN) GO TO 1 

 MWiMUOtMHtLUL 

 IF(MH.LE.H) GO TO 2 

 TR»FK(M* j ) 

 F«(Mtl)»FR(MB+i) 



20 





FH(MK+| )«TH 



TI»F1(M+I) 



FI(M*1)«FI(MR*H 



FI(MH+1).TI 





2 



CONTINUE 



25 





La I 





3 



IF(L.GE.N) GU TO 7 



1STEP»2»L 



tL»L 



UU II Msl.L 



30 





»«3,lU159|6535*FL0/tT(l-M)/E|. 



*H»CU8(A) 



MlPSIN(A) 



DO tt I«M,N,I8TEP 

 J'lH 



35 





IF(ICO.EO.l) GO TO tl 



T««WH»FR(J)»WI«FI(J) 

 TI«*K*FI(J)*WI*FR(J) 

 GO TU 12 





11 



TR«wH»FRCJ)fwI»«I(J) 



110 





TI»WH*FI(J)-WM»H(J) 





12 



FHCJUFBCD-TR 

 FI(J)«FI(l)-TI 

 FR(»»FR(I)*TR 





« 



FI(I)«FI(I)*TI 



45 





LsISTEP 

 GO TU 3 





7 



CONTINUE 



ANSN 



IF(ICO.EQ.l) GO TO 6 



50 





00 5 I«l t N 

 FR(I)=FR(I)/AN 





5 



FICI)«-FI(I)/AN 





6 



HfcTURN 

 tNU 



Figure 3 



;. Lis 



ting of FFT subroutine. 



b. Data Input to Program . 

 f(j) values at 

 At = 1 second 

 intervals 

 (32 values) 



FR = 



(32 values) 



6.000 

 0.590 



-2.600 

 0.562 

 2.000 



-1.390 



-1.400 

 4.238 

 6.000 

 0.590 



-2.600 



5.080 



3.750 



2.184 



0.829 



-1.900 



-2.506 



2.215 



-1.451 



-0.465 



1.445 



2.034 



2.229 



1.391 



0.513 



-0.475 



2.054 



-2.322 



-2.109 



0.257 



1.188 



2.755 



5.438 



6.189 



6.386 



5.080 



3.750 



2.184 



0.829 



-1.900 



-2.506 



2.215 



-1.451 



-0.465 



25 



