中 南 大 学 DSP 技术实验报告 实验名称:快速傅立叶变换(FFT)算法实验 专业班级: 信息0602 学生姓名: 张倩曦 (学号:24) 指导老师: ** 完成日期: 2009年12月2日 中 南 大 学 ·信 息 科 学 与工程学 院 快速傅立叶变换(FFT)算法实验 一.实验目的 1.掌握用窗函数法设计FFT 快速傅里叶的原理和方法; 2.熟悉FFT 快速傅里叶特性; 3.了解各种窗函数对快速傅里叶特性的影响
二.实验设备 PC 兼容机一台,操作系统为Window s2000(或Window s98,Window sXP,以下默认为Window s2000) ,安装Code Composer Studio 2
三.实验原理 1.FFT 的原理和参数生成公式: 公式(1)FFT 运算公式 FFT 并不是一种新的变换,它是离散傅立叶变换(DFT)的一种快速算法
由于我们在计算DFT 时一次复数乘法需用四次实数乘法和二次实数加法;一次复数加法则需二次实数加法
每运算一个X(k)需要4N 次复数乘法及2N+2(N-1)=2(2N-1)次实数加法
所以整个DFT运算总共需要4N^2 次实数乘法和N*2(2N-1)=2N(2N-1)次实数加法
如此一来,计算时乘法次数和加法次数都是和N^2 成正比的,当N 很大时,运算量是可观的,因而需要改进对DFT 的算法减少运算速度
根据傅立叶变换的对称性和周期性,我们可以将DFT 运算中有些项合并
我们先设序列长度为N=2^L,L 为整数
将N=2^L 的序列x(n)(n=0,1,… … ,N-1),按N的奇偶分成两组,也就是说我们将一个N 点的DFT 分解成两个N/2 点的DFT,他们又重新组合成一个如下式所表达的N 点DFT: 一般来说,输入被假定为连续的
当输入为纯粹的实数的时候,我们就可以利用左右对称的特性更好