jwdis.m

此函数通过两点的大地经纬度（Ja,Wa）和（Jb,Wb）计算两点间距离 Dis

function Dis=jwdis(Ja,Wa,Jb,Wb)
% Dis=jwdis(Ja,Wa,Jb,Wb)
% 此函数通过两点的大地经纬度（Ja,Wa）和（Jb,Wb）计算两点间距离 Dis;
% 输入的经纬度值为浮点数，单位为 度，即：xxx.xxx 度;
% 输出的距离为浮点数,单位为 km;

ja=Ja*60;
wa=Wa*60;
jb=Jb*60;
wb=Wb*60;

Dis=sqrt(((ja-jb)*1.855*cos(((Wa-Wb)/2)*pi/180))^2+((wa-wb)*1.849)^2);