交巡警平台分配问题

发布时间:   来源:文档文库   
字号:
如有帮助欢迎下载支持
2015西安航空学院数学建模模拟

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写)B 我们的参赛报名号为(如果赛区设置报名号的话)XXX 所属学校(请填写完整的全名):西安航空学院 参赛队员 (打印并签名 1. 栾天
2. 王辉 3. 李阳 指导教师或指导教师组负责人 (打印并签名
日期: 2015 8 17
1
如有帮助欢迎下载支持
交巡警服务平台的设置与调度
摘要
本文基于交巡警服务平台的设置与调度问题,针对交巡警服务平台管辖范围、警力调度方案、服务平台设置方案及围堵方案建立了动态规划模型、线性规划模型、利用MATLABLINGO等数学软件以及FloydDijkstra等算法解决了上述问题。 针对问题一,在解决交巡警服务平台管辖范围问题时,我们建立了动态规划模型,运用Floyd算法计算每个交巡警服务平台到各个路口的最短距离,并借助MATLAB软件实现了算法,随后我们从中筛选出到达每个交巡警服务平台距离小30米的路口,则连接各路口之间的路线即为A区交巡警服务平台的管辖范围。 在解决警力调度方案问题时,我们建立了以最短路为目标函数的整数规划模型3,采用了求解最短路的Dijkstra算法2,并借助LINGO软件对算法进行了实现,从而得到了对进出该区的13条交通要道实现快速完全封锁的方案。 在解决增加交巡警服务平台个数和具体位置的问题时,我们把握两个原则,是各交巡警服务平台的工作量均衡,二是出警时间尽量控制在3分钟内,综合考虑两个原则,我们找出第一问结论中不包含在平台管辖范围内的路口优先考虑。从而得出需要增设的平台位置和个数。
针对问题二,在解决服务平台设置方案问题时,我们就题目给出的各区的案发率和人口密度在平台分配中所占权重,在spss软件{1}中检测不合理,然后建立多元线性回归模型4spss软件对全市80个平台进行重新分配,解决了平台分配不均匀所导致的资源浪费或资源不足问题。
在解决最佳围堵方案问题时,运用Dijkstra算法计算P点到全市各个路口的最短距离,以十分钟左右为界点确定罪犯逃跑最远范围,只要最远范围附近的平台巡警围堵时间加上罪犯开始逃跑的3分钟<罪犯从P点到达最远范围的时间,即可围堵成功。
总的来说,模型的建立思路清晰、模型简单、假设合理。该模型不仅可解决交巡警服务平台的设置与调度的优化问题,也可给生活中交巡警平台的设置、调度给予参考,可使交巡警在处理警务任务时用较短时间分配最佳救援力量,并选择最优行进路径出警,具有一定的实用性。

关键字:MATLAB LINGO Floyd算法, 整数规划模型
多元线性回归模型, spss软件, Dijkstra算法,

2
如有帮助欢迎下载支持
一, 问题的重述
这是一个对交巡警平台的设置和调度问题,由于警务资源是有限的,如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。 试就某市设置交巡警服务平台的相关情况,建立数学模型分析研究下面的问题: 1)附件1中的附图1给出了该市中心城区A的交通网络和现有的20个交巡警服务平台的设置情况示意图,相关的数据信息见附件2。请为各交巡警服务平台分配管辖范围,使其在所管辖的范围内出现突发事件时,尽量能在3分钟内有交巡警(警车的时速为60km/h)到达事发地。
对于重大突发事件,需要调度全区20个交巡警服务平台的警力资源,对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口,请给出该区交巡警服务平台警力合理的调度方案。
根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况,拟在该区内再增加25个平台,请确定需要增加平台的具体个数和位置。 2)针对全市(主城六区ABCDEF)的具体情况,按照设置交巡警服务平台的原则和任务,分析研究该市现有交巡警服务平台设置方案(参见附件)的合理性。如果有明显不合理,请给出解决方案。
如果该市地点P(第32个节点)处发生了重大刑事案件,在案发3分钟后接到报警,犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯,请给出调度全市交巡警服务平台警力资源的最佳围堵方案。

二, 模型假设及符号说明
1,模型的假设:
1 假设城区所有道路畅通无阻; 2 假设相邻两个交叉路口之间的道路为直线; 3 假设所有事发现场均在城区道路上; 4 假设单一管辖区域内某一时间无并发事件; 5 假设追捕过程警车无故障发生; 6 假设交巡警平台点均位于道路交叉路口 2,符号说明
lmax交巡警从服务平台能够达到的最远距离;
x :在A区突发事件所在地的横坐标;
y :在A区突发事件所在地的纵坐标; xiA区平台的横坐标; yiA区平台的纵坐标; xjA区结点横坐标; yjA区结点纵坐标;
3
如有帮助欢迎下载支持
xij:第i个平台去封锁第j个出口; cij :第i个平台封锁第j个出口所用时间;

三, 问题的分析
问题一,1要求在平台管辖范围内发生事件后,尽量能在3分钟内有交巡警到达。即平台到管辖范围内的所有路口路程在3千米左右。我们根据题目所给数据,借助floyd算法计算出20个交巡警平台到各个路口的最短距离,3千米为判断标准,求解出20个交警服务平台管辖范围(见表二) 2要实现要道快速全封锁,就必须使13个交通要道周围的交巡警服务平台到达出口的距离最短,若到达某些出口时间差异过大,则必然会影响到整个封锁过程,因此排除A区中与出口相距甚远的7个平台。建立整数规划模型对剩余的13个平台到13个出口进行优化分配,得出最佳的封锁方案(见表四)
3)要确定增加交巡警服务平台的个数和位置,我们必须遵循两个原则:一是增加交巡警服务平台以后各个交巡警服务平台的工作量应达到相对均衡,二是增加交巡警服务平台以后各个交巡警服务平台的出警时间应尽量短,(对于那些因距离过长而出警时间过长的地区增加平台),根据(1)所得结论,未包含在平台管辖范围内的路口优先考虑,增设服务平台,使得A区所有路口在平台管辖范围内。 问题二,1按照设置交巡警服务平台的原则和任务,分析研究现有方案的合理性。根据题目所给数据,考虑到现实生活中对平台分配方案的影响因素(地区案发率,地区人口密度),建立多元线性回归模型,运用spss分析现有分配方案中对影响因素所占权重是否合理。假设只考虑以上两种因素,根据spss所得分析数据,改变两种因素对分配方案的影响权重,运用spss分析得到影响因素合理的权重分配,代入多元线性回归模型中求出合理的解决方案。 2)在P点发生重大案件,首先假设(1)嫌疑犯犯案位置固定,2)嫌疑犯逃跑路线各个路线概率均等 3)犯人逃跑速率与警车速率相等,4)搜索逃犯进行封锁应该尽量在最快时间内。运用Dijkstra算法计算P点到全市各个路口的最短距离,以十分钟左右为界点确定罪犯逃跑最远范围,只要最远范围附近的平台巡警围堵时间加上罪犯开始逃跑的3分钟小于罪犯从P点到达最远范围的时间,即可围堵成功。

四, 模型的建立及求解
1,问题一(1)模型的建立及求解:
因为各服务平台管辖范围内出现突发事件时,尽量能在3分钟内有交巡警到达事发地。警车时速60/kmh,因此交巡警从服务平台能够达到的最远距离
3 lmax60km/h*h3km 1
60根据题中所给数据图中1mm=100米,所以3千米对应图中30mm
通过两点之间距离最短的概念入手,又因为道路之间有转弯,所以在所涵
4
如有帮助欢迎下载支持
盖的范围内,可把线路分为两种情况考虑: 直线时:
折线时:xxiyyj22iji230mm2

2i1,2,3...
xxyyjxxyy;

2
j
j建立floyd算法模型求解A区各点之间的最短距离(程序见附录一) 从图的带权邻接矩阵A=[a(i,j]n×n开始,递归地进行n次更新,即由矩D(0=A,按一个公式,构造出矩阵D(1;又用同样地公式由D(1构造出D(2;……;最后又用同样的公式由D(n-1构造出矩阵D(n。矩阵D(nij列元素便是i号顶点到j号顶点的最短路径长度,称D(n为图的距离矩阵。
根据题目所给A区的交通网络与平台设置的示意图,找出到平台路程为3米左右的所有路口,即为此平台的管辖范围,同理求解出20个交警服务平台管辖范围。

基于以上分析,借助floyd算法,用matlab软件求得A区各个路口间的最短距离见下表:
表一 A区各个路口间的最短距离 路口号 1 2 3 4 5
0 0 0 0 0 1
18.86796 0 0 0 0 2 30.5655 21.07724 0 0 0 3 36.96282 40.9176 26.61766 0 0 4 75.9276 72.00694 51.04165 42.02678 0 5


接下来,利用上表计算所得数据,以路程30mm为判断标准,得到A区各交巡警服务平台的管辖范围,其结果如下表所示:

表二 A区各交巡警服务平台管辖范围
服务平台号
1 2 3 4 5 6 7 8 9 10 11 12 13
路程相距3千米的路口
75,76,77,19,79,68,69,70,43,2,44,74,73,72,71,78 44,43,70,3,55,67,68,66,65,69,1,74,73,72,71
65,66,6768,44,2,70,55,54 60,62,63,64,65,66,67,58,57 47,48,49,50,30,51,7,52,53,56 47,48,57,58,59,30,75
30,48,47,5,6,37
31,32,33,34,35,36,379,45,46,47, 33,34,35,36,3216,46,45,47,318,37
10 25,26,27 25 21,22,23,24
5
如有帮助欢迎下载支持
14 15 16 17 18 19 20 14 31
35,36,37,34,45,46,89
40,41,42,43,71,726970,43,2,44,74,73
81,82,83,90,8920,91,8079,19,77,76,78,1,74,73,72 77,76,64,75,79,80,18,83,82,81,84,78,1,74,73,71,69,70,68
85,86,87,88,89,90,82,83,18,81,84
(其中2829,38,39,61,92六个路口未能包括在管辖范围内) 2,问题一(2)模型的建立及求解:
要实现要道快速全封锁,就必须使13个交通要道周围的交巡警服务平台到达出口的距离最短,若到达某些出口时间差异过大,则必然会影响到整个封锁过程。
运用matlab做出A区各个路口的散点图,并标示出封锁口和平台位置如下:
图一 A区个路口散点连线图
(红色点为封锁出口点,蓝色小圈代表交巡警服务平台)
根据以上分析,首先排除A区中与出口相距甚远的7个平台:
1,2,3,17,18,19,20
剩余的13个平台与13个封锁口之间最短距离如下表:
表三 13个平台与封锁口之间的最短距离
12
27.4
14 22.5
79 22.6
82 7
20.14 15.36 10.95 10.98
8.216 8.66 6.202 6.232
4.1421 18.56 14.24 14.24
11.422 20.36 16.04 16.04
13.223 21.76 17.44 17.44
14.624 23.64 19.32 19.29
16.5
6 28 16.19 11.29 11.32 8.56 29 15.49 10.59 10.62 8.01 30 38 11.712.04 6 9.603.16
2 9.373.19
2
0.58 7.5448 7.36 2.46 2.46
1.262 0.35 6.07 5.32 4.88

如有帮助欢迎下载支持
8 9 10 11 12 13 14 15 16
06 21.26 20.112 7.57 3.78
2 9 8 8 8 6 9 8 9.42.6812.614.415.817.710.211.13.786.083.06.68
22 2 82 82 88 62 16 56 7 2 8 8.21.5311.513.314.716.64.934.934.29.76 10.7 7.83
72 2 32 32 32 12 7 2 3
9.2414.115.19.49 7.69 9.09 8.23 4 8 2 13.621.922.89.5 7.7 9.1 3.8 02 2 2 16.720.721.88.58 6.78 6.38 3.5 68 4 4
15.116.012.7 2.7 0.9 0.5 2.38
2 6
10.06.74 3.26 5.06 6.46 8.34 9.16
6
7.9510.812.77.67 9.47 4.75 5.69 8 7 5
15.011.212.20 10 11.8 13.2
8 92 32
10.277 14.707 17.877 14.76 12.572 17.002 20.172 16.06 10.18.8
4 11.34.39
58
4.723.4 8
9.513.17 7 14 17.6 17.17 15.46 9.507 5.097 5.435
20.77 19.06 13.107 8.697 3.6
12.75 8.32 11.0
78 5.85.98 6
11.0 84 16.4.425 1 18.6.758 4
接下来,用剩余的13个平台与13个封锁口建立整数规划模型如下: 我们记派第i个平台封锁第j个出口记为xij
0xij 1i个平台没封锁j出口
i个平台封锁了j出口cij为第i个平台去封锁第j个出口所需的时间。 因为(1)每个平台只能封锁一个出口。 我们得到模型为
minzcijxiji1j11313s.txj113i113ij1(i1,2,3,4...1(j1,2,3,4...
x
ijxij0or1运用lingo软件编程,计算得到A区交巡警服务平台的合理调度方案如下表所示(程序见附录二)

表四 A区交巡警服务平台的合理调度方案 封锁口 交巡警服务平台 封锁口 交巡警服务平台
7
如有帮助欢迎下载支持
12 12 28 15 14 8 29 5 16 16 30 7 21 14 38 9 22 10 48 6 23 13 62 4 24 11 3,问题一(3)求解: 根据(1)中所得结论:
28,29,38,39,6192 六个路口未能包括在交巡警服务平台的管辖范围内,因此我们将增添的5个平台增设在29,38,39,61,92五个路口处,使得A区所有路口都处于交巡警服务平台的管辖范围内。

4,问题二(1)模型的建立及求解:
首先,分析研究现有方案的合理性。根据题目所给数据,考虑到现实生活中对平台分配方案的影响因素:地区案发率,地区人口密度,结合以上两种影响因素与平台分配之间的联系,我们便可以分析题目中交巡警服务平台的合理性: 多元线性回归模型的建立:
若因变量Y与解释变量X1X2X3X4……具有线性关系,它们之间的线性回归模型可表示为:
b1bX+2b2X+k Y=01k+b X+其中为随机扰动项观测值。对于第i个观测值:
Y=b0b1X1i+b2X2i++bkXki+i1,2n
即:
Y11Y12Yn1也即:YXb
假定:
x11x12x1nx21x22x2nbxk101b1xk2b2 2xknnbk8
如有帮助欢迎下载支持
Ei0VariEi2u2Covi,j0Covxj,i0i1,2,n,ni1,2n
ij,i,j=1,2j1,2,k~N0,u2Jn建立案发率,人口密度与平台分布之间的关系
通过题目所给数据分析出,案发率,人口密度在平台分布的设置中所占权重是不同的,各个权重定义为:b1, b2,我们定义区域平台分布数Y的函数为:
Y=b0b1X1+b2X2+
根据题目所给数据建立平台分布影响因素表如下:

表五 平台分布影响因素表
A B C D E F
人口密度所面积 人口 (/权重 占权重
方公里
22 60 27272.73 0.734553103 124.5 0.18 15 59 103 21 2038.83 0.054913054 66.4 0.1 8 4 221 49 2217.19 0.059716933 187.2 0.28 22 5 383 73 1906.01 0.051335732 67.8 0.1 8 4 432 76 1759.26 0.047383225 119.4 0.18 14 4 274 53 1934.31 0.052097953 109.2 0.16 13 4

利用spss软件对题目现有分配进行检测: 题目中给出对全市交巡警服务平台分布为:

表六 现有全市交巡警服务平台分布 区域 A B C D E 平台数 20 8 17 9 15
spss分析数据如下图:

9 F 11
如有帮助欢迎下载支持
2 spss对现有平台分布的检验

根据上图分析数据可得:Sig=0.11>0.05则说明现有平台分布不合理; 极端考虑只有发案率,人口密度两个因素影响分布结果,假设两个因素权重比为1:1,则得到假设分布表如下:

表七 假设交巡警平台分布 区域 A B C D E F 平台数 37 6 13 6 9 7
spss分析数据如下图:

图三 spss分析假设平台分布的合理性
10
如有帮助欢迎下载支持

根据上图分析数据可得:Sig=0.00<0.05则说明假设平台分布合理; 即发案率所占权重为b1=0.465,人口密度所占权重为b2=0.509,常量b0=0.015 因此给出的解决方案如下:将A区平台数增加至37个,B区平台数减少至6个,C区平台数增加至13个,D区平台数减少至6个,E区平台数减少至9个,F区平台数减少至7个。 重新分配各区平台分布如下表:
表八 重新分配各区平台分布表 区域 A B C D E F 平台数 37 6 13 6 9 7
5,问题二(2)模型的建立及求解:
针对P(第32个节点)发生特大刑警事件,对事件进行背景分析: 1,在三分钟后才接到报案; 11
如有帮助欢迎下载支持
2,嫌疑犯逃跑路线各个路线概率均等 3,嫌疑犯逃跑时是驾车逃跑,假设其速率与警车速率相等 4,搜索逃犯进行封锁应该尽量在最快时间内
根据题目所给数据我们经过MatlabDijkstra算法求出P点到图中所有点的最短距离如下表(程序见附录三)
表六 P点到全市各个路口点的最短距离 全市点个数 1 2 3 4 5 P点到此点87.09191 77.92946 57.63896 59.42432 24.69818
的距离
经过MatlabDijkstra算法求出P点到图中所有点的最短时间。 根据最短时间,我们将要在尽量最短的时间内抓到嫌疑犯,转化成在最短时间内控制一个圈的所有结点,让嫌疑犯不能逃出一定区域,我们将这个在一定时间内的时间设在十分钟左右,我们找出在此范围内的结点(除去嫌疑犯3分钟逃出区域内的结点),因为数据量太大,在此我们选出部分结点进行说明,现举出部分可能有用的结点: 表七 部分结点标号
结点标号 28 29 38 40 41 214 215 241 248 370 371 561 将上面信息进行整理得出P点到达各个割点的时间 和由前面程序求出的表格中找出最近平台到达结点的时间 ,当 时说明控制结点成功。最后得出围堵方案如下:
1446224017414805615536501529320371167248174214171241169253
说明:箭头左边为交巡警服务平台所在的路口标号,右边是该平台需要围堵的路口标号。 从而得出,完成合围所花费的时间为7.36分钟。
综合评价
1. 模型的优点
我们所建立的模型思路清晰、假设合理,适用于解决实际问题,通用性、推广性强。
1问题一中首先运用Floyd算法计算出A区交通示意图所有节点之间的最短距离与最短路径,该方法极大地简化了传统算法的计算过程,提高了计算的准确性,为后续问题的求解打下了良好的基础。
针对A区交巡警服务平台的管辖区域的划分,在已得到的各节点之间的最短距离与最短路径的基础上,建立整数规划模型,运用lingo软件求解,快速给出调度方案。并将交巡警服务平台管辖的路口节点以表格的形式给出,清晰明了、简单直观,便于查找与调度。
2问题2中,对于该市现有的交巡警服务平台设置方案的合理性评价,我们兼顾了案发率与人口密度两大影响因素,对所给数据的运用较为全面,因而具有很强的合理性与可行性。 2. 模型的缺点
1由于本题所涉及的数据量较大,部分数据的四舍五入可能给结果带来一定的误差。
2假设量过多,导致模型存在一定的程度的理想化,如现实中无法保12
如有帮助欢迎下载支持
证道路的畅通性;将巡警车与逃犯速度均设为定值且相等,忽略了实际生活所存在的不确定因素,与真实情形尚有一定差距,将或多或少地对实际结果产生影响。
3. 模型的改进
针对于围堵问题,为了使模型具有更强的实用性与合理性,我们对模型做出如下改进:我们假设的是巡警车与罪犯速度相同,但实际中往往是不相等。因此我们将分情况讨论来建立模型。 4. 模型的推广
我们所建立的模型较好地划分了交巡警平台的管辖范围,解决了围堵罪犯的相关问题,有效地改善交巡警平台的工作效率,平衡了各个服务平台的工作量。在经济和科技飞速发展的今天,该模型还可广泛地运用到其他诸如消防救援、安全事故应急等工作中去,具有极高的借鉴意义。

参考文献
[1] blustin spss线性回归分析 815
[2] XIANG__jiangsu Dijkstra算法详解

815
[3] 建模资料 5 运筹学模型 第四,五页 [4] 灰白先知 多元线性回归模型 816



附录一:
floyd算法模型求解A区各点之间的最短距离:
X = A(:,1;
13
如有帮助欢迎下载支持
Y = A(:,2; N = length(X; D = zeros(N,N; for I = 2:N for J = 1:I-1 D(I,J = sqrt((X(I - X(J*(X(I - X(J + (Y(I - Y(J*(Y(I - Y(J; end end D; D1 = D+D'; plot(D; 附录二:
A区交巡警服务平台的合理调度方案:
model: sets: ck/1..13/:; fs/1..13/:; link(ck,fs:x,c; endsets data: c=输入矩阵。 enddata
min=@sum(link(i,j:c*x; @for(ck(i:@sum(fs(j:x(i,j=1; @for(fs(j:@sum(ck(i:x(i,j=1; @for(link:@bin(x; end

附录三:
Dijkstra算法求出P点到图中所有点的最短距离: m=size(B,1; A=zeros(m,m+inf; for i=1:m A(i,i=0; end [m n]=size(C; for i=1:m A(C(i,1,C(i,2=sqrt((B(C(i,1,2-B(C(i,2,2^2+(B(C(i,1,3-B(C(i,2,3^2; A(C(i,2,C(i,1=A(C(i,1,C(i,2; end a=dijkf2(A/10; 14
如有帮助欢迎下载支持
time=a(32,:;

15

本文来源:https://www.2haoxitong.net/k/doc/d4a9615c710abb68a98271fe910ef12d2bf9a91d.html

《交巡警平台分配问题.doc》
将本文的Word文档下载到电脑,方便收藏和打印
推荐度:
点击下载文档

文档为doc格式