如何运用matlab分析解决实际问题?

如何运用matlab分析解决实际问题?

原始M-K检验可参加另一篇博客-MATLAB实现Mann-Kendall趋势/突变检验。本文主要介绍修正M-K趋势检验。

原始M-K检验中的零假设是数据独立且随机排序,然而,数据中正自相关的存在增加了在实际不存在趋势时检测趋势的概率,反之亦然。

原始数据序列图形如下:

绘图代码如下:

M-K趋势检验代码详见另一篇博客-MATLAB实现Mann-Kendall趋势/突变检验。

主函数代码如下:

突变检验图:

代码如下:

Pettitt突变检验代码详见另一篇博客-MATLAB实现Pettitt突变检验/贝叶斯突变检测等。

主函数代码如下:

调用Pettitt函数代码如下:

统计量示意图如下:(第520个值出发生突变,即1998年3月)

原数据示意图如下: