MATLAB中根据两点的经纬度求距离

重点是用matlab计算距离得到的单位：
 

%% 初始化参数
clear all;
close all;
clc;
%% 读取两点的坐标
coordinate_a = [39.85929,116.47070];
coordinate_b = [39.04403,121.63660;
    26.60968,106.73550;
    45.74223,126.62370;
    36.62035,114.49450;
    30.37639,120.16980;
    31.890180,117.35390];
number_b = size(coordinate_b,1)
%% 计算两点间的距离
distance_ab = [];
for i = 1:number_b
    % diatance_a的第一列的单位：度
    distance_ab(i,1) = distance(coordinate_a(1),coordinate_a(2),coordinate_b(i,1),coordinate_b(i,2))
    % diatance_ab_2的第二列单位：千米
    distance_ab(i,2) = distance(coordinate_a(1),coordinate_a(2),coordinate_b(i,1),coordinate_b(i,2))/180*pi*6371;
end

计算出的距离结果为：

其中distance_ab中的第一列的距离单位为度，第二列的距离单位为千米