numpy.random.RandomState.poisson#
方法
- random.RandomState.poisson(lam=1.0, size=None)#
从泊松分布中抽取样本。
泊松分布是大 N 时二项分布的极限。
- 参数:
- lam浮点数或浮点数类型的类数组
固定时间间隔内发生的预期事件数,必须 >= 0。序列必须可广播到请求的大小。
- size整数或整数元组,可选
输出形状。如果给定的形状是,例如
(m, n, k)
,则抽取m * n * k
个样本。如果 size 为None
(默认值),如果lam
是标量,则返回单个值。否则,抽取np.array(lam).size
个样本。
- 返回:
- outndarray 或标量
从参数化的泊松分布中抽取的样本。
另请参阅
random.Generator.poisson
新代码应使用此方法。
注释
泊松分布的概率质量函数 (PMF) 为
\[f(k; \lambda)=\frac{\lambda^k e^{-\lambda}}{k!}\]对于预期间隔为 \(\lambda\) 的事件,泊松分布 \(f(k; \lambda)\) 描述在观察到的间隔 \(\lambda\) 内发生 \(k\) 个事件的概率。
由于输出仅限于 C int64 类型的范围,当 lam 在最大可表示值的 10 个标准差内时,会引发 ValueError。
参考文献
[1]Weisstein, Eric W. “泊松分布。”来自 MathWorld–A Wolfram Web Resource。https://mathworld.wolfram.com/PoissonDistribution.html
[2]维基百科,“泊松分布”, https://en.wikipedia.org/wiki/Poisson_distribution
示例
从分布中抽取样本
>>> import numpy as np >>> s = np.random.poisson(5, 10000)
显示样本的直方图
>>> import matplotlib.pyplot as plt >>> count, bins, ignored = plt.hist(s, 14, density=True) >>> plt.show()
为 lambda 100 和 500 各抽取 100 个值
>>> s = np.random.poisson(lam=(100., 500.), size=(100, 2))