scipy.spatial.distance.

jensenshannon#

scipy.spatial.distance.jensenshannon(p, q, base=None, *, axis=0, keepdims=False)[source]#

计算两个概率数组之间的 Jensen-Shannon 距离（度量）。 这是 Jensen-Shannon 散度的平方根。

两个概率向量 p 和 q 之间的 Jensen-Shannon 距离定义为，

\[\sqrt{\frac{D(p \parallel m) + D(q \parallel m)}{2}}\]

其中 \(m\) 是 \(p\) 和 \(q\) 的逐点均值，而 \(D\) 是 Kullback-Leibler 散度。

如果 p 和 q 的总和不为 1.0，则此例程将对其进行归一化。

参数:
p(N,) array_like

左概率向量

q(N,) array_like

右概率向量

basedouble, 可选

用于计算输出的对数的底，如果没有给出，则例程使用 scipy.stats.entropy 的默认底。

axisint, 可选

计算 Jensen-Shannon 距离的轴。 默认为 0。

在版本 1.7.0 中添加。

keepdimsbool, 可选

如果将其设置为 True，则缩减的轴将作为大小为一的维度保留在结果中。 使用此选项，结果将与输入数组正确广播。 默认为 False。

在版本 1.7.0 中添加。

返回值:
jsdouble 或 ndarray

沿 axis 计算的 p 和 q 之间的 Jensen-Shannon 距离。

备注

在版本 1.2.0 中添加。

示例

>>> from scipy.spatial import distance
>>> import numpy as np
>>> distance.jensenshannon([1.0, 0.0, 0.0], [0.0, 1.0, 0.0], 2.0)
1.0
>>> distance.jensenshannon([1.0, 0.0], [0.5, 0.5])
0.46450140402245893
>>> distance.jensenshannon([1.0, 0.0, 0.0], [1.0, 0.0, 0.0])
0.0
>>> a = np.array([[1, 2, 3, 4],
...               [5, 6, 7, 8],
...               [9, 10, 11, 12]])
>>> b = np.array([[13, 14, 15, 16],
...               [17, 18, 19, 20],
...               [21, 22, 23, 24]])
>>> distance.jensenshannon(a, b, axis=0)
array([0.1954288, 0.1447697, 0.1138377, 0.0927636])
>>> distance.jensenshannon(a, b, axis=1)
array([0.1402339, 0.0399106, 0.0201815])
本页内容