M-K趋势分析法

Mann-Kendall 非参数检验不需要数据遵从一定的分布，不受少数异常值的干扰，因此被广泛应用到水文气象数据的时间序列趋势分析中

  %我自己写的mk突变点分析，已经调试正确
%s是数据序列；sy是起始年份或序号（做图的时候会用到）
function [ufkubk]= mkabrpt（xsy）
sy=load（‘F:\MK\2.txt‘）%这个是年份
x =load（‘F:\MK\1.txt‘）%这个是径流数据
    N=length（x）;  %%F:\毕业论文参考文献\数学方法\mk\闽江水文站M_K\MkUB\月max\竹岐.txt
    ufk=smk（x）;
    y=flipud（x）;
    ubk=smk（y）;
    ubk=-ubk;
    ubk=flipud（ubk）;   
    %绘图
    ey=sy+N-1;
    htx=[sy:ey];
    
    a=ceil（abs（max（[ufk;ubk]））+3）;
    b=floor（min（[ufk;ubk]）-1）;
         
    if a<5
       a=5
    end
    
    if b>-4
       b=-4
    end
    
    figure（2）
    plot（htxufk‘r-‘‘linewidth‘1.5）;
    hold on
    plot（htxubk‘b-‘‘linewidth‘1.5）;
    plot（htx1.96*ones（N1）‘:‘‘linewidth‘1）;
    legend（‘UF统计量‘‘UB统计量‘‘0.05显著水平‘）;
    xlabel（‘年度‘‘FontName‘‘TimesNewRoman‘‘FontSize‘16）;
    ylabel（‘统计量‘‘Font

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1106  2014-08-06 10:18  mk程序\mkabrpt.m

     文件       1994  2015-08-12 22:43  mk程序\mk计算Z值.m

     文件        736  2014-08-06 10:18  mk程序\smk.m

     目录          0  2015-08-12 22:43  mk程序

----------- ---------  ---------- -----  ----

                 3836                    4