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