|
****昨天有位朋友想了解随机数产生
设X是一个在(0,1)间均匀分布的随机变量U(0,1),其采样值为x,则有Y=F-1(X)为满足概率分布F(y)的随机变量.因而产生变量Y的随机数y可以分两步:
(1) 产生U(0,1)的随机数x;
(2) y= F-1(X).
如何做到(1)呢?
一种简单的近似做法为x(I+1)=(ax(i)+c) mod m. c不为0时这种方法称为混合同余法,c为0时称为乘同余法. a和m一般取质数.
还可以通过外部的触发条件来实现.如原了在某一时刻的的空间位置是完全随机的,如单位时间内原子衰变的个数也是服从一定的分布.
|
|