blh坐标转换为xyz matlab代码

以下是将BLH（纬度、经度、海拔）坐标转换为XYZ坐标的Matlab代码：

function [x,y,z] = blh2xyz(latitude, longitude, altitude)
% BLH to XYZ conversion
% Reference: https://en.wikipedia.org/wiki/Geographic_coordinate_conversion#From_geodetic_to_ECEF_coordinates

a = 6378137; % semi-major axis of the Earth
f = 1/298.257223563; % flattening factor of the Earth

% Calculate parameters
b = a*(1-f);
e = sqrt((a^2-b^2)/a^2);
N = a/sqrt(1-e^2*sind(latitude)^2);
cos_lat = cosd(latitude);
sin_lat = sind(latitude);
cos_lon = cosd(longitude);
sin_lon = sind(longitude);

% Convert BLH to XYZ coordinates
x = (N+altitude)*cos_lat*cos_lon;
y = (N+altitude)*cos_lat*sin_lon;
z = ((b^2/a^2)*N+altitude)*sin_lat;
end

其中，`latitude`是纬度，以度为单位，`longitude`是经度，以度为单位，`altitude`是海拔，以米为单位。函数返回的`x`，`y`，`z`是对应的XYZ坐标，以米为单位。