请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4353|回复: 2

请教MATLAB做BFAST突变分析

[复制链接]

新浪微博达人勋

发表于 2018-1-1 19:33:10 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
求大神指导
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-10-25 19:15:59 | 显示全部楼层
请问楼主问题解决了吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-5 06:45:06 | 显示全部楼层
本帖最后由 zhaokg 于 2023-12-29 05:06 编辑

如果有人需要一个贝叶斯算法的话,我开发了一个类似bfast的时间序列趋势,断点/突变分析,和时间序列分解算法叫做BEAST. 这里有更多介绍: http://github.com/zhaokg/Rbeast.

BEAST可以在R和Matlab中调用:

in R
  1. install.packages("Rbeast")
  2. data(Nile)                       #  annual streamflow of the Nile River   
  3. out = beast(Nile, season='none') #  'none': trend-only data without seasonlaity   
  4. print(out)                  
  5. plot(out)
复制代码


in Matlab:
  1. eval(webread('http://b.link/rbeast',weboptions('cert','')))  % install BEAST to your local machine
  2. load('Nile.mat')                                                           % annual streamflow of the Nile River startin from year 1871
  3. out = beast(Nile, 'season', 'none','start', 1871)               % trend-only data without seasonality
  4. printbeast(out)
  5. plotbeast(out)
复制代码


Matlab中另外一个例子
  1. load('googletrend.mat')    % Monthly Google search trend data of 'beach' since Jan 2004
  2. o = beast( beach )         % Apply BEAST to the 'beach time series: beach is a data vector only; the time
  3.                                     % info can be supplied using the start and deltat
  4.                                     % keywords, as in the next commented line
  5. %o = beast( beach,'start',[2004,1],'deltat',1/12 )
  6. printbeast(o)           % print the changepoints detected
  7. plotbeast(o)            % plot the results: o.season.Y and o.trend.Y are the seasonal and trend compoents
复制代码

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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