matlab计算JS散度

计算 JS 散度需要先计算 KL 散度和 JSD 值，其中 KL 散度可以通过以下公式计算：

$KL(P||M) = \sum_{i=1}^{n}p_i log(\frac{p_i}{m_i})$

其中 $P$ 和 $M$ 分别表示两个概率分布，$p_i$ 和 $m_i$ 分别表示这两个分布在第 $i$ 个元素的概率。

然后，可以通过以下公式计算 JSD 值：

$JS(P||M) = \frac{1}{2}(KL(P||M) + KL(M||P))$

现在来看具体的 Matlab 代码：

function jsd = JSD(P, M)
    % 计算 KL 散度
    KL_PM = sum(P .* log2(P./M));
    KL_MP = sum(M .* log2(M./P));
    % 计算 JSD 值
    jsd = 0.5 * (KL_PM + KL_MP);
end

其中 $P$ 和 $M$ 都是概率分布向量，可以根据具体情况进行设置。