python中numpy.random.gamma()函数

python中numpy.random.gamma()函数

根据官网，地址如下：
https://www.w3cschool.cn/doc_numpy_1_11/numpy_1_11-generated-numpy-random-gamma.html?lang=en

#需要导入的库
import numpy as np
import matplotlib.pyplot as plt
import scipy.special as sps

shape,scale=2.,2.
s=np.random.gamma(shape,scale,1000)
count, bins, ignored = plt.hist(s, 50, normed=True)#50：是50个条形图
'''
根据伽马分布的概率密度
伽马分布常被用来模拟电子元件的失效时间，
并且在泊松分布事件之间的等待时间相关的过程中自然产生。
'''
y=bins**(shape-1)*(np.exp(-bins/scale)/(sps.gamma(shape)*scale**shape))
plt.plot(bins,y,linewidth=2,color='r')
plt.show()

运行以上代码如下：


感悟：
“无论哪个行业，老手讲质量，新手拼价格，遇到便宜的不是你捡到宝了，而是被新手练手”这句话是在一个老师里口中说的，觉得很经典就偷偷的记下来，谢谢老师！