2011西京学院数学建模模拟竞赛题目A题硕士论文质量评价问题
硕士论文质量评价模型
硕士论文质量评价模型作者：刘铁来源：《价值工程》2011年第28期摘要：针对2011年西北工业大学校内数学建模竞赛题A题进行建模求解。
建立了合适的评价尺度函数，提出引入修正因子来克服论文评价中的情感因素的办法，完整地解决了题目中规定的问题。
对指导现实的论文质量控制有一定借鉴意义。
Abstract: Aimed at the 2011 Northwestern Polytechnical University campus mathematical modeling contest questions A, establish a model and solution. Give the proper evaluation scale function. In order to overcome the emotional factors in the evaluation of papers, we introduced correction factor into it. It completely solved the problems in the topic, which have a certain reference value for guiding practical paper quality control.关键词：硕士论文；质量评价；回归；修正因子Key words: master thesis；quality evaluation；regression；correct factor中图分类号：C533文献标识码：A 文章编号：1006-4311（2011）２8-0266-021 问题重述某校正开展硕士生质量评价，现搜集到2006、2007、2008年硕士生论文的评阅信息，分别按年存放在相关数据库中。
附件1和附件2中给出2006，2007，2008年各年硕士论文的评阅信息。
2011年全国大学生数学建模竞赛A题论文
城市表层土壤重金属污染分析摘要土壤作为城市环境的重要组成部分，不仅提供人类生存所需的各种营养物质，而且接受来自工业和生活废水、固体废物、农药化肥、及大气降尘等物质的污染．很容易导致金属元素的蓄积，从而造成土壤重金属的污染．本文讨论了某城市表层土壤重金属污染的空间分布分布状况，地区污染程度，以及污染传播特征，污染源等，建立了相应的几何与数学模型或算法，得到了较好的结果，为防治城市表层重金属污染，保护和提高土壤资源和生态环境，提供参考．对于问题一：通过给定数据的相关分析，不考虑地形高低对污染浓度变化的影响，用Matlab 软件编程绘制个重金属元素污染浓度空间分布三维网格图和二维等高线图，综合研究该城市各功能区的空间分布以及污染程度分布．建立了Muller 地积累指数模分析模型：)]/([log 2Bn C Fn ⨯=ℜ，确定污染程度水平分级标准，通过统计计算，分析了各重金属在不同功能区的污染状况及程度．结论是：主干道路区和工业区的重金属元素的污染最严重，其他次之．对于问题二，为说明重金属元素污染的主要原因，采用单因子指数模型和内梅罗综合指数模型进行综合指标评价分析，结合问题一中统计数据进行综合分析，得到个重金属元素在各功能区及城区的综合污染程度指标．污染最严重的功能区是主干道路区，其次按照污染程度从大小的顺序依次为：工业区、生活区、公园绿地区、山区．主干道路区土壤表层重金属元素含量很高，且种类多．根据地区的差异性和元素的特殊性，分析出重金属污染Hg 和Cu 污染是最严重的污染源，且污染最严重的地区在主干道路区和工业区．这些污染主要由于含铅汽油的燃烧、汽车轮胎磨损产生的含锌粉尘、工业污水的排放、生活废水的排放、化肥农药的多度使用、金属矿山的开采．详细情况见正文．对于问题三，为了找出该城区的污染源，在分析出重金属元素的主要传播特征之后，考虑大气空间传播情况，建立了微分方程模型，通过模型求解分析，用其等效的向内（向污染浓度较高的方向搜索）搜索算法，计算确定了重金属元素主要污染源的位置，其中As 较严重的中心污染源坐标分别为：（5291，7349，10）、（12696,3024，27）、（18134、10046、41）、（17814,10707,64）、（27700,11609,165）．这五个污染源主要分布在主干道路区．（5291，5739，10），（12696,3024，27），（17814,10707，64）分布在工业区，其它两种污染源分布在生活区．其余元素的中心污染源见正文．对于问题四，需对前面所建立的模型进行分析与评价并进行模型的优化，在详细分析了前三个问题的求解模型及过程之后，评价出所建立模型的优缺点．在问题三中，重金属元素除了在大气中传播以外，还通过水土流动传播．另外，前几个模型都是静态的，但污染物传播的过程与时间有关，是一个动态的过程．最后建立了一个扩算方程模型进行优化，能够为更好的研究城市地理环境的演变模式做贡献．关键词：重金属污染 地积累指数模型 单因子指数模型 内梅罗综合污染指数 微分方程模型一、问题重述1．1 问题背景随着工业发展和城市化进程的加剧，通过交通运输、工业排放、市政建设和大气沉降等造成城市重金属污染越来越严重．对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究城市不同功能区表层土壤重金属污染特征和污染空间分布性，以便更好的研究人类活动影响下城市地质环境的演变模式．本文就如何应用查证获得的海量数据资料展开城市环境质量评价，研究地质环境的演变模式建立数学模型．附录1列出了采样点的位置、海拔高度及其所属功能区等信息，附录2列出了8种列出了8种主要重金属元素在采样点处的浓度，附录3列出了8种主要重金属元素的背景值．1．2 需要解决的问题有(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度．(2) 通过数据分析，说明重金属污染的主要原因．(3) 分析重金属污染物的传播特征，由此建立数学模型，确定污染源的位置．(4) 分析所建立模型的优缺点，为更好的研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、问题分析该题目一方面通过GPS记录了该城市大量样本点的位置以及所属功能区，再应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据，通过这两个表的数据就大致可以提取出一些对于解决问题的重要信息，另一方面，题目给出了自然区各样本点的重金属元素的背景值，作为重金属污染情况的指标．对于分析研究各个样本点的污染程度至关重要．利用Matlab软件进行三维网格图和等高线图的制作并结合相关的数据统计分析，可以分析该城区不同区域重金属的污染程度．后面利用地积累指数法和内梅尔综合评价指数对数据进行处理，分析污染严重的功能区和重金属．结合图形的分析以及模型的建立综合分析重金属污染物的传播特征．接着对模型进行一定的优化处理，使得处理的结果更准确．三、模型假设1、假设题目所给的数据合理正确．2、该区域的划分是稳定的，不会出现大的变动．3、不考虑观测误差、随机误差和其他外在因素所产生的误差．4、重金属在大气中无穷空间扩散，不受风的影响，其扩散服从热传导定律．5、重金属污染程度连续变化，大气中重金属元素浓度连续变化．6、界限不明显区域有扩大、缩小、消失的过程，穿过大气进入仪器的重金属含量只有浓度大小之分，浓度大小由仪器灵敏度确定．四、变量与符号说明eo lg地积累指数n ()8,7,6,5,4,3,2,1=n 分别表示As,Cd,Cr,Cu,Hg,Ni,Pb,Zn 元素Fn 污染物重金属元素n 的浓度 Bn第n 种重金属元素的背景值上限P 综综合污染指数 n C重金属n 的实测值（ug/g ） max (/)n n C S 重金属污染物中污染指数最大值 (/)n n wg C S重金属污染物中污染指数平均值 n χ 重金属污染物n 的环境质量指数；n α 重金属污染物n 的实测值 n β 重金属污染物的评价标准． Ω 重金属元素通过的平面t 时间 h 海拔高度 V体积五、模型建立与求解针对问题一，首先想到的是用Matlab 软件编程，进行三维网格图、三维曲面图、等高线图和散点图的制作．5．1 问题(1)的分析、模型建立与求解： 5．1．1 问题（1）的分析对于问题一，首先来分析一下， 要给出8种主要重金属在该城区的空间分布， 就必须确定每个重金属元素与他们所对应的地区之间的联系．刚好题目给出了每个样本点的各元素浓度，那么 是不是可以将每种重金属元素含量浓度含量与该目标点所在的功能区建立联系？由此 想到利用Matlab 软件画出每种元素在该城区的三维曲面和空间曲面图．同时 在分析不同区域重金属的污染程度时，考虑到这个污染程度是否可以量化．并且是否能够建立一种模型将这种指标量化．这道问题还要求考虑每个功能区的污染程度， 知道每个功能区的每种重金属污染程度又是不一样的．那 通过什么指标来判断每个功能区的污染程度大小，这也是 为什么用权重作为评价每个功能区的污染程度的指标．5．1．2 问题（1）的模型建立该城区受这八种重金属元素As 、Cd 、Cr 、Cu 、Hg 、Ni 、Pb 、Zn 污染程度不一样．题目提供了每种重金属元素的背景值，那么 怎么样利用这些背景值和每种元素相关的浓度确定不同区域重金属的污染程度？所以 需要找出一种方法来准确的分析该城区内不同区域重金属的污染程度，并且最好能够量化．建立8种主要重金属元素在该城区的空间分布模型如下：引入了一种用于研究沉积物及其他物质中重金属污染程度的该区内不同地域重金属的污染程度的定量指标——地积累指数又称Muller 指数法，Muller 指数法表达式为:)]/([log 2Bn C Fn ⨯=ℜ式中Fn 表示污染物重金属元素n 的浓度；Bn 表示第n 种重金属元素的背景值上限，C 为考虑各地岩石差异可能会引起背景值的变动而取得一系列系数(一般取值为1．5)，用来表征沉积特征、岩石地质及其他影响．Muller 地积累指数评价和分级标准分级标准具体详见表1表1：地积累指数分级标准地积累指数ℜ 分级污染程度105≤ℜ<6及严重污染 54≤ℜ< 5强-及严重污染 43≤ℜ< 4强污染 32≤ℜ< 3中等-强污染 21≤ℜ< 2中等污染 10≤ℜ< 1轻度-中等污染 0≤ℜ 0无污染 该方法指标主要是通过每种重金属元素测得的实际浓度以及他们的相关背景值，计算出每种元素的地积累指数．然后根据上面这张表 就可以判断出每种元素的污染级别，这样就可以对每种元素的污染情况进行分析．然后 再利用Matlab 软件对题目所给数据进行处理，画出相应的网格曲面图和等高线曲线图．这里需要对Matlab 进行编程，首先利用每个样本点的横坐标、纵坐标、海拔高度建立等高线图，程序语句见附录一．通过该图，可以直观的看到该城区各功能区的空间分布．但是这张图不能反映出8种主要元素在城区的污染情况， 需要借助于各种主要元素的浓度．所以 需要再建立一张等高曲线图以及相应的网格曲面图，将主要元素的浓度作为第三坐标，命令语句见附录一．5．1．3 问题（1）的求解过程首先通过Matlab 软件，调用每个样本点的位置相关数据．就是以海拔为第三坐标，并且对每个功能区进行颜色区分，画出该城区每个功能区的二维等高线图．最后把每个样本点显示在图上．得到如下这张图：图一：重金属As空间二维等高线分布图这张图只反映出了该城区各功能区的空间分布，还不能看出每种重金属污染的情况．将每种重金属元素的浓度在图上反应出来，做出该城区重金属污染的二维等高线图．具体程序语句见附录二，得到如下这张图：图二：重金属As分布平面图同时为了对应这张As含量分布平面图，也画出了三维网格曲面图（图三）．图三：重金属As含量分布的空间三维图从空间三维图三中可以看到，有一处的波峰很高说明该处污染情况很严重，有二处处于波峰说明污染情况比较严重的主要有二处，还有一处面积比较广且所处高度稍微低一点这表明该处所受污染情况相对严重且污染的范围较广；同样分析二维等高线图二，图中有一处等高线之间的间距越来越密集且颜色很深表明该处受污染情况很严重，有二处等高线比较密集颜色相对较深表明这二处的污染情况相对严重，还有一处等高线间的距离较密集但是所包围的面积较广说明该处的污染也较严重且污染的面积很广．再结合前面的数据他们中心污染源的坐标分别为：（5291，5739），（12696,3024），（17814,10707）．都是分布在工业区，还有一处污染级别不是特别严重,但是在该处存在着污染源，此处刚好是山林密集区．通过观察图三，会发现刚好有三个点处于波峰，还有个点波峰稍微偏低，但还是能很直观的看出来．再来看一下，Cd这种重金属的城区各功能区的二维等高线图，分布平面图，空间分布图（图四、图五）：图四：重金属Cd空间二维等高线分布图图五：重金属Cd含量空间分布平面图以及相应的三维网格曲面图（图六）：图六：重金属Cd含量空间分布图从空间分布图六中可以看到，污染情况比较严重且面积比较广的主要有一处，还有五处污染也相对严重．以及几处小的污染；同样从二维分布图五可以看出等高曲线所谓面积有一处颜色很深，说明该区域污染情况很严重，同时也观察到又五处等高曲线所围的面积颜色比较深，这说明了这五处区域污染情况相对严重，很明显的是有一处等高曲线所围成的面积比较广且颜色较深，表明了该区域有一处污染情况较严重且污染面积比较广，由此可见不管是从二维还是三维图形进行分析的结果是相吻合的．再结合前面的数据它金属Cd中心污染源的坐标为：（22304,10527）．分布在主干道路区，还有一处污染级别不是特别严重．再观察图三，会发现刚好有三个点处于波峰．如此，通过同样的方法，都能够得到对其它六种种重金属在该城区的空间分布以及污染情况的了解（参见附录三）通过观察每种元素的三维曲面图以及等高曲线图．很容易观察到，每种重金属对该城区都存在或大或小的污染．其中有些地区是存在多种重金属污染，并且污染很严重，通过观察这8张图会发现这六种元素Cd,Cr,Cu,Hg,Ni,Pb 在横坐标在[3000,4000],纵坐标在[3000,6000]这个区域内含量都非常高，大致可以判定这段区域属于重度污染区．下面将题目中所给的数据用excel进行分类处理，得到样本点的地积累指数．然后运用数学统计法得到各种元素污染程度数据分布表，通过这些表就可以确定该城区内不同区域重金属的污染程度．统计该表时，是通过统计每个功能区的总样本点个数，然后通过地积累指数法分别计算出每种样本点的地积累指数，并判断他们的所在的污染级别．然后统计每种污染级别下，各功能区的污染点数占总点数的百分比也就是说的权重，通过该权重就能够分析出每种重金属元素的污染程度大小，以及污染所波及的范围．从而得到每种重金属元素污染最严重的地区．通过Excel对数据运算，得到重金属元素As 污染情况分布表：表二：As污染程度分布数据表下面通过同样的数据处理，得到Cd污染程度数据分布表：表三：Cd污染程度数据分布表其它六种元素的污染程度数据分布表见附录三．表中数值0的意义是在该污染级别下不存在观测的样本点．这是个大样本事件，可以认为该级别污染很轻微，甚至不存在这种级别的污染．而百分比越大，就说明在该污染级别下涉及的样本点比较多，污染波及范围较广．5．1．4问题（1）的结果分析5．1．4．1 As这种重金属污染情况分析由该表可以看出各个区域受As的污染程度，其中一类区即是生活区31．82%无污染，63．64%轻度—中度污染，4．55%为中等污染，无强污染和及严重污染的情况；二类区即是工业区38．89%不受重金属污染，52．78%受轻度—中度污染，5．56%受中等污染，2．78%受中等—强污染；三类区即是山区大多数不受污染，只有15．15%受轻度—中度污染，1．51%受中等污染；四类区即是主干路区47．83%不受污染，50．00%受轻度—中度污染，0．72%受中等污染，1．45%受中等—强污染；五类区即是公园绿地区大多数受轻度—中度污染，25．71%不受污染，2．86%受中等污染．再结合相应的几何图形，会发现在四区存在三个很明显的污染源，在污染源附近会看到，有很多二区的样本点．有个别一区的点，说明这种元素对一区的影响相对来说轻点．所以由分析可知工业区受污染最严重，污染面积达到了61．11%，其次是生活区、主干道路区，生活区污染面积都达到了50%以上，也就是说这三个区有至少一半的土壤受到该元素的不同程度的污染．其余功能区受污染程度就次之．5．1．4．2 Cd这种重金属污染情况分析由该表可以看出各个区域受Cd的污染程度，其中一类区即是生活区29．55%无污染，54．55%轻度—中度污染，13．64%为中等污染，无强污染和及严重污染的情况；二类区即是工业区16．77%不受重金属污染，44．44%受轻度—中度污染，30．56%受中等污染，8．33%受中等—强污染；三类区即是山区大多数不受污染，只有75．76%受轻度—中度污染，21．21%受中等污染；四类区即是主干路区23．91%不受污染，44．2%受轻度—中度污染，26．09%受中等污染，5．07%受中等—强污染；五类区即是公园绿地区大多数受轻度—中度污染，48．57%不受污染，31．43%受轻度-重度污染，11．43%受中等污染，8．57%受中等-强污染．再结合相应的几何图形，会发现在四区存在三个很明显的污染源，在污染源附近会看到，有很多二区的样本点．有个别一区的点，说明这种元素对一区的影响相对来说轻点．所以由分析可知工业区受污染最严重，污染面积达到了61．11%，其次是生活区、主干道路区，生活区污染面积都达到了50%以上，也就是说这三个区有至少一半的土壤受到该元素的不同程度的污染．其余功能区受污染程度就次之．5．1 这六种重金属Cr、Cu、Hg、Ni、Pb、Zn污染情况分析由于重金属含量越多，说明该地区的重金属污染程度越严重．Cr污染最严重的有一处，该中心污染源的坐标为：（3299，6018），所在地区为主干道路区，一定程度上波及到了生活区．一区和四区存在强-及严重污染，一区波及面积达到了52．27%，四区波及面积达到了41．3%，该元素污染最严重的就是生活区．Cu污染最严重的有一处，该中心污染源的坐标为：（2427，3971），所在地区为生活区，一定程度上波及到了工业区和主干道路区．一区和四区存在及严重污染，一区污染波及范围达到了84．09%，四区污染波及范围达到了84．06%，该元素污染最严重的就是生活区和主干道路区．Hg污染最严重的有一处，中心污染源的坐标为：（3299，6018），所在地区为主干道路区，一定程度上波及到了生活区．一区和四区存在及严重污染，一区污染波及范围达到了54．55%，四区污染波及范围达到了50．74%，该元素污染最严重的就是主干道路区．Ni污染最严重的有一处，中心污染源的坐标为：（3299，6018），所在地区为主干道路区，一定程度上波及到了生活区．一区、二区和四区存在及严重污染，一区污染波及范围达到了90．91%，二区污染波及范围达到了94．44%，四区污染波及范围达到了93．48%，该元素污染最严重的就是主干道路区和生活区．Pb污染最严重的有二处，中心污染源的坐标为：（2383，3692）、（5062，4339），所在地区为生活区和主干道路区，一定程度上波及到了工业区．一区和四区存在及严重污染，一区污染波及范围达到了52．73%，四区污染波及范围达到了80．87%，该元素污染最严重的就是主干道路区．Zn污染最严重的有一处，中心污染源的坐标为：（14065，10987），所在地区为主干道路区，一定程度上波及到了工业区．四区存在及严重污染，四区污染波及范围达到了67．39%，该元素污染最严重的就是主干道路区．所以，该城区不同区域重金属污染最严重的区域是主干道路区和工业区，其次是生活区、公园绿地区、山区．5．2 问题（2）的求解：5．2．1问题（2）的分析通过问题一的分析，可粗劣的判断哪几种元素污染比较大，哪个功能区污染比较严重，但是怎么样才能具体到哪个功能区污染最严重，被污染的功能区的土壤哪种重金属污染最严重？所以，针对问题二给出的数据分析，不能简单的进行数据处理．为了使得所寻找出来的原因更有说服力，用两种方法分别进行说明和验证，还要进行综合指标评价．最后确定了最严重的污染地区以及污染最严重的相关元素，根据地区的差异性和元素的特殊性，才能说明重金属污染的主要原因．5．2．2数据的统计分析首先通过数据的处理，建立每个功能区各重金属元素的污染程度样本所占的百分比表．一功能区的相关百分比数据如下：表四：一功能区各重金属污染程度所占百分比在此功能区从总体来看，重金属污染程度处于中等-强污染，其中主要污染来自重金属元素Ni，另外在该区域有少数地方Cu污染及严重．表五： 二功能区各重金属污染程度所占百分比在该功能区重金属Hg 和重金属Ni 的污染极为严重，尤其是在该区域的某些地方．由此可见，在此功能区照成重金属污染的罪魁祸首为重金属元素Hg 和重金属元素Ni ． 通过这两张表， 会发现有些地区之所以污染严重，主要是因为个别元素污染所导致的．所以 要分析重金属污染的原因，就得分析该重金属在该功能区为什么会产生污染．其它三个功能区各重金属污染程度百分比见附录三．通过该附录表 可以看到在该功能区里，重金属污染程度较轻，污染等级集中在轻度污染及以下． 再观察功能区四，重金属污染十分严重，大多数重金属污染元素都集中在在各个功能区，但是在这个功能区，Pb 污染级别比较轻，没有中度甚至以上级别的污染． 再看功能区五，从总体上分析，该地区重金属污染中等、强污染几乎没有，正因如此造成重金属污染的少数种类重金属元素就凸显出来了——Ni 元素和Hg 元素．纵观整体，分析所有的功能区， 很容易发现造成重金属污染的主要重金属元素，他们就是Ni 元素和Hg 元素．知道前面的数据分析理由不充分，只是一个粗劣的判断．为了综合前面处理的数据，准确找出各个功能区污染的主要元素． 需要利用单因子指数法和内梅罗综合污染指数法进行综合评价．5．2．3 单因子指数法和内梅罗综合污染指数法的建立与求解单因子指数法是目前国内土壤重金属的单项污染指数评价方法之一，其计算公式为：n n n βαχ=，式中n χ为重金属污染物n 的环境质量指数；n α为重金属污染物i 的实测值；n β为重金属污染物的评价标准．n χ﹥1表示污染；n χ=1或n χ﹤1表示无污染；且n χ值越大，则污染物越严重．为了更全面的反应各重金属对土壤的不同作用．突出高浓度重金属对环境质量的影响， 采用内梅罗综合污染指数法．其计算公式为:2)/(/22max n wgn n n S C S C P +=）（综,式中max )(n n βα表示重金属污染物种污染指数nn βα的最大值；(/)n n wg C S 表示重金属污染物中污染指数的平均值．土壤污染水平分级标准采用国家土壤环境二级标准．土壤污染综合污染指数分级标准为综合污染指数>3为重污染，2~3为中污染，1~2为轻污染，0．7~1为警戒级，≤ 0．7为安全级．下面为了找到每种元素在该城区的综合污染指数，借助于Matlab 循环计算．编写如下系列命令见附录七．运行程序结果为As 综合污染指数：p=4．0093，分别运行另外几种程序，得到每种重金属元素的综合评价指标，简单结果如下表：。
2011数学建模A题程序及答案参考
程序：clcclearclose all% ZuoBiao.dat和JinShu.dat为导入的两个题目所给数据文件x=D(:,2);y=D(:,3);z=D(:,4);c=D(:,5);xi=linspace(min(x),max(x),100);yi=linspace(min(y),max(y),100);[xi,yi]=meshgrid(xi,yi);zi=griddata(x,y,z,xi,yi);ci=griddata(x,y,c,xi,yi);marker={'*','o','s','^','p'};color={'k','r','m','c','b'};mat={'As','Cd','Cr','Cu','Hg','Ni','Pb','Zn'};str={'等高线','1生活区','2工业区','3山区','4交通区','5公园绿地区'};for j=1:8% 等高线图形figurecontourf(xi,yi,zi,0:10:500);% set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2)% clabel(C,h,[0:10:50,50:50:300,300:100:500]);title(['重金属',mat{j}, ' 在该区域的空间分布图'])xlabel('X')ylabel('Y')colormap hsvcolorbargrid onhold onfor i=1:5loc=c==i;plot(x(loc),y(loc),marker{i},'markerfacecolor',color{i},'MarkerEdgeColor',color{i});endlegend(str,'location','best') for k=1:length(x)text(x(k)-200,y(k)+200,num2str(M(k,j+1)),'fontsize',8); end end % 三维体图figureh=surf(xi,yi,zi); set(h,'cdata',ci); colormap hsvtitle('三维图立体（颜色条表示分类）') xlabel('X') ylabel('Y') colorbar hidden off hold onfor i=1:5 loc=c==i;plot3(x(loc),y(loc),z(loc),marker{i},'markerfacecolor',color{i}); endstr{1}='三维图';legend(str,'location','best')第1章 污染物在水体中迁移模型1.1 流体运动的某些概念 1.1.1 浓度设V ∆是以点(,,)x y z 为中心的微小体积，M ∆是该微小体积内包含的污染物的质量，某时刻t ，点(,,)x y z 的浓度定义为0(,,,)limV MC x y z t V∆→∆=∆(1-1-1)某一时刻污染物质在水域中的分布，一般来说是随着空间位置的变化而变化的。
2011数学建模A题全国一等奖论文(城市重金属污染及地质演变)
Hg 1.41 2 4.197 5 0.226 1 3.672 4 1.645 2
Ni 0.58 1 0.69 1 0.33 1 0.52 1 0.28 1
Pb 1.156 2 1.585 2 0.237 1 1.032 2 0.957 1
Zn 1.78 2 2.01 3 0.087 1 1.812 2 1.081 2
Hg 93.04 642.3 40.95 446.1 109.4
Ni 18.34 19.81 15.45 17.59 14.88
Pb 69.10 93.04 36.55 63.39 60.19
Zn 237.0 277.9 73.2 242.3 145.9
把 Bn , C n 代入到式 （1 ） 后， 得到 5 块区域内各种重金属元素的地积累指数值 （ I geo ） ， 如下表：
图 1：As 元素在各区域的分布 3
从图 1 中可以直观地看出 As 元素主要分布在工业区、生活区与绿地周围，而在交 通区分布较少。在山区大部分地区分布较少，少部分地方有分布。
图 2：Cd 元素在各区域的分布
对于图 2， 可以直观地看出 Cd 元素主要分布在工业区与交通区周围， 次之分布在生 活区与绿地，在山区分布较少。 另外，由其它 6 种重金属元素的分布图（见附件—图像—第 1 问）可以看出 Cr 元 素主要分布在生活区与交通区，在工业区、绿地与山区分布较少。 对于 Cu 元素，发现工业区与交通区分布较多，其它地方分布较少。 对于 Hg 元素，发现只有少部分生活区与交通区分布较多，其它部分基本上没有受 到 Hg 元素的污染。 对于 Ni 元素， 发现在山区、 生活区与工业区分布较多， 在交通区与绿地分布较少。 对于 Pb 元素，发现在交通区分布最多，在部分生活区与部分工业区也有所分布， 而在其它区域都分布较少。 对于 Zn 元素，发现在生活区与工业区分布较多，其它区域分布较少。 4.1.2 地累积指数法评价不同区域重金属污染程度 地累积指数法是一种常用的研究沉积物中重金属污染的定量指标，其计算公式为：
2011数学建模A题论文研究报告
2011数学建模A题论文研究报告1．问题的探究本题是研究某一地区重金属的污染情况，从问题本质来看问题可以大致看做是对待解决问题的初步认识与描述；其次就是对该问题逐步深入进行探究，比如原因等；然后建立问题的数学模型；最后紧接上题，对该模型的探究，譬如可以是模型的修正与推广。
本人认为这样的提问方式一个是循序渐进、不断深入，符合常规，更易于读者理解和思考。
就本题来讲第一题是要求给出不同重金属在该地的空间分布，并建立能表征污染程度的指标，读者看到这样的第一题就会有亲切感，因为问题所需的数据已经给出，只需要对数据进行分析就可，这样就不会让读者感到无从下手；第二题是要求分析污染的原因，引导读者向更深的方向对问题进行探究；第三题就步入正题了，建立重金属污染的传播方式的数学模型，确定金属的污染源；最后是对模型优缺点探究和更加广发的应用。
1.1问题一的探究问题一是得到重要金属的空间分布图并建立重金属的污染指标，我们首先讨论解决该问题的结果是得到该地区不同地理位置的金属含量高低，对于该问题的解决不可避免的就是要运用给出的已知数据，由于无法得知数据的给出情况，这里不再讨论数据的处理方法。
浓度的数据不难得出，本人认为关键之处就是怎样把庞大的数据清晰、简明的标注在相应的地理位置上，便于读者阅读。
本人的想法是建立一个三维的空间分布图，X、Y轴分布表示经度、纬度，Z轴表示该地理位置的重金属浓度，这样在图上就可以标明给出数据，然后通过拟合的方法便可得到一个浓度平面，达到清晰明了的读出金属的浓度的效果。
本方法虽然不失为一个良策，但由于本人能力有限，无法用科学方法得当此图。
在论文中本题的解决采用等值线的方法来描述重金属的空间分布图。
利用三角线性插值的方法可以得到浓度的等值图。
浓度的等值图可以科学清晰的反映金属浓度的空间分布情况，并且等值线的疏密可以体现数据变化的速率，越密说明变化越快，最密点越有可能是污染源；并且可以找出污染严重的区域，结合已知的城市功能分区和地形可以定性分析出该地重金属污染的原因，可以为后题的解答提供有效的依据，和答题方向。
天然肠衣搭配问题 2011年全国大学生数学建模竞赛 A题 优秀论文设计
天然肠衣搭配问题摘要本文针对天然肠衣原料的搭配方案进展设计，充分考虑最优化原如此，在满足搭配方案具体要求同时兼顾效率的情况下，设计线性规划模型，并借助软件Lingo求解出最理想的捆数与搭配方案。
对于题目给出的五个具体要求，我们经过分析之后将其划分优先级，逐层递进地找出答案。
首先我们将条件〔1〕设为最优先条件即对于给定的一批原料，装出的成品捆数越多越好。
在此根底上，条件〔2〕的优先级次之。
对于条件〔3〕和〔4〕，我们经过讨论后认为其意在于放宽较为苛刻的长度与每捆根数要求以符合实际生产。
因而理想情况应是所有捆的根数与长度都恰好满足规格。
当由于给定数据原因使得理想情况不能实现时，再考虑放宽剩余原料的组装长度与根数要求，条件〔3〕与〔4〕的优先级最次。
在建模过程中，我们先对各规格在不考虑〔3〕与〔4〕的情况下进展线性规划，求每种每捆可行搭配方案所能组装出的最大捆数，再将其加和得出各规格的最大捆数。
这种方法在数据量较大的情况下兼顾了准确度与效率。
对上述不能组合的剩余材料我们如此放宽条件。
因条件〔2〕要求最短长度最长的成品数量尽可能多，再结合条件〔4〕中原料可以降级使用的规如此，故我们采用先从规格三的剩余原料考虑，再依次降级并入次级的原料使用考虑搭配。
由于剩余材料数量较少，故可以不必考虑效率问题。
最后满足条件〔5〕将结果求解。
利用上述模型和Lingo软件最后求解出了最大捆数183。
并可以根据原料数量求出具体的搭配方案。
关键词：搭配方案线性规划 Lingo1.问题重述天然肠衣〔以下简称肠衣〕制作加工是我国的一个传统产业，出口量占世界首位。
肠衣经过清洗整理后被分割成长度不等的小段〔原料〕，进入组装工序。
传统的生产方式依靠人工，边丈量原料长度边心算，将原材料按指定根数和总长度组装出成品〔捆〕。
为了提高生产效率，公司计划改变组装工艺，先丈量所有原料，建立一个原料表。
原料按长度分档，通常以为一档，如：3按3米计算，按计算，其余的依此类推。
2011数学建模国赛A题
城市表层土壤重金属污染分析与预测摘要近几十年，由于工业化和城市化的进程加快，土壤重金属污染日益严重。
本文针对重金属污染问题，主要包括四个相关问题，分别为确定重金属元素的空间分布，并分析该城区内不同区域重金属的污染程度，说明重金属污染的主要原因，确定污染源的位置，对模型进行评价与扩展并分析该城市地质环境演变模式。
我们在合理的假设下，先绘制出重金属元素空间分布图，再根据单因子指数与内梅罗综合污染指数分析方法对各功能区污染程度分析，说明重金属污染的主要原因。
接着，在此基础上建立二维污染物空间扩散模型，采用搜索算法在污染较为严重的区域搜索出适量的样本点，再采用多元非线性回归方法对二维扩散模型进行拟合，求解相应的参数。
最后通过对第三问模型优缺点的分析，建立时空扩散模型，更加精确反映城市地质环境的演变情况。
针对问题一，首先用软件Surfer绘制出该城区的空间分布图，得出该城区大致的地形地貌，再将纵坐标改为八种重金属元素的含量，采用克里格插值，得到各种重金属在城区的分布图，得出各重金属在不同功能区的污染物空间分布特点。
再选用单因子指数和内梅罗综合指数评价污染状况，由数据可以反映出工业区和交通区Hg和Cu严重超标，各功能区污染程度分别为：生活区、山区、交通区及公园绿地区受到了轻度污染，工业区受到了重度污染。
针对问题二，对各功能区的污染物含量进行分析，分析推测该城市可能存在规模较大或者数量较多的Hg,Cu,Zn，Pb严重超标的工厂。
为进一步揭示表层土壤重金属污染来源，我们对三个污染相对严重的功能区的超标较严重的重金属元素Hg,Cu,Zn,Pb进行相关性分析，运用SPSS软件计算出各相关系数，揭示出不同功能区内土壤表层重金属污染的不同来源，分析可得，交通区污染主要来源于汽车尾气排排放，工业区污染主要来源于工厂的超标排放。
针对问题三，根据污染物扩散的特征，建立一维污染物扩散模型：Q(x)=Q0eαx−12x2x m，以此推广成污染物的高斯二维空间扩散模型：Q(x,y)=Q0exp[−(x−x o)2−(y−y0)22σ2] 。
2011全国大学生数学建模竞赛A题获奖论文——一篇
城市表层土壤重金属污染分析的数学模型摘要为研究城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式。
本文通过处理和分析已给数据，给出金属的空间分布说明污染程度和主要原因；建立数学模型确定污染源位置；最后收集其他信息讨论城市地质环境的演变模式。
问题一，利用matlab软件作出位置坐标x、y与八种总金属元素浓度的空间分布图；分析采集的重金属元素浓度所在区域的大致情形。
对采集的重金属元素浓度的数据进行分析，并计算单因子和多因子污染指数，根据土壤污染分级标准判断出不同重金属元素在各功能区的污染程度和各功能区的综合污染程度，其中工业区中铜是所有元素在不同功能区中污染程度最严重的，而工业区和交通区的综合污染程度是最严重的。
问题二，首先利用SAS软件对八种重金属元素在五个城区的含量进行主成分分析，得到八种重金属对各功能区的贡献率，可初步推断出工业生产、交通设施和生活垃圾造成重金属污染。
再利用SAS软件对各城区的重金属进行因子分析，进一步判断八种不同重金属污染的原因，如汞污染的原因为工业生产中三废的排放、交通运输业中汽油的燃烧和汽车轮胎磨损产生的粉尘等。
问题三，根据所给数据，分析重金属污染传播特征，即分别是介质的迁移运动、污染物的分散运动、污染物的累积与转化、污染物被环境介质吸收或吸附、污染物的沉淀，然后利用Matlab软件，采用多元纯二次二项式回归分析方法，分别得到每种重金属元素浓度与坐标的回归方程，并根据该方程利用多元函数求极值的方法确定出污染源的可能位置分别为：As（1878.2634，6003.7263，4.5846），Cd(970.5835，3946.7518，6.5891)，Cr（1235.1956，2658.3427，8.5402），Cu（138.4682，6223.4521，3.2461），Hg （1231.5782，2561.5483，5.2478），Ni（12234.2587，5865.1656，23.2461），Pb （2310.68914145.2674，3.2651），Zn（3015.43418642.2365 5.0543）；问题四，基于前三问，分析所建模型的优缺点。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2011西京学院数学建模模拟竞赛题目
A题：硕士论文质量评价问题
我国自1980年建立新的学位制度以来，已初步形成了具有我国特色的研究生招生和培养模式，并且随着社会环境的变化和人才培养的不同要求适时作出调整。
如：入学类型开始有了在职人员和非在职人员、定向培养和非定向培养、直博和提前攻博等类型；招生工作中的计划内定向培养意识逐渐淡化，出现了在职攻读学位等进校不离岗的招生形式等等。
然而,这种多样化的招生和培养方式，也使得研究生论文质量问题日益引起人们关注。
特别是近年来研究生招生规模的迅速扩张,以及研究生在国家发展与社会进步中发挥的越来越大作用，更使研究生论文质量问题成为人们关注的焦点。
如何建立合理的研究生论文质量评价体系，并通过量化的手段找出当前国家在研究生招生过程中存在的问题，进而调整招生政策，改革招生方式，真正吸收综合素质高和研究能力强的优秀学生进入研究生队伍，已成为保证、提高硕士研究生论文质量的第一大关口，是国家乃至个人都十分关注的一项课题。
为全面贯彻科学发展观，落实以质量为核心的发展要求，全面分析和评价我国硕士生质量，制定进一步提高硕士研究生教育质量的政策，需要对硕士生的招生质量、论文质量、培养质量等进行综合评价。
某校正开展硕士生质量评价，现搜集到2006、2007、2008年硕士生论文的评阅信息，分别按年存放在相关数据库中。
附件1和附件2中给出2006，2007，2008年各年硕士论文的评阅信息。
全部存放在Excel表中。
请根据这些信息分析解决以下问题。
1．对2006,2007,2008年各年硕士生论文选题与开题进行总体评价。
包括各专业的评价和各年的总体评价。
2．对2006，2007，2008年各年硕士论文评分的评价。
包括各专业与各年的总体评价。
3．对各专业、各年硕士论文选题开题与论文得分之间的相关性进行分析，你从中得出什么结论？
4．对2006，2007，2008年复审（毕业后的重新评阅）论文的评价。
包括各专业与各年的总体评价。
5．对硕士毕业前后论文的评分结果进行分析，你得出什么样的评论。
说明你的观点与结论。
对此你有什么建议。
附件中数据说明：
表1 硕士论文选题的评分标准
表2 硕士论文开题报告指标
表3 硕士学位论文评分表
表4 等级划分表
附件1 2006_2007_2008年硕士生论文抽查表(Excel表)
数据说明
KT------代表表1中开题的评价结果
XT1,XT2,….,XT6 -----分别代表表2中6项开题的评价结果
A1,A2,….,A7-------分别代表表3中第一位专家的盲审评阅结果
B1,B2,….,B7-------分别代表表3中第二位专家的盲审评阅结果
每篇抽查的硕士论文都被两位专家评审，并有7项打分。
附件2 2006_2007_2008年硕士论文复审抽查表(Excel表)
A1,A2,….,A7-------分别代表表3中专家的盲审评阅结果，复审时每篇论文只有一位专家评审.
送审编号---前两位代表毕业的年，3~4位代表所属学科的序号，5~6位代表该年该学科的送审论文序号。
如080101表示2008年的论文，所属学科是1，在该学科中论文序号为1。