物理分辨率与计算分辨率的探讨.docx

Harbin Institute of Technology 数字信号处理报告 题 目：物理分辨率与计算分辨率的探讨 院 （系） 电子与信息工程学院 学 科 信息与通信工程（51） 学 生 学 号 哈尔滨工业大学 物理分辨率与计算分辨率的探讨 摘要：本文研究了有限长序列在进行DFT变换后的频谱泄漏问题，对于产生这种现象的原理进行了方针和分析，重点介绍了物理分辨率和计算分辨率，明确了二者的定义，分析了他们的区别和联系，并从分辨率的角度解释了有限长序列频谱的一些“奇怪”的现象，最后对减小频率泄漏的方法做了简单说明。 关键词：频谱泄漏 物理分辨率 计算分辨率 分辨率是信号处理课程中一个基本的概念，它包括频率分辨率和时间分辨率，这里我们重点来研究前者。 频率分辨率可从两个方面来定义：第一，某个算法（如谱分析方法、功率谱估计方法等）将原信号中两个靠的很近的谱峰仍然能保持分开的能力，即物理分辨率；第二，在进行DFT时，频率轴上所能得到的最小频率间隔 ，即通常所说的计算分辨率。一般说的频率分辨率是指物理分辨率。 下面针对具体问题做出分析。（文中所有图的横坐标均采用归一化角频率） 1 频谱泄漏问题 DFT是有限长序列的傅里叶变换，有限长数据就相当于原来的无限长数据在时域乘了一个矩形函数（数据截断），时域相乘对应频域卷积，必然会造成得到的频谱与原来的频谱不相同，主要的表现为频谱的扩散，也就是通常所说的频谱泄漏。可以设想如果窗谱是冲击函数，频域卷积之后与原来的频谱将相同，但是冲击函数对应的时域窗为无穷宽的均匀函数，相当于没有进行数据截断，所以，只要窗函数有一定宽度，频谱泄漏是必然要发生的。 但是我们在MATLAB中运行材料中给的程序时发现了一个奇怪的现象，即非整数周期的信号采样的DFT结果会“走样”（频谱泄漏），如图1；而整数周期的，不管我们怎么调整周期个数却不会出现频谱泄漏的问题，如图2、图3和图4，这似乎是与理论分析相矛盾的，这是怎么回事呢？ 图1 非整数周期 图2 一个整数周期 图3 两个整数周期 图4 三个整数周期 2 从分辨率的角度分析 在给出的材料中，从时域和频域两个角度对以上的现象进行了分析，其中时域的分析是很好理解的，频域的分析相对要复杂一下。下面我将结合材料中的频域分析，从分辨率的角度对这个问题进行讨论。 首先，需要明确的是DTFT和DFT的关系。DTFT是离散时间序列的傅里叶变换，把序列映射到连续归一化频率域；DFT是离散时间序列的离散傅里叶变换，把序列映射成离散的频率序列。DTFT是具有物理意义的变换，DFT则是用于近似计算DTFT的工具，而FFT只是DFT的快速算法。最终我们看到的DFT结果是由两部分合成出来的，一部分是具有物理意义的信号的DTFT结果，另一部分是分析手段（加窗，补零等）所带来的误差信号。从材料中，我们已经认识到“DTFT是有限长离散序列的灵魂”。 本文开头简单介绍过物理分辨率和计算分辨率的定义，下面做具体按说明。 物理分辨率是指两个靠的很近的频谱峰值能够分辨的能力，可用来表示。一般来说，在时域抽样率一定的情况下，信号长度越长，即抽样点越大，则物理分辨率越高。有这样的关系 （1） 其中是时域抽样间隔。 需要注意的是，这个是指真正实际的信号长度，抽样点数也是指这个长度上的抽样点数，而不是补零以后的长度或抽样点数。也就是说物理分辨率只取决于时域信号的长度。 计算分辨率是指对于一个点序列做点DFT，所得到的每两根谱线间的距离 （2） 而这里的不再是实际的点数，而是计算DFT时候的点，如果经过补零的话，将是补零以后的点数。 在MATLAB程序中，物理分辨率是实际的分辨率，可是我们看到的都是DFT之后的结果，也就是计算分辨率。所以，当物理分辨率足够高的时候，我们可以适当提高计算分辨率，这样看不到的谱分量就能看到了。但是当时域信号长度不足时，物理分辨率低，即使再怎么提高计算分辨率，也是无济于事的。 这样，所有的理论依据都已经形成，可以用来解释频谱泄漏时的“奇怪现象”了。 整数周期的信号采样的DFT结果真的“不会”出现频谱泄漏的现象吗？ 其实我想原理已经很明白了，频谱泄漏是肯定存在的，只是抽样点少导致计算分辨率低，使我们没有看到，让人产生没有频谱泄漏的错觉。我们假定时域抽样信号x1长度为N（图中为32），在材料给出的程序中，有这样一条语句y1= fft(x1)，做序列的N