f=zeros(31,31);
f(10:20,10:20)=1;
F=fft2(f);
S=abs(F);
imshow(S,[])
Fc=fftshift(f);