多普勒雷达测速方法改进研究

多普勒雷达测速方法改进研究

高玉苓

山南市市场监督管理局 西藏山南  856000

摘要：针对基于多普勒频偏估计的雷达测速方法精度不高的问题，提出采用自相关处理和改进的能量重心法相结合的频偏估计方法。该方法对接收信号进行自相关处理以抑制噪声的影响，并根据能量重心法中谱线选择与矫正精度的关系，对能量重心法进行改进。与常规算法相比，本算法有效提高了多普勒频偏估计的精度和多普勒雷达测速的精度。

关键词：雷达测速；算法研究；多普勒效应；能量重心法

引言

    城市轨道交通系统具有运量大、占地少、单位运量能耗小等诸多优点，成为世界各国交通建设的重点。随着以高铁和地铁为代表的轨道交通系统的快速发展，轨道交通车辆的速度检测作为列车运行控制系统的关键技术之一，是列车安全保障体系的重要组成部分日。实时、准确地获得列车速度是列车安全高效运行的重要保障。随着列车运行速度的不断提升，对列车测速的精度及可靠性也提出了更高的要求比较常用的测速方法有列车转轴转动测速、GPS(global positioning system)导航定位测速和雷达测速转动测速是一种直接接触式的测量方法，结构简单，测量误差较大。GPS导航测速是一种非接触式的测量方法，利用GPS卫星定位实时获得列车的运行速度，受天气及环境影响较大。雷达测速也是一种非接触式测速方法，根据多普勒效应，通过测量多普勒频率差，得到列车的运行速度信息。雷达测速在测速原理、精度及可靠性方面的优势使之更加适合对高速列车进行速度测量，特别是随着微处理器技术的飞速发展，雷达测速技术在列车测速领域具有广阔的应用前景。雷达测速的原理是根据列车的运行与发射信号之间的多普勒频差，解算出列车运行速度信息。因此，高精度频偏估计方法是雷达测速的关键技术之一。频偏估计方法大致可以分为时域估计法和频域估计法。时域法存在较高的信噪比门限，适用于高信噪比条件，而频域法的信噪比适用范围和频偏估计范围均较时域法广，因此得到了广泛的应用日。频域估计法基于信号周期图，通常先利用傅里叶变换获得信号的频谱峰值信息，再利用频谱矫正技术对初始估计的频率进行修正。由于能量泄露和栅栏效应，频率估计值与实际频率值存在较大误差，因此众多学者研究了各种频率校正算法。常见的主要有插值法A和能量重心法。插值法的估计精度受到频率真实谱线位置分布的影响，估计性能不理想。能量重心法是以周期图频率粗估计为基础，根据加窗后信号功率谱的能量重心逼近主瓣重心的特性，对信号频率峰值进行修正。实际中，信号叠加噪声会干扰谱线的正确选择，从而影响能量重心法频率矫正的精度。因此通过优化谱线选择可以有效提高频率估计精度。另外，学者们还把自相关技术应用到频率修正算法中，提高了现有信号频偏估计的精度。在前人研究的基础上，针对现有的基于能量重心法在频率修正方面存在的局限，本文对能量重心法进行改进，并采用自相关处理来抑制噪声，提高频偏估计的精度。

1多普勒雷达测速原理

1.1系统组成

    多普勒测速雷达通过检测多普勒频移完成对目标运动速度的检测，其主要由信号处理单元、基带、射频发射接收单元、天线和接口电路（串口）等构成，系统原理框图如图1所示，测速雷达实物图如图2所示。基带和射频发射接收单元主要产生一定频率的微波信号，并由天线将这个信号定向发射出去，同时对于反射回来的信号进行接收、混频等；信号处理单元主要完成基带信号的处理，以便获得能反映目标真实运动的多普勒频率；接口电路按照使用要求将获得的目标运动速度信息送出。

    为获得稳定可靠的多普勒频率，本系统在一部雷达上集成两个不同角度的多普勒测速模块。通过提取轨道方向的分量进行各自处理，得出最终的结果。即使某一个模块出现了问题也可以保证测速正常，能够有效提高产品的可靠性。

图1 多模块雷达测速原理框图

图2雷达测速实物图

1.2 雷达测速算法原理

雷达测速的原理是利用多普勒效应对目标速度进行测量，如图3所示，当目标与雷达存在相对运动时，会产生多普勒频移，而目标的速度大小可通过多普勒频移计算得到，计算公式为：

              （1）

式中：为多普勒频移；为雷达波的波长；为雷达天线与地面的夹角。从上式可知，提高频率估计精度可以有效提高雷达测速精度。

图3  多普勒效应示意图

2 基于自相关和改进能量重心法的频偏估计算法

设信号模型为单频复正弦信号，对信号进行采样截断后得：      （2）

其中：A为信号幅度；为信号频率为采样频率；信号功率，表示复高斯白噪声，其实部和虚部的均值为0、方差都是，所以复噪声功率为，信噪比(SNR)为E，，对信号进行N点傅里叶变换变换得：

              （3）

其中为采样周期。频谱分辨率为，幅值最大处谱线的索引记为，是最接近信号频率真实值的谱线。由于信号截断引起的能量泄露和采样引起的栅栏效应，不可避免地导致估计精确度受限。利用频谱矫正技术计算修正因子

进行补偿，可以得到矫正后频率估计值为：

     （4）

其中修正因子的实际物理意义是：信号实际频率对应谱线的位置相对于幅度峰值谱线的偏差，并且有：。

2.1 自相关处理

根据信号与噪声之间的相关特性可知，对信号进行自相关预处理，可以有效抑制噪声，提高信噪比，从而提高信号的频率估计精度。对信号进行自相关预处理后的自相关函数可以表示成：

              （5）

其中表示相对时延。从上式可以看出，自相关函数包含了信号的频率和相位信息；，表示零均值相位噪声，其概率密度分布函数为：       （6）

其中表示信噪比。的自相关序列表示为：

       （7）

其中：和分别是正弦分量和余弦分量的自相关函数；和分别是和的互相关函数，是随机噪声，和没有相关性。

当N足够大时，则和的自相关函数分别为：

       （8）

     (9)

近似有：       （10）

所以可以写成：        （11）

综上分析可见，对信号进行自相关预处理，可以有效抑制噪声，提高信号处理的信噪比，从而提高频率估计性能。

2.2  改进的能量重心法

能量重心法是根据信号功率谱的重心逼近主瓣重心的特征性为基础，对频率估计值进行修正的方法。设为信号的功率谱有：

             （12）

其中为自相关参数。

设置为能量重心法计算得到的频率修正值，则有：

              （13）

由能量重心法计算得到的修正值与谱线的选取有关，当真实频率接近量化频点时，估计精度较高，当真实频率接近量化频点中间位置时，估计精度较低。对此文献提出了一种改进算法，记为该算法在峰值谱线左右各增加一条谱线用于计算修正值，方法如下：

    （14）

其中，和，分别为峰值谱线左侧和右侧增加一根谱线计算得到的修正值，为谱线之间幅角相位差的余弦值。

能量重心法的精度不仅与谱线选择有关，还与频偏真实值的分布有关。当真实频率与峰值谱线对应的频率差位于区间[-0.2，0.2]，采用奇数根谱线能量重心法获得的估计精度较高，而位于＜-0.2和＞0.2的区间时，增加一根谱线来获得修正值可以提高频率估计精度。

所以本文对能量重心法进行改进，基于改进的能量重心法的频率估计算法描述如下：

1)对采样后的接收信号进行自相关预处理；

2)利用周期图法得到频率的粗略估计值；

3)采用常规能量重心法计算频率修正值；

4)根据的所在区间获得最终的频率估计结果。当位于[-0.2，0.2]区间上时，则直接采用作为信号频率估计的修正值。当时，如果的符号为正，增加一根峰值谱线右侧的谱线计算得到，的值作为信号频率修正值，如果的符号为负，增加一根峰值谱线左侧的谱线计算得到的值作为频率修正值。

参考文献：

[1]朱东飞，王永皎，杨烨，等.基于光栅阵列的城市轨道列车定位与测速方法.光子学报，2019，48(11)：159-166.

[2]孟凡江，尹逊政.一种基于多传感器融合的冗余测速测距系统设计.铁路计算机应用，2018，27(1)：46-50.

[3]陈文赛、高安全高可靠测速测匝系统的设计与实现口：城市轨道交通研究，2018，21(4)：131-134，147.

[4]张龙，王兴理，马婧，等.脉冲雷达测速方法的改进研究现代雷达，2019，41(10)：74-77.

[5]万灵达，杨晓光.一种基于高精度频率估计算法电子科技，2010，23(10)：79-81.