A=fft2(double(I));
figure,imshow(A);
%取傅里叶变换的模的平方;
sfftI3=fftshift(A); %直流分量移到频谱中心
RR3=real(sfftI3); %取傅立叶变换的实部
II3=imag(sfftI3); %取傅立叶变换的虚部
A3=sqrt(RR3.^2+II3.^2); %计算频谱幅值
A4=(RR3.^2+II3.^2); %取傅里叶变换的模的平方;