24小时热门版块排行榜    

查看: 2303  |  回复: 4

東之寂桃

新虫 (初入文坛)

[求助] 关于REOF的一些问题已有1人参与

本人最近在做REOF分析,在网上有很多REOF的代码都是Fortran的,但是看不懂Fortran的代码,所以想自己用MATLAB弄一个代码,但是对于REOF分析的计算步骤有些疑惑,看了一些书也没找到答案,所以希望在这里得到解答,感激不尽!
       x是m×n的原数据矩阵,若E是经过EOF分析后所获取的空间矩阵,PC是时间系数矩阵,做极大方差旋转之后,a是新因子载荷阵,问题如下:
1)REOF的方差贡献率怎么求?描述:EOF分析结果的方差贡献是通过特征值(λk/∑λ)求得,但是特征值所构成矩阵并没有经过旋转,所以很疑惑如何求REOF各模态的方差贡献率。
2)REOF的时间系数RPC怎么求?
3)旋转后的新因子载荷阵 a 是否就是空间模态?还是说需要做一些其他的处理,例如求出a 的协方差矩阵的特征值与特征向量,然后按照EOF的方法处理得出空间系数和时间系数?
        还望各位虫友不吝赐教!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shen159876

新虫 (小有名气)

【答案】应助回帖

http://www.ncl.ucar.edu/Document/Functions/Built-in/eofunc_varimax.shtml

这是NCL中使用的REOF的函数,你可以搜索了看一下脚本如何写的,应该对你的思路有启发
2楼2018-10-26 15:39:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

燕子小小麦

新虫 (正式写手)

楼主你好请问能分享一下你的程序吗

发自小木虫Android客户端
3楼2019-08-13 21:42:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

young90

金虫 (正式写手)

一直困惑,EOF不就是matlab中的PCA吗

发自小木虫IOS客户端
4楼2019-08-15 09:42:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li123456

新虫 (小有名气)

你好,请问这个问题解决了吗??我用matlab里面的pca算的eof但是用rotatefactors算出reof之后,怎么算方差贡献率呢
5楼2022-04-29 21:49:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 東之寂桃 的主题更新
信息提示
请填处理意见