numpy.random.RandomState.vonmises#

方法

random.RandomState.vonmises(mu, kappa, size=None)#

从 von Mises 分布中抽取样本。

样本从具有指定模式 (mu) 和浓度 (kappa) 的 von Mises 分布中抽取,在区间 [-pi, pi] 上。

von Mises 分布(也称为圆形正态分布)是单位圆上的连续概率分布。它可以被认为是正态分布的圆形类似物。

注意

新代码应使用 vonmises 方法,该方法是 Generator 实例的;请参阅快速入门

参数:
mufloat 或 float 的类数组

分布的模式(“中心”)。

kappafloat 或 float 的类数组

分布的浓度,必须 >=0。

sizeint 或 int 元组,可选

输出形状。如果给定形状为,例如 (m, n, k),则抽取 m * n * k 个样本。如果 size 为 None(默认),如果 mukappa 均为标量,则返回单个值。否则,抽取 np.broadcast(mu, kappa).size 个样本。

返回:
outndarray 或标量

从参数化的 von Mises 分布中抽取的样本。

另请参阅

scipy.stats.vonmises

概率密度函数、分布或累积密度函数等。

random.Generator.vonmises

新代码应使用此方法。

注释

von Mises 分布的概率密度为

\[p(x) = \frac{e^{\kappa cos(x-\mu)}}{2\pi I_0(\kappa)},\]

其中 \(\mu\) 是模式,\(\kappa\) 是浓度,\(I_0(\kappa)\) 是 0 阶修正贝塞尔函数。

von Mises 以 Richard Edler von Mises 的名字命名,他出生于奥匈帝国,也就是现在的乌克兰。他于 1939 年逃往美国,并在哈佛大学担任教授。他研究概率论、空气动力学、流体力学和科学哲学。

参考文献

[1]

Abramowitz, M. 和 Stegun, I. A. (Eds.)。"数学函数手册,包含公式、图表和数学表格,第 9 次印刷",纽约:Dover, 1972。

[2]

von Mises, R.,“概率和统计的数学理论”,纽约:Academic Press,1964 年。

示例

从分布中抽取样本

>>> mu, kappa = 0.0, 4.0 # mean and concentration
>>> 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()
../../../_images/numpy-random-RandomState-vonmises-1.png