clear?all;
close?all;%‘Fig2.22（a）.jpg‘
a=imread（‘C:\Documents?and?Settings\Administrator\桌面\49.bmp‘）;
[nmp]=size（a）;
figuresubplot（121）imshow（a）;title（‘原圖像‘）;
a=double（a）;
b=zeros（19）;
e=zeros（nm）;
f=zeros（nm）;

for?i=2:（n-1）
????for?j=2:（m-1）
????????b（11）=-1+a（i-1j-1）;
????????b（12）=-2+a（i-1j）;
????????b（13）=-3+a（i-1j+1）;
????????b（14）=-4+a（ij-1）;
????????b（15）=-5+a（ij）;
????????b（16）=-6+a（ij+1）;
????????b（17）=-7+a（i+1j-1）;
????????b（18）=-8+a（i+1j）;
????????b（19）=-9+a（i+1j+1）;
????????min=255;
????????for?k=1:9?????????
????????????if?b（1k）????????????????min=b（1k）;
????????????end
????????end
????????e（ij）=min;
????end
end
for?i=2:（n-1）
????for?j=2:（m-1）
????????b（11）=1+e（i-1j-1）;
????????b（12）=2+e（i-1j）;
????????b（1