非常感谢楼主分享,大赞
我想问下,function [Tau,I_sq]=mutual_information(data,tau_max,n)中,I_sq应该是不能小于零的,但是我测试的时候出现了小于零,后来我发现,是不是需要将一下代码改下:
if Ps(i)~=0
H_s=H_s-Ps(i)*log(Ps(i));
elseif Pq(i)~=0
H_q=H_q-Pq(i)*log(Pq(i));
end
改成:
if Ps(i)~=0
H_s=H_s-Ps(i)*log(Ps(i));
end
if Pq(i)~=0
H_q=H_q-Pq(i)*log(Pq(i));
end