立即注册 登录
气象家园 返回首页

wuyejincao的个人空间 http://bbs.06climate.com/?15920 [收藏] [复制] [分享] [RSS]

留言板

facelist doodle涂鸦板

您需要登录后才可以留言 登录 | 立即注册


東之寂桃 2018-8-21 23:13
你好,我最近在做REOF分析,在论坛里看到了你关于REOF的帖子,不过代码都是Fortran的,下载来看了,实在是看不懂,所以想自己用MATLAB弄一个代码,但是对于REOF分析的计算步骤有些疑惑,看了一些书也没找到答案,所以希望你能给与解答,感激不尽!
       若x是m×n的原数据矩阵,E是经过EOF分析后所获取的空间矩阵,PC是时间系数矩阵,做极大方差旋转之后,a是新因子载荷阵,问题如下:
1)REOF的方差贡献率怎么求?描述:EOF分析结果的方差贡献是通过特征值(λk/∑λ)求得,但是特征值所构成矩阵并没有经过旋转,所以很疑惑如何求REOF各模态的方差贡献率。
2)REOF的时间系数RPC怎么求?
3)通过极大方差正交旋转后的新因子载荷阵 a 是否就是空间模态?还是说需要做一些其他的处理?
        还请不吝赐教,或者给一些可供参考的书籍文献。感激不尽!
nunu18 2018-7-11 11:21
wuyejincao: 貌似一般情况下是做标准化的
reof自带标准化处理了,为什么需要做这个呀,标准化不是为了无量纲处理吗?我的理解是量级和单位不同的要素之间分析才需要处理标准化,为什么reof也需要这么处理,而不像和eof一样,可以选择距平处理?
nunu18 2018-7-10 18:50
wuyejincao: 好像程序内部有开关自动选择的,您可以看下程序,具体我有些遗忘了好多年没用了。对于前处理的距平化,标准化,原始场,分别对应分别后的不同含义,你可以参照下 ...
eof我看出来了,是不需要做,程序自带,但是reof没看出来,所以想进一步确定是否需要处理
nunu18 2018-7-9 20:24
楼主,你好!我最近也在用reof做气候的分区研究,我只是想知道用论坛的reof程序运行之前,是否需要对原始资料进行距平化处理?期待您的答复
nunu18 2018-1-22 19:51
楼主都好忙啊
LOSER 2016-10-24 17:17
请教个问题啊,REOF之后,各个模态的解释方差是不是比较接近,不像EOF减小的那么快?
皮特 2015-5-18 15:20
请问楼主是如何计算REOF时间系数的啊?
xiaowan 2015-4-14 13:32
请问时间序列是如何计算出来的,因为旋转的时候只选取了10个模态,那么时间序列也只有10。如果对所有数据旋转,得到的结果又不太对,请问楼主是如何计算REOF时间系数的啊?谢啦
yangwawa 2014-6-24 17:31
楼主您好,我用你在家园里分享的reof程序进行分析处理,我一共要做两组数据,有一个没有问题,出了结果。另外一个就遇到了math错误,你之前说可能是数据有缺测,我的数据是19站,31年的数据,其中有1个站31年的数据都是0,我看程序输出的相关矩阵a就是在这个站是nan,应该就是这个站的问题。但是这个站不是缺测啊, 只是数据为0,我不大会用fortran程序,我看你给的资料预处理程序就是把缺测值变成了0,另外我用matlab做了相关矩阵a就没有问题,但是我又不会用matlab继续编reof程序。所以希望你能帮助我解决这个问题,谢谢!!
敏芬 2013-11-25 15:12
wuyejincao: 对不起,这个问题我也不清楚,你去问问别人吧
好的。谢谢!
敏芬 2013-11-23 10:49
楼主,我最近想把FNL换成CFSR,可是我在ungrib那里出了问题,我想问链接Vtable的时候是链接很多个还是指链接一个?Vtable要做什么修改,可以大概给我说一下怎么做吗?
lemon_tree 2013-9-6 16:10
楼主,你好,你是不是用旋转EOF对中国降水进行了初步分区?
yy456 2012-12-29 11:51
同问,能否留下qq,有问题请教,谢谢
流云飞袖 2012-12-23 18:59
能否留下qq,有问题请教,谢谢
mofangbao 2012-10-12 16:48
当然可以,你先与我联系下,最近可能会有一批新版主,QQ 937062711

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

返回顶部