MK突变检验

clear all
data=load（‘E:\PhD paper\Huai River_Actual evapotranspiration\MK_mutation Test\ETa_AA model.txt‘）;

% Mann-Kandel 突变点分析MATLAB脚本文件，使用前预先导入目标数据，变量名设定为data，再将内容复制至命令框即可。
N=size（data1）;
n=2;
nb=1;
nB=2;
ri=zeros（N1）;
rbi=zeros（N1）;
SFk=zeros（N1）;
SBk=zeros（N1）;
UFk=zeros（N1）;
UBk=zeros（N1）;
while n<=N
    i=1;
    r=zeros（1N）;
    while i        if data（i）<=data（n）
            r（i）=1;
        else
            r（i）=0;
        end
         i=i+1;
    end
    ri（n）=sum（r）;
    SFk（n）=SFk（n-1）+ri（n）;
    UFk（n）=（SFk（n）-n*（n-1）/4）/sqrt（n*（n-1）*（2*n+5）/72）;
    n=n+1;
end
while nb<=N-1
    i=N;
    rb=zeros（1N）;
    while i>nb
        if data（i）<=data（nb）
            rb（i）=1;
        else
            rb（i）=0;
        e