新手做的eof分析，有些小问题，望大家帮忙！

泰达巴神

现有一个经过标准化处理后的资料，格点数m=1769，时间长度n=53，由于m>n，因此做eof时需要进行时空转换。以下为我的程序：***********************************************

X = dlmread('D:、easa.txt');
[m,n]=size(X); %求空间点数和时间
[VR,DR] = eigs(X'*X/n);  %计算特征值和特征向量
   DR = diag(DR); %特征值
   DR = DR(:); %特征值按列排列
   VN = X*VR; %时空转换
   T  = VN'*X; %时间序列
   VN = VN.*repmat(sqrt(DL'),m,1); %特征向量标准化
   T  = T./repmat(sqrt(DL),1,n);  %时间序列标准化
   D100  = DR*100./trace(X*X'/n); %方差贡献率


************************************************
但是计算得到的结果不正确，时间序列画出来的图也有问题。
巧合的是我用论坛中有人分享过的一个叫做“一键eof”的软件的时候得到的时间序列图也是这种正弦曲线，但是用fortran程序做出来的结果就是正确的。
而且还有一个问题：为什么我计算出来的特征值和特征向量都取了6组？我并没有进行任何设置，是默认的吗？
希望大家帮我看看，找出原因，感激不尽！

yangjunhua

楼主很厉害啊

dreamdong

楼主问题解决了吗？

泰达巴神

dreamdong 发表于 2015-12-30 14:44
楼主问题解决了吗？

嗯，EOF已经没问题了

dreamdong

泰达巴神 发表于 2015-12-30 22:05
嗯，EOF已经没问题了

好吧，我最近也在搞这个，如果有问题，还想问问楼主呢。