Manner-Kendall(M-K)---突变检验

原理

去看原文

代码

data=xlsread('C:\Users\Desktop\mk突变检验\data.xls','Sheet1','A2:B52');
x=data(:,1); % 时间列
y=data(:,2); % 数据列
N=length(y);
n=length(y);
Sk=zeros(size(y));
UFk=zeros(size(y));
s=0;
for i=2:n
   for j=1:i
         if y(i)>y(j)
           s=s+1;
         else
           s=s+0;
         end
   end
   ...
end
%...练习数据及完整代码后台回复“MK突变”获取。

结果

非平稳时间序列突变检测 – Bernaola Galvan分割算法

非平稳时间序列突变检测 – Bernaola Galvan分割算法

Pettitt突变检验

Pettitt突变检验

原文有源码，更多内容，请关注地学数据处理分析。