C++求泊松分布

最新推荐文章于 2023-08-07 00:02:21 发布

西蒙情人

最新推荐文章于 2023-08-07 00:02:21 发布

阅读量2.7k

点赞数 10

分类专栏：数学文章标签： c++ 泊松分布

本文链接：https://blog.csdn.net/weixin_44458385/article/details/103265234

版权

数学专栏收录该内容

2 篇文章 0 订阅

订阅专栏

C++实现泊松分布

泊松分布公式

$P(X=k)=\frac{e^{-\lambda}\lambda^{k}}{k!}$
问题：直接用公式求会发生溢出

解决方法

方法一用递推计算泊松分布值

$P(0)=exp(-\lambda)$
$P(k+1)=\frac{\lambda P(k)}{k+1}$

代码大致为：

double result=exp(-lambda);
for(int i=1;i<k;i++)
	result*=(lambda/double(i))

方法二 C++11自带std::poisson_distribution

#include <iostream>
#include <iomanip>
#include <string>
#include <map>
#include <random>
int main()
{
    std::random_device rd;
    std::mt19937 gen(rd());
 
    // if an event occurs 4 times a minute on average
    // how often is it that it occurs n times in one minute?
    std::poisson_distribution<> d(4);
 
    std::map<int, int> hist;
    for(int n=0; n<10000; ++n) {
        ++hist[d(gen)];
    }
    for(auto p : hist) {
        std::cout << p.first <<
                ' ' << std::string(p.second/100, '*') << '\n';
    }
}