caientropy.m

来自「彩票预测 可以帮你预测彩票 自己做的 不是很好 请大家支持」· M 代码 · 共 14 行

function y=caientropy(cai)
[m,n]=size(cai);
h=zeros(30,1);
for i=1:m
    for j=1:n
        if cai(i,j)==0
        else
        h(cai(i,j))=h(cai(i,j))+1;
        end
    end
end
h1=h/(m*n);

y=-sum(h1.*log(h1+(h1==0)));