mk检验程序269770

%最近写论文需要用到MK检验法，网上收集到大量的matlab代码，但是没有一个代码能够
%完全正确运行或者分析信息不全，结合多位网友编写的MK检验法，经过我的改编，顺利得到
%正确的运行结果，谢谢各位网友，希望对有需要的盆友有帮助
% Mann-Kendall突变检测 
% 数据序列y
% 结果序列UFk，UBk2
%--------------------------------------------
%读取excel中的数据，赋给矩阵y
%获取y的样本数
%A为时间和径流数据列
A=xlsread（‘C:\Users\LIJIANCHENG\Desktop\1.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=zeros（size（y））;
% 定义逆序累计量序列Sk2，长度=y，初始值=0
Sk2=zeros（size（y））;
% 定义逆序统计量UBk，长度=y，初始值=0
UBk=zeros（size（y））;
% s归0
s=0;
% 按时间序列逆转样本y
% 也可以使用y2=flipud（y）;或者y2=flipdim（y1）;
for i=1:n
    y2（i）=y（n-i+1）;
end;
% i从2开始，因为根据统计量UBk公式，i=1时，Sk2（1）、E（1）、Var（1）均为0
% 此时UBk无意义，因此公式中，令UBk（1）=0
for i=2:n
   for j=1:i
         if y2（i）>y2（j）
           s=s+1;