第3章5-8 离散傅里叶变换及其快速算法.ppt

第三章 离散傅里叶变换及其快速算法 3. 6 N为合数的FFT算法 例3.2 设xa(t)＝e-t(t≥0)，要求对该信号进行FFT分析。已知频率分辨率为F≤5Hz，信号的最高频率为f0=40Hz。试确定最小记录长度tp、取样周期T和记录长度中的最少取样点数N，并用FFT求出信号的频谱X(k)、振幅谱IX(k)|和相位谱θ(k)。 现将x(n)分解为 用FFT计算循环卷积 令 因此 3.8 线性调频Z变换 在实际应用中，有时只对信号的某一频段感兴趣，或只需计算单位圆上某一段的频谱值。例如，在对窄带信号进行分析时，常希望在窄频带内对频率的取样很密集，以便提高频率分辨率，而在窄频带外不予以考虑。 对于这种情况，如果采用DFT方法，则需要在窄频带内外都增加频域取样点，而窄频带外的计算量是浪费的。此外，有时对非单位圆上的取样感兴趣，例如在语音信号处理中，常常需要知道其Z变换的极点所在处的复频率，这时就需要在这些极点附近的曲线上进行频域取样，这样，就要沿着螺旋线对Z变换取样。 这种沿螺旋线上取样点计算的Z变换，称为线性调频Z变换(Chirp Z Transform，简称CZT)。 下面就来讨论这种变换的原理及其计算方法。 一个长度为N的序列x(n)，其Z变换为 (3.110) 为了使z可以沿着z平面更一般的路径(不只是单位圆)取值，可以沿一段螺旋线对z作等分角取样，这些取样点上的zk表示为 (3.111) 其中M为所要分析的复频谱的点数，不一定等于N。W和A为任意复数，可表示为 (3.112) 得到 (3.114) 取样点zk所在的路径如图3.38所示，图中： (1)A0表示起始取样点z0的矢量长度，通常A0≤1，否则将处于单位圆|z|=1之外。 (2)θ0表示起始取样点z0的矢量的相角，它可以是正值或负值。 (3)φ0表示两相邻取样点矢量之间的角度差。φ0为正时，表示zk的路径沿逆时针方向旋转；φ0为负时，zk的路径沿顺时针方向旋转。 (4)W0表示螺旋线的伸展率。W01时，随着k的增加螺旋线向内盘旋；W01时，则随k的增加螺旋线向外盘旋；W0=1对应于半径为A0的一段弧线，在A0=1时这段弧线是单位圆的一部分。 在 的特殊情况下， 各zk值等间隔均匀分布在单位圆上，这对应于计算序列的DFT。将式(3.111)代入式(3.110)得 (3.115) 直接计算式(3.115)，总共要算M个取样点，需要NM次复数乘法和(N-1)M次复数加法，这与DFT的直接计算类似。当N和M很大时，这个计算量是很大的，因而限制了运算速度。如果将式(3.115)中的因子Wnk的幂nk作如下的变换，则可以将以上运算转换成卷积的形式，并可采用FFT来计算，从而大大提高运算速度。具体来说，将nk写成以下形式： (3.116) 将式(3.116)代入式(3.115)得 令 和 则 (3.119) 从式(3.119)可看出，可以通过计算g(k)与h(k)(此处用变量k代替n)的线性卷积，然后再乘上 来得到X(zk)，即 (3.120) 式(3.120)的计算可用图3.39来表示。 序列g(n)可被看成是一个具有二次相位的复指数序列，这种信号在雷达系统中称为chirp信号(线性调频信号)，线性调频Z变换的名称即来源于此。 归纳起来，CZT的计算步骤如下： (1)选择满足条件L≥(N+M-1)和L=2M的整数L值。 (2)对 添补零取样值，构成长度为L点的序列g(n)。 (3.121) 用FFT计算g(n)的L点DFT，得 (3.122) (3) 按下式构造长为L的序列h(n)， (3.123) 用FFT计算式(3.123)的序列h(n)的L点DFT (4)计算Q(r)＝H(r)·G(r)。 (5)求Q(r)的L点IFFT，得到h(n)与g(n)的循环卷积，即 q(n)的前M个值等于h(n)与g(n)的线性卷积结果，对应于n≥M的值没有意义，不必计算。 (6)最后计算X(zk) 可以看出，CZT算法是非常灵活的。它的输入序列长度N和输出序列长度M可以不相等，且均可为任意数，包括素数；各zk点间的角度间隔φ0可以是任意的，因而频率分辨率可以调整；计算Z变换的取样点的轨迹可以不是圆而是螺旋线；起始点z0可任意选定，也就是可以从任意频率开始对输入数据进行分析。 (a)h(n)是一个长度为N+M-1的有限长序列。 (b)将h(n)延拓成周期为L的周期序列，并取出主值序列。 (c)对g(n)补零值至L点。 (d)q(n)的前M个值，即为h(n)与g(n)的线性卷积。 对频率再进行偶奇分，则得频率的偶数项为 频率的奇数项为 通过类似的推导可得 上面4式所表示的计算都是N/4点的DFT计