




版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、交巡警服務(wù)平臺(tái)的設(shè)置與調(diào)度摘 要對(duì)提供的數(shù)據(jù)進(jìn)行了預(yù)處理后,利用djist算法求出任意兩個(gè)路口節(jié)點(diǎn)之間的最短路程,程序見(jiàn)附錄2,結(jié)果見(jiàn)“a區(qū)節(jié)點(diǎn)間最小路的值.xls”。交巡警服務(wù)平臺(tái)分配管轄范圍的原則是(1)所管轄的節(jié)點(diǎn)到該交巡警服務(wù)平臺(tái)的最近路程應(yīng)盡量小于3000m,(2)各交巡警服務(wù)平臺(tái)管轄范圍內(nèi)的所有節(jié)點(diǎn)發(fā)案率之和應(yīng)盡量接近平均發(fā)案率6.225。由此建立多目標(biāo)規(guī)劃模型。對(duì)該模型利用matlab編程(見(jiàn)附錄1)求出滿(mǎn)意解如圖1所示。對(duì)進(jìn)出該區(qū)的13條交通要道快速封鎖問(wèn)題,我們把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道最短路程之和的最小化作為第一目標(biāo),把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道的最長(zhǎng)
2、路程最小化作為第二目標(biāo),建立多目標(biāo)規(guī)劃模型,利用matlab編程(見(jiàn)附錄4)求解,結(jié)果如表1和圖2所示。根據(jù)a區(qū)各交巡警服務(wù)平臺(tái)的出警時(shí)間方差和工作量方差,考慮在服務(wù)平臺(tái)15附近的a區(qū)出口節(jié)點(diǎn)29增設(shè)一個(gè)服務(wù)平臺(tái),在服務(wù)平臺(tái)7附近的路口節(jié)點(diǎn)48增設(shè)一個(gè)服務(wù)平臺(tái),在路口節(jié)點(diǎn)91增設(shè)一個(gè)服務(wù)平臺(tái)。增設(shè)3個(gè)服務(wù)平臺(tái)后的a區(qū)巡警服務(wù)平臺(tái)分配管轄范圍如圖3所示,合理性有較大改善。對(duì)全市六區(qū)交巡警服務(wù)平臺(tái)設(shè)置合理性評(píng)價(jià)如a區(qū),主要要避免服務(wù)平臺(tái)的出警時(shí)間過(guò)長(zhǎng),不能及時(shí)趕到突發(fā)事件地點(diǎn),并且使各服務(wù)平臺(tái)的工作量盡量均衡。我們計(jì)算了b區(qū)的出警時(shí)間方差和工作量方差如表3所示(程序見(jiàn)附錄6)。建議在99、94、9
3、8處進(jìn)行增設(shè)服務(wù)平臺(tái),另外也可以考慮在95處增加警力。其他各區(qū)的處理方法同b區(qū)。該市p點(diǎn)發(fā)生了重大刑事案件,快速?lài)孪右煞傅恼{(diào)度問(wèn)題,是一個(gè)動(dòng)態(tài)規(guī)劃。我們給出了嫌疑犯逃跑各階段的路口集合,可能的逃跑路線(xiàn)如圖4所示。用matlab編程(程序見(jiàn)附錄7),得到最佳圍堵方案如表4。關(guān)鍵字:動(dòng)態(tài)規(guī)劃 多目標(biāo)規(guī)劃 最短路徑算法 圖論 0-1規(guī)劃一、 問(wèn)題的重述“有困難找警察”,是家喻戶(hù)曉的一句流行語(yǔ)。警察肩負(fù)著刑事執(zhí)法、治安管理、交通管理、服務(wù)群眾四大職能。為了更有效地貫徹實(shí)施這些職能,需要在市區(qū)的一些交通要道和重要部位設(shè)置交巡警服務(wù)平臺(tái)。每個(gè)交巡警服務(wù)平臺(tái)的職能和警力配備基本相同。由于警務(wù)資源是有限的,
4、如何根據(jù)城市的實(shí)際情況與需求合理地設(shè)置交巡警服務(wù)平臺(tái)、分配各平臺(tái)的管轄范圍、調(diào)度警務(wù)資源是警務(wù)部門(mén)面臨的一個(gè)實(shí)際課題。試就某市設(shè)置交巡警服務(wù)平臺(tái)的相關(guān)情況,建立數(shù)學(xué)模型分析研究下面的問(wèn)題:(1)根據(jù)給出的相關(guān)數(shù)據(jù)信息,為各交巡警服務(wù)平臺(tái)分配管轄范圍,使其在所管轄的范圍內(nèi)出現(xiàn)突發(fā)事件時(shí),盡量能在3分鐘內(nèi)有交巡警(警車(chē)的時(shí)速為60km/h)到達(dá)事發(fā)地。對(duì)于重大突發(fā)事件,需要調(diào)度全區(qū)20個(gè)交巡警服務(wù)平臺(tái)的警力資源,對(duì)進(jìn)出該區(qū)的13條交通要道實(shí)現(xiàn)快速全封鎖。實(shí)際中一個(gè)平臺(tái)的警力最多封鎖一個(gè)路口,給出該區(qū)交巡警服務(wù)平臺(tái)警力合理的調(diào)度方案。根據(jù)現(xiàn)有交巡警服務(wù)平臺(tái)的工作量不均衡和有些地方出警時(shí)間過(guò)長(zhǎng)的實(shí)際情
5、況,擬在該區(qū)內(nèi)再增加2至5個(gè)平臺(tái),請(qǐng)確定需要增加平臺(tái)的具體個(gè)數(shù)和位置。(2)針對(duì)全市(主城六區(qū)a,b,c,d,e,f)的具體情況,按照設(shè)置交巡警服務(wù)平臺(tái)的原則和任務(wù),分析研究該市現(xiàn)有交巡警服務(wù)平臺(tái)設(shè)置方案的合理性。如果有明顯不合理,請(qǐng)給出解決方案。如果該市地點(diǎn)某點(diǎn)處發(fā)生了重大刑事案件,在案發(fā)3分鐘后接到報(bào)警,犯罪嫌疑人已駕車(chē)逃跑。為了快速搜捕嫌疑犯,給出調(diào)度全市交巡警服務(wù)平臺(tái)警力資源的最佳圍堵方案。二、 問(wèn)題的分析1、交巡警服務(wù)平臺(tái)分配管轄范圍的基本原則應(yīng)該是(1)所管轄的范圍內(nèi)出現(xiàn)突發(fā)事件時(shí),應(yīng)盡量能在3分鐘內(nèi)有交巡警到達(dá)事發(fā)地,其實(shí)就是所管轄的節(jié)點(diǎn)到該交巡警服務(wù)平臺(tái)的最近路程應(yīng)盡量小于30
6、00m,(2)各交巡警服務(wù)平臺(tái)的工作量應(yīng)盡量均衡,就是各交巡警服務(wù)平臺(tái)管轄范圍內(nèi)的所有節(jié)點(diǎn)發(fā)案率之和應(yīng)盡量接近平均的發(fā)案率。由此可見(jiàn)該問(wèn)題是一個(gè)多目標(biāo)規(guī)劃問(wèn)題。為了求解問(wèn)題,我們首先對(duì)提供的數(shù)據(jù)進(jìn)行一些預(yù)處理,再利用djist算法求出任意兩個(gè)路口節(jié)點(diǎn)之間的最短路程,程序見(jiàn)附錄二,結(jié)果見(jiàn)“a區(qū)節(jié)點(diǎn)間最小路的值.xls”。由于題目提供的節(jié)點(diǎn)坐標(biāo)單位為mm,因此我們將警車(chē)3分鐘走過(guò)的路程根據(jù)比例化為30mm。如該區(qū)發(fā)生重大事件,交巡警對(duì)進(jìn)出該區(qū)的13條交通要道要實(shí)現(xiàn)快速全封鎖,即是求各交巡警服務(wù)平臺(tái)到13條交通要道最小路徑的總體最小化,以滿(mǎn)足盡量實(shí)現(xiàn)快速封鎖的目的。也是一個(gè)多目標(biāo)規(guī)劃問(wèn)題。增設(shè)交巡警
7、服務(wù)平臺(tái)的原則是:(1)對(duì)有些出警時(shí)間過(guò)長(zhǎng)的地方,可以考慮增設(shè)交巡警服務(wù)平臺(tái),以期盡量能在3分鐘內(nèi)有交巡警到達(dá)事發(fā)地。(2)有些工作量過(guò)大的交巡警服務(wù)平臺(tái)附近,可以考慮增設(shè)交巡警服務(wù)平臺(tái),以期工作量應(yīng)盡量均衡。2、對(duì)該市現(xiàn)有交巡警服務(wù)平臺(tái)設(shè)置方案合理性的分析,主要是看目前的方案中是否存在有些服務(wù)平臺(tái)工作量過(guò)大和出警時(shí)間過(guò)長(zhǎng)的情況,如果有明顯不合理,應(yīng)考慮增設(shè)服務(wù)平臺(tái)或?qū)ぷ髁窟^(guò)大的服務(wù)平臺(tái)增加警力。該市p點(diǎn)發(fā)生了重大刑事案件,在案發(fā)3分鐘后接到報(bào)警,犯罪嫌疑人已駕車(chē)逃跑。為了快速?lài)孪右煞?,給出最佳調(diào)度方案,是一個(gè)動(dòng)態(tài)規(guī)劃問(wèn)題。 三、 模型假設(shè)1、任意相鄰兩個(gè)路口之間的道路都是直線(xiàn)。2、增設(shè)服
8、務(wù)平臺(tái)的時(shí),其他服務(wù)平臺(tái)地點(diǎn)不變。3、不考慮交通堵塞,犯罪嫌疑人的車(chē)速與交巡警的車(chē)速相同,均為常量。4、各交巡警服務(wù)平臺(tái)的警力配備基本一致。四、 符號(hào)說(shuō)明:表示從第i個(gè)服務(wù)平臺(tái)到第j個(gè)路口節(jié)點(diǎn)的最短路程。:表示從第i個(gè)服務(wù)平臺(tái)到第j個(gè)路口節(jié)點(diǎn)最短路程的正偏差。:表示第j個(gè)路口節(jié)點(diǎn)是否屬于第i個(gè)服務(wù)平臺(tái),是為1否為0。:表示第i個(gè)服務(wù)平臺(tái)是否負(fù)責(zé)圍堵第j個(gè)交通要道,是為1否為0。:表示第j個(gè)路口節(jié)點(diǎn)的發(fā)案率。:a區(qū)平均每個(gè)服務(wù)平臺(tái)的發(fā)案率 :表示從第k個(gè)服務(wù)平臺(tái)到第j個(gè)交通要道的最短路程。:表示第k個(gè)服務(wù)平臺(tái)是否負(fù)責(zé)封堵第j個(gè)交通要道,是為1否為0。五、 模型的建立和求解5.1 a區(qū)交巡警服務(wù)平
9、臺(tái)分配管轄范圍的優(yōu)化模型根據(jù)問(wèn)題分析,交巡警服務(wù)平臺(tái)分配管轄范圍的原則主要是:(1)所管轄的節(jié)點(diǎn)到該交巡警服務(wù)平臺(tái)的最近路程應(yīng)盡量小于30mm,所以將總的正偏差最小化作為第一目標(biāo)。即最小化 其中:=(2)各交巡警服務(wù)平臺(tái)的工作量應(yīng)盡量均衡,因此將服務(wù)平臺(tái)所管轄的路口節(jié)點(diǎn)的總案發(fā)率的方差最小化作為第二目標(biāo)。即最小化建立如下的多目標(biāo)規(guī)劃模型:對(duì)該模型利用matlab編程(見(jiàn)附錄1)求出滿(mǎn)意解:如圖1所示:圖1 a區(qū)交巡警服務(wù)平臺(tái)分配管轄范圍示意圖5.2 交巡警服務(wù)平臺(tái)警力合理的調(diào)度模型對(duì)于該區(qū)發(fā)生重大事件,交巡警對(duì)13個(gè)交通要道快速全封鎖,根據(jù)問(wèn)題分析,我們把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道最
10、短路程的總和最小化作為第一目標(biāo),把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道的最長(zhǎng)路程最小化作為第二目標(biāo),建立多目標(biāo)規(guī)劃,模型如下: 其中 利用matlab編程(見(jiàn)附錄4)求解,結(jié)果如表1:表1 a區(qū)交巡警服務(wù)平臺(tái)對(duì)13條交通要道全封鎖方案出口節(jié)點(diǎn)標(biāo)號(hào)12222362483029281421241638服務(wù)平臺(tái)標(biāo)號(hào)12111345871514910162 圖2為示意圖,箭頭起點(diǎn)為服務(wù)平臺(tái),終點(diǎn)為封鎖的出口節(jié)點(diǎn)。圖2 a區(qū)交巡警服務(wù)平臺(tái)對(duì)13條交通要道全封鎖方案示意圖5.3 a區(qū)交巡警服務(wù)平臺(tái)增設(shè)方案根據(jù)5.1交巡警管轄區(qū)域的劃分,通過(guò)計(jì)算得出a區(qū)各交巡警服務(wù)平臺(tái)的出警時(shí)間方差和工作量方差(見(jiàn)附錄5)
11、,結(jié)果如下。表2 a區(qū)各交巡警服務(wù)平臺(tái)的出警時(shí)間方差和工作量方差出警時(shí)間方差服務(wù)平臺(tái)標(biāo)號(hào)工作量方差服務(wù)平臺(tái)標(biāo)號(hào)010203060409050100601108012090140100160110190120.0380130.03150140.1440171.6310184.7330194.731816.47164.732036.15205.181346.5426.635141.66710.73171036.181522.807從表1中發(fā)現(xiàn)服務(wù)平臺(tái)15出警時(shí)間過(guò)長(zhǎng),服務(wù)平臺(tái)7出警時(shí)間過(guò)長(zhǎng)且任務(wù)過(guò)重,服務(wù)平臺(tái)20出警時(shí)間較長(zhǎng)且任務(wù)較重,因此考慮在服務(wù)平臺(tái)15附近的a區(qū)出口節(jié)點(diǎn)29增設(shè)一個(gè)服務(wù)平臺(tái),
12、在服務(wù)平臺(tái)7附近的路口節(jié)點(diǎn)48增設(shè)一個(gè)服務(wù)平臺(tái),在18和20之間的路口節(jié)點(diǎn)91增設(shè)一個(gè)服務(wù)平臺(tái)。增設(shè)三個(gè)服務(wù)平臺(tái)后的a區(qū)巡警服務(wù)平臺(tái)分配管轄范圍如圖3所示。圖3 a區(qū)增設(shè)后的巡警服務(wù)平臺(tái)分配管轄范圍示意圖5.4 全市交巡警服務(wù)平臺(tái)設(shè)置合理性評(píng)價(jià)及改進(jìn)方案對(duì)全市六區(qū)交巡警服務(wù)平臺(tái)設(shè)置合理性評(píng)價(jià)的依據(jù)如5.1,主要要避免服務(wù)平臺(tái)的出警時(shí)間過(guò)長(zhǎng),不能及時(shí)趕到突發(fā)事件地點(diǎn),并且使各服務(wù)平臺(tái)的工作量盡量均衡。我們計(jì)算了b區(qū)的出警時(shí)間方差和工作量方差(程序見(jiàn)附錄6)。表3 b區(qū)交巡警服務(wù)平臺(tái)設(shè)置評(píng)價(jià)指標(biāo)出警時(shí)間方差服務(wù)平臺(tái)標(biāo)號(hào)工作量方差服務(wù)平臺(tái)標(biāo)號(hào)09309309509709709909801000100
13、0.25968.15964.419421.48945.2995229.569914.4498從表3中可以看出服務(wù)平臺(tái)99的出警時(shí)間過(guò)長(zhǎng),服務(wù)平臺(tái)94和96的出警時(shí)間較長(zhǎng),服務(wù)平臺(tái)98的工作量很大,服務(wù)平臺(tái)95和94的工作量也較大。因此可以采用5.3的方法在99、94、98處進(jìn)行增設(shè)服務(wù)平臺(tái),另外也可以考慮在95處增加警力。其他各區(qū)的處理方法同b區(qū)。5.5 全市交巡警服務(wù)平臺(tái)警力資源的最佳圍堵嫌疑犯方案對(duì)于該市p點(diǎn)(第32個(gè)節(jié)點(diǎn))處發(fā)生了重大刑事案件,在案發(fā)3分鐘后接到報(bào)警,犯罪嫌疑人已駕車(chē)逃跑。為了快速搜捕嫌疑犯,調(diào)度全市交巡警服務(wù)平臺(tái)警力資源的最佳圍堵方案如下:案發(fā)點(diǎn)為,以為始點(diǎn)的道路之終點(diǎn)
14、為,所有點(diǎn)的集合為稱(chēng)之為狀態(tài),為第一步可能的圍堵的路口集合,同理(集合元素的排序按距p的最短路程升序排列)稱(chēng)之為狀態(tài)為第步可能的圍堵路口集合狀態(tài)轉(zhuǎn)移規(guī)律:1、在第步可能?chē)侣房诩现?,按升序依次判斷犯罪嫌疑人若從此路口逃跑是否能夠及時(shí)堵住(假設(shè)只要交巡警先到路口設(shè)卡,一定能夠堵住,到達(dá)時(shí)立即可以實(shí)施圍堵,設(shè)卡時(shí)間不計(jì))。如果能堵住則繼續(xù)判斷下一個(gè)點(diǎn),如果堵不住,則以該點(diǎn)為始點(diǎn),向外擴(kuò)展一步,即以該點(diǎn)為始點(diǎn)在逃逸方向上找終點(diǎn),將若干個(gè)終點(diǎn)代替該點(diǎn),得到狀態(tài)的圍堵點(diǎn)集合。2、如果某狀態(tài)中所有路口都能堵住,則可輸出最佳調(diào)度方案。通過(guò)編程計(jì)算第一狀態(tài)的路口集合為7,31,33,第二狀態(tài)的路口集合為8,
15、 34,37,47,15,30,第三狀態(tài)的路口集合為5,6,9,10,28,29,36,46,48,237,第四狀態(tài)的路口集合為16,26,35,39,49,45,49,50, 55,59,61,235,236,238,239,247,370,371,嫌疑犯從案發(fā)地點(diǎn)p可能逃跑的路線(xiàn)如圖4所示,不同的線(xiàn)型對(duì)應(yīng)著不同狀態(tài)下可能逃跑的路線(xiàn)及到達(dá)的路口節(jié)點(diǎn)。按照上述思路利用matlab編程計(jì)算的結(jié)果如表4(程序見(jiàn)附錄7):表4 全市交巡警服務(wù)平臺(tái)圍堵嫌疑犯方案圍堵節(jié)點(diǎn)標(biāo)號(hào)569102829364648237服務(wù)平臺(tái)標(biāo)號(hào)56910151711687173p圖4 嫌疑犯從案發(fā)地點(diǎn)p可能逃跑的路線(xiàn)圖六、
16、模型的評(píng)價(jià)與推廣總的來(lái)說(shuō),模型與算法的思路清晰,其結(jié)果可幫助交巡警的調(diào)度和最佳行進(jìn)路徑的選擇,減少了出警時(shí)間,可操作性強(qiáng),有一定的實(shí)際應(yīng)用價(jià)值,模型與算法也有一定的局限性和不足,如個(gè)別地區(qū)的出警時(shí)間超過(guò)3分鐘,也沒(méi)有考慮道路的暢通性等不定因素對(duì)出警時(shí)間的影響等.本模型與算法較好的解決了交巡警平臺(tái)的最優(yōu)選址問(wèn)題,當(dāng)事故發(fā)生時(shí),交巡警可以第一時(shí)間到達(dá)事發(fā)地點(diǎn),有效的改善了交巡警在執(zhí)行任務(wù)中的效率,并且該模型與算法也可運(yùn)用到其他應(yīng)急問(wèn)題中,比如關(guān)于消防救援問(wèn)題、重大生產(chǎn)安全事故應(yīng)急救援問(wèn)題等.七、 參考文獻(xiàn)【1】 趙靜,但琦等編著,數(shù)學(xué)建模與數(shù)學(xué)實(shí)驗(yàn)(第3版), 北京,高等教育出版社,2008?!?
17、】 云舟工作室編著,matlab6數(shù)學(xué)建模基礎(chǔ)教程,北京,人民郵電出版社,2001?!?】 吳翊,吳孟達(dá),成禮智編著,數(shù)學(xué)建模的理論與實(shí)踐,長(zhǎng)沙,國(guó)防科技大學(xué)出版社,1999?!?】 譚永基,俞文吡編著,數(shù)學(xué)模型,上海,復(fù)旦大學(xué)出版社,1997。【5】 王海英等編著 圖論算法及其matlab實(shí)現(xiàn) 北京航空航天大學(xué)出版社 2010【6】 王桂平,王衍,任嘉辰編著 圖論算法理論、實(shí)現(xiàn)及應(yīng)用 北京大學(xué)出版社 2011附件附錄1%a區(qū)交巡警服務(wù)平臺(tái)分配管轄范圍求解clear all;mln = xlsread(d:b,2,a2:b138); %matrix of line between nodesm
18、n = xlsread(d:b,1,a2:e93); %matrix of nodes% create 0-1 matrix of line beteen nodesmln01=zeros(92,92);mlnsize=size(mln);mlnr=mlnsize(1); %row number of mlnfor i=1:mlnr if mln(i,1)=92 & mln(i,2)=92 % only deal with the data in the part a mln01(mln(i,1),mln(i,2)=1; endend% solve the path between the n
19、odesmlnpath=zeros(92,92);m=10000; %set m with a big numberfor i=1:92 for j=1:92 if mln01(i,j)=1 mlnpath(i,j)=sqrt(mn(i,2)-mn(j,2)2+(mn(i,3)-mn(j,3)2); end % set infinite number to the nodes with not connection endendsp=shortpath(mlnpath);for h=1:1000 randservicenode=randperm(20); handlenode=zeros(92
20、,1); %storage the depart result,0 note not servicenode=1:20; for i=1:20 startsn=servicenode(randservicenode(i); handlenode(startsn)=startsn; totalwork=0; clear totalpath startnodea totalpath(1)=0; startnodea(1,:)=startsn,0; startnodecount=1; while totalwork=10 startnode=startnodea(startnodecount,1);
21、 for j=1:mlnr if mln(j,1)=startnode & handlenode(mln(j,2)=0 tempm=startnodea(:,1); tempc=find(tempm=startnode); tempp=mlnpath(startnode,mln(j,2)+totalpath(tempc); if tempp30 totalpath(length(totalpath)+1)=mlnpath(startnode,mln(j,2)+totalpath(tempc); startnodea(size(startnodea,1)+1,:)=mln(j,2),startn
22、ode; handlenode(mln(j,2)=startsn; totalwork=totalwork+mn(mln(j,2),5); end end end startnodecount=startnodecount+1; if size(startnodea,1)sp(uanode,j) minservice=j; mindistance=sp(uanode,j); end end handlenode(uanode)=minservice; uam(i,2)=minservice; uam(i,3)=mindistance; end totalwork=0; largem=0; fo
23、r i=1:20 fservices=find(handlenode=i); work=0; for j=1:length(fservices) work=work+mn(fservices(j),5); end if work6.225 totalwork=totalwork+(work-6.225)2; largem=largem+1; end end finaltotalwork(h)=totalwork; finalhandlenodeh=handlenode; finallargem(h)=largem;enda,v=min(finaltotalwork)finallargem(v)
24、附錄2%求最短路徑函數(shù)function sp=shortpath(a)m=10000;a=a+a;for i=1:length(a) for j=1:length(a) if a(i,j)=0 a(i,j)=m; end endendfor i=1:length(a) pb(1:length(a)=0; pb(i)=1; d(1:length(a)=m; d(i)=0; temp=i; while sum(pb)length(a) tb=find(pb=0); d(tb)=min(d(tb),d(temp)+a(temp,tb); tmpb=find(d(tb)=min(d(tb); temp
25、=tb(tmpb(1); pb(temp)=1; end sp(i,:)=d;end附錄3%尋找a區(qū)出口節(jié)點(diǎn)clear all;mln = xlsread(d:b,2,a2:b928); %matrix of line between nodesmn = xlsread(d:b,1,a2:e93); %matrix of nodes% create 0-1 matrix of line beteen nodesmln01=zeros(92,92);mlnsize=size(mln);mlnr=mlnsize(1); %row number of mlnfor i=1:mlnr if mln(i
26、,1)=92 & mln(i,2)=92 % only deal with the data in the part a mln01(mln(i,1),mln(i,2)=1; endend% solve the path between the nodesmlnpath=zeros(92,92);m=10000; %set m with a big numberfor i=1:92 for j=1:92 if mln01(i,j)=1 mlnpath(i,j)=sqrt(mn(i,2)-mn(j,2)2+(mn(i,3)-mn(j,3)2); end % set infinite number
27、 to the nodes with not connection endendsp=shortpath(mlnpath);j=1;crossnode=zeros(13,1);for i=1:mlnr if (mln(i,1)92) | (mln(i,1)92 & mln(i,2)=92) if mln(i,1)=92 if isempty(find(crossnode=mln(i,1) crossnode(j)=mln(i,1); j=j+1; end else if isempty(find(crossnode=mln(i,2) crossnode(j)=mln(i,2); j=j+1;
28、end end end end附錄4%尋找封鎖交通要道的程序clear all;mln = xlsread(d:b,2,a2:b143); %matrix of line between nodesmn = xlsread(d:b,1,a2:e93); %matrix of nodes% create 0-1 matrix of line beteen nodesmln01=zeros(92,92);mlnsize=size(mln);mlnr=mlnsize(1); %row number of mlnfor i=1:mlnr if mln(i,1)=92 & mln(i,2)=92 % o
29、nly deal with the data in the part a mln01(mln(i,1),mln(i,2)=1; endend% solve the path between the nodesmlnpath=zeros(92,92);m=10000; %set m with a big numberfor i=1:92 for j=1:92 if mln01(i,j)=1 mlnpath(i,j)=sqrt(mn(i,2)-mn(j,2)2+(mn(i,3)-mn(j,3)2); end % set infinite number to the nodes with not c
30、onnection endendsp=shortpath(mlnpath);load out13.matsp13=zeros(13,20);for i=1:length(crossnode) tempsp13(i,:)=sp(crossnode(i),1:20); endsave shortpath.mat sp sp13for i=1:13 for j=1:20 if tempsp13(i,j)=0 tempsp13(i,j)=m; end endendfor h=1:1000 randoutnode=randperm(13); sp13=tempsp13; clear serviceton
31、ode servicetonode=zeros(13,3); % the first column store the out node,second store service node,third store distance servicetonode(:,1)=crossnode; % for i=1:13% if crossnode(randoutnode(i)=20% servicetonode(randoutnode(i),2)=crossnode(randoutnode(i);% servicetonode(randoutnode(i),3)=0; % end% end for
32、 i=1:13 if servicetonode(randoutnode(i),2)=0 a,b=min(sp13(randoutnode(i),:),2); tempm=servicetonode(:,2); while isempty(find(tempm=b) sp13(randoutnode(i),b)=m; a,b=min(sp13(randoutnode(i),:),2); end servicetonode(randoutnode(i),2)=b; servicetonode(randoutnode(i),3)=a; end end sumd=sum(servicetonode(
33、:,3); finalstnh=servicetonode; finalsd(h)=sumd;enda,v=min(finalsd)附錄5%a區(qū)交巡警服務(wù)平臺(tái)設(shè)置合理性指標(biāo)計(jì)算clear all;mln = xlsread(d:b,2,a2:b141); %matrix of line between nodesmn = xlsread(d:b,1,a2:e93); %matrix of nodes% create 0-1 matrix of line beteen nodesmln01=zeros(92,92);mlnsize=size(mln);mlnr=mlnsize(1); %row
34、number of mlnfor i=1:mlnr if mln(i,1)=92 & mln(i,2)30 path=path+(sp(fservices(j),i)-30)2; end end totalpath(i)=path; if work6.225 totalwork(i)=(work-6.225)2; largem=largem+1; endendsorttw,positionw=sort(totalwork)sorttp,positionp=sort(totalpath)附錄6%b區(qū)交巡警服務(wù)平臺(tái)設(shè)置合理性指標(biāo)計(jì)算clear all;load partb.mat % partb.
35、mat sp mn mln sn b131handlenode=minhn;len=length(mn);lensn=length(sn);avework=sum(mn(:,5)/lensn;totalwork=zeros(lensn,1);totalpath=zeros(lensn,1);largem=0;for i=1:lensn fservices=find(handlenode=sn(i); work=0; path=0; for j=1:length(fservices) work=work+mn(fservices(j),5); if sp(fservices(j),i)30 path=path+(sp(fservices(j),i)-30)2; end end totalpath(i)=path; if workavework totalwork(i)=(work-avework)2; largem=large
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 廊坊師范學(xué)院《古代漢語(yǔ)(下)》2023-2024學(xué)年第二學(xué)期期末試卷
- 石家莊幼兒師范高等專(zhuān)科學(xué)?!痘A(chǔ)俄語(yǔ)Ⅰ》2023-2024學(xué)年第一學(xué)期期末試卷
- 2024-2025學(xué)年福建省華安一中、長(zhǎng)泰一中等四校高三語(yǔ)文試題下學(xué)期一模預(yù)考試題含解析
- 上海音樂(lè)學(xué)院《社會(huì)保障概論》2023-2024學(xué)年第二學(xué)期期末試卷
- 湖南食品藥品職業(yè)學(xué)院《現(xiàn)代表面工程技術(shù)》2023-2024學(xué)年第二學(xué)期期末試卷
- 班會(huì)主題:禁毒教育
- 福州市重點(diǎn)中學(xué)2025屆中考英語(yǔ)試題模擬試題含答案
- 2025屆河北省師大附中高三年級(jí)第二學(xué)期期末生物試題試卷含解析
- 漳州輔警考試題庫(kù)2024
- 湖南網(wǎng)絡(luò)工程職業(yè)學(xué)院《中醫(yī)傷科學(xué)》2023-2024學(xué)年第二學(xué)期期末試卷
- 2025屆山西省高三一模地理試題(原卷版+解析版)
- 2024年電信銷(xiāo)售員工年終總結(jié)
- 2025年度執(zhí)業(yè)藥師職務(wù)聘用協(xié)議模板
- Unit3 Weather Part A(說(shuō)課稿)-2023-2024學(xué)年人教PEP版英語(yǔ)四年級(jí)下冊(cè)
- 2-山東工業(yè)技師學(xué)院申報(bào)國(guó)家級(jí)高技能人才培訓(xùn)基地項(xiàng)目申報(bào)書(shū)
- 常用消毒劑的分類(lèi)、配制及使用課件演示幻燈片
- GB 45069-2024懸崖秋千安全技術(shù)要求
- 員工反恐怖協(xié)議
- 2025年高考政治一輪復(fù)習(xí)知識(shí)清單必修四《哲學(xué)與文化》重難點(diǎn)知識(shí)
- 2021年4月17日江蘇事業(yè)單位考試《綜合知識(shí)和能力素質(zhì)》(管理崗客觀題)
- 《臨床技術(shù)操作規(guī)范-放射醫(yī)學(xué)檢查技術(shù)分冊(cè)》
評(píng)論
0/150
提交評(píng)論