从经纬度坐标到欧式距离的转化【matlab实现】

如果给定两点坐标（在平面直角坐标系的情况下），设A（x1,y1），B（x2,y2），很容易可以计算两点距离为（x2-x1）^2+（y2-y1）^2的（1/2）次方,

但给定两点真实坐标，如经纬度，就按照如下方法计算：首先不妨假设点A,点B三维空间向量为：

基于此，我用matlab实现编程：

编程如下：

[location,txt]=xlsread('省会城市经纬度.xlsx','Sheet1','A1:C35');
[m,n]=size(location);
T=zeros(m,m);
for i=[1:1:m]
    for j=[i+1:1:m]
      T(i,j)=distance(location(i,1),location(i,2),location(j,1),location(j,2));
    end
end
for i=[1:1:m]
    for j=[1:1:m]
        T(j,i)=T(i,j);
    end
end
xlswrite('省会城市欧氏距离表.csv',T,'B2:AI35');
xlswrite('省会城市欧氏距离表.csv',txt(:,1),'A1:A35');
xlswrite('省会城市欧氏距离表.csv',txt(:,1)','A1:AI1');