随机数产生和检验

1、关于随机数的产生与检验1第一张，PPT共四十二页，创作于2022年6月第一节 概论1. 意义：由于在统计学的不同技术中需要使用随机数，比如从统计总体中随机抽取样本时，或者在将实验动物随机分配到不同的试验组的过程中，或者在进行蒙特卡罗模拟法计算的时候等等，所以。 2.定义：设随机变量XF(x),则称随机变量X的抽样序列Xi为分布F(x)的随机数。（一）基本概念和定理2第二张，PPT共四十二页，创作于2022年6月定理定理1.1：设 是连续且严格单调上升的分布函数，它的反函数存在，且记为 ，即 2、若随机变量 ，则 的分布函数为 1、若随机变量 的分布函数为 ，则 3第三张，PPT共四十二页，创作

2、于2022年6月证明：设随机变量 的分布函数为 ，当 时，当 时， ；当 时，所以设 的分布函数为 ，则因为 ，对任意 有 。所以 的分布函数为 4第四张，PPT共四十二页，创作于2022年6月定理1.1说明了任意分布的随机数均可由均匀分布 的随机数变换得到。常简称 的随机数为均匀分布随机数。5第五张，PPT共四十二页，创作于2022年6月手工方法：抽签、掷骰子、摇号等；随机数表法：占用内存大，目前已很少使用；物理方法：放射性衰变、电子设备的热噪音、宇宙射线的触发时间等等；不能重复计算；数学方法：使用最广。（二）产生随机数的一般方法：6第六张，PPT共四十二页，创作于2022年6月（三）伪随机

3、数伪随机数：在计算机上用数学方法产生均匀随机数是指按照一定的计算方法而产生的数列，它们具有类似于均匀随机变量的独立抽样序列的性质，这些数既然是依照确定算法产生的，便不可能是真正的随机数，因此常把用数学方法产生的随机数称为伪随机数。伪随机数不可能真随机；需要对产生的伪随机数进行各种检验保证其符合独立性条件且分布为要求的分布；7第七张，PPT共四十二页，创作于2022年6月8第八张，PPT共四十二页，创作于2022年6月均匀随机数的产生：主要有线性同余法（LCG），组合同余法，反馈位移寄存器方法等 第二节 均匀随机数的产生9第九张，PPT共四十二页，创作于2022年6月同余性质：对称性：ab(mo

4、d M)，则ba(mod M).传递性:若ab(mod M)，bc(mod M)，则ac(mod M). （一）同余与线性同余法10第十张，PPT共四十二页，创作于2022年6月性质4：例如：已知1260(mod 16)，M=16，取C=6，a=2，b=10，因为(M,C)=2，则有210(mod 8)，其中M/(M,C)=16/2=8。或者，取C=12， M=16，因为(M,C)=4，则有15(mod 4)，其中M/(M,C)=16/4=4。11第十一张，PPT共四十二页，创作于2022年6月求余运算求余运算的式子A(mod M)定义为：其中 表示求 的整数部分。 12第十二张，PPT共四十

5、二页，创作于2022年6月线性同余法（Linear Congruence Generator，LCG）的递推公式为：13第十三张，PPT共四十二页，创作于2022年6月线性同余法的周期：14第十四张，PPT共四十二页，创作于2022年6月线性同余法产生的序列 一定会重复，因为周期最多只有M个可能取值。15第十五张，PPT共四十二页，创作于2022年6月说明：满周期是T=M时。16满周期第十六张，PPT共四十二页，创作于2022年6月当c0时，下式称为混合同余发生器，当c=0时，称为乘同余发生器，此时当模为素数时，称它为素数模乘同余发生器。 17补充1：混合同余发生器与素数模乘同余发生器第十七张

6、，PPT共四十二页，创作于2022年6月两个常用的混合同余发生器：18第十八张，PPT共四十二页，创作于2022年6月常用的素数模乘同余发生器 ：19第十九张，PPT共四十二页，创作于2022年6月常用的素数模乘同余发生器 ：20第二十张，PPT共四十二页，创作于2022年6月思想： 先用一个随机数发生器产生的随机数列为基础，再用另一个发生器对随机数列进行重新排列得到的新数列作为实际使用的随机数。这种把多个独立的发生器以某种方式组合在一起作为实际使用的随机数，希望能够比任何一个单独的随机数发生器得到周期长、统计性质更优的随机数，即组合发生器。 21补充2：组合发生器 ：第二十一张，PPT共四十

7、二页，创作于2022年6月 Maclaren 和 Marsaglia在1965年提出的著名的组合发生器是组合同余发生器,该算法的具体步骤如下： 22第二十二张，PPT共四十二页，创作于2022年6月2 .用第二个LCG产生一个随机整数 ,要求 ;3. 令 ，然后再用第一个LCG产生一个随机数 ,令 ；置 ;4 .重复23，得随机数列 ，即为组合同余发生器产生的数列。若第一个LCG的模为 ，令 ，则 为均匀随机数。1.用第一个LCG产生 个随机数，一般取 。这 个随机数被顺序地存放在矢量 中。置 ;23步骤：第二十三张，PPT共四十二页，创作于2022年6月检验目的：检验均匀伪随机数符合独立同均

8、匀分布；两种检验方法统计检验：对生成的伪随机数进行假设检验理论检验：从理论上讨论随机数发生器性质统计检验常用近似正态统计量和2统计量以下检验方法一般假设用某发生器生成了均匀分布伪随机数r1，r2，.，rn，来检验这些生成的随机数的各种统计量。 24第三节 随机数检验第二十四张，PPT共四十二页，创作于2022年6月251、特征量检验（参数检验） 第二十五张，PPT共四十二页，创作于2022年6月26第二十六张，PPT共四十二页，创作于2022年6月27注：若卡方值过大，则拒绝原假设（即分布不是均匀的）：检验随机数在（0,1）区间内分布时均匀的（一）卡方检验法：第二十七张，PPT共四十二页，创作

9、于2022年6月（二） Kolmogorov-Smirnov testK-S检验是连续分布的拟合性检验。检验样本的经验分布函数与总体的分布函数间的差异是否显著。28第二十八张，PPT共四十二页，创作于2022年6月29注：R软件检验随机数是否服从某一分布时，可采用这种检验方法。第二十九张，PPT共四十二页，创作于2022年6月3、独立性检验：自相关系数的检验 30 随机数r1,r2,rn中的前后项是否是统计相关性是否是显著的。相关系数反映了数据间的线性相关程度，若独立，则相关系数必为0（反之不一定）。 原假设H0：第三十张，PPT共四十二页，创作于2022年6月31第三十一张，PPT共四十二页

10、，创作于2022年6月R require sample(1:100, 20)#从1到100中无重复抽取20个数；runif(n, min=0, max=1)#产生n个0-1的均匀分布随机数；rnorm (n, mean = 0, sd = 1) #产生n个以0为均值，1为方差的正态分布随机数；rexp： The Exponential Distribution (wiki link) (指数分布，独立随机事件发生的时间间隔)rf：The F Distribution (wiki link) (F分布，两个卡方分布除以各自自由度)rbeta： The Beta Distribution (wik

11、i link)rbinom： The Binomial Distribution (wiki link) (二项分布)rcauchy： The Cauchy Distribution (wiki link) (柯西分布，N阶矩都不存在的分布.)rchisq： The (non-central) Chi-Squared Distribution (wiki link) (卡方分布，正态分布平方的分布)32第三十二张，PPT共四十二页，创作于2022年6月rgamma： The Gamma Distribution (wiki link) (伽玛分布)rpois： The Poisson Dist

12、ribution (wiki link) (泊松分布，单位时间内随机事件发生的次数)rgeom： The Geometric Distribution (wiki link) (几何分布，在第n次伯努利试验中，试验k次才得到第一次成功的机率)rhyper： The Hypergeometric Distribution (wiki link) (超几何分布)rlnorm： The Log Normal Distribution (wiki link) (对数正态分布，正态分布的指数的分布)rlogis： The Logistic Distribution (wiki link) (逻辑分布)r

13、multinom： The Multinomial Distribution (wiki link) (多变量正态分布)rnbinom： The Negative Binomial Distribution (wiki link) (负二项分布)33第三十三张，PPT共四十二页，创作于2022年6月R ks.test分布检验ks.test(data, pnorm, mean(data), sd(data)34第三十四张，PPT共四十二页，创作于2022年6月35第三十五张，PPT共四十二页，创作于2022年6月假设检验Matlab-require%M-file函数f的定义：判断概率函数func

14、tion f=p_judge(A,alpha)% 判别所给数据源在置信率为0.05时的概率分布形式。A=A(:);%数据集A的形式为n1。randperm(n)%产生1到n的均匀分布随机序列a=normrnd(0,1,1,6)%正态分布随机数36第三十六张，PPT共四十二页，创作于2022年6月正态分布mu,sigma=normfit(A);p1=normcdf(A,mu,sigma);H1,s1=kstest(A,A,p1,alpha)n=length(A); if H1=0disp(该数据源服从正态分布。)elsedisp(该数据源不服从正态分布。) end37第三十七张，PPT共四十二页

15、，创作于2022年6月Gama分布phat=gamfit(A,alpha);p2=gamcdf(A,phat(1),phat(2);H2,s2=kstest(A,A,p2,alpha) if H2=0disp(该数据源服从分布。)elsedisp(该数据源不服从分布。)end38第三十八张，PPT共四十二页，创作于2022年6月泊松分布lamda=poissfit(A,alpha);p3=poisscdf(A,lamda);H3,s3=kstest(A,A,p3,alpha)if H3=0disp(该数据源服从泊松分布。) elsedisp(该数据源不服从泊松分布。) end39第三十九张，PPT共四十二页，创作于2022年6月指数分布mu=expfit(A,alpha);p4=expcdf(A,mu);H4,s4=kstest(A,A,p4,alpha)if H4=0disp(该数据源服从指数分布。) elsedisp(该数据源不服从指数分布。