clear
N=100
f=zeros(50,50);
f(15:35,23:28)=1;
figure(1);
imshow(f,'notruesize');
F=fft2(f,N,N);
F2=fftshift(abs(F));
figure(2);
x=1:N;y=1:N;
mesh(x,y,F2(x,y));
colormap(gray);colorbar;
N=100
f=zeros(50,50);
f(15:35,23:28)=1;
figure(1);
imshow(f,'notruesize');
F=fft2(f,N,N);
F2=fftshift(abs(F));
figure(2);
x=1:N;y=1:N;
mesh(x,y,F2(x,y));
colormap(gray);colorbar;