M-K检验法MatlabM文件代码

function pp=MK（x）
l=length（x）;
r=zeros（l1）;
s=zeros（l1）;
E=zeros（l1）;
V=zeros（l1）;
UF=zeros（l1）;
UB=zeros（l1）;
%…………………………………………………………………初始化……………………………………………………………………………

%计算统计量s
for i=1:l
    for j=1:i
        if x（i）>x（j）
           r（i）=r（i）+1;
        end
    end
s（i）=sum（r）;
end
%计算UF
for i=2:l
    E（i）=i*（i-1）/4;
    V（i）=i*（i-1）*（2*i+5）/72;
    UF（i）=（s（i）-E（i））/sqrt（V（i））;
end
UF（1）=0;
%写入文件
        %DLMWRITE（‘MK1.txt‘UF）;
%…………………………………………………………第一部分完成…………………………………………………………………………

%逆序
x=flipud（x）;
%计算统计量s
r=zeros（l1）;
for i=1:l
    for j=1:i
        if x（i）>x（j）
           r（i）=r（i）+1;
        end
    end
s（i）=sum（r）;
end
%计算UB
for i=2:l