【注】学习本篇前最好一定的数学基础,即学习过概率论与数理统计
这个问题是求全概率的,在这里我们给它加上一问求:取得正品且是甲厂生产的概率。
import numpy as np
H = np.array([1/2, 3/10, 2/10]) # 完备事件组概率
C = np.array([9/10, 14/15, 19/20]) # 条件概率
h = np.array(1/2)
c = np.array(9/10)
# 全概率公式
def quangailv(H, C):
return (H * C).sum()
# 贝叶斯公式
def beiyesi(h, c, quangailv):
x = h*c
y = quangailv(H,C)
return x/y
print(quangailv(H, C))
print(beiyesi(h, c, quangailv))
输出结果为: