f=zeros(30,30);
f(5:24,13:17)=1;
imshow(f, 'notruesize');
F=fft2(f,256,256); % 快速傅立葉變換算法只能處矩陣維數為2的冪次，f矩陣不
                      % 是，通過對f矩陣進行零填充來調整   
F2=fftshift(F);      % 一般在計算圖3-1的圖形函數的傅立葉變換時，坐標原點在
                      % 函數圖形的中心位置處，而計算機在對圖像執行傅立葉變換
                      % 時是以圖像的左上角為坐標原點。所以使用函數fftshift進 
                      %行修正，使變換后的直流分量位于圖形的中心；
figure,imshow(log(abs(F2)),[-1 5],'notruesize');