%数值的单位必须是公里
a=6378137.0;
b=6356752.3142;
f=1/298.25722353;
e2=f*(2-f);
wgs84 = wgs84Ellipsoid(‘kilometer’);
%大地坐标系转成地心地固坐标系
[x,y,z] = geodetic2ecef(wgs84,28,112,122);
%地心地固坐标系转成大地坐标系
[lat,lon,alt]=ecef2geodetic(wgs84,x,y,z);
%站心坐标系转地心地固坐标系
[X, Y, Z] = enu2ecef(10000, 20000, 1000, 28, 110, 100, wgs84);
%地心地固坐标系转站心坐标系
[xEast, yNorth, zUp] = ecef2enu(X, Y, Z, 28, 110, 100, wgs84);
%站心坐标系转大地坐标系
[LAT, LON, H] = enu2geodetic(xEast, yNorth, zUp, LAT0, LON0, H0,wgs84);
%大地坐标系转站心坐标系
[xEast, yNorth, zUp] = geodetic2enu(LAT, LON, H, LAT0, LON0, H0,wgs84);