主成分分析分析法.docx

第四节 主成分分析方法地理环境是多要素的复杂系统,在我们进行地理系统分析时,多变量问题是经常会遇到的。变量太多,无疑会增加分析问题的难度与复杂性,而且在许多实际问题中,多个变量之间是具有一定的相关关系的。因此,我们就会很自然地想到,能否在各个变量之间相关关系研究的基础上,用较少的新变量代替原来较多的变量,而且使这些较少的新变量尽可能多地保留原来较多的变量所反映的信息?事实上,这种想法是可以实现的,本节拟介绍的主成分分析方法就是综合处理这种问题的一种强有力的方法。第一节主成分分析方法的原理主成分分析是把原来多个变量化为少数几个综合指标的一种统计分析方法,从数学角度来看,这是一种降维处理技术。假定有n个地理样本,每个样本共有p个变量描述,这样就构成了一个nxp阶的地理数据矩阵:如何从这么多变量的数据中抓住地理事物的内在规律性呢?要解决这一问题,自然要在p维空间中加以考察,这是比较麻烦的。为了克服这一困难,就需要进行降维处理,即用较少的几个综合指标来代替原来较多的变量指标, 而且使这些较少的综合指标既能尽量多地反映原来较多指标所反映的信息, 同时它们之间又是彼此独立的。那么,这些综合指标(即新变量)应如何选取呢?显然,其最简单的形式就是取原来变量指标的线性组合,适当调整组合系数,使新的变量指标之间相互独立且代表性最好。X2,…,Xp,它们的综合指标新变量指标如果记原来的变量指标为Xi,为Xi,X2,…,zm(mep)。贝U坷"】內+G衍++l]p%在(2)式中,系数lj由下列原则来决定:乙与zj(i工j;i,j=1,2,…,m)相互无关;zi是xi,X2,…,Xp的一切线性组合中方差最大者;Z2是与zi不相关的Xi,X2,…,Xp的所有线性组合中方差最大者;;Zm是与Zi,乙, Zm-1都不相关的Xi,X2,…,Xp的所有线性组合中方差最大者。这样决定的新变量指标zi,Z2,…,zm分别称为原变量指标Xi,X2,…,xp的第一,第二,…,第m主成分。其中,乙在总方差中占的比例最大,z2,Z3,…,zm的方差依次递减。在实际问题的分析中,常挑选前几个最大的主成分,这样既减少了变量的数目,又抓住了主要矛盾,简化了变量之间的关系。从以上分析可以看出,找主成分就是确定原来变量 Xj(j=i,2,…,p)在诸主成分Zi(i=i,2,…,m)上的载荷lj(i=i,2,…,mj=i,2,…,p),从数学上容易知道,它们分别是Xi,X2,…,xp的相关矩阵的m个较大的特征值所对应的特征向量。第二节主成分分析的解法主成分分析的计算步骤通过上述主成分分析的基本原理的介绍,我们可以把主成分分析计算步骤归纳如下:计算相关系数矩阵ipJ在公式(3)中,九(i,j=i,2,…,p)为原来变量Xi与Xj的相关系数,其计算公式为因为R是实对称矩阵(即rj=rj),所以只需计算其上三角元素或下三角元素即可。计算特征值与特征向量首先解特征方程丨入I-R|=0求出特征值入「(i=1,2,…,p),并使其按大小顺序排列,即入i>X2》…,》入p>0;然后分别求出对应于特征值入i的特征向量e(i=1,2,…,p)。计算主成分贡献率及累计贡献率P > P主咸分颂献率72”「p),累计贡献率工丫“工化。k-l k-l b・「般取累计贡献率达85-95%的特征值X1,X2,…,Xm所对应的第一,第•…,