import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import beta
# 定义Beta分布的参数
alpha = 2
beta_param = 3
rng = 100
# 生成Beta分布的随机数
random_numbers = beta.rvs(alpha, beta_param, size=1000)*rng
# 绘制直方图
plt.hist(random_numbers, bins=30, density=True, alpha=0.6, color='g')
# 绘制Beta分布的理论曲线
x = np.linspace(0, rng, 1000)
pdf = beta.pdf(x/rng, alpha, beta_param)
plt.plot(x, pdf/rng, 'r-', lw=2)
plt.title('Beta Distribution (α=2, β=5)')
plt.xlabel('Value')
plt.ylabel('Density')
plt.show()

893

被折叠的 条评论
为什么被折叠?



