组-实验四-王栩灿-肖凌峰-张皓.doc

实验四 傅里叶变换及频率滤波 时间：2017.11.3 组员：王栩灿/肖凌峰/张皓 第18组 一、实验目的 1.了解傅里叶变换在图像处理中的应用； 2.掌握频率滤波的原理和特点； 3.掌握傅里叶变换的基本性质； 4.通过实验了解二维频谱的分布特点。 二、实验内容 1. 构造同时包含50Hz和120Hz的正弦信号进行傅里叶变换并分析。 2. 以高斯光束为例（自己选的例子亦可），利用matlab实现二维图像的傅里叶变换，选取不同的抽样间隔，观察抽样间隔对结果的影响并分析。 3. 对某一图像(如lena.jpg） （1）编程实现图像的理想低通、高斯低通、布特沃斯低通滤波，分析不同的截止频率对滤波后图像的影响。 （2）编程实现图像的理想高通、高斯高通、布特沃斯高通滤波，分析不同的截止频率对滤波后图像的影响。 （3）在空域中实现高斯低通滤波，并与频域的结果进行比较。 4. 利用FFT实现二维高斯光束的菲涅尔衍射和夫琅禾费衍射，并分析光束在不同传输距离上的衍射图样（传输距离自定，至少三个数值说明即可）。 利用一幅图像进行菲涅尔衍射和夫琅禾费衍射后会是什么样子？ 三、实验原理 ①快速傅里叶变换 有限长序列可以通过离散傅里叶变换(DFT）将其频域也离散化成有限长序列。但其计算量太大，很难实时地处理问题，因此引出了快速傅里叶变换(FFT). 1965年，Cooley和Tukey提出了计算离散傅里叶变换（DFT）的快速算法，将DFT的运算量减少了几个数量级。从此，对快速傅里叶变换（FFT）算法的研究便不断深入，数字信号处理这门新兴学科也随FFT的出现和发展而迅速发展。根据对序列分解与选取方法的不同而产生了FFT的多种算法，基本算法是基2DIT和基2DIF。FFT在离散傅里叶反变换、线性卷积和线性相关等方面也有重要应用。 快速傅立叶变换（FFT），是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。FFT的基本思想是把原始的N点序列，依次分解成一系列的短序列。充分利用DFT计算式中指数因子 所具有的对称性质和周期性质，进而求出这些短序列相应的DFT并进行适当组合，达到删除重复计算，减少乘法运算和简化结构的目的。此后，在这思想基础上又开发了高基和分裂基等快速算法，随着数字技术的高速发展，1976年出现建立在数论和多项式理论基础上的维诺格勒傅里叶变换算法(WFTA）和素因子傅里叶变换算法。它们的共同特点是，当N是素数时，可以将DFT算转化为求循环卷积，从而更进一步减少乘法次数，提高速度。 ②理想低通/高通滤波器 是一种假想的低通/高通滤波器，其对于高于/低于截止频率的信号完全截止，而对于低于/高于截止频率的信号完全无失真传输。 ③布特沃斯低通/高通滤波器 截止频率位于距原点D0的n阶布特沃斯滤波器（BLPF）的传递函数定义为 与IFPF 不同，BLPF并没有字通过频率和滤除频率之间给出明显截止的急剧不连续性。且空间域的一阶布特沃斯滤波器没有振铃现象。 ④高斯低通/高通滤波器 高斯滤波器的二维形式由下式给出： ⑤Fresnel衍射与Fraunhofer衍射 在光学里，菲涅耳衍射（Fresnel diffraction）指的是光波在近场区域的衍射，即光源或衍射的图样的屏与衍射孔(障碍物)的距离是有限的。菲涅耳衍射积分式可以用来计算光波在近场区域的传播，因法国物理学者奥古斯丁·菲涅耳而命名，是基尔霍夫衍射公式的近似。 夫琅禾费衍射是指把单色点光源放在透镜的焦点上，经过透镜后的单色平行光垂直照射衍射屏时，在屏后面不同距离上会观察到一些衍射现象，其中当屏远离到足够大的距离后，光斑中心出现一个较大的亮斑，外围是一些较弱的明暗相间的同心圆环，此后再往外移动，衍射花样出现稳定分布，中心处总是亮的，只是半径不断扩大而已，这种衍射称为夫琅禾费衍射，又称远场衍射（而当距离较近时为菲涅尔衍射，较近时没有夫琅禾费衍射）。 实验程序与结果 (构造同时包含50Hz和120Hz的正弦信号进行傅里叶变换并分析。 主程序 f1 = fun_GenSinWave(1000,1000,50); %plot(f1); f2 = fun_GenSinWave(1000,1000,120); %figure;plot(f2); f = f1+f2; y = fft(f); figure;plot(y); 正弦波生成函数 function [output] = fun_GenSinWave(fs,N,f) %part one ,generate 50Hz and 100Hz sin signal and do Fourier change fsm = fs;%采样频率 Ns = N;