EGM96计算重力异常的Matlab程序

EGM96计算重力异常的Matlab程序

function EGM96jisuan（vlH）
% 使用这个程序选择EGM96文件即可以计算给定经纬度和高度的重力异常，大地水准面差距和垂线偏差。
% 注意，下面使用了许多已知量，这些量都可以用EGM96canshu这个小程序求得，所以为了提高计算效率，此处全部用了计算出来的已知量。
% 本程序有个很大缺陷就是计算速度很慢，大约需要十几秒。我觉得主要原因是循环比较多，可以考虑以后简化循环，同时提高读取数据的速度。
% 请沈老师指正。——曲腾腾
j2=0.108262982131E-2;
j4=-0.237091120053E-5;
j6=0.608346498882E-8;
j8=-0.142681087920E-10;
j10=0.121439275882E-13;
R=6371000.7900;
gm=0.3986004418e+15;
a=6378137;
e2=0.00669437999013;
geqt=9.7803253359;
k1=0.00193185265246;
f=1/298.257223563;
[filename1pathname1]=uigetfile（‘*.txt‘‘读入EGM96的数据‘）;
fit1=fopen（strcat（pathname1filename1）‘rt‘）;
if（fit1==-1）
    msgbox（‘Input File or Path is not correct‘‘Warning‘‘warn‘）;
    return;
end   
h=1;  %h只是为了控制数组的下标不要为0。 
while  ~feof（fit1）
    line = fgetl（fit1）;
    if ~isstr（line） break end %读到

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    3005546  2008-12-25 01:37  EGM96计算程序\EGM96.txt

     文件       3401  2013-09-29 09:35  EGM96计算程序\EGM96jisuan.asv

     文件       3401  2008-12-25 23:16  EGM96计算程序\EGM96jisuan.m

     文件        101  2013-09-10 20:22  EGM96计算程序\新建文本文档.txt

     目录          0  2013-09-29 09:35  EGM96计算程序

----------- ---------  ---------- -----  ----

              3012449                    5