做eof时碰到了一个很让我困惑的问题！求解答！

泰达巴神

言深深 发表于 2014-11-24 20:28
可以帮你看！但是我能说我不懂matlab嚒

！！！怎么会这样...
那我想问下，我有一个fortran的eof程序，算出来的结果时间序列的趋势是正确的，但数值不对，太大了。时间序列和特征向量要分别除以和乘以对应的特征值，将其标准化。我感觉我程序算出来的特征值太大了

言深深

泰达巴神 发表于 2014-11-24 21:10
！！！怎么会这样...
那我想问下，我有一个fortran的eof程序，算出来的结果时间序列的趋势是正确的，但 ...

特征值大是可能的啊，特征值归一化之后才是解释方差。

liu_senfeng

泰达巴神 发表于 2014-11-24 15:47
X = dlmread('D:\finalcca\sa\ea\easa.txt');
[m,n]=size(X); %求空间点数和时间
[VR,DR] = eigs(X' ...

能否说明下你的X是否已经中心化，是m-by-n还是n-by-m，说详细点嘛，不然其他人想看都看不懂呢

liu_senfeng

泰达巴神 发表于 2014-11-24 15:47
X = dlmread('D:\finalcca\sa\ea\easa.txt');
[m,n]=size(X); %求空间点数和时间
[VR,DR] = eigs(X' ...

我先假设你的 [m,n]=size(X); %求空间点数和时间   跟读取是X是匹配的，则错误的有 第3行 改为  [VR,DR] = eigs(X'*X/m);  第5行改为：DR = DR(:).*m./n;    第8,9行的  DN  是啥？？？ 不对。向量归一化，应该是对每个k的向量逐个除以他的内积，时间系数的求法应该在特征向量归一化之后，T = V_normalized'*X   

泰达巴神

liu_senfeng 发表于 2014-11-25 12:08
我先假设你的 [m,n]=size(X); %求空间点数和时间   跟读取是X是匹配的，则错误的有 第3行 改为  [VR,DR]  ...

第三行和第五行改过以后，就相当于在第三行直接除以n，我本来是[VR,DR] = eigs(X'*X/n);这么写的，但是感觉数值有问题，在网上看了一些教程有的不用除n，就去掉了。按你的修改过后，就和我之前的结果一样了。第8.9行打错了，应该是DR.
我的资料是(61*29,53)，若可能的话我把资料传给你，也想跟你交流一下，谢谢！

泰达巴神

泰达巴神 发表于 2014-11-25 15:50
第三行和第五行改过以后，就相当于在第三行直接除以n，我本来是[VR,DR] = eigs(X'*X/n);这么写的，但是感 ...

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

泰达巴神

泰达巴神 发表于 2014-11-25 15:55
X = dlmread('D:\finalcca\sa\ea\easa.txt');
[m,n]=size(X); %求空间点数和时间
[VR,DR] = eigs(X'*X ...

这是我改过的程序，时间系数画出来还是不对...

泰达巴神

言深深 发表于 2014-11-25 08:40
特征值大是可能的啊，特征值归一化之后才是解释方差。

是呢，我fortran程序里有相应的步骤，但还是不对，呜呜

scwjx520

一直想做EOF，到现在还没开始动手，，赶紧整明白了和我们分享下。

超级努力

谢谢分享~~~~~~~~~~~~~~~