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