高斯正反算批量计算

高斯正反算 通过matlab与Excel正反算进行批量计算。正算时Excel第一列输入纬度，第二列输入经差，sheet2中放椭球参数 反算是Excelsheet1中放坐标，sheet2放椭球参数，sheet3放中央经线

function s=BF（Xaee）
%UNtitleD 输入X坐标，迭代计算大地纬度
%   此处显示详细说明
m（1）=a*（1-ee）;
m（2）=3*ee/2*（m（1））;
m（4）=5*ee/4*（m（2））;
m（6）=7*ee/6*（m（4））;
m（8）=9*ee/8*（m（6））;
A（1）=m（1）+m（2）/2+m（4）*3/8+m（6）*5/16+m（8）*35/128;
A（2）=m（2）/2+m（4）/2+15*m（6）/32+7*m（8）/16;
A（4）=m（4）/8+3*m（6）/16+7*m（8）/32;
A（6）=m（6）/32+m（8）/16;
A（8）=m（8）/128;
BF1=X/A（1）;
BF2=X/A（1）;
while（（BF2-BF1）<（10^（-10）））
    BF1=BF2;
    Fa=-A（2）/2*sin（2*BF1）+A（4）/4*sin（4*BF1）-A（6）/6*sin（6*BF1）+A（8）/8*sin（8*BF1）;
    BF2=（X-Fa）/A（1）;
end
s=BF2;
end

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

     文件        542  2018-11-13 16:02  高斯坐标转换\高斯反算\BF.m

     文件        290  2018-11-09 11:48  高斯坐标转换\高斯反算\D2R.m

     文件        158  2018-11-09 11:50  高斯坐标转换\高斯反算\R2D.m

     文件       1253  2018-11-13 15:52  高斯坐标转换\高斯反算\ZBJS.m

     文件       8461  2018-11-13 16:01  高斯坐标转换\高斯反算\反算的纬度和经差的结果.xlsx

     文件      10917  2018-11-13 15:14  高斯坐标转换\高斯反算\坐标数据.xlsx

     文件        290  2018-11-09 11:48  高斯坐标转换\高斯正算\D2R.m

     文件        886  2018-11-11 16:08  高斯坐标转换\高斯正算\JWDJS.m

     文件        158  2018-11-09 11:50  高斯坐标转换\高斯正算\R2D.m

     文件        435  2018-11-09 12:40  高斯坐标转换\高斯正算\X.m

     文件        435  2018-11-09 12:40  高斯坐标转换\高斯正算\ZWXHC.m

     文件      10320  2018-11-13 15:16  高斯坐标转换\高斯正算\经纬度.xlsx

     目录          0  2018-11-13 16:03  高斯坐标转换\高斯反算

     目录          0  2018-11-13 15:25  高斯坐标转换\高斯正算

     目录          0  2018-11-13 14:43  高斯坐标转换

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

                34145                    15