matlab地心坐标系和GPS坐标系的转换

**

卑微小兰的第一篇：matlab地心坐标系和GPS坐标系的转换，emmm

GPS坐标系和大地坐标系的正反变换

a=6378137;
f=0.00669437999014;%第一偏心率平方
B=40;
L=50;
H=60;
N=a/((1-f*(sind(B)))^0.5)%N为地心到该点的空间距离
X=(N+H)cosd(B)cosd(L);
Y=(N+H)cosd(B)sind(L);
Z=(N(1-f)+H)sind(B);%%我们得到了天线基于大地坐标系的几个坐标XYZ
syms B1;
syms H1;
syms L1;
ep3=atan(Y/X)180/pi-L1;
ep4=Z/sin((B1pi)/180)-N(1-f)-H1;
ep5=atan(((Z(N+H1))/(((X^2+Y2)^0.5)(N(1-f)+H1))))*180/pi-B1;
s=solve(ep3,ep4,ep5,H1,L1,B1);