limentr.m

fs1016w编码标准matlab源码

function H=limentr(Aprob,Aprobtrans)
%求极限熵
H=0;
l=size(Aprob,2);
%for i=1:l
%H=-sum(Aprob(i)*Aprobtrans(:,i).*log(Aprobtrans(:,i)));
%end

for i=1:l
    for j=1:l
        if ~Aprobtrans(i,j)==0
        H=H-Aprob(i)*Aprobtrans(i,j)*log2(Aprobtrans(i,j));
    end
    end
end