任意分布的随机数的产生方法—VC程序实现方法

任意分布的随机数的产生方法—VC程序实现方法2010-11-02 vckbase 戎亚新摘要：

随机数在实际运用中非常之多，如游戏设计，信号处理，通常我们很容易得到平均分布的随机数。但如何根据平均分布的随机数进而产生其它分布的随机数呢？本文提出了一种基于几何直观面积的方法，以正态分布随机数的产生为例讨论了任意分布的随机数的产生方法。

正文：

一、平均分布随机数的产生

大家都知道，随机数在各个方面都有很大的作用，在vc的环境下，为我们提供了库函数rand（）来产生一个随机的整数。该随机数是平均在0~RAND_MAX之间平均分布的，RAND_MAX是一个常量，在VC6.0环境下是这样定义的：

#define RAND_MAX 0x7fff

它是一个short 型数据的最大值，如果要产生一个浮点型的随机数，可以将rand（）/1000.0这样就得到一个0~32.767之间平均分布的随机浮点数。如果要使得范围大一点，那么可以通过产生几个随机数的线性组合来实现任意范围内的平均分布的随机数。例如要产生-1000~1000之间的精度为四位小数的平均分布的随机数可以这样来实现。先产生一个0到10000之间的随机整数。方法如下：

int a = rand（）%10000;

然后保留四位小数产生0~1之间的随机小数：

double b = （double）a/10000.0;

然后通过线性组合就可以实现任意范围内的随机数的产生，要实现-1000~1000内的平均分布的随机数可以这样做：

double dValue = （rand（）%10000）/10000.0*1000-（rand（）%10000）/10000.0*1000;

则dValue就是所要的值。

到现在为止，你或许以为一切工作都已经完成了，其实不然，仔细一看，你会发现有问题的，上面的式子化简后就变为：

double dValue = （rand（）%10000）/10.0-（rand（）%10000）/10.0;

这样一来，产生的随机数范围是正确的，但是精度不正确了，变成了只有一位正确的小数的随机数了，后面三位的小数都是零，显然不是我们要求的，什么原因呢，又怎么办呢。

先找原因，rand（）产生的随机数分辨率为32767，两个就是65534，而经过求余后分辨度还要减小为10000，两个就是20000而要求的分辨率为1000*10000*2=20000000，显然远远不够。下面提供的方法可以实现正确的结果：

double a = （rand（）%10000） * （rand（）%1000）/10000.0;
double b = （rand（）%10000） * （rand（）%1000）/10000.0;
double dValue = a-b;

则dValue就是所要求的结果。在下面的函数中可以实现产生一个在一个区间之内的平均分布的随机数，精度是4位小数。

double AverageRandom（double min,double max）
{
int minInteger = （int）（min*10000）;
int maxInteger = （int）（max*10000）;
int randInteger = rand（）*rand（）;
int diffInteger = maxInteger - minInteger;
int resultInteger = randInteger % diffInteger + minInteger;
return resultInteger/10000.0;
}