function?EGM96jisuan（vlH）
%?使用這個程序選擇EGM96文件即可以計算給定經(jīng)緯度和高度的重力異常，大地水準面差距和垂線偏差。
%?注意，下面使用了許多已知量，這些量都可以用EGM96canshu這個小程序求得，所以為了提高計算效率，此處全部用了計算出來的已知量。
%?本程序有個很大缺陷就是計算速度很慢，大約需要十幾秒。我覺得主要原因是循環(huán)比較多，可以考慮以后簡化循環(huán)，同時提高讀取數(shù)據(jù)的速度。
%?請沈老師指正。——曲騰騰
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的數(shù)據(jù)‘）;
fit1=fopen（strcat（pathname1filename1）‘rt‘）;
if（fit1==-1）
????msgbox（‘Input?File?or?Path?is?not?correct‘‘Warning‘‘warn‘）;
????return;
end???
h=1;??%h只是為了控制數(shù)組的下標不要為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