1、摘要本文就某市的实际情况与需求，合理的建立了有关交巡警服务平台设置与调度的模型，通过图论模型、规划模型以及计算机程序的结合，对题中所述问题进行了求解，获得了比较满意的结果。对于问题一，首先将出警时间的约束转换为距离约束，分别利用朴素的覆盖点集以及微变量逐次调控的方法，得出了该问题的预分配方案以及最后的优化方案。问题一的第二个子问题是匹配问题，我们通过 0-1 变量 match 来标记每个交巡警服务平台是否参与道路的封锁。则当封锁路口时，需要全部 13 个路口全部封锁才能达到目的。警车到达节点所花费的时间应该以最后一个到达对应的节点的警车所需要的时间来决定。于是借助 MATLAB 程序来解决此问
2、题，由上文可得目标函数为： ，得到一个最佳的调度方案，1maxcost(,)(,)nijijmatchij封堵完成的最短时间约为 8 分钟。第三个子问题作为对第一个子问题的优化补充，我们基于前面的结果对各站点的工作量进行尽量的均衡分配，根据再分配的结果，在满足各个服务站点工作量平衡的前提下，得出结论分别要在编号为 42,57,62,90 的四个交点上添加 4 个新的服务站。对于问题二，采取与问题一的第三个子问题相同的原则和任务要求，针对全市现有的服务平台数量进行平台管辖范围的分配，分配原理与问题一的第一个子问题基本类似，只是数据规模的一个扩大问题。而在后来的模型优化过程中，我们引入人口密度的因
3、素，对现有服务平台管辖范围进行重新分配。并利用工作量的均衡性来度量设置方案合理性。最后，对于问题二的最后一个子问题，我们建立了最佳围堵方案模型。考虑到警力资源的限制，不可能完全将该区域的任何一个点都围堵住,这时就需要在原有围堵的基础上改进方案，在那些未被围住的点继续以 3 分钟的圈往外延伸。最后得出合理的围堵方案。关键词：0-1 规划 变量微调 工作量均衡度一、问题重述警察肩负着刑事执法，治安管理，交通管理，服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。如何根据城市的实际情况与需求合理地设置交巡警服务平台，分配各平台的管辖范围，调度警务资
4、源是警务部门面临的一个实际课题。试就某市设置交巡警服务平台的相关情况，建立数学模型分析研究下面的问题：（1）附件 1 中的附图 1 给出了该市中心城区 A 的交通网络和现有的 20 个交巡警服务平台的设置情况示意图。请为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在 3 分钟内有交巡警到达事发地。对于重大突发事件，需要调度全区 20 个交巡警服务平台的警力资源，对进出该区的 13 条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内
5、再增加 2 至 5 个平台，请确定需要增加平台的具体个数和位置。（2）针对全市的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案的合理性。如果有明显不合理，请给出解决方案。如果该市地点 P 处发生了重大刑事案件，在案发 3 分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。二、问题假设（1） 每个交巡警平台的职能与警力的配置相同；（2） 不考虑警车启动与停止，路上受到阻碍所花费的时间；（3） 假设巡警都按最短路径到达各案发路口；（4） 假设犯罪案件都在路口上发生；（5） 道路均为双向；（6） 假设城
6、区内道路无限速；（7） 假设犯罪案件不会在同一时间内发生多起；（8） 交巡警服务平台的节点的犯罪案件的解决不花费时间；（9） 假设犯罪车辆逃跑与警车的追赶速度相同；三、符号说明s1：交巡警管辖距离；v1：警车的平均速度；v2：嫌疑人逃跑的平均速度；t0：警车到达案发路口的时间限制；li：各点到管辖它的交巡警平台的距离；：各点的发案率；Ai：各点的工作量；zi：各节点到管辖它的交巡警服务平台的距离；xi：第 i 个节点的横坐标；yi：第 i 各节点的纵坐标；d：两节点之间的距离；cost（i，j）：i，j 两点的实际最短距离；四、模型预处理（1）交巡警服务平台的管辖范围：该问题要求在道路交点出现
7、突发事件时，交巡警尽量能在 3 分钟内到达事发点。由于警车的时速均衡且为 60km/h，所以可以将时间限制转换为距离限制，由于 ，可求出交巡警在时间限制内管辖范围的最大半径为 3km。为了10svt处理方便，我们使得交巡警服务平台以一整段路为标准来管辖各路段。同时根据附件 2 所给内容，基于各个路口的发案率，假设案件均发生在交叉路口。由此，即可将路段管理转化为对路口的管辖。当突发事件发生时，警车立即出动至所管辖的案发点。（2）交巡警服务平台布置的合理性：该问题主要考虑的是交巡警服务平台的工作量均衡问题，由于 A 区不同地域节点的密度不同，而且发案率也不相同，所以工作量可以表示为 li 与 的乘
8、积，即 。对每个交巡警服务点所管辖的范围以交点为单位，按照工作iAl量均衡的原则去逐个改变交点所属的辖区，最终使各个交巡警平台的工作量达到最优均衡状态。五、问题一的解决问题 1.1由于两节点间的距离公式为： ，借助 C 语言程序22()()ijijdxy（见附录 1） ，可得到各交巡警服务平台（共 20 个，编号 1-20）到各节点（共72 个，编号 21-92）的距离小于 3km 的各节点编号，得到下表交巡警服务平台编号可管辖的节点编号1 42 43 44 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 832 40 4
9、1 42 43 44 64 65 66 67 68 69 70 71 72 73 74 75 76 77 783 38 39 40 43 44 54 55 63 64 65 66 67 68 69 70 75 764 54 55 56 57 58 60 62 63 64 65 66 67 68 75 76 775 30 32 33 46 47 48 49 50 51 52 53 56 57 58 59 616 30 32 47 48 49 50 51 52 53 56 57 58 59 617 30 31 32 33 34 35 46 47 48 498 30 31 32 33 34 35 3
10、6 37 45 46 47 48 49 52 53 569 31 32 33 34 35 36 37 45 4610 无11 21 25 26 27 2812 24 25 2813 21 22 23 2414 2615 3116 32 33 34 35 36 37 45 4617 40 41 42 43 44 67 68 69 70 71 72 73 74 7818 42 68 69 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 9219 63 64 65 66 67 68 69 70 71 72 73 74 75
11、 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 9220 80 81 82 83 84 85 86 87 88 89 90 91表 1对于表 1，有以下四种特殊情况：1.同一个节点被分到不同的交巡警服务平台管辖范围之内；2.有些节点没有被交巡警服务平台所管辖；3.有些交巡警服务平台没有管理任何节点；4.有些距离交巡警服务平台 3km 以内的节点实际路程大于 3km。对于 1，在处理时，先按照每个交巡警服务平台所管辖的节点的多少将其进行升序排列，优先处理管辖的节点较少的平台，每处理完一个平台之后，删除后面平台所管辖的与之重复的节点，然后重新按照每
12、个平台所管辖的交点的多少将其进行升序排列，重复进行判断和删减，直到所有重复的点都被删除。对于 2，将没有被交巡警服务平台所管辖的节点划到离其最近的交巡警服务平台的管辖范围内。对于 3，暂不考虑，在问题 1 的第三小问中解决。对于 4，在解决 1 之前，结合附录 1 中的图 1，将距离某个交巡警服务平台实际路程大于 3km 的节点在表 1 中将其编号删除。根据上述方法把表 1 处理后，得到表 2，其中 28,29,38,39,61,92 没有任何交巡警服务平台可以在规定的半径内管辖它们，所以将这六个点划归到离它们最近的交巡警服务平台管辖。交巡警服务平台编号 可管辖的节点编号1 无2 38 69
13、71 73 74 75 783 28 29 44 54 55 67 68 764 57 58 60 62 63 64 65 665 49 536 47 48 50 51 52 56 597 30 328 无9 无10 无11 26 2712 25 6113 21 22 23 2414 3915 3116 34 35 36 37 45 4617 77 79 8018 无19 无20 81 82 83 84 85 86 87 88 89 90 91 92表二这样分配管辖范围都可以使每个交巡警服务平台的警车在有事故时能够尽快到达事发地点。但忽略了每个服务站工作负担的因素，导致某些工作站工作压力过大，
14、而且有的交巡警服务平台没有工作机会，这个问题将会在第三小问中解决。问题 1.2抽象为一个 20 对 13 的匹配问题来求解，其中对于每个交巡警服务平台是否参与道路的封锁用户 0-1 变量 match 来标记（若第 i 个交巡警服务平台负责第j 个出口处的堵截，则 match(i,j)=1，反之则为 0） 。要求各节点警力到达封锁路口所花费的时间应尽可能短。当封锁路口时，需要全部 13 个路口全部封锁才能达到目的。警车到达节点所花费的时间应该以最后一个到达对应的节点的警车所需要的时间来决定。借助 MATLAB 程序（见附录 2）来解决此问题，由上文可得目标函数为： 1maxcost(,)(,)n
15、ijijmatchij约束条件为：1.每一个服务站最多出动去一个路口（i，j 的组合唯一） ；2.每个通往其他区的路口都要有相应的服务站的警车去封锁；3.各节点警力到达封锁路口所花费的时间应尽可能短。根据上面 C 语言程序的结果，结合附录，最优条件下的匹配方式为：（3,16） ， （4,48） ， （5,30） ， （7,29） ， （10,12） ， （11,24） ， （12,23） ， （13,22） ，（14,21） ， （15,28） ， （16,14） ， （19,38） ， （20,62） 。其中所求的目标函数结果为8015.457。所以结论为：1.2.6.8.9.17.18 号
16、平台的警车队没有任务，3 号平台的警车队去 16 号节点，4 号平台的警车队去 48 号节点，5 号平台的警车队去 30 号节点，7 号平台的警车队去 48 号节点，10 号平台的警车队去 12 号节点，11 号平台的警车队去 24 号节点，12 号平台的警车队去 23 号节点，13 号平台的警车队去22 号节点，14 号平台的警车队去 21 号节点，15 号平台的警车队去 28 号节点，16 号平台的警车队去 14 号节点，19 号平台的警车队去 38 号节点，20 号平台的警车队去 62 号节点；封堵完成的最短时间约为 8 分钟。问题 1.3要解决该问题，首先要解决问题1.1中的工作量不均
17、衡问题，问题1.1中只是在分配管辖任务时对各个节点进行了分配，而现在就要根据不同平台的工作量对A区20个平台的管辖范围实行重分配，总体思想就是在保持较优解得前提下让工作量小的服务平台分担工作量大的服务平台所管辖的路口。分成三步进行：1.对于总工作量较小的服务平台，找出与该平台直接连通的各节点，在3km范围内，将工作量较大的节点划归到现有服务平台的管辖范围中；2.对于工作量较大的的服务平台，将距离服务中心较远的节点划分到与此节点直接联通且工作量较小的服务平台的管辖范围中。3.上两步预处理好后，根据模型预处理的（2）中的公式，算出各交巡警服务平台的工作量，调节各交巡警服务平台所管辖的节点，使各交巡
18、警服务平台的工作量达到最优均衡。经过上述三步的操作，得到表三：交巡警服务平台编号 可管辖的节点编号1 69 77 79 802 71 73 74 75 783 44 55 67 68 764 57 60 62 63 64 65 665 49 53 546 50 51 52 56 58 597 30 32 33 47 48 618 35 36 37 45 469 31 3410 2611 2212 2713 23 24 2514 2115 28 2916 38 39 4017 41 42 43 70 7218 84 87 88 89 90 9119 81 82 8320 85 86 92表三随后
19、，在原来20个服务平台的基础上再添加新的服务平台。添加服务平台的过程分成三步进行：1.找出工作量较大的交巡警服务平台，找出其管辖的节点；2.根据附件1的图1，在这些节点附近找出若干个节点，使得这些节点与其他节点的实际距离小于3km；3.根据问题1.1和问题1.2的步骤，以交巡警服务平台与节点距离3km之内，各交巡警服务平台的工作量尽量均衡为原则，得到新的各交巡警服务平台所管辖的节点。以上三步进行完毕后，可得到共要添加4个交巡警服务平台，且添加的交巡警服务平台的节点的编号为42,57,62,90。六、问题二的解决问题2.1研究交巡警服务平台设置的合理性：本题有两个判别合理性的原则：1.警车能在3
20、min内赶到案发节点；2.交巡警服务平台的工作量均衡度尽量小。对1的合理性判断：全城剩下的5个区，均与问题1.1同理，求出各区中交巡警3min内不能赶到的节点编号，如下表四：区编号 交巡警3min内不能赶到的节点编号A 28 29 38 39 61 92B 122 123 124 151 152 153C 183 199 200 201 202 203 205 206 207 208 209 210 215 238 239 240247 248 251 252 253 257 259 261 262 263 264 268 269 285 286 287 288 299 300 301 302
21、 303 304 312 313 314 315 316 317 318 319D 329 330 331 332 336 337 339 344 362 369 370 371E 387 388 389 390 391 392 393 395 407 408 409 411 412 413 417 418 419 420 438 439 443 445 446 451 452 455 458 459 464 469 471 474 486 487F 505 506 507 508 509 510 512 513 514 515 516 517 518 519 522 523 524 525
22、526 527 529 533 540 541 559 560 561 566 569 574 575 578 580表四由表四可得，共有138个节点交巡警不能在3min内赶到，数量太多。所以按原则1，交巡警服务平台设置是不合理的。对2的合理性的判断：根据模型预处理的（2）中的公式 及问题1中的方法，可以求得各iAl区的各交巡警服务平台的工作量A i如下表五所示：A B C D E F编号 Ai 编号 Ai 编号 Ai 编号 Ai 编号 Ai 编号 Ai1 10.3 93 2.1 166 3.8 320 8.7 372 5.2 475 13.1 2 9.7 94 11.3 167 8.3 32
23、1 12.0 373 4.1 476 13 3 5.6 95 9.5 168 4.7 322 4.4 374 5.5 477 10.7 4 17.1 96 11.5 169 3.4 323 4.2 375 6.1 478 9.5 5 9.7 97 5.6 170 12.9 324 7.9 376 2.6 479 8.7 6 2.5 98 12.1 171 12.4 325 2.2 377 4.2 480 4.7 7 40.4 99 4.3 172 8.3 326 5.1 378 2.6 481 7.2 8 5.0 100 4.5 173 11.5 327 7.6 379 7.4 482 4.4
24、 9 8.2 174 10.1 328 6.7 380 2.5 483 3.3 10 1.6 175 8.7 381 6.2 484 3.8 11 4.6 176 8.1 382 10.3 485 3.3 12 10.3 177 2.2 383 10.013 39.6 178 4.5 384 8.3 14 7.2 179 13.0 385 9.1 15 12.9 180 13.0 386 6.3 16 28.4 181 6.2 17 5.3 182 12.2 18 6.119 3.420 11.5表五由表五可得，交巡警服务平台工作量不均衡，其中工作量最大的是7号，工作量为40.4；最小的是10
25、号，工作量为1.6。所以按原则2，交巡警服务平台设置是不合理的。综上述，现有的交巡警服务平台的设置情况是较不合理的。交巡警服务平台设置的优化方案：在不改变现有交巡警服务平台的位置的情况下，适当增加交巡警服务平台的数量，从而使城区中没有警车不能在3min内赶到案发现场的节点且每个交巡警服务平台的工作量尽量均衡。由上面中的计算结果可知，全城区共有 138 个交巡警服务平台的警车不能在 3min 内赶到案发现场的节点。同样利用求解问题 1 中第三小问的方法与步骤，得到新增加的最小交巡警服务台数目与其位置，从而得到优化。问题 2.2根据要求，在市区P点发生案件3分钟后，警局才接到报警，随而立即采取行动
26、。假设在接到报案那一刻，嫌疑人刚要逃离案发点，此时，假设嫌疑人以远离犯罪现场的逃离方式向四周随机逃跑。由于3分钟内嫌疑人逃跑的最大距离为 。以此可以圈定其逃跑范围，而交巡警在接到报警去围堵嫌疑人时又需20vt要花费3分钟的时间。也就是说从开始逃跑到交巡警到达围堵地点的过程中，嫌疑人共有6分钟的逃跑时间，逃跑的距离为 。在警车出动围堵路口的过程20vt中，嫌疑人在以逃离 的基础上继续逃跑。为了实现成功围堵，要将嫌疑人所有20vt可能经过的路口都堵住。有未被堵住的路口的话，用同样的方法，以那些未被堵住的路口为中心，再去圈定一个三分钟的罪犯逃离范围。此时就有路口三分钟圈重叠的情况。将这些个圈定的路口
27、再以同样的方法调动警力，去尽量封锁嫌疑人可能逃跑所经过的路口。从而实现成功围堵。七、模型评价本文缺点：1.问题解决时的步骤过于繁琐，数据处理的工作量太大，导致结果不精确；2.问题解决时，由于水平有限，把实际问题做了较大的简化，可能导致与实际情况不符；3.本文没有对人口密度的情况进行综合考虑。在城区，在设置交巡警服务平台时应该对居民区附近重点布置，使交巡警能在最快的时间内赶到，以保证城市居民的安全；4.由于时间仓促，没有对问题 2.2 进行具体的求解，只是给出了一种可行的方法。本文优点：所建模型解决了交巡警服务平台的的出警问题，追捕嫌疑人的封堵路口问题。除此之外，模型还可用于消防车的安排问题，事
28、故的救援问题等，对实际情况有较好的参考作用。八、参考文献【1】http:/ C 语言描述，北京：高等教育出版社，2011 年【3】屈婉玲，耿素云，张立昂，离散数学，高等教育出版社，2008 年【4】单锋，朱丽梅，田贺民，数学模型，国防工业出版社，2012 年附录附录 1：#include#include#include#define MAXSIZE 20#define MAXJIEDIAN 92typedef structfloat xMAXSIZE;float yMAXSIZE;float last;pingtai;typedef structfloat xMAXJIEDIAN;float
29、yMAXJIEDIAN;float last;jiedian;void main()int i;int j;int m=0;int t9220=0;pingtai M;float xpt20=413,403,383.5,381,339,335,317,334.5,333,282,247,219,225,280,290,337,415,432,418,444;float ypt20=359,343,351,377.5,376,383,362,353.5,342,325,301,316,270,292,335,328,335,371,374,394;for(i=0;i3acs(m,1)=i;acs
30、(m,2)=asj(i);acs(m,3)=afw(i);m=m+1;endendclear m;clear i;%计算二十个平台距离 13 个出入口的距离for i=1:20for j=1:13Dis2(i,j)=Dis(i,acr(j,2);endendclear i;clear j;%转到 lingo 求解 A 封锁方案%求 A 每个平台的日工作量agzl=zeros(20,1);for i=1:92agzl(afw(1,i),1)=agzl(afw(1,i),1)+jd(i,6);endclear i;time=Dis./10./.*60;%question2（其中读取数据的表格见附件）
