图像频谱图画图——matlab

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);          %显示原图像的频谱