数据信号处理离散傅里叶变换(dft).ppt

* * 　　栅栏效应与频率分辨率是两个不同的概念： 如果在长度为N的序列后补N个零，再进行2N点DFT，使栅栏宽度减半，从而减轻了栅栏效应，得到高密度谱。 但是这种截短后补零的方法不能提高频率分辨率。因为截短已经使频谱变模糊，补零后仅使采样间隔变小，但得到的频谱采样的包络仍是已经变模糊的频谱，所以频率分辨率没有提高。因此，要提高频率分辨率，只有通过增加信号的有效持续时间Tp，来增加采样点数，才能得到高分辨率谱。 * 1 * * * 　　例如，要求计算序列在半径为r的圆上的频谱，那 么N个等间隔采样点为　　　　　，　　　　　 的频谱分量为   令　　　　　　，则 * * (3.4.19) 上式说明，要计算x(n)在半径为r的圆上的N点等间隔频谱分量，可以先对x(n)乘以r－n，再计算N点DFT(FFT)即可得到。若要求x(n)分布在该圆的有限角度内的N点等间隔频谱采样，可以在尾部补M－N个零，仍按(3.4.10)式用DFT分析整个圆上的M点等间隔频谱，最后只取所需角度内的N点等间待间隔采样即可。显然，这种方法的计算量大，效率低。下面要介绍的Chirp－Z变换可使这种谱分析的运算量大大减少。 * * 　　3． Chirp－Z变换 　 设序列x(n)长度为N，要求分析z平面上M点频谱采样值，设分析点zk为  　　　　　　　　　　　　　　　　　　　　（3.4.20） 式中A和W为复数，用极坐标形式表示为 (3.4.21) * * 式中，A0和W0为实数。当k=0时，有  (3.4.22)  由此可见，(3.4.20)式中，A决定谱分析起始点z0的位置，W0的值决定分析路径的盘旋趋势，φ0则表示两个相邻分析点之间的夹角。如果W01， 则随着k增大，分析点zk以φ0为步长向外盘旋；而W01时，向内盘旋, 如图3.4.9所示。 * * 图3.4.9 Chirp－Z变换分析频点分布 * * 将zk代入Z变换公式得到:   利用下面的关系式：   得到： * * 令  (3.4.23)  (3.4.24)  则  (3.4.25) * * (3.4.25)式说明，长度为N的序列x(n)的M点Chirp-Z变换可通过预乘得到y(n)，再计算y(n)与h(n)的卷积v(k)，最后乘以　　　这三个步骤得到。计算方框图如图3.4.10所示。图中，h(n)=　　, 看成一个数字网络的单位脉冲响应，输出v(n)=y(n)*h(n)|n=k=V(k)。 　　当W0=1时， 　　　　，是一个频率随时间线性增长的复指数序列。在雷达系统中，这样的信号称作线性调频信号，并用专用词汇Chirp表示，因此对上述变换起名为线性调频Z变换，简称Chirp-Z变换(CZT)。 * * 图3.4.10 Chirp－Z变换的计算框图 * * 　　下面介绍用DFT(FFT)计算Chirp－Z变换的原理和实现步骤。首先要确定线性卷积的区间。由于序列x(n)的长度为N(即0≤n≤N－1)，因此y(n)的长度也是N，如图3.4.11(a)所示。然而　　　　　是无限长序列，v(n)=y(n)*h(n)必然是无限长序列。而谱分析点数为M，即我们只对0≤k≤M－1区间上的卷积结果感兴趣，因此只要计算出V(k)=v(n)在区间［0 ~M－1］上的M个值就可以了。根据上述要求，只要截取区间［－(N－1)~M－1］上的h(n)就可以了。如图3.4.11(b)所示，这时经线性卷积所得v(n)长度为2N＋M－2。由(3.4.3)式知，y(n) L h(n)是v(n)的周期延拓序列的主值序列，延拓周期为L，即 * * 图3.4.11 Chirp-Z变换中hL(n)的形成 * * 而v(n)的非零区间为［－(N-1)，N＋M－2］。为了用循环卷积代替线性卷积计算出v(n)在［0 ~M－1］区间上的M个序列值，以L为周期