%產(chǎn)生一組[0，1]分布的隨機(jī)數(shù)，采用余數(shù)法
%從文獻(xiàn)中得知下列參數(shù)組較為適用：
%y的初值為1，w=7,模M=10^10(偽隨機(jī)數(shù)周期為5*10^7)
%y的初值為1，w=5^13,模M=2^36(偽隨機(jī)數(shù)周期為2^34,約2*10^10)
%y的初值為1，w=5^17,模M=2^42(偽隨機(jī)數(shù)周期為2^40,約10^12)
%--------------------------------------------------------
%function x=uniform(y,w,M,n)
%n為要產(chǎn)生的隨機(jī)數(shù)個(gè)數(shù)
function x=uniform(y,w,M,n)
for i=1:n
    x(i)=y./M;
    y=rem(w*y,M);%rem()為求余函數(shù)
end