functionCds=GeoCdschange(X,Y,Z) %已知大地直角坐标系求解大地坐标系算法 采用WGS84坐标系 a=6378134; b=6356752.3141; e_1_2=0.00669437999013; e_2_2=0.00673949674227; %X=13194419.145061; %Y=13194419.145061; %Z=18629484.032597;