user_ecef.m

北京到广州铁路路线仿真

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);