C++生成随机数：Beta分布（beta distribution）

double cls_random::randomBeta(
    double alpha,
    double beta)
{
    /*Johnk's beta generator*/
    double u, v;
    double x, y;
    do 
    {
        u=cls_random::randomUniform();
        v=cls_random::randomUniform();
        x=pow(u,1/alpha);
        y=pow(v,1/beta);
    } while (x+y>1);
    return x/(x+y);
}

期望:\operatorname{E}[x] = \frac{\alpha}{\alpha+\beta}\!

方差:V=\frac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}\!

wiki： http://en.wikipedia.org/wiki/Beta_distribution

【推广】 免费学中医，健康全家人
原文地址：https://www.cnblogs.com/yeahgis/p/2592713.html