conv2d.m

采用matlab编写的数字图像恢复程序

  function Ku = conv2d(u,k_hat)

%  Ku = conv2d(u,k_hat)
%
%  Use 2D FFT to compute convolution u*k.  k_hat= fft2(k).

   m = max(size(u));
   m2 = 2*m;
   ze = real(ifft2( ((fft2(u,m2,m2)) .* k_hat), m2,m2 ));
   Ku = ze(1:m,1:m);