Question

我无法保持随机生成的值，这些值通常分布在0和1之间（包括0，不包括1）。我相信算法基本上是正确的，我只是难倒在这里。任何见解都会很棒。

这些是必需的包含文件：

#include <limits.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>

正态分布随机数生成器函数：

float rand_normal(float mean, float stddev) 
{
    static float n2 = 0.0;
    float x, y, r;
    static int n2_cached = 0;
    if (!n2_cached) 
    {
        do
        {
            x = 2.0*rand()/RAND_MAX - 1;
            y = 2.0*rand()/RAND_MAX - 1;
            r = x*x + y*y;
        } while (r==0.0 || r>1.0);
        float d = sqrt(-2.0*log(r)/r);
        float n1 = x*d;
        float result = n1*stddev + mean;
        n2 = y*d;
        n2_cached = 1;
        return result;
    }
    else
    {
        n2_cached = 0;
        return n2*stddev + mean;
    }
}

主要功能仅用于测试目的。

int main()
{
    srand(time(NULL));
    int i;
    float min = 0.5, max = 0.5, r, avg = 0;
    float x, w;
    int n = 10000000;
    for (i=0; i<n; i++)
    {
        r = rand_normal(0.5, 0.09);
        if (r < min)
            min = r;
        else if ( r>max)
            max = r;
        avg += r;
    }
    avg /= (float)n;
    printf("min = %f\nmax = %f\navg = %f\n", min, max, avg);
    return 0;
}

如果有人想知道，这个功能是“植物遗传继承”模拟所必需的。

Answer 1

为什么您希望结果保持在0到1之间？高斯分布具有完全支持，因此无论您正在查看的是什么区间，无论您选择何种均值和方差，总会有一个（可能非常小的）非零概率落在该区间之外。如果你真的想因某种原因将自己限制在[0,1]，那么你可以简单地调用rand_normal，直到你进入那段时间。

另请注意，虽然Box-Müller（您正在使用的算法）易于实现，但这是生成高斯随机变量的最差和最昂贵的方法之一。我所知道的最好和最快的算法是“Ziggurat”方法，其实现可以在

找到

http://www.seehuhn.de/pages/ziggurat

Answer 2

我肯定会创建一个将“rand（）”转换为标准化浮点值的函数。例如：

double
nrand ()
{
  return rand()/(RAND_MAX - 1);
}

此外，以下是一些可能有用的链接：

C中0到1 [0,1]之间的归一化随机数分布

2 个答案: