用于雷达探测的恒虚警检测方法、装置及电子设备

摘要

本发明提供一种用于雷达探测的恒虚警检测方法、装置、电子设备及存储介质，属于雷达信号处理技术领域，所述方法包括：获取至少一个雷达帧的样本数据并生成对应的距离‑多普勒散射中心能量矩阵；针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并计算每个距离单元的目标噪声能量估计系数以及每个距离单元的目标检测门限系数，并计算每个距离单元的参考检测门限值，并基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚警检测。本发明可实现在非均匀杂波分布条件下的目标检测，并且能够达到理想的检测性能。

权利要求

1.一种用于雷达探测的恒虚警检测方法，其特征在于，所述方法包括：获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵；

针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数；

针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数；

在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚警检测。

2.根据权利要求1所述的用于雷达探测的恒虚警检测方法，其特征在于，所述获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵的步骤包括：

将每个雷达帧的样本数据对应的每个距离单元的多普勒维的散射中心按照能量由小到大进行排序，生成该雷达帧对应的距离‑多普勒散射中心能量矩阵 ，表示为：= ；

其中，M表示多普勒维傅立叶变换的点数，N表示距离维傅立叶变换的点数；

其中，第p个距离单元的散射中心能量排序后的向量 表示为：。

3.根据权利要求2所述的用于雷达探测的恒虚警检测方法，其特征在于，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤包括：针对每个雷达帧，以每个距离单元对应的散射中心能量排序中的中值作为该距离单元的噪声能量估计的真值；

其中，设第p个距离单元的噪声能量估计的真值为 ，则 ， 。

4.根据权利要求3所述的用于雷达探测的恒虚警检测方法，其特征在于，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤还包括：针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则。

5.根据权利要求4所述的用于雷达探测的恒虚警检测方法，其特征在于，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤还包括：针对每个雷达帧，设每个距离单元的噪声能量估计系数的初始值为 ，最大值为 ，搜索步长为 ，则第p个距离单元的第 次搜索的噪声能量估计系数为 ；

以   为预设门限，设能量小于所述预设门限的散射中心为 ，其中，， ；

设第p个距离单元的的第 次搜索的噪声能量估计值 ，则;

针对每个距离单元搜索噪声估计系数使得对应的噪声能量估计值与真值之间的差值最小，其中，设第p个距离单元的噪声能量估计系数为 ，则。

6.根据权利要求5所述的用于雷达探测的恒虚警检测方法，其特征在于，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤还包括：针对每个雷达帧，根据每个距离单元的噪声能量估计的真值计算该距离单元的检测门限系数；

其中，设第p个距离单元的多普勒维散射中心能量由小到大排序的预设位置的散射中心能量作为第p个距离单元的噪声能量最大值的估计值 ，并且设第p个距离单元的检测门限系数为 ，则

；

其中，表示固定比例系数。

7.根据权利要求6所述的用于雷达探测的恒虚警检测方法，其特征在于，所述针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数的步骤包括：若所述至少一个雷达帧仅包括一个雷达帧，则将基于所述一个雷达帧计算出的每个距离单元的噪声能量估计系数作为该距离单元对应的目标噪声能量估计系数，并且将基于所述一个雷达帧计算出的每个距离单元的检测门限系数作为该距离单元对应的目标检测门限系数；

若所述至少一个雷达帧包括多个雷达帧，则将不同雷达帧的同一距离单元的噪声能量估计系数的均值作为该距离单元对应的目标噪声能量估计系数，并且将不同雷达帧的同一距离单元的检测门限系数的均值作为该距离单元对应的目标检测门限系数。

8.根据权利要求7所述的用于雷达探测的恒虚警检测方法，其特征在于，所述在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚警检测的步骤包括：针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则；

计算每个距离单元的噪声能量的估计值，其中，设第p个距离单元的目标噪声能量估计系数为 ，并且能量小于 的散射中心为 ， ， ，则第p个距离单元的噪声能量的估计值 被表达为：；

计算每个距离单元的参考检测门限值，其中，第p个距离单元的参考检测门限值为 *， 为第p个距离单元的目标检测门限系数。

9.一种用于雷达探测的恒虚警检测装置，其特征在于，包括：矩阵生成模块，用于获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵；

第一计算模块，用于针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数；

第二计算模块，用于针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数；

恒虚警检测模块，用于在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚警检测。

10.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至8中任一项所述的用于雷达探测的恒虚警检测方法的步骤。

11.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至8中任一项所述的用于雷达探测的恒虚警检测方法的步骤。

说明书

用于雷达探测的恒虚警检测方法、装置及电子设备

技术领域

[0001] 本发明涉及雷达信号处理技术领域，尤其涉及一种用于雷达探测的恒虚警检测方法、装置及电子设备。

背景技术

[0002] 恒虚警检测（Constant False Alarm Rate, CFAR）技术是雷达系统在保持虚警概率恒定条件下对接收机输出的信号与噪声作判别以确定目标信号是否存在的技术。

[0003] 在雷达大量安装使用过程中，面临复杂的应用场景，杂波或噪声能量的分布呈现多样化，经典的恒虚警检测算法（CFAR）仅适用于均匀或者部分非均匀杂波分布条件下的目
标检测，在非均匀杂波条件下无法达到理想的检测性能。

发明内容

[0004] 本发明提供一种用于雷达探测的恒虚警检测方法、装置、电子设备及存储介质，用以解决现有技术中在非均匀杂波条件下无法达到理想的检测性能的问题。

[0005] 本发明提供一种用于雷达探测的恒虚警检测方法，所述方法包括：

[0006] 获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵；

[0007] 针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距
离单元的噪声能量估计系数以及每个距离单元的检测门限系数；

[0008] 针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数；

[0009] 在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并
基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚
警检测。

[0010] 根据本发明所述的用于雷达探测的恒虚警检测方法，所述获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵的
步骤包括：

[0011] 将每个雷达帧的样本数据对应的每个距离单元的多普勒维的散射中心按照能量由小到大进行排序，生成该雷达帧对应的距离‑多普勒散射中心能量矩阵 ，表示为：

[0012] = ；

[0013] 其中，M表示多普勒维傅立叶变换的点数，N表示距离维傅立叶变换的点数；

[0014] 其中，第p个距离单元的散射中心能量排序后的向量 表示为：

[0015] 。

[0016] 根据本发明所述的用于雷达探测的恒虚警检测方法，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根
据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及
每个距离单元的检测门限系数的步骤包括：

[0017] 针对每个雷达帧，以每个距离单元对应的散射中心能量排序中的中值作为该距离单元的噪声能量估计的真值；

[0018] 其中，设第p个距离单元的噪声能量估计的真值为 ，则 ， 。

[0019] 根据本发明所述的用于雷达探测的恒虚警检测方法，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根
据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及
每个距离单元的检测门限系数的步骤还包括：

[0020] 针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则

[0021] 。

[0022] 根据本发明所述的用于雷达探测的恒虚警检测方法，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根
据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及
每个距离单元的检测门限系数的步骤还包括：

[0023] 针对每个雷达帧，设每个距离单元的噪声能量估计系数的初始值为 ，最大值为 ，搜索步长为 ，则第p个距离单元的第 次搜索的噪声能量估计系数为 ；

[0024] 以  为预设门限，设能量小于所述预设门限的散射中心为 ，其中，， ；

[0025] 设第p个距离单元的的第 次搜索的噪声能量估计值 ，则

[0026] ;

[0027] 针对每个距离单元搜索噪声估计系数使得对应的噪声能量估计值与真值之间的差值最小，其中，设第p个距离单元的噪声能量估计系数为 ，则

[0028] 。

[0029] 根据本发明所述的用于雷达探测的恒虚警检测方法，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根
据所述每个距离单元的噪声能量估计的真值，计算每个距离单元的噪声能量估计系数以及
每个距离单元的检测门限系数的步骤还包括：

[0030] 针对每个雷达帧，根据每个距离单元的噪声能量估计的真值计算该距离单元的检测门限系数；

[0031] 其中，设第p个距离单元的多普勒维散射中心能量由小到大排序的预设位置的散射中心能量作为第p个距离单元的噪声能量最大值的估计值 ，并且设第p个距离单元的
检测门限系数为 ，则

[0032] ；

[0033] 其中，表示固定比例系数。

[0034] 根据本发明所述的用于雷达探测的恒虚警检测方法，所述针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应
的目标噪声能量估计系数和目标检测门限系数的步骤包括：

[0035] 若所述至少一个雷达帧仅包括一个雷达帧，则将基于所述一个雷达帧计算出的每个距离单元的噪声能量估计系数作为该距离单元对应的目标噪声能量估计系数，并且将基
于所述一个雷达帧计算出的每个距离单元的检测门限系数作为该距离单元对应的目标检
测门限系数；

[0036] 若所述至少一个雷达帧包括多个雷达帧，则将不同雷达帧的同一距离单元的噪声能量估计系数的均值作为该距离单元对应的目标噪声能量估计系数，并且将不同雷达帧的
同一距离单元的检测门限系数的均值作为该距离单元对应的目标检测门限系数。

[0037] 根据本发明所述的用于雷达探测的恒虚警检测方法，所述在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标
检测门限系数，计算每个距离单元的参考检测门限值，并基于每个距离单元的参考检测门
限值对所述雷达探测数据的多普勒维散射中心进行恒虚警检测的步骤包括：

[0038] 针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则

[0039] ；

[0040] 计算每个距离单元的噪声能量的估计值，其中，设第p个距离单元的目标噪声能量估计系数为 ，并且能量小于 的散射中心为 ， ， ，

[0041] 则第p个距离单元的噪声能量的估计值 被表达为：

[0042] ；

[0043] 计算每个距离单元的参考检测门限值，其中，第p个距离单元的参考检测门限值为* ， 为第p个距离单元的目标检测门限系数。

[0044] 本发明还提供一种用于雷达探测的恒虚警检测装置，包括：

[0045] 矩阵生成模块，用于获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵；

[0046] 第一计算模块，用于针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计
的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数；

[0047] 第二计算模块，用于针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测
门限系数；

[0048] 恒虚警检测模块，用于在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元
的参考检测门限值，并基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒
维散射中心进行恒虚警检测。

[0049] 本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述的用于雷达探
测的恒虚警检测方法的步骤。

[0050] 本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述的用于雷达探测的恒虚警检测方法的步骤。

[0051] 本发明提供的用于雷达探测的恒虚警检测方法、装置、电子设备及存储介质，通过设置对应距离单元的噪声能量估计值和对应距离单元的检测门限系数的乘积作为参考检
测门限值，并以所述参考检测门限值对多普勒维散射中心进行检测，实现在非均匀杂波分
布条件下的目标检测，并且能够达到理想的检测性能。

附图说明

[0052] 为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一
些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些
附图获得其他的附图。

[0053] 图1是本发明提供的用于雷达探测的恒虚警检测方法的流程示意图；

[0054] 图2是本发明提供的确定噪声能量估计的真值的流程示意图；

[0055] 图3是本发明提供的计算噪声能量估计值的流程示意图；

[0056] 图4是本发明提供的计算最佳噪声能量估计系数的流程示意图；

[0057] 图5是本发明提供的计算检测门限系数的流程示意图；

[0058] 图6是本发明提供的计算参考检测门限值的流程示意图；

[0059] 图7（a）图7（d）是本发明提供的应用示例的示意图；~

[0060] 图8是本发明提供的用于雷达探测的恒虚警检测装置的结构示意图；

[0061] 图9是本发明提供的电子设备的结构示意图。

具体实施方式

[0062] 为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，
而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳
动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0063] 本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情
况下可以互换，以便这里描述的实施例能够以除了在这里图示或描述的内容以外的顺序实
施。

[0064] 以下对本发明涉及的技术术语进行描述：

[0065] 雷达目标的散射中心一般位于目标边缘、拐点及联接处，其反映目标精密的几何结构特征，在目标识别中具体有重要的作用和意义。

[0066] 由于雷达工作的电磁环境越来越复杂，充斥着各种杂波、噪声及干扰，如何从回波中提取有效信息，准确、快速地检测到目标至关重要。雷达信号处理常用的方法有脉冲压
缩、积累、动目标、旁瓣相消、数字波束形成、恒虚警检测等。

[0067] 恒虚警（constant false‑alarm，CFAR）是雷达信号处理的重要组成部分。雷达系统通常要求能够在比热噪声更为复杂和不确知的背景环境中检测目标的存在并保持给定
的虚警概率，为此，必须采用自适应门限检测电路。保持信号检测时的虚警概率，使漏检概
率达到最小，或使正确检测概率达到最大。恒虚警处理的基本过程是估计需要检测的单元
中的噪声和干扰电平，并根据估计值设置阈值，然后与检测单元信号进行比较以确定是否
存在目标。

[0068] 本发明提供的用于雷达探测的恒虚警检测方法、装置、电子设备及存储介质，用以解决现有技术中在非均匀杂波条件下无法达到理想的检测性能的问题。通过设置对应距离
单元的噪声能量估计值和对应距离单元的检测门限系数的乘积作为参考检测门限值，并以
所述参考检测门限值对多普勒维散射中心进行检测，实现在非均匀杂波分布条件下的目标
检测，并且能够达到理想的检测性能。

[0069] 下面结合图1‑图9描述本发明的用于雷达探测的恒虚警检测方法、装置、电子设备及存储介质。

[0070] 图1是本发明提供的用于雷达探测的恒虚警检测方法的流程示意图，如图1所示。一种用于雷达探测的恒虚警检测方法，所述方法包括：

[0071] 步骤101，获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵。

[0072] 步骤102，针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计
算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数。

[0073] 步骤103，针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数。

[0074] 步骤104，在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门
限值，并基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进
行恒虚警检测。

[0075] 本发明所述用于雷达探测的恒虚警检测方法的核心是如何准确地估计出目标噪声能量估计系数和目标检测门限系数（即步骤103），以下针对上述步骤101 104进行具体描
~
述。

[0076] 本发明所述用于雷达探测的恒虚警检测方法包括两部分：

[0077] 第一部分是参考检测门限值的自学习模式（即对应上述步骤101 103）；~

[0078] 第二部分是参考检测门限值的实际检测模式（即对应上述步骤104）。

[0079] 通常在雷达初次安装完毕后则进入第一部分所述的参考检测门限值的自学习模式，如下：

[0080] 上述步骤101中，所述雷达信号的样本数据可通过测量获得，通过研究各种目标的多普勒维散射中心能量矩阵可得到目标的散射特性。

[0081] 可选地，上述步骤101中，所述获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵的步骤包括：

[0082] 将每个雷达帧的样本数据对应的每个距离单元的多普勒维的散射中心按照能量由小到大进行排序，生成该雷达帧对应的距离‑多普勒散射中心能量矩阵 ，表示为：

[0083] = ；

[0084] 其中，M表示多普勒维傅立叶变换的点数，N表示距离维傅立叶变换的点数。

[0085] 图2是本发明提供的确定噪声能量估计的真值的流程示意图，如图2所示。上述步骤102中，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距
离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每
个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤包括：

[0086] 步骤201，针对每个雷达帧，以每个距离单元对应的散射中心能量排序中的中值作为该距离单元的噪声能量估计的真值。

[0087] 步骤202，设第p个距离单元噪声能量估计的真值为 ，其中 ， 。

[0088] 上述包含第p个距离单元的多普勒维散射中心能量矩阵 表示为：

[0089] = 。

[0090] 则第p个距离单元的散射中心能量排序后的向量 表示为：

[0091] 。

[0092] 需要说明的是，上述步骤201中可选择每个距离单元对应的多普勒维散射中心能量矩阵 的平均值或中值作为该距离单元的噪声能量估计的真值。但是采取中值的适应
性会更好些，因为在实际雷达系统中有可能会出现较大的野值（比如：同频异步干扰）。

[0093] 图3是本发明提供的计算噪声能量的估计值的流程示意图，如图3所示。上述步骤102中，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离
单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个
距离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤还包括：

[0094] 步骤301，针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则

[0095] 。

[0096] 步骤302，针对每个雷达帧，设每个距离单元的噪声能量估计系数的初始值为 ，最大值为 ，搜索步长为 ，则第p个距离单元的第 次搜索的噪声能量估计系数为 。

[0097] 步骤303，以  为预设门限，设能量小于所述预设门限的散射中心为 ，其中， ， 。

[0098] 步骤304，根据所述散射中心，计算每个距离单元的噪声能量的估计值，其中设第p个距离单元的的第 次搜索的噪声能量估计值 ，则

[0099] 。

[0100] 图4是本发明提供的计算最佳噪声能量估计系数的流程示意图，如图4所示。上述步骤102中，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个
距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算
每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤还包括：

[0101] 步骤401，针对每个距离单元搜索噪声估计系数使得对应的噪声能量估计值与真值之间的差值最小。

[0102] 步骤402，设第p个距离单元的噪声能量估计系数为 ，则

[0103] 。

[0104] 图5是本发明提供的计算检测门限系数的流程示意图，如图5所示。上述步骤102中，所述针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单
元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距
离单元的噪声能量估计系数以及每个距离单元的检测门限系数的步骤还包括：

[0105] 步骤501，针对每个雷达帧，根据每个距离单元的噪声能量估计的真值计算该距离单元的检测门限系数。

[0106] 步骤502，设第p个距离单元的多普勒维散射中心能量由小到大排序的预设位置的散射中心能量作为第p个距离单元的噪声能量最大值的估计值 ，并且设第p个距离单元
的检测门限系数为 ，则

[0107] ；

[0108] 其中，表示固定比例系数。

[0109] 上述中，所述预设位置为约4/5的位置，比如多普勒维为512点时，可选取约4/5的位置即第410个散射中心。

[0110] 上述中，所述 根据实际经验值可选取 。

[0111] 上述步骤103中，所述针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测
门限系数的步骤包括：

[0112] 若所述至少一个雷达帧仅包括一个雷达帧，则将基于所述一个雷达帧计算出的每个距离单元的噪声能量估计系数作为该距离单元对应的目标噪声能量估计系数，并且将基
于所述一个雷达帧计算出的每个距离单元的检测门限系数作为该距离单元对应的目标检
测门限系数；

[0113] 若所述至少一个雷达帧包括多个雷达帧，则将不同雷达帧的同一距离单元的噪声能量估计系数的均值作为该距离单元对应的目标噪声能量估计系数，并且将不同雷达帧的
同一距离单元的检测门限系数的均值作为该距离单元对应的目标检测门限系数。

[0114] 当上述每个距离单元的噪声能量估计系数和检测门限系数计算完毕，并将每个距离单元的目标噪声能量估计系数和每个距离单元的目标检测门限系数进行存储，以便进入
第二部分的实际检测模式时被调用。

[0115] 经过以上步骤的参考检测门限值的自学习模式完成后，进入上述第二部分所述的参考检测门限值的实际检测模式，如下：

[0116] 图6是本发明提供的计算参考检测门限值的流程示意图，如图6所示。上述步骤104中，所述在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应
的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并
基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚
警检测的步骤包括：

[0117] 步骤601，针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则

[0118] 。

[0119] 步骤602，计算每个距离单元的噪声能量的估计值，其中，设第p个距离单元的目标噪声能量估计系数为 ，并且能量小于 的散射中心为 ， ，
。

[0120] 需要说明的是，该步骤中的第p个距离单元的噪声能量估计系数 为上述第一部分的自学习模式中的步骤401计算所得。

[0121] 则第p个距离单元的噪声能量的估计值 被表达为：

[0122] 。

[0123] 由于经过上述第一部分的在自学习模式后，第二部分的实际检测模式是不需要搜索的，所以第p个距离单元的噪声能量的估计值 的没有上标k。

[0124] 步骤603，计算每个距离单元的参考检测门限值，其中，第p个距离单元的参考检测门限值为 * ， 为第p个距离单元的目标检测门限系数。

[0125] 需要说明的是，该步骤中的第p个距离单元的目标检测门限系数 为上述第一部分的自学习模式中的步骤502计算所得。

[0126] 综上所述，本发明所述用于雷达探测的恒虚警检测方法，能有效解决不同场景下检测门限值的调整问题，通过运用自学习的方法实现在非均匀杂波背景下的高效检测，具
有较大的实用价值。

[0127] 以下通过一应用示例对本发明所述用于雷达探测的恒虚警检测方法进行说明。

[0128] 图7（a）图7（d）是本发明提供的应用示例的示意图，如图图7（a）图7（d）所示。~ ~

[0129] 图7（a）表示单帧雷达距离‑多普勒图像：距离512点，多普勒512点。图7（b）表示第241个距离单元对应的多普勒维散射中心图像及初始噪声能量估计系数对应的噪声能量初
始估计值。图7（c）表示噪声能量估计系数第5次迭代后得到的噪声能量估计值，从图中可以
看出，经过多次搜索之后，噪声能量的估计值已接近真值，即噪声能量估计的真值的线与噪
声能量5次搜索后的估计值的线重合。图7（d）表示检测门限值（即本发明所述的参考检测门
限值）。

[0130] 下面对本发明提供的用于雷达探测的恒虚警检测装置进行描述，下文描述的用于雷达探测的恒虚警检测装置与上文描述的用于雷达探测的恒虚警检测方法可相互对应参
照。

[0131] 图8是本发明提供的用于雷达探测的恒虚警检测装置的结构示意图，如图8所示。一种用于雷达探测的恒虚警检测装置800，包括矩阵生成模块810、第一计算模块820、第二
计算模块830、以及恒虚警检测模块840。其中，

[0132] 矩阵生成模块810，用于获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵。

[0133] 第一计算模块820，用于针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量
估计的真值，计算每个距离单元的噪声能量估计系数以及每个距离单元的检测门限系数。

[0134] 第二计算模块830，用于针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检
测门限系数。

[0135] 恒虚警检测模块840，用于在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单
元的参考检测门限值，并基于每个距离单元的参考检测门限值对所述雷达探测数据的多普
勒维散射中心进行恒虚警检测。

[0136] 可选地，所述矩阵生成模块810，还用于：

[0137] 将每个雷达帧的样本数据对应的每个距离单元的多普勒维的散射中心按照能量由小到大进行排序，生成该雷达帧对应的距离‑多普勒散射中心能量矩阵 ，表示为：

[0138] = ；

[0139] 其中，M表示多普勒维傅立叶变换的点数，N表示距离维傅立叶变换的点数；

[0140] 其中，第p个距离单元的散射中心能量排序后的向量 表示为：

[0141] 。

[0142] 可选地，所述第一计算模块820，还用于：

[0143] 针对每个雷达帧，以每个距离单元对应的散射中心能量排序中的中值作为该距离单元的噪声能量估计的真值；

[0144] 其中，设第p个距离单元的噪声能量估计的真值为 ，则 ， 。

[0145] 可选地，所述第一计算模块820，还用于：

[0146] 针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则

[0147] 。

[0148] 可选地，所述第一计算模块820，还用于：

[0149] 针对每个雷达帧，设每个距离单元的噪声能量估计系数的初始值为 ，最大值为 ，搜索步长为 ，则第p个距离单元的第 次搜索的噪声能量估计系数为 ；

[0150] 以  为预设门限，设能量小于所述预设门限的散射中心为 ，其中，， ；

[0151] 设第p个距离单元的的第 次搜索的噪声能量估计值 ，则

[0152] ;

[0153] 针对每个距离单元搜索噪声估计系数使得对应的噪声能量估计值与真值之间的差值最小，其中，设第p个距离单元的噪声能量估计系数为 ，则

[0154] 。

[0155] 可选地，所述第一计算模块820，还用于：

[0156] 针对每个雷达帧，根据每个距离单元的噪声能量估计的真值计算该距离单元的检测门限系数；

[0157] 其中，设第p个距离单元的多普勒维散射中心能量由小到大排序的预设位置的散射中心能量作为第p个距离单元的噪声能量最大值的估计值 ，并且设第p个距离单元的
检测门限系数为 ，则

[0158] ；

[0159] 其中，表示固定比例系数。

[0160] 可选地，所述第二计算模块830，还用于：

[0161] 若所述至少一个雷达帧仅包括一个雷达帧，则将基于所述一个雷达帧计算出的每个距离单元的噪声能量估计系数作为该距离单元对应的目标噪声能量估计系数，并且将基
于所述一个雷达帧计算出的每个距离单元的检测门限系数作为该距离单元对应的目标检
测门限系数；

[0162] 若所述至少一个雷达帧包括多个雷达帧，则将不同雷达帧的同一距离单元的噪声能量估计系数的均值作为该距离单元对应的目标噪声能量估计系数，并且将不同雷达帧的
同一距离单元的检测门限系数的均值作为该距离单元对应的目标检测门限系数。

[0163] 可选地，所述恒虚警检测模块840，还用于：

[0164] 针对每个雷达帧，计算每个距离单元对应的散射中心能量的均值，其中，设第p个距离单元的散射中心能量的均值 ，则

[0165] ；

[0166] 计算每个距离单元的噪声能量的估计值，其中，设第p个距离单元的目标噪声能量估计系数为 ，并且能量小于 的散射中心为 ， ， ，

[0167] 则第p个距离单元的噪声能量的估计值 被表达为：

[0168] ；

[0169] 计算每个距离单元的参考检测门限值，其中，第p个距离单元的参考检测门限值为* ， 为第p个距离单元的目标检测门限系数。

[0170] 图9示例了一种电子设备的实体结构示意图，如图9所示，该电子设备可以包括：处理器(processor)910、通信接口(Communications Interface)920、存储器(memory)930和
通信总线940，其中，处理器910，通信接口920，存储器930通过通信总线940完成相互间的通
信。处理器910可以调用存储器930中的逻辑指令，以执行所述用于雷达探测的恒虚警检测
方法，所述方法包括：

[0171] 获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵；

[0172] 针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距
离单元的噪声能量估计系数以及每个距离单元的检测门限系数；

[0173] 针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数；

[0174] 在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并
基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚
警检测。

[0175] 此外，上述的存储器930中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本
发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以
软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以
使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施
例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（ROM，
Read‑Only Memory）、随机存取存储器（RAM，Random Access Memory）、磁碟或者光盘等各种
可以存储程序代码的介质。

[0176] 另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序
指令被计算机执行时，计算机能够执行上述各方法所提供的所述用于雷达探测的恒虚警检
测方法，所述方法包括：

[0177] 获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵；

[0178] 针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距
离单元的噪声能量估计系数以及每个距离单元的检测门限系数；

[0179] 针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数；

[0180] 在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并
基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚
警检测。

[0181] 又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各提供的所述用于雷达探测的恒虚警检
测方法，所述方法包括：

[0182] 获取至少一个雷达帧的样本数据，并基于所述样本数据生成每个雷达帧对应的距离‑多普勒散射中心能量矩阵；

[0183] 针对每个雷达帧，基于其对应的距离‑多普勒散射中心能量矩阵，确定每个距离单元的噪声能量估计的真值，并根据所述每个距离单元的噪声能量估计的真值，计算每个距
离单元的噪声能量估计系数以及每个距离单元的检测门限系数；

[0184] 针对每个距离单元，基于所有计算出的该距离单元的噪声能量估计系数以及检测门限系数确定该距离单元对应的目标噪声能量估计系数和目标检测门限系数；

[0185] 在雷达对目标物体进行探测的过程中，基于雷达探测数据以及每个距离单元对应的目标噪声能量估计系数和目标检测门限系数，计算每个距离单元的参考检测门限值，并
基于每个距离单元的参考检测门限值对所述雷达探测数据的多普勒维散射中心进行恒虚
警检测。

[0186] 以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单
元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其
中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性
的劳动的情况下，即可以理解并实施。

[0187] 通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上
述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该
计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指
令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行各个实施
例或者实施例的某些部分所述的方法。

[0188] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可
以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；
而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和
范围。