CN108845973A - 一种基于改进Quinn算法的多普勒频率估计方法 - Google Patents

Abstract

本发明公开了一种基于改进Quinn算法的多普勒频率估计方法，在利用奎因Quinn算法进行多普勒频率估计前，对信号进行补零的傅里叶变换FFT运算，使得信号频谱主瓣内的谱线增多，进而利用更多的谱线信息，对峰值位置进行修正，能获得更准确的多普勒频率估计值。该方法包括如下步骤：对N点的扩频信号x(n)进行一倍补零，对一倍补零后的信号做点数为2N的FFT运算，得到频域信号X(k)；选取频域信号X(k)幅度的峰值|X(k₀)|，k₀为峰值|X(k₀)|对应的频域位置坐标；采用Quinn算法进行频域插值，判断Quinn算法插值方向r，其中r＝+1表示插值方向向右进行；r＝‑1表示插值方向向左进行；同时计算Quinn算法的谱线插值搬移的移动量计算多普勒频率估计值其中Δf为FFT运算的鉴频间隔。

Description

一种基于改进Quinn算法的多普勒频率估计方法

技术领域

本发明涉及扩频通信技术领域，具体涉及一种基于改进Quinn算法的多普勒频率估计方法。

背景技术

高动态扩频信号的捕获需要完成对多普勒频率的估计，提高多普勒频率估计精度、减小估计误差可以使跟踪单元快速锁定。其估计精度受限于接收信号的信噪比，信噪比越高，多普勒频率估计精度越高。通常会采用延长相关时间的方式来提高信噪比，但是在高动态低载噪比环境下，延长相关处理时间会带来多普勒频率徙动的问题，不仅不能有效提高信噪比而且存在多普勒频率扩散，无法有效提高多普勒频率估计精度。

目前，扩频信号捕获技术中获得高精度多普勒频率估计值的方法主要有精频捕获方法和频域估计方法两种。

精频捕获方法可分为提高傅里叶变换FFT鉴频精度、改进多普勒频率搜索策略两个方面。提高FFT的鉴频精度往往从增加FFT点数和频谱细化入手。增加FFT点数的方法原理简单直观，但会使捕获时间延长、运算量增大。频谱细化的方法主要有Zoom-FFT、线性调频Z变换：Zoom-FFT通过降低采样频率和频率抽取等方式，将目标频段进行局部放大；线性调频Z变换是在Z变换的基础上，在某一频率附近增加谱线密度，这两种方法实现复杂均较高。搜索策略方面主要是缩小步长搜索方法，即通过多轮多普勒频率搜索，每一轮的频率搜索步进在上一次的基础上减半，来达到提高捕获多普勒频率估计精度的目的，但是该方法捕获时间长、效率低下。

频域估计方法在捕获的粗多普勒频率估计值的基础上，通过相位差分、数值逼近或者频域插值等方法，进一步对小范围内的频率进行处理，获得更高的估计精度，在运算复杂度方面优势明显。其中，基于相位差分的频率估计方法利用相位变化信息对频率进行估计，存在相位模糊的问题，影响频率估计准确性。数值逼近法是利用曲线拟合的方式对sinc形鉴频曲线进行还原，其拟合曲线截断长度的选取难以兼顾拟合误差和运算复杂度性能。频域插值方法有Rife、Mife、和奎因Quinn及其衍生算法等，该类方法利用FFT后主瓣和第一旁瓣的谱线的位置、幅值和相位等信息提高频率估计精度。Rife算法利用了主瓣和第一旁瓣的谱线幅值相对关系进行插值，但是在噪声环境下，当实际多普勒频率与估计多普勒频率之间的差值较小时，即实际信号多普勒频率靠近离散频点时，可能会出现插值方向错误，导致频率估计误差较大。M-Rife算法通过频移处理，避免信号落入偏差较大的区域，能够较好地解决Rife算法插值方向错误的问题，但是该算法存在迭代过程，运算复杂度大。常规Quinn算法在Rife算法基础上，通过引入相位信息来判断插值方向，可解决Rife算法的插值方向错误，且运算复杂度低，但是其只使用了主瓣内两根谱线实部的信息，能达到的频率估计精度有限。所以研究具有高估计精度的Quinn算法及其多普频率估计方法具有实际意义。

目前尚未有一种方法能够解决上述问题。

发明内容

有鉴于此，本发明提供了一种基于改进Quinn算法的多普勒频率估计方法，在利用奎因Quinn算法进行多普勒频率估计前，对信号进行补零的FFT运算，使得信号频谱主瓣内的谱线增多，进而利用更多的谱线信息，对峰值位置进行修正，能获得更准确的多普勒频率估计值。

为达到上述目的，本发明的技术方案提供了一种基于改进Quinn算法的多普勒频率估计方法，包括如下步骤：

步骤(1)对N点的扩频信号x(n)进行一倍补零，补零数量为N，对一倍补零后的信号做点数为2N的傅里叶变换FFT运算，得到频域信号X(k)。

选取频域信号X(k)幅度的峰值|X(k₀)|，k₀为峰值|X(k₀)|对应的频域位置坐标。

步骤(2)采用奎因Quinn算法进行频域插值，判断Quinn算法插值方向r，其中r＝+1表示Quinn算法插值方向向右进行；r＝-1表示Quinn算法插值方向向左进行；同时计算Quinn算法的谱线插值搬移的移动量

判断Quinn算法插值方向r采用如下方法：

使用k₀左右两个频率采样点k₀-2和k₀+2处的谱线X(k₀-2)和X(k₀+2)，计算Quinn算法插值方向判断参数及

其中，Re{·}表示取实部。

则r为：

计算Quinn算法的谱线插值搬移的移动量采用如下方法：

比较k₀左右两边的两个频率采样点k₀+1和k₀-1处的谱线幅值|X(k₀+1)|和|X(k₀-1)|，选取|X(k₀+1)|和|X(k₀-1)|的较大者记为|X(k₁)|。

估计频率和实际信号频率偏差的估计值，即谱线插值搬移的移动量为

步骤(3)多普勒频率估计值为：

其中Δf为FFT运算的鉴频间隔，Δf＝f_s/(2N)，f_s为采样频率。

有益效果：

本方法在常规Quinn算法的基础上，使用补零FFT的方式增加主瓣内谱线数目，进而利用更多谱线信息进行多普勒频率估计，在Quinn算法的插值方向进行判断时，可以使用峰值两个频率采样点k₀-2和k₀+2处的谱线X(k₀-2)和X(k₀+2)对插值的方向做判断；相比已有的精频捕获方法，本发明方法运算复杂度低、耗时短；同时本发明方法采用次峰值进行谱线插值搬移的移动量的计算，由于经过补零，频谱分辨率提高，次峰值位置更加接近峰值位置，因此搬移量计算也就更精确，有效提高了多普勒频率估计精度。

附图说明

图1为本发明实施例提供的基于改进Quinn算法的多普勒频率估计方法流程图；

图2为本发明实施例提供的基于改进Quinn算法的多普勒频率估计方法中步骤(2)展开后的流程图。

具体实施方式

下面结合附图并举实施例，对本发明进行详细描述。

本发明提供了一种基于改进Quinn算法的多普勒频率估计方法，具体流程如图1所示，包括如下步骤：

步骤(1)对N点的扩频信号x(n)进行一倍补零，补零数量为N，对一倍补零后的信号做点数为2N的傅里叶变换FFT运算，得到频域信号X(k)。

选取频域信号X(k)幅度的峰值|X(k₀)|，k₀为峰值|X(k₀)|对应的频域位置坐标。

本发明使用补零FFT的方式增加主瓣内谱线数目，进而利用更多谱线信息进行多普勒频率估计，能够降低运算复杂度，提高多普勒频率的估计精度。

步骤(2)采用奎因Quinn算法进行频域插值，判断Quinn算法插值方向r，其中r＝+1表示Quinn算法插值方向向右进行；r＝-1表示Quinn算法插值方向向左进行；同时计算Quinn算法的谱线插值搬移的移动量

图2示出了针对步骤(2)的展开步骤。

由于步骤(1)的补零操作，可以利用更多谱线信息进行多普勒频率估计，所以在判断Quinn算法插值方向r时，可以使用峰值两个频率采样点k₀-2和k₀+2处的谱线X(k₀-2)和X(k₀+2)对插值的方向做判断。

判断Quinn算法插值方向r采用如下方法，如图2所示：

使用k₀左右两个频率采样点k₀-2和k₀+2处的谱线X(k₀-2)和X(k₀+2)，计算Quinn算法插值方向判断参数及

其中，Re{·}表示取实部。

则r为：

计算Quinn算法的谱线插值搬移的移动量采用如下方法，如图2所示：

比较k₀左右两边的两个频率采样点k₀+1和k₀-1处的谱线幅值|X(k₀+1)|和|X(k₀-1)|，选取|X(k₀+1)|和|X(k₀-1)|的较大者记为|X(k₁)|。

估计频率和实际信号频率偏差的估计值，即谱线插值搬移的移动量为

本发明方法采用次峰值进行谱线插值搬移的移动量的计算，改善了常规Quinn算法偏移量误差大的问题，可有效提高多普勒频率估计精度。

步骤(3)多普勒频率估计值为：

其中Δf为傅里叶变换FFT运算的鉴频间隔，Δf＝f_s/(2N)，f_s为采样频率。

综上，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种基于改进Quinn算法的多普勒频率估计方法，其特征在于，包括：

步骤(1)对N点的扩频信号x(n)进行一倍补零，补零数量为N，对一倍补零后的信号做点数为2N的傅里叶变换FFT运算，得到频域信号X(k)；

选取所述频域信号X(k)幅度的峰值|X(k₀)|，k₀为峰值|X(k₀)|对应的频域位置坐标；

步骤(2)采用奎因Quinn算法进行频域插值，判断Quinn算法插值方向r，其中r＝+1表示Quinn算法插值方向向右进行；r＝-1表示Quinn算法插值方向向左进行；同时计算Quinn算法的谱线插值搬移的移动量

所述判断Quinn算法插值方向r采用如下方法：

使用k₀左右两个频率采样点k₀-2和k₀+2处的谱线X(k₀-2)和X(k₀+2)，计算Quinn算法插值方向判断参数及

其中，Re{·}表示取实部；

则r为：

所述计算Quinn算法的谱线插值搬移的移动量采用如下方法：

比较k₀左右两边的两个频率采样点k₀+1和k₀-1处的谱线幅值|X(k₀+1)|和|X(k₀-1)|，选取|X(k₀+1)|和|X(k₀-1)|的较大者记为|X(k₁)|；

估计频率和实际信号频率偏差的估计值，即谱线插值搬移的移动量为

步骤(3)多普勒频率估计值为：

其中Δf为所述FFT运算的鉴频间隔，Δf＝f_s/(2N)，f_s为采样频率。