Spec-Zone.ru › NumPy 1.19

numpy.random.RandomState.standard_cauchy

метод

RandomState.standard_cauchy(size=None)

Генерирует выборку из стандартного распределения Коши с модой = 0.

Также известно как распределение Лоренца.

Примечание

Новый код должен использовать метод standard_cauchy экземпляра default_rng(); см. random-quick-start.

Параметры
sizeint or tuple of ints, optional

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

Возвращает
samplesndarray or scalar

Сгенерированные выборки.

См. также

Generator.standard_cauchy

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

Примечания

Функция плотности вероятности для полного распределения Коши:

P(x; x_0, \gamma) = \frac{1}{\pi \gamma \bigl[ 1+
(\frac{x-x_0}{\gamma})^2 \bigr] }

а стандартное распределение Коши просто устанавливает x_0=0 и \gamma=1.

Распределение Коши возникает при решении задачи о вынужденном гармоническом осцилляторе, а также описывает уширение спектральных линий. Оно также описывает распределение значений, в которых линия, наклоненная под случайным углом, пересекает ось x.

При изучении гипотез, которые предполагают нормальность, наблюдение за тем, как тесты работают с данными из распределения Коши, является хорошим индикатором их чувствительности к распределению с тяжелыми хвостами, так как распределение Коши очень похоже на распределение Гаусса, но с более тяжелыми хвостами.

Ссылки

1

Справочник NIST/SEMATECH по статистическим методам, «Распределение Коши», https://www.itl.nist.gov/div898/handbook/eda/section3/eda3663.htm

2

Weisstein, Эрик У. «Распределение Коши». Из MathWorld — Wolfram Web Resource. http://mathworld.wolfram.com/CauchyDistribution.html

3

Википедия, «Распределение Коши» https://en.wikipedia.org/wiki/Cauchy_distribution

Примеры

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

>>> import matplotlib.pyplot as plt
>>> s = np.random.standard_cauchy(1000000)
>>> s = s[(s>-25) & (s<25)]  # truncate distribution so it plots well
>>> plt.hist(s, bins=100)
>>> plt.show()
../../../_images/numpy-random-RandomState-standard_cauchy-1.png

© 2005–2020 NumPy Developers
Licensed under the 3-clause BSD License.
https://numpy.org/doc/1.19/reference/random/generated/numpy.random.RandomState.standard_cauchy.html

Spec-Zone.ru

Настройки Оффлайн Что нового Помощь О нас
Spec-Zone .ru
спецификации, руководства, описания, API