numpy.random.vonmises

numpy.random.vonmises(mu, kappa, size=None)

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

Выборка осуществляется из распределения фон Мизеса с заданным режимом (mu) и дисперсией (kappa) на интервале [-π, π].

Распределение фон Мизеса (также известное как круговое нормальное распределение) — это непрерывное распределение вероятностей на единичной окружности. Его можно рассматривать как круговой аналог нормального распределения.

Параметры:	mu : float Режим (центр) распределения. kappa : float Дисперсия распределения, должна быть >=0. size : int или кортеж из целых чисел, необязательно Форма выходного массива. Если заданная форма, например, (m, n, k), то m * n * k выборок генерируется. По умолчанию None, в этом случае возвращается единственное значение.
Возвращает:	samples : скаляр или массив NumPy Сгенерированные выборки, которые находятся в интервале [-π, π].

См. также

scipy.stats.distributions.vonmises

функция плотности вероятности, распределение или кумулятивная функция распределения и т.д.

Примечания

Плотность вероятности распределения фон Мизеса:

где — это режим, — дисперсия, а — модифицированная функция Бесселя нулевого порядка.

Распределение фон Мизеса названо в честь Рихарда Эдлера фон Мизеса, родившегося в Австро-Венгрии, в том, что сейчас Украина. Он бежал в США в 1939 году и стал профессором Гарвардского университета. Он работал в области теории вероятностей, аэродинамики, механики жидкости и философии науки.

Ссылки

[R270]

Abramowitz, M. and Stegun, I. A. (Eds.). “Справочник по математическим функциям с формулами, графиками и математическими таблицами, 9-е издание,” Нью-Йорк: Dover, 1972.

[R271]

von Mises, R., “Математическая теория вероятностей и статистики”, Нью-Йорк: Academic Press, 1964.

Примеры

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

>>> mu, kappa = 0.0, 4.0 # mean and dispersion
>>> s = np.random.vonmises(mu, kappa, 1000)

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

>>> import matplotlib.pyplot as plt
>>> from scipy.special import i0
>>> plt.hist(s, 50, normed=True)
>>> x = np.linspace(-np.pi, np.pi, num=51)
>>> y = np.exp(kappa*np.cos(x-mu))/(2*np.pi*i0(kappa))
>>> plt.plot(x, y, linewidth=2, color='r')
>>> plt.show()

(Исходный код)