runstest.m

来自「伪随机数发生器: “利用随机序列进行间隔控制」· M 代码 · 共 20 行

function [r,zr]=runstest(randomnum)
if size(randomnum(1,:))<=20
    error('输入序列太短');
end
logic=(randomnum>=0.5);
run=abs(diff(logic));
n=length(logic);
n1=sum(logic)+1;
n2=n-n1;
r=sum(run);
u_r=(2*n1*n2)/(n1+n2)+1;
std_r=sqrt((2*n1*n2*(2*n1*n2-n1-n2))/(((n1+n2)^2)*(n1+n2-1)));
zr=(r-u_r)/std_r;
if abs(zr)<sz
    disp('接受独立假设');
else
    disp('拒绝独立假设');
end