python

2024-06-02 14:08| 来源: 网络整理| 查看: 265

一、泊松分布问题：

假设我每天接到骚扰电话的次数服从泊松分布，并且经统计平均每天我会接到20个骚扰电话。请问： 1、我明天接到15个骚扰电话的概率？ 2、我明天接到24个骚扰电话以下的概率（包含24）？

二、泊松分布公式：

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

泊松分布的概率函数为：

$P(X=k)=\frac{\lambda ^{k}}{k!}e^{-\lambda },k=0,1,2,....$

累积概率分布函数为：

$P(X\leq x)=\sum_{k=0}^{x}\frac{\lambda ^{k}e^{-\lambda }}{k!}$

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

均值方差：泊松分布的均值和方差都是 $\lambda$ 。

针对上面的问1这里 $\lambda$ 就是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附近，概率最大，均值两边概率呈递减状态。再者可以看到这结果很像二项分布，也像正态分布，关于两者三个分布的相似，这里不做讨论。

五、查表

查表前面代码中有涉及，其实就是查累积概率。查询 $\lambda$ =100，发生次数小于等于120的概率：

from scipy import stats p = stats.poisson.cdf(120, 100) print(p)

结果：0.9773306709216473

六、随机数生成

生成服从 $\lambda$ =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]

【本文地址】

公司简介

联系我们