c - c中(0,1)的随机数

Question

我想要一个指数分布的变量，这样做似乎非常有效:

double exp_variate(double lambda){
    return -log(my_random_function())/lambda;
}

其中 my_random_function 给出双 x，所以 0 < x < 1 (严格)意味着我永远不会得到 -log(0)或 0 .我如何获得不包括 0 和 1 值的随机数。我敢肯定之前有人会问过这个问题，但我看不到。此外，我对其他生成指数分布值的方法持开放态度。

编辑:> 到 <

Answer 1

大概你的意思是0 < x < 1 ，因为一个数不能同时小于 0 和大于 1。

在那种情况下，我可能会这样做:

double my_random_function() { 
    return (rand() + 1.0) / (RAND_MAX+2.0);
}

C 的随机数生成器通常生成从 0 到 RAND_MAX 的数字。这会将随机数生成器生成的值加一，因此它不能生成 0。除以 RAND_MAX+1 将得到从 1/RAND_MAX 到 (RAND_MAX+1)/(RAND_MAX+1) 的结果，因此后者是 1。因为我们想避免这种情况，所以我们使用 (RAND_MAX+2)，所以我们通过 (RAND_MAX+1)/(RAND_MAX+2) 得到 1.0/(RAND_MAX+2)，这将略大于 0.0 到略小于 1.0。