import numpy as np import matplotlib.pyplot as plt from scipy.stats import beta

设置贝塔分布的参数

alpha = 2 beta = 5

生成贝塔分布的随机样本

samples = beta.rvs(alpha, beta, size=1000)

计算贝塔分布的概率密度函数

x = np.linspace(0, 1, 100) pdf = beta.pdf(x, alpha, beta)

绘制贝塔分布的随机样本直方图和概率密度函数图像

plt.hist(samples, density=True, alpha=0.5, label='Samples') plt.plot(x, pdf, 'r', label='PDF') plt.xlabel('x') plt.ylabel('Probability Density') plt.title('Beta Distribution (alpha=2, beta=5)') plt.legend() plt.show()