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的整數的無重複的隨機排列