第3章 离散傅叶变换（DFT）.pptVIP

　　最后要说明的是，栅栏效应与频率分辨率是不同的两个概念。 如果截取长度为N的一段数据序列，则可以在其后面补N个零，再进行2N点DFT，使栅栏宽度减半，从而减轻了栅栏效应。 但是这种截短后补零的方法不能提高频率分辨率。因为截短已经使频谱变模糊，补零后仅使频域采样间隔变小，但得到的频谱采样的包络仍是原来已经变模糊的频谱，所以频率分辨率没有提高。因此，要提高频率分辨率，就必须对原始信号截取的长度加长（对模拟信号，就是增加采样时间Tp的长度）。 为使频率分辨率提高1倍，即F=5 Hz，要求：     　上面分析了为提高谱分辨率，又保持谱分析范围不变，必须增长记录时间Tp，增加采样点数。应当注意，这种提高谱分辨率的条件是必须满足时域采样定理，即绝对不能保持N不变，通过增大T来增加记录的时间TP。 用快速算法FFT计算时，选用N=1024点。 　　2． 用DFT对序列进行谱分析　 　 我们知道单位圆上的Z变换就是序列的傅里叶变换，即   X(ejω)是ω的连续周期函数。如果对序列x(n)进行N点DFT得到X(k)，则X(k)是在区间［0，2π］上对 X(ejω)的N点等间隔采样，频谱分辨率就是采样间隔2π/N。因此序列的傅里叶变换可利用DFT(即FFT)来计算。 　　对周期为N的周期序列　　，由(2.3.10)式知道，其频谱函数为   其中  由于　　以N为周期，因而X(ejω)也是以2π为周期的离散谱， 每个周期有N条谱线，第k条谱线位于ω=(2π/N)k处，代表 　　的k次谐波分量。而且，谱线的相对大小与 成正比。 由此可见，周期序列的频谱结构可用其离散傅里叶级数系数 　　表示。 　　由DFT的隐含周期性知道，截取　　的主值序列 　　　　　　　，并进行N点DFT，得到：  　　　　　　　　　　　　　　　　　　　　 所以可用X(k)表示　　的频谱结构。 (3.4.16) 　　如果截取长度M等于　　的整数个周期，即M=mN，m为正整数，即 　　　　　　　　　　　　　　   令n=n′＋iN; i=0, 1, …, m－1; n′=0, 1, …, N－1，则 （3.4.17） 因为 所以 　　　　　　　　　　　  由此可见，XM(k)也能表示　　的频谱结构，只是在k=im时，　　　　　　，表示　　的i次谐波谱线，其幅度扩大m倍。而其他k值时，XM(k)=0。当然，X(i)与XM(im)对应点频率是相等的　　　 所以，只要截取　　的整数个周期进行DFT，就可得到它的频谱结构，达到谱分析的目的。 （3.4.18） 　　如果　　的周期预先不知道，可先截取M点进行DFT，即  再将截取长度扩大1倍，截取  比较XM(k)和X2M(k)，如果二者的主谱差别满足分析误差要求，则以XM(k)或X2M(k)近似表示　　的频谱，否则，继续将截取长度加倍，直至前后两次分析所得主谱频率差别满足误差要求。设最后截取长度为iM，则XiM(k0)表示ω=［2π/(iM)］k0点的谱线强度。 　　在很多实际应用中，并非整个单位圆上的频谱都很有意义: 对于窄带信号，往往只希望对信号所在的一段频带进 行谱分析，这时便希望采样能密集地在这段频带内进 行，而带外部分可完全不予考虑。 另外，有时希望采样点不局限于单位圆上。 例如，语音信号处理中，常常需要知道系统极点所对 应的频率，如果极点位置离单位圆较远，则其单位圆 上的频谱就很平滑，如图3.4.8(a)所示，这时很难从中 识别出极点对应的频率。 如果使采样点轨迹沿一条接近这些极点的弧线或圆周 进行，则采样结果将会在极点对应的频率上出现明显 的尖峰，如图3.4.8(b)所示。这样就能准确地测定出极 点频率。 对均匀分布在以原点为圆心的任何圆上的N点 　　　　　　　　　　　　　　 频率采样，可用 DFT(FFT)计算；而沿螺旋弧线采样，则要用线 性调频Z变换(Chirp－Z变换，简称CZT)计算。  ) 1 1 0 e ( π 2 j - = = N k r z k N k ， ， ， ， L 图3.4.8 单位圆与非单位圆谱分析示意图 例如，要求计算序列在半径为r的圆上的频谱，那么N个等间隔采样点为　　　　，　　　　　 ，Zk的频谱分量为  令　　　　　，则 (3.4.19) 上式说明，要计算x(n)在半径为r的圆上的N点等间隔频谱分 量，可以先对x(n)乘以r－n，再计算N点DFT(FFT)即可得到。 若要求