mulmod1.m

来自「关于随机数生成的的源代码」· M 代码 · 共 13 行

function r = MulMOD1(L,alpha,beta,n)
r = zeros(n,1);
x = zeros(n,1);
M = power(2,L);
a = 8*alpha + 5;
x(1) = 4*beta + 1;
r(1) = x(1)/M;

for i=2:n
    y = a*x(i-1);
    x(i) = mod(y,M);
    r(i) = x(i)/M;
end