随机变量模拟
随机变量模拟代码实现
// 导入程序所需包
import numpy as np
import numpy.random as npr
import matplotlib.pyplot as plt
%matplotlib inline
//生成随机数,并将随机数分配到不同的变量中
sample_size=500
rn1=npr.rand(sample_size,3)
rn2=npr.randint(0,10,sample_size)
rn3=npr.sample(size=sample_size)
a=[0,25,50,75,100]
rn4=npr.choice(a,size=sample_size)
//生成可视化试图
fig,((ax1,ax2),(ax3,ax4))=plt.subplots(nrows=2,ncols=2,figsize=(7,7))
ax1.hist(rn1,bins=25,stacked=True)
ax1.set_title(‘rand’)
ax1.set_ylabel(‘frequency’)
ax2.hist(rn2,bins=25)
ax2.set_title(‘randint’)
ax2.grid(True)
ax3.hist(rn3,bins=25)
ax3.set_title(‘sample’)
ax3.set_ylabel(‘frequency’)
ax3.grid(True)
ax4.hist(rn4,bins=25)
ax4.set_title(‘choice’)
ax4.grid(True)