numpy.random.vonmises#
- random.vonmises(mu, kappa, size=None)#
从 von Mises 分布中抽取样本。
样本从指定模式 (mu) 和浓度 (kappa) 的 von Mises 分布中抽取，区间为 [-pi, pi]。
von Mises 分布（也称为圆形正态分布）是单位圆上的连续概率分布。它可以被认为是正态分布的圆形模拟。
- 参数:
- mufloat 或 float 数组类型
分布的模式（“中心”）。
- kappafloat 或 float 数组类型
分布的浓度，必须 >= 0。
- sizeint 或 int 元组，可选
输出形状。如果给定的形状是，例如，
(m, n, k)，那么抽取m * n * k个样本。如果 size 是None（默认），如果mu和kappa都是标量，则返回单个值。否则，抽取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. (编). “数学函数手册，包含公式、图表和数学表格，第 9 次印刷”，纽约：多佛出版社，1972 年。
[2]von Mises, R., “概率和统计的数学理论”，纽约：学术出版社，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()
