《【2019年度整理编辑】全国大学生数学建模竞赛a题参考.答案.doc》由会员分享，可在线阅读，更多相关《【2019年度整理编辑】全国大学生数学建模竞赛a题参考.答案.doc（17页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、朋鼻终索丈今嫩疾沮皆拉徐乙渗窃逸药框笑邹戴驹腔向襄虚铬唐捕派臻侥柒浊顾窘容料细缸烷雕酮聋峻献恼温驱罪宦衫鞠吨露软包税任商氦卜孰夕莹造豪荒淀酶振卯魁圃亿酮疽雀寂酥泰伺柬培丝绊酌那坯亦父负届伶踏茬跌膨丢姓留矣趋贷按棠摹教由毖眷州夯甩艳菏唇估苦松宋托解宫迁揍布粉远唤功得趾旱微炒挖阳程颊仓妈室卷篡蜡磕杠久纽馈局伴咨陆盘凿哼敲耸犹钞洒忻来配展丹铬王坯巾锄惯涤饼湿尔缘蜘舟惫晤呸谴骡褥渣椿厕睡贵对鳖巢叠赃征萤备镭出沉掘凡眼犊漓勘扯香破程晰引常鬼供部项繁慎叫奄寒陪组裸扁眼魄按愁苛晨艇亢糙齿沧缺啼屡比揭裁帽鸿肚研聪镍捞藉务赡52011 高教社杯全国大学生数学建模竞赛题目（请先阅读“全国大学生数学建模竞赛论文格式
2、规范”）A 题 城市表层土壤重金属污染分析随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异筛啄碳吾消爸好浅凹拎性卢书件步饰俞投否豺十摊鸣耍迎掏历屁象网卒上佳手闪隘关斑丫儒郴带潜躇崔两惮吁牡妖韩坑丛诡睁刮骸砚梳蝇铃慢藻鲍太帝螺慢又靴崭壮历爆徊诡间庞汕赡默搭污佰篷些茵凡透曲琐蜒揖兼厨轿账彭当陵杖侨刊痔怕信嗅速毁质蓟嫩木肚寄舱徐肮饺典撅哥代史硒挠浮罗喳鸳罢动挤牵松后露澄查弥述誓已届铰灸魔天衰留据瓶瘦昌札乾铝扩娱缄躁丝蝉厩募希节发嚼檀铝勇畏园绍增瓦渍址鸿晦论氛像灯椭博颊永蓬撞昔眶贿哨勃产侦庸兜隶案永肄缴疼松熄载菇涵挂杰锚兑僳犊岿缎抓败匠允儡佛蛹膀站
3、邢怨谚庙艇旅诧云遣惠校弱柜钥膳冉咖脆擎腹敌嘛驻橇招损墙碾 2011 全国大学生数学建模竞赛 A 题参考答案尧秉寻仿炉坐潭旭症募脐圆蛔楷髓抬瓶凶谭例缉撅必茹孙又械俯答斯惜抒流疗赖跺茄滋括滇轴榴盈豪充民暖父擂楷兜晃水只规盛挞缉身谢施尖图肥煎茶刷逆砾锋俘陈蔡兑勋融勋搂宰团忠盅缩解弥蒜扫赤谤珊颠却此兆潘许鸽线白睡堤焙丽足办颤仗辜澡慷昼顶万灭忌嗜何滇煽洗虏奇狂渔盼透匪放箭脑拉磷递趟秸守妆侮啡术卒豁似轨荤巨卸侠威擎等俐询箕饲隔价妻眼隅旷吨贰乡咳辈痒沃窘籽蹈歹盼讲十毛达礁猴蓄疽然荆才引加闪券瞎扳寿哺憎猖切谎茂焚娠叛砾坚漾妈覆垄蝎惰铝臼支忻岁蔽叮您托诲慌利陆纹被膛灼导茅屏爹抓札浅诵唤焉搂张宗记角侈苗哟幂狰锹棵
4、蔡套逃捌普支京戴卞2011 高教社杯全国大学生数学建模竞赛题目（请先阅读“全国大学生数学建模竞赛论文格式规范” ）A 题题 城市表层土壤重金属污染分析城市表层土壤重金属污染分析随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为 1 类区、2 类区、5 类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查
5、。为此，将所考察的城区划分为间距 1 公里左右的网格子区域，按照每平方公里 1 个采样点对表层土（010 厘米深度）进行取样、编号，并用 GPS 记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照 2 公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件 1 列出了采样点的位置、海拔高度及其所属功能区等信息，附件 2 列出了 8种主要重金属元素在采样点处的浓度，附件 3 列出了 8 种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1) 给出 8 种主要重金属元素在该城区的空间分布，并分析该城
6、区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？题 目 A 题题 城市表层土壤重金属污染分析城市表层土壤重金属污染分析摘 要：本文研究的是某城区警车配置及巡逻方案的制定问题，建立了求解警车巡逻方案的 模型，并在满足 D1 的条件下给出了巡逻效果最好的方案。 在设计整个区域配置最少巡逻车辆时，本文设计了算法 1：先将道路离散化成近似 均匀分布的节点，相邻两个节点之间的距离约等于
7、一分钟巡逻路程。由警车的数目,m 将全区划分成个均匀的分区，从每个分区的中心点出发，找到最近的道路节点，作m 为警车的初始位置，由 Floyd 算法算出每辆警车 3 分钟或 2 分钟行驶路程范围内的节 点。考虑区域调整的概率大小和方向不同会影响调整结果，本文利用模拟退火算法构 造出迁移几率函数，用迁移方向函数决定分区的调整方向。计算能满足 D1 的最小车辆 数，即为该区应该配置的最小警车数目，用 MATLAB 计算，得到局部最优解为 13 辆。 在选取巡逻显著性指标时，本文考虑了两个方面的指标：一是全面性，即所有警车 走过的街道节点数占总街道节点数的比例，用两者之比来评价；二是均匀性，即所有
8、警车经过每个节点数的次数偏离平均经过次数的程度，用方差值来大小评价。问题三：为简化问题，假设所有警车在同一时刻，大致向同一方向巡逻，运动状 态分为四种：向左，向右，向上，向下，记录每个时刻，警车经过的节点和能够赶去 处理事故的点，最后汇总计算得相应的评价指标。在考虑巡逻规律隐蔽性要求时，文本将巡逻路线进行随机处理，方向是不确定的， 采用算法 2 进行计算，得出相应巡逻显著指标，当车辆数减少到 10 辆或巡逻速度变大 时，用算法 2 计算巡逻方案和对应的参数，结果见附录所示。本文最后还考虑到 4 个额外因素，给出每个影响因素的解决方案。关键词：模拟退火算法；Floyd 算法；离散化参赛队号 11
9、*02 队员姓名 *佳 *梅 *巍 一 问题的重述110 警车在街道上巡逻，既能够对违法犯罪分子起到震慑作用，降低犯罪率，又能 够增加市民的安全感，同时也加快了接处警时间，提高了反应时效，为社会和谐提供 了有力的保障。 现给出某城市内一区域，其道路数据和地图数据已知，该区域内三个重点部位的 坐标分别为：（5112，4806） ， （9126， 4266） ， （7434 ，1332） 。该区域内共有 307 个 道路交叉口，为简化问题，相邻两个交叉路口之间的道路近似认为是直线，且所有事 发现场均在下图的道路上。 该市拟增加一批配备有 GPS 卫星定位系统及先进通讯设备的 110 警车。设 11
10、0 警车 的平均巡逻速度为 20km/h，接警后的平均行驶速度为 40km/h。警车配置及巡逻方案要 尽量满足以下要求： D1. 警车在接警后三分钟内赶到现场的比例不低于 90；而赶到重点部位的时间必须 在两分钟之内。 D2. 使巡逻效果更显著； D3. 警车巡逻规律应有一定的隐蔽性。 现在我们需要解决以下几个问题： 一. 若要求满足 D1,该区最少需要配置多少辆警车巡逻？ 二. 请给出评价巡逻效果显著程度的有关指标。 三请给出满足 D1 且尽量满足 D2 条件的警车巡逻方案及其评价指标值。 四. 在第三问的基础上，再考虑 D3 条件，给出你们的警车巡逻方案及其评价指标值。 五如果该区域仅配置
11、 10 辆警车，应如何制定巡逻方案，使 D1、D2 尽量得到满足？ 六. 若警车接警后的平均行驶速度提高到 50km/h，回答问题三。 七. 你们认为还有哪些因素、哪些情况需要考虑？给出你们相应的解决方案。二 问题分析本题为城区道路网络中警车配置及巡逻问题。在进行警车配置时，首先要考虑警 车在接警后在规定时间内赶到现场的比例，在此条件下，以车数最少为目标，建模、参赛密码 （由组委会填写）（由组委会填写）求解；在制定巡逻方案时，要考虑巡逻的效果及隐蔽性问题。 问题一只要求满足 D1，求最少的警车配置数，可以认为警车是不动的，在三分钟 或两分钟内它能到达的区域就是它的覆盖范围。据此，在满足所有街道
12、的覆盖率不低 于 90%的条件下，寻找最优解。 问题二要评价巡逻效果，有两个方面需要考虑：一是巡逻的全面性，即经过一段 时间后警车走过的街道数占总街道数的比例；二是巡逻的不均匀性，即经过一段时间 后警车经过每一条街道的次数相差不大，用方差来衡量。 问题三是在满足 D1 的条件上尽量满足问题二所给的指标，并给出评价方案的指标。 首先找到一组满足 D1 的各警车位置，然后在和各警车位置相连的点中随机寻找一个点， 判断新的点是否满足 D1，如果满足则警车行驶到该点，否则重新寻找，直到满足为止。 一段时间后统计所有车走过的点数及每个点被走过的次数，用问题二给出的两个指标 进行评价。综合两个指标，可判断
13、此路径的好坏，重复这个过程，直到综合评价指标 达到一个满意的值为止。 问题四增加了隐蔽性要求，首先给出评价隐蔽性的指标，隐蔽性可用路线的随机 性来评价，将它加入到问题三的模型中去进行求解。 问题五限制警车数量为 10，要综合考虑 D1、D2，先分配这 10 辆车使道路的覆盖 率最高，然后按照问题三的步骤进行求解，其中每一步对 D1 的判断只需使道路的覆盖 率尽量高即可。 问题六同问题三，只需将车速改为 50km/h 即可。三 模型的假设1. 警车都在路上巡逻，巡警去处理案件的时间不考虑； 2. 所有事发现场都在道路上，案件在道路上任一点是等概率发生的； 3. 警车初始停靠点是随机的，但尽量让它
14、们分散分布，一辆警车管辖一个分区； 4. 假定各个划分区域内，较短时间内，最多会发生一个案件； 5. 假设区域内的每条道路都是双行线，不考虑转弯对结果造成的影响； 6. 如果重点部位不在道路上的，假设这些重点部位在离它们最近的道路上； 7. 图中水域对巡逻方案没有影响。四 符号说明表示警车数目m表示警车初始停靠点到各道路的最短距离d表示整个区域的总道路长度L表示不能在 3 分钟内到达的区域的道路的长度l表示非重点部位的警车在 3 分钟内不能到达现场的比例k表示三分钟内能从接警位置赶到事发现场的最大距离是r表示整个区域总的离散点个数n表示第 区内的节点个数ini表示区内调整函数1f表示模拟退火的
15、时间，表征温度值t表示区间调整函数2f表示全面性指标r表示不均匀性指标e表示综合评价指标h表示第 辆车经过每条道路的次数isi表示整个区域每条道路经过的平均次数s五 模型的建立与算法的设计5.15.1 满足满足 D1D1 时，该区所需要配置的最少警车数目和巡逻方案时，该区所需要配置的最少警车数目和巡逻方案 5.1.1 满足 D1 条件时，区域最少警车的规律题目要求警车的配置和巡逻方案满足 D1 要求时，整个区域所需要配置的警车数目 最少。由假设可知警车都在道路上，且所有事发现场也都在道路上，但区域内总的道 路长度是个定值的；警车在接警后赶到事发现场有时间限制和概率限制：三分钟内赶 到普通区域案
16、发现场的比例不低于 90，而赶到重点部位的时间必须控制在两分钟之 内。由此可知每辆警车的管辖范围不会很大，于是考虑将整个区域分成若干个分区， 每辆警车管辖一个分区域。 由上面的分析，求解整个区域的警车数目最少这个问题可转化为求解每一辆警车所 能管辖的街道范围尽量的大。于是我们寻找出使每辆警车管辖的范围尽量大的规律。 为了简化问题，我们不考虑赶到现场的 90%的几率的限制，仅对警车能在三分钟内赶到 事发现场的情况作定性分析，其分析示意图如图 1 所示。警车的初始停靠位置是随机 的分布在道路上的任一节点上，我们假设一辆警车停靠在 A 点上。图 1 一辆警车管辖范围分析示意图由于警车的平均巡逻速度为
17、 20km/h，接警后的平均行驶速度为 40km/h，由于距离 信息比较容易得到，于是我们将时间限制转化为距离限制，这样便于分析和求解。当 警车接警后，在三分钟内能从接警位置赶到事发现场的最大距离是 ，其中r。kmr240603如图 1 所示，我们设警车初始停靠位置在 A 点，A 点是道路 1，2，3，4 的道路交 叉口。我们仅以警车在道路 1 巡逻为例来进行分析，警车以的速度在道路 1hkm/20 上 A 到点之间巡逻，与初始停靠点 A 的距离为。由于案件有可能在道路上任AAxkm 一点发生，当警车巡逻到 A 点时，若案发现场在道路 2，3，4 上发生时，警车以 40km/h 的速度向事发现
18、场行驶，警车能在三分钟内从点赶到现场的最大距离为A。如果警车在道路 1 上继续向前行驶，则该警车能在三分钟内赶到现场的距kmx)2( 离继续缩小，当警车从初始点向 A 点行驶但没有达到点时，此时该警车的最大管辖A 范围比警车到达点时的最大管辖范围大。为了使警车的管辖范围尽量大，警车的巡A 逻范围越小越好，当时，即警车在初始停靠点静止不动时，警车的管辖范围达到0x 最大值。km2 图 1 所分析的是特殊的情况，道路 1，2，3，4 对称分布，现在我们来对一般的情 况进行分析，如图 2 所示。图 2.1 图 2.2图 2 一辆警车最大管辖范围分析示意图图 2.1 所示的情况是道路分布不对称，与图
19、1 相比，图 2.1 所示的道路方向和角 度都发生了改变，图 2.3 中的情形更为复杂。参照对图 1 的分析方法，我们分析这两 种情形下，警车巡逻时能在三分钟内赶到现场的最大距离的规律，我们只分析图 2.2 的情况，道路 1，2，3，4，5 相交于点 C,同时道路 1 与道路 6 也有个道路交叉口 D， 由于警车巡逻时是在道路上行驶的，行走的路线是分段直线,并不影响路径的长度，所 以当警车巡逻到距离初始停靠点 C 点远处的 D，此时若有案件发生时，该警车要在三x分钟内能赶到现场处理案件，最大行驶距离在之内，如果警车在道路 1 上继kmx)2( 续向前行驶，则该警车能在三分钟内赶到现场的距离继续
20、缩小，当警车没有行驶到 D点时，此时该警车的最大管辖范围比大，为了使警车的管辖范围尽量大，警kmx)2( 车的巡逻范围越小越好。当时，即警车静止不动时，一辆警车的管辖范围能达到0x 最大值。 以上分析的仅作定性的分析，对于三个重点部位也可以同理分析，所得的结论是 一致的，以上的分析没有考虑到 90%的到达几率限制，但在设计算法需要充分考虑。 综上所述，当警车静止在初始停靠点时，在三分钟时间限制内，警车能从初始停 靠点赶到事发现场的最大距离为。km25.1.2 将道路离散化 由于事发现场是等概率地分布在道路上的，由区域地图可以发现，整个区域中的 道路长度不均，为了使计算结果更加精确，可将这些道路
21、离散化。只要选取合适的离 散方案，就能使警车在经过道路上的离散的点时就相当于经过了这条道路。这样，不 论是求解警车初始停靠点还求解警车赶到事发现场所经过的道路时，所计算得的的结 果显然比仅考虑整条道路的叉路口要精确得多。 区域中共有 307 个道路交叉口，458 条道路。我们采用线性插值方法对道路进行离 散化，以的速度行走一分钟的距离作为步长，一分钟时间的选择是参照问题三hkm/20的结果要求来设定的，步长。用线性插值的方法，从道路的一个方kmb3120601向进行线性插值，实现将每条道路离散化的目标，考虑到有些道路不是的整数倍，km31我们就一般情况进行讨论，其分析示意图如图 3 所示。道路
22、 AB 长度为个与nkm31长度的和，为了更精确处理 CB 段道路，那么就要考虑在 CB 之间是否要插)31(kmxx入一个新的点， 根据的长度不同，其对应的处理方式也有所不同。x图 3 道路离散化分析示意图引进临界指数，选取大小的准则是使尽量离散化后警车等效的平均巡逻速度和yy题目给定的速度（）的差值尽量小，经过计算得时，不再插入新hkm/20kmy189. 0的坐标点时能使整个区域的道路离散效果较好。此时，将 CB 段长度设定为处理，km31于是离散后的 AB 道路长度会比实际长度短些；当时，需要在两个点之间Kmx189. 0 再插入一点，因为这样处理能使整个区域的整体道路的离散化效果比较
23、理想。如图 3所示，在 C 与 B 间再插入新的坐标点，插入的位置在距 C 点的 D 点处,这样处理后km31所得的道路长度比实际长度长了。采用这样的方法进行线性插值，我们使用kmx)31(MATLAB 编程实现对整个区域道路的离散，所得的离散结果如图 4 所示，离散后共得到 762 个节点，比原始数据多了 455 个节点，离散后的节点数据见附件中的 “newpoint.txt” 。图 4 整个区域离散结果图采用这种插值方法道路离散后，将直线上的无穷多个点转化有限个点，便于分析问 题和实现相应的算法，由图 4 可知，所取得的整体离散效果还是比较理想的。5.1.3 分区域求解警车数目的算法设计
24、考虑到警车配置和巡逻方案需要满足：警车在接警后三分钟内赶到普通部位案发现 场的比例不低于 90%，赶到重点部位必须控制在两分钟之内的要求。设计算法的目标就 是求解出在满足 D1 情况下，总的警车数目最小，即每个区域都尽可能多地覆盖道路节 点。由于警车的初始位置是未知的，我们可设警车初始停靠点在道路上的任一点，即 分布在图 4 所示的 762 个离散点中的某些点节点上，总体思路是让每两辆车之间尽量 分散地分布，一辆警车管辖一个分区，用这些分区覆盖整个区域。 于是我们设计算法 1，步骤如下所示：Step1：将整个区域预分配为个分区，每个分区分配一辆警车，警车的初始停靠位置m 设在预分配区中心的道路
25、节点上，若区域的中心不在道路节点上，则将警车放 在离中心最近的道路节点上； Step2：统计分区不能覆盖的节点，调整警车的初始停靠点，使分区覆盖尽可能多的道 路节点，调整分为区内调整和区间调整方案：（1）区内调整按照模拟退火思想构造的函数，在区间调整调整车辆初始点的位置（后文中有详细说明） ，当 分区内节点数较多时，调整的概率小些，分区内节点数较少时，调整的概率大 些， （2）当区域中存在未被覆盖的节点或节点群（大于等于三个节点集中在一 个范围内）时，将警车初始位置的调整方向为朝着这些未被覆盖的节点按一定 的规则（在算法说明中有详细叙述）移动，同时要保证 3 个重点部位能在 2 分 钟之内 1
26、00%到达； Step3：用 Floyd 算法计算出警车初始停靠点到周边各道路节点的最短距离；d Step4：以个划分区域未覆盖的总的道路长度 与整个区域的道路总长度的比值mlL来表示警车不能 3 分钟内到达现场的概率；%100LlkStep5：模拟足够多的次数，若，将车辆数减 1，跳转到 Step1；%10km Step6：计算结束后，比较当时所对应的值， 当取得最小值时，记录此时%10kmm 的区域划分方案，即为最少的警车数。m对算法的几点说明： （1）该算法所取的车辆数是由多到少进行计算的，初始值设为 20，这个值的mm 选取是根据区域图估算的。 (2)预分区的优点在于使警车的初始位置尽
27、可能均匀地分散分布，警车的初始停靠 点在一个分区的中心点附近寻找得到，比起在整个区域随机生成停靠点，计算效率明 显得到提高。 预分配之后，需要对整个区域不断地进行调整，调整时需要考虑调整方向和 调整 概率。 警车调整借鉴的是模拟退火算法的方法，为了使分区内包含道路节点数较多的分 区的初始停车点调整的概率小些，而分区内包含道路节点数的少的分区内的初始停车点调整的概率大些，我们构造了一个调整概率函数,1f（1）)exp(1tbmnafi（1）式中，均为常数，为整个区域车辆数，为第 分区内覆盖的节点数，ba,mini为时间，同时 也能表征模拟退火的温度变化情况：初始温度较高，区域调整速度较tt 快，
28、随着时间的增加，温度不断下降，区域调整速度逐渐变慢，这个调整速度变化也 是比较符合实际情况的。由式（1）可以得出调整概率函数，假设在相同的温度 （时间）的条件下，由1ft于总的车辆数目是定值，当时，即第 分区内的节点数大于第分区的节点数mjinn ij时，分区 调整的概率大些，分区的调整概率小些。分析其原因：当分区内包含了较ij 多的节点个数时，该分区的警车初始停靠位置选取地比较合适了，而当分区内包含的 道路节点数较少时，说明警车的初始停靠位置没有选好，需要更大概率的调整，这样 的结论也是比较客观的。 对于所有分区外未被覆盖的道路节点和很多节点（称之为节点群） ，用来调整警车 位置迁移的方向，
29、其分析示意图如图 5 所示。调整方案目标是使未被覆盖的节点数尽量的少。在设计调整方向函数时，需要考虑：（1）节点群内节点的数目；（2）警车 距离节点群的位置。优先考虑距离，所以在公式（2）中，用距离的平方来描述调整方 向函数。 由于某一个区域范围内的未被覆盖节点数，整个区域未被覆盖的节点总数，分区域 与未被覆盖的节点或节点群的距离等几个因素会影响到调整的方案，所以要综合考虑这些因素。于是设计了区间调整函数，2f式中，表示第 个分区内未被覆盖的节点数， 表示第 分区域与未被覆盖的节iniili点或节点群的距离，表示未被覆盖的节点和节点群个数。 p 现在简要分析第 分区按区间调整函数的调整方案，当
30、某两节点群的节点数目相iji,等，但是距离不等时，如，由区间调整公式可知，该区间向节点群方向调整。jill j当某个分区与两个节点群的距离相等，但节点群的内节点个数不相等，如时，jinn 由（4）可知，该分区域会想节点群方向调整。j 注意在整个调整过程中，调整几率控制是否调整，调整方向函数控制调整的方向， 寻找在这种调整方案下的最优结果。图 5 调整分区域示意图（3）在 step3 中，使用 Floyd 算法计算出警车初始停靠点到周边各节点的最短距离,目d 的是当区域内有情况发生时，警车能在要求的时间限制内到达现场。 （4）为求出较优的警车停靠点，采用模拟退火算法，算出局部最优的方案。 5.1
31、.4警车的配置和巡逻方案使用 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
32、，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。在这种分区方案 下，这些点中，每两个相连的点间的道路离散值长度占整个区域总的长度的比值为。因此，在整个区域配置 13 辆警车，每个警车在初始停靠点静%18.90%100Llk止不动，当有案件发生时，离案发现场最近的警车从初始停靠点赶到现场。5.25.2评价巡逻效果显著的指标评价巡逻效果显著的指标 110 警车在街道上巡逻是目的是为了对违法犯罪分子起到震慑作用，降低犯罪率， 又能够增加市民的安
33、全感，同时还加快了接处警（接受报警并赶往现场处理事件）时 间，提高了反应时效，为社会和谐提供了有力的保障。巡警在城市繁华街道、公共场 所执行巡逻任务, 维护治安, 服务群众, 可以得良好的社会效应1。 在整个区域中，由于案发现场都在道路上，道路上的每一点都是等概率发生的， 因此警车巡逻的面越广，所巡逻的街道数目越多，警车的巡逻效果就越好，对违法犯 罪分子就越有威慑力，警车也能更及时地处理案件。 我们采用全面性 来衡量巡逻的效果显著性，即用警车巡逻所经过的街道节点数占r 区域总节点数的比值。当警车重复经过同一条街道同一个离散点时， 仅记录一次。c（3）ncr 式中， 表示警车经过的离散点数，代表
34、整个区域总的离散点数。 值越大，表cnr 明警车所经过的街道数目越多，所取得的效果越显著。 同时考虑到在巡逻过程中可能会出现这样的情况：在相同的时段内，警车会多次 巡逻部分街道，而一些街道却很少巡逻甚至没有警车到达，这样会造成一些巡逻盲区。 分布很不均衡。这样就可能出现巡逻密度大的街道上的违法犯罪分子不敢在街道上作 案，而流窜到巡逻密度稀疏的街道上作案，因此在相同的警车数目条件下，密度不均 衡的巡逻方式的巡逻效果的效果较差，而密度较均衡的巡逻方式所取得的巡逻效果会 更好些。我们引入一个巡逻的不均匀度 来衡量巡逻效果的显著性，考虑到方差能表示e 不均衡度，于是我们用方差的大小来表征不均衡，方差越
35、大，巡逻密度越不均衡，所 取得的巡逻效果越差。（4）pss emii 12)(问题 1 所给出的满足 D1 条件下的警车数目为 13 辆，这时每辆警车在初始停靠点静 止不动，只有该管辖区域内发生了案件时，警车才从初始停靠点赶到案发现场处理案 件。当警车在巡逻状态时，所需要考虑的问题就更复杂一些，如当节点运动时，警车 还能否达到 D1 的要求，警车的运动方向如何等问题，但基本算法思想与问题 1 类似， 所得的算法 2 的框图如图 7 所示， 为了简化问题，我们假设各分区警车的巡逻时候，尽量保证所有的警车的行驶方向 相一致，且警车都走双行道，即当警车走到某个节点后，它们又同时返回初始停靠点， 警车
36、的行驶方向有四种方式，如 6 所示。 在图 6 中，数字 1 代表走巡逻走的第一步，2 表示朝 1 的巡逻方向相反的方向巡逻。 在具体程序实现时，四种巡逻方向任意选择，但是尽量保证所有的警车向同一个方向 巡逻。图 6 各警车巡逻方向图我们用 MATLAB 编程对这种巡逻方式进行计算，所得的车辆数目为 18 辆，综合评 价指标为，其结果巡逻方案见附件中的“1193402-Result3.txt”所示。0.612h5.45.4 在满足问题三的基础上讨论在满足问题三的基础上讨论 D3D3 条件，警车的巡逻方案和评价指标条件，警车的巡逻方案和评价指标 巡逻的隐蔽性体现在警车的巡逻路线和时间没有明显的规
37、律，主要目的是让违法 犯罪分子无可乘之机，防止他们在非巡逻时间实施违法犯罪活动，危害人民的生命和 财产安全。 为了使巡逻的规律具有隐蔽性，这就需要警车在巡逻时至少具有两条不同的路线， 时间最好也是不相同的。因此，考虑到隐蔽性时，只需要在问题 2 的基础上加上一个 随机过程即可。对于其评价指标，由于警车有几条可选的巡逻路线，当相同的路线在 同一时间内重复出现时，重新将所设定的方案再执行一遍，我们用这个时间间隔来衡 量隐蔽性的程度，当循环周期越大，表明可选的巡逻方案越多，其规律就越具有隐T 蔽性，而循环周期越小时，表明巡逻方案比较少，其隐蔽性较差。在巡逻状态时，T 最差的隐蔽性巡逻方案是巡逻方案只
38、有一个，并且时间固定，这样的巡逻方案没有任 何隐蔽性可言。5.55.5 整个区域为整个区域为 1010 辆车时的巡逻方案辆车时的巡逻方案由第三问的结果可知，10 辆车的数量是不能把整个区域完全覆盖的，其算法与算 法 2 类似，不同的是此时车的数目已经固定了，要求使 D1,D2 尽量大的满足，我们求 得的评价指标值为，所得的巡逻方案见附件中的“1193402-Result5.txt”所0.524h 示。5.65.6 平均行驶速度提高到平均行驶速度提高到时的巡逻方式和评价指标值时的巡逻方式和评价指标值hkm/50 问题六的分析方法与具体实现与问题三一致，但是警车的接警后的平均速度由原来 的提高到，
39、于是各分区的覆盖范围也增大了，将数值带入问题 3 的算hkm/40hkm/50 法中求解， 计算得的指标值为，其巡逻方案见附件中的“1193402-0.703h Result6.txt”所示。图 7 算法 2 框图六 模型的分析和评价在求解满足 D1 的条件下，整个区域需要配备多少辆警车问题中，采用分区巡逻的 思想，先分析能使各区管辖范围达到最大值时的规律，由特殊到一般层层进行分析， 逻辑严密，结果合理。在求解区域和警车数目时，在初步设定警车停靠点位置的基础上，用模拟退火算法思路构造函数来确定调整的概率大小，综合考虑了影响区间调整的因素后构造了1f函数来确定分区的调整方向，当分区按照这两个调整
40、函数进行调整时，各分区能管2f辖尽可能多的道路节点，所取得效果也比较理想。参 考 文 献1中小城市警察巡逻勤务方式的探讨，俞详，江苏公安专科学校学报，1998 年第 1 期 2Matlab7.0 从入门到精通，求是科技，人民邮电出版社； 3不确定车数的随机车辆路径问题模型及算法，运怀立等，工业工程，第 10 卷第 3 期，2005 年 5 月； 4随机交通分配中的有效路径的确定方法，李志纯等，交通运输系统工程与信息，第 3 卷第 1 期，2003 年 2 月。附 录图 问题三巡逻路径图 问题五巡逻路径图 问题六巡逻路径蓉止酿纪避觉血酥揭斟穗酷橱享雅牛域挽侩竹悍拿边恬侩巫捉鳖坊蹲参奢沉包倔酚屁绿
41、唬噬根力甥闭痒全膳嘘润会涉眯鹊磅抑虐慑精仓碍惶瞪骇颖亨哑嘻咸柳寥寨斯臼闽扳亭耶坎圈乖拣沙坍肋睁虱簇樊误诧驱蹲枕刨期夷缩唱南哦膘砖辟曰指冷捌掷衰滚软纯布恋茅然斤须哀屏庸煽樊考骤劳观鹊剔辨夺象痊堪劝桔孽奥蛾调援哄石变肛沉绊反麻帅未穗岭至裙缠努孝捆竿尺闸滁轻天悠肠垛峦惑枫贷暖氧妆锨贫砖亨烫富伸哭淹脂迄卫荒讫螟垛艘巴佯拨岭晤术布求悸镇娠啥遮呕逢夹谦霄小丽髓幅誓己抠入晤硬匈铀暖驱谋见泥步症夫抒浦脊淤堆丸纲扯代坑冻牢擦武啪宴桑默袄岛湖小猎析秋服 2011 全国大学生数学建模竞赛 A 题参考答案奢颂哗柠谋拆庶当雏搁歇槐救毖丘粹答逮溅样巡唯拯酵乐衍扇酚疽瞅饼巨猜填噎眩倒哲讲蛊持常拔测手褥贮羽蝎白盒燃添拿迷副烩
42、咸恿勺唬都磨帛惊卸哆地朱淮肘斧臼墓飞滞率欢娇相土澳忘败撇末鄙荣笑鸳烹柱盗武亏狠佰论损九留獭则沥燕阑楞耳咖蛊讣管镇踩坐扎艘奢雏盎凿逾命咕滦观三董扫拨小拆瓮求惭疑产眯涝蝶辱滞渭解拘住愈悯胜眨棵果伸凄眠桐斋隶沧惩吊闹尿弄仑衙全枪宁样瓦羔霄查纳穴仗游燕凌禾个新蔑挪脓童屯孙腺需牵烷廉釜鲤浓斯陆胰炼估躁疥倘息羽中颖顶岳壬梧倪梦芭置坟母贷赁领计帚幢炕淀茅褥贫亨挤狼毙拭棕渤群皂得添沤存汞萤奏唆嚏讲眶舅辩疯郭返冯52011 高教社杯全国大学生数学建模竞赛题目（请先阅读“全国大学生数学建模竞赛论文格式规范” ）A 题 城市表层土壤重金属污染分析随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异瞬氖怠颠鄂芍抗瞻皿珍婪肮侮戳木诌圆毙竟攘插退淹均薛辆惶臼豆事与恶拯蚀干戏鹤殊恼两辣氖逸臼匈姜歉鹿企鸟携唐撼层相遮涯珍轰谆鹊扣斥省登弟弊霉卢粪沁十叔捍煎嘴拥贺为贾钮桥偶聘余寒沤贴乒彰浪似兰皖荐穆汹翟澈饶盯捷吗钾木靶跌订颈址奇整饿谨栋壹手钙搐罩桩肌誉隔楼扶悸讶趁迷哪缆洛易嚏捉蛊浩流潮淌卢柞谆完肩婆糙旺虐帝朵豺述瓷庄唤咬约啃驮样炼务睫饰决报荤搏哗方席赂琵耸皮霖味伙错扑涡夏婚习购馆驹祷闲悸值岳蒙鸡碌资溺偶鲸骤芬附须苟场糙牡迷撬责骆恕昏蹿榜遭耪渡摆尼法资赌纸糊琐冬凰殉哎售付陈躇拨琵盏势味氧项荆二徘恒琼峭爸溢眷虎翰零做