第2章 伪随机数产生.pptVIP

第二章 伪随机数的产生 第二章 伪随机数的产生 一.伪随机数产生的意义 二.产生U(0,1)的乘同余法 三.正态分布N(0,1)的产生 四.逆变法与其它分布随机数的产生 一.伪随机数产生的意义 在GA,SA,TS中都要用到； 在计算机中的固有伪随机数发生器只有U(0,1) 且可重复性不好，没有其他分布； 自己设计的发生器，可控型好、可重复性好，便于仿真比较。 二.产生U(0,1)的乘同余法（1） 乘同余法的计算公式 可产生随机数序列。 问题：怎样设定 和 可以使随机数序列最长？ 二.产生U(0,1)的乘同余法（2） 乘同余法的方法： 若 的整数，当 {x}满足以下条件时，可以达到最大周期 (序列长度) 为3(Mod8)或5 (Mod8)的数 ； 为奇数，一般取为1。 二.产生U(0,1)的乘同余法（3） 乘同余法举例说明： = =16 = 3 =1,3,9,11,1,3,9,11 = 5 =1,5,9,13,1,5,9,13 = 3 =2,6,2,6… 可得整数序列 ，要想获得U(0,1)，见下面 二.产生U(0,1)的乘同余法（4） 产生U(0,1)步骤： ； 令 。 二.产生U(0,1)的乘同余法（5） 产生U(0,1)举例说明： = =16 =3，x0=1 =1/16, 3/16, 9/16, 11/16, 1/16, 3/16, 9/16, 11/16 … =3，x0=2 =2/16, 6/16, 2/16, 6/16… 二.产生U(0,1)的乘同余法（6） 优秀编程举例： A=65539 M=2147483648=231 Xi=AXi IF(Xi0) Xi=Xi+M(M对应于计算机中最大整数) 三.正态分布N(0,1)的产生（1） 正态分布可以用多个U(0,1)来近似，若 是独立同分布， 较大，则 近 似正态分布，且满足 及 则 三.正态分布N(0,1)的产生（2） 令： 一般n取12则： 其中： （详见下页） 三.正态分布N(0,1)的产生（3） 注： 四.逆变法与其它分布随机数的产生（1） 逆变法 四.逆变法与其它分布随机数的产生（2） 是分布函数， ，如何产生 ？ 设 ， 是随机变量 产生 ， 是 分布函数 逆变法的目的：产生 分布的随机数 四.逆变法与其它分布随机数的产生（3） 逆变法的步骤： 已知 ，或由 求 即 ，令 推导 产生 用 得到 四.逆变法与其它分布随机数的产生（4） 负指数分布的产生 负指数函数的密度函数： 四.逆变法与其它分布随机数的产生（5） 负指数函数的分布函数的产生过程： ① 令 ② ③ 产生 则 ④ 即 四.逆变法与其它分布随机数的产生（6） 产生 是负指数分布的。 四.逆变法与其它分布随机数的产生（6） 思考1： 爱尔朗（Erlang）分布是m个负指数分布的和。设 为负指数分布，则 为Erlang分布。试设计其随机数的产生。 四.逆变法与其它分布随机数的产生（6） Answer： 产生 四.逆变法与其它分布随机数的产生（6） 思考2： 试用逆变法设计分布如图所示密度函数 f(x) 的伪随机数序列的产生方法。 * * 序列 满足以下关系式： 常数 取模(除以M后的余数) 或 取整 98.7 0 1 0 1 分布函数 1 0 1 密度函数 = 1 , 0≤x≤1 0 , 其它 1