Matlab模拟均匀球体剖面重力异常正演 上传版本.zip

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(x,gg) xlabel('') ylabel('') title('均匀球体剖面重力异常') %数据生成文本% y=[x;gg]; fid=fopen('均匀球体剖面重力异常.txt','wt'); %wt以文本格式写入 fprintf(fid,'%4.2f %.8e\n',y); fclose(fid);