第三章随机数的产生与检验演示文稿

第三章随机数的产生与检验演示文稿1现在是1页\一共有37页\编辑于星期一优选第三章随机数的产生与检验现在是2页\一共有37页\编辑于星期一定理定理1.1：设是连续且严格单调上升的分布函数，它的反函数存在，且记为，即2、若随机变量，则的分布函数为1、若随机变量的分布函数为，则3现在是3页\一共有37页\编辑于星期一证明：设随机变量的分布函数为，当时，当时，；当时，所以设的分布函数为，则因为，对任意有。所以的分布函数为4现在是4页\一共有37页\编辑于星期一定理1.1说明了任意分布的随机数均可由均匀分布的随机数变换得到。常简称的随机数为均匀分布随机数。5现在是5页\一共有37页\编辑于星期一手工方法：抽签、掷骰子、摇号等；随机数表法：占用内存大，目前已很少使用；物理方法：放射性衰变、电子设备的热噪音、宇宙射线的触发时间等等；不能重复计算；数学方法：使用最广。（二）产生随机数的一般方法：6现在是6页\一共有37页\编辑于星期一（三）伪随机数伪随机数：在计算机上用数学方法产生均匀随机数是指按照一定的计算方法而产生的数列，它们具有类似于均匀随机变量的独立抽样序列的性质，这些数既然是依照确定算法产生的，便不可能是真正的随机数，因此常把用数学方法产生的随机数称为伪随机数。伪随机数不可能真随机；需要对产生的伪随机数进行各种检验保证其符合独立性条件且分布为要求的分布；7现在是7页\一共有37页\编辑于星期一8现在是8页\一共有37页\编辑于星期一均匀随机数的产生：主要有线性同余法（LCG），组合同余法，反馈位移寄存器方法等

第二节均匀随机数的产生9现在是9页\一共有37页\编辑于星期一同余性质：对称性：a≡b(modM)，则b≡a(modM).传递性:若a≡b(modM)，b≡c(modM)，则a≡c(modM).

（一）同余与线性同余法10现在是10页\一共有37页\编辑于星期一性质4：例如：已知12≡60(mod16)，M=16，取C=6，a=2，b=10，因为(M,C)=2，则有2≡10(mod8)，其中M/(M,C)=16/2=8。或者，取C=12，M=16，因为(M,C)=4，则有1≡5(mod4)，其中M/(M,C)=16/4=4。11现在是11页\一共有37页\编辑于星期一求余运算求余运算的式子A(modM)定义为：其中表示求的整数部分。12现在是12页\一共有37页\编辑于星期一线性同余法（LinearCongruenceGenerator，LCG）的递推公式为：13现在是13页\一共有37页\编辑于星期一线性同余法的周期：14现在是14页\一共有37页\编辑于星期一线性同余法产生的序列一定会重复，因为周期最多只有M个可能取值。15现在是15页\一共有37页\编辑于星期一说明：满周期是T=M时。16满周期现在是16页\一共有37页\编辑于星期一当c≠0时，下式称为混合同余发生器，当c=0时，称为乘同余发生器，此时当模为素数时，称它为素数模乘同余发生器。

17补充1：混合同余发生器与素数模乘同余发生器现在是17页\一共有37页\编辑于星期一两个常用的混合同余发生器：18现在是18页\一共有37页\编辑于星期一常用的素数模乘同余发生器：19现在是19页\一共有37页\编辑于星期一常用的素数模乘同余发生器：20现在是20页\一共有37页\编辑于星期一思想：

先用一个随机数发生器产生的随机数列为基础，再用另一个发生器对随机数列进行重新排列得到的新数列作为实际使用的随机数。这种把多个独立的发生器以某种方式组合在一起作为实际使用的随机数，希望能够比任何一个单独的随机数发生器得到周期长、统计性质更优的随机数，即组合发生器。

21补充2：组合发生器：现在是21页\一共有37页\编辑于星期一

Maclaren和Marsaglia在1965年提出的著名的组合发生器是组合同余发生器,该算法的具体步骤如下：22现在是22页\一共有37页\编辑于星期一2.用第二个LCG产生一个随机整数,要求;3.令，然后再用第一个LCG产生一个随机数,令；置;4.重复2~3，得随机数列，即为组合同余发生器产生的数列。若第一个LCG的模为，令，则为均匀随机数。1.用第一个LCG产生个随机数，一般取。这个随机数被顺序地存放在矢量中。置;23步骤：现在是23页\一共有37页\编辑于星期一检验目的：检验均匀伪随机数符合独立同均匀分布；两种检验方法统计检验：对生成的伪随机数进行假设检验理论检验：从理论上讨论随机数发生器性质统计检验常用近似正态统计量和ⅹ2统计量以下检验方法一般假设用某发生器生成了均匀分布伪随机数r1，r2，...，rn，来检验这些生成的随机数的各种统计量。24第三节随机数检验现在是24页\一共有37页\编辑于星期一251、特征量检验（参数检验）现在是25页\一共有37页\编辑于星期一26现在是26页\一共有37页\编辑于星期一27注：若卡方值过大，则拒绝原假设（即分布不是均匀的）：检验随机数在（0,1）区间内分布时均匀的（一）卡方检验法：现在是27页\一共有37页\编辑于星期一（二）Kolmogorov-SmirnovtestK-S检验是连续分布的拟合性检验。检验样本的经验分布函数与总体的分布函数间的差异是否显著。28现在是28页\一共有37页\编辑于星期一29注：R软件检验随机数是否服从某一分布时，可采用这种检验方法。现在是29页\一共有37页\编辑于星期一3、独立性检验：自相关系数的检验30随机数r1,r2,……,rn中的前后项是否是统计相关性是否是显著的。相关系数反映了数据间的线性相关程度，若独立，则相关系数必为0（反之不一定）。

原假设H0：现在是30页\一共有37页\编辑于星期一31现在是31页\一共有37页\编辑于星期一Rrequiresample(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：TheExponentialDistribution(wikilink)(指数分布，独立随机事件发生的时间间隔)rf：TheFDistribution(wikilink)(F分布，两个卡方分布除以各自自由度)rbeta：TheBetaDistribution(wikilink)rbinom：TheBinomialDistribution(wikilink)(二项分布)rcauchy：TheCauchyDistribution(wikilink)(柯西分布，N阶矩都不存在的分布...)rchisq：The(non-central)Chi-SquaredDistribution(wikilink)(卡方分布，正态分布平方的分布)32现在是32页\一共有37页\编辑于星期一rgamma：TheGammaDistribution(wikilink)(伽玛分布)rpois：ThePoissonDistribution(wikilink)(泊松分布，单位时间内随机事件发生的次数)rgeom：TheGeometricDistribution(wikilink)(几何分布，在第n次伯努利试验中，试验k次才得到第一次成功的机率)rhyper：TheHypergeometricDistribution(wikilink)(超几何分布)rlnorm：TheLogNormalDistribution(wikilink)(对数正态分布，正态分布的指数的分布)rlogis：TheLogisticDistribution(wikilink)(逻辑分布)rmultinom：TheMultinomialDistribution(wikilink)(多变量正态分布)rnbinom：TheNegativeBinomialDistribution(wikilink)(负二项分布)

33现在是33页\一共有37页\编辑于星期一Rks.test分布检验ks.test(data,"pnorm",mean(data),sd(data))34现在是34页\一共有37页\编辑于星期一35现在是35页\一共有37页\编辑于星期一假设检验Matlab-require%M-file函数f的定义：判断概率函数functionf=p_judge(A,alpha)%判别所给数据源在置信率为0.05时的概率分布形式。

A=A(:);%数据集A的形式为n×1。

randperm(n)%产生1到n的均匀分布随机序列a=normr