MATLAB 傅里叶变换

function [ output_args ] = Untitled2( input_args )
I1=imread('lena.bmp'); 
%读入原图像文件
I2=imread('Peppers.bmp');  
%读入原图像文件
subplot(3,2,1);
imshow(I1); 
title('lena原图')
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(3,2,2);
imshow(A1);          
title('lena傅里叶变换')
subplot(3,2,3);
imshow(I2);          
title('Peppers原图')
fftI2=fft2(I2);                       
%二维离散傅立叶变换
sfftI2=fftshift(fftI2);             
%直流分量移到频谱中心
RR2=real(sfftI2);                    
%取傅立叶变换的实部
II2=imag(sfftI2);                  
%取傅立叶变换的虚部
A2=sqrt(RR2.^2+II2.^2);             
%计算频谱幅值
A2=(A2-min(min(A2)))/(max(max(A2))-min(min(A2)))*225;
%归一化
subplot(3,2,4);
imshow(A2);          
title('Peppers傅里叶变换')

end