wavelet.for

来自「Fortran编的波动方程正演模拟,一维的数值算法」· FOR 代码 · 共 19 行

      SUBROUTINE WAVELET(WT,N,FR,DT)
	DIMENSION WT(N)
	REAL WT
	REAL FR,DT
	PI=ATAN(1.0)*4.0
	N2=N/2
	DO I=1,N2
		A1=PI*FR*DT*(I)
		A2=A1*A1
		B1=EXP(-A2)
		WT(N2+1+I)=(1.0-2.0*A2)*B1
	END DO
	DO I=1,N2
		WT(I)=WT(N-I+1)
	END DO
	WT(N2+1)=1.0
	WRITE(*,*)WT
	RETURN
	END