其他
100种分析思维模型之:大数定理
import numpy as np
import scipy.stats as st
import matplotlib as mpl
import matplotlib.pyplot as plt
# 防止中文乱码
mpl.rcParams['font.sans-serif']=[u'simHei']
mpl.rcParams['axes.unicode_minus']=False
# 定义次数和范围
n = 1000
x = np.arange(1, n+1, 1)
# 定义数组:出现正面的比例
y = []
# 循环计算
for i in x:
data = st.bernoulli.rvs(size=i, p=0.5) # 伯努利分布数据
y.append(sum(data)/len(data)) # 追加出现正面的比例
# 设置图表的大小
plt.figure(figsize=(18.8, 8))
# 作图设置
plt.xlabel('抛硬币次数', fontsize=20)
plt.ylabel('出现正面的比例', fontsize=20)
plt.tick_params(axis='both', which='major', labelsize=16)
# 作图并展示
plt.plot(x,y)
plt.plot([n,0.5], [0.5,0.5], 'r') # 均值参考线
plt.show()
↑ 点击此卡片关注 ↑