bash函数$ RANDOM应该具有统一分布吗?

时间:2013-04-11 01:57:57

标签: linux bash statistics

据我所知,bash函数$ RANDOM会在一个范围内生成随机整数,但是,这些数字是否应该遵循(或近似)一个统一的离散分布?

2 个答案:

答案 0 :(得分:5)

我刚打印$ RANDOM一百万次,把它变成直方图,用gnu​​meric查看它,图表显示非常正态的分布!

for n in `seq 1 1000000`; do echo $RANDOM ; done > random.txt
gawk '{b=int($1/100);a[b]++};END{for (n in a) {print n","a[n]}}' random.txt > hist.csv
gnumeric hist.csv

因此,如果您想要近似线性分布,请使用$(($ RANDOM%$ MAXIMUM))并且不要使用大于16383的$ MAXIMUM或8192安全。如果你想要非常大的数字,你可以多次连接$ RANDOM%1000,只要你处理前导零。

如果你想要正态分布,请使用$(($ RANGE * $ RANDOM / 32767 + $ MINIMUM)),并记住这只是整数数学。

答案 1 :(得分:3)

Bash document实际上并没有这么说:

  

RANDOM

     

每次引用此参数时,都会生成0到32767之间的随机整数。   为此变量赋值会使随机数生成器播种。

读到这一点,我肯定假设它意图是线性的;恕我直言,因为它不是别的。

但是看一下bash源代码,$RANDOM的实现旨在产生一个线性分布(这来自bash 4.2源代码):

/* The random number seed.  You can change this by setting RANDOM. */
static unsigned long rseed = 1;
static int last_random_value;
static int seeded_subshell = 0;

/* A linear congruential random number generator based on the example
   one in the ANSI C standard.  This one isn't very good, but a more
   complicated one is overkill. */

/* Returns a pseudo-random number between 0 and 32767. */
static int
brand ()
{
  /* From "Random number generators: good ones are hard to find",
     Park and Miller, Communications of the ACM, vol. 31, no. 10,
     October 1988, p. 1195. filtered through FreeBSD */
  long h, l;

  /* Can't seed with 0. */
  if (rseed == 0)
    rseed = 123459876;
  h = rseed / 127773;
  l = rseed % 127773;
  rseed = 16807 * l - 2836 * h;
#if 0
  if (rseed < 0)
    rseed += 0x7fffffff;
#endif
  return ((unsigned int)(rseed & 32767));       /* was % 32768 */
}

正如评论所暗示的那样,如果您想要好的随机数,请使用其他内容。