一种在复平面内用图形分析进行信号频谱分析的方法.docVIP

一种在复平面内用图形分析进行信号频谱分析的方法 信号处理中，经常会碰到数字频率、模拟频率这样的概念。如何理解它们之间的关系这个问题一直困扰着许多人。特别是对于数字频率，或者只有一个模糊的概念，或者只简单的知道与2有关。本文给出了一种在复平面内用图形分析的方法进行信号频谱分析的方法，详细地分析了数字频率与模拟频率之间的关系，同时很好的解释了频谱混叠现象。 谈到频率，首先映入脑海的是一串周期性的信号，比如正弦波、方波等等。而在生活中，与频率有关的事务很多，像建筑工连续的敲打声、每周晨跑的次数、走路步伐的快慢等等，生活中无处不在。这里所谈的频率是指信号处理中的信号频率，模拟信号的频率是很直观的，如国内家用电器用到的220V、50Hz交流电。这里将先从离散傅里叶变换入手，对数字频率与模拟频率的关系给出更准确地解释，并对奈奎斯特采样定理和频谱混叠现象作另一番解析，这将对从事信号处理特别是数字信号处理工作的工程师们有一定的帮助。 在复平面图中分析离散傅里叶变换[1] 对于数字频率，很多电子设计工程师只有一个模糊的概念，或者说只是简单知道它与有关系。而对数字频率和模拟频率的关系并不太了解，或者只是简单的知道他们存在下面这样的关系 (1) 式中表示模拟频率，表示数字频率，表示采样频率。 为了说明数字频率与模拟频率之间的关系，引入序列的离散傅里叶变换[2]： 　　　 　 (2) 将上式展开成矩阵形式[3]，即 为旋转因子。对取模值及功率谱[4]，上面的矩阵在复平面内分布如下图1。图中只显示方向，不代表大小。 图1 功率谱在单位圆上的分布 例如取k=1，就是图1中N个复向量叠加，k=2，就是图2中N个复向量叠加，以此类推。图中向量的模值为的大小，为方便在图中不予表示，而只是显示在复平面内的方向。 (3) (4) 数字频率与模拟频率之间的关系 有了以上分析，再来理解数字频率与模拟频率之间的关系。 当为实序列方波信号时，则的功率谱在图1的圆上均匀分布。下图4是方波信号的功率谱大小的分布图。注意图1为图2的俯视图。 图4 序列的功率谱在单位圆周上的立体分布图 根据傅里叶变换的共轭对称性[5]，与一样大，即图1 的权值是关于实轴对称的，这在图4中也可以看出来。 有了以上分析，我们再看式(1)，并将它变成如下形式： (5) 式5中K是比值。为分析二者的关系，将数字频率置于图4中来分析，并作出的位置分布图5，可得出两点结论: (1) 可在圆周上任意位置，然而由于圆的周期性，的实际处理范围限制在； (2) 由于傅里叶变换的共轭对称性，当对未知信号进行离散傅立叶变换后，得到当频率分布在实轴下方时，则它与关于实轴对称分布的频率采样效果是一样的。因此有实际意义的范围就只在范围内了。 图5 功率谱 这样一来，比值K的范围是。因此对任何模拟频率为的未知信号以采样率采样后，对我们有意义的范围是，即。到这里，我们可以从离散傅立叶变换过程中看到了奈奎斯特采样定理是怎么得来的了。 可见经过离散傅立叶变换后得到的是怎样反过来得到模拟频率的，所以在一定意义上说是合理的，而是不合适的。因为当时，就会发生频谱混淆现象，所以就不适合了。而当时还会发生频谱重叠现象。 图形法分析频谱混叠 下面举例分析频谱混叠现象。频谱混叠应该包含两层含义：第一层含义是频谱混淆；第二层含义是频谱重叠。先来看频谱混淆。例如用采样频率胃Hz的对信号对频率为Hz的信号进行采样，这时，这样一来应该是，分布于实轴下方，如图4(右)所示。根据以上分析，这与对称分布在实轴上方的发生了混淆。因此实际采样得到的频率应该是，对应的模拟频率为(80-50H)Hz =30Hz，这与一致。用时域信号的图形来表示如图6(左)所示，实线表示正弦波频率为50Hz，虚线表示正弦波频率为30Hz，线杆表示用80Hz的采样点。 图6 用80Hz采样50Hz频率信号出现的频谱混淆现象 再来看频谱重叠现象，即分析采样频率在条件下会出现的情况。这时，由知道，得到的数字频率的范围是，这在圆周上表现为周期为的周期性，即又回到了实轴以上的平面，与频率发生了重叠，频谱重叠这一词就是这样得来的。例如当为1KHz，信号频率1.2KHz，采样后的频率为(1.2-1)KHz=0.2KHz，这与是一致的，如图7所示：实线表示模拟信号频率为1.2KHz，虚线表示频谱重叠形成的频率为0.2KHz的信号。 图用1000Hz采样1200Hz频率信号出现的频谱重叠现象 继续旋转增大，它又回到了实轴以下的平面，这时既会出现频谱重叠现象又会出