Matlab利用距平变量对温度场进行主成分分析

作业记录
题目：东海 30度 N，125度E， 30N,130度E，30度N，135度E 三网格 1976.12-1977.11的月平均水温如下
试将温度场进行主成分分析

解答：

1. 利用距平变量计算主成分思路
   主成分分析是一种通过降维技术把多个变量化为少数几个主成分（新的综合变量）的多元统计方法。这些主成分能够反映原始变量的大部分信息，通常表示为原始变量的线性组合，且要求这些主成分各自所包含的信息互不重叠，即主成分之间互不相关。求取主成分的具体步骤：
   第一步，整理资料，根据需要进行中心化或者标准化处理，建立资料矩阵X_(m×n)，有m个原始变量（空间点），n个样本；
   第二步，求取X(m×n)的m阶协方差矩阵S；
   第三步，求矩阵S的特征值和特征向量。其中特征值从大到小排列λ1≥λ2≥⋯≥λm，对应特征向量构成矩阵从左向右排列B=(b1 b2 ⋯ bm)；
   第四步，计算主成分的方差贡献以及累积方差贡献，根据贡献值大小取前若干个主成分，计算主成分（新的综合变量）的数值。

2. 程序实现

clear all;clc
%--------利用距平变量计算主成分---------%
%%%%%原始数据 %%%%%
lat=30;
lon=[125,130,135];
t=1:12;
sst=[17.7,15.0,14.4,13.2,15.0,17.0,22.5,27.9,28.0,27.2,24.6,22.2;
     21.1,21.1,19.6,21.4,22.2,24.2,25.8,28.8,28.5,27.5,26.1,24.3;
     21.4,20.0,18.9,19.7,21.5,23.4,25.0,28.3,28.5,27.8,25.9,24.3];
%%%%%距平资料阵的计算%%%%
SST(1,:)=sst(1,:)-mean(sst(1,:));%sst中各个分量扣除自身的平均值
SST(2,:)=sst(2,:)-mean(sst(2,:));
SST(3,:)=sst(3,:)-mean(sst(3,:));

S=SST*SST';%协方差阵
[L,namta]=eig(S);%求矩阵S的特征值L和特征向量namta
d=diag(namta);%取出特征值矩阵列向量（提取出每一主成分的贡献率）
namta_des=sort(d,'descend'); %将特征值矩阵列向量按照贡献率的大小descend降序排列
v=fliplr(L);       %将特征向量L按照特征值大小从上向下排列
%%%%第k个主成分，k=1,2,3%%%%%
  for k=1:3
       perc(k)=namta_des(k)/sum(d); %第k个主成分的方差贡献率
  end
pc1=v(:,1)'*SST;%%第一主成分
pc2=v(:,2)'*SST;%%第二主成分
pc3=v(:,3)'*SST;%%第三主成分

lon=125:5:135;
t=1:12;
[T,Lon] = meshgrid(t,lon);
figure(1);
contourf(T,Lon,SST,'ShowText','on')
colorbar 
hold on
axis tight;
xlabel('t=12','FontSize',14);ylabel('lon','FontSize',14);
title('距平资料阵SST','FontSize',20); 

figure(2)
subplot(2,1,1)
plot(v(:,1),lon,'*-')
hold on
plot(v(:,2),lon,'o-')
hold on 
plot(v(:,3),lon,'.-')
hold on
legend('l1','l2','l3');
axis tight;
xlabel('lk(k=1,2,3)','FontSize',14);ylabel('lon','FontSize',14);
title('第k主成分的系数向量','FontSize',20);
hold on

subplot(2,1,2)
plot(t,pc1,'*-')
hold on
plot(t,pc2,'o-')
hold on 
plot(t,pc3,'.-')
hold on
legend('pc1','pc2','pc3');
axis tight;
xlabel('t=12','FontSize',14);ylabel('yk(k=1,2,3)','FontSize',14);
title('第k主成分的时间序列(PC)','FontSize',20);
hold on

结果如下：
把X绘成填色和等值线图如下所示：

前三个主成分的系数向量（3行1列）以及时间序列（1行12列）图如下所示：

各个主成分的方差贡献率如下：

程序一如既往幼稚简单，说不定还有很多错误，如有赐教，不甚感激~
see you next time~ bye-bye~
立冬快乐，盼诸君安~ Over！