主成分分析分析法

《主成分分析分析法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、第四节主成分分析方法地理环境是多要素的复杂系统，在我们进行地理系统分析时，多变量问题是经常会遇到的。变量太多，无疑会增加分析问题的难度与复杂性，而且在许多实际问题中，多个变量之间是具有一定的相关关系的。因此，我们就会很自然地想到，能否在各个变量之间相关关系研究的基础上，用较少的新变量代替原来较多的变量，而且使这些较少的新变量尽可能多地保留原來较多的变量所反映的信息？事实上，这种想法是可以实现的，木节拟介绍的主成分分析方法就是综合处理这种问题的一种强冇力的方法。第一节主成分分析方法的原理主成分分析是把原来多个变量化为少数几个综合指标的一种

2、统计分析方法，从数学角度來看，这是一种降维处理技术。假定有n个地理样本，每个样本共有P个变量描述，这样就构成了一个nXp阶的地理数据矩阵：如何从这么多变量的数据屮抓住地理事物的内在规律性呢？要解决这一问题，自然要在p维空间中加以考察，这是比较麻烦的。为了克服这一困难，就需要进行降维处理，即用较少的几个综合指标来代替原来较多的变量指标，而且使这些较少的综合指标既能尽量多地反映原來较多指标所反映的信息，同时它们之间又是彼此独立的。那么，这些综合指标(即新变量)应如何选取呢？显然，其最简单的形式就是取原来变量指标的线性组合，适当调整组合系数，

3、使新的变量指标之间相互独立且代表性最好。如果记原来的变量指标为Xi,x2,・・・，Xp,它们的综合指标新变量指标为Xi,X2,…，zm(mWp)。贝ij在（2）式屮，系数1订由下列原则来决定:(l)zi与Zj(iHj；i,j=l,2,…，m)相互无关;（2）z】是x”X2,…，Xp的一切线性组合中方差最大者；Z2是与Zi不相关的X1,X2,…,Xp的所有线性组合中方差最大者;；Z■是与Z

4、,Z2,Zl都不相关的X”X?,…，Xp的所冇线性组合中方差最大者。这样决定的新变量指标Zi,Z2,…，zm分别称为原变量指标X”X2,…，xP的第一

5、，第二，…，第m主成分。其中，zi在总方差屮占的比例最大，Z2,z3,…,的方差依次递减。在实际问题的分析中，常挑选前儿个最大的主成分，这样既减少了变量的数口，又抓住了主要矛盾，简化了变量之间的关系。从以上分析可以看出，找主成分就是确定原來变量x」（j二1,2,…，p）在诸主成分Zi（i=l,2,…，m）上的载荷lij（i=l,2,…,m；j=l,2,…，p）,从数学上容易知道，它们分别是X"X2,・・・，x,的相关矩阵的m个较大的特征值所对应的特征向量。第二节主成分分析的解法主成分分析的计算步骤通过上述主成分分析的基本原理的介绍，我们

6、可以把主成分分析计算步骤归纳如下：（1）计算相关系数矩阵C3)kJ.在公式（3）中，ru（i,计算公式为,j=l,2,…，p）为原来变量Xi与&的相关系数，其（4）因为R是实对称矩阵（即甘门），所以只需计算其上三角元素或卜三角元素即可。（2）计算特征值与特征向量首先解特征方程丨入1-RI二0求出特征值入「（i二1,2,…，p）,并使其按大小顺序排列，即22…、2入卩20；然后分别求出对应于特征值入i的特征向量ei（i二1,2,…，p）o（3）计算主成分贡献率及累计贡献率主咖I则率彳戈汗压叫乙胁文?r■龙sk-lUlr-一般取累计贡献率达

7、85-95%的特征值X,入2,・・・，♦所对应的第一，第二，，第m（mWp）个主成分。（4）计算主成分载荷p（%叫）=Jy金敬afc=l>2>■■■>p）Cs）由此可以进一步计算主成分得分：第三节主成分分析应用实例主成分分析实例对于某区域地貌-水文系统，其57个流域盆地的九项地理要素：XI为流域盆地总高度（m）x2为流域盆地山口的海拔高度（m）,x：j为流域盆地周长（m）,X；为河道总长度（km）,xs为河表2-14某57个流域盆地地理要素数据1TOOM901.1042.481302.TBS31.8302189144602.1B64.3

8、M306.8339T.038992666261.6002.0ft)383.04221.126461641002.T606.32011T4.84490.198661309901.1422.080326.10026.T2B0tsn8M00.19&10.210fB4.290M.961T221084168.T»16.000654.6003S.665861fi10401.3001.290133.60022.2109119263B88.44T30.0)0atfl.(DO29.12261016«6.LT4LL.38382<0T923.31196086

9、202.8906.810623.6B02T.24712860MSOT.«DT.190304.90011.02413L23T699T2.O4B2.999282.T3)29.61914669TOO4.12022