图像二维FFT变换matlab仿真,对比加入窗函数之后的图像频谱和相位，含仿真操作录像

clc; clear; close all; warning off; addpath(genpath(pwd)); %-------------------------------------------------------------------------% %实验要求：用函数fft2对图像进行二维离散Fourier变换。用图像的形式分别显示清晰的 %频谱幅度和相位。对图像用窗函数加权，观察其Fourier变换，解释加窗前后的变化。 %对图像作Fourier变换和反变换。将图像频谱的幅度和相位部分分离。分别对幅度部分和 %相位部分求Fourier反变换，观察得到的空域结果并解释所观察到的现象。 %-------------------------------------------------------------------------% f=imread('test.jpg'); imshow(f); title('原始图像'); f=double(f); F=fft2(f); fd=abs(fftshift(F)); fdf=ifft2(fd);%对幅度傅里叶反变换 xw=angle(F); xwf=ifft2(exp(j*xw)); figure; subplot(2,2,1); imshow(log(abs(fftshift(F))),[])%显示频谱幅度 title('图像频谱幅度'); subplot(2,2,2); imshow(angle(fftshift(F)),[]) title('图像相位'); %--------------------------- subplot(2,2,3); imshow(log(1+abs(fdf)),[]); title('图像频谱幅度的逆变换'); subplot(2,2,4); imshow(xwf,[]); title('图像相位的逆变换'); %----------------------------------- [M,N]=size(f); H1=hamming(M); H2=hamming(N); HH=H1*H2'; figure, subplot(2,2,1); imshow(HH,[]); title('生成一个圆形的汉明窗'); g=f.*HH; subplot(2,2,2); imshow(uint8(g)); title('加汉明窗后的图像'); G=fft2(g); subplot(2,2,3); imshow(log(abs(fftshift(G))),[])%显示加窗后的频谱幅度 title('加窗后的频谱幅度'); subplot(2,2,4); imshow(angle(fftshift(G)),[])%显示频谱幅度 title('加窗后的图像相位'); %----------------------------------------------------