傅里叶变换

%采用傅里叶快速变换生成频谱图
%------方法一
f=imread('lv.jpg');
f=rgb2gray(f);
r=f;
p=fft2(f);
p=fftshift(p);
%p=abs(p);
imshow(log(1+abs(p)),[]);


%-------------
figure;
j=fft2(r);
j=fftshift(j);
j=abs(j/256);
imshow(uint8(j),[]);