主成分分析分析法.doc

第四节主成分分析方法地理环境是多要素的复杂系统,在我们进行地理系统分析时,多变量问题是经常会遇到的。变量太多,无疑会增加分析问题的难度与复杂性,而且在许多实际问题中,多个变量之间是具有一定的相关关系的。因此,我们就会很自然地想到,能否在各个变量之间相关关系研究的基础上,用较少的新变量代替原来较多的变量,而且使这些较少的新变量尽可能多地保留原来较多的变量所反映的信息?事实上,这种想法是可以实现的,本节拟介绍的主成分分析方法就是综合处理这种问题的一种强有力的方法。第一节主成分分析方法的原理主成分分析是把原来多个变量化为少数几个综合指标的一种统计分析方法, 从数学角度来看,这是一种降维处理技术。假定有 n个地理样本,每个样本共有 p个变量描述,这样就构成了一个 n×p阶的地理数据矩阵: 如何从这么多变量的数据中抓住地理事物的内在规律性呢?要解决这一问题,自然要在 p维空间中加以考察,这是比较麻烦的。为了克服这一困难,就需要进行降维处理,即用较少的几个综合指标来代替原来较多的变量指标,而且使这些较少的综合指标既能尽量多地反映原来较多指标所反映的信息,同时它们之间又是彼此独立的。那么,这些综合指标(即新变量)应如何选取呢?显然,其最简单的形式就是取原来变量指标的线性组合,适当调整组合系数,使新的变量指标之间相互独立且代表性最好。如果记原来的变量指标为 x 1,x 2,…,x p,它们的综合指标——新变量指标为x 1,x 2,…,zm(m≤p)。则在( 2)式中,系数 l ij由下列原则来决定: (1)z i与z j(i≠j;i,j=1 ,2,…,m)相互无关; (2)z 1是x 1,x 2,…,x p的一切线性组合中方差最大者; z 2是与 z 1不相关的 x 1,x 2,…,x p的所有线性组合中方差最大者;……;z m是与 z 1,z 2,…… z m-1都不相关的 x 1,x 2,…,x p的所有线性组合中方差最大者。这样决定的新变量指标 z 1,z 2,…,zm分别称为原变量指标 x,x 2,…,x p 的第一,第二,…,第m主成分。其中,z 1在总方差中占的比例最大,z 2,z 3,…, z m的方差依次递减。在实际问题的分析中,常挑选前几个最大的主成分,这样既减少了变量的数目,又抓住了主要矛盾,简化了变量之间的关系。从以上分析可以看出,找主成分就是确定原来变量 x j(j=1 ,2,…,p)在诸主成分 z i(i=1 ,2,…,m)上的载荷 l ij(i=1 ,2,…,m;j=1 ,2,…,p),从数学上容易知道,它们分别是 x 1,x 2,…,x p的相关矩阵的 m个较大的特征值所对应的特征向量。第二节主成分分析的解法主成分分析的计算步骤通过上述主成分分析的基本原理的介绍,我们可以把主成分分析计算步骤归纳如下: (1)计算相关系数矩阵在公式( 3)中, r ij(i,j=1 ,2,…,p)为原来变量 x i与x j的相关系数,其计算公式为因为 R是实对称矩阵(即r ij=r ji),所以只需计算其上三角元素或下三角元素即可。(2)计算特征值与特征向量首先解特征方程| λI-R |=0求出特征值λ(i=1 ,2,…,p),并使其按大小顺序排列,即λ 1≥λ 2≥…,≥λ p≥0;然后分别求出对应于特征值λ i的特征向量 e i(i=1 ,2,…,p