已知四个岛礁的wgs-84坐标（包括海拔），求这四个岛礁之间的距离。用matlab实现

在MATLAB中计算两个WGS-84坐标点之间的距离通常使用Haversine公式，因为考虑到地球是个球体。以下是如何使用该公式来计算四个岛礁之间距离的一个简单示例：

首先，假设你有四个岛礁的经纬度（经度`lat1`, `lon1`, `lat2`, `lon2`）和对应的海拔`alt1`, `alt2`（在海平面高度）。注意，MATLAB中的纬度是从南极开始逆时针记为正，经度从本初子午线向西为正。

% 岛礁坐标 (以度为单位)
lat1 = [your_latitude_1; your_latitude_2; your_latitude_3; your_latitude_4];
lon1 = [your_longitude_1; your_longitude_2; your_longitude_3; your_longitude_4];
alt1 = [your_altitude_1; your_altitude_2; your_altitude_3; your_altitude_4];

% 将纬度和经度转换为弧度
lat_rad = deg2rad(lat1);
lon_rad = deg2rad(lon1);

% Haversine公式计算两点间的距离 (忽略海拔影响)
a = sin((lat_rad(2,:) - lat_rad(1,:))./2).^2 + cos(lat_rad(1,:)) * cos(lat_rad(2,:)) .* sin((lon_rad(2,:) - lon_rad(1,:))./2).^2;
c = 2 * atan2(sqrt(a), sqrt(1-a));
dist = earth_radius * c; % earth_radius是一个常数，地球平均半径（约为6371公里）

% 添加海拔差异的影响（如果需要）
if ~isempty(alt1)
    alt_diff = alt1(2:end) - alt1(1:end-1); % 仅考虑相邻岛屿的高度差
    dist = dist + hypot(dist, alt_diff); % 使用勾股定理
end

% 结果存储在一个矩阵中
dist_matrix = dist';

请注意，这段代码假定了地球是一个完美的球体（使用的是平均半径`earth_radius`），实际情况中可能需要考虑椭球形状的影响，特别是对于大尺度的距离。另外，海拔差的处理依赖于具体的应用场景，上述代码仅仅给出了一个基本的近似方法。