交巡警论文数学建模.docx

1、交巡警论文数学建模2011高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为

2、（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：交巡警服务平台的设置与调度摘要在城市中警务资源分配存在不均衡的现象，和警力调度存在困难，本文建立模型解决交巡警服务平台合理调度的问题，和平均分配问

3、题。对于问题一：第一小问，首先通过路程=时间*速度，将3分钟转换为3千米，用matlab编程对A区的各个节点标号，画图并用Floyd算法求出任意两点间的路程。把得到的数据在excel中筛选出离20个交巡警服务平台小于3km的节点，并按就近原则将各个节点划分给各个交巡警服务平台，对于大于3km分配给离他最近的平台。第二小问，对进出该区的13条要道进行快速的封锁合理调度方案的问题，将问题转化成13个平台封锁13个要道的问题，并且这13个平台到这13个节点的时间最短。引入0-1规划建立模型，再用lingo求解得到，得到不同的平台封锁不同的进出口的调度方案。第三小问，对增加交警服务平台的问题，主要对现

4、有平台没办法在3分钟内到达的节点处增加，首先用Flody算法算出这些点的3km范围内的所有节点，为了使增加后的服务平台尽可能均衡，将增加的平台设在距离现有平台最远的节点上。对于问题二：第一小问，本题提出服务面积，服务人口，和案发率三个指标对影响交巡警服务平台的设置。对三个指标按公式量化后给出一个综合评价值，并判断出该市的平台设置明显不合理。为了改善和提高合理性，需要在D区增设平台。第二小问，通过MATLAB画出A区图和罪犯在3分钟内可能到达的面积区域，并画出全市的交通路线图，分析罪犯可能的逃跑路线，根据围堵范围最小，和围堵的面积最小的原则，做出最优的围堵方案。关键词：Flody算法 0-1规划

5、 lingo 调度方案一、问题重述根据某市的交巡警服务平台的相关情况，建立数学模型分析研究下面的问题：1）根据出的该市A区的交通网络图，和现有的20个交巡警服务平台的设置情况示意图，为各个交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。当发生重大事件的时候，需要全区20个交巡警服务平台，对进出该区的13条要道进行快速的封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。根据现有交巡警服务平台分配，出现警力剩余，和警力不足的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台

6、的具体个数和位置。2）针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案（参见附件）的合理性。如果有明显不合理，请给出解决方案。如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。二、问题分析针对问题一：1) 首先通过速度将时间转换为距离，对A区的各个节点标号。通过软件求出任意两点间的路程，并筛选出离20个交巡警服务平台小于3km的节点，并按就近原则将各个节点划分给各个交巡警服务平台。2) 对进出该区

7、的13条要道进行快速的封锁合理调度方案的问题，将问题转化成13个平台封锁13个要道的问题，并且这13个平台到这13个节点的时间最短。引入0-1规划建立模型，再用lingo求解得到，得到不同的平台封锁不同的进出口的调度方案。3) 对增加交警服务平台的问题，主要对现有平台没办法在3分钟内到达的节点的3km范围内增加，首先用Flody算法算出这些点的3km范围内的所有节点，通过就近原则，将增加的平台设在距离最近的节点上。针对问题二：1) 对于评价全市现有交巡警服务平台设置方案合理性的问题，本文对服务面积，服务人口，和案发率三个指标进行归一化处理，用层次分析法确定各指标的权值，得到对全市各区平台设置进

8、行综合评价，并判断现有方案的不合理性，并给出改善的方法。2) 对围堵A区32节点发生案件疑犯的问题，画出全市的交通路线图，和报案时可能到达的位置。分析罪犯可能的逃跑路线，做出最优的围堵方案。三、模型假设1) 假设所有的道路均为双向通道，2) 假设没有发生任何的交通拥堵问题3) 假设所有警车和罪犯逃跑速度均为行驶4) 假设所有案件发生在路口的节点上四、符号说明 任意两节点间的距离、 任意两节点的横坐标 任意两节点的纵坐标 第i个交巡警平台相对于第j个节点的封锁 第i个交巡警平台到第j个节点的最短路程 警车与罪犯的行驶速度 三种指标的综合结果 分别为服务面积、发案率、服务人口三个指标的权重 对各区

9、的三种指标的平均数 接到报警后罪犯逃至某一节点的时间 接到报警后交巡警到达封锁点的时间 罪犯从32号节点逃跑到某一节点的距离 交巡警从服务平台到达封锁节点的距离五、模型的建立与求解5.1问题一模型的建立：5.1.1针对该市A区的交通网络图，和现有的20个交巡警服务平台的设置情况示意图，为各个交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地的问题，首先可以算出平台管辖范围为3km。用matlab对各点标号【1】可得：图1图1：A区各个节点的标号在给各个节点标完号后用公式（5-1-1）算出两个节点间的路线的长度，并用fl

10、oyd算法【2】，算出任意两点间的最短路程。 公式（5-1-1）用excel对数据进行筛选得到离20个交巡警服务平台小于3km的节点，为保证能在最短的时间到达事发地，对节点离两个或两个以上的交巡警服务平台都小于3km的节点，按就近原则分配得到表1表1：各交巡警服务平台管辖的各个路口节点交巡警服务平台管辖的各个节点16768697173747576782404344723970354556566457606263645495051525356585967306132474883346931343545101126271225132122232414152829163637381741421880

11、81828319777920848586878889909192其中28，29，38，39，61，92六个点，任何一个交巡警都不能在3分钟内到达，因此将这六个点分给离他各自最近的服务平台。由表1，可得各个交巡警服务平台管辖节点，因此可在原该市A区的交通网络图中画出各个交巡警服务平台的管辖区域，见图2。图2：各个交巡警管辖区域5.1.2对进出该区的13条要道进行快速的封锁合理调度方案的问题。引入0-1变量，若第i个平台对第j个路口节点封锁记为1，若第i个平台对第j个路口节点不封锁记为0。设为第i个交巡警平台到第j个节点的最短路程,V表示警车的速度。则可以建立以下函数。约束条件：利用lingo编程

12、求最优解【3】可得下表2结果。表2：封锁方案交巡警服务平台标号出入A区路口节点路程2383.98214620.355482.47587300.583182910.49319161.532510227.70711243.80531212013230.51414015284.7158162110.0063由表可以看出全部封锁所用时间为10.4931分钟，封锁过程如图3所示。图3：封锁线路图5.1.3对于警力不均衡需要增加交巡警服务平台的问题。从表一中可以看出，28、29、38、39、61、92六个点，任何一个交巡警都不能在3分钟内到达，所以需要增加平台使这六个点在3分钟内有警车到达。在5.1.1中

13、的Flody算法算出任意两点间的距离，通过筛选，得到可以增加平台的点如表3表3：距离没有被覆盖的点小于3km的点没有被覆盖的点282938可以增加平台的点28、2928、2938、39、40没有被覆盖的点396192可以增加平台的点38、39、4048、6187、88、89、90、91、92为了使管辖的区域范围尽可能均衡，所以使增加的交巡警服务平台离原有的平台尽可能的远，于是给出一种增加平台的方案，增加平台的标号为：29、39、61、92。增加后的管辖区域为表4：增加平台后的管辖范围交巡警服务平台管辖的各个节点116768697173747576782240434472703354556566

14、44576062636455495051525356585966773032474888334699313435451010111126271212251313212223241414151516163637171741421818808182831919777920208485868788899091292829393839616192925.2问题二5.2.1判断现有交巡警服务平台的合理性和改进。按照要求，交巡警平台的设置主要遵循以下原则： 1、警情主导警务原则：根据管区道路交通流量、拥堵状况、治安复杂情况、发案量高低，科学确定平台管控区域； 2、快速处警原则：城区接警后确保快速到达现场(

15、本文要求三分钟内到达)；3、方便与安全原则：按照醒目、规范，方便群众和确保安全的原则，科学设置平台。 平台设置在遵循上述三大原则的基础上，应当结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况，在充分考虑现有警力和财力并确保安全的条件下，科学确定平台的数量和具体位置。结合本题提出服务面积，服务人口，和案发率三个指标对影响交巡警服务平台的设置。对三个指标按公式量化后给出一个综合评价值，判断交巡警服务平台的合理性。 （5-2-1） 表5：各区平均服务人数单位：万人区域平台数人口平均服务人口A20603B8212.625C17492.882D9738.111E15765.06

16、6F11534.811表6：各区平均服务面积单位：平方公里区域平台数面积平均管辖面积A20221.1B810312.875C1722113D938342.556E1543228.8F1127424.909表6：各区平均案发率区域平台数案发率平均案发率A20124.56.225B866.48.3C17187.211.012D967387.533E15119.47.96F11109.29.97根据三个指标的重要性用层次分析法和归一化对三个指标首先赋予它的相对尺度，服务面积的相对尺度为4服务人口的相对尺度为2案发率的相对尺度为3以此得出它们的成对比较矩阵：用MATLAB求出对比较矩阵的权重得分别为0

17、.6、0.2、0.2。带入公式（5-2-1）中得到一个综合指标。见表9.表9：各区综合指标区域综合指标A3.2302B 8.9713C10.0879D 23.2249E16.5775F15.4616平均值12.9256由表9可以明显看出D区的分配明显最不合理，同时说明，该市现有交巡警服务平台设置方案明显不合理，改进方案，D区的平台少面积大人口多，现有警力明显不足，所以应在D区增设交巡警服务平台。5.2.2 对P点发生案件疑犯的围堵方案通过MATLAB画出A区图和罪犯3分钟内能到达的面积区域。如下图4：六、模型评价与改进缺点：在分配管辖范围时本文只假设案件发生在节点位置处，具有一定片面性。优点：

18、在研究管辖范围时，采用最短路径算法，确保了在案发时，交巡警平台会以最快的速度赶到案发地。图4由图可得：罪犯从32号节点可逃至28、29、30、38、48、62号节点，通过这些节点可以逃至C区、D区或F区，因此，先通过MATLAB画出全市路线图4和节点与节点间的距离表得出模型与结果（图5、表10）：图5表9交巡警出警路线交巡警起点通往交巡警终点节点编号X轴坐标Y轴坐标X轴坐标Y轴坐标节点编号2403377.5381381624381343388.5330.540167167399131366369169210396250350239171284409298.5378235173295382296

19、37223732010134317034237032191355174340371480417312372307561六、模型的评价与改进缺点：1) 本文假设案件发生在节点位置具有片面性。2) 对于分配封锁路线和围堵的问题存在一些主观性，科学性较低。3) 模型也有一定的局限性，本文中是在理想化的条件下，如道路通畅，出警车辆和人员配备理想化。4) 5.12中对封锁路线可以进行改进，让16号巡警台往14号节点，14号巡警台到21号，7号巡警台往29号节点，8把号巡警台往30号节点。这样就可以缩短封锁时间。优点：1) 本模型较好的解决了交巡警的出警和调度问题，对现有交巡警服务平台设置方案的不合理提出

20、解决方案。2) 本文对交巡警的封锁路线和围堵罪犯的方案用路线图表示出来，能更加形象和直观的看出出警方法。 七、模型的推广对罪犯的拦截模型可以运用到其他领域，比如火警的出警，医院的急救，不安全事故的快速救援。八、参考文献1 XX知道：设置交巡警服务平台的原则和任务http: /2王莲芬, 许树柏，层次分析法引论，中国人民大学出版社, 1990附录：1MATLAB画各节点程序： Q=xlsread(E:数学建模交巡警服务平台2 (2).xls,全市交通路口节点数据);x1=Q(1:92,2);y1=Q(1:92,3);for i=1:92 plot(x1,y1,.); hold on j=num2

21、str(i); j= ,j; text(x1(i),y1(i),j);end hold on T1=x1,y1; T2=load (E:数学建模交巡警服务平台T2.txt); j=Q(1:20,2); k=Q(1:20,3); e=219 280 337 251 234 225 212 243 246 314 371 315 381; f=316 292 328 277 271 265 290 328 337 367 330 374 381; for i=1:140 x3=T1(T2(i,1),1); y3=T1(T2(i,1),2); x4=T1(T2(i,2),1); y4=T1(T2(i

22、,2),2); a=x3,x4;b=y3,y4; plot(x3,y3,B.,x4,y4,B.,j,k,RO,a,b,e,f,R*); hold on endT2数据1 751 782 443 453 654 394 635 495 506 597 327 478 98 479 3510 3411 2211 2612 2514 2115 715 3116 1416 3817 4017 4217 8118 8118 8319 7920 8621 2222 1323 1324 1324 2525 1126 2726 1027 1228 2928 1529 3030 730 4831 3231 34

23、32 3333 3433 834 935 4536 3536 3736 1636 3937 738 3938 4139 4040 241 1741 9242 4343 243 7244 345 4646 846 5547 4847 647 548 6149 5049 5350 5151 5251 5952 5653 5253 5454 5554 6355 356 5757 5857 6057 458 5960 6261 6062 462 8563 6464 6564 7665 6666 6766 7667 4467 6868 6968 7569 7069 7169 170 270 4371 7

24、271 7472 7373 7473 1874 174 8075 7676 7777 7877 1978 7979 8080 1881 8282 8382 9083 8484 8585 2086 8786 8887 8887 9288 8988 9189 2089 8489 9090 9191 922 通过MATLAB用floyd算法算任意两点间的最短路程的程序dist=1000*ones(92);for i=1:92 dist(i,i)=0;endD=xlsread(E:数学建模交巡警服务平台2 (2).xls,全市交通路口节点数据);no=D(1:92,1); x=D(1:92,2); y

25、=D(1:92,3); E=xlsread(E:数学建模交巡警服务平台2 (2).xls,全市交通路口的路线);sp=E(:,1);ep=E(:,2);for k=1:length(sp) if sp(k)=92&ep(k)=92 x1=x(sp(k); y1=y(sp(k); x2=x(ep(k); y2=y(ep(k); t=sqrt(x1-x2)2+(y1-y2)2)/10; dist(sp(k),ep(k)=t; dist(ep(k),sp(k)=t; endendfor k=1:92 for i=1:92 for j=1:92 if dist(i,k)+dist(k,j)dist(i,j) dist(i,j)=dist(i,k)+dist(k