数字信号处理第2版张小虹4数字信号处理4课件教学.pptVIP

复乘数减少。 运算结构比基2、基4都复杂。 Nlog2N = NM →(N/3)log2N = NM/3 与基2FFT算法相比 与基4FFT算法相比 (3N/8)log2N =3 NM/8 →(N/3)log2N = NM/3 基2FFT复数乘法次数为Nlog2N = NM，分裂基的复乘数 就以(N/3)log2N = NM/3计，二者相比分裂基的复乘数也 减少了2/3。即使与基4FFT的复数乘法次数 (3N/8) log2N = 3NM/8相比，复乘数亦有减少。再考虑 几个特殊因子的影响，实际复数乘法次数会更少。 研究统计表明，分裂基算法最接近理论上所需乘法次 数的最小值。 在利用DFT作频率分析时，会遇到以下几种希望频率分 辨率得到改善的情况。 4.7 线性调频z变换算法 1、求短序列DFT时，把较短序列的N点均匀分布在单 位圆上，得到的频率分辨率2?/N是很低的。采用补零 的方法增加点数，可以提高频谱密度，但N的增加意味 着计算量的增加，是低效的算法。 2、实际问题中往往会遇到包括两个或两个以上分段 的频谱，由于对各分段的兴趣不同，所以对各分段的 之间128点窄带信号频谱进行分析，按常规方法是计 算1024点的离散傅里叶变换，然后取出需要128点 频谱采样值。 率要求并不一样。例如我们需要对?= ?? /8~? /8 3、由第2章频响函数的几何作图得知，极点离单位圆 单位圆较远时，往往很难确定极点所对应的频率。 越近，谐振峰越明显。而如图4.7-1所示，当极点离 很难确定峰值频率 jIm[z] Re[z] B A |X(ej?)| 0 ? |X(ej?)| 0 ? ? 的尖峰，有利于谐振峰频率的识别。 极点的弧线进行，那么在极点对应的频率上会出现明显 如果不是在单位圆上采样，而是沿着如图所示靠近这些 的问题，人们提出了线性调频z变换算法。 总之，为了在不增加计算量的情况下，解决频率分辨率 4.7.1 CZT定义 A=1为半径的圆上。 按(4.7-1)式频谱的计算必然是从z平面的实轴开始，以 设x(n)为已知时宽为N的有限长序列，其z变换为 式中：z= esT = e?Tej?T = Ae j? A = e?T是实数，?=?T是数字频率。 Z[x(n)] x(n) z?n n=0 = X(z) = ? N?1 如果希望得到任意起始点和以螺线规律变化的z值，可设 k=0,1,2,…,M?1 式中A0、 W0为正实数。将k=0,1,2,…,M?1代入上式， 得到 A=A0ej?0 W=W0e ?j?0 zk=AW?k (4.7-2) =A0ej?0W0 ejk?0 ?k z0 =A0ej?0 ； z1=A0 W0 ej(?0+ ?0) ； ?1 z2=A0 W0 ej(?0+2?0) ； ?2 …… zM?1=A0 W0 ej[(?0+(M?1)?0)] ； ?(M?1) 式中?0是A的起始角，由A可以确定频谱采样的起点。 W0是螺线的伸展率， W01，随着k增加，螺线内旋； W01，随着k增加，螺线外旋。 ?0是z平面上相邻 zk 之间的夹角。 zM?1=A0 W0 ej[(?0+(M?1)?0)] ； ?(M?1) k=0,1,2,…,M?1 (4.7-3) 将 zk 代入z变换公式，得到CZT的定义： CZT[x(n)] = X(zk) x(n) zk n=0 = ? N?1 ?n n=0 = ? N?1 x(n) A?nW nk (M不需要等于N) 图4.7-2 Re[z] jIm[z] ?0 zM?1 (M-1)?0 螺线采样如图4.7-2所示。 z0 =A0ej?0 ； z1=A0 W0 ej(?0+ ?0) ； ?1 z2=A0 W0 ej(?0+2?0) ； ?2 …… A0 z0 A0W0 ?1 z1 A0W0 ?2 z2 ?0 ?0 每计算一个 X (zk) 值需要有N次复数乘法，一共有M个 X(zk)，所以得到直接计算CZT的计算量为： 复数乘法 mF =N·M 当N与M较大时，可以利用FFT减少运算量，提高运算 速度。 4.7.2、利用FFT的CZT算法 1、思路 利用恒等式 n=0 ? N?1 x(n) A?nW nk X(zk) = n=0 ? N?1 x(n)A?nW = [k2+n2?(k?n)2] 2 1 (4.7-4) n=0 ? N?1 x(n) A?nW W k2 2 1 n2 2 1 =W (k?n)2 2 1 ? [k2+n2?(k?n)2] 2 1 nk= g(n) = x(n) A?nW n2 2 1 令 ， n2 2 1 ? h(n) = W 则 (4.7-5)