matlab经纬度距离计算及数据分类合并程序

A=xlsread('F:\BaiduYunDownload\我的资源\大创\CBD\house.xlsx','a1:d1707'); %读入地铁站附近小区经纬度数据,第一列小区编号，第3列地铁站编号，第4列距离 A1=xlsread('F:\BaiduYunDownload\我的资源\大创\CBD\小区去0+7042.xls','a1:c7042'); %第一列房价，第二三列经纬度 for i=1:7042 for j=1:1562 if A(j,1)==i A2(j,:)=A1(i,:); end end end xlswrite('F:\BaiduYunDownload\我的资源\大创\CBD\house_price_loc.xlsx',A2); %*****************************A2为6号线沿线小区房价及经纬度******************************************* %***************************************************************************************************** B=[116.368122,39.921484;116.323066,39.989956;116.458594,39.910288]; %金融街，中关村，CBD for i= 1: 1562 for j=1:3 DIS(i,j)=dis(A2(i,2),A2(i,3),B(j,1),B(j,2)); %DIS为各小区距三大中心的距离 end [C(i,1),C(i,2)] = min(DIS(i,:)); %C为各小区距某一中心的距离。将最小距离放在1，将离什么近放在2 end xlswrite('F:\BaiduYunDownload\我的资源\大创\CBD\house_center_dis.xlsx',C) xlswrite('F:\BaiduYunDownload\我的资源\大创\CBD\house_3center_dis.xlsx',DIS)