matlab高斯正反算程序6,[转载][转]大地主题解算——高斯平均引数法正反算Matlab源码...

再次强烈鄙视whu出版社！！！！！！

%-----------------------------------------------GaussGeodeticDirect.m-----------------------------------------

%高斯平均引数法大地主题正算

�, L2, A1, B1, L1, A2 单位 rad

%S 单位 m

function [B2, L2, A2] = GaussGeodeticDirect(B1, L1, A1, S, eps,

RefEllipsoid)

a = RefEllipsoid(1);

b = RefEllipsoid(2);

c = RefEllipsoid(3);

alpha = RefEllipsoid(4);

e_2 = RefEllipsoid(5);

e_2_ = RefEllipsoid(6); %e'_2

e = e_2^0.5;

dB0 = S*cos(A1)/(a*( 1-e_2 )*( 1-e_2*sin(B1)*sin(B1) )^-1.5);

dL0 = S*sin(A1)*sec(B1)/(a*( 1-e_2*sin(B1)*sin(B1) )^-0.5);

dA0 = dL0 * sin(B1);

Bm = 0;

Lm = 0;

Am = 0;

dB1 = 0;

dL1 = 0;

dA1 = 0;

while(1)

Bm = B1