二维卷积定理

Question

我试图通过MATLAB验证2D问题的卷积定理：在空间域中使用滤波器进行卷积相当于在频域中与滤波器相乘。我写了以下代码。在第2步之后，我得到一个模糊的图像。但是在第3步之后，我没有得到与预期相同的模糊图像，代码的输出图如下所示。我的假设出了什么问题？

% 1- Read image and create a domain frequqncy window
im = imread('High_Detail.jpg');
im = imresize(im,[100,100]);
im = rgb2gray(im);
im = im2double(im);
window_fft_shifted = zeros([100,100]);
window_fft_shifted(50-11:50+11,50-11:50+11) = 1;

% 2- FFT then Multiplying it by window then get IFFT, a blurred image expected
im_fft = fft2(im);
im_fft_shifted = fftshift(im_fft);
im_fft_shifted_W = im_fft_shifted .* window_fft_shifted;
im_result_1 = ifft2(ifftshift(im_fft_shifted_W));

% 3- Get IFFT of window, convolve im with it, the same blurred image is expected
window = ifft2(ifftshift(window_fft_shifted));
im_result_2 = conv2(im, window);

% 4- Plot
figure;
subplot(2,2,1); imshow(im);
subplot(2,2,3); imshow(im_result_1);
subplot(2,2,4); imshow(im_result_2);

Answer 1

该行：

window = ifft2(ifftshift(window_fft_shifted));

应包括ifftshift以补偿MATLAB反对称信号表示：

window = ifftshift(ifft2(ifftshift(window_fft_shifted)));