一种实序列FFT的新算法.pdfVIP

维普资讯 第 21卷 第 2期 东北重型机械学院学报 Vo_【2lNo2 Jun 1997 — 种实序列FFT的新算法 塑 Tlvql1．7 —丽 T^／9／／l／ 摘 要 针对实际工业应用中信号多为实数的特点，综合PFA，CTA，WFTA等算法的优 分关霎键类词号餮至TN7137翠快速算法1I共轭对一法并1=。[二T{榴I】够勺L 0 引言 自从 1％5年Cooley—Tukey提出离散傅立叶变换算法…以来，人们又相继研究提 出了许多快速算法，如Splix—Radix算法2【1，温纳格雷德(wFrA)算法 3【1、素因子算法 (PFA)【 公因子算法 (CFA)等，这些算法在减少运算量和提高运算效率方面都有很多 优 ．然而，由于实际工业应用中信号多为实数序列，而以上算法都以复数形式输入， 因此存在一定冗余计算 、1985年黄俊钦和刘整扯提出实数序列FFT的一种快速算法I6】， 刺用实序列频谱的复共轭对称陛，使基2递归算法中每个子序列的谱值只需计算一半， 但其在节省乘法次数同时增加了加法次数 ．本文结合了PFA，CFA和WFTA算法的优 点，并充分利用了实数序列输入的特点，提出了一—种实序列FFT的快速算法，经比较该 算法有更少的运算量 、 l 理论基础 一 维DFT变换为二维或多维DFT，可以获得许多计算上的优越性．根据—维输入数 据个数 有CFA和PFA两类变换形式 ． 当N=N1v‘2，(Nl，Ⅳ2)=l时．可将一维DFT映 捞 二维DFT．N点输入序列 的DFT为 x()=∑ 月) ， ，Ⅳ=e—j2n／ 刺用下标变换式={Ⅳ2町 ”1+NI ) k=( l+Ⅳlk2)Ⅳ可得=维傅立 叶变换形式： l,k2)=∑l∑x(n2) I 其中， 《·)表示mad ，一=0，l，…，一1；f=l，2；ⅣlⅣ =lrood ，^r2Ⅳ三-=l mod Nt． 实现方法为：先将Ⅳ点输入数据x(nl，2)分为N1组，每组有 个点．对每组数据完 成Ⅳ2点DFT，得到 ”1，2)，送回原址；然后按 2点分组，得到v‘2组x(n1，2)，每组完 ■ 成Ⅳl点DFT再分刷送回原址，即得到骺 的 1，2)．即由Ⅳl组Ⅳ2点DET和 组Ⅳl 1995年4月19日收到 维普资讯 蔓!塑 生苎 二 型!!!堕 !!! 点DFT嵌套组合而成 ．这就是PFA算法 ． 当Lv=N1N2，(Ⅳl，Ⅳ2)≠1时，利用下标变换式 ”=”l+Nin~，七=,v2l+ 2，可得： x(kl2)=∑I∑ l，2)诒2 l[ ] 1。0 L O · 其中 n；,kT0，L…，一1；f_L，2．式中右边第二个括号内的因子称为附加旋转因子 (TF)． 它可以归入对”l的求和，也可归入对 2的求和 ．因此．仍然可以实现由二个—维变换的 嵌套 ．只是由于TF的存在．不可改变嵌套次序，这就是CFA算法 ． CFA的实现方法为：先对N个输入数据x(nl，n2)按r~1分组．完成Ⅳl组Ⅳ2点DFT, 得到X(ttDkg；再按 分组，得到 组v‘】点DFT输入值，对每一个值乘一项 毋 项，作为新的输入值，最后作 组 点DFT-得到x(kl, ．即Jv闺 【点DFT运算， Ⅳ次与} ：的相乘运算，以及Ⅳ2组Nt,