在數位電腦上用數學方法產生的、統計意義下具有在區間(0,1)上均勻總體簡單子樣性質的數值序列{rn|n=1,2,…;0≤rn≤1}。
用蒙特卡羅法模擬求解一個實際問題,要用到各種不同分佈的隨機變數、隨機向量量和隨機過程η的抽樣序列{ηn|n=1,2,…},稱它們為隨機數。如常用的二項分佈隨機數、均勻分佈隨機數、二維正態分佈隨機數等,其中最基本、最重要的是區間(0,1)上均勻分佈的隨機數。因此,如何在計算機上產生偽隨機數備受重視。
在一臺b)進制(如二進制或十進制)、尾數字長為k位的計算機上,不考慮符號和階碼,可以表示bk個不同的數,即0,1,2,…,bk-1。在數學計算機上產生偽隨機數,就是選取m個整數x1,x2,…,xm作為初值和一個適於遞推計算的數學公式
![](/img3/10666.gif)
![](/img3/10667.gif)
![](/img3/10668.gif)
選用不同的遞推計算公式
![](/img3/10669.gif)
![](/img3/10670.gif)
![](/img3/10671.gif)
![](/img3/10672.gif)
![](/img3/10673.gif)
在一臺尾數字長為k位的二進制計算機上,取模M=2k,乘同餘法的遞推計算同餘式為
![](/img3/10674.gif)
![](/img3/10675.gif)
![](/img3/10675.gif)
![](/img3/10676.gif)
對產生的偽隨機數,要經過一定的理論分析和各種統計檢驗,以檢查得到的序列是否具有在區間(0,1)上均勻總體簡單子樣所應具有的各種統計性質,如分佈的均勻性、取值的隨機性、前後的獨立性和分段序列統計性質的一致性等。進行上述統計性質檢驗的方法很多,常用的有參數檢驗、均勻性檢驗、獨立性檢驗、連檢驗和各種不同的組合規律性檢驗等。
有瞭偽隨機數{rn},利用各種不同的抽樣算法,如直接抽樣、變換抽樣、舍選抽樣、復合抽樣等,就可以產生模擬計算中需要的各種不同分佈的隨機數。
設隨機變量η的分佈函數F(x)連續,逆函數F-1(y)存在,則R=F(η)為區間(0,1)上均勻分佈的隨機變量。利用這一原理,從隨機數{rn}出發,就可以直接得到η的抽樣序列
![](/img3/10677.gif)
![](/img3/10678.gif)
![](/img3/10679.gif)
![](/img3/10680.gif)
在概率統計的理論研究和實際應用中,經常遇到具有密度函數
![](/img3/10681.gif)
![](/img3/10682.gif)
上述抽樣算法,要用到對數、開方、正弦和餘弦等算法,速度較慢。在計算機上,靈活多變、計算量省的舍選抽樣和復合抽樣更經常的用來產生所要的各種不同分佈的隨機數。以隨機餘弦
![](/img3/10683.gif)
![](/img3/10684.gif)
![](/img3/10685.gif)
![](/img3/10686.gif)
![](/img3/10687.gif)
參考書目
中國科學院計算中心概率統計組編著:《概率統計計算》,科學出版社,北京,1979。
徐鐘濟編著:《蒙特卡羅方法》,上海科學技術出版社,上海,1985。