Mann-Kendall突变检测

% Mann-Kendall突变检测 % 数据序列y % 结果序列UFk，UBk2 %-------------------------------------------- %读取excel中的数据，赋给矩阵y %获取y的样本数 %A为时间和径流数据列

% Mann-Kendall突变检测 
% 数据序列y
% 结果序列UFk，UBk2
%--------------------------------------------
%读取excel中的数据，赋给矩阵y
%获取y的样本数
%A为时间和径流数据列
A=xlsread（‘宜昌.xlsx‘）;
x=A（:1）;%时间序列
y=A（:2）;%径流数据列
N=length（y）;
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）;
end;
% ------------------------------正序列计算end
% 逆序列计算---------------------------------
% 构造逆序列y2，长度=y，初始值=0
y2=z

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     109058  2019-04-04 15:41  宜昌.xlsx

     文件       2922  2017-03-13 17:10  mktest.m

----------- ---------  ---------- -----  ----

               111980                    2