曲面上重力异常的计算之matlab程序
葛粲
在三维重力异常正演问题中,最常见的就是在直角坐标系下将异常体划分为一个个均匀密度的立方体,计算各个立方体的在测点产生的重力异常,然后累加起来。这种方法是最直接和显而易见的方法,但是在编程计算时往往效率很低,主要是涉及到多层循环(5层循环),既要对空间中的异常体循环(3层循环),又要对观测平面循环(2层循环)。利用matlab强大的矩阵计算功能,我们可以一次性获得一个密度异常立方体在整个观测平面上产生的观测异常。这样只需要对密度异常体进行循环并累加就可以完成三维重力异常正演,大大提高了计算效率。
一个密度异常立方体在观测点P产生的重力异常公式推导如下:
matlab子程序:
function [ gs ] = SurfaceGravityByPrism( xx,yy,zz,x,