高斯随机数产生原理及代码_笔记.pdfVIP

高斯随机数产生原理及代码_笔记--第1页

一:随机数的产生

C++中不提供random函数,但是提供了rand函数,产生0~RAND_MAX之间的整数,但严格意义上来

讲生成的只是伪随机数(pseudo-randomintegralnumber).生成随机数时需要指定一个种子，如果

在程序内循环，那么下一次生成随机数时调用上一次的结果作为种子。但如果分两次执行程序，那么由于

种子相同，生成的“随机数”也是相同的。rand()函数不接受参数，默认以1为种子（即起始值）。若

随机数生成器总是以相同的种子开始，所以形成的伪随机数列也相同，失去了随机意义。（但这样便于程

序调试）.

C++中另一函数srand（），可以指定不同的数（无符号整数变元）为种子。但是如果种子相同，伪

随机数列也相同。一个办法是让用户输入种子，但是仍然不理想。比较理想的是用变化的数，比如时间来

作为随机数生成器的种子。time的值每时每刻都不同。所以种子不同，所以，产生的随机数也不同。

rand函数产生随机数的方法:

1如果要产生0~10的10个整数，可以表达为：

intN=rand()%11;

这样，N的值就是一个0~10的随机数，如果要产生1~10，则是这样：

intN=1+rand()%11;

总结来说，可以表示为：

a+rand()%(b-a+1)

其中的a是起始值，(b-a+1)是整数的范围。

若要0~1的小数，则可以先取得0~10的整数，然后均除以10即可得到随机到十分位的10个随机小数，

若要得到随机到百分位的随机小数，则需要先得到0~100的10个整数，然后均除以100，其它情况依

此类推。当要求的精度较高的时候，使用RAND_MAX作为分母。如果要求左闭右开的话，分母设置为

RAND_MAX+1即可。精度要求高时形式如下：

x=a+((rand()%RAND_MAX)/(double)RAND_MAX)*(b-a);/*xbelong[a,b]*/

通常rand()产生的随机数在每次运行的时候都是与上一次相同的，这是有意这样设计的，是为了便

于程序的调试。若要产生每次不同的随机数，可以使用srand(seed)函数进行随机化，随着seed的不

同，就能够产生不同的随机数。还可以包含time.h头文件，然后使用srand(time(0))来使用当前时

间使随机数发生器随机化，这样就可以保证每两次运行时可以得到不同的随机数序列(只要两次运行的间隔

超过1秒)。后一种方法性能更优！

在工程应用时，我们一般将系统当前时间(Unix时间)作为种子，这样生成的随机数更接近于实际意

义上的随机数。给一下例程如下：

#includeiostream

#includetime.h

usingnamespacestd;

voidmain(){

inti=0;

intcount[6]={0};

srand(time(0));/*根据系统时间设置随机数种子,