均匀球体重力异常正演模拟Matlab代码

均匀球体重力异常正演模拟Matlab代码

function gg=qtmx（）
%      均匀球体重力异常     %
%质心坐标（x0y0z0）
%xyz为采样点
%球心埋深D（m），球的半径R（m），剩余密度p（kg/m^3）
%剩余质量M=4/3*pi*R^3*p
%重力异常g=G*M*D/（x^2+y^2+D^2）^（3/2）  （m/s^2）
%G=6.67e-11（m^3/kg.s^2） 万有引力常数
%pi  圆周率
%数据保存在‘均匀球体重力异常.txt’

clear
clc

%球体参数%
r=100;        %球的半径
pd=1*10^3;    %剩余密度
x0=0;         %球心坐标（x0y0z0）
y0=0;
D=200;        %球心埋深z0

%采样区间%
x=（-2000:50:2000）;
y=（-2000:50:2000）;
z=0;

%常数%
G=6.67e-11;

%计算异常%
[x1y1]=meshgrid（xy）;   %生成网线节点矩阵
gg=G*（（4/3）*pi*r^3*pd）*（D-z）./（（（x1-x0）.^2+（y1-y0）.^2+（D-z）^2）.^（3/2））*10^5;%单位mGal

%成图%
figure（1）%图1
mesh（x1y1gg）%三维
xlabel（‘‘）
ylabel（‘‘）
title（‘均匀球体重力异常‘）
figure（2）
contourf（x1y1gg）%二维
title（‘均匀球体重力异常‘）

%数据生成文本%
%t=[x1（:）‘
%   y1（:）‘
%  gg（:）‘];
%fid=fopen（‘均匀球体重力异常.txt‘‘wt‘）; %wt以文本格式写入
%fprintf（fid‘%4.2f %4.2f %.2e\n‘t）;
%fclose（fid）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1006  2015-12-01 21:23  qtmx.m