爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: 简单幸福

[程序设计] REOF 的matlab代码

  [复制链接]

新浪微博达人勋

发表于 2013-9-6 12:55:36 | 显示全部楼层
谢谢楼主的分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-9-16 16:58:17 | 显示全部楼层
谢谢lz~ 我来验证一下~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-12-18 19:03:14 | 显示全部楼层
感谢分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-4 20:10:30 | 显示全部楼层
能不能用rotatefactors函数?默认就是用varimax旋转的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-4 20:10:50 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-4 20:17:25 | 显示全部楼层
本帖最后由 zosuper 于 2014-1-4 20:18 编辑

在[p, nc] = size(x);这一句之前给x一个值,例如eof分解出来空间场的前4个场,然后计算旋转后的结果。
与直接用[xr,T] = rotatefactors(x,'method','varimax','Normalize','off');计算出来的一样呀!
也就是说如果有了x,直接调用这个函数就可以实现帖子给出的功能了。xr就是旋转后的结果。如果y是eof分解得到的空间场,用前4个特征向量的话,[yr,T] = rotatefactors(y(:,1:4),'method','varimax','Normalize','off');即可
密码修改失败请联系微信:mofangbao
回复 支持 2 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-1-9 19:13:55 | 显示全部楼层

和fortran算出来的不太一样 。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-10 09:12:45 | 显示全部楼层
Brick_peng 发表于 2014-1-9 19:13
和fortran算出来的不太一样 。。。

趋势一样就可以了呀
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-18 17:41:59 | 显示全部楼层
好东西啊,真心感谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-8 22:46:23 | 显示全部楼层
楼主~~你给的这个代码是用来求REOF的最后一步极大方差正交旋转过程的吧~~前面的步骤还得处理好对吗
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

快速回复 返回顶部 返回列表