基于雷达的杂波速度谱图的建立方法

摘 要： 在雷达信号处理中，为了对低速运动杂波进行有效的抑制，研究了一种杂波速度谱图的建立方法。此杂波速度谱图的建立在FPGA中实现，通过对雷达实际回波数据在FPGA中的处理得到运动杂波速度图。实验结果表明，该方法设计出的杂波速度谱图与仿真结果一致，可以有效地抑制静止和慢速运动的杂波。并且由于超大规模可编程器件的高效处理能力，使其工程实现方便、灵活，具有很强的实用性。

杂波和干扰的有效抑制已经成为雷达信号处理的一个主要研究方向。现代雷达面临的杂波环境相当复杂，普通的动目标对消系统主要是对固定的地物杂波有一个很好的抑制效果，但对像气象、海浪、箔条等具有一定速度的运动杂波而言，杂波谱的中心会偏离零多普勒频率，如不采取有效措施，将无法很好地消去这类杂波[1-2]。因此考虑到杂波的运动性、区域性及时变性，为自适应地抑制低速运动杂波，建立动态杂波速度谱图是一个很好的选择。

现有的文献对杂波轮廓图和杂波幅度图研究较多[3-6]。杂波轮廓图可以用来选择正常支路和MTD(或MTI)相参处理支路，在杂波轮廓图内的数据采用相参信号处理后输出，在杂波轮廓图之外的部分采取正常支路输出，从而提高弱小目标在弱杂波或无杂波条件下的检测能力[6]。杂波幅度图可以用在恒虚警检测中。参考文献[7]和参考文献[8]对三维立体杂波图的形成与应用进行了研究。而目前对杂波速度图的研究相对较少。

1 自适应杂波抑制技术

对于运动杂波，在计算得到其多普勒频率之后，抑制杂波的方法一般有两种[9]：(1)对回波信号u(t)进行运动杂波谱中心补偿，将运动杂波谱中心移动到零频，再用凹口位于零频的MTI滤波器抑制运动杂波；(2)直接采用凹口位于fd处的MTI滤波器来抑制运动杂波，而凹口于fd处的MTI滤波器权系数可预先存储在一个滤波器权系数库中，如图1所示。

不管是对运动杂波进行多普勒中心补偿，还是利用权库法，都需要计算出运动杂波谱中心，即运动杂波的速度。因此在这里设计了一种求杂波速度谱图的方法，统计缓慢移动杂波的速度，将杂波的多普勒中心频率存储，利用FPGA的高效处理能力，通过软件编程的方法实现杂波速度谱图的建立，工程实现方便。

2 速度谱图的实现
求杂波速度的方法有多种，一种是经典的谱估计FFT法，速度精度与FFT的点数有关。在信号重频已知的情况下，FFT点数越高速度分辨精度越高。还可以用arctan求相位差的方法求杂波速度，在求解反正切函数的运算中，传统的方法有查表法、多项式近似法、查表与多项式结合法、逐位法等[10]。然而，这些方法在速度和精度上达不到要求并且硬件实现困难。

为了方便在FPGA中实现，适合FPGA中的操作流程，本设计选择互相关法进行多普勒速度的求解。