利用matlab实现，wgs-84坐标与ecef坐标转换

在MATLAB中，可以使用这些公式来实现WGS-84坐标系（大地坐标，经纬度和海拔）到地球中心直角坐标系（ECEF，即地心笛卡尔坐标系）的转换。以下是一个简单的示例：

% 定义WGS-84椭球参数
a = 6378137; % 半径赤道处（单位：米）
b = a * sqrt(1 - pow((a^2 / pow(EarthEquatorialRadius, 2)), 2)); % 半径极点处
f = (a - b) / a; % 扰动率

% 给定的经纬度和海拔
lat = deg2rad(your_latitude); % 将纬度从度转换为弧度
lon = deg2rad(your_longitude); % 同上，将经度转换为弧度
height = your_elevation; % 海拔（单位：米）

% 计算辅助变量
N = a / sqrt(1 - f^2 * sin(lat)^2);
enu_x = (N + height) .* cos(lat) .* cos(lon);
enu_y = (N + height) .* cos(lat) .* sin(lon);
enu_z = (pow(a, 2) / N + height) .* sin(lat);

% ECEF坐标
ecef_x = enu_x;
ecef_y = enu_y;
ecef_z = ecef_z;

% 显示结果
disp(['ECEF coordinates: ', num2str(ecef_x), ', ', num2str(ecef_y), ', ', num2str(ecef_z)]);

请确保替换`your_latitude`, `your_longitude`, 和 `your_elevation`为实际的经纬度和海拔值。