import numpy as np
import matplotlib.pyplot as plt
import math
#Show Normal Distribution/正态分布展示
def demo1(sample_number):
mu ,sigma = 0, 1
np.random.seed(0)
s = np.random.normal(mu, sigma, sample_number)
plt.hist(s, bins=100, density = True)
plt.xlabel(‘z’)
plt.ylabel(‘density’)
plt.title(‘The Normal Standard Distribution’)
plt.show()
if name == ‘main’:
demo1(1000000)
#Proof of the Center Limit Theory/中心极限定理证明
u = 0
sig = math.sqrt(0.5) # δ
t = 10000
a = np.zeros(10000)
for i in range(t):
a += np.random.uniform(-5, 5, 10000)
a = (a - u)/(sig/math.sqrt(t))
plt.hist(a, bins=100, color=‘g’, alpha=0.75)
plt.grid()
plt.show()
以上程序参考网上资料, 谢谢。
Thanks to the contributor on the website.