numpy.random.vonmises

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

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

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

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

Примечание

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

Параметры

mufloat or array_like of floats

Мода (“центр”) распределения.

kappafloat or array_like of floats

Дисперсия распределения, должна быть ≥0.

sizeint or tuple of ints, optional

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

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

outndarray or scalar

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

См. также

scipy.stats.vonmises

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

Generator.vonmises

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

Примечания

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

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

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

Ссылки

1

Абрамович, М. и Стегун, И. А. (ред.). «Справочник по математическим функциям с формулами, графиками и математическими таблицами, 9-е издание», Нью-Йорк: Dover, 1972.

2

фон Мизес, Р., «Математическая теория вероятности и статистики», Нью-Йорк: 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, density=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()