課程設(shè)計(jì)交巡警服務(wù)平臺(tái)的設(shè)置與調(diào)度

1、交巡警服務(wù)平臺(tái)的設(shè)置與調(diào)度摘 要對(duì)提供的數(shù)據(jù)進(jìn)行了預(yù)處理后，利用DjiST算法求出任意兩個(gè)路口節(jié)點(diǎn)之間的最短路程，程序見附錄2，結(jié)果見“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編程（見附錄1）求出滿意解如圖1所示。對(duì)進(jìn)出該區(qū)的13條交通要道快速封鎖問題，我們把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道最短路程之和的最小化作為第一目標(biāo)，把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道的最長(zhǎng)

2、路程最小化作為第二目標(biāo)，建立多目標(biāo)規(guī)劃模型，利用matlab編程（見附錄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í)間過長(zhǎng)，不能及時(shí)趕到突發(fā)事件地點(diǎn)，并且使各服務(wù)平臺(tái)的工作量盡量均衡。我們計(jì)算了B區(qū)的出警時(shí)間方差和工作量方差如表3所示（程序見附錄6）。建議在99、94、9

3、8處進(jìn)行增設(shè)服務(wù)平臺(tái)，另外也可以考慮在95處增加警力。其他各區(qū)的處理方法同B區(qū)。該市P點(diǎn)發(fā)生了重大刑事案件，快速圍堵嫌疑犯的調(diào)度問題，是一個(gè)動(dòng)態(tài)規(guī)劃。我們給出了嫌疑犯逃跑各階段的路口集合，可能的逃跑路線如圖4所示。用matlab編程（程序見附錄7），得到最佳圍堵方案如表4。關(guān)鍵字：動(dòng)態(tài)規(guī)劃 多目標(biāo)規(guī)劃 最短路徑算法 圖論 0-1規(guī)劃一、 問題的重述“有困難找警察”，是家喻戶曉的一句流行語。警察肩負(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ù)部門面臨的一個(gè)實(shí)際課題。試就某市設(shè)置交巡警服務(wù)平臺(tái)的相關(guān)情況，建立數(shù)學(xué)模型分析研究下面的問題：（1）根據(jù)給出的相關(guān)數(shù)據(jù)信息，為各交巡警服務(wù)平臺(tái)分配管轄范圍，使其在所管轄的范圍內(nèi)出現(xiàn)突發(fā)事件時(shí)，盡量能在3分鐘內(nèi)有交巡警（警車的時(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í)間過長(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(xiàn)）的具體情況，按照設(shè)置交巡警服務(wù)平臺(tái)的原則和任務(wù)，分析研究該市現(xiàn)有交巡警服務(wù)平臺(tái)設(shè)置方案的合理性。如果有明顯不合理，請(qǐng)給出解決方案。如果該市地點(diǎn)某點(diǎn)處發(fā)生了重大刑事案件，在案發(fā)3分鐘后接到報(bào)警，犯罪嫌疑人已駕車逃跑。為了快速搜捕嫌疑犯，給出調(diào)度全市交巡警服務(wù)平臺(tái)警力資源的最佳圍堵方案。二、 問題的分析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ā)案率。由此可見該問題是一個(gè)多目標(biāo)規(guī)劃問題。為了求解問題，我們首先對(duì)提供的數(shù)據(jù)進(jìn)行一些預(yù)處理，再利用DjiST算法求出任意兩個(gè)路口節(jié)點(diǎn)之間的最短路程，程序見附錄二，結(jié)果見“A區(qū)節(jié)點(diǎn)間最小路的值.xls”。由于題目提供的節(jié)點(diǎn)坐標(biāo)單位為mm，因此我們將警車3分鐘走過的路程根據(jù)比例化為30mm。如該區(qū)發(fā)生重大事件，交巡警對(duì)進(jìn)出該區(qū)的13條交通要道要實(shí)現(xiàn)快速全封鎖，即是求各交巡警服務(wù)平臺(tái)到13條交通要道最小路徑的總體最小化，以滿足盡量實(shí)現(xiàn)快速封鎖的目的。也是一個(gè)多目標(biāo)規(guī)劃問題。增設(shè)交巡警

7、服務(wù)平臺(tái)的原則是：（1）對(duì)有些出警時(shí)間過長(zhǎng)的地方，可以考慮增設(shè)交巡警服務(wù)平臺(tái)，以期盡量能在3分鐘內(nèi)有交巡警到達(dá)事發(fā)地。（2）有些工作量過大的交巡警服務(wù)平臺(tái)附近，可以考慮增設(shè)交巡警服務(wù)平臺(tái)，以期工作量應(yīng)盡量均衡。2、對(duì)該市現(xiàn)有交巡警服務(wù)平臺(tái)設(shè)置方案合理性的分析，主要是看目前的方案中是否存在有些服務(wù)平臺(tái)工作量過大和出警時(shí)間過長(zhǎng)的情況，如果有明顯不合理，應(yīng)考慮增設(shè)服務(wù)平臺(tái)或?qū)ぷ髁窟^大的服務(wù)平臺(tái)增加警力。該市P點(diǎn)發(fā)生了重大刑事案件，在案發(fā)3分鐘后接到報(bào)警，犯罪嫌疑人已駕車逃跑。為了快速圍堵嫌疑犯，給出最佳調(diào)度方案，是一個(gè)動(dòng)態(tài)規(guī)劃問題。 三、 模型假設(shè)1、任意相鄰兩個(gè)路口之間的道路都是直線。2、增設(shè)服

8、務(wù)平臺(tái)的時(shí)，其他服務(wù)平臺(tái)地點(diǎn)不變。3、不考慮交通堵塞，犯罪嫌疑人的車速與交巡警的車速相同，均為常量。4、各交巡警服務(wù)平臺(tái)的警力配備基本一致。四、 符號(hào)說明：表示從第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ù)平臺(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編程（見附錄1）求出滿意解：如圖1所示：圖1 A區(qū)交巡警服務(wù)平臺(tái)分配管轄范圍示意圖5.2 交巡警服務(wù)平臺(tái)警力合理的調(diào)度模型對(duì)于該區(qū)發(fā)生重大事件，交巡警對(duì)13個(gè)交通要道快速全封鎖，根據(jù)問題分析，我們把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道最

10、短路程的總和最小化作為第一目標(biāo)，把參與封堵各交巡警服務(wù)平臺(tái)到各出口要道的最長(zhǎng)路程最小化作為第二目標(biāo)，建立多目標(biāo)規(guī)劃，模型如下： 其中 利用matlab編程（見附錄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ū)域的劃分，通過計(jì)算得出A區(qū)各交巡警服務(wù)平臺(tái)的出警時(shí)間方差和工作量方差（見附錄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í)間過長(zhǎng)，服務(wù)平臺(tái)7出警時(shí)間過長(zhǎng)且任務(wù)過重，服務(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í)間過長(zhǎng)，不能及時(shí)趕到突發(fā)事件地點(diǎn)，并且使各服務(wù)平臺(tái)的工作量盡量均衡。我們計(jì)算了B區(qū)的出警時(shí)間方差和工作量方差（程序見附錄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í)間過長(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)警，犯罪嫌疑人已駕車逃跑。為了快速搜捕嫌疑犯，調(diào)度全市交巡警服務(wù)平臺(tái)警力資源的最佳圍堵方案如下：案發(fā)點(diǎn)為，以為始點(diǎn)的道路之終點(diǎn)

14、為，所有點(diǎn)的集合為稱之為狀態(tài)，為第一步可能的圍堵的路口集合，同理（集合元素的排序按距P的最短路程升序排列）稱之為狀態(tài)為第步可能的圍堵路口集合狀態(tài)轉(zhuǎn)移規(guī)律：1、在第步可能圍堵路口集合中，按升序依次判斷犯罪嫌疑人若從此路口逃跑是否能夠及時(shí)堵?。僭O(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)度方案。通過編程計(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可能逃跑的路線如圖4所示，不同的線型對(duì)應(yīng)著不同狀態(tài)下可能逃跑的路線及到達(dá)的路口節(jié)點(diǎn)。按照上述思路利用matlab編程計(jì)算的結(jié)果如表4（程序見附錄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可能逃跑的路線圖六、

16、模型的評(píng)價(jià)與推廣總的來說，模型與算法的思路清晰，其結(jié)果可幫助交巡警的調(diào)度和最佳行進(jìn)路徑的選擇，減少了出警時(shí)間，可操作性強(qiáng)，有一定的實(shí)際應(yīng)用價(jià)值，模型與算法也有一定的局限性和不足，如個(gè)別地區(qū)的出警時(shí)間超過3分鐘，也沒有考慮道路的暢通性等不定因素對(duì)出警時(shí)間的影響等.本模型與算法較好的解決了交巡警平臺(tái)的最優(yōu)選址問題，當(dāng)事故發(fā)生時(shí)，交巡警可以第一時(shí)間到達(dá)事發(fā)地點(diǎn)，有效的改善了交巡警在執(zhí)行任務(wù)中的效率，并且該模型與算法也可運(yùn)用到其他應(yīng)急問題中，比如關(guān)于消防救援問題、重大生產(chǎn)安全事故應(yīng)急救援問題等.七、 參考文獻(xiàn)【1】 趙靜,但琦等編著，數(shù)學(xué)建模與數(shù)學(xué)實(shí)驗(yàn)（第3版）， 北京，高等教育出版社，2008。【2

17、】 云舟工作室編著，MATLAB6數(shù)學(xué)建?；A(chǔ)教程，北京，人民郵電出版社，2001。【3】 吳翊，吳孟達(dá)，成禮智編著，數(shù)學(xué)建模的理論與實(shí)踐，長(zhǎng)沙，國防科技大學(xué)出版社，1999?！?】 譚永基，俞文吡編著，數(shù)學(xué)模型，上海，復(fù)旦大學(xué)出版社，1997?！?】 王海英等編著 圖論算法及其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