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
MATLAB 傅里叶变换
最新推荐文章于 2023-10-07 20:46:48 发布