一维正态分布随机数序列产生的方法

2021-09-06 13:27:57 浏览次数：

一维正态分布随机数序列产生的方法【摘要】正态分布在数理统计中具有基础性的作用，因此产生高质量的正态分布有重要的意义。现在我们总结了一些方法来介绍一维正态分布随机数序列的生法中心极限定理，Hasiting 有理逼近法，统计工具箱，反函数法，舍选法，R软件及一维正态随机数的检验。【关键词】正态分布；一维；随机数。正态分布的定义若连续随机变量的概率密度为 然而，得到服从正态分布的随机数的基本思想是先得到服从均匀分布的随机数，再将服从均匀分布的随机数转变为服从正态分布。一维正态分布随机数序列产生的方法一、用matlab统计工具箱产生R normrndMU，SIGMA （生成均值为 MU，标准差为 SIGMA 的正态随机数）R normrndMU，SIGMA,m（生成1 m个正态随机数）R normrndMU，SIGMA,m,n（生成m行n列的 m n 个正态随机数）例如 1 R normrnd0,1,1 5生成5个正态0,1 随机数 2 R normrnd1 2 3;4 5 6,0.1,2,3 生成期望依次为1,2,3;4,5,6,方差为 0.1 的 2 3 个正态随机数。二利用分布函数的反函数1、设连续型随机变量Y的概率函数为 fx,需产生给定分布的随机数.算法（1）产生n个RND 随机数r1，r2，，rn； （2）所得yi,i1,2,,n 即所求.2、基本原理设随机变量Y的分布函数Fy是连续函数，而且随机变量XU0,1，令ZF1X。则Z与Y有相同分布。

　　证明 FZz PF1X z PXFzGFz Fz 因Gx是随机变量X 的分布函数若Y的概率密度为 fy，由YF1X可得对给出定的0,1上均匀分布随机数ri，则具有给定分布的随机数 yi 可由方程 解出。

　　若要得到分布函数为Fx的随机变量Y。可令，其中u是服从均匀分布的随机变量，有 因而，对于任意的分布函数，只要求出它的反函数，就可以由服从均匀分布的随机变量实例来生成服从该分布函数的随机变量实例。现在来看正态分布的分布函数，对于，其分布函数为 显然，要想求其反函数是相当困难的，同时要想编程实现也很复杂。可见，用此种方法来生成服从正态分布的随机变量实例并不可取。三利用中心极限定理1、2、产生服从N,2的算法步骤（1）产生n 个RND 随机数r1，r2，，rn；（2）3 计算 yx ，y 是服从 N,2 分布的随机数。3、原理分析设1，2，，n是n个相互独立的随机变量,且iU0,1,i 1,2,,n,由中心极限定理知 ，渐近服从正态分布N0,l 。注意我们现在已经能产生0，1均匀分布的随机数了，那么我们可以利用这个定理来产生标准正态分布的随机数。现在我们产生n个0，1均匀分布随机数，我们有 为方便起见，我们特别选 n 12，则 这样我们很方便地就把标准正态分布随机数计算出来了。

　　4、中心极限定理说明，无穷多个具有相同分布的随机变量之和是服从标准正态分布的。相应的，如果多个随机变量均为0,1区间的均匀分布，则生成标准正态分布的式子为其中，k 为0,1区间均匀分布的随机变量的个数。k 值可以根据计算精度来选取，若取 k12，则。此时算得的将不超过-6,6的范围，这对于一般可靠度的计算已达足够精度。四使用Box Muller方法先证明令，则令，则有。接下来再来得出Box Muller方法设为一对相互独立的服从正态分布的随机变量。则有概率密度函数令,其中，则有分布函数令如果服从均匀分布，则的分布函数即为。最后，可以用代替，令为，其中，，得从而，只需要有两个服从均匀分布的随机变量，就能通过公式来得到一个服从正态分布的随机变量。用Box Muller方法来生成服从正态分布的随机数是十分快捷方便的。五、Hasiting 有理逼近法这是一种计算速度快，也能满足一定精度的算法。我们可以构造分布函数反函数的近似逼近公式，来产生标准正态分布的随机数。其计算公式为这里 ，系数为 a0 2.515517 b1 1.432788 a1 0.802853 b2 0.189269 a2 0.010328 b3 0.001308六、舍选法1、基本思想实质上是从许多RND随机数中选出一部分,使之成为具有给定分布的随机数。2、设随机变量X的概率密度函数为fx，存在实数ab，使PaXb1。算法步骤 1 选取常数，使fx1，xa,b； 2 产生两个RND 随机数r1 、r2，令y abar1 ； 3 若r2fy，则令xy；否则剔除 r1和r2,重返步骤2，重复循环,产生的随机数x1，x2，，xN的分布由概率函数 fx 确定。3、舍选法算法原理分析设PaZb1，Z的概率密度为fz,（1）选常数，使fz1，za，b；（2）随机变量X1，X2相互独立XiU0,1,令Y1abaX1Ua,b;（3）若X2fY1，则令X Y1，否则剔除X1，X2重复到2；则随机变量X的分布与Z相同。注可选取有限区间a1,b1，使得 （是很小的正数）例如，取 a13,b13，有 在区间a1,b1上应用舍选法,不会出现较大的系统误差。七、其他抽样方法直接抽样法、变换抽样值序抽样 【参考文献】【】魏宗舒；概率论与数理统计M。北京高等教育出版社,1995【】马砚儒.经验数学期望及性质。内蒙古民族大学学报,2002 【】程维虎，杨振海；舍选法几何解释及曲边梯形概率密度随机数生产算法J。数理统计与管理，2006,25【】杨振海，程维虎统计模拟。数理统计与管理，【】杨自强,魏公毅.产生伪随机数的若干新方法.数值计算机应用,2001【6】王萍,许海洋.一种新的随机数组合发生器的研究 计算机技术与发展，2006