用MATLAB实现将GPS的BLH坐标转换为直角坐标

% 以下是MATLAB代码实现将GPS的BLH坐标转换为直角坐标
% 输入经度、纬度、海拔高度
% 输出直角坐标系下的x、y、z坐标

function [x, y, z] = BLH2XYZ(lon, lat, h)
% 地球椭球体参数
a = 6378137; % 长半轴
b = 6356752.3142; % 短半轴
f = (a - b) / a; % 扁率
e2 = (a^2 - b^2) / a^2; % 第一偏心率的平方
% 弧度制转换为度数制
lon = lon / 180 * pi;
lat = lat / 180 * pi;
% 计算卯酉圈曲率半径
N = a / sqrt(1 - e2 * sin(lat)^2);
% 计算直角坐标系下的x、y、z坐标
x = (N + h) * cos(lat) * cos(lon);
y = (N + h) * cos(lat) * sin(lon);
z = (N * (1 - e2) + h) * sin(lat);
end