numpy.random.weibull

random.weibull(a, size=None)

Генерация выборок из распределения Вейбулла.

Генерация выборок из однопараметрического распределения Вейбулла с заданным параметром формы a.

\[X = (-ln(U))^{1/a}\]

Здесь U — случайная величина, равномерно распределённая на интервале (0,1].

Более распространённое двухпараметрическое распределение Вейбулла, включающее параметр масштаба \(\lambda\), определяется как \(X = \lambda(-ln(U))^{1/a}\).

Примечание

Новый код должен использовать метод weibull экземпляра Generator; см. Быстрый старт.

Параметры:

afloat или array_like из float

Параметр формы распределения. Должен быть неотрицательным.

sizeint или кортеж из int, необязательно

Форма выходных данных. Если заданная форма, например, (m, n, k), то m * n * k выборок генерируется. Если size имеет значение None (по умолчанию), то возвращается одно значение, если a является скалярным. В противном случае генерируется np.array(a).size выборок.

Возвращаемое значение:

outndarray или скаляр

Выбранные значения из распределения Вейбулла.

См. также

scipy.stats.weibull_max

scipy.stats.weibull_min

scipy.stats.genextreme

gumbel

random.Generator.weibull

которую следует использовать в новом коде.

Примечания

Распределение Вейбулла (или распределение экстремальных значений типа III для наименьших значений, распределение Вейбулла типа III или распределение Розина—Раммлера) относится к классу обобщённых распределений экстремальных значений (ОЭЗ), используемых при моделировании задач экстремальных значений. Этот класс включает распределения Гумбеля и Фреше.

Плотность вероятности для распределения Вейбулла:

\[p(x) = \frac{a} {\lambda}(\frac{x}{\lambda})^{a-1}e^{-(x/\lambda)^a},\]

где \(a\) — параметр формы, а \(\lambda\) — параметр масштаба.

Функция имеет свой пик (моду) в \(\lambda(\frac{a-1}{a})^{1/a}\).

Когда a = 1, распределение Вейбулла сводится к экспоненциальному распределению.

Ссылки

[1]

Вейбулл В., Королевский технологический институт, Стокгольм, 1939 г. «Статистическая теория прочности материалов», Известия Шведской инженерной академии № 151, 1939 г., Литографическая типография Генерального штаба, Стокгольм.

[2]

Вейбулл В., «Статистическая функция распределения широкого применения», Журнал прикладной механики ASME, статья 1951 г.

[3]

Википедия, «Распределение Вейбулла», https://en.wikipedia.org/wiki/Weibull_distribution

Примеры

Генерация выборок из распределения:

>>> a = 5. # shape
>>> s = np.random.weibull(a, 1000)

Отображение гистограммы выборок вместе с функцией плотности вероятности:

>>> import matplotlib.pyplot as plt
>>> x = np.arange(1,100.)/50.
>>> def weib(x,n,a):
...     return (a / n) * (x / n)**(a - 1) * np.exp(-(x / n)**a)

>>> count, bins, ignored = plt.hist(np.random.weibull(5.,1000))
>>> x = np.arange(1,100.)/50.
>>> scale = count.max()/weib(x, 1., 5.).max()
>>> plt.plot(x, weib(x, 1., 5.)*scale)
>>> plt.show()