图像频谱图画图——matlab

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