%?Mann-Kendall突變檢測?
%?數(shù)據(jù)序列y
%?結(jié)果序列UFk，UBk2
%--------------------------------------------
%讀取excel中的數(shù)據(jù)，賦給矩陣y
%獲取y的樣本數(shù)
%A為時間和徑流數(shù)據(jù)列
A=xlsread（‘宜昌.xlsx‘）;
x=A（:1）;%時間序列
y=A（:2）;%徑流數(shù)據(jù)列
N=length（y）;
n=length（y）;
%?正序列計算---------------------------------
%?定義累計量序列Sk，長度=y，初始值=0
Sk=zeros（size（y））;
%?定義統(tǒng)計量UFk，長度=y，初始值=0
UFk=zeros（size（y））;
%?定義Sk序列元素s
s?=?0;
%?i從2開始，因為根據(jù)統(tǒng)計量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
%?逆序列計算---------------------------------
%?構(gòu)造逆序列y2，長度=y，初始值=0
y2=z

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

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

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

-----------?---------??----------?-----??----

???????????????111980????????????????????2