用 matlab 把WGS84 转ECEF坐标系

在 MATLAB 中，将 WGS84 经纬度坐标转换到地球椭球体（如 WGS84 ellipsoid）上的 ECEF （Earth-Centered, Earth-Fixed，地心、静止参考系）坐标通常涉及地理坐标系统的转换。这需要一些数学计算，包括球面坐标到平面直角坐标的转换以及椭球坐标到笛卡尔坐标的转换。

以下是大致步骤：

1. **导入经纬度数据**：
使用 `geodetic2enu` 函数，它接受 WGS84 地理坐标 (lat, lon) 和海拔高度作为输入。

[lon, lat, height] = geodetic2enu(lat_deg, lon_deg, 0); % lat_deg 和 lon_deg 是以度为单位的纬度和经度

2. **假设高度为0**，因为 ECEF 参考系是相对于地面的。如果考虑海拔高度，需要对结果进行平移。

3. **从球面坐标到平面直角坐标**：
在 WGS84 球模型下，先将纬度转换为弧度，然后使用球面三角学公式得到 X 和 Y 分量。

R = 6378137; % WGS84 半径（米）
lat_rad = deg2rad(lat);
x = R * cos(lat_rad) * cos(lon_rad);
y = R * cos(lat_rad) * sin(lon_rad);

4. **得到 Z 分量**（垂直距离）：
这是在 WGS84 海拔上直接加上 `height` 的值。

z = height;

5. **组合 XYZ 坐标**：
最终的 ECEF 坐标就是 `(x, y, z)`。

注意：这个过程假设了 WGS84 模型是一个完美的椭球体，实际上在实际应用中可能会引入更复杂的校正和精度问题。如果你需要更高精度的转换，可能需要使用专门的地理信息库，如 MATLAB 的 `geodetic toolbox` 或者开源库 `proj4j`。