地球椭球参数计算

%地理坐标 Latitude—纬度 Longitude—经度 CenterLongitude 中央子午线经度（单位：度）%
Latitude=30.4777183851287;
Longitude=114.339694365886;
CenterLongitude =114;

B=Latitude*pi/180;
L=Longitude*pi/180;
L0=CenterLongitude*pi/180;

 

%椭球参数 a-地球长半轴　b-短半轴%
a=6378137;
b=6356752.31424518;
%椭球参数 f-扁率 e-椭圆的第一偏心率 e1-椭圆的第二偏心率 %
f=(a-b)/a;
e=sqrt(a*a-b*b)/a;
e1=sqrt(a*a-b*b)/b;
%椭球参数 W-第一基本纬度函数 V-第二基本纬度函数%
W=sqrt(1-e*e*sin(B)*sin(B));
V=sqrt(1+e1*e1*cos(B)*cos(B));

%椭球参数 M-子午圈曲率半径 N-卯酉圈曲率半径%
M=a*(1-e*e)/(W*W*W);
N=M*V*V;

%平均曲率半径R 取过地面某点的所有方向的平均值来作为这个球体的半径近似计算  %
R=sqrt(M*N);