[数学]金属污染分析.ppt

生活区：内梅罗综合指数为3.17，属于重度污染区，其Cu、Zn的污染级别为重度，其他6种元素的污染级别为中度； 工业区：内梅罗综合指数为13.53，属于重度污染区，Hg、Cu为重度污染，Cu和Zn为中度污染，Cr和Ni为无污染，As、Cd和Pb为轻度污染。 从表中可以看出，该城区内生活区、工业区、主干道路区属于重污染区，公园绿地区属于中等污染区，山区属于轻度污染区。 为了更形象地体现不同区域重金属的污染程度，将表6中的数据做成直方图如下： 各区域污染指数图 5.3 重金属污染的主要原因 5.3.1聚类分析模型建立 聚类分析又称群分析，它是研究（样品或指标）分类问题的一种多元统计方法，所谓类，通俗地说，就是指相似元素的集合。严格的数学定义是较麻烦的，在不同问题中类的定义是不同的。聚类分析起源于分类学，在考古的分类学中，人们主要依靠经验和专业知识来实现分类。随着生产技术和科学的发展，人类的认识不断加深，分类越来越细，要求也越来越高，有时光凭经验和专业知识是不能进行确切分类的，往往需要定性和定量分析结合起来去分类，于是数学工具逐渐被引进分类学中，形成了数值分类学。后来随着多元分析的引进，聚类分析又逐渐从数值分类学中分离出来而形成一个相对独立的分支。。聚类分析内容 非常丰富，有系统聚类法、有序样品聚类法、动态聚类法、模糊聚类法、图论聚类法、聚类预报法等。本文主要介绍常用的系统聚类法。聚类分析有两种：一种是对样本的分类，称为Q型，另一种是对变量（指标）的分类，称为R型。本文考虑八种金属元素在五个区域上的污染浓度变化的不同，是对金属元素进行聚类，因此应采用Q型聚类方法。Q型聚类分析的主要作用：⒈可以综合利用多个变量的信息对样本进行分析。⒉分类结果直观，聚类谱系图清楚地表现数值分类结果。⒊聚类分析所得到的结果比传统分类方法更细致、全面、合理。 经过这种标准化所得的新数据，各要素的极大值为1，极小值为0，其余的数值均在0与1之间。 （2）定义样品之间距离，计算样品两两距离，得一距离阵记为 ，开始每个样品自成一类，显然这时 。 （3）找出 的非对角线最小元素，设为 ，则将 和 合并成一个新类，记为 ，即 。 我们采取系统聚类方法中的最短距离法。 最短距离法聚类的步骤如下（1）首先对数据进行极差的标准化 （4）给出计算新类与其它类的距离公式： 将 中第p、q行及p、q列用上面公式并成一个新行新列，新行新列对应 ，所得到的矩阵记为 。 （5）对 重复上述对 的（2）、（3）两步得 ；如此下去，直到所有的元素并成一类为止。如果某一步 中非对角线最小的元素不止一个，则对应这些最小元素的类可以同时合并。 5.4.2 聚类分析模型求解 八种金属元素在五个区域的平均浓度，如表9所示 元素 1生活区 2工业区 3山区 4主干道路区 5公园绿地区 As (μg/g) 6.27 7.25 4.04 5.71 6.26 Cd (ng/g) 289.96 393.11 152.32 360.01 280.54 Cr (μg/g) 69.02 53.41 38.96 58.05 43.64 Cu (μg/g) 49.40 127.54 17.32 62.21 30.19 Hg (ng/g) 93.04 642.36 40.96 446.82 114.99 Ni (μg/g) 18.34 19.81 15.45 17.62 15.29 Pb (μg/g) 69.11 93.04 36.56 63.53 60.71 Zn (μg/g) 237.01 277.93 73.29 242.85 154.24 平均浓度表 用excel对数据进行极差标准化处理，得到 结果，如表10所示： 元素 1生活区 2工业区 3山区 4主干道路区 5公园绿地区 As (μg/g) 0.00 0.01 0.00 0.00 0.00 Cd (ng/g) 0.45 0.61 0.23 0.56 0.43 Cr (μg/g) 0.10 0.08 0.05 0.08 0.06 Cu (μg/g) 0.07 0.19 0.02 0.09 0.04 Hg (ng/g) 0.14 1.00 0.06 0.69 0.17 Ni (μg/g) 0.02 0.02 0.02 0.02 0.02 Pb (μg/g) 0.10 0.14 0.05 0.09