高教社杯全国大学生数学建模竞赛.docx
- 文档编号:25781498
- 上传时间:2023-06-14
- 格式:DOCX
- 页数:31
- 大小:97.77KB
高教社杯全国大学生数学建模竞赛.docx
《高教社杯全国大学生数学建模竞赛.docx》由会员分享,可在线阅读,更多相关《高教社杯全国大学生数学建模竞赛.docx(31页珍藏版)》请在冰豆网上搜索。
高教社杯全国大学生数学建模竞赛
2011高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们将受到严肃处理。
我们参赛选择的题号是(从A/B/C/D中选择一项填写):
B
我们的参赛报名号为(如果赛区设置报名号的话):
所属学校(请填写完整的全名):
山西大学
参赛队员(打印并签名):
1.田彩星
2.刘艳红
3.程卓
指导教师或指导教师组负责人(打印并签名):
日期:
2011年9月12日
赛区评阅编号(由赛区组委会评阅前进行编号):
2011高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
赛区评阅记录(可供赛区评阅时使用):
评
阅
人
评
分
备
注
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
交巡警服务平台的设置与调度
摘要
本文研究的是城市交巡警的分配及围堵方案的制定问题,结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,在现有警力的情况下,尽可能科学地合理地安排各个交巡警平台的管辖范围,从而达到效率最大化目的。
文章针对分配各交巡警服务平台的管辖范围的问题,建立0-1整数规划模型,采用“就近原则”给出求解模型的多项式时间算法。
针对考虑工作量均衡要求的平台辖区分配问题和调度警务资源问题,提出相应的整数规划模型,并将其转化为图论问题,通过二分搜索,多次调用最大匹配算法求解模型,得到最优解。
针对合理设置交巡警服务平台问题,综合考虑时间要求和工作量均衡要求,给出度量平台设置合理程度的指标,并通过穷举法给出A区增加平台的方案。
根据提出的指标,分析研究该市现有交巡警服务平台设置方案。
对该市发生的重大刑事案件,也给出一个围堵嫌犯的方案。
关键词:
整数规划二部图最大匹配交巡警服务平台最短距离
1
一、问题的背景
面对各种突发事件,即使在科技高度发达的今天,也有显得束手无策的时候,近十年来,我国科技带动生产力不断发展,国家经济实力不断增强,然而另一方安全生产形势却相当严峻,每年因各类生产事故造成大量的人员伤亡、经济损失。
尤其是在一些大目标点,作为人类经济、文化、政治、科技信息的中心,由于其“人口集中、建筑集中、生产集中、财富集中”的特点,一旦发生重大事故,将会引起相当惨重的损失。
为了保障安全生产、预防各类事故。
我国正在各省(市)目标点逐步设立交巡警平台。
这一警种拥有包括枪支在内的“高精尖”装备,代替过去的交警和巡警。
交巡警平台是将刑事执法、治安管理、交通管理、服务群众四大职能有机融合的新型防控体系。
在人流量极大、治安状况比较复杂、交通持续比较混乱的事故多发带产生强大的司法制衡力、社会治安的驾驭力、打击罪犯的冲击力。
保证在事故发生的第一时间赶到现场。
大力的减少了社会上各种混乱行为的发生。
使居民的生命财产安全得以保障。
二、问题的重述
警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。
为了更有效地贯彻实施这些职能,需要在市区的一些交通要道和重要部位设置交巡警服务平台。
每个交巡警服务平台的职能和警力配备基本相同。
由于警务资源是有限的,如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。
现我们就某一市路口节点和交巡警服务平台已知的条件下,根据交巡警服务平台设置原则和任务,合理的分配交巡警服务平台,根据实际情况建立数学模型分析研究下面的问题:
问题一:
附录中给出了该市中心城区A的交通网络和现有的20个交巡警服务平台的设置情况示意图,相关的数据信息见附录。
请为各交巡警服务平台分配管辖范围,使其在所管辖的范围内出现突发事件时,尽量能在3分钟内有交巡警(警车的时速为60km/h)到达事发地。
问题二:
对于重大突发事件,需要调度全区20个交巡警服务平台的警力资源,对进出该区的13条交通要道实现快速全封锁。
实际中一个平台的警力最多封锁一个路口,请给出该区交巡警服务平台警力合理的调度方案。
问题三:
根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况,拟在该区内再增加2至5个平台,请确定需要增加平台的具体个数和位置。
问题四:
针对全市(主城六区A,B,C,D,E,F)的具体情况,按照设置交巡警服务平台的原则和任务,分析研究该市现有交巡警服务平台设置方案(参见附录)的合理性。
如果有明显不合理,请给出解决方案。
问题五:
如果该市地点P(第32个节点)处发生了重大刑事案件,在案发3分钟后接到报警,犯罪嫌疑人已驾车逃跑。
为了快速搜捕嫌疑犯,请给出调度全市交巡警服务平台警力源的最佳围堵方案。
2
三、模型的假设
1.假设题目所给数据真实可靠;
2.交巡警到达出事地块路口即为到达出事地点;
3.交巡警到达每条街道的路口即为到达此街道;
4.在任何车况、路况下,在任何时刻警车的时速为60km/h,且是匀速行驶;
5.罪犯逃逸的时速与警车的时速相等。
四、模型的建立与求解
准备工作
首先构造边赋权无向图G如下:
以A区所有路口的集合作为图G的顶点集合,两个顶点相邻当且仅当对应的两个路口是同一条街道的两端;每条边上的权为对应街道的长度,即其两端路口
和
的欧氏距离
(程序见附件)。
然后在图G上运用Floyd算法可求出A区92个路口两两之间的最短距离矩阵D和最短路径矩阵R(见附件)。
问题一
1.模型I的建立
问题一要求为各交巡警服务平台分配管辖范围,使其在所管辖的范围内出现突发事件时,尽量能在3分钟内有交巡警(警车的时速为60km/h)到达事发地。
对A区的每个路口
和每个交巡警服务平台
,用0-1变量
表示平台
对路口
的管辖关系。
当路口
由平台
管辖时,
;否则,
。
要求每个路口都由一个平台管辖,即有
目标是在尽量短时间内有交巡警到达事发地,因此,对应的目标函数为
其中,
表示由平台
到达路口
的最短时间。
由假设4,警车的时速恒为60km/h,因此时间与距离成正比(在此,时间(分钟)
10=距离(千米)。
)这说明上述的目标函数等价于
3
其中,
表示由平台
到达路口
的距离。
综上,我们得到问题一的0-1整数规划模型P1:
2.模型I的求解
模型I虽然是个整数规划,但存在一个好的算法去求解它。
在此,我们利用“就近原则”为各路口分配所属交巡警服务平台。
也就是根据A区92个路口两两之间的最短距离矩阵D,对每个路口
,计算它与所有20个交巡警服务平台的距离。
如果在这20个交巡警服务平台,路口
与平台
的距离最短,则路口
就由平台
管辖(即,此时有
)。
由此“就近原则”所得的交巡警服务平台管辖路口的方案为:
Aorder=[12345678910
11121314151617181920
1313131312111115157
9789916161622
171722987755
5553354554
7444331112
1211111911918
18181820202020202020
2020],
这里Aorder是一个92维向量,它的第
个元素Aorder(
)表示第
个路口由第Aorder(
)个平台管辖。
比如:
Aorder(21)=13表示第21个路口由第13个平台管辖。
(见附件)
3解法的正确性和复杂度
3.1解法的正确性
设由上述“就近原则”所得的问题P1的解为
。
假设
为P1的任意一个解。
设解
的值为
。
在解
中,必有管辖
路口的平台,设为
。
因此有
。
由“就近原则”,离
路口最近的平台为
。
因此,
,解法的正确性得证。
3.2解法的复杂度
4
由该解法解有
个路口、
个平台的问题,只需检查每个路口到每个平台的距离矩阵找出相应的最小值,其复杂度为
。
3.3解的分析
在上面由“就近原则”所给的分配方案Aorder中,从15号平台到29号路口的时间最长,为5.70053分钟。
大于3分钟的巡查路线有从15好平台到28号路口,从16号平台到38号路口,从7号平台到61路口,从20号平台到92号路口。
由2.3.1中的讨论知,在当前平台的设置方案下,任何平台管辖范围的分配方案都不能保证当出现突发事件时,在3分钟内有交巡警到达事发地。
4.问题一的模型II
模型I仅仅考虑了时间的因素,而未考虑到交巡警服务平台工作量的均匀程度,因此在模型I的基础上我们建立了模型II。
该模型同时考虑到时间和工作量共同作用的影响。
可将
路口的发案率
看为交巡警服务平台管辖该路口时所需承担的工作量。
因此,
是平均工作量。
为使交巡警服务平台的工作量均衡,可增加约束条件:
其中,
,
为给定的两个较小的整数。
由此得模型P2:
5.问题一的模型III
5.1模型III的建立
如果在模型P2中不考虑发案率,认为每个路口的工作量都是相等的,并进一步取
,就得到下面的模型P3:
5
5.2模型III的求解
现要将A区的92个路口分配给20个交巡警服务平台。
对给定的正数K,构造二部图G(K)如下:
将每个交巡警服务平台看成5个顶点,共有100个顶点组成一部;将每个路口看成一个顶点,共92个顶点组成一部。
两个不同部中两点连边当且仅当对应的平台和路口的距离不大于K。
利用匈牙利算法求出G(K)的最大匹配M。
对属于M的每条边ij,令
。
如果M的边数=92,则说明
满足P3的约束条件,从而P3的最优值不大于K;如果M的边数<92,则说明
不满足P3的约束条件,从而P3的最优值大于K。
由此思想,我们通过二分搜索确定临界的K值,也就是P3的最优值(见附件)。
5.3模型III的解
由上述算法求得的P3的解(即,交巡警服务平台管辖路口的方案)为:
Aorder=[2455816151010111212141415163221141414131312131515151010101616779943424887168696995867667553332171717171819201181920118192011841920118192017],
这里Aorder是一个92维向量,它的第
个元素Aorder(
)表示第
个路口由第Aorder(
)个平台管辖。
比如:
Aorder(21)=14表示第21个路口由第14个平台管辖。
5.4模型III的解的分析
在由“就近原则”获得的分配方案中(见2.2节),一个平台到由它管辖的路口的最长时间为5.70053分钟。
在模型III中的分配方案中,一个平台到由它管辖的路口的最长时间为6.50592分钟。
这说明在时间指标上“就近原则”获得的分配方案优于模型III中的方案。
对一个给定的分配管辖范围的方案中,用
表示交巡警服务平台工作量均衡程度,其中
表示
平台的工作量,
为平均工作量(详见问题一的模型II中所述)。
在由“就近原则”获得的分配方案中,
。
在模型III中的分配方案中,
。
这说明在考虑工作量的均衡时,模型III中的分配方案优于由“就近原则”获得的分配方案。
6
问题二:
问题二需调度全区的20个交巡警服务平台对进出该区的13条交通要道实现快速全封锁,且要求一个平台的警力最多封锁一个路口。
同问题一的模型III类似,该问题的0-1规划模型P4为:
与问题一的模型III的解法类似,对正数K,构造一部为20个平台、另一部为13条要道的二部图G(K)。
如果M的边数=13,则P4的最优值不大于K;如果M的边数<13,则P4的最优值大于K。
由此思想,我们通过二分搜索,多次调用求二部图最大匹配的匈牙利算法,确定临界的K值,也就是P4的最优值。
求解程序见附录。
其解如下:
路口节点
12
14
16
21
22
23
24
28
29
30
38
48
62
分管交巡警平台
10
16
2
11
13
14
12
15
7
5
1
4
18
所需时间(分钟)
7.5866
6.7417
7.3881
5.0723
0.9055
6.4733
3.5916
4.7518
8.0155
3.1829
5.8809
7.3959
6.7344
由上表我们可以得到20个平台将13个交通要道进行封锁,需要的最长时间是从7号平台到29号路口,这需要8.0155分钟。
7
问题三:
为了解决问题三,需要在该区增加2—5个平台,而对于每个平台的设置方案,根据其所管辖的路口,总可以得到其所对应的管辖区域,即对应于每个平台都有一个分片方案。
在整体设置中,设每个平台所管辖的路口发案率为
进而可以算出工作量
,最后得到工作量的标准差
,其中
是每个平台工作量的期望;而根据第一问的求解过程,我们可以得到每个平台到各个路口的最长时间
。
设
是比例系数
,
是每个方案的得分。
通过穷举法,我们得到在不同比例系数和增加不同平台个数下得分最低的方案的得分情况。
下表给出了相应的数据。
比如下表第4行第4列元素25.8930表示比例系数λ=0.3时,在所有增加2个平台的方案中分数
最低的方案对应的分数为25.8930。
比例
G
台数
λ=0.1
λ=0.2
λ=0.3
λ=0.4
λ=0.5
λ=0.6
λ=0.7
λ=0.8
λ=0.9
0
17.0917
21.5266
25.9614
30.3962
34.8311
39.2659
43.7007
48.1356
52.5704
1
17.7747
22.1337
26.4926
30.8516
35.2105
39.5695
43.9284
48.2874
52.6463
2
17.0038
21.4484
25.8930
30.3376
34.7822
39.2268
43.6714
48.1160
52.5606
3
17.2783
21.6924
26.1065
30.5206
34.9347
39.3488
43.7629
48.1770
52.5911
4
17.9076
22.2518
26.5960
30.9401
35.2843
39.6285
43.9727
48.3169
52.6611
5
18.6242
22.8888
27.1533
31.4179
35.6825
39.9470
44.2116
48.4761
52.7407
8
通过对数据的分析,我们可以得最优到当增加2个平台时,交巡警服务平台的工作量不均衡和有些地方出警时间过长的情况可以得到最优缓解。
当λ=0.2时,最优增加平台的方案为在28,30两个路口处增加平台。
问题四
根据上题A区中分配管辖范围的“就近原则”的方法,用可分别得到B,C,D,E,F区的距离矩阵以及各区交巡警服务平台的分配管辖范围,为了使交巡警服务平台设置的更加合理,我们设置了标准差,平均工作量,最长时间,平均时间以及利用权重对各区评定分数,将以上五个方面作为评判分配合理性的标准,拟以下表格加以对A,B,C,D,E,F六区进行比较:
区
工作量标准差
最长时间
总时间
得分数
(权重:
W)
W=0.1
W=0.3
W=0.5
W=0.7
W=0.9
A
12.6569
57.0053
1.0322e+003
17.0917
25.9614
34.8311
43.7007
52.5704
B
10.4048
44.7031
1.1146e+003
13.8346
20.6943
27.5540
34.4136
41.2733
C
22.6768
68.6054
3.5962e+003
27.2697
36.4554
45.6411
54.8269
64.0126
D
10.1054
160.6282
1.3030e+003
25.1577
55.2623
85.3668
115.4713
145.5759
E
15.4957
191.0506
2.6137e+003
33.0512
68.1621
103.2731
138.3841
173.4951
F
16.3426
84.7985
2.6249e+003
23.1882
30.0338
43.7250
57.4161
77.9529
表一
由表格中的数据分析可得:
该市E区的工作量的标准差最大,得分最高,说明其工作量的分配不均衡,平台设置最不合理。
我们分析造成这种情况的原因可能是由于E区人口多,城区面积大。
问题五:
若该市地点P(第32个节点)处发生了重大刑事案件,且在案发后犯罪嫌疑人立即驾车逃跑,三分钟后接到报警。
设嫌疑人逃跑时的时速与警车的时速相同,均为60Km/h,每个交巡警服务平台可以同时派出两个警队。
为了快速搜捕,调度全市警察进行围堵,由假设我们建立如下解决方案:
11、12、13、16号交巡警服务平台原地围堵,4号交巡警服务平台围堵标号为62的出入口,F区的480号交巡警服务平台围堵第561号路口节点,15号交巡警服务平台围堵第29个出入口,D区的351号交巡警服务平台围堵第371号路口节点,C区的173号交巡警服务平台围堵第235、237号路口节点。
9
由问题一我们已求得20个交巡警服务平台对A区的十三个交通要道实现快速封锁需要8分钟,在案发3分钟后接到报警,此时嫌疑人已逃跑11分钟。
首先我们找出P点到各交通要道的时间,其中大于11分钟的出入路口有12、21、22、23、24五个出入口,那么12、11、13三个路口节点的交巡警服务平台可原地围堵;大于3分钟的出入路口有14、16、28、29、38、62六个出入口。
从P到62号出入口所用时间为9.1分钟,而离62号出入口最近的4号交巡警服务平台到62号出入口所需时间为0.35分钟,所以在接到报警后由4号交巡警服务平台围堵第62号出入口;从P到16号出入口所用时间为3.分钟,所以接到报警后16、14号交巡警服务平台原地围堵,那么嫌疑人不会从12、14、16、21、22、23、24号出入口逃出;在对16号出入口进行围堵之后,嫌疑人只能从经由39号路口节点到达38号出入口,此时所需时间为6.49分钟,离其最近的17号交巡警服务平台到38号出入口需4.76分钟,所以嫌疑人进入F区,经计算F区的480号交巡警服务平台到达561号路口节点所需时间为4.5分钟,嫌疑人进入F区需8.8分钟,所以480号交巡警服务平台围堵第561号路口节点;从15号交巡警服务平台到第29号出入口所需时间5.7分钟,从P到第29号出入口所需时间9.2分钟,所以15号交巡警服务平台围堵第29号路口节点;对于28号交巡警服务平台,此时嫌疑人已进入D区离其最近的第351号交巡警服务平台围堵第371号路口节点;小于3分钟的出入口有30、48两个,那么在接到报警之前嫌疑人已经进入C区,经计算由C区的173号交巡警服务平台围堵第235、237号路口节点。
以上方案我们假设每个交巡警服务平台可以同时派出两个警队,按照设置交巡警服务平台的原则和任务,每个警局只能围堵一个出入口,所以对进入C区的犯罪嫌疑人我们可以将所有的出入C区的路口封锁。
若将该问题进行推广,在任何一处案件发生,罪犯逃离时,在接到报案的同时,各平台警员都可在第一时间赶赴交通要道(该市各区间的出入口)进行围堵。
七、模型的评价与推广
1.模型的评价
(1)上述所有模型的建立都是在假设的前提下,没有考虑城区人口,城区的面积,车速的变化以及道路交通状况等因素造成的影响。
(2)在上述模型中我们均以交巡警平台到各个路口的时间最短作为目标函数,将工作量尽可能上均衡这一要求放在了约束条件中,事实我们也可以将工作量的标准差最小作为目标函数,将交巡警平台到各个路口的时间最短作为约束条件。
(3)对于0-1整数规划我们没有用传统的求解规划问题的方法来求解,而是转化为了图论中的二部图最大匹配问题来求解。
由于求解二部图最大匹配有好的算法,因此此法可推广用于求解更大规模的问题。
遗憾的是在将模型P2转化为P3时丢失了其中的发案率因素,而仅仅考虑了路口的数目这一因素。
2.模型的推广
本模型较好的解决了交巡警平台的最优选址问题,当案件发生时,交巡警可以第一时间到达事发地点,有效的改善了交巡警在执行任务中的效率,在经济迅猛发展的今天,城市加速扩张,人口迅速增长,交巡警平台的设置是平安城市的最好保障.该模型也可运用到其他最优选址问题中去,比如关于消防救援工作最优路径问题、重大生产安全事故应急救援问题、公共交通的最优路径问题等.同时也可利用该模型算法拓展模型在其他领域的适用范围.
10
参考文献
[1]:
钱湔.运筹学[M].北京:
科学出版社,2000.
[2]:
薛定宇,陈阳泉.初等运用数学效果的matlab求解[M].北京:
清华大学出版社,2004.
[3]:
石辛民,郝正清.基于matlab的适用数值计算[M].北京:
清华大学出版,北京交通大学出版社,2006.
11
附录:
问题一的程序
function[Aorder]=Asortorder(A)
%用“就近原则”求辖区分配方案。
%A是警局到路口的距离矩阵,A(i,j)表示从i警局到j路口的距离。
%Aorder是分配方案向量,Aorder(i)=j表示第i个路口归第j个警局管。
m=size(A,2);
Aorder=zeros(1,m);
n=size(A,1);
fori=1:
m
Aorder(i)=1;
forj=1:
n
ifA(j,i) Aorder(i)=j; end end end end floyd算法: function[D,R]=floyd(A) %floyd算法,由图的邻接矩阵
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高教 全国大学生 数学 建模 竞赛