利用Madah内建的快速傅氏变换函数可以方便地在波数域计算重磁异常导数。介绍了基于Madab语言的波数域求导的算法,给出了程序源代码,讨论了一些有助于提高计算精度的编程技巧。通过模型试验和数据分析,发现在计算垂向导数时波数域求导算法的精度比傅氏级数的精度有明显改善;而水平导数的计算,2种方法的精度相当。在某区钾盐勘探中,用该方法处理高精度重力剖面数据,取得了较好效果。
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 1 D 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 calcula- tion, 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.