二维离散傅立叶变换.doc

图的二维傅立叶变换二维傅立叶变换建立输入图，在6464的黑色图矩阵的中心建立1616的白色矩形图点阵，形成图文件。对输入图进行二维傅立叶变换，将原始图及变换图都显示于屏幕上。调整输入图中白色矩形的，再进行变换，将原始图及变换图都显示于屏幕上比较变换结果调整输入图中白色矩形的尺寸40，44），再进行变换，将原始图及变换图都显示于屏幕上比较变换结果。的中心建立1616的白色矩形图点阵，形成图文件。 对输入图进行二维傅立叶变换傅立叶原始图及变换图都显示于屏幕上比较变换结果。调整输入图中白色矩形的调整输入图中白色矩形的尺寸40，44），再进行变换，输出图像 （2）图像的二维离散傅立叶 2 3 1 2 3 图1.1图像的二维离散傅立叶原始图及变换图都显示于屏幕上 %原始图象 f=zeros(64,64);%输入64*64的黑色图像矩阵 f(25:40,25:40)=1;%建立16*16的白色矩行图像点阵 figure(1); subplot(231),imshow(f); title(原始图像)%显示原图像 F=fft2(f);%傅立叶变换 subplot(232) imshow(abs(F));title(傅里叶变换图像);%显示傅里叶变换图像 F2=fftshift(abs(F));%频谱中心化 subplot(233); imshow(abs(F2));title(中心化傅里叶频谱图);%显示中心化傅里叶频谱图 x=1:64; y=1:64; subplot(234); mesh(abs(real(F)));title(三维频谱图);%显示三维频谱图 subplot(235) mesh(x,y,F2(x,y)); title(FFT) 2、调整输入图中白色矩形的 %原始图象 f=zeros(64,64);%输入64*64的黑色图像矩阵 f(47:63,47:63)=1;%建立16*16的白色矩行图像点阵 figure(1); subplot(231),imshow(f); title(原始图像)%显示原图像 F=fft2(f);%傅立叶变换 subplot(232) imshow(abs(F));title(傅里叶变换图像);%显示傅里叶变换图像 F2=fftshift(abs(F));%频谱中心化 subplot(233); imshow(abs(F2));title(中心化傅里叶频谱图);%显示中心化傅里叶频谱图 x=1:64; y=1:64; subplot(234); mesh(abs(real(F)));title(三维频谱图);%显示三维频谱图 subplot(235) mesh(x,y,F2(x,y)); title(FFT) 整输入图中白色矩形的尺寸40，44），再进行变换的程序 40×40 clear %原始图象 f=zeros(64,64);%输入64*64的黑色图像矩阵 f(13:52,13:52)=1;%建立16*16的白色矩行图像点阵 figure(1); subplot(231),imshow(f); title(原始图像)%显示原图像 F=fft2(f);%傅立叶变换 subplot(232) imshow(abs(F));title(傅里叶变换图像);%显示傅里叶变换图像 F2=fftshift(abs(F));%频谱中心化 subplot(233); imshow(abs(F2));title(中心化傅里叶频谱图);%显示中心化傅里叶频谱图 x=1:64; y=1:64; subplot(234); mesh(abs(real(F)));title(三维频谱图);%显示三维频谱图 subplot(235) mesh(x,y,F2(x,y)); title(FFT) 4×4 clear %原始图象 f=zeros(64,64);%输入64*64的黑色图像矩阵 f(13:52,13:52)=1;%建立16*16的白色矩行图像点阵 figure(1); subplot(231),imshow(f); title(原始图像)%显示原图像 F=fft2(f);%傅立叶变换 subplot(232) imshow(abs(F));title(傅里叶变换图像);%显示傅里叶变换图像 F2=fftshift(abs(F));%频谱中心化 subplot(233); imshow(abs(F2));title(中心化傅里叶频谱图);%显示中