I =imread('C:\Users\wangd\Desktop\in000155.jpg'); %读入原图像文件
I1 = rgb2gray(I);
subplot(1,2,1);imshow(I1); %显示原图像
fftI1=fft2(I1); %二维离散傅立叶变换
sfftI1=fftshift(fftI1); %直流分量移到频谱中心
RR1=real(sfftI1); %取傅立叶变换的实部
II1=imag(sfftI1); %取傅立叶变换的虚部
A1=sqrt(RR1.^2+II1.^2); %计算频谱幅值
A1=(A1-min(min(A1)))/(max(max(A1))-min(min(A1)))*225;%归一化
subplot(1,2,2);imshow(A1); %显示原图像的频谱
