MATLAB:计算数组坐标外的LatLon距离

function distance = calculateDistance(lat1, lon1, lat2, lon2)
    R = 6371; % 地球半径，单位为千米

    lat1 = deg2rad(lat1);
    lon1 = deg2rad(lon1);
    lat2 = deg2rad(lat2);
    lon2 = deg2rad(lon2);

    dlat = lat2 - lat1;
    dlon = lon2 - lon1;

    a = sin(dlat/2)^2 + cos(lat1) * cos(lat2) * sin(dlon/2)^2;
    c = 2 * atan2(sqrt(a), sqrt(1-a));

    distance = R * c; % 最终距离，单位为千米
end

lat1 = 40.7128; % 纽约市的纬度
lon1 = -74.0060; % 纽约市的经度
lat2 = 34.0522; % 洛杉矶的纬度
lon2 = -118.2437; % 洛杉矶的经度

distance = calculateDistance(lat1, lon1, lat2, lon2);
disp(distance); % 输出结果