熵的MATLAB源程序

关于图像处理的 首先我把一个图像生成一个矩阵，然后要把这其中的数字算它的统计特性，根据熵的公式：算出信源熵 （首先算出矩阵中数字的出现概率，然后对概率取对数再乘以概率 对所以可能的概率都这样计算 最后连加取反） 我要的是MATLAB程序，希望谁帮忙 谢谢 急急

%^^^^^^^^^^^^^^^^^^ Get Entropy Of Image Y ^^^^^^^^^^^^^
p = imhist(Y(:));
% remove zero entries in p
p(p==0) = [];
% normalize p so that sum(p) is one.
p = p ./ numel(Y);
H(1,i) = -sum(p.*log2(p));

求助两次了。。抱歉啊，我不熟悉matlab，没有用过，帮你顶起，看朋友们谁熟悉吧。