[原创]Matlab气候突变检测程序合集（更新了）

sfhstcn2

漠小北 发表于 2019-2-15 15:10
前辈加楼主：

         您好！

不敢当不敢当。是这样子，是关于“突变检验”和“突变点”含义理解的问题，“突变点”并不是“序列趋势变化点f'(x)=0甚至f''(x)=0”，而是指的“序列出现了某种明显均值或者方差变化的地方”并且这种变化通过了检验的显著水平，因而“突变点”可能不是一点，甚至可能是一段。Pettitt检验你给的这个序列结果显示是2005年，P=0.0042，显著的；你做的t检验结果中2005年这点突变也是显著的。当然t检验结果显示后面那一段都是显著的，而Pettitt检验只有2005年这一点，这是由于两者算法本身的差异，Pettitt检验只给出了检验所计算的统计量中最大值所在处，并不能说明其它点就不是显著突变点。

漠小北

sfhstcn2 发表于 2019-2-15 23:55
不敢当不敢当。是这样子，是关于“突变检验”和“突变点”含义理解的问题，“突变点”并不是“序列趋势变 ...

哇，谢谢前辈解答，那如果我是想找突变点，就是变化趋势方向发生转变（由增加变为减少或者由减少变为增加）的那个点，请问有什么检验可以做吗？再次感谢您啦。。。

sfhstcn2

漠小北 发表于 2019-2-17 13:21
哇，谢谢前辈解答，那如果我是想找突变点，就是变化趋势方向发生转变（由增加变为减少或者由减少变为增加 ...

可以对序列做个低通滤波，然后做一阶差分（相当于求一阶导），找到一阶差分后序列的等于或者最接近于0的点（f'(x)=0极值点）

漠小北

sfhstcn2 发表于 2019-2-17 13:58
可以对序列做个低通滤波，然后做一阶差分（相当于求一阶导），找到一阶差分后序列的等于或者最接近于0的 ...

好的，非常感谢前辈的解答{:eb500:}{:eb500:}

luobetty

谢谢

小小南瓜

楼主，这是我用MK突变检验做出来的图，显示突变点在2013年左右，但是我用pettitt出来的结果是2009年，且2009年不显著。我做累积距平发现序列在2000年、2003年、2009年有较大的增长幅度，这样看好像2009年是突变点更合理一点。MK突变检验有什么适用条件吗，为什么出来几种方法出来的结果不一致，我处理了好几组数据都发现用累积距平观察出来的突变点和MK分析出来的不一样，这是什么原因呐……希望楼主帮忙解答一下啦

小小南瓜

这个是累积距平图

zyhsdau

新人想下载还没有金钱，正在想办法搞金钱

sfhstcn2

小小南瓜 发表于 2019-3-7 16:54
楼主，这是我用MK突变检验做出来的图，显示突变点在2013年左右，但是我用pettitt出来的结果是2009年，且200 ...

MK最适合检验的就是在某个时刻突然出现持续性增加或减少的序列，MK检验挺严格的，感觉你这个序列在09年后增加持续得还不够久，不过还是检测出了2013年，可能即使09年平均值虽然开始变化，但还不能确定这时就开始增长了。建议还是MK检验的结果有说服力些。

zyhsdau

好东西好东西