随机事件的模拟
一、实验目的
会使用随机数命令rnd,进而求服从各种分布的随机数列(样本)。熟悉随机事件的模拟。
二、预备知识
随机数产生器
所有分布的随机数的产生方法都始于均匀分布随机数。一旦具备了均匀分布随机数产生器,其他分布
的随机数都可使用直接法、反转(inversion)法或拒绝(rejection)法获得。
(1) 随机数产生的基本方法
1 直接法(direct)
直接法源于分布的定义。假设产生二项分布随机数,可以认为其随机数就是在n次抛硬币之后,某一
面出现的次数(每次抛掷时,此面出现的概率为p)。如果产生了n个均匀分布随机数,数出大于p的次数m,
那么结果m就是参数为n和p的二项分布的随机数。
2 反演法(inversion)
反演方法的理论基础在于均匀分布与其他连续分布之间的关系。假设P为一连续分布,其逆为F ;U是
一个均匀分布随机数,则 (U)服从F分布。因此,可将某种分布的逆函数作用于均匀分布的随机数,获得
这种分布的随机数。遗憾的是,这种方法通常并不是最有效的。
3 拒绝法(rejection)
对于某些分布,其函数形式使得运用直接法和反演法来产生随机数比较困难或费时较多。在这种情况
下,拒绝法或许能很好地解决这一问题。
若要产生概率密度函数为f的某种分布的随机数,若采用拒绝法,首先需找到另一分布密度函数g和一
个常数c,并满足以下条件:
f(x)<=cg(x),
然后通过以下步骤进行:
1. 产生概率密度为g(x)的G分布的随机数x
2. 设参数 ;
3. 产生一个均匀分布随机数u
4. 如果r× <1,返回x
5.否则重复第1步至第3步。
为提高效率,产生G分布随机数的方法要简单,而且常数c的值要小。迭代的期望值为c。
(2) 产生随机数的通用函数
在工具箱中,提供了通用的随机数产生因数random和特定分布的随机数产生函数(以rnd结尾)。可以
直接调用这些函数来获得所需的随机数,而不必经过(1)中所述的过程。
random( )
功能:产生可选分布的随机数。
格式: y=random(‘name’,A1,A2,A3,m,n)
说明: random函数产生统计工具箱中任一分布的随机数。‘name’为相应分布