mk突变分析 matlab程序

A=xlsread（‘d:/2.xlsx‘）;
x=A（:1）;
y=A（:2）;
N=length（y）;
n=length（y）;
Sk=zeros（size（y））;
UFk=zeros（size（y））;
s=0;
for i=2:n
    for j=1:i
        if y（i）>y（j）
            s=s+1;
        else
            s=s+0;
        end;
    end
    Sk（i）=s;
    E=i*（i-1）/4;
    Var=i*（i-1）*（2*i+5）/72;
    UFk（i）=（Sk（i）-E）/sqrt（Var）;
end;
y2=zeros（size（y））;
Sk2=zeros（size（y））;
UBk=zeros（size（y））;
s=0;
for i=1:n
    y2（i）=y（n-i+1）;
end;
for i=2:n
    for j=1:i
        if y2（i）>y2（j）
            s=s+1;
        else
            s=s+0;
        end;
    end;
    Sk2（i）=s;
    E=i*（i-1）/4;
    Var=i*（i-1）*（2*i+5）/72;
    U