randu1.m

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

function randmtx=randU1(row,col,seed,k1,k2)
if nargin<3
    seed=1973272912;
end
m=2^31-1;
a=630360016;
randmtx(1,1)=mod(a*seed,m);
for i=2:row*col
    randmtx(1,i)=mod((a*randmtx(1,i-1)),m);
end
randmtx=reshape(randmtx,row,col);
if nargin<4
    randmtx=randmtx/m;
    return
end
randmtx=randmtx/m;
randmtx=randmtx*(k2-k1)+k1;