entropy.m

%entropy源程序

function h=entropy(x,n)
%entropy计算一个矩阵的熵的一阶估计
%返回n个像素的矩阵x的一阶估计(若n忽略，则默认为256)。
%这个估计是假设x中每个元素发生的概率是独立的。

error(nargchk(1,2,nargin));%检查输入参数
if nargin<2%对n设置默认
    n=256;
end

x=double(x);%把输入参数变成double型
xh=hist(x(:),n);%计算n个bin直方图
xh=xh/sum(xh(:));%计算概率

i=find(xh);
h=-sum(xh(i).*log2(xh(i)));%计算熵