資源簡(jiǎn)介
M-K檢驗(yàn)法Matlab函數(shù)式M文件代碼,含計(jì)算、繪圖部分,能很好符合魏鳳英老師的例子
代碼片段和文件信息
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);
%…………………………………………………………………初始化……………………………………………………………………………
%計(jì)算統(tǒng)計(jì)量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
%計(jì)算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;
%寫(xiě)入文件
????????%DLMWRITE(‘MK1.txt‘UF);
%…………………………………………………………第一部分完成…………………………………………………………………………
%逆序
x=flipud(x);
%計(jì)算統(tǒng)計(jì)量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
%計(jì)算UB
for?i=2:l
????
評(píng)論
共有 條評(píng)論