mk突变点检测_MK检验突变分析 matlab

% Mann-Kendall突变检测

% 数据序列y

% 结果序列UFk，UBk2

%读取excel中的数据，赋给矩阵y

%获取y的样本数

%A为时间和降水数据列

x=降水(:,1);%时间序列

y=降水(:,2);%降水数据列

N=length(x);

n=length(y);

% 正序列计算---------------------------------

% 定义累计量序列Sk，长度=y，初始值=0

Sk=zeros(size(y));

% 定义统计量UFk，长度=y，初始值=0

UFk=zeros(size(y));

% 定义Sk序列元素s

s = 0;

% i从2开始，因为根据统计量UFk公式，i=1时，Sk(1)、E(1)、Var(1)均为0

% 此时UFk无意义，因此公式中，令UFk(1)=0

for i=2:n

for j=1:i

if y(i)>y(j)

s=s+1;

else

s=s+0;

end;

end;

Sk(i)=s;

E=i*(i-1)/4; % Sk(i)的均值

Var=i*(i-1)*(2*i+5)/72; % Sk(i)的方差

UFk(i)=(Sk(i)-E)/sqrt(Var);