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：先将道路离散化成近似均匀分布的节点，相邻两个节点之间的距离约等于一分钟巡逻路程。由警车的数目 ,m将全区划分成 个均匀的分区，
7、从每个分区的中心点出发，找到最近的道路节点，作m为警车的初始位置，由 Floyd 算法算出每辆警车 3 分钟或 2 分钟行驶路程范围内的节点。考虑区域调整的概率大小和方向不同会影响调整结果，本文利用模拟退火算法构造出迁移几率函数，用迁移方向函数决定分区的调整方向。计算能满足 D1 的最小车辆数，即为该区应该配置的最小警车数目，用 MATLAB 计算，得到局部最优解为 13 辆。在选取巡逻显著性指标时，本文考虑了两个方面的指标：一是全面性，即所有警车走过的街道节点数占总街道节点数的比例，用两者之比来评价；二是均匀性，即所有警车经过每个节点数的次数偏离平均经过次数的程度，用方差值来大小评价。问题三
8、：为简化问题，假设所有警车在同一时刻，大致向同一方向巡逻，运动状态分为四种：向左，向右，向上，向下，记录每个时刻，警车经过的节点和能够赶去处理事故的点，最后汇总计算得相应的评价指标。在考虑巡逻规律隐蔽性要求时，文本将巡逻路线进行随机处理，方向是不确定的，采用算法 2 进行计算，得出相应巡逻显著指标，当车辆数减少到 10 辆或巡逻速度变大时，用算法 2 计算巡逻方案和对应的参数，结果见附录所示。本文最后还考虑到 4 个额外因素，给出每个影响因素的解决方案。关键词：模拟退火算法；Floyd 算法；离散化参赛队号 11*02 队员姓名 *佳 *梅 *巍 一 问题的重述110 警车在街道上巡逻，既能够
9、对违法犯罪分子起到震慑作用，降低犯罪率，又能够增加市民的安全感，同时也加快了接处警时间，提高了反应时效，为社会和谐提供了有力的保障。现给出某城市内一区域，其道路数据和地图数据已知，该区域内三个重点部位的坐标分别为：（5112，4806） ， （9126， 4266） ， （7434 ，1332） 。该区域内共有 307 个道路交叉口，为简化问题，相邻两个交叉路口之间的道路近似认为是直线，且所有事发现场均在下图的道路上。该市拟增加一批配备有 GPS 卫星定位系统及先进通讯设备的 110 警车。设 110 警车的平均巡逻速度为 20km/h，接警后的平均行驶速度为 40km/h。警车配置及巡逻方案
10、要尽量满足以下要求：D1. 警车在接警后三分钟内赶到现场的比例不低于 90；而赶到重点部位的时间必须在两分钟之内。D2. 使巡逻效果更显著； D3. 警车巡逻规律应有一定的隐蔽性。现在我们需要解决以下几个问题：一. 若要求满足 D1,该区最少需要配置多少辆警车巡逻？二. 请给出评价巡逻效果显著程度的有关指标。三请给出满足 D1 且尽量满足 D2 条件的警车巡逻方案及其评价指标值。四. 在第三问的基础上，再考虑 D3 条件，给出你们的警车巡逻方案及其评价指标值。五如果该区域仅配置 10 辆警车，应如何制定巡逻方案，使 D1、D2 尽量得到满足？ 六. 若警车接警后的平均行驶速度提高到 50km/
11、h，回答问题三。七. 你们认为还有哪些因素、哪些情况需要考虑？给出你们相应的解决方案。二 问题分析本题为城区道路网络中警车配置及巡逻问题。在进行警车配置时，首先要考虑警车在接警后在规定时间内赶到现场的比例，在此条件下，以车数最少为目标，建模、参赛密码 （由组委会填写）求解；在制定巡逻方案时，要考虑巡逻的效果及隐蔽性问题。问题一只要求满足 D1，求最少的警车配置数，可以认为警车是不动的，在三分钟或两分钟内它能到达的区域就是它的覆盖范围。据此，在满足所有街道的覆盖率不低于 90%的条件下，寻找最优解。问题二要评价巡逻效果，有两个方面需要考虑：一是巡逻的全面性，即经过一段时间后警车走过的街道数占总街
12、道数的比例；二是巡逻的不均匀性，即经过一段时间后警车经过每一条街道的次数相差不大，用方差来衡量。问题三是在满足 D1 的条件上尽量满足问题二所给的指标，并给出评价方案的指标。首先找到一组满足 D1 的各警车位置，然后在和各警车位置相连的点中随机寻找一个点，判断新的点是否满足 D1，如果满足则警车行驶到该点，否则重新寻找，直到满足为止。一段时间后统计所有车走过的点数及每个点被走过的次数，用问题二给出的两个指标进行评价。综合两个指标，可判断此路径的好坏，重复这个过程，直到综合评价指标达到一个满意的值为止。问题四增加了隐蔽性要求，首先给出评价隐蔽性的指标，隐蔽性可用路线的随机性来评价，将它加入到问题
13、三的模型中去进行求解。问题五限制警车数量为 10，要综合考虑 D1、D2，先分配这 10 辆车使道路的覆盖率最高，然后按照问题三的步骤进行求解，其中每一步对 D1 的判断只需使道路的覆盖率尽量高即可。问题六同问题三，只需将车速改为 50km/h 即可。三 模型的假设1. 警车都在路上巡逻，巡警去处理案件的时间不考虑；2. 所有事发现场都在道路上，案件在道路上任一点是等概率发生的；3. 警车初始停靠点是随机的，但尽量让它们分散分布，一辆警车管辖一个分区；4. 假定各个划分区域内，较短时间内，最多会发生一个案件；5. 假设区域内的每条道路都是双行线，不考虑转弯对结果造成的影响；6. 如果重点部位不
14、在道路上的，假设这些重点部位在离它们最近的道路上；7. 图中水域对巡逻方案没有影响。四 符号说明表示警车数目m表示警车初始停靠点到各道路的最短距离d表示整个区域的总道路长度L表示不能在 3 分钟内到达的区域的道路的长度l表示非重点部位的警车在 3 分钟内不能到达现场的比例k表示三分钟内能从接警位置赶到事发现场的最大距离是r表示整个区域总的离散点个数n表示第 区内的节点个数i i表示区内调整函数1f表示模拟退火的时间，表征温度值t表示区间调整函数2f表示全面性指标r表示不均匀性指标e表示综合评价指标h表示第 辆车经过每条道路的次数isi表示整个区域每条道路经过的平均次数五 模型的建立与算法的设计
15、5.1 满足 D1 时，该区所需要配置的最少警车数目和巡逻方案5.1.1 满足 D1 条件时，区域最少警车的规律题目要求警车的配置和巡逻方案满足 D1 要求时，整个区域所需要配置的警车数目最少。由假设可知警车都在道路上，且所有事发现场也都在道路上，但区域内总的道路长度是个定值的；警车在接警后赶到事发现场有时间限制和概率限制：三分钟内赶到普通区域案发现场的比例不低于 90，而赶到重点部位的时间必须控制在两分钟之内。由此可知每辆警车的管辖范围不会很大，于是考虑将整个区域分成若干个分区，每辆警车管辖一个分区域。由上面的分析，求解整个区域的警车数目最少这个问题可转化为求解每一辆警车所能管辖的街道范围尽
16、量的大。于是我们寻找出使每辆警车管辖的范围尽量大的规律。为了简化问题，我们不考虑赶到现场的 90%的几率的限制，仅对警车能在三分钟内赶到事发现场的情况作定性分析，其分析示意图如图 1 所示。警车的初始停靠位置是随机的分布在道路上的任一节点上，我们假设一辆警车停靠在 A 点上。图 1 一辆警车管辖范围分析示意图由于警车的平均巡逻速度为 20km/h，接警后的平均行驶速度为 40km/h，由于距离信息比较容易得到，于是我们将时间限制转化为距离限制，这样便于分析和求解。当警车接警后，在三分钟内能从接警位置赶到事发现场的最大距离是 ，其中r。kmr24063如图 1 所示，我们设警车初始停靠位置在 A
17、 点，A 点是道路 1，2，3，4 的道路交叉口。我们仅以警车在道路 1 巡逻为例来进行分析，警车以 的速度在道路 1hkm/0上 A 到 点之间巡逻， 与初始停靠点 A 的距离为 。由于案件有可能在道路上任 Ax一点发生，当警车巡逻到 A 点时，若案发现场在道路 2，3，4 上发生时，警车以40km/h 的速度向事发现场行驶，警车能在三分钟内从 点赶到现场的最大距离为。如果警车在道路 1 上继续向前行驶，则该警车能在三分钟内赶到现场的距kmx)2(离继续缩小，当警车从初始点向 A 点行驶但没有达到 点时，此时该警车的最大管辖范围比警车到达 点时的最大管辖范围大。为了使警车的管辖范围尽量大，警
18、车的巡逻范围越小越好，当 时，即警车在初始停靠点静止不动时，警车的管辖范围达到0x最大值 。k2图 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 也有个道路交
19、叉口 D， 由于警车巡逻时是在道路上行驶的，行走的路线是分段直线,并不影响路径的长度，所以当警车巡逻到距离初始停靠点 C 点 远处的 D，此时若有案件发生时，该警车要在三x分钟内能赶到现场处理案件，最大行驶距离在 之内，如果警车在道路 1 上继kmx)2(续向前行驶，则该警车能在三分钟内赶到现场的距离继续缩小，当警车没有行驶到 D点时，此时该警车的最大管辖范围比 大，为了使警车的管辖范围尽量大，警)(车的巡逻范围越小越好。当 时，即警车静止不动时，一辆警车的管辖范围能达到0x最大值。以上分析的仅作定性的分析，对于三个重点部位也可以同理分析，所得的结论是一致的，以上的分析没有考虑到 90%的到达
20、几率限制，但在设计算法需要充分考虑。综上所述，当警车静止在初始停靠点时，在三分钟时间限制内，警车能从初始停靠点赶到事发现场的最大距离为 。km25.1.2 将道路离散化由于事发现场是等概率地分布在道路上的，由区域地图可以发现，整个区域中的道路长度不均，为了使计算结果更加精确，可将这些道路离散化。只要选取合适的离散方案，就能使警车在经过道路上的离散的点时就相当于经过了这条道路。这样，不论是求解警车初始停靠点还求解警车赶到事发现场所经过的道路时，所计算得的的结果显然比仅考虑整条道路的叉路口要精确得多。区域中共有 307 个道路交叉口，458 条道路。我们采用线性插值方法对道路进行离散化，以 的速度
21、行走一分钟的距离作为步长，一分钟时间的选择是参照问题三hkm/20的结果要求来设定的，步长 。用线性插值的方法，从道路的一个方kmb31206向进行线性插值，实现将每条道路离散化的目标，考虑到有些道路不是 的整数倍，km31我们就一般情况进行讨论，其分析示意图如图 3 所示。道路 AB 长度为 个 与n长度的和，为了更精确处理 CB 段道路，那么就要考虑在 CB 之间是否要插)31(kmx入一个新的点， 根据 的长度不同，其对应的处理方式也有所不同。x图 3 道路离散化分析示意图引进临界指数 ，选取 大小的准则是使尽量离散化后警车等效的平均巡逻速度和y题目给定的速度（ ）的差值尽量小，经过计算
22、得 时，不再插入新hkm/20 kmy189.0的坐标点时能使整个区域的道路离散效果较好。此时，将 CB 段长度设定为 处理，k31于是离散后的 AB 道路长度会比实际长度短些；当 时，需要在两个点之间Kx.再插入一点，因为这样处理能使整个区域的整体道路的离散化效果比较理想。如图 3所示，在 C 与 B 间再插入新的坐标点，插入的位置在距 C 点 的 D 点处,这样处理后km31所得的道路长度比实际长度长了 。采用这样的方法进行线性插值，我们使用kmx)31(MATLAB 编程实现对整个区域道路的离散，所得的离散结果如图 4 所示，离散后共得到762 个节点，比原始数据多了 455 个节点，离
23、散后的节点数据见附件中的“newpoint.txt”。图 4 整个区域离散结果图采用这种插值方法道路离散后，将直线上的无穷多个点转化有限个点，便于分析问题和实现相应的算法，由图 4 可知，所取得的整体离散效果还是比较理想的。5.1.3 分区域求解警车数目的算法设计考虑到警车配置和巡逻方案需要满足：警车在接警后三分钟内赶到普通部位案发现场的比例不低于 90%，赶到重点部位必须控制在两分钟之内的要求。设计算法的目标就是求解出在满足 D1 情况下，总的警车数目最小，即每个区域都尽可能多地覆盖道路节点。由于警车的初始位置是未知的，我们可设警车初始停靠点在道路上的任一点，即分布在图 4 所示的 762
24、个离散点中的某些点节点上，总体思路是让每两辆车之间尽量分散地分布，一辆警车管辖一个分区，用这些分区覆盖整个区域。于是我们设计算法 1，步骤如下所示：Step1：将整个区域预分配为 个分区，每个分区分配一辆警车，警车的初始停靠位置m设在预分配区中心的道路节点上，若区域的中心不在道路节点上，则将警车放在离中心最近的道路节点上；Step2：统计分区不能覆盖的节点，调整警车的初始停靠点，使分区覆盖尽可能多的道路节点，调整分为区内调整和区间调整方案：（1）区内调整按照模拟退火思想构造的函数，在区间调整调整车辆初始点的位置（后文中有详细说明） ，当分区内节点数较多时，调整的概率小些，分区内节点数较少时，调
25、整的概率大些， （2）当区域中存在未被覆盖的节点或节点群（大于等于三个节点集中在一个范围内）时，将警车初始位置的调整方向为朝着这些未被覆盖的节点按一定的规则（在算法说明中有详细叙述）移动，同时要保证 3 个重点部位能在 2 分钟之内 100%到达；Step3：用 Floyd 算法计算出警车初始停靠点到周边各道路节点的最短距离 ；dStep4：以 个划分区域未覆盖的总的道路长度 与整个区域的道路总长度 的比值ml L来表示警车不能 3 分钟内到达现场的概率；%10LlkStep5：模拟足够多的次数，若 ，将车辆数 减 1，跳转到 Step1；10kmStep6：计算结束后，比较当 时所对应的 值
26、， 当 取得最小值时，记录此时的区域划分方案， 即为最少的警车数。对算法的几点说明：（1）该算法所取的车辆数 是由多到少进行计算的， 初始值设为 20，这个值的m选取是根据区域图估算的。(2)预分区的优点在于使警车的初始位置尽可能均匀地分散分布，警车的初始停靠点在一个分区的中心点附近寻找得到，比起在整个区域随机生成停靠点，计算效率明显得到提高。 预分配之后，需要对整个区域不断地进行调整，调整时需要考虑调整方向和 调整概率。警车调整借鉴的是模拟退火算法的方法，为了使分区内包含道路节点数较多的分区的初始停车点调整的概率小些，而分区内包含道路节点数的少的分区内的初始停车点调整的概率大些，我们构造了一
27、个调整概率函数 ,1f（1）)exp(1tbmnafi（1）式中， 均为常数， 为整个区域车辆数， 为第 分区内覆盖的节点数，ba, in为时间，同时 也能表征模拟退火的温度变化情况：初始温度较高，区域调整速度较tt快，随着时间的增加，温度不断下降，区域调整速度逐渐变慢，这个调整速度变化也是比较符合实际情况的。由式（1）可以得出调整概率函数 ，假设在相同的温度 （时间）的条件下，由1f t于总的车辆数目 是定值，当 时，即第 分区内的节点数大于第 分区的节点数mjini j时，分区 调整的概率大些，分区 的调整概率小些。分析其原因：当分区内包含了较i多的节点个数时，该分区的警车初始停靠位置选取
28、地比较合适了，而当分区内包含的道路节点数较少时，说明警车的初始停靠位置没有选好，需要更大概率的调整，这样的结论也是比较客观的。对于所有分区外未被覆盖的道路节点和很多节点（称之为节点群） ，用来调整警车位置迁移的方向，其分析示意图如图 5 所示。调整方案目标是使未被覆盖的节点数尽量的少。在设计调整方向函数时，需要考虑：（1）节点群内节点的数目；（2）警车距离节点群的位置。优先考虑距离，所以在公式（2）中，用距离的平方来描述调整方向函数。由于某一个区域范围内的未被覆盖节点数，整个区域未被覆盖的节点总数，分区域与未被覆盖的节点或节点群的距离等几个因素会影响到调整的方案，所以要综合考虑这些因素。于是设
29、计了区间调整函数 ，2f式中， 表示第 个分区内未被覆盖的节点数， 表示第 分区域与未被覆盖的节ini ili点或节点群的距离， 表示未被覆盖的节点和节点群个数。 p现在简要分析第 分区按区间调整函数的调整方案，当某两节点群 的节点数目相i ji,等，但是距离不等时，如 ，由区间调整公式可知，该区间向节点群 方向调整。jil当某个分区与两个节点群的距离相等，但节点群的内节点个数不相等，如 时，jin由（4）可知，该分区域会想节点群 方向调整。j注意在整个调整过程中，调整几率控制是否调整，调整方向函数控制调整的方向，寻找在这种调整方案下的最优结果。图 5 调整分区域示意图（3）在 step3 中，使用 Floyd 算法计算出警车初始停靠点到周边各节点的最短距离 ,目d的是当区域内有情况发生时，警车能在要求的时间限制内到达现场。（4）为求出较优的警车停靠点，采用模拟退火算法，算出局部最优的方案。5.1.4 警车的配置和巡逻方案使用 MATLAB 编程实现算法 1 得到，整个区域配备 13 辆警车，这些警车静止在初始停靠点时，能满足 D1 要求。警车的初始停靠位置分别为道路交叉节点6，25，30，37，82，84，110，111，126，214，253，258，278 处。每个警车所管辖的交叉点（原始的交叉节点）如图 6 所示，求解的分区结果见附录所示。
