正确生成随机指数值c++

Question

我想生成属于指数分布的随机数。这是我写的

    int size = atoi(argv[2]);
    double *values = (double*)malloc(sizeof(double)*size);

    double gamma = atof(argv[1]);
    if(gamma<=0.0){
        cout<<"Insert gamma"<<endl;
        return 0;
    }


    for(int i=0; i<size; i++){
        values[i]=0;

    }

    srand ( time(NULL) );
    for(int i=0; i<size; i++){
        x = ((double) rand() / (RAND_MAX));
        //cout << random <<endl;
        value=(log(1.0-x)/(-gamma));
        //count each value
        values[value]=values[value]+1.0;
    }

但它们并没有涵盖所有矢量的大小。它们或多或少地覆盖了矢量的 10%，其他字段都是 0，并且由于在我需要进行线性插值之后我想减少那些 'empty space' 以便至少有一个值数组的每个单元格，我该怎么做？例如，我有一个 100000 的向量，只有前 60 个字段填充了值，因此从 60 到 999999 的单元格都是 0，当我进行线性回归时，它们会对公式产生负面影响。

Answer 1

好的，我看到了错误

您产生了 size 次事件。您确实需要更多事件来填充直方图

PS

填充 bin #n（n 在 [0...size）范围内）的概率由表达式

给出

prob = exp(-gamma*n) - exp(-gamma*(n+1))

其中 gamma 等于 0.01，比如说 n 大约 1000 会给你大约 4*10^-7 的概率。因此，即使要在这个容器中获得一个事件，您也需要采样大约 250 万次

PPS

并且使用库指数抽样虽然总体上很好，但不会给您带来任何好处，因为据我所知您的抽样还可以

正确生成随机指数值c++

Generate random exponential value correctly c++

c++

random

linear-regression

exponential-distribution