如果size(I，3)>1
g = im 2 double(RGB 2 gray(I))；
其他
g = im2 double(I)；
結束
%還可以在找譜之前加壹個窗口，這裏省略窗口。
f = FFT shift(ff T2(G))；
imF = log 10(ABS(F)+1)；%振幅譜取對數。
imshow(imF，[])；