projects / pdp8.git / blame
| Commit | Line | Data |
|---|---|---|
| 81e70d48 PH |
1 | C -+-+-+-+-+ \ e H A E P R P . F T \ e -+-+-+-+-+\r |
| 2 | C\r | |
| 3 | C PREP FOR THE FFT\r | |
| 4 | C\r | |
| 5 | SUBROUTINE PRPFFT\r | |
| 6 | INCLUDE HAEBUF.FI\r | |
| 7 | INCLUDE HAEGSA.FI\r | |
| 8 | INCLUDE HAEPTI.FI\r | |
| 9 | INCLUDE HAETTY.FI\r | |
| 10 | INCLUDE HABRK.FI\r | |
| 11 | C\r | |
| 12 | INTEGER LOG2\r | |
| 13 | EXTERNAL LOG2\r | |
| 14 | \fC\r | |
| 15 | INTEGER I,J,K,DA(3),BAND,LIMIT @ DO LOOP COUNTERS & LIMITS\r | |
| 16 | REAL TEMP,CPUTIM\r | |
| 17 | C\r | |
| 18 | C\r | |
| 19 | 10 CONTINUE\r | |
| 20 | C\r | |
| 21 | C HERE WE TEST FOR CORRECT FREQUENCY LIMITS OF THE EEG BAND (1-5)\r | |
| 22 | DO 12 BAND=1,5\r | |
| 23 | DO 12 LIMIT=1,2\r | |
| 24 | IF (FREQU(BAND,LIMIT).LE.SAMRAT/2) GOTO 12\r | |
| 25 | WRITE (TTO,101) BAND,SAMRAT\r | |
| 26 | FREQU(BAND,LIMIT)=COND(LIMIT.EQ.1,SAMRAT/2.-.1,SAMRAT/2.)\r | |
| 27 | 12 CONTINUE\r | |
| 28 | C\r | |
| 29 | IF (OPTION.NE.12 .OR. ENDS.LT.SPAN) CALL ASKHIM(4) @ CO* DOESN'T ASK SCILLY QUESTIONS IF THERE ARE CORRECT PARAMETERS\r | |
| 30 | C\r | |
| 31 | C TEST FOR VALID INPUT PARAMETERS\r | |
| 32 | C\r | |
| 33 | IF BREAK(11) RETURN @ USER GETS RID OF PROGRAM\r | |
| 34 | INCR=SAMRAT*SPAN @ NUMBER OF POINTS FOR EACH FFT --> INCR\r | |
| 35 | IF (INCR.LE.2048) GOTO 22 @ MORE THAN 2048 POINTS TO SCAN?\r | |
| 36 | WRITE (TTO,6) SPAN @ TOO LESS MEMORY\r | |
| 37 | GOTO 10\r | |
| 38 | 22 CONTINUE\r | |
| 39 | IF (SPAN.LE.ENDS) GOTO 23 @ LESS THAN ONE PASS?\r | |
| 40 | WRITE (TTO,8) SPAN,ENDS\r | |
| 41 | GOTO 10\r | |
| 42 | 23 CONTINUE\r | |
| 43 | IF (ENDS.LE.SAMCNT/SAMRAT) GOTO 24 @ IF THE END DOES NOT EXCEED THE SAMPLE COUNT THEN WE GO 24 ELSE\r | |
| 44 | I=SAMCNT/SAMRAT @ WE COMPUTE THE SAMPLING TIME IN SECONDS\r | |
| 45 | WRITE (TTO,2) ENDS,I @ IN ORDER TO PRINT AN DIAGNOSTIC MSG\r | |
| 46 | ENDS=SAMCNT/SAMRAT @ AND THE WE BUMP THE SAMPLING TIME INTO ENDS\r | |
| 47 | 24 CONTINUE\r | |
| 48 | EXPON=LOG2(INCR) @ TEST FOR POWER OF 2\r | |
| 49 | IF (EXPON.GT.0) GOTO 20 @ IF INCR IS A POWER OF TWO THEN GOTO 20\r | |
| 50 | WRITE (TTO,3) SAMRAT,SPAN @ NO POWER OF TWO\r | |
| 51 | GOTO 10 @ READ NEXT VALUE\r | |
| 52 | 20 TEMP=AMOD(FLOAT(ENDS-BEGIN),FLOAT(SPAN)) @ COMPUTE SECONDS TO TRANSFORM\r | |
| 53 | IF (TEMP.LT.0.1) GOTO 30 @ ARE THERE SOME SECONDS LEFT?\r | |
| 54 | TEMP=AMOD(FLOAT(SAMCNT)/FLOAT(SAMRAT),FLOAT(SPAN))\r | |
| 55 | WRITE (TTO,4) TEMP @ TELL THE USER THERE IS SOMETHING REMAINING\r | |
| 56 | 30 ENDS=BEGIN+(ENDS-BEGIN)/SPAN*SPAN @ COMPUTE THE END OF THE SPECTRUM TO SCAN\r | |
| 57 | C\r | |
| 58 | I=NCOND(OVRLAP,2*(ENDS-BEGIN)/SPAN,(ENDS-BEGIN)/SPAN) @ IF THE INTERMEDIATE RESULTS FIT ONTO UNIT 5 THEN\r | |
| 59 | I=I/4 @ FOUR DATA WORDS PER RECORD\r | |
| 60 | I=I+STSCAN @ COMPUTE NUMBER OF RECORDS NESESCARY\r | |
| 61 | IF (MAXBL5.GE.I) GOTO 40 @ RETURN ELSE PRINT AN ERROR MSG\r | |
| 62 | WRITE (TTO,9) I\r | |
| 63 | OPTION=0 @ PREVENT LOOPING IF CO*\r | |
| 64 | GOTO 10\r | |
| 65 | C\r | |
| 66 | C\r | |
| 67 | 40 FREQU(1,2)=AMIN1(SAMRAT/2.,FREQU(1,2)) @ HIGHEST FREQUENCY OF THE INPUT SIGNAL --> UPPER LIMIT OF THE FIRST BAND\r | |
| 68 | IF BREAK(11) RETURN\r | |
| 69 | CALL ADATE(DA) @ ASK FRTS FOR THE CURRENT DATE AND INSERT IT INTO THE HEADER\r | |
| 70 | WRITE (3,1) SETNUM,CHANEL,LABEL,DA\r | |
| 71 | WRITE (3,100) ((FREQU(I,J),J=1,2),I=1,5)\r | |
| 72 | RETURN\r | |
| 73 | C\r | |
| 74 | C\r | |
| 75 | 1 FORMAT ('1SATZ:',I3,' KANAL',I2,' :',10A6,2(A2,1H/),A2)\r | |
| 76 | 100 FORMAT (1H0,15X,5(F6.1,' -',F6.1,5X))\r | |
| 77 | 101 FORMAT (' FUER EEG BAND',I2,' REICHT DIE MESSRATE VON',I4,\r | |
| 78 | * ' HZ NICHT AUS. BANDGRENZE VERRINGERT.')\r | |
| 79 | 2 FORMAT (' RECHNUNG SOLL MIT',I5,' SEC UEBER DAS DATENENDE',\r | |
| 80 | * I5,' SEC HINAUSGEHEN.')\r | |
| 81 | 3 FORMAT (' DAS PRODUKT AUS MESSRATE',I3,' [HZ] UND SPANNE',\r | |
| 82 | * I3,' [SEK] IST KEINE ZWEIERPOTENZ.')\r | |
| 83 | 4 FORMAT (' AM ENDE DES SPEKTRUMS BLEIBEN',F4.1,' SEKUNDEN',\r | |
| 84 | * ' UNBERUECKSICHTIGT.')\r | |
| 85 | 6 FORMAT (' MEHR ALS 2048 PUNKTE FUER EINE FFT:',I7)\r | |
| 86 | 8 FORMAT (' SPANNE MIT',I3,' SEK. GROESSER ALS DIE MESSZEIT',I3\r | |
| 87 | * ,' SEKUNDEN.')\r | |
| 88 | 9 FORMAT (' FILE 5 FUER DIE ZWISCHENWERTE ZU KLEIN:',I5,\r | |
| 89 | * ' SPANNE VERGROESSERN.')\r | |
| 90 | END\r | |
| 91 | \1a\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0 |