It is convenient to calculate derivative in wave number space using the inbuilt FFT function in Matlab. This paper describes a matlab algorithm for derivative calculation of 1D gravity and magnetic anomaly in wave number space, gives code of program and deals with some programming skills for precision improvement. Through model test and data analysis, it is found that the precision of the FFT method is much higher than that of the Fourier Series method in vertical derivative calculation. In horizontal derivative calculation, however, the precisions of the two methods are almost the same. During kainite exploration in a certain basin, satisfactory result was obtained in the processing of gravity profile data using this method.
肖锋, 孟令顺, 吴燕冈. 在波数域计算一维重磁异常导数的Matlab语言算法[J]. 物探与化探, 2008, 32(3): 316-320.
XIAO Feng, MENG Ling-shun, WU Yan-gang. A MATLAB ALGORITHM FOR DERIVATIVE CALCULATION OF 1D GRAVITY AND MAGNETIC ANOMALY IN WAVE NUMBER SPACE. Geophysical and Geochemical Exploration, 2008, 32(3): 316-320.