%已知一個30×30大小的二值圖像，在圖像中間有個長為5高為20的白色區域，其它為黑色。
%對這幅圖進行傅立葉變換分析（主要應用FFT算法）
f = zeros(30,30);
f(5:24,13:17) = 1;
%定義圖像數組
imshow(f,'notruesize')

F = fft2(f);
%二維傅立葉變換（FFT算法）
mesh(fftshift(abs(F)));
%繪制頻譜圖

F2 = fftshift(log(abs(F)));
imshow(F2,[-1 5],'notruesize'); 
%顯示頻譜圖像，頻譜的零頻率系數被移到頻譜中間
colormap(jet); colorbar

%在上面的變換前的矩陣沒有被填充，下面比較一下填充矩陣填充后的情況。
F = fft2(f,256,256);
%在變換前f被用0填充成256×256的矩陣，變換后的矩陣大小也是256×256
imshow(fftshift(log(abs(F))),[-1 5]); 
colormap(jet); colorbar