function?[Zbeta]=trendMK（xdt）
%Z?Mann-Kendall?統計量
%beta?傾斜度
%x?序列
%dt?數據的時間間隔
%
%%?序列長度檢驗
n=length（x）;
if?n<11
????message=‘數據點應多于10個‘;
????warning（message）
end
%%?計算統計量
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
%%?給出正態分布0.05處的值
normal=1.9

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

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

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

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