CN111999715B - 异质杂波下目标知识辅助自适应融合检测方法 - Google Patents

Abstract

本发明公开了一种异质杂波下目标知识辅助自适应融合检测方法，属于宽带雷达信号处理领域。针对实际异质杂波环境中不同距离单元间杂波协方差矩阵结构非均匀特点，构建易于数学处理的合理先验分布，将杂波先验信息与辅助数据中所含异质杂波信息进行高效融合，构建易于数学处理的异质杂波协方差矩阵结构概率密度函数的几何平均估计方法，通过考虑距离扩展目标不同散射点之间是否统计独立的具体情况，分别求解距离扩展目标未知复幅度的最大似然估计表达式，构建具有闭合形式的距离扩展目标知识辅助自适应融合检测器，提高了宽带雷达对异质杂波环境的适应能力，提升了复杂电磁环境下宽带雷达对弱小目标的探测能力，具有推广应用价值。

Description

异质杂波下目标知识辅助自适应融合检测方法

技术领域

本发明隶属于宽带雷达信号处理领域，具体涉及一种异质杂波下目标知识辅助自适应融合检测方法。

背景技术

不同于低分辨率窄带雷达，宽带雷达带宽大、距离分辨率高，在抗干扰、反侦察、精确探测及成像、高精度跟踪、目标识别等方面具有明显的优势，在现代军事和民用领域引起了广泛重视，已经成为现代雷达发展的一个重要方向。窄带雷达的距离分辨单元一般远大于常见目标的几何尺寸，目标回波信号只占据一个距离分辨单元，实际目标常被作为“点目标”进行处理。而宽带雷达目标的回波信号不仅仅只占据一个距离分辨单元，而是分布在不同的径向距离分辨单元中，呈现为“一维距离像”，形成“距离扩展目标”。随着宽带雷达的广泛应用，距离扩展目标检测问题正受到越来越多的关注，成为近年来雷达信号处理界的热点和难点问题之一。

一方面，由于宽带雷达观测到的目标回波分布在多个径向距离单元上，若仍采用窄带雷达的点目标检测方法，针对单个距离单元对回波信号进行目标检测，并利用邻近的距离单元采样进行背景杂波统计特性估计，则距离扩展目标强散射点的能量会泄漏到邻近距离单元中形成“信号污染”现象，并进一步对单个待检测距离单元的目标信号构成遮蔽效应，导致无法检测到目标，由于未能充分利用分布在多个距离单元中目标回波的全部能量，这不但不能体现宽带雷达的优势，反而会降低对距离扩展目标的检测能力。另一方面，宽带雷达目标自适应检测面临目标自然环境复杂多变、电磁干扰等多种因素的影响，杂波统计特性不再满足独立同分布的均匀环境假设条件，这一问题在宽带雷达目标检测场景下尤为突出，导致现有距离扩展目标检测方法难以取得理想的检测效果。现有杂波模型主要包括均匀环境、部分均匀环境、非高斯杂波环境、异质杂波环境等。其中，部分均匀环境下距离扩展目标检测模型中，假设多个待检测距离单元的数据(也称主数据)与只含纯杂波的参考距离单元数据(也称辅助数据)中的杂波分量具有相同的协方差矩阵结构但功率水平不同；而以复合高斯分布建模的非高斯杂波环境中，假设主数据和辅助数据中的杂波分量具有相同的协方差矩阵结构但不同距离单元间杂波功率水平不同。前述两种杂波模型中，均假定杂波分量具有相同的协方差矩阵结构，但在宽带雷达探测面临的复杂海杂波、地物杂波等实际环境下，由于各种干扰因素的存在，不同距离单元间杂波协方差矩阵结构的均匀性会被进一步破坏，前述均匀、部分均匀和非高斯杂波模型中杂波分量具有相同协方差矩阵结构的假设条件将无法满足。此时需采用异质杂波模型进行环境建模，异质杂波环境下不同距离单元间杂波协方差矩阵结构相似但不相同，需通过合理的先验分布将辅助数据和主数据中的杂波协方差矩阵结构联系起来。

在宽带雷达探测场景下，针对实际异质杂波环境中不同距离单元间杂波协方差矩阵结构非均匀特点，如何利用合理先验分布，基于辅助数据中所含异质杂波信息来准确估计主数据中多个距离单元中的杂波统计信息，在满足估计准确性的条件下构建便于数学处理的估计方法，进而构建具有闭合形式的距离扩展目标检测统计量，在增强检测器恒虚警率(CFAR)特性的前提下，进一步提高宽带雷达对异质杂波环境的适应能力，改善目标检测性能，是提升复杂电磁环境下宽带雷达探测能力的关键，也是目前急需解决的难题之一。

发明内容

在宽带雷达探测场景下，针对实际异质杂波环境中不同距离单元间杂波协方差矩阵结构非均匀特点，如何构建易于数学处理的合理先验分布，将杂波先验信息与辅助数据中所含异质杂波信息进行高效融合，构建易于数学处理的协方差矩阵结构估计方法，在满足估计准确性的条件下为宽带雷达距离扩展目标检测器CFAR设计提供扎实基础，进而构建具有闭合形式的距离扩展目标知识辅助自适应融合检测器，通过考虑距离扩展目标不同散射点之间是否统计独立的具体情况，进一步提高宽带雷达对异质杂波环境的适应能力，改善异质杂波下的宽带雷达目标检测性能，提升复杂电磁环境下宽带雷达对弱小目标的探测能力。

本发明所述异质杂波下目标知识辅助自适应融合检测方法包括以下技术措施：

步骤1从P个待检测距离单元获取主数据X，从与待检测距离单元临近的K个纯杂波距离单元获取K个辅助数据y_k(k＝1,2,…K)，利用逆Wishart分布对杂波协方差矩阵结构M进行先验分布建模，进而基于K个辅助数据分别构建M的K个后验概率密度函数，将K个后验概率密度函数进行几何平均，构建异质杂波协方差矩阵结构概率密度函数的几何平均估计方法，获得异质杂波协方差矩阵结构概率密度函数的几何平均估计

具体步骤包括：

由P个待检测距离单元回波复幅值构成P个待检测数据向量x_m(m＝1,2,…,P)，进而构建N×P维的主数据X＝[x₁,x₂,…,x_P]，即x_m表示第m个主数据分量，P为大于1的自然数；以P个待检测距离单元为中心，在其前后分别连续取一定数量的不包含目标的距离单元回波复幅值，构成K个只含纯杂波的辅助数据y_k(k＝1,2,...K)，其中，x_m和y_k均为N×1维的复向量，N表示雷达接收阵元数与相干处理脉冲数的乘积；

主数据和辅助数据中的杂波分量服从条件复高斯分布，但其杂波协方差矩阵结构M为随机矩阵，且服从自由度为L、均值为LR₀的逆Wishart分布，其中，M和R₀均为N×N维的Hermitian复数矩阵，M的先验分布概率密度函数(PDF)可表示为

其中，exp[·]表示指数函数，函数tr(·)和det(·)分别表示对矩阵求迹和求行列式，函数I[L,(LR₀)^-1]定义为：

上式中，Γ(·)表示Gamma函数。

令N×N维的复数矩阵R_k(k＝1,2,…K)表示为

其中，上标“H”表示共轭转置。

给定第k个辅助数据y_k时矩阵M的条件PDF(即以辅助数据y_k为条件的M的后验PDF)可以表示为：

其中，f(y_kM)表示已知M时辅助数据y_k的条件复高斯分布PDF，f(y_k)表示y_k的PDF，可由下式计算

结合杂波先验分布信息和K个辅助数据y_k(k＝1,2,...K)进行杂波协方差矩阵结构估计，鉴于前述多个函数和后续检测统计量设计中多涉及乘除运算，基于K个条件PDF f(My_k)(k＝1,2,…K)，采用几何平均方法估计出矩阵M的PDF为：

上式为M的PDF几何平均估计表达式，可见其具有闭型形式，为后续设计闭合形式的距离扩展目标知识辅助自适应融合检测统计量奠定了基础。

步骤2基于异质杂波协方差矩阵结构概率密度函数的几何平均估计

在有目标和无目标两种假设下求解主数据X的概率密度函数，根据广义似然比检验准则，分别求解距离扩展目标P个散射点未知复幅度的最大似然估计，针对距离扩展目标不同散射点之间是否统计独立的两种情况，分别构建异质杂波下距离扩展目标知识辅助自适应融合检测的检测统计量λ；具体步骤包括：

情况1：不考虑距离扩展目标不同散射点间的独立性，即针对距离扩展目标不同散射点之间非统计独立的情况；

在无目标的H₀假设下，基于矩阵M的PDF的几何平均估计

主数据X的PDF可表示为

在有目标的H₁假设下，基于矩阵M的PDF的几何平均估计

主数据X的PDF可表示为

其中，上述两式中，f_X(X|M,H₀)和f_X(X|M,H₁)分别表示H₀和H₁两种假设下主数据X在M已知时的复高斯分布条件PDF；1×P维的复幅度向量a表示距离扩展目标的未知复幅度，具体包含P个散射点的未知复幅度；p表示已知的空时导向矢量，是一个N×1维的单位向量，根据雷达系统工作参数确定；

根据广义似然比检验准则，距离扩展目标检测统计量可初步表示为

令N×N维的复矩阵C为：

复幅度向量a的最大似然估计可表示为

将式(11)的

代入式(9)中替换其中的未知量a，并对式(9)中Λ₁(X)的表达式进行等价数学变换，可得异质杂波下距离扩展目标知识辅助自适应融合检测统计量λ₁，即

其中，I_N表示N×N维的单位矩阵；

在不考虑距离扩展目标不同散射点间独立性的情况下，可令

λ＝λ₁ (13)

情况2：针对距离扩展目标不同散射点间统计独立的情况；

在无目标的H₀假设下，基于矩阵M的PDF的几何平均估计

主数据X的第m个分量x_m的PDF可表示为

在有目标的H₁假设下，基于矩阵M的PDF的几何平均估计

主数据X的第m个分量x_m的PDF可表示为

其中，上述两式中，

和

分别表示H₀和H₁两种假设下x_m在M已知时的复高斯分布条件PDF；1×P维的复幅度向量a可表示为a＝[a₁,a₂,…,a_P]，即包含了P个散射点的未知复幅度a_m(m＝1,2,…,P)，其中，复标量a_m(m＝1,2,…,P)表示距离扩展目标第m个散射点的未知复幅度；

根据广义似然比检验准则，距离扩展目标检测统计量可表示为

距离扩展目标第m个散射点的未知复幅度a_m(m＝1,2,…,P)的最大似然估计可表示为

将式(17)的

代入式(16)中替换其中的未知量a_m(m＝1,2,…,P)，并对式(16)中Λ₂(X)的表达式进行等价数学变换，可得散射点独立时异质杂波下距离扩展目标知识辅助自适应融合检测统计量λ₂，即

上式中，|·|表示取复数的模值。

针对距离扩展目标不同散射点间统计独立的情况，可令

λ＝λ₂ (19)

步骤3为保持检测方法的CFAR特性，根据预设的虚警概率设置检测门限T；将主数据对应的检测统计量λ与检测门限T进行比较，若λ≥T，则判定当前待检测距离单元存在距离扩展目标，主数据不作为后续其他待检测距离单元的辅助数据；反之若λ＜T，则判定当前待检测距离单元不存在距离扩展目标，主数据作为后续其他待检测距离单元的辅助数据。

与背景技术相比，本发明的有益效果是：1)针对实际异质杂波环境中不同距离单元间杂波协方差矩阵结构非均匀特点，构建了易于数学处理的杂波协方差矩阵结构合理先验分布；2)将杂波先验信息与辅助数据中所含异质杂波信息进行高效的几何平均融合，构建了易于数学处理的协方差矩阵结构估计方法，在满足估计准确性的条件下为距离扩展目标检测器CFAR设计奠定扎实基础；3)针对距离扩展目标不同散射点之间独立与否，分别求解距离扩展目标未知复幅度的最大似然估计表达式，利用复幅度简洁有效的表达式，分别构建了异质杂波下距离扩展目标知识辅助自适应融合检测器，均具有闭合形式的表达式，且在保持CFAR特性的同时降低了检测统计量的计算复杂度，提高了宽带雷达对异质杂波环境的适应能力，改善了异质杂波下的目标检测性能，提升了复杂电磁环境下宽带雷达对弱小目标的探测能力；4)本发明方法适用于部分非宽带雷达探测情形，例如，使用低/中分辨率雷达检测大目标或检测以相同速度运动的空间邻近点目标群(舰艇编队、飞机编队、车辆编队等情况)，具有很好的应用前景。

附图说明

图1是本发明所提出的异质杂波下目标知识辅助自适应融合检测方法的功能模块图。图1中，1.中间矩阵计算模块，2.杂波协方差矩阵结构概率密度函数的几何平均估计模块，3.距离扩展目标知识辅助自适应融合检测器构建模块，4.检测判决模块。

具体实施方式

下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明做出的任何修改和改变，都落入本发明的保护范围。

为验证本发明所述方法的有效性，本具体实施方式给出了二个实施例，第一个实施例针对宽带雷达对地探测环境，第二个实施例针对宽带雷达对海探测环境。

实施例1：

参照说明书附图1，实施例1的具体实施方式分为以下几个步骤：

步骤A1利用对地探测宽带雷达，对待检测地域周围的无目标范围进行雷达照射，获得与待检测距离单元邻近的不包含目标的距离单元回波复幅值，构成K个只含纯地杂波的辅助数据y_k(k＝1,2,…K)，将辅助数据送至中间矩阵计算模块(1)；在中间矩阵计算模块(1)中，依据式(3)计算矩阵R_k(k＝1,2,…K)，并将矩阵R_k(k＝1,2,...K)送至杂波协方差矩阵结构概率密度函数的几何平均估计模块(2)；在杂波协方差矩阵结构概率密度函数的几何平均估计模块(2)中，依据式(6)计算杂波协方差矩阵结构M概率密度函数的几何平均估计

并将

送至距离扩展目标知识辅助自适应融合检测器构建模块(3)；

值得注意的是，步骤A1中，利用逆Wishart分布对杂波协方差矩阵结构进行先验分布建模，充分考虑了实际建筑物、林地、草地等不同地杂波环境导致的异质杂波情况，针对异质杂波下不同距离单元间杂波协方差矩阵结构非均匀特点，构建的地杂波协方差矩阵结构先验分布易于数学处理；另外，通过将杂波先验信息与辅助数据中所含异质杂波信息进行高效的几何平均融合，所获得的协方差矩阵结构估计表达式具有闭型形式，便于后续数学处理，在满足估计准确性的条件下，为对地探测宽带雷达距离扩展目标检测器CFAR设计奠定扎实基础。

步骤A2在距离扩展目标知识辅助自适应融合检测器构建模块(3)中，依据式(12)和式(13)计算异质地杂波下距离扩展目标知识辅助自适应融合检测的检测统计量λ，并将λ送至检测判决模块(4)；

值得注意的是，针对地面目标不同散射点间相关性较强的情况，本实施例中考虑距离扩展目标不同散射点间的非独立性，进而采用复幅度向量a的最大似然估计表达式(11)进行检测器设计，可充分兼顾不同散射点间相关性的影响，提高特定场景下的检测性能；另外，在步骤A2中，所构建的复杂地面环境异质杂波下距离扩展目标知识辅助自适应融合检测器，在保持CFAR特性的同时，可进一步提高宽带对地雷达在异质地杂波环境下的适应能力，改善复杂地面环境异质杂波下的距离扩展目标检测性能，提升复杂电磁环境下宽带雷达对弱小地面目标的探测能力。本实施例方法同样适用于使用低/中分辨率雷达检测以相同速度运动的空间邻近车辆点目标群，具有很好的应用前景。

步骤A3根据预设的虚警概率设置检测门限T：具体为，设定虚警概率为P_fa，根据蒙特卡洛方法，依据前期积累的100/P_fa个实测地杂波数据计算检测阈值T。进一步，检测统计量λ与检测门限T进行比较，若λ≥T，则判定当前待检测距离单元存在距离扩展目标，主数据X不作为后续其他待检测距离单元的辅助数据；反之若λ＜T，则判定当前待检测距离单元不存在距离扩展目标，主数据X作为后续其他待检测距离单元的辅助数据。

实施例2：

参照说明书附图1，实施例2的具体实施方式分为以下几个步骤：

步骤B1利用对海探测宽带雷达，对待检测海域周围的无目标范围进行雷达照射，获得与待检测距离单元邻近的不包含目标的距离单元回波复幅值，构成K个只含纯海杂波的辅助数据y_k(k＝1,2,...K)，将辅助数据送至中间矩阵计算模块(1)；在中间矩阵计算模块(1)中，依据式(3)计算矩阵R_k(k＝1,2,...K)，并将矩阵R_k(k＝1,2,...K)送至杂波协方差矩阵结构概率密度函数的几何平均估计模块(2)；在杂波协方差矩阵结构概率密度函数的几何平均估计模块(2)中，依据式(6)计算杂波协方差矩阵结构M概率密度函数的几何平均估计

并将

送至距离扩展目标知识辅助自适应融合检测器构建模块(3)；

值得注意的是，步骤B1中，利用逆Wishart分布对杂波协方差矩阵结构进行先验分布建模，充分考虑了不同海况海洋环境中杂波尖峰等导致的异质杂波情况，针对海洋环境异质杂波下不同距离单元间杂波协方差矩阵结构非均匀特点，构建的海杂波协方差矩阵结构先验分布易于数学处理；另外，通过将杂波先验信息与辅助数据中所含异质杂波信息进行高效的几何平均融合，所获得的协方差矩阵结构估计表达式具有闭型形式，便于后续数学处理，在满足估计准确性的条件下，为对海宽带雷达距离扩展目标检测器CFAR设计奠定扎实基础。

步骤B2在距离扩展目标知识辅助自适应融合检测器构建模块(3)中，依据式(18)和式(19)计算异质海杂波下距离扩展目标知识辅助自适应融合检测的检测统计量λ，并将λ送至检测判决模块(4)；

值得注意的是，针对海面目标不同散射点间相互独立的情况，本实施例中采用单个距离单元散射点复幅度a_m(m＝1,2,…,P)的最大似然估计表达式(17)进行检测器设计，可充分利用不同散射点间的独立性信息，提高特定场景下的检测性能；另外，在步骤B2中，所构建的海洋环境异质杂波下距离扩展目标知识辅助自适应融合检测器，在保持CFAR特性的同时，可进一步提高宽带对海雷达对异质海杂波环境的适应能力，改善海洋环境异质杂波下的距离扩展目标检测性能，提升复杂电磁环境下宽带雷达对弱小海面目标探测能力。

步骤B3根据预设的虚警概率设置检测门限T：具体为，设定虚警概率为P_fa，根据蒙特卡洛方法，依据前期积累的100/P_fa个实测数据计算检测阈值T；考虑到海杂波获取难度大，若实际获得的纯海杂波实测数据量Z少于100/P_fa，则缺少的100/P_fa-Z个杂波数据可利用海杂波仿真模型进行仿真获得，其中的模型参数根据已获得的纯海杂波实测数据进行合理估计设定。进一步将检测统计量λ与检测门限T进行比较，若λ≥T，则判定当前待检测距离单元存在距离扩展目标，主数据X不作为后续其他待检测距离单元的辅助数据；反之若λ＜T，则判定当前待检测距离单元不存在距离扩展目标，主数据X作为后续其他待检测距离单元的辅助数据。

Claims (4)

1.异质杂波下目标知识辅助自适应融合检测方法，其特征在于，包括以下步骤：

步骤1针对P个待检测距离单元的主数据X和K个辅助数据，利用逆Wishart分布对其中的杂波协方差矩阵结构进行先验分布建模，进而基于K个辅助数据分别构建杂波协方差矩阵结构的K个后验概率密度函数，将杂波协方差矩阵结构的K个后验概率密度函数进行几何平均，构建异质杂波协方差矩阵结构概率密度函数的几何平均估计方法，获得异质杂波协方差矩阵结构概率密度函数的几何平均估计；

步骤2基于异质杂波协方差矩阵结构概率密度函数的几何平均估计，在有目标和无目标两种假设下，求解主数据X的概率密度函数，根据广义似然比检验准则，分别求解距离扩展目标P个散射点未知复幅度的最大似然估计，针对距离扩展目标不同散射点之间是否统计独立的两种情况，分别构建异质杂波下目标知识辅助自适应融合检测的检测统计量；

步骤3根据预设的虚警概率设置检测门限T；将主数据对应的检测统计量λ与T进行比较，若λ≥T，则判定当前待检测距离单元存在距离扩展目标，X不作为后续其他待检测距离单元的辅助数据；反之若λ＜T，则判定当前待检测距离单元不存在距离扩展目标，X作为后续其他待检测距离单元的辅助数据。

2.根据权利要求1所述的异质杂波下目标知识辅助自适应融合检测方法，其特征在于，在所述步骤1中：

基于杂波协方差矩阵结构M的K个后验概率密度函数，采用几何平均方法估计M的概率密度函数，得到异质杂波协方差矩阵结构概率密度函数的几何平均估计

为

其中，exp[·]表示指数函数，函数tr(·)和det(·)分别表示对矩阵求迹和求行列式，Γ(·)表示Gamma函数，N表示雷达接收阵元数与相干处理脉冲数的乘积，矩阵R_k，k＝1,2,...K，具体表示为

上标“H”表示共轭转置，L表示逆Wishart分布的自由度，R₀为N×N维的Hermitian复数矩阵，K个辅助数据y_k，k＝1,2,...K，均为N×1维的复向量。

3.根据权利要求1所述的异质杂波下目标知识辅助自适应融合检测方法，其特征在于，在所述步骤2中：

针对距离扩展目标不同散射点之间非统计独立的情况，得到异质杂波下目标知识辅助自适应融合检测的检测统计量λ为

其中，I_N表示N×N维的单位矩阵；p表示已知的空时导向矢量，是一个N×1维的单位向量，根据雷达系统工作参数确定；det(·)表示对矩阵求行列式；矩阵R_k，k＝1,2,...K，表示为

上标“H”表示共轭转置，L表示逆Wishart分布的自由度，R₀为N×N维的Hermitian复数矩阵，K个辅助数据y_k，k＝1,2,...K，均为N×1维的复向量。

4.根据权利要求1所述的异质杂波下目标知识辅助自适应融合检测方法，其特征在于，在所述步骤2中：

针对距离扩展目标不同散射点之间统计独立的情况，得到异质杂波下目标知识辅助自适应融合检测的检测统计量λ为

其中，N×1维复向量x_m，m＝1,2,…,P，表示第m个主数据分量，|·|表示取复数的模值；p表示已知的空时导向矢量，是一个N×1维的单位向量，根据雷达系统工作参数确定；上标“H”表示共轭转置，L表示逆Wishart分布的自由度，R₀为N×N维的Hermitian复数矩阵，K个辅助数据y_k，k＝1,2,...K，均为N×1维的复向量。

Images