二基于主成分分析方法确认重金属污染成因模型

1、二 基于主成分分析方法确认重金属污染成因模型2.1 主成分分析法基本原理主成分分析法是指标筛选最常用的方法之一，它利用降维的思想，根据实际需要，从中可取几个较少的综合指标尽可能多地反映原来指标的信息，把多指标转化为少数几个综合指标，减少变量的个数。故这几个综合指数能包含原信息量的80%以上即可，从而筛选出主要因子。一个主成分不足以代表原来的n个变量，因此需要寻找第二个乃至第三、第四主成分，第二个主成分不应该再包含第一个主成分的信息，统计上的描述就是让这两个主成分的协方差为零，几何上就是这两个主成分的方向正交。具体确定各个主成分的方法如下：设 表示第i 个主成分，i = 1,2,n，可设其中对每

2、一个i ，均有且使得的值达到最大；不仅垂直于，而且使的值达到最大；同时垂直于和，并使)的值达到最大。1）对原始数据进行标准化处理表1 原始数据标准化处理Table 1 Standardization of initial data样本指标x1x2xn1Y11Y12Y1n2Y21Y22Y2n3Y31Y32Y3nmYm1Ym2Ymn标准化为：第j个指标的样本均值：第j个指标的样本方差：经标准化处理后可得标准化矩阵：2)计算相关系数矩阵R计算标准化后的每两个指标间的相关系数，得到相关系数矩阵，即个指标的协方差矩阵。相关系数矩阵：式中 ， ， 是第i 个指标与第j 个指标的相关系数,是一个实对称的相关

3、矩阵3）计算特征值与特征向量：根据前面得到的相关矩阵，由下式可求出相关矩阵的特征值将上式改写成矩阵的形式为就可得到个非负特征值，从而得到对应于特征根的个单位化特征向量，构成一个正交矩阵，记为，则其中中的为第个主分量，为第个分量。4）计算主成分对于个样本中的第个样本，根据，则得到个主成分如下对于全部的个样本，则有即：，整理得：式中 样本主成分，标准化样本。5）主成分选择为了合理选择少数几个主成分来有效地描述原来个指标所构成的一组样本，要引入主成分贡献率的概念及其计算方法。若为相关矩阵的第个特征根，则= 第个主成分的贡献率； = 前个主成分的累计贡献率。 样本前个主成分的累计贡献率表明了前个主成分

4、能够反映原样本信息量的程度。当其达到一定水平时（），说明采用前个主成分来描述原样本所包含的信息量已经可以达到要求。2.2 模型建立与求解 土壤中重金属主成分分析模型基于附表数据cumcm2011A附件_数据，对每个区域的8类不同重金属污染程度进行绘图，其结果如下，图2.1 城区5类区域8种不同重金属污染情况 由图2.1容易发现，生活区Cd、Zn和Hg污染程度较高，分别占到总体污染比例的0.33、0.27和0.20；工业区Hg、Cd成为污染的主要原因，分别占到该区污染比例的0.67、0.16；山区Cd、Zn和Pb造成该区污染的主要原因，分别占总体污染比例的0.35、0.22和0.16；交通区Hg

5、、Cd和Zn，分别占0.43、0.25、0.12；公园绿化区Cu、Cd和Zn是造成污染的重要原因，分别占到总体污染比例的0.38、0.29和0.14.另外，基于原始数据，得到5类不同区域的相关系数矩阵如下表2.2所示表2.2 5类不同区域8种不同重金属污染参数相关系数表Table 2 the correlation coefficient of 8 different heavy mental in 5 areas元素AsCdCrCuHgNiPbZnAs1Cd0.38051Cr0.23850.34921Cu0.53120.49870.37591Hg0.29340.36710.15050.197

6、91Ni0.60530.28260.52710.43410.21141Pb0.45020.80180.41570.50520.34030.30041Zn-001710.34640.41420.23790.24240.33430.32761元素AsCdCrCuHgNiPbZnAs1Cd0.32861Cr0.37960.54101Cu0.15290.56650.91971Hg0.18130.53320.90210.98351Ni0.68970.48870.69830.50280.47911Pb0.39470.82920.67540.66970.61250.57761Zn0.51770.75360.

7、6951062170.59040.63410.73881元素AsCdCrCuHgNiPbZnAs1Cd-.29091Cr0.11320.06651Cu0.52660.08950.36361Hg0.07500.2463-.00590.50511Ni0.07790.04880.94520.3576-.04481Pb-.20540.76580.10730.12230.22620.02801Zn-.17610.60590.62730.25220.16960.62870.58981元素AsCdCrCuHgNiPbZnAs1Cd0.12151Cr0.13940.37301Cu0.09200.32360.8

8、9451Hg-.00410.21110.01190.03211Ni0.22790.35050.86950.28620.03981Pb0.05990.61470.42770.50530.26580.39581Zn0.18800.29420.39520.43160.11840.50310.48211元素AsCdCrCuHgNiPbZnAs1Cd0.35821Cr0.68940.56401Cu0.10690.50030.35731Hg0.17560.05420.02260.13631Ni0.69110.43270.73950.2667-.04751Pb0.26500.59760.39740.7563

9、0.38890.16761Zn0.28530.71170.50900.52110.06310.29780.74761表2.3 不同区域重金属前3个主成分贡献率Table 3 the contribution rate of top 3 heavy mental in different area项目AsCdCrCuHgNiPbZn生活区主成分10.3520.4120.3380.3830.2590.3610.4220.263主成分2-.6070.1610.220-.2310.122-.2380.1060.649主成分3-.009-.4030.4750.0240.4210.505-.3350.25

10、7工业区主成分10.2260.3430.3990.3790.3690.3350.3740.375主成分20.6750.066-.183-.412-.4090.3750.0430.168主成分30.283-.6280.3240.1680.2280.341-.412-.229山区主成分1-.0050.3440.4650.2960.1860.4230.3470.519主成分20.468-.4750.3360.329-.0520.350-.441-.121主成分30.3840.113-.2930.4940.600-.3290.139-.147交通区主成分10.1210.3200.4510.4670.0

11、870.4580.3630.332主成分2-.1150.403-.280-.2250.655-.2820.4200.102主成分30.9410.011-.151-.1970.022-.023-.0570.218公园区主成分10.3210.4100.4090.3430.1030.3360.3400.404主成分2-.4480.081-.3350.3600.278-.4900.4300.225主成分30.340-.226-.003-.1470.8650.0110.094-.235 土壤中重金属的聚类分析聚类树状图可以形象地反映元素间的远近和贴近程度，因而可以有效的揭示元素间的相互联系。由于本题中所有数据均是在