基--2按频率抽取的FFT算法Decimation-in-Frequency(DIF)讲稿.ppt

作业 P200页，第8、9题 第七节 分裂基FFT算法 一、发展史 自从基2快速算法出现以来，人们仍在不断寻求更快的算法。基4FFT算法就比最初的基2FFT算法更快。 从理论上讲，用较大的基数还可以进一步减少运算次数，但要以程序（或硬件）变得更为复杂为代价。甚至得不偿失。 1984年，法国的杜梅尔(P.Dohamel)和霍尔曼(H.Hollmann)将基2分解和基4分解糅合在一起，提出了分裂基FFT算法。其运算量比前几种算法都有所减少，运算流图却与基2FFT很接近，运算程序也很短。它是目前一种实用的高效新算法。 二、分裂基FFT算法原理 结论1 N/2点 DFT N/4 点 DFT . . . N/4 点 DFT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . -j 结论2 结论3 3.直接利用FFT流图方法的注意点 （1）FFT与IFFT连接应用时，注意输入输出序列的排列顺序，即应注意是自然顺序还是倒序。 （2）FFT和IFFT共用同一个程序时，也应注意利用FFT算法输入输出的排列顺序，即应注意自然顺序还是倒位序 （3）当把频率抽取FFT流图用于IDFT时，应改称时间抽取IFFT流图。 （4）当把时间抽取FFT流图用于IDFT时，应改称频率抽取IFFT流图。 作业 用C语言完成N=128点的IDIT，IDIF。 第六节线性调频Z变换 一、引入 以上提出FFT算法，可以很快地求出全部DFT值。即求出有限长序列x(n)的z变换X(z)在单位园上N个等间隔抽样点zk处的抽样值。它要求N为高度复合数。即N可以分解成一些因子的乘积。例N=2L 实际上：（1）也许对其它围线上z变换取样发生兴趣。如语音处理中，常常需要知道某一围线z变换的极点所处的复频率。 （2）只需要计算单位圆上某一段的频谱。如窄带信号，希望在窄带频率内频率抽样能够非常密集，提高分辨率，带外则不考虑。 （3）若N是大素数时，不能加以分解，又如何有效计算这种序列DFT。例N=311，若用基2则须补N=28=512点，要补211个零点。 二、问题提出 由上面三个问题提出： 为了提高DFT的灵活性，须用新的方法。 线性调频z变换(CZT)就是适用这种更为一般情况下，由x(n)求X(zk)的快速变换 CZT：来自于雷达专业的专用词汇。 三、算法原理1.定义 Z 变 换 采 用 螺 线 抽 样, 可 适 用 于 更 一 般 情 况 下 由 x(n) 求X(zk) 的 快 速 算 法, 这 种 变 换 称 为 线 性 调 频 Z 变 换 ( 简 称 CZT ). 2.CZT公式推导1 为适应z可以沿平面内更一般的路径取值，故： 沿z平面上的一段螺线作等分角的抽样，则z的取样点Zk可表示为： 已 知 x(n) ，0≤n≤N-1,则它的z变换是： 其中M：表示欲分析的复频谱的点数。M不一定等于N。A， 都为任意复数。 2.CZT公式推导2 2.CZT公式推导3 3.用CZT求解DFT的流图 4.CZT变换各点的值 5、图形 6、说明1 （1）A为起始样点位置 6、说明2 （2）zk是z平面一段螺线上的等分角上某一采样点。 6、说明3 6、说明4 6、说明5 7、CZT的实现步骤1 7、CZT的实现步骤2 7、CZT的实现步骤3 7、CZT的实现步骤4 7、CZT的实现步骤5 7、CZT的实现步骤6 7、CZT的实现步骤7 8、CZT变换运算流程图 9、CZT运算量的估算1 9、CZT运算量的估算2 10、CZT运算量与直接运算量比较 当M、N足够小时，直接算法运算量少。 但M、N值比较大时（大于50），CZT算法比直接算法的运算量少得多。 例M=50，N=50，N*M=2500次 而CZT1600次。 11、CZT算法的特点 与标准FFT算法相比，CZT算法有以下特点： （1）输入序列长N及输出序列长M不需要相等。 （2）N及M不必是高度合成数，二者均可为素数。 （3）Zk的角间隔 是任意的，其频率分辨率也是任意的。 （4）周线不必是z平面上的圆，在语音分析中螺旋周线具有某些优点。 （5）起始点z0可任意选定，因此可以从任意频率上开始对输入数据进行窄带高分辨率的分析。 （6）若A=1,M=N,可用CZT来计算DFT，即使N为素数时，也可以。 总之，CZT算法具有很大的灵活性，在某种意义上说，它是一个一般化的DFT。 12、CZT变换的应用1 （1）利用CZT变换计算DFT。 12、CZT变换的应用2 （2）对信号的频谱进行细化分析。其中对窄带信号频谱或对部分感兴趣的频谱进行细化分析。 这样CZT只对感兴趣的频