一种基于自相关函数和DTFT的频率估计方法

本发明涉及通信，尤其涉及一种基于自相关函数和dtft的频率估计方法。

背景技术：

1、加性白噪声背景中的正弦信号频率估计是通信领域的一个重要的研究课题。例如利用正弦信号频率估计能够减小或消除多普勒效应的影响。多普勒频移是通信中经常遇到的问题，多普勒效应不仅会降低码元同步精度，而且会带来载波频率的偏移。所以对载波频偏的准确估计是解决通信问题的关键步骤。目前，在移动通信系统和卫星通信系统中，都需要先对多普勒效应导致的载波频移进行准确估计，再结合估计结果进行频率校正。因此对正弦信号频率估计的研究有着重要的理论意义和重要价值。

2、时域估计算法和频域估计算法是现有正弦信号频率估计的两大类。时域频率估计算法有最大似然算法、自相关算法、线性预测算法以及最小二乘算法等。频域估计算法是利用离散傅里叶变换进行正弦信号频率估计，该类方法与时域估计方法相比，复杂度相对较低，在低信噪比时的性能更好。

3、dft频率估计方法，物理意义明确，具有较高的信噪比增益和对参数不敏感等优点，因此已经成为正弦信号频率估计领域一种广受关注与研究的方法。近年来许多基于dft的频率估计方法为了提高估计精度进行额外的步骤从而提升了计算复杂度。这些算法在估计精度和计算复杂度方面不能达到很好的平衡。当前基于dft的频率估计算法都需要先进行fft运算，fft运算占据了算法运算量的大部分，对正弦信号频率估计进行实时处理增加了障碍。

技术实现思路

1、本发明主要解决现有技术的计算复杂度高和估计精度低的技术问题，提出一种基于自相关函数和dtft的频率估计方法，在较大降低运算量的同时，提高了估计精度。

2、本发明提供了一种基于自相关函数和dtft的频率估计方法，包括以下过程：

3、步骤1，采集待处理的复正弦信号，获得n个样本的观测序列；

4、步骤2，利用自相关函数进行频率估计，得到初始频率fr；

5、步骤3，在初始频率附近选择三个抽样点，并计算这三个抽样点的dtft抽样值，根据三个抽样点的dtft抽样值求出相对频率偏差δ，得到精估计频率fδ；

6、步骤4，利用步骤2得到的初始频率和步骤3得到的精估计频率得到复正弦信号频率估计。

7、进一步的，在步骤4之后，还包括：

8、步骤5，利用迭代运算，得到新的精估计频率更新复正弦信号频率估计。

9、进一步的，在步骤1中，加性白噪声背景下的复正弦信号的观测序列可表示为：

10、x(n)＝s(n)+z(n), n＝0,1,2,…,n-1 (1)

11、

12、其中，z(n)为加性高斯白噪声，a为信号的幅度、f0为信号频率、θ0为初相、fs为采样频率、n为采样点数，j表示虚数单位。

13、进一步的，步骤2包括如下过程：

14、在相位随机变化的情况下，观测序列估计复正弦信号频率归结为求使下式最大似然方程达到最大的初始频率fr；

15、

16、其中，δt表示采样间隔，x(i-1)表示复正弦信号，i表示采样点；

17、定义非归一化的相关函数，相关函数表示为：

18、

19、其中，vk为等效噪声，k表示加权系数；

20、在初始频率与采样间隔乘积frδt较小的情况下，相关函数r(k)的相位随加权系数k的增加近似线性增加，对于一定频率的信号，利用相关函数r(k)估计初始频率fr时加权系数k越大估计误差越小；因此arg{e[r(k)]}＝2πfrkδt；

21、在采样点数n较大，初始频率与采样间隔乘积frδt较小的情况下，有：

22、

23、整理化简得到初始频率fr的近似估计值为：

24、

25、其中，表示估计出来的初始频率，m表示采样点数；

26、用m项相关函数估计频率，有：

27、

28、为了便于计算和实时处理，省略加权系数，此时频率为：

29、

30、m等于1可以得到初始频率的公式：

31、

32、按照公式(9)，计算得到初始频率

33、进一步的，在步骤3中，抽样点的选择规则：初始频率附近任意位置以及距离该任意点左右间隔为长度阈值的两点，分别为第一抽样点、第二抽样点和第三抽样点；其中，长度阈值为

34、三个抽样点的dtft抽样值，分别为第一抽样值sg、第二抽样值sg-h、第三抽样值sg+h；其中，g为偏离初始频率的离散频率值；

35、根据第一抽样值sg、第二抽样值sg-h、第三抽样值sg+h求出相对频率偏差δ，得到精估计频率fδ。

36、进一步的，步骤3包括如下步骤3.1至步骤3.3：

37、步骤3.1，在0≤f＜fs时，计算序列单一频率复正弦信号的信号模型s(n)在任意频率f处的频谱s(f)：

38、

39、其中，n为采样点数，本发明用δ∈[-0.5,0.5]表示相对频率偏差，信号频率表示为

40、

41、步骤3.2，把和s(n)代入s(f)，a表示与初始频率估计值的距离间隔，可以得到任意位置的dtft谱线sa的表达式：

42、

43、步骤3.3，令a分别等于g，g-h，g+h时，得到第一抽样值sg、第二抽样值sg-h、第三抽样值sg+h的表达式，具体步骤如下：

44、当n的值远远大于π(δ-a)时，

45、将代入dtft谱线sa的表达式得到新的dtft谱线sa表达式：

46、

47、当a分别等于g，g-h，g+h时，可以得到第一抽样值sg、第二抽样值sg-h、第三抽样值sg+h的表达式：

48、

49、

50、

51、经过一些数学运算可以得到：

52、

53、整理得到相对频率偏差δ的估计表达式为：

54、

55、其中，表示估计出来的相对频率偏差；

56、精估计频率fδ的估计表达式为：

57、

58、其中，表示估计出来的精估计频率；

59、利用公式(18)计算得到精估计频率

60、本发明提供的一种基于自相关函数和dtft的频率估计方法，与现有技术相比具有以下优点：

61、1、首先采用信号采样序列的相关函数得到频率的初始频率，然后采用位于初始频率值附近的三个dtft抽样值得到频率的精估计频率，并利用初始频率和精估计频率得到复正弦信号频率估计；还可以进一步利用迭代运算，得到精度更高的复正弦信号频率估计。

62、2、本发明与同类型方法比较，本发明的算法复杂度最低，估计精度最高，而通信领域中的应用对算法实时性和频率估计精度的要求一般都很高。相比于其他方法，本发明更加适用于通信领域。

63、3、本发明使用的三个dtft抽样点的位置是任意的，因此本发明更具有一般性。能够更加灵活方便得解决通信领域中的频率估计问题。