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

desertsky

暖气团 发表于 2018-3-9 20:45
楼主，为什么我做的Pettitt检测出来的突变点在最后一年，并且还通过了0.5的显著性检验

我也是呢 ：（

desertsky

sfhstcn2 发表于 2018-3-14 13:11
看了下，这个方法只寻找统计量“秩和”的最大值作为突变点，对于一些趋势变化不明显的序列确实是有可能出 ...

感谢楼主的无私奉献，有关pettitt的检验，发现用如下程序是可以的，但用您的程序是出现了如前两位所说的问题，即最后位是突变点而且通过了显著性检验，请明察，菜鸟小白，多拍砖 ：）
https://www.mathworks.com/matlabcentral/fileexchange/60973-pettitt-change-point-test-for-univariate-time-series-data

青风海魅

sfhstcn2 发表于 2018-7-25 23:50
这就可能涉及到MK检验图的解读问题了，趋势变化是看UF线，UF线的转折点才是我们通常理解的趋势变化点(f'( ...

M-K检验等属均值突变检验，建议采用滑动F检验或者改进的有序聚类分析等均方差突变检验方法，综合诊断突变点

青风海魅

sfhstcn2 发表于 2018-7-25 23:50
这就可能涉及到MK检验图的解读问题了，趋势变化是看UF线，UF线的转折点才是我们通常理解的趋势变化点(f'( ...

M-K检验等属均值突变检验，建议采用滑动F检验或者改进的有序聚类分析等均方差突变检验方法，综合诊断突变点

青风海魅

sfhstcn2 发表于 2018-7-25 23:50
这就可能涉及到MK检验图的解读问题了，趋势变化是看UF线，UF线的转折点才是我们通常理解的趋势变化点(f'( ...

M-K检验等属均值突变检验，建议采用滑动F检验或者改进的有序聚类分析等均方差突变检验方法，综合诊断突变点

青风海魅

zzcsunmoonstar 发表于 2018-7-25 21:04
构造的时间序列是一个很明显的先增后减突变，不知道为什么检测不出来？

M-K检验等属均值突变检验，建议采用滑动F检验或者改进的有序聚类分析等均方差突变检验方法，综合诊断突变点

sfhstcn2

青风海魅 发表于 2018-8-10 21:55
M-K检验等属均值突变检验，建议采用滑动F检验或者改进的有序聚类分析等均方差突变检验方法，综合诊断突变 ...

好建议，顶你

liuyg123

您好，您这个程序包内关于mk检验，函数内有一段s(k)=s(k-1)+r;请问为什么不是直接s(k)=r;呢？秩序列不是前面有多少大于他的数就加多少个1吗？为什么还要再加前一个秩序列的数呢

liuyg123

liuyg123 发表于 2018-8-29 11:41
您好，您这个程序包内关于mk检验，函数内有一段s(k)=s(k-1)+r;请问为什么不是直接s(k)=r;呢？秩序列不是前 ...

已经想明白了

sunjiajia

好不要用呢，为了积分找了好多mk检验