联合直方图，联合熵和互信息

Question

我有问题。我实现了任何向量的熵计算代码。我需要计算两个向量之间的互信息。我可以使用这个公式：H（x）+ H（y）-H（x，y）或H（x）-H（x | y）。我附上了我的代码。我的数据没有整数值（值为double）。我想采取两个不同向量的联合直方图。

for i = 1:14
    for j = 0:1023

        histogram =hist(data(i,j+1:Windows+j),bin);
        pdf = histogram./sum(histogram);

        Temp = -sum(pdf .* log2(pdf));

        entropyData(i,j+1) = Temp;
    end
end

如果我的代码适用于计算每个向量的熵，我的问题是如何计算联合直方图或pdf来找到联合熵？