计算JS散度MATLAB

JS散度是一种用于度量两个概率分布之间距离的方法，可以通过以下公式进行计算：

JS(P,Q) = 0.5 * KL(P,M) + 0.5 * KL(Q,M)

其中，P和Q是两个概率分布，M是它们的平均分布，KL(P,M)和KL(Q,M)是它们与平均分布之间的KL散度。

下面是MATLAB代码实现：

function jsd = jsdiv(P, Q)
    M = 0.5 * (P + Q);
    jsd = 0.5 * kl_div(P, M) + 0.5 * kl_div(Q, M);
end

function kl = kl_div(P, Q)
    kl = sum(P .* log2(P ./ Q));
end

其中，P和Q是两个概率分布向量，可以通过MATLAB的hist函数计算得到。函数jsdiv返回P和Q的JS散度值，函数kl_div返回P和Q之间的KL散度值。