fft算法心得(共9篇).docxVIP

fft算法心得(共9篇) 　　电子信息工程综合课程设计报告书　　DSP课程设计报告　　题目：学院：　　FFT快速算法C程序计算机与信息工程学院　　专业班级：09级电子信息工程2班小组成员：指导教师：时间：　　刘洋XX　　杨龙　　XX年03月21日　　目录　　1.设计目的....错误！未定义书签。　　设计目的.....................................................................错误！未定义书签。使用设备.....................................................................错误！未定义书签。　　2.设计任务与要求错误！未定义书签。　　3.原理与分析..4.实验步骤....5.实验程序....6.心得体会....错误！未定义书签。错误！未定义书签。错误！未定义书签。错误！未定义书签。　　1设计目的　　设计目的　　1．掌握用窗函数法设计FFT快速傅里叶的原理和方法；2．熟悉FFT快速傅里叶特性；　　3．了解各种窗函数对快速傅里叶特性的影响。　　使用设备　　PC兼容机一台，操作系统为WindowsXX(或Windows98，WindowsXP，以下默认为WindowsXX)，安装TC软件。　　2设计任务与要求　　1按FFT流程设计程序；　　2实现程序的抽样输入倒位序输出自然序3可以观察到各次谐波　　3原理与分析　　1．FFT的原理和参数生成公式　　rkk　　x(k)??x1(r)WN?WN?x2(r)WNrk?X1(k)?WNkX2(k)　　r?0　　2　　r?0　　2　　N　　?12　　N?12　　公式FFT运算公式　　FFT并不是一种新的变换，它是离散傅立叶变换的一种快速算法。由于我们在计算DFT时一次复数乘法需用四次实数乘法和二次实数加法；一次复数加法则需二次实数加法。每运算一个X需要4N次复数乘法及2N+2=2次实数加法。所以整个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：　　一般来说，输入被假定为连续的。当输入为纯粹的实数的时候，我们就可以利用左右对称的特性更好的计算DFT。　　我们称这样的RFFT优化算法是包装算法：首先2N点实数的连续输入称为“进包”。其次N点的FFT被连续被运行。最后作为结果产生的N点的合成输出是“打开”成为最初的与DFT相符合的2N点输入。　　使用这战略，我们可以划分FFT的大小，它有一半花费在包装输入O的操作和打开输出上。这样的RFFT算法和一般的FFT算法同样迅速，计算速度几乎都达到了两次DFT的连续输入。下列一部分将描述更多的在TMS320C54x上算法和运行的细节。　　2．程序流程图　　4实验步骤　　1．实验准备-启动TC。　　2．打开工程，编写程序；　　3．编译连接，并查看userscreen；　　5实验程序　　#include#include　　//***********定义结构体*************//structcmplex{　　floatzhen;floatjia;};　　//**********复乘****************//　　structcmplexls(structcmplexa1,structcmplexa2){　　structcmplexc;　　=**;=*+*;return(c);}　　voidmain()　　{inti,k,N=8,nv,nm,j=0,m,le,le1,ip;floatpai=,tmp,sum,*a;structcmplexu,w,t,x[8];　　//***********抽样函数***************//for(i=0;i.......................................................................................3　　DFT的Matlab实现......................................................