交巡警服务平台的设置及调度-等奖论文

1、-PAGE . z.第 PAGE 2页，共 NUMPAGES 1页2012年数学建模大赛全国一等奖论文-PAGE . z.交巡警效劳平台的设置与调度摘要本文结合*城市实际情况对交巡警效劳平台的设置和调度进展了深入研究，解决了以下问题：借助于Warshall-Floyd算法得出了区任意两点间的最短路，并按照距离最近原则将各路口分派给相应的平台，得到各平台的管辖围见表1，其中，除6个节点外，其余86个节点的出警时间均小于3分钟。建立二部图的最大匹配模型解决了13条要道快速全封锁问题，得最短封锁时间约8分1秒，各平台警力调度方案如下：效劳台号4571011封锁路口4830291221以出警时间不超过

2、3分钟为首要准则分析得出需增加4个效劳平台，通过计算机搜索比较了所有可能的72种方案后，按照工作量均方差最小原则确定出新增平台位置分别为28、39、48、87号路口，此时，工作量均方差取得最小值2.3703。在引入影响巡警效劳平台设置合理性的3个指标根底上，建立熵权模糊评判模型，对平台设置合理性进展判决，得出现有平台设置不合理，其中区和区尤为明显，针对其工作量大且3km平台覆盖率低的情况提出了解决方案。证明了关于围堵的一个结论，提出了一端围堵法，确定出了为实现围堵所需要封锁的随时间变化而变化的路口集合，并将其与全城所有效劳平台构成动态二部图，根据匈牙利算法得出了在此方法下的最短围堵时间为10.

3、79分钟，需调用37个平台警力，具体围堵方案如下：效劳平台17166167168169170171172封锁路口92248252175254178182213关键词 Warshall-Floyd算法 二部图 匈牙利算法 模糊评判 -. z.一 问题的重述警察肩负着刑事执法、治安管理、交通管理、效劳群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警效劳平台。每个交巡警效劳平台的职能和警力配备根本一样。由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警效劳平台、分配各平台的管辖围、调度警务资源是警务部门面临的一个实际课题。试就*市设置交巡警效劳

4、平台的相关情况，建立数学模型分析研究下面的问题：1附件1中的附图1给出了该市中心城区A的交通网络和现有的20个交巡警效劳平台的设置情况示意图，相关的数据信息见附件2。请为各交巡警效劳平台分配管辖围，使其在所管辖的围出现突发事件时，尽量能在3分钟有交巡警警车的时速为60km/h到达事发地。对于重大突发事件，需要调度全区20个交巡警效劳平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警效劳平台警力合理的调度方案。根据现有交巡警效劳平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区再增加2至5个平台，请确定需要增加平台的具体个数

5、和位置。2针对全市主城六区A，B，C，D，E，F的具体情况，按照设置交巡警效劳平台的原则和任务，分析研究该市现有交巡警效劳平台设置方案参见附件的合理性。如果有明显不合理，请给出解决方案。如果该市地点P第32个节点处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警效劳平台警力资源的最正确围堵方案。二 模型的分析此题主要研究交巡警管辖围分配、设置与调度问题。能否及时响应一个或者多个路口的效劳请求，是此问题所迫切的要求。由于交巡警平台设置及效劳要求都在路口，自然可将问题抽象成图论模型。管辖围确实立首先依赖于各节点间的最短路，这可借助于图论中的

6、最短路算法Warshall-Floyd得以解决，进而以最快出警时间为原则进展分配。对于多个路口的同时请求效劳，需要尽快的分配各个不同平台警力到相应路口，由于要求每个平台至多效劳一个路口，将问题转化为由平台和路口所构成的二部图的匹配问题。当需要围堵犯罪嫌疑人时，首要问题是确定需要封锁的依时间变化的路口集合，进而可将其转化为动态二部图的匹配问题。在考虑已有交巡警设置方案合理性时，先结合全市的具体情况寻找决定六区各交巡警效劳平台设置合理性的一些指标，采用熵权的模糊评判方法，对平台设置合理性进展判决，并针对判决结果对设置方案进展合理的调控。三 模型假设犯罪嫌疑人和警车速度均为60km/h；效劳平台接警

7、后即可立即出警；一个效劳平台的警力最多封锁一个路口；交巡警效劳平台均设在路口；相邻节点间的道路为直线段。四 符号说明：初始距离矩阵表示路线的距离，假设路口之间无直达路线，取：最短距离矩阵：=1,2,3，20表示区的20个交巡警效劳平台：以为顶点划分为边集的二部图：二部图的匹配：嫌疑人或者警车的移动速度五 模型的建立与求解5.1 基于Warshall-Floyd算法的最短距离分配 根本思想运用Floyd算法计算城区中任意两点间的最短路程，并将每个路口交给距离最近的交巡警效劳平台管辖。 最近距离分配算法Step1：由全市交通路口的路线及路口节点坐标数据，由假设5根据勾股定理计算出初始距离矩阵程序见

8、附录1。Step2：然后依据Warshall-Floyd算法得出任意两个路口之间的最短距离矩阵程序见附录2，记其中的前20行为。Step3：对的每一列取最小值，并记录最小值大于3km的数值，设第列的最小值由第行取得，则将路口交由第个效劳平台管辖假设有两行均取得最小值则任取其一即可程序见附录3。 分配结果及分析表1 A区交巡警平台管辖围表巡警平台12345管辖的路口1、67、68、69、71、73、74、75、76、782、39、40、43、44、70、723、54、55、65、66 4、57、60、62、63、645、49、50、 51、52、53、56、58、59巡警平台678910管辖的路

9、口67、30、32、34、47、48、61、8、33、469、31、34、35、4510巡警平台1112131415管辖的路口11、26、2712、2513、21 、22 、23、 241415、28、29巡警平台1617181920管辖的路口16、36、37、3817、41、4218、80、81、82、8319、77、7920、84、85、86、87、88、89、90、91、92由于1-20号设置了巡警平台，因此由自己管辖， 28、29、38、39、61、92号路口与最近的巡警平台的距离均大于3km，分别为4.75、5.70、3.41、3.68、4.19、3.60单位：km，即无法在3min

10、到达。能在三分钟之能到达的路口节点占总结点数的93.5%。5.2 基于二部图的快速全封锁方案由假设2可知，一个巡警平台的警力最多封锁一个路口，要实现快速全封锁，就是要使13条交通要道在最短时间全部由20个巡警平台中的*13个平台一一封锁封锁时间以最后一个路口被封锁的时间计。 基于二部图的快速全封锁方案的思想对于*个时间，建立一个二部图，其中分别表示13个要道与A区的20个效劳平台。边表示平台可在时间到达要道，即。使用匈牙利算法得到的最大匹配。如果饱和，说明可从得到一个全封锁。否则，就增大时间，重新循环上述操作。 基于二部图的快速全封锁方案的实现设20个巡警平台分别到达13个交通要道的时间从小到

11、大依次是。设，由二部图边的构造知。故假设可在时间封锁即找到饱和的匹配，则一定可在时间封锁。即有下述结论：结论1：最短封锁时间一定是*个巡警平台到达*个交通要道的时间。基于上述结论，封堵时间依次取，直到对应的二部图存在饱和的匹配，即可找到最短全封锁时间。Step1：表示13个交通要道距离20个巡警平台的最短距离矩阵，记20个巡警平台分别到达13个交通要道的时间从小到大依次是.Step2：依次取,确定二部图，记其对应的邻接矩阵为，其中Step3：使用匈牙利算法得到的最大匹配，假设饱和，则当前为最短全封锁时间，算法终止，否则转Step2。 最快封锁方案及其检验由附录程序4，得最小全封锁时间=8分1秒

12、，此时对应的二部图具有一个饱和13个要道的匹配由此矩阵，得全封锁方案如下矩阵的1-13行分别表示12、14、16、21、22、23、24、28、29、30、38、48、62号路口，1-20列分别表示1-20号巡警效劳台.表2 A区对13条要道的全封锁方案效劳台号12345710封锁路口38166248302912距离5.887.394.397.403.188.027.59效劳台号111213141516封锁路口212224232814距离5.076.882.396.474.756.74经检验上表第3行，效劳台到各个对应封锁路口的距离均不超过8.02km约合8分1秒,说明所得结果是正确的。5.3

13、 A区交巡警效劳平台调整方案 交巡警平台各指标定义效劳平台工作量：指该效劳平台管辖的各路口节点日均报案率之和；*地的出警时间：指管辖该地的巡警从交巡警平台到达该地的时间。 总体调整思路交巡警效劳平台增设原则：以对辖区所有节点的快速响应为首要考虑因素，统筹兼顾各效劳平台工作量的均衡性。即平台设置的首要目标是增加尽可能少的效劳平台2-5个以实现对该区所有节点3分钟以的全覆盖，并且在新增平台后，全区所有路口节点的管辖权按就近原则重新分配后各平台工作量尽可能均衡。Step1：新增平台数目确定筛选出出警时间长于3分钟即距离管辖平台大于3千米的路口节点，称其为偏远节点。表3 长于三分钟的偏远节点与巡警平台

14、的距离偏远节点282938396192出警距离km4.75.73.43.64.13.6上述6个出警时间超过3分钟的节点可分为4组如图1所示，从最短距离矩阵中可知，不同组类节点间距离均大于6km，故新增一个节点至多能使一类节点的出警时间降到3分钟以，因此至少需要新增4个效劳平台。图1 四组偏远路口另一方面，从4组中各取一个节点作为效劳平台即可满足对这6个节点3分钟的全覆盖。初步方案可定为在28、 38、 61、92位置设立效劳平台，从新分配管辖围后工作量的均方差为3.1005。Step2：平台位置确实定总体思想：新增四个平台，在满足所有节点出警时间均不超过3分钟的前提下，使得调整后的工作量具有尽

15、可能小的均方差。记为这6个节点28、29;38、39;61;92的3km以邻域节点，因为一个效劳平台最多只能覆盖一个组类，故新增的4个节点应分别取自以及。由最短距离矩阵得，故共有种新增方案。用计算机对这些方案进展比较，选出使各效劳平台工作量方差最小的新增位置，分别为28、39、 48、87号，按5.1重新分配后工作量均方差为2.3703。图2 两种方案下不同平台的工作量图注释：方案一为根据28,38,61,92共4个点选择的分配方案结果 方案二为根据28，39，48，87共4个点选择的分配方案结果由图可知方案二的各平台工作量较方案一有明显收敛，即各平台工作量较均衡，工作负荷很大9以上和很小2以

16、下的平台数量都明显减少，但图局部平台的工作量并未改变，这是由以下两方面原因造成的：一是由于问题节点的位置原本就比较偏远，其周边的节点有限，因此新增效劳点后，能产生的影响有限；二是由于要优先满足3分钟覆盖这一限制条件，因此新增效劳台的位置被限定在一定的区域，无法全局安排，因此，只能产生局部影响。5.4基于熵权的交巡警效劳平台设置的模糊综合评价模型根本思想 先确定影响主城六区各交巡警效劳平台设置合理性的3个影响指标各区效劳平台平均的工作量、各区效劳平台节点平均覆盖率、各区效劳平台效劳人口密度，再采用熵权的模糊评判方法，对平台设置合理性进展判决，并针对判决结果对设置方案进展合理的调控。确定影响合理性

17、因素 根据设置交巡警效劳平台的原则和任务，结合现有效劳平台工作量的不均衡和局部地方出警时间过长情况，可以确定出影响交巡警效劳平台设置合理性的3项指标：效劳平台的工作量、效劳平台节点覆盖率、平台效劳人口密度。 针对全市各区的具体情况，基于全局的考虑，先将各指标具体阐述如下： 1各区效劳平台的平均工作量2各区效劳平台节点覆盖率 基于第一问Warshall- Floyd算法计算出区中各路口与最近的巡警平台的距离均大于3km的效劳平台个数为6个的方法，求出全城582个节点间的最短路矩阵，通过Matlab编程依次可以算出剩余五区各个路口与最近的巡警平台的距离均大于3km的效劳平台个数，记为平台覆盖不到的

18、路口数，则各区效劳平台节点覆盖率为： 3区效劳平台效劳人口密度 模型算法 Step1：建立影响交巡警效劳平台设置合理性的因素域 Step2：建立评判集 Step3：在影响交巡警效劳平台设置合理性的因素域与评判集之间进展隶属度分析，建立模糊关系矩阵 矩阵中表示因素域中第个因素对于等级域中第个等级的隶属度。 Step4：确定，其为3个因素对交巡警效劳平台合理性指标的权重，并满足 Step5：求出，评价结果属于中最大值对应的合理性等级。 建立评价集 效劳平台设置合理性大小是相对且模糊的，不可能定性描述，是属于模糊集理论。则根据模糊数学理论，可以将3个指标分为6个等级，分别为非常合理，较合理，合理，不

19、是很合理，不合理，明显不合理。 熵权法权重向量确实定 1判断矩阵权值 要比较3个指标、对于效劳平台设置的影响，根据上述指标计算公式，通过E*cel和Matlab编程求解得各区效劳平台平均的工作量、各区效劳平台节点平均覆盖率以及各区效劳平台效劳人口密度万人/效劳平台如下表：表4 6城区3项指标的统计结果区域各效劳平台平均工作量各效劳平台节点平均覆盖率各效劳平台效劳人密度A6.22593.75%3B8.391.78%2.625C11.01169.48%2.882D7.53376.92%8.111E7.9668.93%5.067F9.92767.59%4.818将上述6城区3项指标归一化为判断矩阵：

20、 23个评价指标的熵值 根据熵权法的公式中熵的定义公式计算出3个评价指标的熵值为： 模型的求解 交巡警效劳平台合理性模糊综合评价等于与两个矩阵的乘积， 即 Matlab编程见附录5求解得： 结合评价指标的熵值： 根据最大隶属度原则,0.195最大，所对应的是明显不合理，故全市六区的巡交警效劳平台设置不合理。解决方案 结合表4，发现区效劳平台的平均工作量为11.011发案率，是六个区中工作量最大的，而区效劳平台对道路节点的平均覆盖率为69.48%，又是六个区中非常低的，效劳人口还较少，说明区平台设置很不合理，而区效劳平台平均工作量也很大，节点平均覆盖率最低，区平台设置也不合理。即： 1从效劳节点

21、的覆盖率出发 根据前面求得六区各路口与最近的巡警平台的距离均大于3km的效劳平台个数和对应的节点，得出、区交巡警平台覆盖不到3km局部路口如下：表5 C、F交巡警平台未覆盖路口城区 该城区覆盖不到的节点名称覆盖不到总节点数 183，199，200，201，202，203，205，206，20747 486，487，505，506，507，508，509，510，51235故建议在上述孤立节点附近增设效劳平台，以增大效劳节点的覆盖率。 2从效劳平台的工作量出发 从全市六区交通网络与平台设置的示意图看出、区是在郊区，同时、两区的工作量又很大，故建议在原有效劳平台上增加值班巡警人数根底上，在所辖城区

22、边上安排机动巡逻车，进展机动巡逻。5.5基于全动态二部图的围堵方案要实现对犯罪嫌疑人的围堵，最完美的设想是将犯罪嫌疑人刚刚经过的节点和正在前往的节点分别封锁，从而将犯罪嫌疑人限定在两个有连线且已经封锁的节点间。 一种尝试的封堵方案：两端封堵法由最短距离矩阵可以确定从32节点出城所需的最短时间，确定封堵时间上限设，记表示距离32号节点车程不超过的所有路口，表示所有与中任一节点相邻的路口集合，易知，假设能在时刻封锁，则是一个可行的围堵时间。Step1：依次此取，Step2：建立二部图,其中表示全城所有的效劳平台集合，其中的边表示平台可在时间到达路口，即与的车程不超过Step3：求的最大匹配，假设饱

23、和，算法完毕，否则转Step1。结果：执行此算法未找到饱和匹配，说明此方案下不能完成围堵任务。 改进的封堵方案：一端封堵法与两端封堵法相对应，一端封堵法在于封锁，相应于两端封锁法，此方法只需封锁更少的节点。结论二： 假设交巡警可在时间全封锁，则嫌疑人无法逃离本市。证明：反证法。假设嫌疑人可以逃离本市，设其逃离本市经过的路口序列为：，其中表示本城17个出城口之一。因，显然有。记，注意到，故。由的最小性，知。因相邻路口，故，从而有. 另一方面，因包含了从32节点车程不超过的所有路口以及，故嫌疑人到达的时间大于. 这是一个矛盾，因为在时刻已被*巡警封锁，证毕。算法实现：将两端封锁法中的改为即可程序见

24、附录6。结果分析：当时，算法找到一组围堵方案如下表所示：表6 交巡警平台要封锁的路口号及其之间的距离效劳平台123451011121314封锁路口8485901861932224471459487平台与路口的距离km4.096.647.857.2210.27.713.816.45.059.64效劳平台1617166167168169170171172173封锁路口52192248252175254178182213221平台与路口的距离km10.75.488.713.814.982.226.1784.899.7效劳平台174175178179320321372475476477封锁路口2742

25、12275277370371491568530535平台与路口的距离km6.996.493.17.818.9810.18.014.655.63效劳平台478480481482483484485封锁路口544554555563528565567平台与路口的距离km2.138.145.638.2610.69.065.1检验：由表中的平台与路口的距离均不超过10.8km，结果合理。六 模型的评价与改进6.1 模型的评价1.模型的优点1在设置新的交巡警平台时，在优先考虑到出警时间的同时也统筹兼顾到各平台工作量的平衡性；2运用到了模糊评判模型，将模糊的且不易定性描述的合理性大小进展量化评判。3将封锁固定

26、路口及围堵嫌疑人的问题统一转化为依时间变化的二部图的匹配问题，借助于匈牙利算法，得以简单高效解决。2.模型的缺点过于强调出警时间3km覆盖，使得交巡警效劳平台的工作量不太均衡。6.2 模型的改进放松对出警时间的要求，建立工作量和出警时间的多目标规划，以使得管辖方案的配置更为合理。参考文献1王文波，数学建模及其根底知识详解M，：大学.2振航，数学建模，M：中国人民大学，2004.3明哲，金俊，石端银，图论及其算法M, ：机械工业，2010.4卫国，MATLAB程序设计教程M，：中国水利水电，2005.5启源，金星，数学建模案例选集M，：高等教育，2006.-. z.附录1. A区初始距离矩阵Ma

27、lLab程序function M=initdisM()load JDL* % 载入节点和路线的数据M=inf(92,92);for i=1:92 M(i,i)=0;endm,n=size(L*);for i=1:m start=L*(i,1); endd=L*(i,2); if(start=92&endd=92) M(start,endd)=sqrt(JD(start,1)-JD(endd,1)2+(JD(start,2)-JD(endd,2)2); M(endd,start)=M(start,endd); endend2最短距离矩阵functionD,R=floydwarshall()D=i

28、nitdisM();n=length(D);for(i=1:n) for(j=1:n) R(i,j)=j; endendfor(k=1:n) for(i=1:n) for(j=1:n) if(D(i,k)+D(k,j)D(i,j) D(i,j)=D(i,k)+D(k,j); R(i,j)=k; end; end; endend3.管辖围的分配load JDL*D,R=floydwarshall()D20_92=D(1:20,1:92);Min=min(D20_92)a=zeros(92,1);for i=1:92 a(i)=find(D20_92(:,i)=Min(i);enda %分配方案b

29、=zeros(20,1);for i=1:20 ctrlby_i=find(a=i); fananlv=JD(:,3); b(i)=sum(fananlv(ctrlby_i);endb %每个平台的工作量4.1匈牙利最大匹配算法%二部图的最大匹配算法匈牙利算法%A为二部图的矩阵表示%返回值M为最大匹配function M=ma*match(A)m n=size(A);M=zeros(m,n);y=zeros(1,n);%求一个极大初始分配for i=1:m for j=1:n if(A(i,j)&y(j) M(i,j)=1; y(j)=1; break; end endendwhile(1)

30、*=zeros(1,m); %0表示未标记 y=zeros(1,n); for i=1:m if(any(M(i,:) %*i是非饱和的 *(i)=-(n+1); %标记，-表示未扫描, y共有n个 end end while(1) %尝试寻找M增广链 flag=0; for i=1:m if(*(i)0) %标记但未扫描 *(i)=-*(i); %正号表已扫描 for j=1:n if(A(i,j)&y(j)=0&M(i,j)=0) y(j)=-i; flag=1;%出现新标记的y end end end end if(flag=0) break;end flag=0; for j=1:n

31、if(y(j)0&any(M(:,j) %Breakthrough:找到增广链。存在一个标记且非饱和的yj flag=1; k=y(j); M(k,j)=1; while(*(k)=n+1) %倒退求M增广链，修改M M(k,*(k)=0; M(y(*(k),*(k)=1; k=y(*(k); end break; end end if(flag=0) break;end %Non-Breakthrough end4.2基于二部图匹配的13条交通要道的封锁方案m=13;n=20;FF=1:20;%20个效劳平台YD=12 14 16 21 22 23 24 28 29 30 38 48 62;

32、%13个要道D=floydwarshall();L=D(YD,FF);l=sort(L(:);for i=1:length(l) a=l(i); A=(L=l(i)+eps) %建立二部图，其中边表示在了平台与要道距离不超过了l(i) Q=ma*match(A)if (sum(Q(:)=length(YD) % 饱和13个路口 break; %找到最小的距离时间 endenda5.模糊评判程序function E*ample8_6A=1/12,1/6,1/3,2/3,1/6,1/12;R=0.12 0.16 0.22 0.15 0.16 0.19; 0.20 0.20 0.15 0.16 0.15 0.14; 0.11 0.10 0.11 0.31 0.19 0.18 ;fuzzy_zhpj(3,A,R) %调用综合评判函数end%functionB=fuzzy_zhpj(model,A,R) %模糊综合评判B=;m,s1=size(A);s2,n=size(R);if(s1=s2) disp(A的列不等于R的行);else if(model=1) %主因素决定型 for(i=1:m) for(j=1:n) B(i,j)=0; for(k=1:s1) *=0; if(A(i,k)R(k,j) *=A(i,k)