function?[Zbeta]=trendMK（xdt）
%Z?Mann-Kendall?統(tǒng)計(jì)量
%beta?傾斜度
%x?序列
%dt?數(shù)據(jù)的時(shí)間間隔
%
%%?序列長(zhǎng)度檢驗(yàn)
n=length（x）;
if?n<11
????message=‘?dāng)?shù)據(jù)點(diǎn)應(yīng)多于10個(gè)‘;
????warning（message）
end
%%?計(jì)算統(tǒng)計(jì)量
S=0;
for?i=1:n-1
????for?j=i+1:n
????????S=S+sign（x（j）-x（i））;
????end
end
var_S=（n*（n-1）*（2*n+5））/18;
if?S==0
????Z=0;
elseif?S>0
????Z=（S-1）/sqrt（var_S）;
else
????Z=（S+1）/sqrt（var_S）;
end
%%?給出正態(tài)分布0.05處的值
normal=1.9

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????914??2010-09-01?10:00??trendMK.m

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

??????????????????914????????????????????1