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

设置贝塔分布的参数

a = 2 b = 5

生成贝塔分布的随机样本

samples = beta.rvs(a, b, size=1000)

绘制贝塔分布的概率密度函数

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

plt.plot(x, pdf, 'r-', lw=2, label='Beta(' + str(a) + ', ' + str(b) + ') PDF') plt.hist(samples, bins=30, density=True, alpha=0.6, color='g', label='Samples')

plt.xlabel('x') plt.ylabel('Probability Density') plt.title('Beta Distribution Example') plt.legend() plt.show()