气象统计应用15多个变量的主成分幻灯片课件.pptVIP

气象统计应用 赵传湖 海洋环境学院 chzhao@ouc.edu.cn 多个变量的主成分 主成分的表示 研究对象为一气象要素场，场中有P个空间点，样本容量为n，记p个空间点上要素为x1,x2,…,xp，其观测值为xki(k=1,2,…,p;i=1,2,…,n)由这p个变量组成新变量 若y满足方差极大的要求，则y称为原p个变量的主成分（主分量），若进一步有 则有 主成分的导出 要求y具有极大的方差 又 则 即 在 条件下的极值问题 则转化为 求极值的问题 即有 整理得 要使v有非零解，必须有 问题转化为求矩阵S的特征值和特征向量 主成分的计算步骤 根据资料阵计算协方差阵S或相关阵R 求其特征值和特征向量 计算给定的累积贡献，确定特征向量和主成分的个数 计算主成分 P148~149例子 1.利用距平资料计算主成分 x x = -0.5600 -0.2400 -0.8300 1.5400 0.1500 0.6600 1.0400 0.2500 0.0600 1.3400 0.4500 -0.5300 -0.0600 -0.6400 0.4600 0.4400 -1.0400 -1.8300 -0.3600 -0.1400 0.5600 0.1400 -0.8400 0.8600 0.4400 -0.5400 0.4600 0.5400 0.2500 0.6600 0.0400 0.4500 -0.8300 -2.5600 -1.4400 -1.1300 -2.0600 -0.5400 0.5600 -1.8600 -0.4400 0.2600 -0.3600 -0.2400 -0.6300 1.3400 1.1500 0.3600 -0.7600 -0.0400 -0.1300 0.6400 0.9500 -0.1300 0.0400 -0.6400 -0.4300 1.1400 1.1500 -1.5300 -0.0600 1.2500 1.0600 0.6400 -0.4400 0.4600 -0.4600 -1.3400 -1.1300 1.3400 1.6500 1.4600 0.6400 0.6500 1.6600 0.1400 1.3500 -0.1300 -0.7600 -0.7400 0.1600 -1.2600 0.7500 -0.5300 0.3400 0.4500 0.8600 -0.6600 -1.5400 -0.8300 （1）原始数据 （2）计算协方差阵 s=cov(x,1); s = 1.0144 0.4563 0.1934 0.4563 0.7299 0.2566 0.1934 0.2566 0.7118 （3）计算特征值和特征向量 [v,d]=eig(s); v = 0.4697 0.5015 0.7266 -0.8157 -0.0683 0.5744 0.3377 -0.8625 0.3770 d = 0.3610 0 0 0 0.6197 0 0 0 1.4755 （4）特征值排序 [newd,i]=sort(diag(d),descend); newd = 1.4755 0.6197 0.3610 I = 3 2 1 特征向量排序 newv=v(:,I); newv = 0.7266 0.5015 0.4697 0.5744 -0.0683 -0.8157 0.3770 -0.8625 0.3377 （5）计算累积方差贡献 G=newd/sum(newd);%计算方差贡献 Gs=0; for i=1:length(G) Gs(i)=sum(G(1:i));%计算累积方差贡献 if Gs(i)0.8 break;