python 验证大数定理

大数定理简单来说是指实验次数越多 事件A发生的频数会接近事件A的期望值。或者说 事件A发生的频率接近于他的概率。

import numpy as np #导入numpy数学包
import matplotlib.pyplot as plt#导入pyolot画图工具包
balls=np.array([0,1,1,2,3])#0代表白球,1代表黑球,2代表蓝球,3代表红球 两个黑球 其他球各一个
N=2000 #试验次数
fre=np.zeros(shape=N) #频率结果数组
def get_frequence(all_random_number): #计算取出黑球频率
    n=len(all_random_number)#实验总次数
    frequence=(np.sum(all_random_number==1))
    return frequence/n#返回频率

for j in range(1,N+1):#循环 增加试验次数取球 每次试验次数+1
    result=np.zeros(shape=j)#结果数组
    for i in range(j):#取球实验 取J次球
        data=np.random.choice(balls,1)#随机取一个球
        result[i]=np.sum(data)#choice取出来是一个数组,要用sum转成整数
    fre[i]=get_frequence(result)#每次实验取i次球得到的概率放入频率结果数组
plt.rcParams['font.family'] = ['sans-serif']#让 plt显示汉字
plt.rcParams['font.sans-serif'] = ['SimHei']#让 plt显示汉字
plt.plot(range(1,N+1),fre)#绘制结果
plt.grid(ls='--',c='y')#设定坐标
plt.title('伯努利大数定律演示----抽取黑球的频率与抽取次数的关系')#标题
plt.xlabel('实验次数')#X轴标号
plt.ylabel('黑球出现频率')#Y轴标号
plt.yticks(ticks=np.arange(0.1,1,0.1))#设置Y轴间隔
plt.show()#显示图片