MATLAB写一个ECEF转WGS-84的代码

MATLAB是一种强大的数值计算和数据可视化工具，用于编写地球坐标转换的代码时，可以利用其内置的数学函数和地理坐标处理库。以下是使用MATLAB编写从地球惯性坐标系(ECEF)到世界大地坐标系统(WGS-84)的简单代码示例，这个过程通常涉及到椭球参数和转换矩阵：

% WGS-84常数
a = 6378137; % 卫星椭球半径（平均）
f = 1 / 298.257223563; % 反照率

% ECEF坐标转换为WGS-84经纬度
function [lat, lon, h] = ecefToWgs84(xyz)
    % 地理坐标转换参数
    e2 = (1 - f) / (1 + f);
    e = sqrt(1 - e2);
    
    N = a / sqrt(1 - e^2 * sind(lat)^2); % 纬度正高
    h = xyz(3); % 高度
    x = xyz(1) / N;
    y = xyz(2) / N;
    z = xyz(3) / N;
    
    lat_rad = atan2(sqrt(1 - e^2 * z^2), z);
    lambda = atan2(y, x);
    
    % 转换为度
    lat = deg(lat_rad);
    lon = deg(lambda);
end

% 示例：假设xyz是一个3x1的ECEF坐标向量
xyz_example = [10000; 20000; 30000]; % 东经10000米，北纬20000米，海拔30000米
[lat_example, lon_example, h_example] = ecefToWgs84(xyz_example);