标准库rand()函数的缺陷以及Blitz++随机数生成的简介

简介: 所以在这里我将简单介绍一下如何使用Blitz++库中的随机数生成类。不过在这里我不能够证明Blitz++随机数算法相对于标准库有什么优越性。Blitz++的源代码是开放的,你可以完全了解它的随机数算法的实现。

当我们需要在某个任务中使用随机数,通常我们习惯于使用标准库的rand函数。像这样:srand(time(0)); // 时间种子

        rand() % MAX_RAND 

标准库的rand函数使用线性同余算法,是生成速度相当快的一种随机数生成算法。在多数情况下也确实能满足我们的要求,但是对于一些特殊目的应用这个算法生成的随机数是不行的,比如某些加密算法,蒙特卡罗积分等(在.NET中创建随机密码的加密安全随机数就不能使用Random类的线性同余随机数,而要使用System.Security.Cryptography命名空间中的相关随机数生成类)。

这个线性同余算法的实现可以在很多书籍中找到。下面我给出一个The C Programming Langurage 2nd中的一个实现,这也是普遍使用的标准库随机数算法的实现:

   unsigned long int next =1;

 

   /* rand:  return pseudo-random integer on 0..32767 */

    int rand(void)

   {

       next = next *1103515245+12345;

       return ( unsigned int )(next/65536)%32768;

   }

 

   /* srand:  set seed for rand() */

   void srand(unsignedint seed)

   {

       next = seed;

   }

 

这个实现的问题在于rand函数return行中的那个32768,在标准库中这个数字定义为RAND_MAX宏,在VisualC++和Mingw32编译器的stdlib.h头文件(或者cstdlib)中你都可以发现RAND_MAX的值为32768。也就是说这个算法的随机数分布在0--RAND_MAX中,而在一般编译器中就是0--32768。假设你的算法需要的是300000多个的随机数,那么使用rand函数会产生重负次数近30次的随机数!

所以在这里我将简单介绍一下如何使用Blitz++库中的随机数生成类。不过在这里我不能够证明Blitz++随机数算法相对于标准库有什么优越性。Blitz++的源代码是开放的,你可以完全了解它的随机数算法的实现

在Blitz++中随机数类组织在ranlib namespace中,使用方法也非常简单,seed成员函数种入种子后,再用random成员函数就可以了。Blitz++中包括了产生各种概率分布情况的随机数,列举如下:均匀分布(Uniform),正态分布(Normal),指数分布(Exponential),Beta分布,Gamma分布,Χ方分布,F分布,离散均匀分布。具体地可以参考Blitz++文档的第九章。本文将会演示正态分布的随机数类。

				NormalUnit<>()   
				标准正态分布,
				μ
				 = 0, 
				σ
				 = 1;
				 
		
				Normal<>(T mean, T standardDeviation) 
				正态分布,N(
				μ
				, 
				σ
				^2
				)
				,其中mean就是
				μ
				,standardDeviation
				就是
				σ
				
		

Blitz++ 中随机数类的seed是共享的,调用一个类的seed后,其他类也同样种入了相同的种子。对于种子的来源,Blitz++文档中有这样一个注意事项:

Note: you may be tempted to seed the random number generator from a static initializer. Don't do it! Due to an oddity of C++, there is no guarantee on the order of static initialization when templates are involved. Hence, you may seed the RNG before its constructor is invoked, in which case your program will crash. If you don't know what a static initializer is, don't worry -- you're safe!

(幸运的是我也不太清楚static initializer具体指什么,:-) )

1 ,先来看一个标准正态分布的例子,根据3 σ 法则(正态分布的随机数几乎全部落在( μ -3 σ , μ +3 σ ) 中),这些随机数应该大部分落在(-3, 3)之间。

下面的程序产生10000个正态分布随机数,然后导入Matlab生成图形,横坐标是随机数的序数,纵坐标是随机数值。很显然,大部分点在3 σ区间内。在区间 ( μ - σ , μ + σ)中数据点最密。

#include <random/uniform.h>

#include <random/normal.h>

#include <iostream>

#include <fstream>

#include <ctime>

#include <cstdlib>

using namespace std;

using namespace ranlib;

 

int main()

{

    const size_t MAXGEN =10000;

    NormalUnit<float> rnd;

 

    rnd.seed(time(0));

 

    ofstream file("c:\\file.txt");

 

    // generator Normal distribution random number

    for (size_t i=0; i<MAXGEN;++i)

    {

        file << rnd.random()<<" ";

    }

}

2  下面是一个服从于 μ= 10,σ= 2的正态分布,根据3σ法则,大部分点应该落在(4,16)之间。

#include <random/uniform.h>

#include <random/normal.h>

#include <iostream>

#include <fstream>

#include <ctime>

#include <cstdlib>

using namespace std;

using namespace ranlib;

 

int main()

{

    const size_t MAXGEN =1000000;

    const int mu =10;

    const int sigma =2;

    Normal < float > rnd(mu,sigma);

 

    rnd.seed(time(0));

 

    ofstream file("c:\\file.txt");

 

    // generator Normal distribution random number

    for (size_t i=0; i<MAXGEN;++i)

    {

        file << rnd.random()<<" ";

    }

}

 

3,生成前述正态分布的钟型曲线(PDF)。

这里产生1M的float随机数,然后乘以10转型为整数,统计在区间0-170之间随机数出现的概率,这样再导入Matlab生成图形就是一个近似的正态分布的钟型曲线了。

 

 

#include <random/uniform.h>

#include <random/normal.h>

#include <iostream>

#include <fstream>

#include <ctime>

#include <cstdlib>

using namespace std;

using namespace ranlib;

 

int main()

{

    const size_t MAXGEN =1000000;

    const int mu =10;

    const int sigma =2;

    Normal < float > rnd(mu,sigma);

 

    rnd.seed(time(0));

 

    ofstream file("c:\\file.txt");

 

    float *rndArray =newfloat[MAXGEN];

    // generator Normal distribution random number

    for (size_t i=0; i<MAXGEN;++i)

    {

        rndArray[i]= rnd.random();

    }

 

    int *rndConvertIntArray =newint[MAXGEN];

    int multiple =10;

    // convert float random number to integer

    for (size_t i=0; i<MAXGEN;++i)

    {

        rndConvertIntArray[i]=int(rndArray[i]* multiple);

    }

 

    const size_t PDFSIZE =(mu +3* sigma)* multiple;

    const size_t redundancy =10;

    int *pdf =newint[PDFSIZE+redundancy];

    for (size_t i=0; i<PDFSIZE+redundancy;++i)

    {

        pdf[i]=0;

    }

 

    // generator PDF(Probability Distribution Function), Normal distribution is a "bell" shape

    for (size_t i=0; i<MAXGEN;++i)

    {

         if (rndConvertIntArray[i]<= PDFSIZE+redundancy-1)

        {

            ++pdf[rndConvertIntArray[i]];

        }

    }

 

    for (size_t i=0; i<PDFSIZE+redundancy;++i)

    {

        file << pdf[i]  <<" ";

    }

 

    delete [] rndArray;

    delete [] rndConvertIntArray;

    delete [] pdf;

}

 

钟型曲线,当然不是特别光滑 ( ), 大部分随机数都出现在( 4  16 )区间内。

 

总结, Blitz++ 的随机数生成类应该说是值得信任的。

目录
相关文章
|
5月前
Verilog基础:$random系统函数的使用
Verilog基础:$random系统函数的使用
105 0
|
5月前
|
资源调度 算法 JavaScript
Python基础专题 - 超级详细的 Random(随机)原理解析与编程实践
Python基础专题 - 超级详细的 Random(随机)原理解析与编程实践
395 0
|
4月前
|
Python
python之数值计算、math库、随机数
python之数值计算、math库、随机数
|
5月前
|
资源调度 Python
|
5月前
|
Python
Python中的math和cmath模块:数学运算的得力助手
Python作为一种功能强大的编程语言,提供了丰富的数学运算功能。其中,math和cmath模块就是Python中用于数学运算的重要工具。math模块提供了基本的数学函数和常量,适用于实数运算;而cmath模块则提供了对复数运算的支持,使得Python在数学计算和工程应用中更加灵活和强大。
|
5月前
|
Rust 数据安全/隐私保护
rust每日一库 rand 生成随机数
rust每日一库 rand 生成随机数
154 0
|
Java
如何使用Java实现随机数生成器
在Java编程中,需要生成随机数的情况非常常见。本文将介绍如何使用Java中提供的相关类和方法来实现随机数生成器。
93 1
|
数据挖掘 Python
Python基础 | 你想要的随机数生成都在这里
Python基础 | 你想要的随机数生成都在这里
74 0
|
Go
Go语言 rand(随机数)包
Go语言 rand(随机数)包
157 0
|
编解码 安全 算法
随机数探秘|如果python不用random怎样生成随机数?
随机数探秘|如果python不用random怎样生成随机数?
502 0