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