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

为了计算WGS-84坐标下的四个岛礁之间的相对距离，我们需要先确定每个点的三维坐标（经度、纬度和海拔）。在MATLAB中，我们可以手动编写一个函数来执行球面三角法[^3]，因为不调用库函数意味着我们不能直接使用内置的经纬度距离计算工具。

下面是一个简单的示例，展示如何手动计算两个点之间的大圆距离：

function distance = sphericalDistance(lat1, lon1, lat2, lon2)
    % 转换成弧度
    lat1_rad = deg2rad(lat1);
    lon1_rad = deg2rad(lon1);
    lat2_rad = deg2rad(lat2);
    lon2_rad = deg2rad(lon2);

    % 经度差（弧度）
    deltaLon = lon2_rad - lon1_rad;

    % 利用Haversine公式计算球面距离
    a = sin(deltaLat/2)^2 + cos(lat1_rad).*cos(lat2_rad).*sin(deltaLon/2)^2;
    c = 2*atan2(sqrt(a), sqrt(1-a));
    R = 6371; % 地球平均半径（km）
    distance = R * c; % 单位为公里
end

% 假设有四个岛礁的坐标
latitudes = [your_latitude_1, your_latitude_2, your_latitude_3, your_latitude_4];
longitudes = [your_longitude_1, your_longitude_2, your_longitude_3, your_longitude_4];
altitudes = [your_altitude_1, your_altitude_2, your_altitude_3, your_altitude_4];

% 计算所有对的距离
distances = zeros(length(latitudes), length(latitudes));
for i = 1:length(latitudes)
    for j = i+1:length(latitudes)
        distances(i,j) = sphericalDistance(latitudes(i), longitudes(i), latitudes(j), longitudes(j));
        distances(j,i) = distances(i,j); % 对称性
    end
end

请注意，这段代码假设所有输入的纬度和经度都是以度为单位，且忽略了海拔对距离的影响。如果要考虑到海拔，可以稍微修改Haversine公式中的计算部分，但那通常会更复杂一些。