《信息隐藏技术实验教程》MATLAB源码

%独立性检验的统计函数 function result=freq(rando,rande) [rowo,colo]=size(rando); [rowe,cole]=size(rande); result=zeros(11,11); for i=1:colo k=1; if rando(1,i)<=0.1 & rando(1,i)>0 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,i)>0.6 & rande(1,i)<=0.7 result(k,7)=result(k,7)+1; end if rande(1,i)>0.7 & rande(1,i)<=0.8 result(k,8)=result(k,8)+1; end if rande(1,i)>0.8 & rande(1,i)<=0.9 result(k,9)=result(k,9)+1; end if rande(1,i)>0.9 & rande(1,i)<1 result(k,10)=result(k,10)+1; end end k=2; if rando(1,i)<=0.2 & rando(1,i)>0.1 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,i)>0.6 & rande(1,i)<=0.7 result(k,7)=result(k,7)+1; end if rande(1,i)>0.7 & rande(1,i)<=0.8 result(k,8)=result(k,8)+1; end if rande(1,i)>0.8 & rande(1,i)<=0.9 result(k,9)=result(k,9)+1; end if rande(1,i)>0.9 & rande(1,i)<1 result(k,10)=result(k,10)+1; end end k=3; if rando(1,i)<=0.3 & rando(1,i)>0.2 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,i)>0.6 & rande(1,i)<=0.7 result(k,7)=result(k,7)+1; end if rande(1,i)>0.7 & rande(1,i)<=0.8 result(k,8)=result(k,8)+1; end if rande(1,i)>0.8 & rande(1,i)<=0.9 result(k,9)=result(k,9)+1; end if rande(1,i)>0.9 & rande(1,i)<1 result(k,10)=result(k,10)+1; end end k=4; if rando(1,i)<=0.4 & rando(1,i)>0.3 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,i)>0.6 & rande(1,i)<=0.7 result(k,7)=result(k,7)+1; end if rande(1,i)>0.7 & rande(1,i)<=0.8 result(k,8)=result(k,8)+1; end if rande(1,i)>0.8 & rande(1,i)<=0.9 result(k,9)=result(k,9)+1; end if rande(1,i)>0.9 & rande(1,i)<1 result(k,10)=result(k,10)+1; end end k=5; if rando(1,i)<=0.5 & rando(1,i)>0.4 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,i)>0.6 & rande(1,i)<=0.7 result(k,7)=result(k,7)+1; end if rande(1,i)>0.7 & rande(1,i)<=0.8 result(k,8)=result(k,8)+1; end if rande(1,i)>0.8 & rande(1,i)<=0.9 result(k,9)=result(k,9)+1; end if rande(1,i)>0.9 & rande(1,i)<1 result(k,10)=result(k,10)+1; end end k=6; if rando(1,i)<=0.6 & rando(1,i)>0.5 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,i)>0.6 & rande(1,i)<=0.7 result(k,7)=result(k,7)+1; end if rande(1,i)>0.7 & rande(1,i)<=0.8 result(k,8)=result(k,8)+1; end if rande(1,i)>0.8 & rande(1,i)<=0.9 result(k,9)=result(k,9)+1; end if rande(1,i)>0.9 & rande(1,i)<1 result(k,10)=result(k,10)+1; end end k=7; if rando(1,i)<=0.7 & rando(1,i)>0.6 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,i)>0.6 & rande(1,i)<=0.7 result(k,7)=result(k,7)+1; end if rande(1,i)>0.7 & rande(1,i)<=0.8 result(k,8)=result(k,8)+1; end if rande(1,i)>0.8 & rande(1,i)<=0.9 result(k,9)=result(k,9)+1; end if rande(1,i)>0.9 & rande(1,i)<1 result(k,10)=result(k,10)+1; end end k=8; if rando(1,i)<=0.8 & rando(1,i)>0.7 if rande(1,i)>0 & rande(1,i)<=0.1 result(k,1)=result(k,1)+1; end if rande(1,i)>0.1 & rande(1,i)<=0.2 result(k,2)=result(k,2)+1; end if rande(1,i)>0.2 & rande(1,i)<=0.3 result(k,3)=result(k,3)+1; end if rande(1,i)>0.3 & rande(1,i)<=0.4 result(k,4)=result(k,4)+1; end if rande(1,i)>0.4 & rande(1,i)<=0.5 result(k,5)=result(k,5)+1; end if rande(1,i)>0.5 & rande(1,i)<=0.6 result(k,6)=result(k,6)+1; end if rande(1,