一种雷达恒虚警检测门限自适应设置方法与流程

本发明属于雷达目标检测技术领域，涉及雷达恒虚警(CFAR)检测新方法；通过该恒虚警检测门限自适应设置方法，实现了雷达恒虚警检测门限的精细调整，减小了抑制局部强杂波带来的检测损失。

背景技术：

通常，雷达目标检测为了抑制虚警，大多采用恒虚警检测方法，该方法按照一定准则，根据噪声和杂波的强弱得到环境估值，并对该估值乘以一个系数作为最终的检测门限。该检测方法可根据噪声和杂波强弱来调整检测门限，达到虚警概率基本恒定的效果。

恒虚警检测器原理框图如图1所示，是统一模式下的恒虚警检测器的原理框图，信号经过平方率检波或者线性检波后按一定规则产生检测门限，不同类型的恒虚警检测器参考单元的选择和门限的产生方式不同，y为待检单元，两侧为保护单元，防止目标本身对门限值的影响，Z为根据临近检测单元信号幅度得到的环境估值，c为门限系数，用来调节门限大小。

最早的均值类恒虚警检测是单元平均恒虚警(CA-CFAR)检测，但是为了降低杂波非均匀分布对检测性能的影响，后来又出现了最大选择恒虚警(GO-CFAR)检测、最小选择恒虚警(SO-CFAR)检测和加权单元平均恒虚警(WCA-CFAR)检测。但是上述恒虚警检测方式在均匀背景下检测效果比较理想，对于局部存在强杂波干扰的环境，为了抑制局部强杂波带来的高虚警，需要将门限系数α提高，以此增大检测门限，降低虚警。但是该处理方式导致检测概率的降低，对于杂波较弱的区域来说，带来了不必要的检测损失。

为此本发明提出一种自适应的CFAR检测门限设置方法，该方法可对不同区域设置不同的检测门限系数，从而在强杂波区域设置较高的检测门限，已达到降低虚警的目的，同时在弱杂波区域设置较低的检测门限，以提高该区域的检测概率。从雷达整个检测区域考虑，雷达的检测性能得到较大改善。

技术实现要素：

要解决的技术问题

为了避免现有技术的不足之处，本发明提出一种雷达恒虚警检测门限自适应设置方法。

技术方案

一种雷达恒虚警检测门限自适应设置方法，其特征在于步骤如下：

步骤1：将雷达探测区域划分为若干个子区域；

步骤2：对子区域进行量化，量化时可采用计算子区域的水平投影面积或者空间体积；

步骤3：计算各子区域内的点迹密度L，计算时使用子区域点数除以该子区域对应的量化值；

步骤4：设置点迹密度门限上限M，设置点迹密度门限下限N，其中M≥N；所述点迹密度门限上限为根据雷达虚警概率计算出最大点迹数量，然后用该值除以整个检测区域的量化值；所述的点迹密度门限下限为降低雷达虚警概率，根据降低后的虚警概率计算出最大点迹数量，然后用该值除以整个检测区域的量化值；

步骤5：将各子区域内的点迹密度L与点迹密度门限值进行比较：当L≥M，则增大该区域内恒虚警检测门限系数c，提高恒虚警检测门限；当L≤N，则减小该区域内恒虚警检测门限系数c；

步骤6：设置最大恒虚警门限系数值c_max，当恒虚警门限系数值大于等于c_max时，不再增加；设置最小恒虚警门限系数值c_min，当恒虚警门限系数值小于等于c_min时，不再减小。

步骤1中按照方位、俯仰、距离三个维度根据需要任意组合进行划分。

步骤4中可设置多个点迹密度门限上限或点迹密度门限下限。

有益效果

本发明与现有技术相比具有以下优点：

1、对CFAR检测门限进行精细调节；

2、降低了强杂波区域的虚警概率；

3、提高弱杂波区域的检测概率。

附图说明

图1恒虚警检测器原理框图

图2雷达自适应CFAR门限设置方法处理流程图

图3雷达自适应CFAR门限设置方法算法流程图

具体实施方式

现结合实施例、附图对本发明作进一步描述：

本发明涉及的雷达恒虚警检测门限自适应设置方法为：

1)将雷达探测区域按照一定规则划分为若干子区域；

2)按照一定大小(一定面积或者一定体积)对雷达探测空间区域(水平投影面积或者空间体积)进行量化。

3)统计各子区域内的点迹(信号处理机检测到的过门限信号)密度L，该密度计算时使用子区域点数除以该子区域对应的量化值；

4)设置点迹密度门限上限M和下限N，满足M≥N，根据雷达虚警概率计算出最大点迹数量，然后用该值除以整个检测区域的量化值，其结果作为门限上限，门限上限大于0。将雷达虚警概率降低一定程度，根据降低后的虚警概率计算出最大点迹数量，然后用该值除以整个检测区域的量化值，其结果作为门限下限，门限下限大于等于0；

5)将各子区域内的点迹密度L与点迹密度门限值进行比较，当L≥M，则增大该区域内恒虚警检测门限系数c，提高恒虚警检测门限，当L≤N，则减小该区域内恒虚警检测门限系数c；

上述操作在雷达工作过程中同步进行，根据当前扫描周期内点迹密度与点迹密度门限值的对比结果，调整下一个扫描周期内的恒虚警检测门限系数。

对于上述第一步，划分子区域包括但不限于在方位、俯仰、距离三个维度根据需要任意组合进行划分，例如，按照方位角度划分，每10°划分为一个子区域，若雷达探测为360°扫描，则将检测区域划分为36个子区域；在此基础之上再按照距离进行划分，将雷达探测量程按照“远中近”等分为三段，则整个检测区域被划分为108个子区域；

对于上述第二步，对检测区域进行量化，包括但不限于根据每个子区域的投影面积、空间体积进行归一；

对于上述第四步，包括但不限于设置一或多个上限值和下限值；可将虚警概率提高不同程度，利用提高后的虚警概率计算其它上限值，可将虚警概率适当降低不同程度，利用降低后的虚警概率计算其它下限值；

对于上述第五步，恒虚警门限系数增大条件可设置为点迹密度大于等于上限值或者点迹密度大于上限值；恒虚警门限系数减小条件可设置为点迹密度小于等于下限值或者点迹密度小于下限值；若设置多个上限值或者下限值，则可根据不同的门限值调节恒虚警门限系数增大或者较小的速度。

以某对空情报雷达为例，雷达探测范围：俯仰0～60°，方位0～360°，距离180km；

a)将雷达探测区域按照方位角每10°等间隔划分为36个子区域；将雷达探测量程按照“远、中、近”等分为3段，每段距离60km；俯仰不做划分，则整个检测区域被划分为108个子区域；

b)计算各子区域的水平面投影面积，将0～60km距离段的一个子区域水平投影面积作为单位面积，量化为1，则60～120km距离段上每个子区域的水平投影面积量化为为3；120～180km距离段上每个子区域的水平投影面积量化为5，整个检测区域对应的量化值为324；

c)按公式计算每个子区域内的点迹密度：

子区域点迹密度＝子区域内点迹总数/子区域水平面投影面积量化；

d)假设根据虚警概率计算出的最大点迹数量为600，设置第一点迹密度上限M1＝600/324，M1取2，将虚警概率提高3倍，则对应的最大点迹数量为1800，设置第二点迹密度上限M2为6，将虚警概率降低为原来的1/2，则对应的最大点迹数量为300，设置点迹密度下限点击密度下限N为1，满足M2≥M1≥N；

e)若当前扫描周期内某一子区域点迹密度大于M2，则提高恒虚警门限系数值加a，对应恒虚警检测门限提高1dB；单位面积内点迹数量大于M1小于M2，则将恒虚警门限系数值加b，对应恒虚警检测门限提高0.5dB；若单位面积内点迹数量小于N，则将恒虚警门限系数值减c，对应恒虚警检测门限降低0.5dB；否则，恒虚警门限系数值不变；

f)设置最大恒虚警门限系数值c_max，该值设置为纯噪声背景下，系统可接受的子区域最小检测概率对应的恒虚警检测门限系数，当恒虚警门限系数值大于等于c_max时，不再增加，该条件优先于e项所列条件；设置最小恒虚警门限系数值c_min，该值设置为纯噪声背景下，系统可接受最大虚警概率对应的检测门限系数。当恒虚警门限系数值小于等于c_min时，不再减小，该条件优先于e项所列条件。