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

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

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

clear
clc

%球的参数%
d=200;      %球心埋深D
r=100;       %球的半径R
p=2*10^3;        %剩余密度p
G=6.67e-11;

%计算异常%
x=（-1000:10:1000）;   
m=（4/3）*pi*r^3*p;
n=x.^2+d^2;
gg=G*m*d./（n.^（3/2））*10^5;  %单位mGal

%成图%
plot（xgg）
xlabel（‘‘）
ylabel（‘‘）
title（‘均匀球体剖面重力异常‘）

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         711  2015-12-01 21:10  qtmxpm.m