生成指数分布的随机数

Question

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

double ran_expo(lambda){
    double u;

    u = rand() / (RAND_MAX + 1.0);

    return -log(1- u) / lambda;
}

我正在使用它(我没有在这里放置代码的其他部分)来生成指数分布的随机数lambda = 0.05.(概率密度函数是lambda * exp(-lambda * x)).

但是,我总是得到非常小的数字,例如 0.000041,或类似的东西1.#INF00(这是什么？).

实际上,对于指数分布lambda = 0.05,生成的数字应该非常大,即大多数都不是30.我的结果非常奇怪.另外,精度不尽如人意,仅为10 ^( - 6).我试过了long double,但它仍然是这样的.

我在windows下使用DEV C++.

Answer 1

这是因为您没有声明类型lambda,更正了结果在您寻求的范围内.如果未声明,旧编译器将认为它是int.

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

double ran_expo(double lambda){
    double u;
    u = rand() / (RAND_MAX + 1.0);
    return -log(1- u) / lambda;
}

int main(void)
{
    int i;
    srand((unsigned)time(NULL));
    for (i=0; i<20; i++)
        printf("%f\n", ran_expo(0.05));
    return 0;
}

节目输出:

0.025040
16.582459
4.296027
33.079902
17.589123
13.073084
8.624299
45.254803
34.611211
27.454302
3.825699
39.168172
24.790600
14.411160
7.247698
0.301951
1.917010
9.065004
3.187146
3.627885