在Python中生成Gamma分布的随机数

在Python中，我们可以使用numpy.random模块的gamma函数来生成Gamma分布的随机数。Gamma分布是一种连续概率分布，经常用于描述等待时间、生存时间等。它的概率密度函数为：

f(x) = x^(k-1) * e^(-x/θ) / (θ^k * Γ(k))

其中，k为形状参数，θ为尺度参数，Γ为伽玛函数。

以下是一个例子，展示如何在Python中生成Gamma分布的随机数：

import numpy as np
import matplotlib.pyplot as plt

# 设置参数
k = 2    # 形状参数
theta = 1    # 尺度参数

# 生成随机数
random_numbers = np.random.gamma(k, theta, 1000)

# 绘制直方图
plt.hist(random_numbers, bins=30, density=True, alpha=0.6, color='g', edgecolor='black')

# 绘制概率密度函数曲线
x = np.linspace(0, 10, 1000)
y = x**(k-1) * np.exp(-x/theta) / (theta**k * np.math.gamma(k))
plt.plot(x, y, 'r-', label='gamma pdf')

# 显示图例和标题
plt.legend()
plt.title('Gamma Distribution')

# 显示图像
plt.show()

在上面的例子中，我们使用了numpy.random.gamma函数生成1000个随机数，参数k设置为2，参数theta设置为1。然后，使用matplotlib.pyplot模块绘制了随机数的直方图，并画出了Gamma分布的概率密度函数曲线。

你可以根据自己的需求，调整k和theta的值，生成不同形状和尺度的Gamma分布的随机数。同时，你也可以调整bins的值，改变直方图的精度。