混合同余法產(chǎn)生均勻分布隨機數(shù)產(chǎn)生方法總結

1、混合同余法產(chǎn)生均勻分布隨機數(shù)產(chǎn)生方法 總結主要學習混合同余法產(chǎn)生各種分布的隨機數(shù)的方法，見參考文獻1, 2，重點參考2。其中要注意混合同余法產(chǎn)生隨機數(shù)的參數(shù)的選取。1混合同余法產(chǎn)生均勻分布的隨機數(shù)1.1混合同余法通過同余運算生成偽隨機數(shù)的方法稱為同余法，常用的同余法包括加同余法、 乘同余法、混合同余法、除同余法。其中乘同余法和混合同余法的性能更好，有 速度快、內存省、周期長、統(tǒng)計特性好等優(yōu)點?；旌贤喾ㄊ荓ehmer在1951年提出的，其迭代公式為 ：Xn =mod(AXn-1 +C, M)(1.1)(1.2)Yn=Xn/M公式(1.1)、(1.2)中，mod表示求余函數(shù)，A,C, M均為正整

2、數(shù)。其中M是模 數(shù)，A是乘數(shù)，C是增量，Xo為初始值(0空X。乞M)，當C=0時，稱此算法為 乘同余法；若C = 0，則稱算法為混合同余法，當C取不為零的適當數(shù)值時，有 一些優(yōu)點，但優(yōu)點并不突出，故常取 C=0。Xn是在(0, M )內服從均勻分布的隨 機變量，Y,則是在(0, 1)內服從均勻分布的隨機變量。式中X。, A,C, M的取值并 不是隨意的，模M大小是發(fā)生器周期長短的主要標志，常見有 M為素數(shù)，取A 為M的原根，則周期T=M -1。試驗統(tǒng)計表明，用以下參數(shù)進行混合同余法產(chǎn)生 的隨機序列的統(tǒng)計特性較好：Xn=mod(314159269X n-1+453806245,231)(1.3)

3、(1.4)Xn=mod(515Xn-1+1,235)735Xn=mod(2 Xn-1+1, 2 )(1.5)Xn二mod(2045X n-1+1, 220)(1.6)107Xn=mod(7X n-1,10 ),Xo=1,T=5 10(1.7)Xn=mod(5 13Xn-1,236),X 0 = 1,T=23 2 1010(1.8)Xn=mod(5 17Xn-1,212),Xo = 1,T=2 401012(1.9)31Xn=mod(32719X n-1, 2 -1), X。任意(1.10)31Xn=mod(16807X n-1, 2 -1), X0任意(1.11)31Xn=mod(122070

4、3125X n-1, 2 -1),X0任意(1.12)15Xn=mod(9869X n-1 +6925, 2 -1)(1.13)3131在式(1.10)(1.12)中 T =2 -2 , 16807、32719、1220703125 都是 2 -1 的原根?；旌贤喾óa(chǎn)生的隨機序列具有以下特點：Xn重復周期較小，由于Xn取值在(0, M )內，其周期TEM ,T受X°,A,C,M的值得影響，在編程實現(xiàn)時，浮點運算也會對T產(chǎn)生影響用此方法產(chǎn)生的隨機序列，在一個周期內任意兩個隨機數(shù)不可能相等，這往往與實際情況不相符經(jīng)Hull和Dobell證明，只有X0,A,C,M滿足以下一些關系才能實現(xiàn)

5、周期最 大化，即T=M，條件如下：C與M互質(或互素，即它們的最大公約數(shù)為1)設q為某一質數(shù)，M分別能被q和4整除，且(A-1 )能被q和4整除產(chǎn)生具有最大周期的偽隨機序列的混合同余法算法為：kXn =mod(4a+1)Xn-1(2c+1), 2 )(1.14)Yn=Xn/2k(1.15)k由于M=2 ,k-2時，M只有一個素數(shù)因子 2，且4也是M的因子，此時 A=4a+1，正好滿足了 T=M的第二個條件；而此時C=2c+1剛好與M互質，即滿足T = M的第一個條件1.2改進的混合同余法改進的混合同余法的迭代公式如下：Xn =mod(AXn-i n, M)(1.16)Yn=Xn/M(1.17)

6、改進的混合同余法具有以下特點：比混合同余法產(chǎn)生的周期長，T _M允許某個偽隨機數(shù)重復發(fā)生，且重復發(fā)生的次數(shù)為T / M 偽隨機序列的周期T一般與初始值Xo的選取無關，只有極個別的情況除外1.3原根相關知識1.3.1歐拉函數(shù)在數(shù)論，對正整數(shù)n,歐拉函數(shù)是少于或等于n的數(shù)中與n互質的數(shù)的數(shù)目。 此函數(shù)以其首名研究者歐拉命名，它又稱為Euler's totient function、©函數(shù)、歐拉商數(shù)等。例如© (8)=4因為1,3,5,7均和8互質。1.3.2原根定義1 設m > 1，(a, m) = 1，則使ra 三 1(mod m)(1.18)成立的最小的正整數(shù)

7、r，稱為a對模m的指數(shù)，記為m(a)，在不致誤會的情況下， 簡記為(a)0由Euler定理，當r = (m)時式(1)成立，因此，恒有 m(a) (m)。若 a 三 b (mod m)，(a, m) = 1，則顯然有 亦(a) = 6m(b)。定義2 若、m(a) = (m)，則稱a是模m的原根。例如，當m = 7時，因為12321 三 2，22三 4，23三 1 (mod 7)，所以、7(2) = 3。又因為12345631 三 3, 32三 2, 33三 6, 34三 4, 35三 5, 36三 1 (mod 7),所以刀(3) = 6 = (7), 3是模7的原根。以后，在談到a對模m的指數(shù)時，總假定m > 1, (a, m)=參考文