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';} }

總結

以上是生活随笔為你收集整理的C++求泊松分布的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。