trendMK 秩和检验MATLAB实现

利用MATLAB 实现Mann Kendall非参数检验，秩和检验

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