matlab 经纬度 东北天 直角坐标系转换代码

要在matlab中进行经纬度到东北天直角坐标系的转换，你可以使用以下代码：

function [x, y, z] = convertToENU(latitude, longitude, altitude, refLat, refLong, refAlt)
    % 地球半径
    R = 6371000;
    
    % 弧度转换
    lat1 = deg2rad(latitude);
    lon1 = deg2rad(longitude);
    lat0 = deg2rad(refLat);
    lon0 = deg2rad(refLong);
    
    % 计算东北天直角坐标系中的位置
    cosLat = cos(lat1);
    sinLat = sin(lat1);
    cosLon = cos(lon1);
    sinLon = sin(lon1);
    
    cosLat0 = cos(lat0);
    sinLat0 = sin(lat0);
    cosLon0 = cos(lon0);
    sinLon0 = sin(lon0);
    
    e = (R * cosLat) / sqrt(1 - 0.00669437999014 * (sinLat * sinLat));
    
    N = R / sqrt(1 - 0.00669437999014 * (sinLat * sinLat));
    
    x = (e + altitude) * cosLat * sin(lon1 - lon0);
    y = (e + altitude) * (cosLat0 * sinLat - sinLat0 * cosLat * cos(lon1 - lon0));
    z = (N + altitude) * sin(lat1 - lat0);
    
end

这个函数接受经度（latitude）、纬度（longitude）和高度（altitude）作为输入参数，以及参考点的经度（refLong）、纬度（refLat）和高度（refAlt）。然后根据经纬度和地球的半径，计算出对应的东北天直角坐标系中的位置，并返回x、y和z的值。