Mann-mandall趋势检验 （Pettitt and Mk）

Mann-mandall趋势检验，Pettitt突变点检验

function [ mk ] =MK（ x ）
%UNtitleD Summary of this function goes here
%   Detailed explanation goes here
X=x;
n=length（X）;
S=0;
    for k=1:（n-1）
        for j=（k+1）:n
            S=S+sign（X（j）-X（k））;
        end   
    end
 VarS=n*（n-1）*（2*n+5）/18;
    if S>0
        Z=（S+1）/sqrt（VarS）;
    elseif S==0
        Z=0;
    else
        Z=（S-1）/sqrt（VarS）;
    end
 
Qi=zeros（n*（n-1）/21）;
counter=1;
for k=1:（n-1）
    for j=（k+1）:n
        Qi（counter）=（X（j）-X（k））/（j-k）;
        counter=counter+1;
    end
end
Alpha=1-normcdf（Z01）;
disp（‘趋势的显著水平为‘）;
disp（Alpha）;

Q=median（Qi）;
disp（‘整体趋势的变化速率为‘）;
disp（Q）;
mk=[ZAlphaQ];
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         433  2014-05-03 12:45  Pettitt.asv
     文件         503  2014-05-04 17:40  Pettitt.m
     文件         503  2014-05-06 10:15  Pettitt2.asv
     文件         537  2014-05-06 11:04  Pettitt2.m
     文件         678  2014-05-06 11:43  MK.asv
     文件         687  2014-05-06 12:55  MK.m