確定罪犯藏身地點(diǎn)問(wèn)題(已修改)
2025-06-19 21:26 本頁(yè)面
【正文】 4．確定罪犯藏身地點(diǎn)問(wèn)題已知罪犯系列犯罪的地點(diǎn)坐標(biāo)如下：PlaceDimensionLongitude153176。4939。N1176。3239。W253176。4839。N1176。3139。W353176。5039。N1176。2939。W453176。4839。N1176。4539。W553176。4939。N1176。3139。W653176。2539。N2176。1539。W753176。4739。N1176。4639。W853176。3939。N1176。4639。W953176。2739。N2176。1339。W1053176。4239。N1176。5239。W1153176。4739。N1176。4539。W1253176。4839。N1176。4039。W1353176。4939。N1176。3439。W（1）將角度制的坐標(biāo)轉(zhuǎn)換為弧度制坐標(biāo)；（2）將系列犯罪的地點(diǎn)（13個(gè)）的球面坐標(biāo)轉(zhuǎn)換為平面坐標(biāo)；并畫(huà)出這些點(diǎn)的圖形（要求帶適當(dāng)寬度的網(wǎng)格，*10^3，列寬h=*10^3）。（3）求這13個(gè)點(diǎn)之間距離的最大值和最小值；分別求這13個(gè)點(diǎn)相鄰點(diǎn)橫坐標(biāo)和縱坐標(biāo)之差絕對(duì)值的最大值和最小值；4．對(duì)所得13個(gè)點(diǎn)的每個(gè)點(diǎn)，在其周?chē)?6個(gè)點(diǎn)（可疑點(diǎn)），其中格寬為d=500。將這16個(gè)點(diǎn)與其余的12個(gè)犯罪的地點(diǎn)進(jìn)行比較，若與某犯罪的地點(diǎn)距離小于等于d，則去掉該可疑點(diǎn)。求出所有剩余的點(diǎn)（構(gòu)成可疑點(diǎn)集合）。（5）對(duì)于d=500:500:60000，重復(fù)上述過(guò)程，對(duì)于不同的d，求可疑點(diǎn)的個(gè)數(shù)，畫(huà)出d與可疑點(diǎn)個(gè)數(shù)的圖形，并確定d為何值使，可疑點(diǎn)個(gè)數(shù)最小。（6）對(duì)上述所得的d，類似問(wèn)題4，求可疑點(diǎn)YD。（7）對(duì)于YD中每個(gè)點(diǎn)（如第j個(gè)點(diǎn)），計(jì)算該點(diǎn)到問(wèn)題2中的13個(gè)點(diǎn)的距離，建立該點(diǎn)的衰減函數(shù)：，求的最大值及YD中對(duì)應(yīng)的點(diǎn)——罪犯所在的可能性最大的點(diǎn)，并將該點(diǎn)加在問(wèn)題2的圖形上。（8）對(duì)于d= (21:30)*500，重復(fù)上述過(guò)程，對(duì)于不同的d，求罪犯所在的可能性最大的點(diǎn)，并將這些點(diǎn)加在問(wèn)題2的圖形上。四 編輯m文件 disiti1function 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,。53,48,。53,50,。53,48,。53,49,。53,25,。53,47,。53,39,。53,27,。53,42,。53,47,。53,48,。53,49,]結(jié)果：disiqi1(weidu) 輸入經(jīng)度為：jingdu=[1,32,。1,31,。1,29,。1,45,。1,31,。2,15,。1,46,。1,46,。2,13,。1,52,。1,45,。1,40,。 1,34,]結(jié)果：disiti1(jingdu) 以地心為原點(diǎn)，赤道平面為投影平面，其中a是坐標(biāo)點(diǎn)的緯度，b是坐標(biāo)點(diǎn)的經(jīng)度。地球半徑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))。endx1=(min(d(:,1))1250): 1250:(max(d(:,1))+1250)。y1=(min(d(:,2))2500): 2500:(max(d(:,2))+2500)。[X1,Y1] = meshgrid(x1,y1)。hold on。plot(X1,Y1,39。b39。)。plot(X139。,Y139。,39。b39。)。plot(d(:,1),d(:,2),39。o39。)hold off。vpa(d)輸入：a =[ ]b =[ ]運(yùn)行m文件函數(shù)：disiti2(a,b)結(jié)果為：ans = [ , ][ , ][ , ][ , ][ ,
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1