利用混合乘同余法产生(0,1)均匀分布随机数.docVIP

概念 混合同余法： 用一常数与已产生的伪随机数之积的同余式做递推公式,产生均匀伪随机数的方法。 即： 第一步：递推式 其中： M=2^k,k2整数；A=1(mod4),即 为正整数,初值为非负整数。 第二步： 可以证明， 是周期为2k的伪随机数。 混合同余法（0，1）均匀分布的随机数 L=60,取k=2,c=1,n=5,l=60,=1 公式原理： 见上第一第二步。 程序： clear clc k=8; c=1; n=5; x=1; M=2^k; A=2^n+1; store=[ ]; for i=1:60 ee=x/M; store(i)=ee; x=mod(A*x+c,M); end plot(1:60,store) title(混合同余法产生随机数) store= 0.00390625000000000 0.132812500000000 0.386718750000000 0.765625000000000 0.269531250000000 0.898437500000000 0.652343750000000 0.531250000000000 0.535156250000000 0.664062500000000 0.917968750000000 0.296875000000000 0.800781250000000 0.4