function s=rand01(p,m,n)
%產(chǎn)生一個(p,1-p)的0－1隨機變量
%輸入?yún)?shù)： 
%       p:  0－1分布中0的概率
%       m,n:產(chǎn)生的隨機變量樣本個數(shù)m×n   
%輸出： 產(chǎn)生的隨機變量樣本矢量
x = rand(m,n);
s = ( sign( x-p+eps )+1 )/2;