EOF程序的相关问题

学者

各位大神，这个EOF程序里载荷向量所对应的时间系数是哪个，多谢多谢！

[A,B,C]=xlsread('SST.xls');
xx=A;
[m,b]=size(xx);
x=xx';
for i=1:b
    x(i,:)=(x(i,:)-mean(x(i,:)));
end
ss=x;
r=ss*ss'/m;
[v,d]=eig(r);
pc=v'*xx';
D=rot90(d,2);
diagonal=diag(D);
eof=fliplr(v);
pc=flipud(pc);
%检查eof的正交性
ae=eof*eof';
%检查pc的正交性
ap=pc*pc';
%回复到x的距平值
ac=eof*pc;
%以上说明计算数据正确
% 求各模态的方差贡献率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sum_d=sum(diagonal);                 % 总的特征值
e1=(diagonal(1,1)/sum_d)*100%;       %90.2268%
e2=(diagonal(2,1)/sum_d)*100%;       %6.0384%
e3=(diagonal(3,1)/sum_d)*100%;       %2.1832%
e4=(diagonal(4,1)/sum_d)*100%;       %0.4853%
e5=(diagonal(5,1)/sum_d)*100%;       %0.3441%
e6=(diagonal(6,1)/sum_d)*100%;
e7=(diagonal(7,1)/sum_d)*100%;
e8=(diagonal(8,1)/sum_d)*100%;
e9=(diagonal(9,1)/sum_d)*100%;
e10=(diagonal(10,1)/sum_d)*100%;
pc1=pc(1,:);pc2=pc(2,:);pc3=pc(3,:);
f1=eof(:,1);   
eof1=reshape(f1,23,1);     
  f2=eof(:,2);   
  eof2=reshape(f2,23,1);
  f3=eof(:,3);   
  eof3=reshape(f3,23,1);

言深深

没有有matlab分析过，不过程序看下来，似乎应该是pc吧，空间模态是eof?二者乘积回复到原来的场。

叶、空づ

言深深 发表于 2016-9-3 20:50
没有有matlab分析过，不过程序看下来，似乎应该是pc吧，空间模态是eof?二者乘积回复到原来的场。

进行EOF分解的数据一定要是三维数组吗？一维表示经度，一维表示纬度，一维表示时间这样吗？

学者

言深深 发表于 2016-9-3 20:50
没有有matlab分析过，不过程序看下来，似乎应该是pc吧，空间模态是eof?二者乘积回复到原来的场。

是的，我这个程序导入的是二维数，出来的结果用surfer出图可以吧

言深深

学者 发表于 2016-9-4 10:29
是的，我这个程序导入的是二维数，出来的结果用surfer出图可以吧

不清楚呢出来的数据形式，surfer的要求是经纬度，数据，txt格式的即可。

2daidai

谢谢楼主

fcb61723062

谢谢楼主分享

curtainsky

学习了学习了