producerand.m

来自「共上载5个源码」· M 代码 · 共 14 行

function prorand=producerand(hang,lie,seed,m,a,c)
if nargin<6  c=453806245;      end
if nargin<5  a=314159269;      end
if nargin<4  m=2^31;         end
if nargin<3  seed=2^10 ;      end
if nargin<2  lie=1;             end
if nargin<1  hang=1;           end

x(1)=mod((a*seed+c),m);
for i=2:hang*lie
     x(i)=mod((a*x(i-1)+c),m);
end
x=x/m;
prorand=reshape(x,hang,lie);