一种基于背景预筛选的自适应CFAR目标检测方法

本发明属于雷达目标检测技术领域，涉及一种基于背景预筛选的自适应cfar目标检测方法，适用于在多目标、杂波边缘、均匀环境情况下自适应地检测目标。

背景技术：

雷达系统中常用恒虚警率((constantfalsealarmrate，cfar))方法对回波信号进行目标检测，目前常用的五种cfar检测器为：单元平均cfar检测器(ca-cfar)、最大选择cfar检测器(go-cfar)、最小选择cfar检测器(so-cfar)、有序统计量cfar检测器(os-cfar)、自适应cfar检测器，包含分布式cfar检测器(vi-cfar)、非均匀杂波估计cfar检测器(hce-cfar)以及其他的各种新型cfar。ca-cfar技术适用于局部均匀杂波中存在单个小目标的情况，在均匀杂波背景下取得了接近最优的检测性能，在杂波边缘和多目标情况下检测性能严重下降。针对这个问题，出现了os-cfar，适合大目标和多个目标情况，改善了多目标背景下的性能严重下降问题，但是在均匀杂波背景下的检测性能严重下降，并且在杂波边缘情况下的虚警尖峰控制能力差。针对杂波边缘情况，出现了go-cfar，go-cfar在杂波边缘环境有较好的检测性能，而在均匀杂波环境下，检测性能反而下降。so-cfar检测器是一种选小cfar检测器，在多目标环境中具有很好的检测性能。由于在杂波边缘环境中得到的检测阈值偏低，会比ca-cfar检测器带来更多的虚警。为了克服cfar的上述缺点，出现了自适应cfar检测器。vi-cfar是ca-cfar、os-cfar、go-cfar、so-cfar四个cfar检测器的组合或变体，可以根据杂波背景自适应地选择cfar检测器，但是在多目标环境下会呈现检测性能的下降。公开号为cn111337894a的专利申请“一种智能参考单元平均恒虚警检测方法”，将两种cfar检测的输出结果进行融合，并建立剔除信息矩阵提高目标检测效果，但该方法实现较为耗时。hce-cfar检测器是一种适用于非均匀背景下的自适应cfar方法，该方法首先找出不同背景下的临界单元，判断检测单元位于哪个背景区域，然后对该区域进行背景功率评估，最后进行cfar检测。该方法在杂波边缘检测较好，但如果参考窗分布均匀，hce-cfar算法同样只采用一部分参考窗进行背景功率估计，这样会损失一部分检测性能且增加算法复杂度去寻找杂波临界点。为了寻找一种可在多目标、杂波边缘、均匀环境下均具有良好检测性能的自适应cfar检测器。公开号为cn111929679a的专利申请“一种自适应加权截断统计恒虚警检测方法”，根据参考单元中的数据确定截断深度，比较参考单元数据与截断深度并进行迭代更新直至参考单元所有数据小于截断深度方可进行背景功率估计，由于涉及迭代计算，该方法效率不高。

因此，亟需一种能够在多目标、杂波边缘、均匀环境情况下自适应地检测目标的新方法。

技术实现要素：

有鉴于此，本发明的目的在于提供一种基于背景预筛选的自适应cfar目标检测方法，直接利用均值和方差进行潜在目标检测，利用背景功率均值代替潜在目标进行背景功率估计使得多目标情况下检测效率更高；另外，还在杂波边缘检测前增加了背景均匀性评估，提高hce-cfar在均匀背景下的检测效果；最终可在多目标、杂波边缘以及均匀环境下自适应地检测目标。

为达到上述目的，本发明提供如下技术方案：

一种基于背景预筛选的自适应cfar目标检测方法，首先对参考窗口进行参数估计剔除背景中的疑似目标，避免目标遮蔽效应；其次对背景均匀性进行评估，确保背景均匀时背景功率估计能使用到更多的参考窗口；对不均匀的背景检测杂波边缘，确定只保留属于本检测单元所在背景的参考窗口进行背景功率估计，提高算法在边缘环境的检测效果，避免由于背景功率估计的不准确导致边缘处出现漏警或虚警；最后再用得到的参考窗口利用ca-cfar对检测单元进行判决。具体包括以下步骤：

s1：接收雷达回波信号，并对信号进行平方律检波处理；

s2：设置检测单元左右参考窗的大小，将平方律检波后的数据载入滑窗；

s3：通过下式计算参考窗口的均值μb及方差σb²，检测出潜在目标点用均值μb替代；

s4：通过最大似然法计算似然相关函数值和l0，估计背景参考窗口的均匀性，若背景均匀，则用原参考窗口作为背景进行步骤s6；若背景不均匀，则进行步骤s5；其中l0为背景均匀时的似然相关函数值，为杂波边缘临界单元所处位置的似然相关函数值；

s5：利用步骤s4中的计算杂波临界值，将目标一侧的参考窗口选出进行步骤s6；

s6：进行背景功率估计；

s7：将背景功率估计值与标称因子相乘得到检测门限；

s8：将检测单元数据与检测门限进行比较从而做出判决。

进一步，步骤s1中，通过匹配滤波器获取雷达回波信号各个距离单元的i、q两路数据，将i、q两路数据进行平方求和即平方律检波处理，此时可将处理后的数据看作符合瑞利分布模型进行后续的检测操作。

进一步，步骤s2中，设置检测单元左右参考窗的大小，设置左右参考窗长度相等，为n/2，设置检测单元两侧的若干个距离单元b1···bn为背景窗口，将步骤一中经过平方律检波处理后的数据载入检测单元背景窗口。

进一步，步骤s3中，均值μb及方差σb²的计算公式为：

其中，bi表示参考窗口中的第i个参考单元的数据，n表示检测单元s的参考窗口中共有n个参考单元，μb和σb²分别为检测单元s的参考窗口的均值及方差。

更进一步，步骤s3中，对潜在目标的检测用以下方式判定：

bi-μb＞k·σb²,i＝1,...,n

其中，k为常数。

进一步，步骤s4中，设m为两种背景分布的临界单元，则m-1位于背景区域1，m+1位于背景区域2。设检测单元s所处位置为s，当s＜m用背景区域1的功率估计λ1²进行检测，参考单元数为m。当s＞m用背景区域2的功率估计λ2²进行检测，参考单元数为n-m。

背景均匀时的似然相关函数为：

其中，为经过步骤s3处理后的参考窗口的功率，表示为：其中xi为当前参考窗口中的第i个参考单元的数据，n为剔除疑似目标后参考单元的数量；

背景均匀性判断具体为：

若则杂波背景具有一致性；若则背景杂波不均匀。

更进一步，步骤s4中，杂波边缘临界单元所处位置的似然相关函数为：

其中，和分别为杂波临界单元位置分割后背景区域1和背景区域2的背景功率，l0为背景均匀时的似然相关函数值，为经过步骤s3处理后的参考窗口的功率。

进一步，步骤s5中，杂波临界值的计算公式为：

其中，是对对背景区域1和背景区域2分界临界点的估计值；估计值求解后，比较检测目标所处位置s与的大小关系。如果则目标处于背景区域1，使用杂波功率估计λ1²计算检测门限xt。如果则目标处于背景区域2，使用功率估计λ2²计算检测门限xt。

进一步，步骤s6中，计算经过预筛选后剩余的数据的和值，具体过程为：

其中，xi(i＝1...j)为第i个参考单元数据，j为预筛选后参考窗口中检测单元的数目，z为参考单元的功率估计值，背景均匀时可直接使用步骤四中已经计算出的作为背景功率。

进一步，步骤s7中，自适应检测门限为：

xt＝t*z

其中，z为经过预筛选后剩余的数据的和值；t为门限因子，在杂波符合高斯分布的假设下为：

t＝pfa^-1/2j-1

其中，pfa的值为人为设置的虚警率，通常取10^-3。

进一步，步骤s8中，将检测单元数据与检测门限进行比较：

其中，h1和h0分别表示有目标和无目标。

本发明的有益效果在于：

(1)本发明方法进行了潜在目标剔除，可以更好地在多目标的情况下对目标进行检测，避免目标遮蔽效应。

(2)本发明方法对背景的均匀性进行评估，可以保证在均匀环境下得到更准确的背景功率估计。

(3)本发明方法对杂波临界单元进行检测，可以使得目标检测方法在杂波环境有更好的检测效果。

本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。

附图说明

为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作优选的详细描述，其中：

图1为基于背景预筛选的自适应cfar目标检测方法系统框图；

图2为各方法在单目标环境中的检测结果对比图及局部放大图；

图3为各方法在多目标环境中的检测结果对比图及局部放大图；

图4为各方法在杂波边缘环境中的检测结果对比图及局部放大图；

图5为各方法在复杂环境中的检测结果对比图及局部放大图。

具体实施方式

以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。

请参阅图1～图5，图1为本发明提供的一种基于背景预筛选的自适应cfar目标检测方法，通过多次对背景进行预筛选，使得目标检测方法在多目标以及不均匀情况下具有自适应性，并尽量保证检测的效率。具体步骤为：

步骤一：接收雷达回波信号，并对信号进行平方律检波处理；

通过匹配滤波器获取雷达回波信号各个距离单元的i、q两路数据，将i、q两路数据进行平方求和即平方律检波处理，将处理后的数据看作符合瑞利分布模型进行后续的检测操作。

步骤二：设置检测单元左右参考窗的大小，将平方律检波后的数据载入滑窗；

设置检测单元左右参考窗的大小，左右均n＝10个参考单元，假设检测单元两侧的若干个距离单元b1···bn为背景窗口，将步骤一中经过平方律检波处理后的数据载入检测单元背景窗口。

步骤三：计算参考窗口的均值μb及方差σb²，通过比较剔除背景中的潜在目标，检测出的潜在目标点的数据用均值μb替代；

其中均值μb及方差σb²的计算方式为：

步骤三中对潜在目标的检测用以下方式判定：

bi-μb＞k·σb²,i＝1,...,n.(3)

本实施例中k取4。

步骤四：通过最大似然法计算似然相关函数值和l0，估计背景参考窗口的均匀性，若背景均匀，则用原参考窗口作为背景进行步骤六，若背景不均匀，则进行步骤五；

设m为两种背景分布的临界单元，则m-1位于背景区域1，m+1位于背景区域2。设检测单元s所处位置为s，当s＜m用背景区域1的功率估计λ1²进行检测，参考单元数为m。当s＞m用背景区域2的功率估计λ2²进行检测，参考单元数为n-m。背景均匀和杂波临界点的似然相关函数分别为：

其中，分别为：

其中xi为当前参考窗口中的第i个参考单元的数据，为杂波边缘临界单元所处位置的似然相关函数值，和分别为杂波临界单元位置分割后背景区域1和背景区域2的背景功率，l0为背景均匀时的似然相关函数值，为经过步骤三处理后的参考窗口的功率。

将l0与进行对比，若则杂波背景具有一致性，则用原参考窗口作为背景进行步骤六。若则背景杂波不均匀，则进行步骤五。

步骤五：利用步骤四中的计算杂波临界值，将目标一侧的参考窗口选出进行步骤六；

由于步骤四中因此背景区域1和背景区域2的临界点可以表示为：

由此可知，是对背景区域1和背景区域2分界临界点的估计值。估计值求解后，比较检测目标所处位置s与位置的大小关系。如果则目标处于背景区域1，使用杂波功率估计λ1²计算检测门限xt。如果则目标处于背景区域2，使用功率估计λ2²计算检测门限xt。

步骤六：进行背景功率估计；

计算经过预筛选后剩余的数据的和值，具体过程为：

其中xi(i＝1...j)为第i个参考单元数据，j为预筛选后参考窗口中检测单元的数目，z为参考单元的功率估计值，背景均匀时可直接使用步骤四中已经计算出的作为背景功率。

步骤七：将背景功率与标称因子相乘得到检测门限；

自适应检测门限为:

xt＝t*z(11)

其中t为门限因子，在杂波符合高斯分布的假设下为：

t＝pfa^-1/2j-1(12)

步骤八：将检测单元数据与检测门限进行比较从而做出判决。

其中h1和h0分别表示该检测单元有目标或无目标。

通过以下仿真实验对本发明效果作进一步验证说明。

(一)实验条件

实验运行平台：matlabr2016a、intelcorei5-7400处理器、8.00gb内存、3.00ghz主频。

(二)仿真实验内容

分别在单目标环境、多目标环境、杂波边缘环境、同时存在多目标及杂波边缘的复杂环境下进行目标检测，对比方法为ca-cfar、so-cfar、os-cfar。三个对比方法的参数设置均为：检测单元前后沿各有3个保护单元和10个参考单元。本方法(在仿真图中用pre-cfar表示)前后窗各设置10个参考单元，不设置保护单元，参数k取4。

单目标环境下的对比实验中，snr＝20db，目标处于距离单元50处，信噪比为35db。从仿真对比图2中可看出单目标环境下各方法均能正确检测出目标，且不会产生虚警。

多目标环境下的对比实验中，snr＝20db，目标处于距离单元50处及距离单元58处，信噪比分别为35db和40db。从仿真对比图3中可看出多目标环境下ca-cfar不能同时把两个目标检测出来，本方法和另外两个方法均能检测出目标，但是so-cfar在某些杂波处产生了虚警。

杂波边缘环境下的对比实验中，前100个距离单元snr＝20db，后100个距离单元snr＝30db，目标处于距离单元95处，信噪比为40db。从仿真对比图4中可看出杂波边缘环境下ca-cfar不能同时把两个目标检测出来，本方法和so-cfar能够准确地在边缘将目标检测出来，但是so-cfar依旧会在某些杂波处产生虚警。

复杂环境下的对比实验中，前100个距离单元snr＝20db，后100个距离单元snr＝30db，目标分别处于距离单元50、58、95处，信噪比分别为35db、40db、40db。从仿真对比图5中可看出复杂环境下ca-cfar未在多目标背景下检测出目标，而so-cfar和os-cfar虽然准确检测处理紧挨着的两个目标以及杂波边缘的目标，但是同时在杂波处产生了虚警。

综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。

最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。