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