2023年12月9日发(作者:2017云南会考数学试卷)
2011 高教社杯全国大学生数学建模竞赛题目
(请先阅读“全国大学生数学建模竞赛论规范式格文”)
A题城市表层土壤重金属污染分析
随着城市经济的快速发展和城市人口的不断增 ,加人类活动对城市环境质量的影响
日显突出。对城市土壤地质环境异常的查证,以及如何应用查证获得的海量数据开料资
展城市环境质量评价,研究人类活动影响下城市地质环境的演变模式,日益成为人们关
注的焦点。
按照功能划分,城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区
等,分别记为 1 类区、 2 类区、⋯ ⋯ 、 5 类区,不同的区域环境受人类活动影响的程度
不同。
现对某城市城区土壤地质环境进行调查。为此,将所考察的城区划分为间距
左右的网格子区域,按照每平方公里
编号,并用 GPS 记录采样点的位置。应用专门仪器测试分析,获得了每个样所含本的
多种化学元素的浓度数据。另一方面,按照 2 公里的间距在那些远离人群及工业活动的
2 列出了 8 种
自然区取样,将其作为该城区表层土壤中元素的背景值。
附件 1 列出了采样点的位置、海拔高度及其所属功能区等信息,附件
主要重金属元素在采样点处的浓度,附件
现要求你们通过数学建模来完成以下任 :务(1) 给出 8 种主要重金属元素在该城区的空间分布,并分析该城区内不同区域重金
属的污染程度。
(2) 通过数据分析,说明重金属污染的主要原 。因(3) 分析重金属污染物的传播特征,由此建立模型,确定污染源的位置。
(4) 分析你所建立模型的优缺点,为更好地研究城市地质环境的演变模式,还应收
集什么信息?有了这些信息,如何建立模型解决问?题3 列出了 8 种主要重金属元素的背景值。
1 公里
1 个采样点对表层土( 0~10 厘米深度)进行取样、 题 目 A 题 城市表层土壤重金属污染分析
摘 要:
本文研究的是某城区警车配置及巡逻方案的制定问题,建立了求解警车巡逻方案的
模型,并在满足 D1的条件下给出了巡逻效果最好的方案。
在设计整个区域配置最少巡逻车辆时,本文设计了算法 1:先将道路离散化成近似
均匀分布的节点,相邻两个节点之间的距离约等于一分钟巡逻路程。由警车的数目 m ,
将全区划分成 m个均匀的分区,从每个分区的中心点出发,找到最近的道路节点,作为
警车的初始位置,由 Floyd 算法算出每辆警车 3 分钟或 2 分钟行驶路程范围内的节点。
考虑区域调整的概率大小和方向不同会影响调整结果, 本文利用模拟退火算法构造出迁
移几率函数,用迁移方向函数决定分区的调整方向。计算能满足 D1 的最小车辆数,即
为该区应该配置的最小警车数目,用 MATLAB计算,得到局部最优解为 13 辆。
在选取巡逻显著性指标时,本文考虑了两个方面的指标:一是全面性,即所有警车
走过的街道节点数占总街道节点数的比例,用两者之比来评价;二是均匀性,即所有警
车经过每个节点数的次数偏离平均经过次数的程度,用方差值来大小评价。
问题三:为简化问题,假设所有警车在同一时刻,大致向同一方向巡逻,运动状态
分为四种:向左,向右,向上,向下,记录每个时刻,警车经过的节点和能够赶去处理
事故的点,最后汇总计算得相应的评价指标。
在考虑巡逻规律隐蔽性要求时,文本将巡逻路线进行随机处理,方向是不确定的, 采用算法 2 进行计算,得出相应巡逻显著指标,当车辆数减少到 10 辆或巡逻速度变大
时,用算法 2 计算巡逻方案和对应的参数,结果见附录所示。
本文最后还考虑到 4 个额外因素,给出每个影响因素的解决方案。
关键词:模拟退火算法; Floyd 算法;离散化
参赛队号 11***02
参赛密码
(由组委会填写)
队员姓名 * 佳 ** 梅 * 巍
一 问题的重述
110警车在街道上巡逻,既能够对违法犯罪分子起到震慑作用,降低犯罪率,又能够增加市民的安全感,同时也加快了接处警时间,提高了反应时效,为社会和谐提供了
有力的保障。
现给出某城市内一区域, 其道路数据和地图数据已知, 该区域内三个重点部位的坐
标分别为:(5112,4806),(9126, 4266 ),(7434 ,1332)。该区域内共有 307个道路
交叉口,为简化问题,相邻两个交叉路口之间的道路近似认为是直线,且所有事发现场
均在下图的道路上。
该市拟增加一批配备有 GPS卫星定位系统及先进通讯设备的 110 警车。设 110 警车
的平均巡逻速度 为 20km/h,接警后的平均行驶速度为 40km/h。警车配置及巡逻方案要
尽量满足以下要求:
D1. 警车在接警后三分钟内赶到现场的比例不低于 90%;而赶到重点部位的时间必须在
两分钟之内。
D2. 使巡逻效果更显著;
D3. 警车巡逻规律应有一定的隐蔽性。
现在我们需要解决以下几个问题:
一. 若要求满足 D1,该区最少需要配置多少辆警车巡逻?
二. 请给出评价巡逻效果显著程度的有关指标。
三.请给出满足 D1且尽量满足 D2条件的警车巡逻方案及其评价指标值。
四. 在第三问的基础上,再考虑 D3条件,给出你们的警车巡逻方案及其评价指标值。
五.如果该区域仅配置 10 辆警车,应如何制定巡逻方案,使 D1、D2尽量得到满足?
六. 若警车接警后的平均行驶速度提高到 50km/h,回答问题三。
七. 你们认为还有哪些因素、哪些情况需要考虑?给出你们相应的解决方案。
二 问题分析
本题为城区道路网络中警车配置及巡逻问题。在进行警车配置时,首先要考虑警车
在接警后在规定时间内赶到现场的比例, 在此条件下, 以车数最少为目标, 建模、求解;
在制定巡逻方案时,要考虑巡逻的效果及隐蔽性问题。
问题一只要求满足 D1,求最少的警车配置数,可以认为警车是不动的,在三分钟
或两分钟内它能到达的区域就是它的覆盖范围。据此,在满足所有街道的覆盖率不低于
90%的条件下,寻找最优解。
问题二要评价巡逻效果,有两个方面需要考虑:一是巡逻的全面性,即经过一段时
间后警车走过的街道数占总街道数的比例;二是巡逻的不均匀性,即经过一段时间后警
车经过每一条街道的次数相差不大,用方差来衡量。
问题三是在满足 D1 的条件上尽量满足问题二所给的指标, 并给出评价方案的指标。
首先找到一组满足 D1 的各警车位置,然后在和各警车位置相连的点中随机寻找一个点,
判断新的点是否满足 D1,如果满足则警车行驶到该点,否则重新寻找,直到满足为止。
一段时间后统计所有车走过的点数及每个点被走过的次数, 用问题二给出的两个指标进
行评价。综合两个指标,可判断此路径的好坏,重复这个过程,直到综合评价指标达到
一个满意的值为止。
问题四增加了隐蔽性要求,首先给出评价隐蔽性的指标,隐蔽性可用路线的随机性
来评价,将它加入到问题三的模型中去进行求解。
问题五限制警车数量为 10,要综合考虑 D1、D2,先分配这 10 辆车使道路的覆盖
率最高,然后按照问题三的步骤进行求解,其中每一步对 D1 的判断只需使道路的覆盖率尽量高即可。
问题六同问题三,只需将车速改为 50km/h 即可。
三 模型的假设
1.
2.
3.
警车都在路上巡逻,巡警去处理案件的时间不考虑;
所有事发现场都在道路上,案件在道路上任一点是等概率发生的;
警车初始停靠点是随机的,但尽量让它们分散分布,一辆警车管辖一个分区;
4. 假定各个划分区域内,较短时间内,最多会发生一个案件;
5. 假设区域内的每条道路都是双行线,不考虑转弯对结果造成的影响;
6. 如果重点部位不在道路上的,假设这些重点部位在离它们最近的道路上;7. 图中水域对巡逻方案没有影响。
四 符号说明
m 表示警车数目
d
表示警车初始停靠点到各道路的最短距离
L
表示整个区域的总道路长度
l
表示不能在 3 分钟内到达的区域的道路的长度
k
表示非重点部位的警车在 3 分钟内不能到达现场的比例
r 表示三分钟内能从接警位置赶到事发现场的最大距离是
n 表示整个区域总的离散点个数
n 表示第
i
区内的节点个数
i
f 表示区内调整函数
1
t 表示模拟退火的时间,表征温度值
f 表示区间调整函数
2
r 表示全面性指标
e 表示不均匀性指标
h
表示综合评价指标
s 表示第
i
辆车经过每条道路的次数
i
s 表示整个区域每条道路经过的平均次数
五 模型的建立与算法的设计 8. 满足 D1时,该区所需要配置的最少警车数目和巡逻方案
5.1.1 满足 D1条件时,区域最少警车的规律
题目要求警车的配置和巡逻方案满足 D1要求时,整个区域所需要配置的警车数目
最少。由假设可知警车都在道路上,且所有事发现场也都在道路上,但区域内总的道路
长度是个定值的;警车在接警后赶到事发现场有时间限制和概率限制:三分钟内赶到普
通区域案发现场的比例不低于 90%,而赶到重点部位的时间必须控制在两分钟之内。 由
此可知每辆警车的管辖范围不会很大,于是考虑将整个区域分成若干个分区,每辆警车
管辖一个分区域。
由上面的分析,求解整个区域的警车数目最少这个问题可转化为求解每一辆警车所
能管辖的街道范围尽量的大。于是我们寻找出使每辆警车管辖的范围尽量大的规律。为
了简化问题,我们不考虑赶到现场的 90%的几率的限制,仅对警车能在三分钟内赶到事
发现场的情况作定性分析,其分析示意图如图 1 所示。警车的初始停靠位置是随机的分
布在道路上的任一节点上,我们假设一辆警车停靠在 A点上。
图1 一辆警车管辖范围分析示意图
由于警车的平均巡逻速度为 20km/h,接警后的平均行驶速度为 40km/h,由于距离
信息比较容易得到,于是我们将时间限制转化为距离限制,这样便于分析和求解。当警
车接警后,在 三分 钟内 能从 接警位置赶到事发 现场 的最 大距 离是 r , 其中
。
3
r 40 2km
60
如图 1 所示,我们设警车初始停靠位置在 A点,A点是道路 1,2,3,4 的道路交叉
口。我们仅以警车在道路 1 巡逻为例来进行分析,警车以
20km / h
的速度在道路 1 上 A
\' \'
到
A
点之间巡逻,
A
与初始停靠点 A 的距离为
xkm
。由于案件有可能在道路上任一点
发生,当警车巡逻到 A 点时,若案发现场在道路 2,3,4 上发生时,警车以 40km/h 的
速度向事发现场行驶,警车能在三分钟内从
\'
A
点赶到现场的最大距离为
(2 x) km
。如
果警车在道路 1 上继续向前行驶,则该警车能在三分钟内赶到现场的距离继续缩小,当\'
\'
警车从初始点向 A点行驶但没有达到
A
点时,此时该警车的最大管辖范围比警车到达
A
点时的最大管辖范围大。为了使警车的管辖范围尽量大,警车的巡逻范围越小越好,当
x 0时,即警车在初始停靠点静止不动时,警车的管辖范围达到最大值
2km
。
图 1 所分析的是特殊的情况,道路 1,2,3,4 对称分布,现在我们来对一般的情
况进行分析,如图 2 所示。
图2.1
图2
图 2.2
一辆警车最大管辖范围分析示意图
图2.1 所示的情况是道路分布不对称,与图 1 相比,图 2.1 所示的道路方向和角度
都发生了改变,图 2.3 中的情形更为复杂。参照对图 1的分析方法,我们分析这两种情
形下,警车巡逻时能在三分钟内赶到现场的最大距离的规律, 我们只分析图 2.2 的情况,
道路 1,2,3,4,5 相交于点 C,同时道路 1 与道路 6 也有个道路交叉口 D, 由于警车
巡逻时是在道路上行驶的,行走的路线是分段直线 , 并不影响路径的长度,所以当警车
巡逻到距离初始停靠点 C点x远处的 D,此时若有案件发生时,该警车要在三分钟内能
赶到现场处理案件, 最大行驶距离在
(2 x) km
之内,如果警车在道路 1 上继续向前行驶,
则该警车能在三分钟内赶到现场的距离继续缩小,当警车没有行驶到 D点时,此时该警
车的最大管辖范围比
(2 x) km
大,为了使警车的管辖范围尽量大, 警车的巡逻范围越小
越好。当
x 0时,即警车静止不动时,一辆警车的管辖范围能达到最大值。
以上分析的仅作定性的分析,对于三个重点部位也可以同理分析,所得的结论是一
致的,以上的分析没有考虑到 90%的到达几率限制,但在设计算法需要充分考虑。
综上所述,当警车静止在初始停靠点时,在三分钟时间限制内,警车能从初始停靠
点赶到事发现场的最大距离为
2km
。
9. 将道路离散化
由于事发现场是等概率地分布在道路上的,由区域地图可以发现,整个区域中的道
路长度不均,为了使计算结果更加精确,可将这些道路离散化。只要选取合适的离散方
案,就能使警车在经过道路上的离散的点时就相当于经过了这条道路。这样,不论是求
解警车初始停靠点还求解警车赶到事发现场所经过的道路时, 所计算得的的结果显然比
仅考虑整条道路的叉路口要精确得多。
区域中共有 307 个道路交叉口,458条道路。我们采用线性插值方法对道路进行离散
化,以
20km / h
的速度行走一分钟的距离作为步长,一分钟时间的选择是参照问题三的结果要求来设定的,步长
b
1
20
60
。用线性插值的方法,从道路的一个方向进
1
km
3
的整数倍, 我们
1
行线性插值, 实现将每条道路离散化的目标, 考虑到有些道路不是
km
3
1
1
与
)
就一般情况进行讨论, 其分析示意图如图 3 所示。道路 AB长度为 n 个
x( x km
km
3
3
长度的和,为了更精确处理 CB段道路,那么就要考虑在 CB之间是否要插入一个新的点,
根据 x的长度不同,其对应的处理方式也有所不同。
图 3 道路离散化分析示意图
引进临界指数
y
,选取
y
大小的准则是使尽量离散化后警车等效的平均巡逻速度和题
目给定的速度(
20km / h
)的差值尽量小,经过计算得
y 0.189km
时,不再插入新的坐
1
标点时能使整个区域的道路离散效果较好。此时,将
km
CB段长度设定为
3
处理,于是
离散后的 AB道路长度会比实际长度短些; 当
x 0.189 Km
时,需要在两个点之间再插入
一点,因为这样处理能使整个区域的整体道路的离散化效果比较理想。如图 3 所示,在
1
的D点处, 这样处理后所得的C与 B间再插入新的坐标点,插入的位置在距 C点
道
km
3
1
路长度比实际长度长了
(
3
x )km
。采用这样的方法进行线性插值,我们使用 MATLAB编
程实现对整个区域道路的离散, 所得的离散结果如图 4 所示,离散后共得到 762个节点,
比原始数据多了 455 个节点,离散后的节点数据见附件中的“ ”。图 4 整个区域离散结果图
采用这种插值方法道路离散后,将直线上的无穷多个点转化有限个点,便于分析问
题和实现相应的算法,由图 4 可知,所取得的整体离散效果还是比较理想的。
10. 分区域求解警车数目的算法设计
考虑到警车配置和巡逻方案需要满足:警车在接警后三分钟内赶到普通部位案发现 场的比例不低于 90%,赶到重点部位必须控制在两分钟之内的要求。设计算法的目标就
是求解出在满足 D1 情况下,总的警车数目最小,即每个区域都尽可能多地覆盖道路节
点。由于警车的初始位置是未知的,我们可设警车初始停靠点在道路上的任一点,即分
布在图 4 所示的 762 个离散点中的某些点节点上, 总体思路是让每两辆车之间尽量分散
地分布,一辆警车管辖一个分区,用这些分区覆盖整个区域。
于是我们设计算法 1,步骤如下所示:
Step1:将整个区域预分配为 m 个分区,每个分区分配一辆警车,警车的初始停靠位置
设在预分配区中心的道路节点上,若区域的中心不在道路节点上,则将警车放
在离中心最近的道路节点上;
Step2:统计分区不能覆盖的节点,调整警车的初始停靠点,使分区覆盖尽可能多的道
路节点,调整分为区内调整和区间调整方案: (1)区内调整按照模拟退火思想
构造的函数,在区间调整调整车辆初始点的位置(后文中有详细说明) ,当分区
内节点数较多时, 调整的概率小些, 分区内节点数较少时, 调整的概率大些,(2)
当区域中存在未被覆盖的节点或节点群 (大于等于三个节点集中在一个范围内)
时,将警车初始位置的调整方向为朝着这些未被覆盖的节点按一定的规则(在算法说明中有详细叙述)移动,同时要保证 3 个重点部位能在 2 分钟之内 100%
到达;
Step3:用 Floyd 算法计算出警车初始停靠点到周边各道路节点的最短距离
d
;
Step4:以 m 个划分区域未覆盖的总的道路长度
l
与整个区域的道路总长度
L
的比值
l
k
L
100%
来表示警车不能 3 分钟内到达现场的概率;
Step5:模拟足够多的次数,若
k 10%
,将车辆数 m 减1,跳转到 Step1;
Step6:计算结束后,比较当
k 10%
时所对应的 m 值, 当m 取得最小值时,记录此时
的区域划分方案, m即为最少的警车数。
对算法的几点说明:
(1)该算法所取的车辆数 m 是由多到少进行计算的, m初始值设为 20,这个值的
选取是根据区域图估算的。
(2) 预分区的优点在于使警车的初始位置尽可能均匀地分散分布, 警车的初始停靠点
在一个分区的中心点附近寻找得到,比起在整个区域随机生成停靠点,计算效率明显得
到提高。
预分配之后,需要对整个区域不断地进行调整,调整时需要考虑调整方向和 调整概
率。
警车调整借鉴的是模拟退火算法的方法, 为了使分区内包含道路节点数较多的分区
的初始停车点调整的概率小些, 而分区内包含道路节点数的少的分区内的初始停车点调
整的概率大些,我们构造了一个调整概率函数
f ,
1
f
1
bm n
a exp(
i
)
t
(1)
(1)式中,
a, b均为常数, m 为整个区域车辆数,
n 为第i
分区内覆盖的节点数, t
i
为时间,同时 t 也能表征模拟退火的温度变化情况:初始温度较高,区域调整速度较快,
随着时间的增加,温度不断下降,区域调整速度逐渐变慢,这个调整速度变化也是比较
符合实际情况的。
由式(1)可以得出调整概率函数
f ,假设在相同的温度 t (时间)的条件下,由于
1
总的车辆数目 m是定值,当 ni
nj
时,即第i
分区内的节点数大于第
j
分区的节点数时,
分区i
调整的概率大些,分区
j
的调整概率小些。分析其原因:当分区内包含了较多的
节点个数时,该分区的警车初始停靠位置选取地比较合适了,而当分区内包含的道路节
点数较少时,说明警车的初始停靠位置没有选好,需要更大概率的调整,这样的结论也
是比较客观的。
对于所有分区外未被覆盖的道路节点和很多节点(称之为节点群) ,用来调整警车位
置迁移的方向,其分析示意图如图 5 所示。调整方案目标是使未被覆盖的节点数尽量的
少。在设计调整方向函数时,需要考虑: (1)节点群内节点的数目; (2)警车距离节点
群的位置。优先考虑距离,所以在公式( 2)中,用距离的平方来描述调整方向函数。
由于某一个区域范围内的未被覆盖节点数,整个区域未被覆盖的节点总数,分区域与未被覆盖的节点或节点群的距离等几个因素会影响到调整的方案,
些因素。于是设计了区间调整函数
f ,
2
所以要综合考虑这
式中,
n 表示第
i
个分区内未被覆盖的节点数, li
表示第
i
分区域与未被覆盖的节
i
点或节点群的距离,
p
表示未被覆盖的节点和节点群个数。
现在简要分析第
i
分区按区间调整函数的调整方案,当某两节点群
等,但是距离不等时,如 li
i, j
的节点数目相
j
方向调整。 l
j
,由区间调整公式可知,该区间向节点群
当某个分区与两个节点群的距离相等,但节点群的内节点个数不相等,如
(4)可知,该分区域会想节点群
j
方向调整。
ni
n 时,由
j
注意在整个调整过程中,调整几率控制是否调整,调整方向函数控制调整的方向,
寻找在这种调整方案下的最优结果。
图 5 调整分区域示意图
(3)在 step3 中,使用 Floyd 算法计算出警车初始停靠点到周边各节点的最短距离
目的是当区域内有情况发生时,警车能在要求的时间限制内到达现场。
(4)为求出较优的警车停靠点,采用模拟退火算法,算出局部最优的方案。
11. 警车的配置和巡逻方案
d
,
使用 MATLAB编程实现算法 1 得到,整个区域配备 13 辆警车,这些警车静止在初始 停靠点时,能满足 D1要求。警车的初始停靠位置分别为道路交叉节点 6,25,30,37,
82,84,110,111,126,214,253,258,278 处。每个警车所管辖的交叉点(原始的 交
叉 节 点 ) 如 图 6 所 示 , 求 解 的 分 区 结 果 见 附 录 所 示 。图 6 满足 D1条件下的区分划分图
13个分区共覆盖了 252个交叉点,另外的 55 个原始交叉点没有被这些分区域覆盖:
137,138,151,159,167,168,170,174,175,186,188,189,211,215,226,
242,255,260,261,262,263,267,270,271,272,275,282,283 ,284,287,
288,289,292,296,297,299,304,305,307。在这种分区方案下,这些点中,每
两 个 相 连 的 点 间 的 道 路 离 散 值 长 度 占 整 个 区 域 总 的 长 度 的 比 值 为
l
k
L
100% 90.18%
。因此,在整个区域配置 13 辆警车,每个警车在初始停靠点静
止不动,当有案件发生时,离案发现场最近的警车从初始停靠点赶到现场。
12. 评价巡逻效果显著的指标
110 警车在街道上巡逻是目的是为了对违法犯罪分子起到震慑作用,降低犯罪率,
又能够增加市民的安全感, 同时还加快了接处警 (接受报警并赶往现场处理事件) 时间,
提高了反应时效,为社会和谐提供了有力的保障。巡警在城市繁华街道、公共场所执行
[1]
巡逻任务, 维护治安, 服务群众 , 可以得良好的社会效应
。
在整个区域中,由于案发现场都在道路上,道路上的每一点都是等概率发生的,因
此警车巡逻的面越广,所巡逻的街道数目越多,警车的巡逻效果就越好,对违法犯罪分
子就越有威慑力,警车也能更及时地处理案件。
我们采用全面性 r 来衡量巡逻的效果显著性,即用警车巡逻所经过的街道节点数占
区域总节点数的比值。当警车重复经过同一条街道同一个离散点时, c仅记录一次。c
r
n
(3)
式中, c 表示警车经过的离散点数, n 代表整个区域总的离散点数。 r 值越大,表
明警车所经过的街道数目越多,所取得的效果越显著。
同时考虑到在巡逻过程中可能会出现这样的情况:在相同的时段内,警车会多次巡
逻部分街道,而一些街道却很少巡逻甚至没有警车到达,这样会造成一些巡逻盲区。分
布很不均衡。这样就可能出现巡逻密度大的街道上的违法犯罪分子不敢在街道上作案,
而流窜到巡逻密度稀疏的街道上作案,因此在相同的警车数目条件下,密度不均衡的巡
逻方式的巡逻效果的效果较差,而密度较均衡的巡逻方式所取得的巡逻效果会更好些。
我们引入一个巡逻的不均匀度 e来衡量巡逻效果的显著性, 考虑到方差能表示不均衡度,
于是我们用方差的大小来表征不均衡,方差越大,巡逻密度越不均衡,所取得的巡逻效
果越差。
m
( s
i
2
s)
i 1
(4)
e
p
问题 1 所给出的满足 D1 条件下的警车数目为 13 辆,这时每辆警车在初始停靠点静
止不动,只有该管辖区域内发生了案件时, 警车才从初始停靠点赶到案发现场处理案件。
当警车在巡逻状态时,所需要考虑的问题就更复杂一些,如当节点运动时,警车还能否
达到 D1的要求,警车的运动方向如何等问题,但基本算法思想与问题 1 类似,所得的
算法 2 的框图如图 7 所示,
为了简化问题,我们假设各分区警车的巡逻时候,尽量保证所有的警车的行驶方向
相一致,且警车都走双行道,即当警车走到某个节点后,它们又同时返回初始停靠点,
警车的行驶方向有四种方式,如 6 所示。
在图 6 中,数字 1 代表走巡逻走的第一步, 2 表示朝 1 的巡逻方向相反的方向巡逻。
在具体程序实现时,四种巡逻方向任意选择,但是尽量保证所有的警车向同一个方向巡
逻。图 6 各警车巡逻方向图
我们用 MATLAB编程对这种巡逻方式进行计算,所得的车辆数目为 18 辆,综合评
价指标为
h 0.612,其结果巡逻方案见附件中的“ ”所示。
13. 在满足问题三的基础上讨论 D3条件,警车的巡逻方案和评价指标
巡逻的隐蔽性体现在警车的巡逻路线和时间没有明显的规律, 主要目的是让违法犯
罪分子无可乘之机,防止他们在非巡逻时间实施违法犯罪活动,危害人民的生命和财产
安全。
为了使巡逻的规律具有隐蔽性,这就需要警车在巡逻时至少具有两条不同的路线,
时间最好也是不相同的。因此,考虑到隐蔽性时,只需要在问题 2 的基础上加上一个随
机过程即可。对于其评价指标,由于警车有几条可选的巡逻路线,当相同的路线在同一
时间内重复出现时,重新将所设定的方案再执行一遍,我们用这个时间间隔来衡量隐蔽
性的程度,当循环周期
T
越大,表明可选的巡逻方案越多,其规律就越具有隐蔽性,而
循环周期
T
越小时,表明巡逻方案比较少,其隐蔽性较差。在巡逻状态时,最差的隐蔽
性巡逻方案是巡逻方案只有一个, 并且时间固定,这样的巡逻方案没有任何隐蔽性可言。
14. 整个区域为 10 辆车时的巡逻方案
由第三问的结果可知, 10辆车的数量是不能把整个区域完全覆盖的, 其算法与算法
2 类似,不同的是此时车的数目已经固定了,要求使 D1,D2尽量大的满足,我们求得的
评价指标值为
h 0.524,所得的巡逻方案见附件中的“ ”所示。
15. 平均行驶速度提高到
50km / h
时的巡逻方式和评价指标值
问题六的分析方法与具体实现与问题三一致,但是警车的接警后的平均速度由原来
的40km / h提高到
50km / h
,于是各分区的覆盖范围也增大了, 将数值带入问题 3 的算法
中求解, 计算得的指标值为
h 0.703
,其巡逻方案见附件中的 “ ”
所示。图7 算法 2 框图六 模型的分析和评价
在求解满足 D1 的条件下,整个区域需要配备多少辆警车问题中,采用分区巡逻的
思想,先分析能使各区管辖范围达到最大值时的规律,由特殊到一般层层进行分析,逻
辑严密,结果合理。
在求解区域和警车数目时,在初步设定警车停靠点位置的基础上,用模拟退火算法
思路构造函数
f 来确定调整的概率大小, 综合考虑了影响区间调整的因素后构造了
1
f
2
函
数来确定分区的调整方向,当分区按照这两个调整函数进行调整时,各分区能管辖尽可
能多的道路节点,所取得效果也比较理想。
参 考 文 献
[1] 中小城市警察巡逻勤务方式的探讨,俞详,江苏公安专科学校学报, 1998 年第 1 期
[2]Matlab7.0 从入门到精通,求是科技,人民邮电出版社;
[3] 不确定车数的随机车辆路径问题模型及算法, 运怀立等,工业工程, 第 10 卷第 3 期,
2005 年 5 月;
[4] 随机交通分配中的有效路径的确定方法, 李志纯等, 交通运输系统工程与信息, 第 3
卷第 1 期,2003 年 2 月。
附 录图 问题三巡逻路径
图 问题五巡逻路径图 问题六巡逻路径
更多推荐
警车,巡逻,道路,调整
发布评论