基于预白化MK趋势和突变分析MATLAB代码

Jude · 发表于 2018-10-17 10:39:16

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

MK趋势突变检验要求时间序列是独立的，那么对于有自相关性时间序列，应该做预白化处理，用预白化处理后的数据进行分析。本人参考了一篇用一阶自相关系数进行预白化的文章。用MKdata中的数据，先用MATLAB中的autocorr函数进行自相关分析，发现，存在自相关性，如图所示。按照文献中公式（4）方法去除自相关性结果如另一张图所示。
现将代码和文献上传，压缩包中prewhitening_MK为主程序，可以输出突变交点图像（如图所示）以及MK趋势检验统计Z值；UFK是参考网上的MK突变分析代码；S是用于MK趋势分析计算统计值Z值的；prewhiten是预白化程序。
MK突变结果示意的横纵坐标大小有点问题，是因为，我用了比较大号的字体，源程序中可以改小。
现奉源码，以飨读者。

江湖小鱼儿 · 发表于 2018-10-17 11:37:49

非常感谢楼主分享，最近正在学习MK相关的改进方法，希望能跟楼主请教

看云淡风清 · 发表于 2018-10-17 11:54:44

感谢楼主分享，非常有用！

Jude · 发表于 2018-10-19 10:43:50

江湖小鱼儿发表于 2018-10-17 11:37
非常感谢楼主分享，最近正在学习MK相关的改进方法，希望能跟楼主请教

互勉互勉，这个是我去年做的情况，本想投文章，但是被退稿了，所以就发论坛里供大家讨论讨论吧。

sciren · 发表于 2018-11-4 18:57:19

谢谢楼主感谢

我是一只付球球 · 发表于 2018-11-17 22:42:44

感谢楼主！

我才不是马大哈 · 发表于 2018-11-19 21:48:24

非常感谢分享

sciren · 发表于 2018-12-6 17:54:09

大神为什么我把MKdata,dat换成30行一列的就说我没定义mkdata啊

Jude · 发表于 2018-12-7 13:49:26

sciren 发表于 2018-12-6 17:54
大神为什么我把MKdata,dat换成30行一列的就说我没定义mkdata啊

这个问题好抽象，我也不知怎么回答了

sciren · 发表于 2018-12-8 14:42:06

Jude 发表于 2018-12-7 13:49
这个问题好抽象，我也不知怎么回答了

大神我昨天重新弄了可以了当时更改的时候不知道为啥就是跑不出来

		自动登录	找回密码
密码			立即注册

[源程序] 基于预白化MK趋势和突变分析MATLAB代码

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