[M, N] = size(I);
F = fft2(I);
sigma = 10;
H = lpfilter('gauss', M, N, sigma);
G = H.*F;
g = real(ifft2(G));
imshow(g, []);