1 REM *** ANALYSIS OF A WAVE PROFILE *** 



2 REM *** BY FOURIER SERIES APPROXIMATION *** 



3 REM THIS PROGRAM MUST BE RUN FROM A POGITION A STATIO\ Plt^FORE 



4 REM A MAXIMUM AMPLITUDE. CALCULATE THE NUMBER OF STATIONS TO BE 



5 REM MEASURED (S) AND THE NUMBER OF POINTS TO BE USED IN THE 



6 REM PROGRAM DEPENDING UPON FREQUENCY (Dl). SI>D1+15 



7 REM T:PERI0D,D1:N0. OF PTS. USED, K U ORDER SOUGHT 



8 REM Q:MAX VOLTAGE TO TRIGGER, Sr NO. OF STATIONS EXPECTED 



9 REM S1:N0. OF DATA PTS.,T5:TEST NUMBER, C 1 =CONVERSION( FT. - V) 



10 DIM A[101,B[10],C[50, 10],F[50],E[3,50],G[50,5«],R(4],Z[3] 



11 READ T,D1,K1,Q,S,S1,T5,C1^ 



12 G0SU3 1000 



13 PRINT 



lA PRINT " *****=***** TEST NUNBER" , T5 , "**********" 



15 PRINT 



16 PRIMT " CORRECTION FACTOR r",Cl 



17 PRINT 



18 FOR NUl TO S 



19 CALL (1 ,D,F) 



20 IF D >3 THEN 22 



21 GOTO 19 



22 CALL (1,D,F) 



23 IF D <: Q THEN 25 



24 GOTO 22 



25 FOR N2: 1 TO SI 



26 CALL (1,D,F) 



27 LET G[Nl,N2]rD 



28 NEXT N2 



31 FOR X: 1 TO Dl+1 



55 LET F[X]: G[N1,X+15] 



35 NEXT X 



40 GOSUB 75 



45 NEXT Nl 



50 GOSUB 380 



60 GOSUB 735 



70 STOP 



75 REM SUBPROGRAM TO CALCULATE MAXIMUM AMPLITUDE USING FOURIER 



76 REM SERIES AT EACH STATION. 

 80 FOR X= 1 TO Dl+1 



85 LET F[X]: Gt N1,X+15 ] 



90 NEXT X 



91 LET E:D1+1 



92 LET G:K1 



100 LET F[l ]: (F[l >t-FtEl)/2 



102 LET F[E]=F[1 ] 



1 1 1 LET A4=1/2*F[ 1 ] 



1 12 LET H:T/D1 



1 15 FOR J:2 TO Dl 



1 16 LET A4:A4-t-F[ J] 



1 18 NEXT J 



120 LET A: 1/(2*3. 14I6)*H*(A4+l/2*F[Dl + n) 



35 



