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突变检验
原文有源码,更多内容,请关注地学数据处理分析。