基于DFT相位的正弦波频率和初相的估计方法

基于DFT相位的正弦波频率和初相的估计方法摘要：利用分段DFT频谱的相位差消除了初相对频率估计的影响且避免了相位模

糊问题。给出了频率和初相估计的均方根误差计算公式。理论分析和Monte Carlo模拟结果显示频率估计均方根误差接近Cramer-Rao(CR)下限，初相估计均方根误差略高于CR下限的2倍。阈值信噪比远远低于基于时域瞬时相位的频率和初相估计方法。在信噪比为6dB，采样点数为1024的情况下，频率估计均方根误差约为DFT频率分辨率的1%，初相估计均方根误差约为2度。该方法已用于FMCW液位测量雷达并取得1mm距离测量精度。

关键词：频率估计；相位估计；雷达测距；DFT

Abstract：Use block of phase DFT spectrum to eliminate the influence of relative frequency

early estimates and avoid the phase fuzzy problem. Given the frequency and initial phase estimate of the RMS error calculation formula. The theory analysis and Monte Carlo simulation results indicate that the frequency estimation RMS error close to Cramer-Rao (CR) floor, the RMS error estimation is slightly higher than the lower limit CR 2 times. Threshold value signal-to-noise ratio is much lower than the instantaneous phase frequency based on the time domain and initial phase estimation. In the signal-to-noise ratio is 6 dB, sampling points for 1024, frequency estimation RMS error of about 1% of the DFT frequency resolution, early phase estimate RMS error of about 2 degrees. The method has been used to level measurement radar and obtain FMCW 1 mm distance measurement accuracy.

Keywords：Frequency estimation、Phase estimates、Radar ranging、DFT

1 引言

正弦波信号频率的估计是通信、雷达、声纳以及电子对抗等领域信号处理中的一个重要问题。基于参数模型的谱估计、最大熵谱估计等方法具有频率分辨率高的优点，但对于长序列，其运算量大，不利于实时处理。而基于DFT的谱分析方法，可采用快速算法，即FFT，因此运算速度快，特别适合于实时信号处理。但DFT的频率分辨力和频率估计精度取决于信号的测量时间长度，信号测量时间过长不但给实时处理带来困难，而且在一些应用中信号的持续时间是有限的，不能任意延长，使得DFT的频率分辨力和估计精度受到限制，因此一般只利用DFT实现频率的粗测[1,2]，文献[3]提出了利用数值搜索方法提高DFT频率估计精度。由于栅栏效应DFT频谱在主瓣之内有2条谱线，利用这2条谱线的幅度可以实现频率插值以及提高频率测量精度[4]。文献[5]提出了直接在时域采用最小二乘线性回归的方法，利用瞬时相位估计信号频率和初相。为了避免直接测量瞬时相位的整周模糊问题，文献[6]提出了相位查分方法。在时域测量瞬时相位估计信号频率和初相的缺点是需要较高的信噪比。文献[7]提出了一种利用DFT频谱的相位和频率插值的综合算法提

高频率测量精度的方法。本文提出一种利用信号DFT 频谱相位提高频率测量精度的新方法，利用分段DFT 消除了相位测量中的整周模糊问题。理论分析和计算机模拟结果表明，该方法的信噪比阈值比基于时域瞬时相位的频率估计方法低很多，在很低信噪比情况下，仍能得到很高的频率和初相估计精度，频率估计均方根误差接近CR 下限。该方法还可以估计信号的初相，初相估计均方根误差高于CR 下限的2倍。

2 FFT 相位差法估计频率和初相原理

设观测信号为单一频率复正弦波信号

00()exp[(2)]s t a j f t πφ=⋅+ (1)

其中，a 、0f 、0φ分别为信号的幅度、频率和初相。对上述信号进行采样，设信号的记录时间长度分为两个长度相同的序列，1()s n 对应前/2N 点，则采样序列可记为

00()exp[(2)],0,1,2,s t a j f t n πφ=⋅+=…,N-1 (2)

2()s n 对应后/2N 点，则

100()exp[(2)],0,1,2,s t a j f t n πφ=⋅+=…,N/2-1 (3)

210()()exp(),0,1,2,s n s n j f T n π==…,N/2-1 (4)

分别对1()s n 和2()s n 进行/2N 点DFT ，得到离散频谱

1()exp(),0,1,2,k k S k A j k φ==…,(N/2-1) (5)

210()()exp(),0,1,2,S k S k j f T k π==…,(N/2-1) (6)

其中k A 和k φ分别为1()S k 的幅度项与相位项

00sin[(/2)]sin[2(/2)/]

k a k f T A k f T N ππ⋅-=- (7) k 00(12/)(/2)N f T k ϕφπ=+-- (8)

根据式(6)知，1()S k 与2()S k 的幅度项完全一样。由式(7),幅度最大处值对应的离散频率为

00[/2]k f T = ([x]表示取最接近x 的整数)。利用DFT 的最大谱线粗测频率为0

ˆk f k f =∆，2/f T ∆=为DFT 的频率分辨率。从式(8)可见，DFT 最大谱线的相位包含信号频率与DFT 最大谱线位置的偏差信息，但由于0φ未知，不能直接利用DFT 的相位来估计频率。用1ϕ和2ϕ分别表示1()S k 和2()S k 在最大谱线处的相位，则两者的差值为

21002f T k ϕϕϕππ∆=-=- (9)