确定罪犯藏身地点问题2.docVIP

4．确定罪犯藏身地点问题 已知罪犯系列犯罪的地点坐标如下： Place Dimension Longitude 1 53°497.73N 1°3230.83W 2 53°4828.74N 1°3158.44W 3 53°5028.40N 1°2932.95W 4 53°4838.71N 1°4549.71W 5 53°494.08N 1°3155.91W 6 53°2557.38N 2°153.14W 7 53°4759.94N 1°4619.78W 8 53°3915.67N 1°4646.64W 9 53°2740.55N 2°1326.73W 10 53°4241.97N 1°5226.05W 11 53°4721.45N 1°4547.33W 12 53°4830.40N 1°4016.68W 13 53°4920.17N 1°3440.65W （1）将角度制的坐标转换为弧度制坐标； （2）将系列犯罪的地点（13个）的球面坐标转换为平面坐标；并画出这些点的图形（要求带适当宽度的网格，行宽1.25*10^3，列宽h=2.5*10^3）。 （3）求这13个点之间距离的最大值和最小值；分别求这13个点相邻点横坐标和纵坐标之差绝对值的最大值和最小值； 4．对所得13个点的每个点，在其周围找16个点（可疑点），其中格宽为d=500。将这16个点与其余的12个犯罪的地点进行比较，若与某犯罪的地点距离小于等于d，则去掉该可疑点。求出所有剩余的点（构成可疑点集合）。 （5）对于d=500:500:60000，重复上述过程，对于不同的d，求可疑点的个数，画出d与可疑点个数的图形，并确定d为何值使，可疑点个数最小。 （6）对上述所得的d，类似问题4，求可疑点YD。 （7）对于YD中每个点（如第j个点），计算该点到问题2中的13个点的距离，建立该点的衰减函数：，求的最大值及YD中对应的点——罪犯所在的可能性最大的点，并将该点加在问题2的图形上。 （8）对于d= (21:30)*500，重复上述过程，对于不同的d，求罪犯所在的可能性最大的点，并将这些点加在问题2的图形上。 四 4.1 编辑m文件 disiti1 function radian=disiti1(x) d = size(x,1); for i=1:d; a=deg2rad(x(i,1)); b=deg2rad(x(i,2)/60); c=deg2rad(x(i,3)/3600); radian(i,1)=a+b+c; end 输入纬度: weidu=[53,49,7.73;53,48,28.74;53,50,28.40;53,48,38.71;53,49,4.08;53,25,57.38;53,47,59.94;53,39,15.67;53,27,40.55;53,42,41.97;53,47,21.45;53,48,30.40;53,49,20.17] 结果： disiqi1(weidu) 0.9393 0.9391 0.9397 0.9392 0.9393 0.9326 0.9390 0.9364 0.9331 0.9374 0.9388 0.9391 0.9394 输入经度为： jingdu=[1,32,30.83;1,31,58.44;1,29,32.95;1,45,49.71;1,31,55.91;2,15,3.14;1,46,19.78;1,46,46.64;2,13,26.73;1,52,26.05;1,45,47.33;1,40,16.68; 1,34,40.65] 结果： disiti1(jingdu) 0.0269 0.0268 0.0260 0.0308 0.0267 0.0393 0.0309 0.0311 0.0388 0.0327 0.0308 0.0292 0.0275 4.2 以地心为原点，赤道平面为投影平面，其中a是坐标点的纬度，b是坐标点的经度。地球半径R为6380000米。 x = Rcos(a)cos(b)y = Rcos(a)sin(b) 新建M文件： function plane = disiti2(a,b) c = size(a,1); r = 6400000; for i=1:c d(i,1) = r*cos(a(i))*cos(b(i)); d(i,2) = r*cos(a(i))*sin(b(i)); end x1=(min(d(:,1))-1250): 1250:(max(d(:,