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/archive/2012/07/15/2592713.html

欢迎关注

微信关注下方公众号,第一时间获取干货硬货;公众号内回复【pdf】免费获取数百本计算机经典书籍

原创文章受到原创版权保护。转载请注明出处:https://www.ccppcoding.com/archives/55300

非原创文章文中已经注明原地址,如有侵权,联系删除

关注公众号【高性能架构探索】,第一时间获取最新文章

转载文章受原作者版权保护。转载请注明原作者出处!

(0)
上一篇 2023年2月9日 上午6:24
下一篇 2023年2月9日 上午6:24

相关推荐