FFT算法仿真与性能分析.docx

有什么问题可以联系我哈,qq:1040426813都是过来人,我很乐意帮助学弟学妹滴~~~西南科技大学课 程 设 计 报 告课程名称：通信专业方向设计 设计名称：题目3：FFT算法仿真与性能分析姓 名：王帅 学 号:班 级：通信0804 指导教师：秦明伟起止日期：2011.12.12-2011.12.26西南科技大学信息工程学院制课程设计任务书学生班级：通信0804学生姓名：王帅学号计名称：题目3：FFT算法仿真与性能分析起止日期：2011.12.12-2011.12.26指导教师：秦明伟设计要求：（1）分析不同基（基2、基4）下1024点FFT算法的实现结构（蝶形单元）；（2）分析不同数据点长度下（1024、512、256、128）基2FFT算法的实现结构（蝶形单元）；（3）给出不同数据点长度下（1024、512、256、128）信号的FFT结果对比分析（谱）；课程设计学生日志时间设计内容2011.12.12分析第三题,大致构思整体思路;2011.12.13上网查阅有关知识;在图书馆,借阅相关资料;2011.12.14温习《数字信号处理》中有关FFT的知识；2011.12.15全面理解FFT的原理；2011.12.16分析基2下1024点FFT算法的实现结构;2011.12.17分析基4下1024点FFT算法的实现结构;2011.12.18分析1024点基2FFT算法的实现结构；2011.12.19分析512点基2FFT算法的实现结构；2011.12.20分析256点基2FFT算法的实现结构；2011.12.21分析128点基2FFT算法的实现结构；2011.12.22编写FFT的Matab实现程序,分析不同点的频谱;2011.12.23撰写设计报告；2011.12.24撰写设计报告；2011.12.25整理思路,温习相关知识,细心准备答辩；2011.12.26答辩课程设计考勤表周星期一星期二星期三星期四星期五课程设计评语表指导教师评语：成绩：指导教师：年月日FFT算法仿真与性能分析分析研究的目的和意义掌握FFT的基本原理；学习Matlab的相关使用，学会用Matlab解决实际问题；培养自己全面、独立思考的能力。设计原理FFT算法概述传统直接计算DFT有很多问题。设复序列x(n) 长度为N点，其DFT为：其中（ k=0，，…，N-1），整个N点DFT运算共需要实数乘法次数4 N2次，实数加法次数：N×2(2N-1)= 2N(2N-1)次。比如当N==1024点时，DFT运算时间= ==1048576次，即一百多万次的复乘运算！对于实时性很强的信号处理(如雷达信号处理)来讲，迫切需要改进DFT的计算方法，以减少运算量。如石油勘探，24道记录，每道波形记录长度5秒，若每秒抽样500点，则●每道总抽样点数=500*5=2500点；● 24道总抽样点数=24*2500=6万点；● DFT运算时间=N2=(60000)2=36*108次；显然，直接计算DFT所带来的运算量是十分庞大的，那如何去解决呢？思路是：由于DFT的运算量与成正比，如果一个大点数的DFT能分解为若干小点数DFT的组合，显然可以减少运算工作量！比如，将一个N（=2m）点的长序列的DFT分解为N/2个2点的短序列的DFT，则运算量由N2变为N/2×22，对于大N，可显著减少运算量，提高运算效率。例如10242(1024/2)×22；另外，还可利用系数的以下特性对DFT进行分解：对称性周期性可约性（4）还有两条， 。根据这两点提高计算效率的方法，提出了多种不同的算法。1965年,库利(J.W.Cooley)和图基(J.W.Tukey)在Mathematics of Computation杂志发表论文《机器计算傅立叶级数的一种算法》，首先提出FFT算法,这就是库利——图基算法（时间抽取基2FFT算法）。对于N点DFT,仅需(N/2)log2N 次复数乘法运算.例如N=1024=210 时，需要(1024/2)log2 210 =512*10=5120次。5120/1048576=4.88% ,速度提高20倍在此基础上又发展出许多改进算法如:桑德——图基算法（频率抽取基2FFT算法）；WFTA算法（小N素数组合算法）；CZT算法（线性调频Z变换算法）；ZFFT算法（局部频谱细化算法）等等。（本设计研究是以按时间抽取DIT为例）基2下FFT算法原理设序列点数 N = 2L（L 为整数，若不满足，则补零）。N=2的整数幂的FFT算法称基-2FFT算法。将序列x(n)按n的奇偶分成两组：则x(n)的DFT:其中,，式中，X1(k)和X2(k)分别是x1(r)和x2(r)的N/2的DFT。另外，式中k的取值范围是：0，