数据标准化 将地球离散为25km*25km的单元格 子午线半径近似为40009km，赤道半径近似为40076km 因此，近似的，地球被划分成1603.04个（长）*1600.36个（宽）的2,565,441.0944个正方形 每个单元格对应经度0.2245733107096516618424992514223度（取小数点后六位），不妨设为a 对应纬度0.11247469319403134294783673673423度（取小数点后六位），设为b 建立平面直角坐标系并且令原点为（-90°，-90°） 之后执行如下操作： 得到用户每次登陆位置信息，可以进行统计得到用户家的位置
home=zeros(58227,3); count=zeros(1604+5,1604+5); max=0; j=-1; for i=1:4747281 if Data(i,1)==j && Data(i,2)>0 && Data(i,3)>0 if Data(i,2)<1604 && Data(i,3)<1604 count(Data(i,2),Data(i,3))=count(Data(i,2),Data(i,3))+1; if count(Data(i,2),Data(i,3))>max max=count(Data(i,2),Data(i,3)); max_x=Data(i,2); max_y=Data(i,3); end else count(1604,1604)=count(1604,1604)+1; if count(1604,1604)>max max=count(1604,1604); max_x=Data(1604,2); max_y=Data(1604,3); end end else count=zeros(1604,1604); if j~=-1 home(j+1,1)=j; %因为是从0开始的 home(j+1,2)=max_x; home(j+1,3)=max_y; end j=Data(i,1); max=0; if Data(i,2)>0 && Data(i,3)>0&& Data(i,2)<1604 && Data(i,3)<1604 count(Data(i,2),Data(i,3))=1; end if Data(i,2)>0 && Data(i,3)>0&&Data(i,2)<1604 && Data(i,3)<1604&& count(Data(i,2),Data(i,3))>max max=count(Data(i,2),Data(i,3)); max_x=Data(i,2); max_y=Data(i,3); end end end根据家的位置求distance
Distance=zeros(4747281,1); for i=1:4747280 a=(abs(Data(i,2)-home(Data(i,1)+1,2)))^2; b=(abs(Data(i,3)-home(Data(i,1)+1,3)))^2; Distance(i,1)=2500*sqrt(a+b); end使用ksdensity函数做出概率密度曲线
