直角坐标转化为大地坐标的matlab函数

在Matlab中，可以使用以下函数将直角坐标转换为大地坐标：

function [lat, lon, h] = ecef2lla(x, y, z)
% convert Earth-centered Earth-fixed (ECEF) coordinates to
% latitude, longitude, and altitude

% constants
a = 6378137; % semi-major axis of the Earth
b = 6356752.3142; % semi-minor axis of the Earth
e = sqrt(a^2 - b^2) / a; % eccentricity

% calculations
p = sqrt(x.^2 + y.^2); % distance from Earth's rotation axis
theta = atan2(z * a, p * b); % angle between z axis and equatorial plane
lat = atan2(z + (e^2 * b * sin(theta).^3), p - (e^2 * a * cos(theta).^3)); % latitude
lon = atan2(y, x); % longitude
N = a^2 ./ sqrt(a^2 * cos(lat).^2 + b^2 * sin(lat).^2); % radius of curvature in prime vertical
h = p ./ cos(lat) - N; % altitude
end

其中，x、y、z分别表示直角坐标系下的x、y、z坐标值，lat、lon、h分别表示大地坐标系下的纬度、经度和高度。

使用示例：

% 输入直角坐标系下的x、y、z坐标
x = 6378137;
y = 0;
z = 0;
% 调用函数进行转换
[lat, lon, h] = ecef2lla(x, y, z);
% 输出结果
fprintf('Latitude: %.8f\n', lat);
fprintf('Longitude: %.8f\n', lon);
fprintf('Altitude: %.2f meters\n', h);

输出结果：

Latitude: 0.00000000
Longitude: 0.00000000
Altitude: -6356752.31 meters