快速傅立叶变换算法及其实现技术.pdfVIP

第23卷第3期 邢台职业技术学院学报 V01．23No．3 of 2006年6月 JournalXingtaiPolytechnicCollege Jun．2006 快速傅立叶变换算法及其实现技术 杨永昌，杨端，王凯 (桂林航天工业高等专科学校，广西桂林541004) 摘要：介绍FFT算法的基本原理并比较深入地分析其特点，详细介绍FFT算法的DSP实现，总 结了按时间抽取的傅立叶快速变换在数字信号处理器件上的实现方法。 关键词：快速傅立叶变换；数字信号处理器件；蝶形运算 中图分类号：0174．22 文献标识码：A 一、引言 傅利叶变换是一种将信号从时域变换到频域的变换形式，是声学、语音、电信和信号处理等领域中 一种重要的分析工具。离散傅利叶变换(DPT)是连续傅利叶变换在离散系统中的表示形式，由于DFT的 快速傅利叶变换(FFT)算法，它是快速计算DFT的一种高效方法，可以明显地降低运算量，大大地提高 DFT的运算速度，运算时间缩短一至两个数量级，从而使DFT在实际应用中得到了广泛的应用。 法和一次加法，而且还提供专门的FFT指令(如位倒序寻址、实现FFT算法所必需的比特反转等)，使得 算法的DSP实现。 二、FFT算法的基本原理及其特点 FFT是一种高效实现离散傅立叶变换(DFl‘)的快速算法，是数字信号处理最为重要的工具之一。 1．FFT算法的基本原理 设X(n)是一个长度为M的有限长序列，则x(n)的N点离散傅里叶变换(DFT)为： ．N。．一-I ． x(尼)=：x(力)阿0芦，k=0，l，…Ⅳ一l (1) n=O (1)式中，形Ⅳ=e吖撕川，称为旋转因子或蝶形因子。 对于N值来说，直接计算它的DFT所需的计算量是很大的。不难发现旋转因子％=e吖撕川具有对称 性形‘_Ⅳ=一矽v“川2和周期性W‘Ⅳ=形Ⅳ“～。而FFT的基本思想是利用巩的周期性和对称性，使长序 率的目的。 FFT)两大类。DITFFT算法是时域内将每一级输入序列依次按奇／偶分成2 FFT)矛U按频率抽取的FFT(DIF 个短序列进行计算，而DIFFFT算法是频域内将每一级输入序列按奇／偶分成2个短序列进行计算。两者 的区别是旋转囚子出现的位置不同，但算法是一样的。在DITFFT算法中，旋转因子蝶出现在输入端， 而DIFFFT算法中它出现在输出端。在此，主要分析DIFFFT算法。… 2．时间抽取(DIT)的FFT算法特点分析 收稿日期：2006—02—15 作者简介：杨永昌(1979_一)，四川通江人，桂林航天工业高等专科学校电子工程系，助教。 36 万方数据 邢台职业技术学院学报 2006年第3期 位置上的附属节点值，需要第m列的P和q位置上的复数节点值，如果将E+I和以+I(g)分别存放在原 来存放在以(p)和瓦(g)的同一存储寄存器内，则为完成整个计算，实际上只需要一列存储N个复数 的寄存器，这种计算方式称为“同址”，为了使计算可按以上讨论同址地进行，计算之前应将输入数据进 行倒位序存储。 ^厂 ^厂 ． 当N=2L时，每级有；个蝶形结，第一级的；个蝶形结都是同一种蝶形运算，系数为％”。第二 ‘’ Z Z ‘1r ^， 级的芸个蝶形结共有两种蝶形运算，一种系数为％o，另一种系数为巩M4，每种各有尝个蝶形结。 Z —4 依此类推、第L级的芸个蝶形结，