python生成泊松分布随机数_python-介绍泊松分布(poisson分布)

一、泊松分布问题：

假设我每天接到骚扰电话的次数服从泊松分布，并且经统计平均每天我会接到20个骚扰电话。

请问：

1、我明天接到15个骚扰电话的概率？

2、我明天接到24个骚扰电话以下的概率(包含24)？

二、泊松分布公式：

首先要清楚，泊松分布是离散的，也就是说我接到骚扰电话次数必须是整数，要么就是15次，要么就是16次...,而不会是15.5次。

泊松分布的概率函数为：

累积概率分布函数为：

其实就是概率函数的累加。

均值方差：泊松分布的均值和方差都是

。

针对上面的问1这里

就是20，k就是15。

三、解决问题

方法一：根据公式编写代码

import math

def p_possion(k, m):

kjie = 1 #k!

for i in range(1, k+1):

kjie*=i

pk = math.pow(m, k)/kjie*math.e**(-m)

return pk

p = p_possion(15, 20)

print("接到15个骚扰电话的概率：",p)

p=0

for i in range(0,25):

p+=p_possion(i, 20)

print("接到24个骚扰电话以下的概率：",p)

结果：

接到15个骚扰电话的概率： 0.0516488535317584

接到24个骚扰电话以下的概率： 0.8432273781737633

方法二：用scipy包中的scipy.stats.poisson

from scipy import stats

p = stats.poisson.pmf(15, 20)

print("接到15个骚扰电话的概率：",p)

p = stats.poisson.cdf(24, 20)

print("接到24个骚扰电话以下的概率：",p)

四、泊松概率及累积概率分布

仍以骚扰电话为例，计算接到电话0~50次的概率：

from scipy import stats

import matplotlib.pyplot as plt

plt.rcParams['font.sans-serif']=['SimHei']#用来正常显示中文标签

plt.rcParams['axes.unicode_minus']=False#用来正常显示负号

X=range(0,51)

Y=[]

for k in X:

p = stats.poisson.pmf(k, 20)

Y.append(p)

plt.bar(X, Y, color="red")

plt.xlabel("次数")

plt.ylabel("概率")

plt.title("接到骚扰电话次数及概率")

plt.show()

结果：

首先可以看出，在均值20附近，概率最大，均值两边概率呈递减状态。再者可以看到这结果很像二项分布，也像正态分布，关于两者三个分布的相似，这里不做讨论。

五、查表

查表前面代码中有涉及，其实就是查累积概率。查询

=100，发生次数小于等于120的概率：

from scipy import stats

p = stats.poisson.cdf(120, 100)

print(p)

结果：0.9773306709216473

六、随机数生成

生成服从

=50的泊松分布随机数100个：

from scipy import stats

# 设置random_state时，每次生成的随机数一样。不设置或为None时，多次生成的随机数不一样

sample = stats.poisson.rvs(mu=50, size=100, random_state=3)

print(sample)

结果：

[51 45 60 40 34 53 54 45 45 49 51 46 48 61 47 53 47 48 45 49 52 45 43 50

50 54 54 47 47 46 36 72 54 55 52 37 42 41 54 54 55 58 53 53 51 43 58 38

63 50 44 53 48 43 53 45 67 37 51 42 54 47 59 55 54 55 55 46 60 43 54 45

59 44 58 45 51 58 56 47 54 33 55 50 58 49 60 37 51 43 50 52 52 45 42 44

49 54 52 48]

author：蓝何忠

email：[email protected]