matlab语言 lla2ecef,Matlab的 lla2ecef函数的C语言实现

下面是一个简单的C语言实现版本：

#include <math.h>

// 地球半径
#define EARTH_RADIUS 6378137.0

void lla2ecef(double lat, double lon, double alt, double* x, double* y, double* z) {
    double cos_lat = cos(lat);
    double sin_lat = sin(lat);
    double cos_lon = cos(lon);
    double sin_lon = sin(lon);
    double N = EARTH_RADIUS / sqrt(1.0 - 0.00669437999014 * sin_lat * sin_lat);

    *x = (N + alt) * cos_lat * cos_lon;
    *y = (N + alt) * cos_lat * sin_lon;
    *z = (N * (1.0 - 0.00669437999014) + alt) * sin_lat;
}

该函数的参数包括：

- `lat`：纬度，单位为弧度
- `lon`：经度，单位为弧度
- `alt`：高度，单位为米
- `x`：输出参数，ECEF坐标系中的x分量，单位为米
- `y`：输出参数，ECEF坐标系中的y分量，单位为米
- `z`：输出参数，ECEF坐标系中的z分量，单位为米

该函数计算给定的地心经纬度坐标系中的点的ECEF坐标系中的坐标。注意，这里的地心经纬度坐标系并不是标准的WGS84坐标系，而是经过一些简化处理的版本。