cumcm2011A附件_数据

承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题.我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出.我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性.如有违反竞赛规则的行为,我们将受到严肃处理.我们参赛选择的题号是(从A/B/C/D中选择一项填写): B我们的参赛报名号为(如果赛区设置报名号的话):所属学校(请填写完整的全名):参赛队员(打印并签名) :1.2.3.指导教师或指导教师组负责人(打印并签名):日期: 2011 年 9 月 11 日赛区评阅编号(由赛区组委会评阅前进行编号):编号专用页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):交巡警服务平台的设置与调度摘要由于警务资源是有限的,所以根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是有关部门面临的一个实际课题.本文着力于通过所给资料,寻找最优化的交巡台设置与调度方案. 按照设置交巡警服务平台的原则和任务,我们首先对问题1用Floyd算法,提出最佳的交巡警服务平台管辖区域划分方案,缩短了出警时间,平衡了工作量,然后采用回溯法,给出了应对突发事件的警力比较合理调度方案；对于问题2,我们将其归结为全局的配置问题,首先用优化后的Floyd 算法对该市现有六城区的交巡警服务平台设置进行改进,其次以时间最短、围堵区域最小为原则,提出了应对重大刑事案件的最佳围堵方案.对于问题1,本文将最短时间问题转化为单向最短路径问题.我们没有运用经典的求最短距的Dijkstra算法,采取时间复杂度更简便的Floyd算法,应用Matlab编程,以出警时间最短为原则,将72个交通节点分配给20个交巡警服务平台；对于出现突发事件,本文采用回溯法,以最节省警力、实现全区封锁联动时间（即封锁路口最长时间）最短为目标,成功的实现了应对突发事件时警力的合理调度；对于某些交巡警服务平台工作量大、出警时间过长等问题,本文利用Mathematica对附表2中的数据进行分析,整理分析A区各节点事故发生率后,利用图论的相关知识,提出应增设4个服务平台,基本实现警力的最优配置.最后,借助于Matlab和Mathematica软件,对附件中所提供的数据进行了筛选,去除异常数据,对残缺数据进行适当补充,并从中随机抽取了3组数据(每组8个采样)对理论结果进行了数据模拟,结果显示,理论结果与数据模拟结果吻合良好.而对于问题2,我们对附件中所提供的A,B,C,D,E,F六城区的数据进行了整合与分析,并做出了直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,并结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,在充分考虑现有警力和财力并确保安全的条件下,科学分析现有平台的数量和具体位置的合理性.数据显示C区和F区的事故发生率较高、交巡警服务平台工作量高于全市平均水平、交巡警服务平台平均每天出警时间过长,针对以上问题我们再次利用均衡二分法,并考虑区域边界处的设点拥挤问题,提出了在C区增设5个交巡平台、F区增设1个交巡平台.对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),确定三分钟后犯罪嫌疑人逃逸的可能覆盖范围,从而利用回溯法的思想采用Matlab编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的穷举算法,快速地形成围堵区域,并实现了围堵区域最小的目的.实现了资源调配问题的优化决策.考虑到该城市未来发展规划,只需对本文所建模型进行适当改进即可,在此不进行详细解答.关键词最短路径 Floyd算法回溯法穷举法优化决策目录交巡警服务平台的设置与调度 (1)摘要 (1)1.问题重述 (1)2.问题分析 (1)2.1对于问题一的分析 (1)2.2对问题二的分析 (1)3.模型假设 (2)4.定义与符号说明 (2)5.模型的建立与求解 (2)5.1 问题一的模型 (2)5.1.1 模型建立 (2)5.1.2 模型求解 (3)5.2 问题二的模型 (8)5.2.1 模型建立 (8)5.2.2 模型求解 (9)7.模型的评价与推广 (11)8. 附件 (12)附件1:用Floyd算法分配个服务平台管辖区域 (12)附件2:邻接矩阵的matlab实现程序 (23)附件3:围堵方案的java实现程序 (30)附件4:全区的交巡警平台有效覆盖范围(有效代表三分钟内可以到达) (31)附件5:用Mathmatica求数据均值与方差 (31)附件6:输入任意两点的坐标,输出两点间距离 (32)附件7:A区各线路距离 (32)1.问题重述“有困难找警察”,是家喻户晓的一句流行语.警察肩负着刑事执法、治安管理、交通管理、服务群众四大职责.为了更有效地贯彻实施这些职能,需要在市区的一些交通要道、人员密集区和重要部位设置交巡警服务平台.每个交巡警服务平台的职能和警力配备基本相同.由于警务资源的有限性,根据城市的实际情况与需求,合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题.本文着力于寻找最优化的设置与调度方案.问题1要求合理分配交巡警服务平台的管辖范围,使其在所管辖的范围内出现突发事件时,尽量能在3分钟内有交巡警(警车的时速为60km/h)到达事发地；对于重大突发事件,给出该区交巡警服务平台警力合理的调度方案,尽快封锁道路；拟在该区内再增加2至5个平台,以减少出警时间、平均工作量,确定需要增加平台的具体个数和位置.问题2要求分析研究该市现有交巡警服务平台设置方案的合理性并给出解决方案；如果该市地点P(第32个节点)处发生了重大刑事案件,在案发3分钟后接到报警,犯罪嫌疑人已驾车逃跑.为了快速搜捕嫌疑犯,给出调度全市交巡警服务平台警力资源的最佳围堵方案.2.问题分析本题所要解决的是A区以及全市的安巡警服务平台设置与调度问题,根据现实生活状况,我们首先要考虑的是警力资源的限制,即要使得所布置的警力尽可能的少.其次是在交巡台数量最少的情况下,力求警员到达现场的时间在3分钟以内,解决突发状况.2.1对于问题一的分析该市中心城区A的交通网络有92个节点和20个交巡警服务平台,要求当突发事件发生时,尽量能在3分钟内有交巡警到达事发地,已知警车的时速为V=60km/h,我们将最短时间转化为最短路问题,应用Floyd算法,求解出A区距离每一节点最近的交巡台,即将该节点分配给该交巡台.对于重大突发事件,要实现对进出该区的13条交通要道进行快速封锁,即需调度交巡台尽快到达13个节点,重复Floyd算法,找出最近交巡台,即可找出调配方案.但需注意的是,有的出入口本来就有交巡台,但为了达最优化,需进行重新分配,故应用回溯法,找到调度方案.现有交巡台工作量不均衡和有些地方出警时间过长,统计A区各个交巡台案发率,计算均值与方差,在案发率较高地带增设交巡台,平衡工作量,尽量缩短出警时间.2.2对问题二的分析对于问题二,是对问题一的进一步改进与推广,在遵循警情主导警务原则,快速出警原则与方便与安全原则,结合辖区地域特征、人口分布和治安状况等实际情况,充分考虑现有警力和财力并确保安全的条件下,设置交巡平台,重复上一问的做法,评估交巡平台的合理性.对于改进方案,应考虑城区内部工作量,城区之间的联系以及城市边界的警力调度.对于突发状况的围堵方案,应在最短时间内对可能逃逸区域进行合围,最小范围内缩小包围圈.3.模型假设1.假设题中所给数据均真实可靠.2.出警时道路恒畅通(无交通事故、交通堵塞等发生),警车行驶正常,警车及肇事车辆行驶时均以60km/h匀速行驶,转弯处不需要花费时间.3.事故均发生在路口节点,两节点连线上认为没有事故发生.4.每条线路行驶都是双向的.5.考虑肇事车辆在P点向各个方向逃逸的概率相等.6.在整个行驶中,车辆只在主要干道行驶.7.发生事故时,忽略反应调度时间.4.定义与符号说明m任意两个标志点i与j之间的距离ijm标志点间的距离组成的距离矩阵n标志点的邻接矩阵n邻接矩阵的元素ijD相邻标志点间的距离矩阵D相邻标志点i与j间的距离ijW标志点的权值矩阵d标志点间的最短距离矩阵d标志点i与j之间的最短距离ijv肇事车辆逃逸速度15.模型的建立与求解5.1 问题一的模型5.1.1 模型建立此问是关于最短路径的模型分析及MATLAB的实现A区道路状况及交巡台的设置如图1所示.本文应用Floyd算法,通过构造距离矩阵,依次找出距离每一节点最近的交巡台,使得有事故发生时,交巡警在最短时间内到达事故现场,以此为依据分配管辖区域.如果道路不通时,认为两端节点的距离为无穷.图1 A 区各节点及服务平台示意图当有重大突发事件时,要对进出该区的13条交通要道进行快速封锁,固定13个出入口,应用回溯法,找到距离节点最近的交巡平台.封锁时间决定于最后到达节点的时间,由于一个平台的警力最多封锁一个路口,至少需调动13个平台的警力.为达到工作量的均衡和出警时间尽可能的短,需进行优化决策.考虑每一节点案发率的不同,在A 区增设2到5个平台,使得每一平台的工作量均衡,平均出警时间大体相同.5.1.2 模型求解首先我们可以根据题中所给的各个标志点的坐标,用matlab 计算出任意两点之间的直线距离,得到92*92的距离矩阵:1111n n nn m m m m m ⎛⎫ ⎪= ⎪ ⎪⎝⎭根据题中的分布图,我们可以得到各标志点的邻接矩1111n n nn n n n n n ⎛⎫ ⎪= ⎪ ⎪⎝⎭即如果两个点相邻,则邻接矩阵中相对应的元素的值为1,否则为0；例如:3和44这两个点相邻,那么3,4444,3==n n .根据Floyd 算法,我们是要求出任意两节点之间的距离,所以我们需要得到相邻两个结点的直线距离.我们可以利用距离矩阵的元素ij m 与ij n 的点乘积得到相邻标志点间的距离矩阵:1111.*n n nn D D D m n D D ⎛⎫ ⎪== ⎪ ⎪⎝⎭对于D 中不相邻点间距离0改为无穷大(Inf)从而得到节点与节点间的权值矩阵:1111n n nn W W W W W ⎛⎫ ⎪= ⎪ ⎪⎝⎭ 即如果15和10之间不相邻,也即不能直接到达,那么D 中的010,15=D 和015,10=D 都将变成10,15W 和15,10W 等于无穷大(Inf),否则则等于D 中相应元素的数据.运用Floyd 算法求出任意两点间最短距离,得到最短距离矩阵d :1111n n nn d d d d d ⎛⎫ ⎪= ⎪ ⎪⎝⎭由Floyd 算法,运行MATLAB 程序,可统计出距离每一节点最近的交巡台的位置,MATLAB 运行结果如表1所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.交巡台—节点 距离 交巡台—节点 距离13—21 27.0831 4—57 18.681513—22 9.0554 6—58 23.841413—23 5.0000 6—59 16.031213—24 23.8537 4—60 17.924012—25 17.8885 4—(61) 52.105511—26 9.0000 4—62 3.500011—27 16.4330 4—63 10.308715—(28) 47.5184 4—64 9.363215—(29) 57.0052 3—65 15.23987—30 5.8310 3—66 18.40129—31 20.5572 1—67 14.91587—32 11.4018 1—75—68 10.79278—33 8.2765 1—69 5.00009—34 5.0249 2—70 8.60239—35 4.2426 1—74—71 11.265016—36 6.0828 2—72 16.403116—37 11.1818 18—73 19.723116—(38) 34.0588 1—74 6.26502—(39) 36.8219 1—75 6.26502—40 19.1442 1—76 9.800517—418.5000 19—779.848917—429.84891—78 6.40312—438.000019—79 4.47212—449.846818—808.06239—4510.9508 18—81 6.70828—469.300518—8210.79357—4712.806218—83 5.38527—4812.902120—8411.75225—49 5.000020—85 4.47215—508.485320—86 3.60505—5112.893220—8714.65115—5217.194420—8812.94645—5311.708220—8914.75223—5422.708918—9019.52563—5512.659020—9116.00605—5621.437020—(92)36.0060表1 该市A区指定节点到交巡警服务平台最短距离由上表可初步确定A区20个交巡台的管辖范围,如表2所示.带括号的节点为发生事故时任意交巡台都不能在三分钟内赶到节点.交巡台序号辖区内节点辖区内案发率交巡台序号辖区内节点辖区内案发率1 67 68 69 71 74 75 76 78 9.42 40 43 44 70 72 39 9.7 3 54 55 65 66 5.6 4 57 60 62 63 64 6.6 5 49 50 51 52 53 56 7.7 6 58 59 4.5 7 30 32 47 48 61 9 8 33 46 5 9 31 34 35 45 8.2 10 1.6 11 26 27 4.6 12 25 4 13 21 22 23 24 8.5 14 2.5 15 (28) (29) 4.8 16 36 37 (38) 5 17 41 42 5.3 18 73 80 81 82 83 719 77 79 3.4 20 84 85 86 87 88 89 9190 (92)11.5表2 该市A区交巡警服务平台所管辖交叉路口清单图2 A区各交巡台管辖区域示意图需要说明的是,同一条路整体归一个交巡台管理.当有重大突发事件时,固定13个进出A区的节点,运用回溯法,结合上表,找到距离节点最近的交巡台,以此来达到总体时间的最短,我们一共可以得到四个方案,在这个过程中可以发现,有些交巡台要避免去最近的节点封锁而去较远的节点,以此来节省警力.具体封锁方案如表3、表4所示.最短调度时间均为8.0155.方案一:交巡台过程出入口节点号2 路径40→3938 时间 3.9822min4 路径62 时间0.3500min6 路径47→4830 时间 3.1829min7 路径3029 时间8.0154min8 路径4748 时间 3.0995min9 路径35→3616 时间 1.5083min10 路径26→2712 时间7.5863min11 路径22 时间 3.2696min12 路径25 24时间 3.5916min13 路径23 时间0.5000min14 路径21 时间 3.2649min15 路径28 时间 4.7518min16 路径14 时间 6.7417min表3 A区突发事件封锁方案一方案二方案三方案四路口标号平台号路口标号平台号路口标号平台号1213121312101416142314161661691662114211121142210221022122311231423132412241224112815281528152972972973083063093819381738148548948862206220622表4 A区突发事件封锁方案二、三、四在对交巡台均衡工作量,加快出警时间方面,综合各节点的案发率、交巡台到其辖区内任一节点的路程进行综合评估,做出优化决策.在案发率较高地带增设交巡台,以缓解周围交巡台的工作压力,为达均衡工作量的目的,将32号节点从7号交巡台归到8号交巡台,44号节点从2号交巡台归到3号交巡台,39号节点从2号交巡台归到16号交巡台,47号节点从6号交巡台归到7号交巡台,61号节点从7号交巡台归到4号交巡台.这样,A区每交巡台平均每天处理案件数从6.1950件、方差6.8289降到每天处理5.1917件、方差2.2182,极大的协调了工作量.对于个别节点的重新划分,会增加出警时间,但在总体上平均每天的出警时间大大缩短了.综上考虑,共增设4个交巡台,重新分配的结果如表4所示.交巡台序号负责区域内的节点管辖区域内的案发率平均每天出警时间1 69 71 74 75 78 6.6 3.50692 40 43 70 72 7.2 6.67363 54 55 44 5.2 4.93144 57 60 61 62 63 6.4 7.70155 49 50 51 52 53 56 7.7 5.94566 58 59 47 6.1 6.39497 30 48 5.9 3.52588 33 46 5.0 2.27489 35 45 4.9 2.127110 1.611 26 27 4.6 2.394612 25 4.0 2.862213 23 24 5.7 3.823914 2.515 (28) (29) 4.8 14.158016 36 37 (38) (39) 6.4 10.056217 41 42 5.3 2.568918 73 80 81 83 5.9 3.843819 77 79 3.4 1.145720 85 86 87 (92) 6.4 5.535421 22 2.8 2.523931 32 34 4.9 4.396266 64 65 67 68 76 5.1 2.665590 82 84 88 89 91 6.2 3.2171表5 优化后的A区交巡台管辖区域示意图图3 A区增设平台示意图图中方块所示节点即为增设平台处.5.2 问题二的模型5.2.1 模型建立对于问题2,对附件中所提供的A,B,C,D,E,F六城区的数据进行整合,做出直观的图表.遵循警情主导警务原则、快速出警原则、方便与安全原则,结合辖区地域特征、人口分布、交通状况、治安状况和未来城市发展规划等实际情况,充分考虑现有警力和财力并确保安全,科学分析现有平台的数量和具体位置的合理性.对于该市地点P(第32个节点)处发生了重大刑事案件的围堵问题,本文将其归结为资源调配问题.本文合理假设了犯罪嫌疑人的车行驶速度(分三种情况考虑:等于警车速度,警车速度的二倍,警车速度的一半),并确定三分钟后犯罪嫌疑人的车行驶的最远距离,从而利用回溯法的思想采用Matlab 编程确定犯罪嫌疑人的车的所有可能位置.以时间最短、围堵区域最小为原则,采用改进的双层Floyd 算法,快速地形成围堵区域,并使围堵区域尽可能的小.5.2.2 模型求解全市整体状况如表5所示,数据显示C 区和F 区的事故发生率较高、交巡警服务平台工作量高于全市平均水平且交巡警服务平台平均每天出警时间过长,针对以上问题本文再次利用问题1的Floyd 算法,并考虑区域边界处的设点拥挤问题,本文提出了在C 区增加5个服务平台、在F 区增加1个服务平台. 全市六个城区 城区面积城区人口 平台数 平均人口 全区案发率 各区平台案发率均值A 22 60 20 2.727 124.5 6.625B 103 21 8 0.204 66.4 8.3C 221 49 17 0.223 187.2 11.012D 383 73 9 0.191 67.8 7.533E 432 76 15 0.176 119.4 7.96F 274 53 11 0.193 109.2 9.927 均值 53.3333 13.3333 0.619 112.4167 8.5595表6 全市整体状况图4 全市增设交巡台位置示意图(方块所示区域)对于P 点发生重大刑事案件,动用全市警力进行围堵,我们希望使得包围圈尽可能的小,由于犯罪嫌疑人的车速度未知,我们分以下三种情况进行考虑:1)当犯罪嫌疑人的车速与警车速度同,即h km v /601=.运用穷举法,对肇事车辆可能的逃逸路线进行分析,以3分钟路程为半径,找到肇事车辆逃逸的覆盖范围,如图5所示,其中实线表示可能路径,在此范围内有8、9、10、15号共4个交巡平台,保证这4个平台警力不动,组成第一范围包围圈.图5 肇事车逃逸3分钟内覆盖区域示意图进一步分析可能的逃逸路线,调度16号交巡台到36号节点,2号交巡台到3号节点,3号交巡台到55号节点,6号交巡台到47号节点,组成第二组半包围,保证对A 区的封锁.若肇事车辆经36号节点逃往16号节点,则会与16号交巡台在途中相遇.对于从32号节点经7号节点逃逸到30号和47号节点,存在从A 区逃往其他城区的可能,需调动其他城区交巡台的支援.将C 区119号交巡台调度到237号节点,将D 区320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点,至此,在全市范围内实现全面封锁.2)当犯罪嫌疑人的车速比警车车速小,即h km v /601<,我们令h km v /301=方法同1),寻找分钟逃逸范围内所覆盖的全部交巡台,经过整合分析,保持7号、8号、9号、15号共4个交巡台原地封锁,10号交巡台到34号节点封锁,6号交巡台到47号节点封锁,16号交巡台到36号节点封锁,3号交巡台经55号节点到46号节点进行封锁,2号交巡台经3号节点到45号节点封锁,在此过程中,10号、2号和3号交巡台会在途中与肇事车辆相遇.3)当犯罪嫌疑人的车速比警车车速大,即h km v /601>,我们令h km v /1201=由于肇事车辆逃逸速度较快,可能会逃逸到C 区和F 区,故需调动C 区和F 区警力进行围堵.A 区将20号交巡台调到62号节点,16号交巡台调到36号节点,2号交巡台经40号到39号节点,17号交巡台调到41号节点,15号、10号、4号、3号、5号、7号、8号、9号交巡台原地封锁,其余交巡台向其邻近的路口节点进行增援.经过分析,肇事车辆可能由28号、48号、30号进入C 区及A 、D 两区的交汇地带,或由16号节点逃逸到F 区,在此,对C 区、D 区、F 区交巡台进行如下调配,实现全市封锁:C 区:240号交巡台调到239号节点,170号交巡台调到225节点,167号交巡台调到259节点.D 区:320号交巡台调度到371号节点,321号交巡台经368号、369号节点到370号节点.F 区:477号交巡台调度到501号节点,518号交巡台调到521号节点,478号节点调到527号节点,484号节点到571号节点.7.模型的评价与推广本文避免了时间复杂度较复杂的Dijkstra 算法,选用Floyd 算法,在求最短路径上提高了效率,代码编写简单.模型的建立思路清晰,遵循可操作性、科学性、可比性原则,该模型建立出了在较理想状态下交巡警平台的最优设置,减少出警时间,均衡工作量,提高工作效率,在遇突发事件时,可尽快实现道路封锁,给生活中交巡警平台的设立予参考,具有一定的实际应用价值,也可以应用于其他适用区域.模型的运算由矩阵、向量的运算组成,易于用数学软件求解和验证.本模型较好的解决了交巡警平台的最优选址问题,当事故发生时,交巡警可以第一时间到达事发地点,有效的改善了交巡警在执行任务中的效率,在经济迅猛发展的今天,城市加速扩张,人口迅速增长,交巡警平台的设置是平安城市的最好保障.该模型也可运用到其他最优选址问题中去,比如关于消防救援工作最优路径问题、重大生产安全事故应急救援问题、公共交通的最优路径问题等. 同时也可利用该模型算法拓展模型在其他领域的适用范围.该模型也有一定的局限性,如现实中不能时刻都保证道路的畅通性.既不能保证出警的时间总是维持在3分钟之内.忽略了实际地形对于车速的影响以及实际生活中存在的不定因素.参考文献[1] [徐孝凯,王凤禄],《数据结构简明教程》第二版,北京:清华大学出版社,2005年4月1日[2] [李建中,骆吉洲],《华章数学译丛》第二版,北京:机械工业出版社,2002年6月[3] [陈庆华等],《组合最优化技术及其应用》第1版,北京:国防科技大学出版社,1989年8月[4] [W.T.Tutte],《Graph Theory 》,英国:Cambridge University Press,2001年3月1日8.附件附件1:用Floyd算法分配个服务平台管辖区域[area_a _x,area_a _y]=find (location_all _daolu<=92);road_index _a=[area_a _x,area_a _y];road_index _aa1=find (road_index _a (:,2)==1);a2=find (road_index _a (:,2)==2);A=road_index _a (a1,1);B=road_index _a (a2,1);[c,ia,ib] = intersect (A,B); csize (c)for i=1:140tt=c (i);uu=location_all _daolu (tt,:);uu1=uu (1);uu2=uu (2);vv1=location_a _zuobiao (uu1,:);vv2=location_a _zuobiao (uu2,:);ww1=[vv1 (1),vv2 (1)];ww2=[vv1 (2),vv2 (2)];line (ww1,ww2)endfor i=1:140tt=c (i);uu=location_all _daolu (tt,:);uu1=uu (1);uu2=uu (2);vv1=location_a _zuobiao (uu1,:);vv2=location_a _zuobiao (uu2,:);ww1=[vv1 (1),vv2 (1)];ww2=[vv1 (2),vv2 (2)];line (ww1,ww2,'k')end% ??? Error using ==> line% String argument is an unknown option.for i=1:140tt=c (i);uu=location_all _daolu (tt,:);uu1=uu (1);uu2=uu (2);vv1=location_a _zuobiao (uu1,:);vv2=location_a _zuobiao (uu2,:);ww1=[vv1 (1),vv2 (1)];ww2=[vv1 (2),vv2 (2)];line (ww1,ww2,'Color',[.8 .8 .8])endsave data_b _problem% A区节点间的邻接矩阵load data_b _problem;matric_lingjie=zeros (92,92);[xx,yy]=find (matric_lingjie==0);matric_lingjie (xx,yy)=inf;for i=1:92matric_lingjie (i,i)=0;endfor i=1:140tt=c (i);uu=location_all _daolu (tt,:);uu1=uu (1);uu2=uu (2); % 端点序号vv1=location_a _zuobiao (uu1,:); % 第一个端点坐标vv2=location_a _zuobiao (uu2,:); % 第二个端点坐标% 计算端点间距离distance=sqrt ((vv1 (1)-vv2 (1))^2+(vv1 (2)-vv2 (2))^2);matric_lingjie (uu1,uu2)=distance;matric_lingjie (uu2,uu1)=distance; % 赋值给邻接矩阵end[D,R]=floyd (matric_lingjie);matric_fenkuai=D (1:20,:);for i=1:20for j=1:92if matric_fenkuai (i,j)>30matric_fenkuai (i,j)=0;endendendti=zeros (1,92);ti (1)=text (location_a _zuobiao (1,1),location_a _zuobiao (1,2)+1.5,'1'); ti (2)=text (location_a _zuobiao (2,1),location_a _zuobiao (2,2)+1.5,'2'); ti (3)=text (location_a _zuobiao (3,1),location_a _zuobiao (3,2)+1.5,'3'); ti (4)=text (location_a _zuobiao (4,1),location_a _zuobiao (4,2)+1.5,'4');ti (6)=text (location_a _zuobiao (6,1),location_a _zuobiao (6,2)+1.5,'6'); ti (7)=text (location_a _zuobiao (7,1),location_a _zuobiao (7,2)+1.5,'7'); ti (8)=text (location_a _zuobiao (8,1),location_a _zuobiao (8,2)+1.5,'8'); ti (9)=text (location_a _zuobiao (9,1),location_a _zuobiao (9,2)+1.5,'9'); ti (10)=text (location_a _zuobiao (10,1),location_a _zuobiao (10,2)+1.5,'10'); ti (11)=text (location_a _zuobiao (11,1),location_a _zuobiao (11,2)+1.5,'11'); ti (12)=text (location_a _zuobiao (12,1),location_a _zuobiao (12,2)+1.5,'12'); ti (13)=text (location_a _zuobiao (13,1),location_a _zuobiao (13,2)+1.5,'13'); ti (14)=text (location_a _zuobiao (14,1),location_a _zuobiao (14,2)+1.5,'14'); ti (15)=text (location_a _zuobiao (15,1),location_a _zuobiao (15,2)+1.5,'15'); ti (16)=text (location_a _zuobiao (16,1),location_a _zuobiao (16,2)+1.5,'16'); ti (17)=text (location_a _zuobiao (17,1),location_a _zuobiao (17,2)+1.5,'17'); ti (18)=text (location_a _zuobiao (18,1),location_a _zuobiao (18,2)+1.5,'18'); ti (19)=text (location_a _zuobiao (19,1),location_a _zuobiao (19,2)+1.5,'19'); ti (20)=text (location_a _zuobiao (20,1),location_a _zuobiao (20,2)+1.5,'20'); ti (21)=text (location_a _zuobiao (21,1),location_a _zuobiao (21,2)+1.5,'21'); ti (22)=text (location_a _zuobiao (22,1),location_a _zuobiao (22,2)+1.5,'22'); ti (23)=text (location_a _zuobiao (23,1),location_a _zuobiao (23,2)+1.5,'23'); ti (24)=text (location_a _zuobiao (24,1),location_a _zuobiao (24,2)+1.5,'24'); ti (25)=text (location_a _zuobiao (25,1),location_a _zuobiao (25,2)+1.5,'25'); ti (26)=text (location_a _zuobiao (26,1),location_a _zuobiao (26,2)+1.5,'26'); ti (27)=text (location_a _zuobiao (27,1),location_a _zuobiao (27,2)+1.5,'27'); ti (28)=text (location_a _zuobiao (28,1),location_a _zuobiao (28,2)+1.5,'28'); ti (29)=text (location_a _zuobiao (29,1),location_a _zuobiao (29,2)+1.5,'29'); ti (30)=text (location_a _zuobiao (30,1),location_a _zuobiao (30,2)+1.5,'30'); ti (31)=text (location_a _zuobiao (31,1),location_a _zuobiao (31,2)+1.5,'31'); ti (32)=text (location_a _zuobiao (32,1),location_a _zuobiao (32,2)+1.5,'32'); ti (33)=text (location_a _zuobiao (33,1),location_a _zuobiao (33,2)+1.5,'33'); ti (34)=text (location_a _zuobiao (34,1),location_a _zuobiao (34,2)+1.5,'34'); ti (35)=text (location_a _zuobiao (35,1),location_a _zuobiao (35,2)+1.5,'35'); ti (36)=text (location_a _zuobiao (36,1),location_a _zuobiao (36,2)+1.5,'36'); ti (37)=text (location_a _zuobiao (37,1),location_a _zuobiao (37,2)+1.5,'37'); ti (38)=text (location_a _zuobiao (38,1),location_a _zuobiao (38,2)+1.5,'38'); ti (39)=text (location_a _zuobiao (39,1),location_a _zuobiao (39,2)+1.5,'39'); ti (40)=text (location_a _zuobiao (40,1),location_a _zuobiao (40,2)+1.5,'40'); ti (41)=text (location_a _zuobiao (41,1),location_a _zuobiao (41,2)+1.5,'41'); ti (42)=text (location_a _zuobiao (42,1),location_a _zuobiao (42,2)+1.5,'42'); ti (43)=text (location_a _zuobiao (43,1),location_a _zuobiao (43,2)+1.5,'43'); ti (44)=text (location_a _zuobiao (44,1),location_a _zuobiao (44,2)+1.5,'44'); ti (45)=text (location_a _zuobiao (45,1),location_a _zuobiao (45,2)+1.5,'45'); ti (46)=text (location_a _zuobiao (46,1),location_a _zuobiao (46,2)+1.5,'46'); ti (47)=text (location_a _zuobiao (47,1),location_a _zuobiao (47,2)+1.5,'47'); ti (48)=text (location_a _zuobiao (48,1),location_a _zuobiao (48,2)+1.5,'48');ti (50)=text (location_a _zuobiao (50,1),location_a _zuobiao (50,2)+1.5,'50'); ti (51)=text (location_a _zuobiao (51,1),location_a _zuobiao (51,2)+1.5,'51'); ti (52)=text (location_a _zuobiao (52,1),location_a _zuobiao (52,2)+1.5,'52'); ti (53)=text (location_a _zuobiao (53,1),location_a _zuobiao (53,2)+1.5,'53'); ti (54)=text (location_a _zuobiao (54,1),location_a _zuobiao (54,2)+1.5,'54'); ti (55)=text (location_a _zuobiao (55,1),location_a _zuobiao (55,2)+1.5,'55'); ti (56)=text (location_a _zuobiao (56,1),location_a _zuobiao (56,2)+1.5,'56'); ti (57)=text (location_a _zuobiao (57,1),location_a _zuobiao (57,2)+1.5,'57'); ti (58)=text (location_a _zuobiao (58,1),location_a _zuobiao (58,2)+1.5,'58'); ti (59)=text (location_a _zuobiao (59,1),location_a _zuobiao (59,2)+1.5,'59'); ti (60)=text (location_a _zuobiao (60,1),location_a _zuobiao (60,2)+1.5,'60'); ti (61)=text (location_a _zuobiao (61,1),location_a _zuobiao (61,2)+1.5,'61'); ti (62)=text (location_a _zuobiao (62,1),location_a _zuobiao (62,2)+1.5,'62'); ti (63)=text (location_a _zuobiao (63,1),location_a _zuobiao (63,2)+1.5,'63'); ti (64)=text (location_a _zuobiao (64,1),location_a _zuobiao (64,2)+1.5,'64'); ti (65)=text (location_a _zuobiao (65,1),location_a _zuobiao (65,2)+1.5,'65'); ti (66)=text (location_a _zuobiao (66,1),location_a _zuobiao (66,2)+1.5,'66'); ti (67)=text (location_a _zuobiao (67,1),location_a _zuobiao (67,2)+1.5,'67'); ti (68)=text (location_a _zuobiao (68,1),location_a _zuobiao (68,2)+1.5,'68'); ti (69)=text (location_a _zuobiao (69,1),location_a _zuobiao (69,2)+1.5,'69'); ti (70)=text (location_a _zuobiao (70,1),location_a _zuobiao (70,2)+1.5,'70'); ti (71)=text (location_a _zuobiao (71,1),location_a _zuobiao (71,2)+1.5,'71'); ti (72)=text (location_a _zuobiao (72,1),location_a _zuobiao (72,2)+1.5,'72'); ti (73)=text (location_a _zuobiao (73,1),location_a _zuobiao (73,2)+1.5,'73'); ti (74)=text (location_a _zuobiao (74,1),location_a _zuobiao (74,2)+1.5,'74'); ti (75)=text (location_a _zuobiao (75,1),location_a _zuobiao (75,2)+1.5,'75'); ti (76)=text (location_a _zuobiao (76,1),location_a _zuobiao (76,2)+1.5,'76'); ti (77)=text (location_a _zuobiao (77,1),location_a _zuobiao (77,2)+1.5,'77'); ti (78)=text (location_a _zuobiao (78,1),location_a _zuobiao (78,2)+1.5,'78'); ti (79)=text (location_a _zuobiao (79,1),location_a _zuobiao (79,2)+1.5,'79'); ti (80)=text (location_a _zuobiao (80,1),location_a _zuobiao (80,2)+1.5,'80'); ti (81)=text (location_a _zuobiao (81,1),location_a _zuobiao (81,2)+1.5,'81'); ti (82)=text (location_a _zuobiao (82,1),location_a _zuobiao (82,2)+1.5,'82'); ti (83)=text (location_a _zuobiao (83,1),location_a _zuobiao (83,2)+1.5,'83'); ti (84)=text (location_a _zuobiao (84,1),location_a _zuobiao (84,2)+1.5,'84'); ti (85)=text (location_a _zuobiao (85,1),location_a _zuobiao (85,2)+1.5,'85'); ti (86)=text (location_a _zuobiao (86,1),location_a _zuobiao (86,2)+1.5,'86'); ti (87)=text (location_a _zuobiao (87,1),location_a _zuobiao (87,2)+1.5,'87'); ti (88)=text (location_a _zuobiao (88,1),location_a _zuobiao (88,2)+1.5,'88'); ti (89)=text (location_a _zuobiao (89,1),location_a _zuobiao (89,2)+1.5,'89'); ti (90)=text (location_a _zuobiao (90,1),location_a _zuobiao (90,2)+1.5,'90'); ti (91)=text (location_a _zuobiao (91,1),location_a _zuobiao (91,2)+1.5,'91'); ti (92)=text (location_a _zuobiao (92,1),location_a _zuobiao (92,2)+1.5,'92');fenpei=[1 1 1 671 681 691 711 731 741 751 761 782 22 392 402 432 442 702 723 33 543 553 653 664 44 574 604 624 634 645 55 495 505 515 525 535 565 585 596 67 77 307 327 477 618 8 8 338 469 9 9 31 9 34 9 359 4510 1011 11 11 2611 2712 1212 2513 13 13 21 13 22 13 2313 2414 1415 15 15 2815 2916 16 16 36 16 3716 3817 17 17 4117 4218 18 18 80 18 81 18 8218 8319 19 19 7719 7920 20 20 84 20 85 20 8620 8820 8920 9020 9120 92];for i=1:92if fenpei (i,1)==1plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b');hold on;elseif fenpei (i,1)==2plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g');hold on;elseif fenpei (i,1)==3plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r');hold on;elseif fenpei (i,1)==4plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c');hold on;elseif fenpei (i,1)==5plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m');hold on;elseif fenpei (i,1)==6plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'y');hold on;elseif fenpei (i,1)==7plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'k');hold on;elseif fenpei (i,1)==8plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b+');hold on;elseif fenpei (i,1)==9plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bo');hold on;elseif fenpei (i,1)==10plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'bs');hold on;elseif fenpei (i,1)==11plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'g+');hold on;elseif fenpei (i,1)==12plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'go');hold on;elseif fenpei (i,1)==13plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'gs');hold on;elseif fenpei (i,1)==14plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'r+');hold on;elseif fenpei (i,1)==15plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'ro');hold on;elseif fenpei (i,1)==16plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'rs');hold on;elseif fenpei (i,1)==17plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'c+');hold on;elseif fenpei (i,1)==18plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'co');hold on;elseif fenpei (i,1)==19plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'cs');hold on;elseplot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'m+');hold on;endendfor i=1:92plot (location_a _zuobiao (i,1),location_a _zuobiao (i,2),'b');hold on;endaxis ([190,460,250,410]);for i=1:140tt=c (i);uu=location_all _daolu (tt,:);uu1=uu (1);uu2=uu (2);vv1=location_a _zuobiao (uu1,:);vv2=location_a _zuobiao (uu2,:);ww1=[vv1 (1),vv2 (1)];ww2=[vv1 (2),vv2 (2)];line (ww1,ww2,'Color',[.5 .5 .5])endti=zeros (1,92);ti (2)=text (location_a _zuobiao (2,1),location_a _zuobiao (2,2)+1.5,'2'); ti (3)=text (location_a _zuobiao (3,1),location_a _zuobiao (3,2)+1.5,'3'); ti (4)=text (location_a _zuobiao (4,1),location_a _zuobiao (4,2)+1.5,'4'); ti (5)=text (location_a _zuobiao (5,1),location_a _zuobiao (5,2)+1.5,'5'); ti (6)=text (location_a _zuobiao (6,1),location_a _zuobiao (6,2)+1.5,'6'); ti (7)=text (location_a _zuobiao (7,1),location_a _zuobiao (7,2)+1.5,'7'); ti (8)=text (location_a _zuobiao (8,1),location_a _zuobiao (8,2)+1.5,'8'); ti (9)=text (location_a _zuobiao (9,1),location_a _zuobiao (9,2)+1.5,'9'); ti (10)=text (location_a _zuobiao (10,1),location_a _zuobiao (10,2)+1.5,'10'); ti (11)=text (location_a _zuobiao (11,1),location_a _zuobiao (11,2)+1.5,'11'); ti (12)=text (location_a _zuobiao (12,1),location_a _zuobiao (12,2)+1.5,'12'); ti (13)=text (location_a _zuobiao (13,1),location_a _zuobiao (13,2)+1.5,'13'); ti (14)=text (location_a _zuobiao (14,1),location_a _zuobiao (14,2)+1.5,'14'); ti (15)=text (location_a _zuobiao (15,1),location_a _zuobiao (15,2)+1.5,'15'); ti (16)=text (location_a _zuobiao (16,1),location_a _zuobiao (16,2)+1.5,'16'); ti (17)=text (location_a _zuobiao (17,1),location_a _zuobiao (17,2)+1.5,'17'); ti (18)=text (location_a _zuobiao (18,1),location_a _zuobiao (18,2)+1.5,'18'); ti (19)=text (location_a _zuobiao (19,1),location_a _zuobiao (19,2)+1.5,'19'); ti (20)=text (location_a _zuobiao (20,1),location_a _zuobiao (20,2)+1.5,'20'); ti (21)=text (location_a _zuobiao (21,1),location_a _zuobiao (21,2)+1.5,'21'); ti (22)=text (location_a _zuobiao (22,1),location_a _zuobiao (22,2)+1.5,'22'); ti (23)=text (location_a _zuobiao (23,1),location_a _zuobiao (23,2)+1.5,'23'); ti (24)=text (location_a _zuobiao (24,1),location_a _zuobiao (24,2)+1.5,'24'); ti (25)=text (location_a _zuobiao (25,1),location_a _zuobiao (25,2)+1.5,'25'); ti (26)=text (location_a _zuobiao (26,1),location_a _zuobiao (26,2)+1.5,'26'); ti (27)=text (location_a _zuobiao (27,1),location_a _zuobiao (27,2)+1.5,'27'); ti (28)=text (location_a _zuobiao (28,1),location_a _zuobiao (28,2)+1.5,'28'); ti (29)=text (location_a _zuobiao (29,1),location_a _zuobiao (29,2)+1.5,'29'); ti (30)=text (location_a _zuobiao (30,1),location_a _zuobiao (30,2)+1.5,'30'); ti (31)=text (location_a _zuobiao (31,1),location_a _zuobiao (31,2)+1.5,'31'); ti (32)=text (location_a _zuobiao (32,1),location_a _zuobiao (32,2)+1.5,'32'); ti (33)=text (location_a _zuobiao (33,1),location_a _zuobiao (33,2)+1.5,'33'); ti (34)=text (location_a _zuobiao (34,1),location_a _zuobiao (34,2)+1.5,'34'); ti (35)=text (location_a _zuobiao (35,1),location_a _zuobiao (35,2)+1.5,'35'); ti (36)=text (location_a _zuobiao (36,1),location_a _zuobiao (36,2)+1.5,'36'); ti (37)=text (location_a _zuobiao (37,1),location_a _zuobiao (37,2)+1.5,'37'); ti (38)=text (location_a _zuobiao (38,1),location_a _zuobiao (38,2)+1.5,'38'); ti (39)=text (location_a _zuobiao (39,1),location_a _zuobiao (39,2)+1.5,'39'); ti (40)=text (location_a _zuobiao (40,1),location_a _zuobiao (40,2)+1.5,'40'); ti (41)=text (location_a _zuobiao (41,1),location_a _zuobiao (41,2)+1.5,'41'); ti (42)=text (location_a _zuobiao (42,1),location_a _zuobiao (42,2)+1.5,'42'); ti (43)=text (location_a _zuobiao (43,1),location_a _zuobiao (43,2)+1.5,'43'); ti (44)=text (location_a _zuobiao (44,1),location_a _zuobiao (44,2)+1.5,'44');。

城市表层土壤重金属污染分析模型摘要：本模型研究人类活动对某城市不同功能区表层土壤重金属污染的影响程度及其主要原因。

通过计算各个功能区土壤重金属含量的单项污染指数，综合污染指数以及各重金属元素间的相关系数，以国家土壤污染等级划分标准来给出各功能区土壤的污染程度等级；并对各功能区的土壤污染情况作了对比。

对于问题（1）在综合分析了附件1，2，3所给出的数据后，应用Matlab 软件得出了8种主要重金属元素在该城区的空间分布（地势图），并在地势图上标出各重金属元素在对数值(某重金属实际的浓度/(1.5*相应的背景值))取对数后的值不小于1但小于2 （该值若在此范围内说明该地被污染）城区中的位置，由此我们可以直观分析该城区内不同区域重金属的污染程度，对其污染程度进行了等级的划定，同时给出了两个不同的污染程度评价模型，得出：a:每个功能区的重金属单项污染指数i p 从大到小依次顺序；b:每种重金属单项污染指数i p 在各个功能区从大到小依次顺序；c:各个功能区的综合污染指数；d:对单项污染指数和综合污染指数依据标准分别给出了污染程度分级表。

对于问题（2）是找出重金属污染的主要原因，我们利用处理问题（1）时得出的某项重金属的单项污染指数与地质累积指数，找出每个功能区的主要的污染源，在此基础上分析各功能区污染的原因，并给出合理分析结果。

同时根据内梅罗(N.L.Nemerow)污染指数考虑综合污染程度，得出工业区与交通区综合污染程度最高，其中Pb-Cd ，Ni-Cr 相关性显著，来自同一污染源的概率较大，具有较强的伴生关系，由此给出关于污染源由水流，空气，地渗扩散到其它功能区而导致其他功能区受污染的合理原因。

在解决问题（3）时，我们通过Matlab 软件作出8种重金属元素在城区内的浓度变化，粗略得出各重金属元素在土壤中传播快慢，Pb-Cd ，Ni-Cr 相关性显著，具有较强的伴生关系，因此在传播过程中有协同传播的特点，在污染源传播模型建立时对这种协同传播特性予以考查。

§4 一、名词解释
名词解释与符号说明
1．奖项等级：只比赛成绩划分的不同奖项，如一等奖，二等奖，三等奖，成功参 赛奖。 2．获奖比例：学校参加比赛获得某一奖项的队伍数量占所有队伍数量的比例。 3．比赛成绩：参赛队伍获得的比赛卷面成绩。 4．规模成绩：每个学校组织参赛的规模，主要包括组织参赛的队伍数量和参赛队 伍的获奖情况两个方面的因素。 5．综合实力：学校的综合实力主要是一个学校组织参赛的规模和比赛获得的奖项 状况决定的，所以学校的实力是比赛成绩与规模成绩的总和。
§3
零；
模型的假设
1.在安徽赛区的排名中， 假设专科组和本科组的记分标准一样， 不做另外分组处理； 2.假设如果一个学校那一年没有参赛， 则该年获得各个等级奖项的参赛队伍数记为 3.如果一个学校在某个奖项等级获奖空缺，也将参赛队伍记为零； 4.每年的考试难度没有差别； 5.每个同学的学习能力基本不变，并且发挥其真实水平； 6.影响学生成绩的因素主要有真实成绩与进步程度； 7.每个学生处于相同的考试环境中； 8.所给的数据时学校的真实考试成绩，没有作弊问题的影响。
n
Cj
i Wi
a
P X T N S
m ji
i
mn
优化因子 为某高校该年奖项的平均得分 为某高校该年所获奖项的总得分 为某高校该年每个学校参赛的队伍 为某高校该年所有参赛队伍的总数
§5
解。 一、问题一的分析与求解
模型的建立与求解
从所要解决的问题和对问题所做的假设出发， 分别对三个问题进行详细的分析与求
4
安徽科技学院 安徽理工大学 安徽绿海商务职业学院 安徽农业大学 安徽三联学院 安徽商贸职业技术学院 安徽师范大学 安徽新华学院 安徽新闻出版职业技术学院 安庆师范学院 蚌埠学院 亳州师范高等专科学校 巢湖学院 池州学院 滁州学院 阜阳师范学院 阜阳师范学院信息工程学院 合肥工业大学 合肥师范学院 合肥学院 河海大学文天学院 淮北师范大学 淮北师范大学信息学院 淮南联合大学 淮南师范学院 黄山学院 江淮学院 解放军电子工程学院 解放军陆军军官学院 六安职业技术学院 马鞍山师范高等专科学校 桐城师范高等专科学校 铜陵学院 皖西学院 芜湖信息技术职业学院 宿州学院 中国科学技术大学 ②年综合规模评比

城市表层土壤重金属污染分析摘要城市表层土壤重金属污染严重影响人的生活和生产，在已知319个样本的坐标以及各重金属污染物浓度情况下，按照各问题的要求，我们建立相应的模型来对具体问题进行求解。

首先为了了解该城市的地形，我们使用了MATLAB软件绘制各分区分布地形图,根据区域化变量理论，运用克里金插值法对未采样点的区域化变量的值进行线性无偏最优估计，绘制出各污染元素的拟合分布图。

在分析不同区域受污染程度时，以题中所给的元素背景值以及国家土壤环境质量准则为基础制定污染等级评判标准，据此进行单因子指数评价以及模糊数学评价分析，得出每个区域内主要污染元素以及污染程度。

接下来，为了分析污染的原因，我们将单因子指数的结果与聚类分析算法的结果进行比较，得出生活区内铬、铜、镉、锌含量较高；工业区内的主要污染元素为镉、铬、汞、镍、铅、锌；山区各重金属含量均较低；交通区内镉、铜、汞污染较重；在公园绿地区内，镉、汞、铅、锌的污染较为严重。

更进一步，为了确定每一种重金属元素的来源区域，我们结合重金属污染物在土壤内易沉积，不易传播的特性以及当地海拔高度，通过主成分分析法和相关分析法计算出污染物的所在区域，并与surfer的绘图点进行比较，发现两者有很好的吻合度，找到了Cd的污染主要来自工业区和交通区；Cr污染主要来自交通区和生活区；Cu污染主要来自交通区和工业区；Hg污染主要来自交通区；Ni污染主要来自交通区；Pb污染主要来自生活区和工业区；Zn污染主要来自交通区，生活区，工业区。

对于污染物传播特性的分析，我们首先用SPSS对各区域内各污染因子进行正态检验，然后对于正态分布和非正态分布的因子我们各给出一种函数描述其传播特性。

对于问题四，我们提出建立城市污染物数据库，记录城市污染物随时间和空间的变化情况，有利于对未知点进行时间和空间的预测。

关键字：重金属污染、克里金插值法、多因素分析法、聚类分析、主成分分析法一、问题重述城市是人类社会不断发展的产物,是人类重要的生活环境,也是人类生产生活的主要场所. 随着经济和社会发展速度的不断加快,城市环境问题,尤其是重金属作为一个威胁人类生存和城市发展的重要的城市环境污染物,引起了学者、专家们的关注.随着工农业发展和城市化进程的加剧, 通过农业灌溉、工业排放、交通运输、市政建设和大气沉降等造成城市土壤重金属的污染日趋严重。