Abstract:The procedure that the Gaussian quadrature rule is introduced into Fourier transformation calculation is shown in rigorous mathematical deduction based on the Shift sampling DFT (SFT) theory.The conclusion is that a Fourier integral can be approached with high precision by the weighted sum of a series of SFT with different shift parameters.The weigh is a half of the Gaussian quadrature coefficient,and the sampling parameter is a half of the Gaussian quadrature node plus 0.5.This conclusion provides a rigorous theoretical basis for improving the precision of forward modeling of potential fields (gravity and magnetic) in wave-number domain.Since the sufficient conditions of both the SFT theory and the Gaussian quadrature theory are limited real functions,the Gaussian FFT algorithm based on the conclusion are applicable in forward and inverse Fourier transforms of any limited real functions.
柴玉璞, 万海珍. 移样离散傅里叶变换在重磁勘探中的应用[J]. 石油地球物理勘探, 2020, 55(6): 1358-1363.
CHAI Yupu, WAN Haizhen. Shift sampling DFT(SFT) theory and its application in gravity and magnetic prospecting. Oil Geophysical Prospecting, 2020, 55(6): 1358-1363.
Chai Y P and Hinze W J.Gravity inversion of an interface above which the density contrast varies exponentially with depth[J].Geophysics,1988,53(6):837-845.
[2]
Wu L Y and Tian G.High-precision Fourier forward modeling of potential fields[J].Geophysics,2014,79(5):G59-G68.
[3]
柴玉璞.Fourier变换数值计算的偏移抽样理论[J].中国科学(E辑),1996,26(5):68-74.CHAI Yupu.Shift sampling theory of fourier transform computation[J].Science in China(E),1997,40(1):21-27.