20第二十讲：CTZ及分裂基FFT.pptVIP

以上提出FFT算法，可以很快地求出全部DFT值。即求出有限长序列x(n)的z变换X(z)在单位园上N个等间隔抽样点z 处的抽样值。它要求N为高度复合数。即N可以分解成一些因子的乘积。例N=2L 实际上：（1）也许对其它围线上z变换取样发生兴趣。如语音处理中，常常需要知道某一围线z变换的极点处的复频率。 （2）只需要计算单位圆上某一段的频谱。如窄带信号，希望在窄带频率内频率抽样能够非常密集，提高分辨率，带外则不考虑 （3）若N是大素数时，不能加以分解，又如何有效计算这种序列DFT。例N=311，若用基2则须补N=28=512点，要补201个零点。 由上面三个问题提出： 为了提高DFT的灵活性，须用新的方法。 线性调频z变换(CZT)就是适用这种更为一般情况下，由x(n)求X(zk)的快速变换 CZT：来自于雷达专业的专用词汇。 二、算法原理 Z 变 换 采 用 螺 线 抽 样, 可 适 用 于 更 一 般 情 况 下 由 x(n) 求X(zk) 的 快 速 算 法, 这 种 变 换 称 为 线 性 调 频 Z 变 换 ( 简 称 CZT ). 为适应z可以沿平面内更一般的路径取值，故： 沿z平面上的一段螺线作等分角的抽样，则z的取样点Z 可表示为： 3.用CZT求解DFT的流图 4.CZT变换各点的值 5、图形 6、说明 （1）A为起始样点位置 三、CZT的实现步骤 8、CZT变换运算流程图 四、CZT运算量的估算 9 、CZT运算量与直接运算量比较 五、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。 六、CZT变换的应用 （1）利用CZT变换计算DFT。 （2）对信号的频谱进行细化分析。其中对窄带信号频谱或对部分感兴趣的频谱进行细化分析。 （3）求解z变换X(z)的零、极点。 用于语音信号处理过程中。 具体方法：利用不同半径同心圆，进行等间隔的采样。 作业 * 第二十讲 线性调频Z变换 六、CZT变换的应用 五、CZT算法的特点 四、CZT运算量的估算 二、算法原理 一、CTZ的引入 三、CZT的实现步骤 一、CTZ的引入 k 1.定义 已 知 x(n) ，0≤n≤N-1,则它的z变换是： 其中M：表示欲分析的复频谱的点数。M不一定等于N。A， 都为任意复数。 k 2.CZT公式推导 （2）zk是z平面一段螺线上的等分角上某一采样点。 1 2 3 4 4 3 2 3 2 1 2 3 2 3 4 6 11 20 30 16 7 26 32 11 20 30 11 20 30 当M、N足够小时，直接算法运算量少。 但M、N值比较大时（大于50），CZT算法比直接算法的运算量少得多。 例M=50，N=50，N*M=2500次 而CZT1600次。 这样CZT只对感兴趣的频率区段进行采样。计算量小很多，有利于实时处理。或在保证实时处理的情况下，可大提高频率分辨率。