%最近寫論文需要用到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;