python (SimPy) 生成服从 erlang 分布的随机数

Question

我正在使用 Python(主要是 SimPy 包，但它与我认为的问题无关)，对一些系统进行建模并运行模拟。为此，我需要生成服从分布的随机数。到目前为止，我通过导入随机数(例如来自随机导入 *)并使用 expovariate 或 normalvariate 方法对一些分布(如指数分布和正态分布)做得很好。但是，我无法随机找到任何生成遵循 Erlang 分布的数字的方法。所以:

1. 有没有我忽略的方法？
2. 我是否必须导入其他一些库？
3. 我可以做一些解决方法吗？ (我认为我可以使用指数分布来生成随机的“Erlang”数字，但我不确定如何做。一段代码可能对我有帮助。

提前致谢!

Answer 1

Erlang 分布是 Gamma 分布的一个特例，它以 numpy.random.gamma (reference) 的形式存在。只需为 k(“形状”)参数使用整数值。另请参阅 scipy.stats.gamma，了解 PDF、CDF 等的功能。