


H = fftshift(hpfilter('gaussian', 500, 500,10));
surf(H(1:10:500, 1:10:500));
axis([0 50 0 50 0 1]);
colormap(gray);
grid off; axis off;