创新互联Python教程:超几何分布python代码实现

超几何分布属于离散型随机变量的概率分布问题,随机变量可以取有限个值,在每取一个值时可以求出一个概率,此时求解的方法就是采用古典概型公式。

产品抽样检查中经常遇到一类实际问题,假定在N件产品中有M件不合格品,即不合格率  。

在产品中随机抽n件做检查,发现k件不合格品的概率为   ,k=0,1,2,...,min{n,M}。

相关推荐:《python基础教程》

亦可写作 

(与上式不同的是M可为任意实数,而C表示的组合数M为非负整数)

  为古典概型的组合形式,a为下限,b为上限,此时我们称随机变量X服从超几何分布(hypergeometric distribution)。

需要注意的是:

(1)超几何分布的模型是不放回抽样。

(2)超几何分布中的参数是M,N,n,上述超几何分布记作X~H(n,N,M)。

python实例代码:

import numpy as np
import matplotlib.pyplot as plt
# 超几何分布 hypergeometric(ngood, nbad, nsample, size=None) 好的总数、坏的总数、每次采样数、试验次数
# np.random.hypergeometric(10,20,5,size=4)
# 在一个口袋中装有30个球,其中有10个红球,其余为白球,这些球除颜色外完全相同。游戏者一次从中摸出5个球。
摸到至少4个红球就中一等奖,那么获一等奖的概率是多少?
s = np.random.hypergeometric(10,20,5,size=1000000)
p = sum(s>=4)/1000000.
print(p)
fig = plt.figure(figsize=(8,6))
a1 = fig.add_subplot(2,2,1)
a1.hist(s ,bins=20,color='k',alpha=0.3)
plt.show()

新闻标题:创新互联Python教程:超几何分布python代码实现
网站路径:http://www.shufengxianlan.com/qtweb/news46/318496.html

网站建设、网络推广公司-创新互联,是专注品牌与效果的网站制作,网络营销seo公司;服务项目有等

广告

声明:本网站发布的内容(图片、视频和文字)以用户投稿、用户转载内容为主,如果涉及侵权请尽快告知,我们将会在第一时间删除。文章观点不代表本网站立场,如需处理请联系客服。电话:028-86922220;邮箱:631063699@qq.com。内容未经允许不得转载,或转载时需注明来源: 创新互联