《数学建模论文城市交通巡警平台的设置与调度.doc》由会员分享，可在线阅读，更多相关《数学建模论文城市交通巡警平台的设置与调度.doc（38页珍藏版）》请在三一办公上搜索。
1、城市交通巡警平台的设置与调度摘 要由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。本文要解决的就是某市设置交巡警服务平台设置方案，以及如何处理在确保突发事件问题。对于第一问，根据附件中的各点的坐标和图中所给的各标志点之间的相邻关系，我们求得任意两个相邻标志点的直线距离，根据附件中的全市交通路口的路线做出了邻接矩阵，再用Floyd算法求得任意两点间的最短距离。在此基础上，为了确定需要增加平台的具体个数和位置，采用主成分分析法。应用迪杰斯特拉（Dijkstra）算法进行搜索得到了该区交巡警服务平台警力合理的
2、调度方案。对于第二问，给出了设置交巡警服务平台的可量化的原则和任务，对现有方案进行评价然后进行优化；案发地点在A区，题目没有给出逃犯的车速，这里要处理好，怎样叫实现了围堵也是需要考虑的问题。关键字：邻接矩阵、距离矩阵、整数线性规划、主成分分析、surfer作图一问题的重述警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。每个交巡警服务平台的职能和警力配备基本相同。根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源。就某市设置交巡警服务平台的相关情况，建立数学模型分析研究
3、下面的问题：（1）为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。（2）针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设
4、置方案（参见附件）的合理性。如果有明显不合理，请给出解决方案。如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。二、问题的分析问题一中有三个小问题，分别讨论在现有巡警台不变的情况下，确定出每个巡警台的控制范围，要求在三分钟之内尽可能到达；当有案件发生时，各交巡警按预定的路线到达指定路口封锁该路口，要求我们给出各节点接到指示时他们的行车路线；根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。根
5、据给出的地图和其他数据，运用matlab软件使用Dijkstra算法以及floyd算法，确定出了最短路径，从而可以计算得出每个巡警台所能控制的范围。不仅仅要考虑运行路线的最短和优化性，还要考虑时间尽可能较少的优化。问题二三基本假设1.不考虑巡警在实际工作中所出现的故障而导致延误追捕。2.假设各站点的警力量是平均一致且为一固定值（巡警台人数高峰期和低潮期的平均值为单一均值）。3.在整个路途中，通过各种通讯工具，走的路程都是最短路程。4.不考虑巡警车在行驶过程中出现的塞车、抛锚等耽误时间的情况。5.不考虑警员所消耗的时间。7.在整个路途中，转弯处不需要花费时间8.由于题目没有给出逃犯的车速，假设逃
6、犯的速度应该不大于警车的时速；（1）I、要求各交巡警服务平台在其分配的管辖范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地，也即要求服务平台的管辖范围为一个圆域，半径为3公里。作图的代码见附录1，得到图1。图1 管辖范围计算每个平台到所覆盖的节点的距离。如果一个节点同时被多个平台覆盖，那么选取离其最近的平台。该区交巡警服务平台警力的调度方案如下：表1平台1平台到节点间的距离及范围（备注：单位：百米）管辖范围节点6716.19节点6812.07节点695节点7115.3节点7213节点7318.36节点746.26节点759.3节点7612.84节点786.4
7、平台2平台到节点间的距离及范围（备注：单位：百米）管辖范围节点3936.82节点4019.14节点438节点449.49节点708.6平台3平台到节点间的距离及范围（备注：单位：百米）管辖范围节点5422.71节点5512.66节点6515.24节点6618.4平台4平台到节点间的距离及范围（备注：单位：百米）管辖范围节点5718.68节点6017.39节点6152.1节点623.5节点6310.31节点6419.37平台5平台到节点间的距离及范围（备注：单位：百米）管辖范围节点4714.56节点495节点508.49节点5112.3节点5216.6节点5311.71节点5620.84节点58
8、23.86节点5913.22平台6空缺平台7平台到节点间的距离及范围（备注：单位：百米）管辖范围节点305.83节点3211.4节点4812.9平台8平台到节点间的距离及范围（备注：单位：百米）管辖范围节点338.28节点469.3平台9平台到节点间的距离及范围（备注：单位：百米）管辖范围节点3120.55节点345.02节点354.24节点455.02平台10空缺平台11平台到节点间的距离及范围（备注：单位：百米）管辖范围节点269节点2716.43平台12平台到节点间的距离及范围（备注：单位：百米）管辖范围节点2517.89平台13平台到节点间的距离及范围（备注：单位：百米）管辖范围节点2
9、127.09节点229.06节点235节点2423.85平台14空缺平台15平台到节点间的距离及范围（备注：单位：百米）管辖范围节点2842.52节点2957.01平台16平台到节点间的距离及范围（备注：单位：百米）管辖范围节点366.08节点3711.18节点3834.06平台17平台到节点间的距离及范围（备注：单位：百米）管辖范围节点418.5节点429.85平台18平台到节点间的距离及范围（备注：单位：百米）管辖范围节点808.06节点816.71节点8211.73节点835.39平台19平台到节点间的距离及范围（备注：单位：百米）管辖范围节点779.85节点794.47平台20平台到节
10、点间的距离及范围（备注：单位：百米）管辖范围节点8411.75节点854.47节点863.61节点8714.66节点8813.52节点899.49节点9013.03节点9115.99节点9236.04II、为了得到该区交巡警服务平台警力合理的调度方案，首先我们计算A区92个节点彼此之间的实际距离1、首先我们可以根据题中所给的各个标志点的坐标，用matlab计算出任意两点之间的直线距离，得到92*92的距离矩阵m：代码见附录22、根据题中交通路口的路线，我们可以得到各标志点的邻接矩阵：，即如果两个点相邻，则邻接矩阵中相对应的元素的值为1，否则为0；例如：1和2这两个点相邻，那么。代码见附录33、
11、根据Floyd算法，我们是要求出各标志点任意两两之间的实际交通距离，所以我们需要得到A区相邻两个标志点的沿公路的交通距离。我们可以利用距离矩阵的元素与的点乘积得到相邻标志点间的距离矩阵：4、 我们可以将中不相邻点间距离0改为无穷大（Inf）从而得到标志点与标志点间的权值矩阵： ，即如果1和5之间不相邻，也即不能直接到达，那么D中的和都将变成和等于无穷大（Inf），否则则等于D中相应元素的数据。5、运用Floyd算法求出任意两点间最短距离，得到最短距离矩阵： 算法见附录4，矩阵见附录5。对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁，也即要
12、封锁13条交通要道在A区内的交通路口节点，标号为12、14、16、21、22、23、24、28、29、30、38、48。该问题转化为整数线性规划问题：Min +约束：i、j、k、l、m、n、p、q、r、s、t、u20的整数且互不相同；这里应用迪杰斯特拉（Dijkstra）算法进行搜索，算法描述如下：（1）假设用带权的邻接矩阵arcs 来表示带权有向图S为1-20整数的集合，它的初始状态为20个元素的集合，初值为0。i、j、k、l、m、n、p、q、r、s、t、u这12个未知数各自按先后顺序实现从1-20的循环，每一次选择了一个1-20之间的整数，则将其从S的集合中删除。（2）D表示总共的cost
13、，每一次对未知数选择了一个整数v后，则加上以该未知数为列数，v为行数的d矩阵中相应的元素值。（3）如果后面算的的D值比前面所算得的D值小，那么选取小的那个D值（4）在操作（1）、（2）中共循环2020!/8!次。结论：i=12,j=14,k=16,l=14,m=15,n=13,p=12,q=11,r=10,s=7,t=3,u=6表2 封锁分配图节点标号121416212223242829303848对应平台1014169151312117836Cost63.6400104.685.045.0026.9327.2980.1530.6124.4421.93MiniunTotal cost469.6
14、3III、为了确定需要增加平台的具体个数和位置，这里我们采用主成分分析。 计算相关系数矩阵我们首先计算未设置为服务平台的72个节点如果设置为服务平台，那么相应的考虑了距离和发案率的工作量。从矩阵d后72列和后72行抽出得到7272的矩阵，该矩阵中的每一个元素乘以相应行数对应的发案率得到矩阵这样既考虑了平台相应的出警时间，又考虑了相应的发案率。由矩阵X可以计算出相关系数矩阵 在上式中，rij（i，j=1，2，p）为原变量的与之间的相关系数，矩阵X的第i列和第j列，其计算公式为 ，为矩阵X第i列和第j列元素的平均值，n=72 计算特征值与特征向量首先解特征方程，通常用雅可比法（Jacobi）求出特
15、征值，并使其按大小顺序排列，即0。 计算主成分贡献率及累计贡献率主成分的贡献率为累计贡献率为(i=1,2,)选取累计贡献率达8595%的特征值建议增加的平台个数为3个，分别是22、23、24。主成分分析matlab代码见附录6 （2）I、根据犯罪率和人口密度与巡警台的正比关系，可以得出，人口密度越大，犯罪率越高的地方，更应该增加巡警台的设置。首先形象地表示各区每个标点的犯罪率高低以及人口密度，从而更好地得出结论A区各结点犯罪率的标示图 B区各结点犯罪率的标示图 C区各结点犯罪率的标示图D区各结点犯罪率的标示图 E区各结点犯罪率的标示图 F区各结点犯罪率的标示图 由于犯罪率和人口密度与巡警台的正
16、比关系以及图中所表示的情况可以得出，当前设置的巡警台存在不合理性，更改的结果是：A区情况第一问已经给出，B区警力配置基本合理，不需再多做调度。C区需要增加一个，增加在273路口处。D区中应该新增8个巡警台。E区中需要增加4个巡警台F区中新增4个巡警台较为合理。 城区单位面积警力利用率各辖区管辖人口出口分配量警力分配合理度A城区0.9160B城区0.08212C城区0.08495D城区0.02734E城区0.03762F城区0.04534II、对全市交通路口的路线表建立有向图D(V,A)在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑，按假设8逃犯最多逃出3公里。 对于追捕逃犯问题，针对案发后犯罪嫌
17、疑人去向不明，我们采用圈套式方法，利用动态规划进行分析，找出人力，物力及时间达到一个平衡点。根据第一题的第二小问，我们可以计算出来，A区13个交通要道出口的每个封锁时间为t1，t2，t3，t4，t5，t6，t7，t8，t9，t10，t11，t12，t13，及用时最长的路口时间为T1和用时最短的路口的时间为T2。同时，找到从P出A区最短的线路（见图P）事实上，经过计算得出，犯罪嫌疑人只有可能在两个区中，即A区和C区，我们先考虑犯罪嫌疑人跑出A区到C区的情况一：犯罪嫌疑人由P-节点30，大约需要1.8分钟，也就是说犯罪嫌疑人在3分钟之后已经离开A区，进入C区，所以此时我们应该考虑C区巡警台的围捕问
18、题。经计算可以得出，出动173,174 号平台的警力封锁216,299号节点即可。情况二：犯罪嫌疑人还在A区，可供他选择也就是两个方向，第一小方面是往左边逃跑（如情况二图一），也就只有三种可能出项的情况，通过计算可以得出，巡警台15封锁28号路口，10平台封锁26路口，14平台封锁14路口即可。另一方面是往右边逃跑（如情况二图二），通过计算得出，2,3,4号巡警台往最近的路口处进行封堵就可以达到围捕成功。（图P）五、模型的评价与改进优点：1.、对于给出的大量的数据，我们首先对数据进行处理，将其转化为有用的数据，如：附件中的各点的坐标转化为点与点之间的距离。2、模型建立的思路简单清晰，并且可以得
19、到很好的效果。3、模型的假设很符合实际生活，以致模型可以很好的运用于相对应的实际生活中。不足：1、在第一问计算需要增加平台的具体个数和位置时，只是按照经验选取了累计贡献率达8595%的特征值。2、在第(2)问中，提出的指标体系没有理论根基。五、参考文献六、附录附录1：clc,clear all;load data1.txt; load data2.txt;c=413 403 383.5 381 339 335 317 334.5 333 282 247 219 225 280 290 337 415 432 418 444 ;d=359 343 351 377.5 376 383 362 35
20、3.5 342 325 301 316 270 292 335 328 335 371 374 394 ;for i=1:140 x1=data1(data2(i,1),1);y1=data1(data2(i,1),2); x2=data1(data2(i,2),1);y2=data1(data2(i,2),2); a=x1,x2;b=y1,y2; plot(x1,y1,*,c,d,r*,a,b); hold on endfor j = 1:20 cen=c(j), d(j) %圓心 r= 30%半徑 theta=linspace(0,2*pi); x=r*cos(theta)+cen(1);
21、 y=r*sin(theta)+cen(2); plot(x,y) axis equal end这里data1.txt为A区交通网络中路口节点的横纵坐标文件，data2.txt为A区交通网中连接两路口节点路线的起点标号和终点标号文件。附录2：dist(x,y) 求任意两点距离function m=dist(a,b)% a为横坐标，b为纵坐标,且均默认为列向量；for i=1:size(a,1) for j=1:size(b,1) m(i,j)=(a(i)-a(j)2+(b(i)-b(j)2)0.5; endend其中a为A区交通网络中92个路口节点的横坐标，b为A区交通网络中92个路口节点的纵
22、坐标。附录3load sjsxX.txtload sjsxY.txtsjX=sjsxXsjY=sjsxYN=92for i=1:Nfor j=1:NJZ1(i,j)=0;endendfor k=1:143JZ1(sjX(k),sjY(k)=1;End其中sjsxX.txt为A区交通网中连接两路口节点路线的起点标号文件，sjsxY.txt为A区交通网中连接两路口节点路线的终点标号文件附录4%floyd.m%采用floyd算法计算A区中每对顶点最短路%d是矩离矩阵%r是路由矩阵function d,r=floyd(a) n=size(a,1); d=a; for i=1:n for j=1:n r
23、(i,j)=j; end end r; for 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)=r(i,k); end end end %k % d %r end附录5矩阵dDistance = Columns 1 through 10 0 18.8680 30.5655 36.9628 75.9276 81.6088 96.0469 78.6924 81.7863 135.3403 18.8680 0 21.0772 40.9176 72.0069 78.8923 88.0738 69
24、.3001 70.0071 122.3315 30.5655 21.0772 0 26.6177 51.0416 58.1055 67.4036 49.0637 51.2957 104.7771 36.9628 40.9176 26.6177 0 42.0268 46.3276 65.8502 52.3283 59.7013 112.0591 75.9276 72.0069 51.0416 42.0268 0 8.0623 26.0768 22.9456 34.5254 76.4853 81.6088 78.8923 58.1055 46.3276 8.0623 0 27.6586 29.50
25、42 41.0488 78.5684 96.0469 88.0738 67.4036 65.8502 26.0768 27.6586 0 19.4551 25.6125 50.9313 78.6924 69.3001 49.0637 52.3283 22.9456 29.5042 19.4551 0 11.5974 59.7369 81.7863 70.0071 51.2957 59.7013 34.5254 41.0488 25.6125 11.5974 0 53.7587 135.3403 122.3315 104.7771 112.0591 76.4853 78.5684 50.9313
26、 59.7369 53.7587 0 175.8408 161.5549 145.3694 154.2992 118.6971 120.2830 92.8493 102.0417 95.2733 42.4382 198.7083 185.9704 168.1822 173.2808 134.1641 133.9589 108.2589 121.4352 116.9273 63.6396 208.0024 192.3876 177.9979 189.4525 155.6663 157.6991 130.1076 137.7044 129.7998 79.2086 148.9228 133.154
27、0 119.1354 132.3301 102.6499 106.3297 79.1770 82.1736 72.8629 33.0606 125.3196 113.2828 94.8591 100.4353 63.8905 65.7951 38.1838 48.1923 43.5660 12.8062 82.0792 67.6831 51.8773 66.2288 48.0416 55.0364 39.4462 25.6223 14.5602 55.0818 24.0832 14.4222 35.3306 54.4266 86.3539 93.2952 101.6514 82.5984 82
28、.2982 133.3754 22.4722 40.3113 52.4619 51.4125 93.1343 97.7394 115.3516 99.0581 103.1601 156.8949 15.8114 34.4384 41.4638 37.1652 79.0253 83.4865 101.7104 85.9796 90.8240 144.5579 46.7547 65.4370 74.2243 65.1249 106.5317 109.5536 130.9695 116.7497 122.5765 176.0824 181.5709 165.7106 151.7638 164.317
29、5 132.4575 135.2479 107.6151 113.2453 104.6375 57.1402 199.4618 183.6981 169.5590 181.5248 148.4924 150.8145 123.1666 130.0250 121.8277 72.2496 210.1904 194.3399 180.3282 192.3337 159.1132 161.3196 133.6899 140.7924 132.6386 82.7587 212.5135 198.2171 182.0254 190.3083 153.3786 154.2012 127.3146 137.
30、9801 131.7004 78.2624 195.1333 181.1767 164.6003 172.4014 135.3514 136.2094 109.2886 120.0771 114.0175 60.4152 167.3708 152.8823 136.9535 146.5512 111.8660 113.8640 86.2670 94.4378 87.2353 35.3836 170.9247 156.9243 140.4066 148.8036 112.8373 114.3208 86.9382 96.5000 90.0125 36.7865 172.8034 160.7016
31、 142.3701 146.6092 107.3313 107.1868 81.4371 94.9868 91.0824 39.1152 168.4429 157.1146 138.2109 140.9441 100.8464 100.1848 75.2728 90.0250 87.1436 37.9473 99.3227 92.1792 71.3180 67.8178 26.5707 26.4008 5.8310 24.5459 31.4006 52.8015 98.3260 88.3629 68.5000 71.1214 34.6554 37.7359 11.1803 19.6596 20
32、.1246 42.0119 87.0919 77.9295 57.6390 59.4243 24.6982 29.4109 11.4018 8.6313 14.7648 53.2541 86.4696 76.3217 56.5088 60.5991 28.6356 33.9559 15.6205 8.2765 10.0000 51.4782 86.5867 75.0017 56.1471 63.5138 35.2598 41.1005 22.3886 12.7769 5.0249 49.2164 79.5550 67.1193 48.9923 59.2220 37.1214 44.0114 2
33、9.8329 14.5774 4.2426 55.7853 80.9568 67.6018 50.4505 62.5879 42.1070 49.0102 33.8378 19.5576 8.5440 54.7449 85.4400 72.4431 54.8840 65.6220 41.7732 48.1664 30.4138 18.8282 7.2801 50.0100 51.0392 34.5398 24.4387 48.5412 56.0357 64.0703 62.7694 43.4108 39.8497 89.1403 49.3964 33.5261 21.9146 45.6098
34、53.6004 61.6117 61.2944 41.8629 39.0512 89.3588 37.5832 19.1442 21.1009 47.5946 67.2347 74.9567 78.1313 58.6941 56.6789 106.6419 31.5634 17.4428 36.1732 58.3095 86.8116 94.1077 100.1312 80.7976 79.3363 129.0242 16.1555 16.0312 36.1836 50.6582 86.1626 92.6121 103.5761 85.0323 86.0233 138.3112 16.1245
35、 8.0000 28.6400 45.7193 79.2023 85.8836 95.9010 77.2172 78.0064 130.2498 23.0217 9.4868 11.6297 34.0771 62.6498 69.6419 78.6448 59.9708 61.1310 113.9517 73.0068 61.0082 42.4647 52.7376 34.1321 41.5933 32.0156 13.7295 9.0000 62.3618 71.8471 61.2046 41.6083 48.9004 28.1603 35.6931 28.6531 9.3005 10.81
36、67 64.2573 88.9550 83.2166 62.1550 56.2694 14.5602 14.8661 12.8062 20.7966 31.0483 63.7024 99.1413 93.3006 72.2582 66.0927 24.0832 21.9317 12.1655 28.2931 36.7151 59.0762 72.1803 67.5426 46.5108 39.3859 5.0000 13.0384 26.9258 19.9625 31.3209 76.2168 71.7844 69.8928 49.4292 36.2802 8.4853 10.0499 34.
37、4093 30.3727 41.7612 84.9588 67.9890 66.1551 45.7739 32.6382 10.5119 13.7295 36.5308 30.4138 41.5030 86.6170 64.5600 62.1289 41.6203 30.0042 12.0416 17.0880 37.1618 28.7141 39.3573 86.4002 65.7647 60.8358 39.8026 34.0771 11.4018 19.1050 31.7805 20.5548 30.8869 79.3221 43.1856 38.5876 18.0624 18.2003
38、 33.6155 40.3113 53.0094 36.7491 42.5441 95.8541 42.4264 33.5261 12.6590 26.4622 39.4081 46.8615 54.7449 36.5034 39.5601 93.3006 60.8769 57.9828 37.4066 27.2259 15.1327 21.0238 38.8973 28.2931 38.2753 87.0919 55.2472 56.2161 37.5832 18.6815 24.8646 28.0045 50.3612 40.6602 50.4009 99.3340 62.6099 63.
39、6553 44.7018 25.8118 21.0950 22.3607 47.1699 40.3547 51.0000 97.3088 66.1287 65.0000 44.9138 30.3356 13.4164 16.0312 39.4462 32.9317 43.8634 89.4986 52.6972 56.4004 39.7398 15.9452 32.3110 34.3657 58.1378 48.7904 58.4123 107.4151 85.9069 85.6037 65.4847 49.2164 19.4165 12.0000 37.5899 41.5030 53.0377 87.8009 38.8330 43.9090 30.1040 3.5000 42.2966 46.0435 66.7608 54.0231 61.8466 113.7409 27.2029 34.1760 25.1446 10.3078 52.0096 56.5685 75.1332 60.4525 66.7308 119.9208 22.1359 25.4951 17.2409 15.9138 53.9351 59.4811 75.1066 58.8430 63.6946 117.3925
