说明:利用Matlab内建的快速傅氏变换函数可以方便地在波数域计算重磁异常导数。介绍了基于Matlab语言的波数域求导的算法,给出了程序源代码,讨论了一些有助于提高计算精度的编程技巧。通过模型试验和数据分析, 发现在计算垂向导数时波数域求导算法的精度比傅氏级数的精度有明显改善;而水平导数的计算,2种方法的...
说明:地球物理学重力勘探中,用最小二乘法对重力异常进行正反演的程序
说明:重力异常反演,《MATLAB程序设计及在地球物理中的应用》一书中的代码,手动敲的。
说明:一维二度重力异常向上及向下延拓都需要在已知剖面上取值,而且延拓高度应为取值点距的整数倍,因此称为等间距延拓。下面应用Matlab程序实现铅垂台阶向上、向下延拓10m的重力异常与原始异常的对比。 计算铅垂台阶的延拓异常,首先要建立原始异常的函数。铅垂台阶重力异常的表达式为: ∆g=Gσ[πH-...
说明:利用离散余弦变换求取重力或磁异常的水平与垂直导数,简单实用。