191816. lajstromszámú szabadalom • Eljárás változó amplitudójú vagy fázisú szeizmikus kijelzési görbe előállítására szeizmikus jelenséget jellemző változóról felvett szabálytalan lefutású vagy szinuszos jellegű mérési görbék alapján, különösen geológiai kutatásoknál

r»o onoon ooooo

1

191.816

2

RST (I)»REAL(TRACE(I))

CST (I)»AIMAG(TRACE(I))

20 CONTINUE

DO 30 1*2,LH RST(I)*2.0»RST(I)

CST(I)* 2.OSCST(I)

30 CONTINUE

DO 40 I-LH1.LT RST(I)*0.0 CST(I)«0.0 40 CONTINUE DO 45 I-l.LT

TRACE(I)*CMPLX(RST(I),CST (I))

45 CONTINUE

CALL FFT(LT. TRACE. 1.0)

DO 47 1*1.LT

RST (I)«REAL(TRACE(I))

CST(I)* AIMAG(TRACE(I))

47 CONTINUE

XXXXXXXXXXXXXXXXXlXXXXXKXXXXXKXKXXXXXtXXKXXXXX * CALCULATE PHASE VALUES I

XXXXXXXXXXKXXSXXXKXXXXXXXXXKXXXXXXXXXKXXXXXtXX

DO 50 I-l.LT

IF(RST(I).EQ.O.AND.CST(I).EQ.O) CO TO 55 ENV(I)•AT AN2(CST(I) , RST(I))

55 IF(RST(I).EQ.O.AND.CST(I).EQ.O) ENV(I)«0.0 50 CONTINUE

DO 65 I-l.LT

WRITE(6. 60) RST(I). CST(I). ENV(I),PUT(I)

65 CONTINUE 60 FORMAT(10X.4F10.3)

RETURN

END

SUBROUTINE FFT(LX. CX.SIGNI)

LX-NUMBER OF SAMPLE POINTS-MUST BE A POWER OF 2

CX«INPUT AND OUTPUT DATA

S-i<m A/Ja

hWtSlt "f/ÍA BtCrVf'-* FOR rejiftQfib ntWtF&ftf SIGNI IS DESIGNATED IN MAIN PROGRAM COMPLEX CX (LX). CAF.G. CEXP. CW. CTEMP J«1

SC«1./FLOAT(LX)

DO 30 I«l,LX lF(I.GT.J) GO TO 10 CTEMP«CX(J)

CX(J)»CX(I)

CX(I)«CTEMP 10 M* LX/2

20 IF(J.LE.M) GO TO 30 J« J-M M-M/2

18