numpy.random.RandomState.exponential#

方法

random.RandomState.exponential(scale=1.0, size=None)#

從指數分佈中抽取樣本。

其機率密度函數為

\[f(x; \frac{1}{\beta}) = \frac{1}{\beta} \exp(-\frac{x}{\beta}),\]

對於 x > 0,其他地方為 0。\(\beta\) 是尺度參數,它是速率參數 \(\lambda = 1/\beta\) 的倒數。速率參數是指數分佈的另一種廣泛使用的參數化 [3]

指數分佈是幾何分佈的連續類比。它描述了許多常見情況,例如多次暴雨中測量的雨滴大小 [1],或維基百科的頁面請求之間的時間 [2]

注意

新程式碼應使用 exponential 方法,此方法屬於 Generator 實例;請參閱快速入門

參數:
scalefloat 或 float 的 array_like

尺度參數 \(\beta = 1/\lambda\)。必須是非負數。

sizeint 或 int 的 tuple,選用

輸出形狀。如果給定的形狀是,例如 (m, n, k),則會抽取 m * n * k 個樣本。如果 size 是 None (預設值),如果 scale 是純量,則傳回單一值。否則,會抽取 np.array(scale).size 個樣本。

傳回:
outndarray 或 純量

從參數化的指數分佈中抽取的樣本。

另請參閱

random.Generator.exponential

新程式碼應使用此方法。

參考文獻

[1]

Peyton Z. Peebles Jr., “Probability, Random Variables and Random Signal Principles”, 4th ed, 2001, p. 57.

[2]

Wikipedia, “Poisson process”, https://en.wikipedia.org/wiki/Poisson_process

[3]

Wikipedia, “Exponential distribution”, https://en.wikipedia.org/wiki/Exponential_distribution

範例

真實世界範例:假設一家公司有 10000 名客戶支援專員,且客戶來電之間的平均時間為 4 分鐘。

>>> n = 10000
>>> time_between_calls = np.random.default_rng().exponential(scale=4, size=n)

客戶在接下來 4 到 5 分鐘內來電的機率是多少?

>>> x = ((time_between_calls < 5).sum())/n 
>>> y = ((time_between_calls < 4).sum())/n
>>> x-y
0.08 # may vary