信号检测研究报告.docx

信号检测与估计研究报告——基于DFT相位的正弦波频率和初相的高精度估计方法声明：本文摘自电子学报2001年9月第9期，作者是齐国清，贾欣乐摘要：利用分段DFT频谱的相位差消除了初相对频率估计的影响且避免了相位模糊问题。给出了频率和初相估计的均方根误差计算公式。理论分析和Monte Carlo模拟结果显示频率估计均方根误差接近Cramer-Rao(CR)下限，初相估计均方根误差略高于CR下限的2倍。阈值信噪比远远低于基于时域瞬时相位的频率和初相估计方法。在信噪比为6dB，采样点数为1024的情况下，频率估计均方根误差约为DFT频率分辨率的1%，初相估计均方根误差约为2度。该方法已用于FMCW液位测量雷达并取得1mm距离测量精度。关键词：频率估计；相位估计；雷达测距；DFT1 引言 正弦波信号频率的估计是通信、雷达、声纳以及电子对抗等领域信号处理中的一个重要问题。基于参数模型的谱估计、最大熵谱估计等方法具有频率分辨率高的优点，但对于长序列，其运算量大，不利于实时处理。而基于DFT的谱分析方法，可采用快速算法，即FFT，因此运算速度快，特别适合于实时信号处理。但DFT的频率分辨力和频率估计精度取决于信号的测量时间长度，信号测量时间过长不但给实时处理带来困难，而且在一些应用中信号的持续时间是有限的，不能任意延长，使得DFT的频率分辨力和估计精度受到限制，因此一般只利用DFT实现频率的粗测[1,2]，文献[3]提出了利用数值搜索方法提高DFT频率估计精度。由于栅栏效应DFT频谱在主瓣之内有2条谱线，利用这2条谱线的幅度可以实现频率插值以及提高频率测量精度[4]。文献[5]提出了直接在时域采用最小二乘线性回归的方法，利用瞬时相位估计信号频率和初相。为了避免直接测量瞬时相位的整周模糊问题，文献[6]提出了相位查分方法。在时域测量瞬时相位估计信号频率和初相的缺点是需要较高的信噪比。文献[7]提出了一种利用DFT频谱的相位和频率插值的综合算法提高频率测量精度的方法。本文提出一种利用信号DFT频谱相位提高频率测量精度的新方法，利用分段DFT消除了相位测量中的整周模糊问题。理论分析和计算机模拟结果表明，该方法的信噪比阈值比基于时域瞬时相位的频率估计方法低很多，在很低信噪比情况下，仍能得到很高的频率和初相估计精度，频率估计均方根误差接近CR下限。该方法还可以估计信号的初相，初相估计均方根误差高于CR下限的2倍。2 FFT相位差法估计频率和初相原理设观测信号为单一频率复正弦波信号(1)其中，、、分别为信号的幅度、频率和初相。对上述信号进行采样，设信号的记录时间长度分为两个长度相同的序列，对应前点，则采样序列可记为 (2)对应后点，则(3)(4)分别对和进行点DFT，得到离散频谱(5) (6)其中和分别为的幅度项与相位项(7)(8)根据式(6)知，与的幅度项完全一样。由式(7),幅度最大处值对应的离散频率为([x]表示取最接近x的整数)。利用DFT的最大谱线粗测频率为，为DFT的频率分辨率。从式(8)可见，DFT最大谱线的相位包含信号频率与DFT最大谱线位置的偏差信息，但由于未知，不能直接利用DFT的相位来估计频率。用和分别表示和在最大谱线处的相位，则两者的差值为(9)当在范围内变化时，由式(9)知，在到之间变化。因此可以利用可以对与DFT最大谱线对应的频率的偏差进行估计 (10)也可以定义与T无关的相对频率偏差，当在范围内变化时，在范围内变化，可以利用直接得到的估计值。然后按下式计算的估值(11)由于首先利用DFT的最大谱线位置对频率进行了粗测，DFT相位差只用来估计信号频率在两条谱线之间的位置，因而相位测量不会出现模糊。将式(9)代入式(8),得 (12)于是初相估计为 (13)当N较大时，上式可近似为 (14)3 噪声对测量精度的影响 在加性白噪声背景下，观测信号可表示为。其中为式(1)所表示的纯信号，为复白噪声，其均值为0，功率谱密度为。前N/2点采样序列记为。设系统的等效带宽为，则的功率(方差)为。采样后的信噪比为。白噪声为平稳随机过程，不满足Fourier变换得绝对可积条件，不能对其进行Fourier变换，因此一般只分析其功率谱密度。但功率谱密度不包含相位信息，无法分析噪声对相位测量的影响。对于采样后的白噪声序列，可将其DFT变换看作是若干个随机变量的线性组合，每项DFT系数仍为随机变量，所以噪声序列的DFT仍为随机序列。因此我们可以定义噪声序列在概率意义上的N/2点DFT变换(15)上式只在概率意义上成立，因为对于的不同次实现，的值是随机变化的。和分别表示的幅度和相位，均为随机量。可以通过上式分析的统计特性。当为高斯白噪声序列时，也服从高斯分布，而且对于同一次DFT变换，不同的k，或同一个离散频率k不同次的DFT变