我有问题。我实现了任何向量的熵计算代码。我需要计算两个向量之间的互信息。我可以使用这个公式: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来找到联合熵?