shannon计算的matlab函数代码

clear  all
close
w=imread（‘woman.bmp‘）; 
% subplot（121）; 
imshow（w）; 
figure;imhist（w）;
%w=double（w）;
[mn]=size（w）;  L=256;  t=w;    N=m*n;
HH=zeros（LL）;    %改进的二维直方图，横坐标为像素f，均值为g，纵坐标为abs（f-g）
for i=1:m
    for j=1:n
        f=w（ij）;
        if（i==1 && j==1） g=（w（ij+1）+w（i+1j））/2; end
        if（i==1 && j==n） g=（w（ij-1）+w（i+1j））/2; end
        if（i==m && j==1） g=（w（i-1j）+w（ij+1））/2; end
        if（i==m && j==n） g=（w（ij-1）+w（i-1j））/2; end
        
        if（i==1 && j>1 && j        if（i==m && j>1 && j