|
[求助]
求助,关于EOF,REOF分析及作图
我的目的:对某一区域若干年降水特征进行分析
前期工作:因无气象学以及matlab基础,搜了很多资料,学习修改了一些代码,但基础问题仍不明白。作图可能会用到arcgis之类软件,这部分内容可以求助课题组其他人,但EOF和REOF的工作还是要我自己完成的。
求助内容:1,分析数据的原始格式是什么(行为时间,列为地点?时间是以日月还是以年为单位?)。
2,我的EOF代码是否有问题?
%% 特征向量及贡献率
[EOF,E] = eig(data*data');
EOF = fliplr(EOF); % 特征向量
E = fliplr(flipud(E)); % 特征根
PC = EOF*data;% 时间函数?时间场?主成分?
PC1 = PC(:,1)/sqrt(E(1,1)); % xdata*xdata’的第一特征向量
PC2 = PC(:,2)/sqrt(E(2,2)); % xdata*xdata’的第二特征向量
PC3 = PC(:,3)/sqrt(E(3,3)); % xdata*xdata’的第三特征向量
diagonal=diag(E);
sum_d=sum(diagonal); % 总的特征值
e1=(diagonal(1,1)/sum_d)*100;
e2=(diagonal(2,1)/sum_d)*100;
e3=(diagonal(3,1)/sum_d)*100;
%% REOF
REOF = varimax(EOF,E);% ?
3,不同模态的提取,以及作图?
真诚求助,祝大家多发paper |
|