user_ecef.m

GPS/GALILEO/BEIDOU兼容接收机自主完好性检测程序

function [xu,yu,zu]=user_ecef(long,lat,h)

%*********************************
% ** %经纬度坐标转化为地心地固坐标 
% **********************************


Re=6378137;
f=1/298.257;
trans=2*pi/360;
RN=Re./sqrt(((1-f).^2).*(sin(trans*lat))^2+(cos(trans*lat)).^2);
xu=(RN+h)*cos(trans*lat).*cos(trans*long);
yu=(RN+h)*cos(trans*lat).*sin(trans*long);
zu=(RN*((1-f)^2)+h)*sin(trans*lat);