clear
I=imread（‘dip3.bmp‘）;
I=rgb2gray（I）;
subplot（121）;
imshow（I）;
k=0.0025;
[mn]=size（I）;
%spectrum=zeros（mn）;
%H=zeros（mn）;
for?u=1:m
??????for?v=1:n
???????H（uv）=exp（（-1）*k*（（u^2+v^2）^（double（5/6））））
??????end
end

F=fft2（I）;
g=imnoise（uint8（ifft