Matlab產生隨機數

Matlab(mathworks.com) 隨機數生成方法:

第一種方法是用 random 語句,其一般形式爲
                     y = random('分佈的英文名',A1,A2,A3,m,n),
表示生成 m 行 n 列的 m × n 個參數爲 ( A1 , A2 , A3 ) 的該分佈的隨機數。例如:
(1) R = random('Normal',0,1,2,4): 生成期望爲 0,標準差爲 1 的(2 行 4 列)2× 4 個正態隨機數
(2) R = random('Poisson',1:6,1,6):  依次生成參數爲 1 到 6 的(1 行 6 列)6 個 Poisson 隨機數


第二種方法是針對特殊的分佈的語句:
一. 幾何分佈隨機數  (下面的 P,m 都可以是矩陣)
   R = geornd(P)   (生成參數爲 P 的幾何隨機數)
   R = geornd(P,m)  (生成參數爲 P 的 × m 個幾何隨機數)
                                     1
   R = geornd(P,m,n)  (生成參數爲 P 的 m 行 n 列的 m × n 個幾何隨機數)
    例如
(1)  R = geornd(1./ 2.^(1:6)) ( 生成參數依次爲 1/2,1/2^2,到 1/2^6 的 6 個幾何隨機數)
(2)  R = geornd(0.01,[1 5]) (生成參數爲 0.01 的(1行5列)5 個幾何隨機數).


二.Beta 分佈隨機數
R = betarnd(A,B)  (生成參數爲 A,B 的 Beta 隨機數)
R = betarnd(A,B,m)  (生成 × m 個數爲 A,B 的 Beta 隨機數)
                           1
R = betarnd(A,B,m,n)  (生成 m 行 n 列的 m × n 個數爲 A,B 的 Beta 隨機數).


三.正態隨機數
R = normrnd(MU,SIGMA)  (生成均值爲 MU,標準差爲 SIGMA 的正態隨機數)
R = normrnd(MU,SIGMA,m)  (生成 1× m 個正態隨機數)
                                      
R = normrnd(MU,SIGMA,m,n) (生成 m 行 n 列的 m × n 個正態隨機數)
   例如
(1) R = normrnd(0,1,[1 5])   生成 5 個正態(0,1) 隨機數
                                 
(2) R = normrnd([1 2 3;4 5 6],0.1,2,3)  生成期望依次爲[1,2,3;4,5,6], 方差爲 0.1 的 2× 3 個正態隨機數.


四.二項隨機數:類似地有
R = binornd(N,P)  R = binornd(N,P,m)   R = binornd(N,p,m,n)
   例如
   n = 10:10:60;   r1 = binornd(n,1./n)  或 r2 = binornd(n,1./n,[1 6]) (都生成參數分別爲
   1          1   ), L, ( 60, ) 的6個二項隨機數.
(10,
    10          60


五.自由度爲 V 的 χ 2 隨機數:
R = chi2rnd(V)    R = chi2rnd(V    R = chi2rnd(V
                                     ,m)             ,m,n)


六.期望爲 MU 的指數隨機數(即 Exp                      隨機數):
                                       1
                                       MU
R = exprnd(MU)   R = exprnd(MU,m)  R = exprnd(MU,m,n)


七.自由度爲 V1, V2 的 F 分佈隨機數:
   R = frnd(V1,V2)   R = frnd(V1, V2,m)  R = frnd(V1,V2,m,n)


八. Γ ( A, λ ) 隨機數:
   R = gamrnd(A,lambda)  R = gamrnd(A,lambda,m)  R = gamrnd(A,lambda,m,n)


九.超幾何分佈隨機數:
   R = hygernd(N,K,M)   R = hygernd(N,K,M,m)  R = hygernd(N,K,M,m,n)

表示N個產品中有K個次品,抽取M次湊到R個次品


十.對數正態分佈隨機數
   R = lognrnd(MU,SIGMA)  R = lognrnd(MU,SIGMA,m)  R = lognrnd(MU,SIGMA,m,n)


十一.負二項隨機數:
   R = nbinrnd(r,p)   R = nbinrnd(r,p,m)   R = nbinrnd(r,p,m,n)


十二.Poisson 隨機數:
   R = poissrnd(lambda)   R = poissrnd(lambda,m)  R = poissrnd(lambda,m,n)
    例如,以下 3 種表達有相同的含義:lambda = 2;  R = poissrnd(lambda,1,10)
(或 R = poissrnd(lambda,[1 10])  或 R = poissrnd(lambda(ones(1,10)))


十三.Rayleigh 隨機數:
   R = raylrnd(B)    R = raylrnd(B,m)   R = raylrnd(B,m,n)


十四.V 個自由度的 t 分佈的隨機數:
   R = trnd(V)    R = trnd(V,m)   R = trnd(V,m,n)


                                              42
十五.離散的均勻隨機數:
R = unidrnd(N)   R = unidrnd(N,m)  R = unidrnd(N,m,n)


十六.[A,B] 上均勻隨機數
R = unifrnd(A,B)   R = unifrnd(A,B,m)  R = unifrnd(A,B,m,n)
例如 unifrnd(0,1:6)與 unifrnd(0,1:6,[1 6]) 都依次生成[0,1] 到[0,6]的6個均勻隨機數.:


十七.Weibull 隨機數
R = weibrnd(A,B)   R = weibrnd(A,B,m)  R = weibrnd(A,B,m,n)

十八.randperm(n)產生1到n的整數的無重複的隨機排列


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章