
原始M-K检验可参加另一篇博客-MATLAB实现Mann-Kendall趋势/突变检验。本文主要介绍修正M-K趋势检验。
原始M-K检验中的零假设是数据独立且随机排序,然而,数据中正自相关的存在增加了在实际不存在趋势时检测趋势的概率,反之亦然。
原始数据序列图形如下:
绘图代码如下:
M-K趋势检验代码详见另一篇博客-MATLAB实现Mann-Kendall趋势/突变检验。
主函数代码如下:
突变检验图:
代码如下:
Pettitt突变检验代码详见另一篇博客-MATLAB实现Pettitt突变检验/贝叶斯突变检测等。
主函数代码如下:
调用Pettitt函数代码如下:
统计量示意图如下:(第520个值出发生突变,即1998年3月)
原数据示意图如下:
