FFT算法DSP实现(DIT基2FFT运算流图).pdfVIP

FFT 算法的研究与DSP 实现 1 2 王靖琰 ，王靖斌 1 中国科学院上海应用物理研究所，上海 (201800) 2 山东农业大学信息工程与技术学院，泰安 (271018) E-mail：wjycas@ 摘 要：本文介绍了DIT 基2FFT 、DIF 基2FFT 、分裂基FFT 和混合基FFT 算法的原理， 并在TMS320C6713B 浮点型DSP 上实现了算法，最后对各种算法进行了对比。 关键词：时域抽取基2FFT ，频域抽取基2FFT ，分裂基FFT ，混合基FFT ，TMS320C6713B 1．引言 傅里叶变换是一种将信号从时域变换到频域的变换形式。是声学、图像、电信和信号处 理等领域中一种重要的分析工具[1] 。其中，离散傅里叶变换(DFT)更是数字信号处理领域不 可缺少的工具之一，特别是快速傅里叶变换(FFT) 的出现使得DFT在实际应用中得到了广泛 的应用。FFT算法就是不断地把长序列的DFT逐次分解为几个短序列的DFT ，并利用的周期 性和对称性来减少DFT 的运算次数这一基本原理来加快运算速度的。 一般高基数FFT算法要比低基数FFT算法运算量要小，计算速度要快，但需要具有更强 大的数据并行处理能力的器件支持且编程难度也增大了。本文分别在浮点型 DSPTMS320C6713B上实现了DIT 、DIF基2FFT算法，分裂基FFT算法和混合基FFT算法及其 逆变换，并对语音信号进行变换和逆变换的实验，对比分析了各种算法的时间复杂度与准确 度。 2 ．各种FFT算法原理介绍与比较 现有的FFT算法包括针对N=2L 点的基2FFT 、基4FFT 、分裂基FFT 以及针对N不等于2L 的混合基FFT 。其中每种算法又可以分为时域抽取算法（DIT ）和频域抽取算法（DIF ）。下 面分别介绍它们各自的特点并对它们进行比较，导出各自的使用场合。 2.1 DIT 和 IDF 基 2 FFT 算法 L [2] 对于 N=2 点 DFT 可按照如下方法做时间抽取 ： N N N N −1 −1 −1 −1 2 2 2 2 X (k ) x (2r )W 2rk =+ x (2r +1)W (2r +1)k x (2r )W rk =+W k x (2r +1)W rk (1) ∑ N ∑ N ∑ N N ∑ N r 0 r 0 r 0 2 r 0 2 N N 这样就把一个 N 点 DFT 转换为 个基本碟形运算和 2 个 点 DFT 运算，如时域抽取 2 2 基 2FFT 一级分解框图所示： -1- x(2r) N/2 点 X(k)