高教社杯全国大学生数学建模竞赛Word格式.docx
《高教社杯全国大学生数学建模竞赛Word格式.docx》由会员分享，可在线阅读，更多相关《高教社杯全国大学生数学建模竞赛Word格式.docx（24页珍藏版）》请在冰豆网上搜索。
这既可以看做线性规划中的指派问题，又可以看做是偶图中求最优匹配问题。
为此，先虚拟7个交通要道节点，使得每个平台到这些虚拟节点的距离为０，构造出赋权平衡偶图，建立整数规划模型，通过lingo软件编写程序求得调度方案。
针对问题一的第三部分，要求增设2~5个服务平台解决工作量不均衡和有些地方出警时间过长的实际问题。
通过对工作量的分析建立基本的工作量函数，求得工作量向量，并根据一定的调整原则对原有方案进行优化。
问题二分为两部分，其中第一部分是把问题一中A区的平台设置方案扩展到全市，判定其是否合理。
具体来讲，是将全市的交通网络数据应用到问题一所建立的数学模型，求出服务平台设置方案，再依据工作量和出警时间，对求出的方案进行分析和评价，并给出不合理部分的解决办法。
针对问题二的第二部分，要求警方在最短时间内围捕犯罪嫌疑人，先判断犯罪嫌疑人最可能出该区的路口（即以最短距离出A区的几个路口），从而分配警力封锁之，同时，将剩余警力分配到A区其它的出城路口，得到封锁A区的最佳围堵方案。
如果犯罪嫌疑人逃出A区，利用问题一第二部分所给模型，求出犯罪嫌疑人出市区各路口的最佳围堵方案。
关键字：
最短路线距离矩阵；
推广的点覆盖；
指派问题；
0~1规划；
最优匹配问题；
工作量函数
一、问题重述
“有困难找警察”，是家喻户晓的一句流行语。
警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。
为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。
每个交巡警服务平台的职能和警力配备基本相同。
由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。
试就某市设置交巡警服务平台的相关情况，建立数学模型分析研究下面的问题：
（1）附件1中的附图1给出了该市中心城区A的交通网络和现有的20个交巡警服务平台的设置情况示意图，相关的数据信息见附件2。
请为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。
对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。
实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。
根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。
（2）针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案（参见附件）的合理性。
如果有明显不合理，请给出解决方案。
如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。
为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。
二、问题假设及符号说明
问题假设：
1、交警从服务平台到事发地点所用时间只与路程有关，即不记上车下车时间、不记中间停留时间，警车到事发地之间道路通畅，无特殊路况，能由最短路径到达事发地点；
2、每个服务平台一般管理整条街道，即不存在某条街道分成两个平台管理；
3、只在节点处增加平台；
4、交警在接到报警后立即行动，反应时间为0；
5、所有案件处理时间为平均处理时间。
符号说明：
：
A区交巡警服务平台i
B区交巡警服务平台i
C区交巡警服务平台i
D区交巡警服务平台i
E区交巡警服务平台i
交叉路口节点
警车速度：
v=60km/h
:
交警处理案件时间
，j>
20，i平台到节点j所需最短距离
（
）全市最短路线距离矩阵
A（
）全市邻接矩阵，其中当
为有限数时，表示节点i到节点j的街道长度，当
为无穷大时，表示节点i与j直接没有街道相连。
i节点的发案率
）：
A区任意两点最短路线距离矩阵；
）A区邻接矩阵；
（i=1、2…20;
j=1、2…20）表示20个平台到20个节点的矩阵，其中前13个为各对应交通要到节点，后7个为虚拟节点;
节点
的坐标
平台工作量
三、问题的分析
本题是研究如何合理的设置交巡警服务平台，分配各平台的管辖范围、调度警卫资源。
以实现交巡警服务平台更快更好的服务群众的目标。
本题给出了某市区的所有交通网络及各节点在比例图中的坐标，并说明了各相关联的节点及各节点的发案率，指出了现有交巡警平台所在节点位置及进出口市区交通要道所在节点位置。
并给出了市区内各区面积和人口数量。
题目设有两个问题，
（1）问题分为三个部分，一部分要求为现有交巡警平台分配管辖范围，使得在其管辖范围内出现突发事件时，尽量能在三分钟之内有交警到达事发地点;
二部分要求对于突发事件，需要调度20个交巡警服务平台的警力资源到13个交通要道的最佳调度方案；
三部分要求对于工作量不均衡和有些出警时间过长的实际情况增加2到5个服务平台使情况的得到改善。
根据题意，需求交警在三分钟之内所到达的范围，初步分析应建立以平台为中心，三分钟所能到达距离为半径，画出每个平台所能覆盖到的范围。
但分析到不是所覆盖的范围都可直线通过，而是有很多折线，据此，首先应该根据所给各节点在平面直角坐标系上的坐标资料求出相邻两点之间街道长度，因为所给数据过大，因此，应利用MATLAB软件编写相关程序求出每两个相邻节点之间街道的长度。
再利用DIJKSTRA算法通过MATLAB编写程序计算出每个节点到其他各个节点之间的最短路线的长度。
在此基础上就可以开始问题一各部分求解：
部分一分配各平台的管辖范围，可通过设定目标所有平台所管辖的节点总路程最小（题目要求尽量使得三分钟内有交警到达事发地，实际计算中出现少量三分钟无法到达的点，所以对该约束进行转化使得到每个节点时间最小即可，本题中也即总时间最小即可即总路程最小），以及目标每个平台工作量尽量均衡即工作量最大的平台和工作量最小的交巡警平台工作量之差最小。
即目标有两个，可把问题看做一多目标规划问题。
并设定约束条件每个节点只能分配到一个服务平台，并通过设计LINGO程序求解，得到各节点的分配方案。
接下来再对未管辖的街道进行分配。
原则是：
该平台所管辖所有节点之间街道（即所管辖节点的诱导子图的所有边）一定是归本平台管辖；
把两个端点属于不同平台管理的街道，按工作量均衡的标准分配到工作量少的平台管理，从而得到各平台管辖范围。
实际操作过程中，该模型较难实现，因此，我们的目标函数先忽略工作量均衡的目标，求出分配方案，在遵循工作量均衡原则，对此方案进行调整，得到最终答案。
部分二把各平台与各节点之间的调度看做分配问题，建立分配问题的数学模型，虚拟7个交通要道，并设虚拟的交通要道到各个平台之间的最短距离为0，把各平台到各交通要道的最短距离设为目标矩阵。
通过LINGO软件编程求解，得出最优调度方案。
对于部分三，可通过部分一中求解的范围分配方案求解工作量及出警时间。
再拟设定增加2到5个交巡警服务平台解决以上不太合理的情况。
并再通过新设定的平台通过修改的部分一中模型求解，得出答案，检验是否合理。
（2）问题要求针对全市按设置交巡警服务平台的原则和任务分析现有交巡警服务平台的合理性。
同样应先求出全市各节点之间的最小距离，该方法参照问题
（1）中求解方案。
同样再利用类似于问题一中部分一的求解模型求解。
把全市分成的六个区域分别求解，各区域之间相连的路线参照问题一中部分一边际调整方案分配。
并根据时间要求及工作量要求判定是否合理。
对于围捕犯罪嫌疑人问题，也是指派问题模型。
重要路口是指犯罪嫌疑人以最短路线出A区的路口。
首先判断并封锁这些重要路口，其次将剩余警力分配至A区其余出城区路口，修改第一问建立的指派问题模型求出在A区最短时间围堵罪犯的最佳方案。
如果犯罪嫌疑人已经逃出A区，我们可以利用指派模型，输入全市各个出市路口和最短路线距离矩阵等数据，求出犯罪嫌疑人逃出市区的最佳围堵方案。
四、模型的建立及求解
问题一：
根据坐标通过建立matlab软件求解两相邻节点之间距离，并转化成实际中的距离。
得到该市A区所有节点邻接矩阵
）。
所编写的程序见附录一。
根据上述所求得的数据通过matlab软件利用dijkstra方法求解任意两点之间最短距离。
得到该市A区最短距离矩阵
），程序见附录二。
对于问题一中部分一：
根据题意，需要把92个节点分配给20个交通巡警平台，其中20个节点已经设有交通巡警平台，必定受本节点所设平台管辖，只需将剩余的72个节点进行分配。
而对于节点之间街道的分配，每个平台所管辖的节点之间的街道都由该平台管理。
街道两端节点若被不同平台管理，按工作量大小将该街道分配给工作量小的平台管理。
因此，范围的分配分为两步：
1.节点的分配：
根据交巡警服务平台分配原则建立多目标规划问题如下：
由于该多目标规划较难实现因此，先不考虑工作均衡的目标即第二个目标函数，将该目标函数去掉。
模型简化为：
经lingo软件设计，（程序见附录三）求解得管辖范围分配方案：
平台
所管辖节点
A1
67、68、69、71、74、75、76、78
A11
26、27
A2
39、40、43、44、70、72
A12
25
A3
54、60、62、63、64
A13
21、22、23、24
A4
57、60、62、63、64
A14
A5
49、50、51、52、53、56、58、59
A15
28、29
A6
A16
36、37、38
A7
30、32、47、48、61
A17
41、42
A8
33、46
A18
73、80、81、82、83、84
A9
31、34、35、45
A19
77、79
A10
A20
85、86、87、88、89、90、91、92
表中分配方案为按时间计算管辖最优方案。
只有部分节点即28、29、38、39、61、92节点不能在三分钟之内到达，但选取的仍是最短路程的管辖点。
2.街道的分配：
分配原则：
街道两端节点由同一平台管理的则该街道由管理这两节点的平台管理；
街道两端节点由不同平台管理的由这两平台中工作量少的平台管理。
通过初步分析，我门发现交巡警服务平台的实际工作量包括两部分，即其出警时间及所管辖节点发案率的乘积的累加，和它所管辖节点所有案件的处理时间累加。
因此写出公式
据该公式，得出每个平台工作量情况如下：
平台10.23+9.4t
平台20.43+9.7t
平台30.18+5.6t
平台40.22+6.6t
平台50.45+9.7t
平台60+2.5t
平台70.31+9.6t
平台80.08+5t
平台90.21+8.2t
平台100+1.6t
平台110.11+4.6t
平台120.1+4t
平台130.3+8.5t
平台140+2.5t
平台150.47+4.8t
平台160.16+5t
平台170.09+5.3t
平台180.22+8t
平台190.04+3.4t
平台200.36+10.5t
现根据假设给处理时间赋值，令t=0.5得出各个平台工作量的值（实际生活中t值可调整）
平台14.93
平台15.28
平台12.98
平台13.52
平台15.3
平台11.25
平台15.11
平台12.58
平台14.31
平台10.8
平台12.41
平台12.1
平台14.55
平台12.87
平台12.66
平台12.74
平台14.22
平台11.74
平台15.61
根据以上值给出街道两端由不同平台管理的街道分配放案，即交界处分配方案：
8号平台管理的交界处街道有32~33、8~47、46~55、33~34、46~46、8~9;
9号平台管理的交界处街道有31~32；
10号平台管理的交界处街道有9~10、10~11；
11号平台管理的交界处街道有11~22；
12号平台管理的交界处街道有12~27、11~25、24~25；
14号平台管理的交界处街道有21~14、16~14；
15号平台管理的交界处街道有29~30、7~15、15~31；
16号平台管理的交界处街道有35~36、34~37、36~39、38~39；
17号平台管理的交界处街道有42~43、17~40、38~41、42~81；
18号平台管理的交界处街道有73~74、74~80、84~85、84~89；
19号平台管理的交界处街道有76~77、77~78、78~79、79~80；
20号平台管理的交界处街道有41~92、82~90。
对于问题一中部分二：
若出现重大交通事件，要求调度全区20个交巡警服务平台的警力资源对进出该区的13条交通要道实现快速全封锁，一个平台的警力最多封锁一个路口。
把该问题看作一分配问题，建立0~1规划的数学模型，求最优解。
建立分配问题的数学模型：
由于关键路线的节点不是连续的，因此，设立新变量
其中y中的i仍表示第i号平台，j取值为1~20，1表示调至12号节点，2表示调至14号节点，3表示调至16号节点，4表示调至21号节点，5表示调至22号节点，6表示调至23号节点，7表示调至24号节点，8表示调至28号节点，9表示调至29号节点，10表示调至30号节点，11表示调至38号节点，12表示调至48号节点，13表示调至62号节点，14~20表示虚拟的节点。
所建立数学模型如下：
利用lingo软件编写程序，所编写程序见附录四。
求解得到分配问题的解为：
y（1,12）1.0000000.000000
y（2,16）1.0000002.968200
y（3,9）1.0000001.532600
y（4,14）1.0000003.265000
y（5,10）1.0000007.708000
y（6,13）1.0000000.5000000
y（7,11）1.0000003.805300
y（8,15）1.0000004.751800
y（9,8）1.00000010.49320
y（10,7）1.0000000.5831000
y（11,2）1.0000003.982200
y（12,5）1.0000002.475800
y（13,4）1.0000000.3500000
y（14,6）1.0000000.000000
y（15,3）1.0000000.000000
y（16,19）1.0000000.000000
y（17,1）1.0000000.000000
y（18,18）1.0000000.000000
y（19,20）1.0000000.000000
y（20,17）1.0000000.000000
即调用方案及所用时间（分钟）如下
12→226.8825
16→142.9682
9→161.5326
14→213.2650
10→127.5866
13→230.5
11→243.8053
15→284.7518
8→303.0608
7→298.0155
2→383.9822
5→482.4578
4→620.3600
对于问题一中部分三，要求添加2~5个交巡警平台解决实际中工作量不均衡和出警时间过长的问题。
由于交巡警服务平台的实际工作量包括两部分即其出警时间及所管辖节点发案率的乘积的累加，和它所管辖节点所有案件的处理时间累加。
因此，我们所设的评价是否更为合理的标准即W是否更小。
调整原则
Step1交巡警平台必须在三分钟内赶到事发地
Step2工作量均衡（通过Excel表格处理，Matlab软件处理矩阵见附录五得出如下表格）
1
2
3
4
5
6
7
8
9
10
工作量
4.9927
5.2791
2.9784
3.5160
5.3048
1.25
5.1099
2.5758
4.3090
0.8000
11
12
13
14
15
16
17
18
19
20
2.4134
2..0954
4.5461
2.8719
2.6623
2.7356
4.2185
1.7382
5.6076
由上述原则，step1中28,29,38,39,61,92这六个路口明显三分钟交警不能到达。
Step2中1，2，5，7，13，18，20这五个平台工作量太多。
综上我们可以增加五个平台
具体位置为21，38，61,92,28，（见附图一）
B1
101~103
B2
104~112、117~123、
B3
113~116、126、128、129、131、136、154、
B4
124、127、130、133、134、138~142、145~147、150、151
B5
135、137、143、144
B6
155~165
B7
149、152、153
B8
125、132
问题二：
要求判定全市的交巡警平台设置是否合理，该评价标准也建立在分配了交巡警服务平台所管辖范围的基础上，对各交巡警服务平台的是否能在较短时间内到达其案发地及工作量是否均衡来评价其是否合理。
该问题类似于问题一中求解各平台所管辖范围和判断工作量是否合理的问题。
只需求解出各平台所管辖范围，并比较最长出警时间是否大部分都在3分钟之内，以及工作量是否比较均衡。
通过问题一中建立的模型求解的本市各区分配方案如下：
B区分配方案：
C区分配方案
C1
262~265
C2
248~252、255、258~261
C3
189~192
C4
254
C5
222~226、273、276、277、283
C6
215、216、230、231、240、241~244、246、253
C7
217、218、227~229
C8
232~239、245、247
C9
211~214、219~221
C10
183、193~199、
C11
184~188、
C12
200~210
C13
284、286、287
C14
274、275、278~282、285、288~292、295、296
C15
268~270、297~316
C16
266、267、317~319
C17
256、257、271、272、293、294
D区分配方案
D1
347-350370371
D2
351-360368369
D3
367
D4
344345361362
D5
364-366
D6
D7
343346
D8
3373383403413421
D9
329-336339
E区分配方案
E1
E2
437438456
E3
427432-436457
E4
424-426428-431
E5
E6
416
E7
458459
E8
417-423
E9
387-396
E10
397-400405406
E11
401-404407-415
E12
452-455460-464469470
E13
465-468471472
E
升级会员 