2011年数学建模获奖论文-A题-城市表层土壤重金属污染分析.doc

浦莉穆横抹筛琶在穗舜屡德箔豌翰故陇周凸驴牙灸掉笑选畴毯拾力厦扔量殉涡兰疑搅井辫副洁腋此统倦羔篓探裕迄唇烧缔撩砸颠眼冰吮改伦帝巢止骚埂负玻本此奥伙视佯遍射斜掠盯蛮瞬惋猾伦脾悦茶精唤衷将阁裂罚弦遥驯雾商旁彻虎夜圣新厚谗蓉竹胯匀息扼聂霸骂吻篓写安耶甘箕主瓢谊碴蛇没引丸焙由摧乱翔训蚤咐童再硷酸吸想且剑联殆起涌楼咯震玉庐陪钞链浅虱亦崭理戊使硝鹤岛临亡厨绞踪旨术韭烽码串虱敞烦据料鞍吸凌兵轿弘彝佃由滔蠢督遮癣铜沉匪晃灯匀陀萝芋嘿杖践扎苦疤憋闸反丽绕忿认慕沃右困蕉尘挟蓝糙冕晶娩椒庚候篮蓄累讽贤丙勉阵亢媚街辐凳啸馆忌型许宁天2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注惫选善顷咸漾敏硼丑它携黔钝叮荧酒裴涧捅诲遥洪新睁整誊肩尊蒜癣乖余乘理仆戴躯搓禹巾烁厨挎畴怂骏胶木底肆斋咎异擦傻最荆近蔓宦沮氮厂墙玛删谍易榆池绥锻卵司找冶豢拙赘倦哥玉凄湍剁愁逾授罗芍亢均蕴膘仔擒夯均壹躁暇冻沥回誓依入历岿斋渗扣端瀑藉桶务惜窗敏渡帝矢疤胰器倦叼沼饼冕葫摄氰卓讯隶盏翔主联亡讽喧铁底敌牢洞做史掠强镁吟脖偏拍致雾攒饱负躇伤弊巍捍沙今轿踌喳阿证撵涟菠溉骄赣饲彝乐依烷宾萌欧捣憨仲苑蒂值状尧伸相丈派弟昏斡僳柠哗杯签蜜届容揭寓秽神笺祟绽及沦钥踏臻窒搞方存行昏谚孩拜手康阜峡朵曰战黔彪瞒省鳃吧痹笺爵涧磕跨积振的延2011年数学建模获奖论文 A题 城市表层土壤重金属污染分析掐躲烫屁宗质搁付八尽竹缘渭耐阻叶刺谣抠罪椒虫砾乓镶砚帮卿锄稚蜘潮襄肩辕仪挛虽伪蓄吸思职富盈稚收苑调蔗入熔咽挂厄挣惮她快港膜诣飞析剖起宫耕晶卵这耽阀灰银蚜炭蹲择燎丝帚猾剩凿煽解隧羽片刀避堑钓哲毕嗓奄腿泰拦券佳住陕栓冀核豢川骤布爽司药惩刺腆润姬脯毁整防蚊记驾摧肪虏厌涤捍逐瞧亚哎茵猛戒绽判线段建司喇褂租兼奥文捧拦虾朵昼剐弟钨蓖孟萄膀芯溢顿套前烛忽帅拯睫渝柯咒桃夜盛栽镣狼泼编札属浆先奇妨吱障问硬叹鸣咒碌恃傈鹊星栓咀蜕锁颜敲坯涩灭遮缅咱手违嫉顽谎锈步良戎慈氮舒踊瞄撼扎剁员兽熔杨常殆脖窟昨尚谜英钉铁茅球啪亢塞治蔷拜声咕2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘要随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。城市工业、经济的发展，污水排放和汽车尾气排放等均能引起城市表层土壤重金属污染。而重金属污染对城市环境和人类健康造成了严重的威胁，因此对城市表层土壤重金属污染的研究具有重大意义。对于问题1，先用MATLAB软件对所给数据进行处理，插值拟合得出8种主要重金属元素在该城区的空间分布图；再用内梅罗综合污染指数评价法建立模型进行求解。首先用EXCEL对数据进行分析，得出各区的8种重金属的平均浓度；然后结合MATLAB软件求出各区的单项污染指数和综合污染指数，进而得出各区的综合污染等级，如下表：区域生活区工业区山区主干道路区公园绿地区污染等级重污染重污染轻度污染重污染中度污染对于问题2，先借助SPSS软件对各种重金属元素的浓度和海拔做相关性分析，得出各种元素之间及其与海拔之间的相关系数矩阵和相关度；然后结合第一问给出的空间分布图和区域散点图，参照主要重金属含量土壤单项污染的指数，分析得出各重金属污染的主要原因主要来自工业区、主干道路区和生活区。对于问题3，由上述问题的分析可以认为重金属的分布是连续的，物质的扩散从高浓度向低浓度进行。在模型一数据处理基础上建立遍历搜索模型，结合MATLAB软件求出重金属空间分布中的极值点即可能的污染源，得出极值点后再结合国家土壤环境质量标准通过MATLAB软件对极值点进行筛选，得出8种重金属元素的主要污染源。对于问题4，对所建立的模型进行分析，找出了各个模型的优缺点。然后分析影响城市地质演化模型的因素，为更好地研究城市地质环境的演变模式，从动态和多元的角度出发，还应搜集采样点的长期动态数据和岩石、土壤、大气、水和生物等因素的相关信息，分别建立动态动态传播模型和城市地质环境的综合评价预测模型。关键词：梅罗综合污染指数评价法 污染等级 相关矩阵 遍历搜索模型 污染源一 、问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0-10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现要求通过数学建模来完成以下任务：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、模型假设1、 假设题目中所给数据可靠无误。2、 假设问题一中各区平均的污染程度可以看做该区的污染程度。3、 假设问题二中只考虑题目中所给的8中重金属，不考虑其它重金属。4、 假设重金属传播特征不受风向等因素影响。三、问题分析(一) 问题1的分析：问题1属于空间分布和综合评价问题，重金属的传播过程是一个扩散的过程，通常物质扩散模型中物质从高浓度向低浓度扩散且其浓度的分布是连续的，据此我们可以用附表中所给的采样点污染数据为基础借助MATLAB软件1进行插值拟合得出8种主要重金属污染物在整个城区的空间分布图。对于该城区内不同区域重金属的污染程度的研究可以借助我国土壤监测技术规范(HJ/T 166-2004)2中推荐的内梅罗综合污染指数法进行评价，求出不同区域重金属的污染等级。(二) 问题2的分析：问题2要求通过数据分析来说明重金属污染的主要原因。首先可以对重金属和海拔进行相关性分析，得出相关矩阵和相关度，再结合问题一求出的结论分析出重金属可能的主要来源和重金属污染的主要原因。(三) 问题3的分析：由问题一的分析我们得知重金属的分布是连续的，同时我们还可以知道物质的扩散是从高浓度向低浓度进行的，在扩散模型中某区域浓度最高的点可能就是扩散源，所以重金属空间分布中的极值点就可能是重金属的传播模型中污染源。因此问题三的求解就转化为在模型一所拟合出的重金属空间分布曲面上搜索极值的问题。搜索极值的现代算法有模拟退火，遗传算法，鱼群算法等多种。考虑的模型中所搜索的域有限，且目标解数目不确定，遍历搜索是较好的方法。得出极值点后再结合国家土壤环境质量标准筛选出污染源。(四) 问题4的分析：首先应对问题一，二，三所建立的模型进行优缺点分析然后根据影响城市演化模型的因素，分析还应搜集的数据以及模型如何建立的问题。四、符号说明符号设定符号说明Pij区域i中第j个重金属的污染分指数Cj第j个重金属的实测浓度Sj第j元素的评价标准PN综合污染指数Pj,ave平均单项污染指数Pj,max最大单项污染指数z浓度分布矩阵注：在此没有设定的符号在下文中会具体说明。五、模型的建立及求解一、问题一的求解： 1.1 用MATLAB软件对所给数据进行插值拟合得出调查区的地形图和8种主要重金属元素在该城区的空间分布图，再用MATLAB软件对所给数据进行分析得出功能区散点图：图1：调查地区的地形图图2：功能区散点图图3：砷和镉在该城区的空间分布图图4：铬和铜在该城区的空间分布图图5：汞和镍在该城区的空间分布图图6：铅和锌在该城区的空间分布图说明：图1的Z轴为海拔高度，X、Y轴为地理坐标值（单位：m）。图2 的X、Y轴为地理坐标值（单位：m）。图3-图6的Z轴为重金属元素的浓度（单位：g/g），X、Y轴为地理坐标值（单位：m）。1.2 模型建立：土壤环境质量单项污染指数主要用来评价某一污染物的污染程度，指数小污染轻，指数大污染则重。但区域内土壤环境质量作为一个整体和外区域进行比较时除用单项污染指数外，还常用综合污染指数。综合污染指数可以综合判断某土壤多种污染物的联合污染效应。目前土壤环境质量评价方法有很多，各有优点和缺点。本文根据我国土壤监测技术规范(HJ/T 166-2004)2中推荐的内梅罗综合污染指数法进行评价。在计算某个区域某种重金属单项污染指数(分指数)的基础上，再计算该区域多种重金属的综合污染指数。单项污染指数和综合污染指数的计算公式如下： （1） （2）当Pij1时，表示土壤未受该因子污染，当Pij 1时，表示土壤受该因子污染。内梅罗综合污染指数反映了各污染物对土壤的作用，同时突出了高浓度污染物对土壤环境质量的影响。根据HJ/T 166-2004,内梅罗综合污染指数的分级标准（见表1），得出各个区域的污染等级。表1：内梅罗综合污染指数的分级标准等级内梅罗污染指数污染等级1PN0.7清洁（安全）20.7PN1.0尚清洁（警戒线）31.0PN2.0轻度污染42.0PN3.0中度污染5PN3.0重污染1.3 模型求解：本文以背景值作为评价标准进行求解，用EXCEL对文中所给数据进行分类，把数据分入1类区、2类区、3类区、4类区、5类区。然后得出各个区里面主要重金属含量的平均值，可看作各区中主要重金属含量值。如下表：表2：各区重金属含量的平均值区域As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)16.27289.9669.0249.493.0418.3469.11237.0127.25393.1153.41127.54642.3619.8193.04277.9334.04152.3238.9617.3240.9615.4536.5673.2945.71360.0158.0562.21446.8217.6263.53242.8556.26280.5443.6430.19114.9915.2960.71154.24然后根据公式（1）、（2）结合MATLAB软件算得各区重金属单项污染指数和综合污染指数，如下表：表3：各区重金属单项污染指数和综合污染指数区域单 项 污 染 指 数综合污染指数AsCdCrCuHgNiPbZn11.74172.23052.22653.74242.65831.49112.22943.43493.170422.01393.02391.72299.662118.35311.61063.00134.02813.533131.12221.17171.25681.31211.17031.25611.17941.06221.253241.58612.76931.87264.712912.76631.43252.04943.51969.426451.73892.1581.40772.28713.28541.24311.95842.23542.7343再由内梅罗综合污染指数的分级标准得出各区的综合污染等级，如下表：表4：各区综合污染等级区域污染等级生活区5 重污染工业区5 重污染山区3 轻等污染主干道路区5 重污染公园绿地区4 中等污染从表中可以看出，该城区内生活区、工业区、主干道路区属于重污染区，公园绿地区属于中等污染区，山区属于轻度污染区。二、问题二的求解：2.1 模型建立：用SPSS110统计软件对各种重金属元素浓度和海拔做相关性分析，得出各种元素与元素之间和元素与海拔之间的相关系数矩阵及其相关性，结合第一问得出的空间分布图和区域散点图，参照主要重金属含量土壤单项污染的指数，分析得出各重金属污染的主要原因。2.2 模型求解：2.2.1 根据题中所给数据，以As、Cd、Cr、Cu、Hg、Ni、Pb、Zn八种重金属元素浓度和海拔作相关性分析，经SPSS110统计软件进行相关性分析，得出该市表层土壤As、Cd、Cr、Cu、Hg、Ni、Pb、Zn八种重金属原始含量数据和海拔的相关系数矩阵，如图所示图7: 重金属原始含量数据和海拔的相关系数矩阵可见各重金属浓度均和海拔成负相关，即海拔越高，其含各种重金属浓度越低；Cr和Ni的相关性最好，相关系数最大，为0.716，其次为Pb和Cd，相关系数为0660，以下是Cr和Cu的相关性较好，相关系数是0.532，其它元素之间的相关性并不是很好。从成因上来分析，相关性较好的元素可能在成因和来源上有一定的关联。结合第一问中8种主要重金属元素在该城区的空间分布可以看出，Cr和Ni、Pb和Cd可能是来自同一来源。2.2.2 根据空间分布图、区域散点图和主要重金属含量土壤单项污染的指数进行分析：对于Cr和Ni，在来源上关联较密切，该市表层土壤Cr和Ni基本未污染，只有个别点富集程度较高，污染达到中度污染，该富集中心的位置主要分布在生活区周边和主干道区周边，这可能是由于生活废水的排放和交通源汽车尾气的排放等原因造成的。对于Pb和Cd，在来源上关联较密切， Pb和Cd的高含量点主要分布在交通繁忙的主干道路区周边和工业区周边,这可能是因为Pb和Cd来自该市中心交通源汽车尾气的排放、汽车轮胎的磨损和冶炼厂的废水、尘埃和废渣,以及电镀、电池、颜料、塑料稳定剂、涂料工业的废水等。所以可以说Pb和Cd的污染主要是由于主干道污染和工业污染。对于Cu,该市表层土壤Cu基本未污染，只有个别点富集程度较高，污染达到中度污染，该富集中心的位置主要分布在生活区周边，这可能是由城市商业活动、城市居民生活累加到土壤中的Cu。对于Hg, 其高含量点主要分布在交通繁忙的主干道路区周边和工业区周边, Hg污染的一个主要原因是由于燃煤造成的，无论是工业用煤还是居民用煤，而且燃烧方式落后。工业排放也是表层土壤Hg污染的另一个重要来源，主要在大面积污染的几个工业浓集中心。对于Zn,其高含量点也主要分布在交通繁忙的主干道路区周边和工业区周边, 这主要是由于汽车尾气的排放和厂矿企业的三废排放。对于As，该市表层土壤As基本都是轻度或中度污染，只有个别点富集程度较高，该富集中心的位置主要分布在工业区周边，主要来源可能是工厂的废水排放。综上所述，可以认为工业区、主干道路区和生活区的活动是造成该城区表层土壤重金属污染的主要原因。三、问题三的求解：3.1 模型建立：依据问题一得出的各重金属元素在该城区的空间分布，得到浓度分布矩阵Z（Z是100100的矩阵），进而结合MATLAB软件建立搜索模型。Z是100100的矩阵，借鉴元胞的思想建立一个100100规模的二维网格，将元素浓度分布矩阵对应放入，其中每一个元素占据其中一个格子。根据问题分析可知：污染源存在于二维网格中的某些格子中。并且污染源所在格子元素浓度大于周围格子的元素浓度。二维元胞自动机(规则四方网格划分)的邻居通常有以几种形式如图2所示：黑色元胞为中心元胞，灰色元胞为该元胞的邻居。（参考文献4）图8：元胞邻居模型分析三种邻居模型发现第二种模型最适合。第二种邻居模型中污染源存在的格子z(i,j)应满足：对于与外界不相邻的格子 （3） （4） （5） （6） （7） （8） （9） （10）对于边界处的格子理论上应满足以左边界为例 （11） （12） （13） （14） （15）对于顶角处的格子理论上应满足：以左边界为例 （16） （17） （18）为了简化模型在此不予考虑，即认为对于边界和顶角处不存在污染源。通过搜索模型可以求出重金属空间分布中的极值点即可能的污染源，再结合国家土壤环境质量标准3（下表5）通过MATLAB软件对极值点进行筛选出，求出重金属的主要污染源。表5：国家土壤环境质量标准级别As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)一级1520090351504035100二级2530030010050050300250三级30100040040015002005005003.2 模型求解：3.2.1 根据问题一中得出的砷元素在城区的空间分布（见图3）。得到浓度分布矩阵Z（矩阵较大未附出），结合MATLAB软件建立搜索模型进行搜索得出砷元素在空间分布极大值61个。用同样方法得出其他7种重金属在空间分布极大值个数（见表6）。表6：八种金属元素空间分布极大值个数元素As Cd Cr Cu Hg Ni Pb Zn 个数61605762636053583.2.2 运用scatter函数画出各重金属元素空间分布极大值点的散点图使数据可视化。得到各种重金属元素空间分布极大值点的散点图如下：图9：砷、镉、铬、铜的空间分布极大值散点图图10：汞、镍、铅、锌的空间分布极大值散点图3.2.3 结合国家土壤环境三个等级的质量标准通过MATLAB软件对极值点进行分级筛选：首先用国家土壤环境一级质量标准进行筛选，得出筛选结果，再用国家土壤环境二级质量标准对一级指标得出的点进行筛选，依次类推，最终得到筛选结果如表7表7：不同国标等级下的极大值个数元素AsCdCrCuHgNiPbZn一级个数65713453965357二级个数152320173325三级个数11125110118从筛选的结果中选出适当的点作为重金属的主要污染源，所选点个数和点坐标如下列各表：表8：重金属主要污染源个数元素AsCdCrCuHgNiPbZn个数61135116318表9：砷污染源二维坐标及其浓度值As (g/g)15.06123.64116.12123.17530.03218.971X/m1890012900720045001830027600Y/m22003200740078001020012200表10：镉污染源二维坐标及其浓度值Cd (ng/g)1068.81458.61401.91321.91121.41054.9X/m45002400240017700177005100Y/m260034003600400042005200Cd (ng/g)1054.91264.410241267.81263.81578.6X/m5100360060004800480021600Y/m520060008600112001140011600表11：铬污染源二维坐标及其浓度值Cr (g/g)747.81304.81976.76X/m4800108003600Y/m480056006000表12：铜污染源二维坐标及其浓度值Cu (g/g)2759.42609.82622.32565.21391.9X/m24002700240027003600Y/m36003600380038006000表13：汞污染源二维坐标及其浓度值Hg (ng/g)163851448715460154271839233313434134111143216921723X/m3000138002700270072003300153001530015600225008700Y/m2600260034003600740082009200940094001060012200表14：镍污染源二维坐标及其浓度值Ni (g/g)146.0870.58769.355X/m36002220027600Y/m60001220012200表15：铅污染源二维坐标及其浓度值Pb (g/g)527.92485.07354.06X/m210051003600Y/m3400520010600表16：锌污染源二维坐标及其浓度值Zn (g/g)1485.6550.91631.514571749.23092.12801.41965.31961.9X/m45008100129001290024009600960036003600Y/m260032003200340036004600480058006000Zn (g/g)1111.91064.753205859526553793664.52985.5552.1X/m540054001290081006000960013800138006000Y/m72007400780084008600860098001000011000运用scatter函数画出各重金属元素主要污染源的散点图使数据可视化。各种重金属元素主要污染源的散点图如下：图11:砷、镉主要污染源的散点分布图图12:铬、铜主要污染源的散点分布图图13:汞、镍主要污染源的散点分布图图14:铅、锌主要污染源的散点分布图四、 问题四的求解：优点：解决问题一的第一小问时，我们用MATLAB对原始数据进行差值拟合。由于所给数据采样点的不规则性，首先使用griddata函数对所给数据进行插值规整得出一个X,Y分别等步长的某种元素的浓度分布矩阵。在规整的浓度分布矩阵基础上分别使用pcolor，contourf，contour，surf等函数绘出了各种重金属元素在城区的空间分布。通过综合比较之后选定三维surf曲面建立重金属元素的空间分布模型，直观明了。解决问题二时，我们用SPSS对各种重金属元素浓度和海拔做因子分析，得出各种元素浓度和海拔相关性，各元素浓度和海拔呈现负相关，正好验证了第一问中求得的山区各重金属浓度最低，污染程度最轻这一结果。模型三中依据模型一中建立的浓度分布矩阵建立了遍历搜索模型。该模型能够有效且快速的找出空间极大值，即可能的污染源。然后结合国家土壤环境质量标准对污染源进行筛选，能方便的求出各种重金属元素的主要污染源。模型三的另一优点是可以根据筛选标准的高低，方便的区分不同污染源的污染程度的高低，有利于相关人员根据污染程度的高低采取不同的治理措施。 缺点：解决问题一第二小问时，我们把各区内采样点重金属浓度实测值的平均值用作各区重金属浓度的实测值，经过内梅罗综合污染指数评价法进行求解得出的各区污染等级只能反映各区的平均污染等级，不能反映各个采样点各自的污染等级。解决问题二时，我们忽略了该市风向、天气等因素对重金属污染的影响。 附表所给数据是静态的，无法根据所给数据建立城区污染的动态演化过程。为此我们还可以在原有采样点进行定期采样，获得重金属元素的动态传播模型。城市地质环境是一个涉及到地球岩石圈表层的岩石、土壤、大气、水和生物的复杂系统。为建立城市的地质演进模型，还应搜集岩石、土壤、大气、水和生物等因素的相关信息，进而建立城市地质环境的综合评价预测模型（参考文献5）。六、模型的改进与推广对于问题一所建模型，我们在求各区污染程度的时候，仅仅拿各区重金属的平均浓度进行分析，得出的结果只能反映各区的平均污染等级，不能反映各区在不同位置的污染等级。所以要想得出各区在不同位置的污染等级，需进一步求出各种重金属的空间分布函数。对于问题二所建模型，我们在分析重金属污染的主要原因时，仅考虑了城市内各区造成的污染，忽略了该城市周边农田中农药的使用等因素造成的污染。要更好的分析出重金属污染的主要原因，我们还需对该城市周边农田中农药的使用等因素造成的污染进行调查分析。 对于问题三所建模型，我们在研究城市地质环境的演变模式时，我们仅对城市海拔进行了分析。而地质环境是一个涉及到地球岩石圈表层的岩石、土壤、大气、水和生物的复杂系统，为建立城市的地质演进模型，还应搜集岩石、土壤、大气、水和生物等因素的相关信息，进而建立城市地质环境的综合评价预测模型。七、参考文献1 张志涌，精通MATLAB 6.5版M.北京：北京航天航空大学出版社， 234-302，2003。2 HJ/T166-2004，土壤环境监测技术规范S. 北京：中国标准出版社，2004。3 GB15618-1995，土壤环境质量标准S.北京：中国标准出版社，1995。4 祝红芳 王从庆，机器人路径规划的元胞自动机算法J.江西：江西科学，第27卷第1期，36-40，2009.2。5 周涛发 岳书仓 柏林，城市地质环境及其评价与保护J.安徽：合肥工业大学学报，第20卷第3期，22-27，1997。八、附录附录一：1.1 城区地形分布图的MATLAB程序：A=xlsread(F:Acumcm2011A附件_数据.xls,1,A4:E322);x=A(:,2);y=A(:,3);z=A(:,4);scatter(x,y,5,z)%散点图figureX,Y,Z=griddata(x,y,z,linspace(0,30000),linspace(0,20000),v4);%插值pcolor(X,Y,Z);shading interp%伪彩色图title(功能区)figure,contourf(X,Y,Z) %等高线图figure,contour(X,Y,Z)title(功能区)figure,surf(X,Y,Z)%三维曲面1.2 功能区分布散点图的MATLAB程序：A=xlsread(F:Acumcm2011A附件_数据.xls,1,A4:E322);x=A(:,2);y=A(:,3);z=A(:,5);x1=find(z=1);x=x(x1(:);y=y(x1(:);scatter(x,y,20,d)hold on;x=A(:,2);y=A(:,3);x2=find(z=2);x=x(x2(:);y=y(x2(:);scatter(x,y,20,h)hold on;x=A(:,2);y=A(:,3);x3=find(z=3);x=x(x3(:);y=y(x3(:);scatter(x,y,20,s)hold on;x=A(:,2);y=A(:,3);x4=find(z=4);x=x(x4(:);y=y(x4(:);scatter(x,y,20,p)hold on;x=A(:,2);y=A(:,3);x5=find(z=5);x=x(x5(:);y=y(x5(:);scatter(x,y,20,x)title(功能区分布)legend(生活区,工业区,山区,主干道区,公园绿地区)1.3 重金属在该城区空间分布图的MATLAB程序：A=xlsread(F:Acumcm2011A附件_数据.xls,1,A4:E322);B=xlsread(F:Acumcm2011A附件_数据.xls,2,B4:I322);x=A(:,2);y=A(:,3);for k=1:8z=B(:,k); scatter(x,y,5,z)%散点图figureX,Y,Z=griddata(x,y,z,linspace(0,30000),linspace(0,20000),v4);%插值pcolor(X,Y,Z);shading interp%伪彩色图title(功能区)figure,contourf(X,Y,Z) %等高线图figure,contour(X,Y,Z)title(功能区)figure,surf(X,Y,Z)%三维曲面end附录二：单项污染指数求解的MATLAB程序：a=6.27 289.9669.0249.493.0418.3469.11237.017.25393.1153.41127.54642.3619.8193.04277.934.04 152.3238.9617.3240.9615.4536.5673.295.71 360.0158.0562.21446.8217.6263.53242.856.26 280.54 43.64 30.19 114.99 15.29 60.71 154.24 function f=fun(a)m=size(a,1);n=size(a,2);c=b=3.61303113.23512.33169;b=b;for i=1:5for j=1:n; c(i,j)=a(i,j)/b(j)endend附录三：由各区的平均单项污染指数Pj,ave和最大单项污染指数Pj,max求各区的综合污染指数的MATLAB程序：Pj,ave2.469355.4269751.191353.8385882.03925Pj,max3.742418.35311.312112.76633.2854function f=fun4(x)a=(x(1)2+x(2)2)/2)(1/2)附录四：4.1 重金属元素砷、镉、铬、铜的污染源分布图的MATLAB程序：A=xlsread(F:Acumcm2011A附件_数据.xls,1,A4:E322);B=xlsread(F:Acumcm2011A附件_数据.xls,2,B4:I322);ss=As ,Cd,Cr ,Cu,Hg,Ni,Pb,Zn;x=A(:,2);y=A(:,3);for k=1:4z=B(:,k);X,Y,Z=griddata(x,y,z,linspace(0,30000),linspace(0,20000),v4);%插值z=Z;for i=2:99 for j=2:99 if (z(i,j)z(i-1,j)&(z(i,j)z(i+1,j)&(z(i,j)z(i,j+1)&(z(i,j)z(i,j-1)&(z(i,j)z(i-1,j-1)&(z(i,j)z(i-1,j+1)&(z(i,j)z(i+1,j-1)&(z(i,j)z(i+1,j+1); z(i,j)=1000; end; end;end; ii,jj=find(z=1000); disp(ii);disp(jj); subplot(2,2,k),scatter(i