实验截图:
- 快速傅里叶变换结果图:
实验代码:
img = rgb2gray(imread('erciyuan.jpg'));
img1 = fft2(img);
img1 = fftshift(img1);
img1 = log(abs(img1)+1);
subplot(1,2,1);
imshow(img);
title('原图');
subplot(1,2,2);
imshow(img1,[]);
title('快速傅里叶变换结果');
结果分析:
- 使用matlab自带的fft2()和fftshift()函数可以实现简单的快速傅里叶变换,fft2()可以返回矩阵的二维傅里叶变换,fftshift()通过将零频分量移动到数组中心,重新排列傅里叶变换。