博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
boost random library的使用
阅读量:4941 次
发布时间:2019-06-11

本文共 1743 字,大约阅读时间需要 5 分钟。

 
生成满足一定分布的随机数,是统计模拟、系统仿真等应用中最基本的要求。matlab中提供了函数可以生成各种常见分布的随机数,c++使用boost random库也可以很容易实现。

一、例子

boost random库的文档提供了一个例子,模拟掷色子。投掷一个均匀的色子,六个面每个面出现的概率应该是相等的,也就是说,投掷出的点数应该是服从{ 1 2 3 4 5 6 } 上的均匀分布的。

#include <boost/random.hpp>
#include <boost/random.hpp>
#include <iostream>
int main()
{
   boost::mt19937gen;                                     
   boost::uniform_int<>dist(1,6);
   boost::variate_generator<boost::mt19937&,boost::uniform_int<>>die(gen,dist);
    for( int i=0; i<10; ++i) {
       std::cout<<die() <<std::endl;
   }
}

 

上面的例子虽然简单,但也可以看出要生成满足一定分布的伪随机数的一般步骤:

首先,我们需要一个伪随机数发生器(generator)。语句

boost::mt19937 gen;

声明了一个mt19937类型的伪随机数发生器gen。boost一共提供了17种类型伪随机发生器供使用者选择,它们是:

minstd_rand,  rand48, lrand48, ecuyer1988, kreutzer1986, hellekalek1995, mt11213b, mt19937, lagged_fibonacci 607, lagged_fibonacci1279, lagged_fibonacci2281, lagged_fibonacci3217, lagged_fibonacci4423, lagged_fibonacci9689, lagged_fibonacci19937, lagged_fibonacci23209, lagged_fibonacci44497。

总之,名字是一个比一个丑啊。到底选用何种发生器,需要使用者在性能、质量和内存之间进行权衡,一般的应用的话,boost文档推荐使用mt19937就够了。也可以选择多种发生器,使结果相互印证,从而更可信。

其次,需要一个分布函数。由伪随机数发生器生成的伪随机数一般是一个很大区间上的均匀分布(如[1,2^99999999])。而要产生我们想要的分布的随机数(本例中,我们想要的分布为 {1 2 3 4 5 6 }上的均匀分布),还要对它们做某些变换。

boost::uniform_int<> dist(1, 6);

就是干这个事的。

boost random提供了如下的分布函数:

uniform_smallint ------ 一个整数集上的离散均匀分布

uniform_int ------ 一个整数集上的离散均匀分布

uniform_01 ------ [0, 1)上的连续均匀分布

uniform_real ------ [ min, max)上的连续均匀分布

bernoulli_distribution ------ 伯努利分布

geometric_distribution ------ 几何分布

triangle_distribution ------ 三角分布

exponential_distribution ------ 指数分布

normal_distribution ------ 正态分布

lognormal_distribution ------ 对数正态分布

uniform_on_sphere ------ 球面上的均匀分布

 

 

最后,把以上两个结合生成我们想要的随机数发生器:

boost::variate_generator<boost::mt19937&, boost::uniform_int<> > die(gen, dist);

转载于:https://www.cnblogs.com/wjgaas/p/3848979.html

你可能感兴趣的文章
《第一行代码》学习笔记2-Android开发特色
查看>>
bzoj3396 [Usaco2009 Jan]Total flow 水流
查看>>
20165231 2017-2018-2 《Java程序设计》第3周学习总结
查看>>
(180905)如何通过梯度下降法降低损失----Google机器学习速成课程笔记
查看>>
(响应式PC端媒体查询)电脑屏幕分辨率尺寸大全
查看>>
Crystal Reports拉报表报错:Error detected by database DLL
查看>>
Java获取新浪微博cookies
查看>>
面试题总结
查看>>
【BZOJ1095】捉迷藏(动态点分治)
查看>>
Table Basics [转载]
查看>>
Logback 学习笔记
查看>>
并查集
查看>>
11、组件注册-使用FactoryBean注册组件
查看>>
nyoj_95_众数问题_map练习
查看>>
For循环
查看>>
020-安装centos6.5后的生命历程
查看>>
面试介绍项目经验(转)
查看>>
创建并设置ASP.NET的会话状态服务器(SQL State Server)
查看>>
<metro>Google的验证
查看>>
SQL中NUMERIC和DECIMAL的区别
查看>>