matlab数字图像处理 有约束最小平方恢复

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