python实现概率质量函数Poisson

Poisson分布Python实现

最新推荐文章于 2024-10-10 17:34:13 发布

saber艾瑞莉娅

最新推荐文章于 2024-10-10 17:34:13 发布

阅读量2.2k

点赞数 1

CC 4.0 BY-SA版权

分类专栏：概率论文章标签：概率质量函数 python

本文链接：https://blog.youkuaiyun.com/qq_36627030/article/details/79394220

概率论专栏收录该内容

1 篇文章

订阅专栏

本文介绍了使用Python实现Poisson概率质量函数的两种方法。第一种方法利用了scipy库中的poisson模块，第二种方法则通过自定义函数实现。同时展示了不同参数下Poisson分布的概率质量函数图像。

用python实现概率质量函数Possion

原理如下：

方法一：用scipy库的poisson

import scipy.stats as stats

plt.figure(figsize=(12.5,4))
a = np.arange(16)
poi = stats.poisson
lambda_ = [1.5,4.25]
colors = ['#348ABD','#A60628']

plt.bar(a,poi.pmf(a,lambda_[0])#第一个参数是K，第二个参数是lambda的值

,color = colors[0],label="$\lambda = %.lf$"%lambda_[0],

alpha=0.60,edgecolor=colors[0],lw="3")

plt.bar(a,poi.pmf(a,lambda_[1]),color = colors[1],label="$\lambda = %.lf$"%lambda_[1],
alpha=0.60,edgecolor=colors[1],lw="3")

第二种方法是我自己写了Possion函数，结果跟上面一样：

def poissons(k,lambdas):
    p = np.zeros((1,len(k)))
    for i in range(len(k)):
        if i == 0:
            p[0,i] = math.pow(lambdas,k[i]) * math.exp(-lambdas)
        else:
            sum_k = 1
            a = k[i].astype('float64')
            while a != 0:
                sum_k = sum_k * a
                a = a - 1
            p[0,i] = math.pow(lambdas,k[i]) * math.exp(-lambdas) / sum_k
    return p