经纬高转xyz坐标 matlab的代码

以下是将经纬高转换为xyz坐标的MATLAB代码：

function [x, y, z] = llh2xyz(lat, lon, h)
% Convert geodetic coordinates (latitude, longitude, height) to Cartesian
% coordinates (x, y, z).
% Input:
%   lat: geodetic latitude in degrees
%   lon: geodetic longitude in degrees
%   h: height above the ellipsoid in meters
% Output:
%   x: Cartesian x-coordinate in meters
%   y: Cartesian y-coordinate in meters
%   z: Cartesian z-coordinate in meters

% Constants for WGS84 ellipsoid
a = 6378137; % semimajor axis in meters
f = 1/298.257223563; % flattening

% Convert latitude and longitude to radians
lat = deg2rad(lat);
lon = deg2rad(lon);

% Compute the radius of curvature in the prime vertical
N = a ./ sqrt(1 - (2*f - f^2)*sin(lat).^2);

% Compute Cartesian coordinates
x = (N + h) .* cos(lat) .* cos(lon);
y = (N + h) .* cos(lat) .* sin(lon);
z = ((1 - f)^2 * N + h) .* sin(lat);
end

使用示例：

lat = 39.9; % latitude in degrees
lon = 116.4; % longitude in degrees
h = 50; % height above the ellipsoid in meters
[x, y, z] = llh2xyz(lat, lon, h);

此代码将返回北京市中心的WGS84坐标系中的xyz坐标。