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