mk.m

用于mann-kendall趋势分析

function data=mk(v)

n=length(v);
for i=1:n-1
    k(i)=0;
    for j=i+1:n
        temp=v(j)-v(i);
        if temp>0
            k(i)=k(i)+1;
        elseif temp<0
            k(i)=k(i)-1;
        end
    end
end
data=sum(k);