CN111999716B - 基于杂波先验信息的目标自适应融合检测方法 - Google Patents

基于杂波先验信息的目标自适应融合检测方法 Download PDF

Info

Publication number
CN111999716B
CN111999716B CN202010911710.5A CN202010911710A CN111999716B CN 111999716 B CN111999716 B CN 111999716B CN 202010911710 A CN202010911710 A CN 202010911710A CN 111999716 B CN111999716 B CN 111999716B
Authority
CN
China
Prior art keywords
clutter
target
distance
covariance matrix
heterogeneous
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010911710.5A
Other languages
English (en)
Other versions
CN111999716A (zh
Inventor
简涛
谢梓铿
卢仁伟
刘瑜
周坚毅
何友
刘书岩
但波
张健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Naval Aeronautical University
Original Assignee
Naval Aeronautical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Naval Aeronautical University filed Critical Naval Aeronautical University
Priority to CN202010911710.5A priority Critical patent/CN111999716B/zh
Publication of CN111999716A publication Critical patent/CN111999716A/zh
Application granted granted Critical
Publication of CN111999716B publication Critical patent/CN111999716B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414Discriminating targets with respect to background clutter

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于杂波先验信息的目标自适应融合检测方法,属于雷达信号处理领域。针对实际异质杂波环境下杂波协方差矩阵结构非均匀特点,构建易于数学处理的合理先验分布,将杂波先验信息与辅助数据中所含异质杂波信息进行算术平均融合,在满足估计准确性的条件下降低了估计算法的计算复杂度,针对距离扩展目标不同散射点之间独立与否,通过求和前去耦操作,设计简单易行的距离扩展目标散射点复幅度简化最大似然估计方法,分别构建了基于杂波先验信息的距离扩展目标自适应融合检测器,提高了宽带雷达对异质杂波环境的适应能力,改善了异质杂波下的目标检测性能,提升复杂电磁环境下宽带雷达对弱小目标的探测能力,具有推广应用价值。

Description

基于杂波先验信息的目标自适应融合检测方法
技术领域
本发明隶属于宽带雷达信号处理领域,具体涉及一种基于杂波先验信息的目标自适应融合检测方法。
背景技术
不同于低分辨率窄带雷达,宽带雷达带宽大、距离分辨率高,在抗干扰、反侦察、精确探测及成像、高精度跟踪、目标识别等方面具有明显的优势,在现代军事和民用领域引起了广泛重视,已经成为现代雷达发展的一个重要方向。窄带雷达的距离分辨单元一般远大于常见目标的几何尺寸,目标回波信号只占据一个距离分辨单元,实际目标常被作为“点目标”进行处理。而宽带雷达目标的回波信号不仅仅只占据一个距离分辨单元,而是分布在不同的径向距离分辨单元中,呈现为“一维距离像”,形成“距离扩展目标”。随着宽带雷达的广泛应用,距离扩展目标检测问题正受到越来越多的关注,成为近年来雷达信号处理界的热点和难点问题之一。
一方面,由于宽带雷达观测到的目标回波分布在多个径向距离单元上,若仍采用窄带雷达的点目标检测方法,针对单个距离单元对回波信号进行目标检测,并利用邻近的距离单元采样进行背景杂波统计特性估计,则距离扩展目标强散射点的能量会泄漏到邻近距离单元中形成“信号污染”现象,并进一步对单个待检测距离单元的目标信号构成遮蔽效应,导致无法检测到目标,由于未能充分利用分布在多个距离单元中目标回波的全部能量,这不但不能体现宽带雷达的优势,反而会降低对距离扩展目标的检测能力。另一方面,宽带雷达目标自适应检测面临目标自然环境复杂多变、电磁干扰等多种因素的影响,杂波统计特性不再满足独立同分布的均匀环境假设条件,这一问题在宽带雷达目标检测场景下尤为突出,导致现有距离扩展目标检测方法难以取得理想的检测效果。现有杂波模型主要包括均匀环境、部分均匀环境、非高斯杂波环境、异质杂波环境等。其中,部分均匀环境下距离扩展目标检测模型中,假设多个待检测距离单元的数据(也称主数据)与只含纯杂波的参考距离单元数据(也称辅助数据)中的杂波分量具有相同的协方差矩阵结构但功率水平不同;而以复合高斯分布建模的非高斯杂波环境中,假设主数据和辅助数据中的杂波分量具有相同的协方差矩阵结构但不同距离单元间杂波功率水平不同。前述两种杂波模型中,均假定杂波分量具有相同的协方差矩阵结构,但在宽带雷达探测面临的复杂海杂波、地物杂波等实际环境下,由于各种干扰因素的存在,不同距离单元间杂波协方差矩阵结构的均匀性会被进一步破坏,前述均匀、部分均匀和非高斯杂波模型中杂波分量具有相同协方差矩阵结构的假设条件将无法满足。此时需采用异质杂波模型进行环境建模,异质杂波环境下不同距离单元间杂波协方差矩阵结构相似但不相同,需通过合理的先验分布将辅助数据和主数据中的杂波协方差矩阵结构联系起来。
在宽带雷达探测场景下,针对实际异质杂波环境中不同距离单元间杂波协方差矩阵结构非均匀特点,如何利用合理先验分布,基于辅助数据中所含异质杂波信息来准确估计主数据中多个距离单元中的杂波统计信息,在兼顾主数据杂波协方差矩阵结构估计准确性的同时,建立简单易行的距离扩展目标多散射点幅度估计方法,进而构建具有闭合形式的距离扩展目标检测统计量,在增强检测器恒虚警率(CFAR)特性的前提下,进一步提高宽带雷达对异质杂波环境的适应能力,改善目标检测性能,是提升复杂电磁环境下宽带雷达探测能力的关键,也是目前急需解决的难题之一。
发明内容
在宽带雷达探测场景下,针对实际异质杂波环境中不同距离单元间杂波协方差矩阵结构非均匀特点,如何构建易于数学处理的合理先验分布,将杂波先验信息与辅助数据中所含异质杂波信息进行合理融合,构建简洁有效的协方差矩阵结构估计方法,在兼顾主数据杂波协方差矩阵结构估计准确性和算法运算量的同时为宽带雷达距离扩展目标检测器CFAR设计提供扎实基础,进而设计简单易行的距离扩展目标多散射点幅度估计方法,构建具有闭合形式的距离扩展目标知识辅助自适应融合检测器,在控制算法运算量的同时,通过考虑距离扩展目标不同散射点之间是否统计独立的具体情况,进一步提高宽带雷达对异质杂波环境的适应能力,改善异质杂波下的宽带雷达目标检测性能,提升复杂电磁环境下宽带雷达对弱小目标的探测能力。
本发明所述基于杂波先验信息的目标自适应融合检测方法包括以下技术措施:
步骤1从P个待检测距离单元获取主数据X,从与待检测距离单元临近的纯杂波距离单元获取K个辅助数据,利用逆Wishart分布对杂波协方差矩阵结构M进行先验分布建模,进而基于K个辅助数据分别构建M的K个后验概率密度函数,将K个后验概率密度函数进行算术平均,构建异质杂波协方差矩阵结构概率密度函数的算术平均估计方法,获得异质杂波协方差矩阵结构概率密度函数的算术平均估计
Figure BDA0002663534820000021
具体步骤包括:
由P个待检测距离单元回波复幅值构成P个待检测数据向量xm(m=1,2,…,P),进而构建N×P维的主数据X=[x1,x2,...,xP],即xm表示第m个主数据分量,P为大于1的自然数;以P个待检测距离单元为中心,在其前后分别连续取一定数量的不包含目标的距离单元回波复幅值,构成K个只含纯杂波的辅助数据yk(k=1,2,...K),其中,xm和yk均为N×1维的复向量,N表示雷达接收阵元数与相干处理脉冲数的乘积;
主数据和辅助数据中的杂波分量服从条件复高斯分布,但其杂波协方差矩阵结构M为随机矩阵,且服从自由度为L、均值为LR0的逆Wishart分布,其中,M和R0均为N×N维的Hermitian复数矩阵,M的先验分布概率密度函数(PDF)可表示为
Figure BDA0002663534820000022
其中,exp[·]表示指数函数,函数tr(·)和det(·)分别表示对矩阵求迹和求行列式,函数I[L,(LR0)-1]定义为:
Figure BDA0002663534820000031
上式中,Γ(·)表示Gamma函数。
令N×N维的复数矩阵Rk(k=1,2,...K)表示为
Figure BDA0002663534820000032
其中,上标“H”表示共轭转置。
给定第k个辅助数据yk时矩阵M的条件PDF(即以辅助数据yk为条件的M的后验PDF)可以表示为:
Figure BDA0002663534820000033
其中,f(yk|M)表示已知M时辅助数据yk的条件复高斯分布PDF,f(yk)表示yk的PDF,可由下式计算
Figure BDA0002663534820000034
结合杂波先验分布信息和K个辅助数据yk(k=1,2,...K)进行杂波协方差矩阵结构估计,直观而简洁的方式是进行算术平均,基于K个条件PDF f(M|yk)(k=1,2,...K),采用算术平均方法估计出矩阵M的PDF为:
Figure BDA0002663534820000035
上式为M的PDF算术平均估计表达式,可见其具有闭型形式且运算简洁,为后续设计闭合形式的距离扩展目标自适应融合检测统计量奠定了基础。
步骤2基于异质杂波协方差矩阵结构概率密度函数的算术平均估计,在有目标和无目标两种假设下求解主数据X的概率密度函数,根据简化广义似然比检验准则,忽略求和运算对距离扩展目标的P个散射点未知复幅度的耦合影响,采用求和前去耦,直接在每次求和运算前对P个散射点未知复幅度进行极值求解,融合杂波协方差矩阵结构先验分布信息与单个辅助数据中的杂波信息,逐一求解距离扩展目标散射点未知复幅度的K个简化最大似然估计,针对距离扩展目标不同散射点之间是否统计独立的两种情况,分别构建异质杂波环境下基于杂波先验信息的距离扩展目标自适应融合检测的检测统计量;具体步骤包括:
情况1:不考虑距离扩展目标不同散射点间的独立性,即针对距离扩展目标不同散射点之间非统计独立的情况;
在无目标的H0假设下,基于矩阵M的PDF的算术平均估计
Figure BDA0002663534820000041
主数据X的PDF可表示为
Figure BDA0002663534820000042
在有目标的H1假设下,基于矩阵M的PDF的算术平均估计
Figure BDA0002663534820000043
主数据X的PDF可表示为
Figure BDA0002663534820000044
其中,上述两式中,fX(X|M,H0)和fX(X|M,H1)分别表示H0和H1两种假设下主数据X在M已知时的复高斯分布条件PDF;1×P维的复幅度向量a表示距离扩展目标的未知复幅度,具体包含P个散射点的未知复幅度;p表示已知的空时导向矢量,是一个N×1维的单位向量,根据雷达系统工作参数确定;
根据广义似然比检验准则,距离扩展目标检测统计量可初步表示为
Figure BDA0002663534820000045
其中,IN表示N×N维的单位矩阵;
从式(9)可知,由于存在迭代求和运算,难以获得距离扩展目标未知复幅度向量a的最大似然估计,因此需对广义似然比检验准则进行简化修正,求解a的其他合理估计表达式;若忽略求和运算对a的耦合影响,即采用求和前去耦,直接在每次求和运算之前对a进行极值求解,可得a的K个简化最大似然(SML)估计为:
Figure BDA0002663534820000046
将式(10)的
Figure BDA0002663534820000047
代入式(9)中替换各求和项中的未知量a,并对式(9)中Λ1(X)的表达式进行等价数学变换,针对异质杂波环境,可得基于杂波先验信息的距离扩展目标自适应融合检测统计量λ1,即
Figure BDA0002663534820000051
在不考虑距离扩展目标不同散射点间独立性的情况下,可令
λ=λ1 (12)
情况2:针对距离扩展目标不同散射点间统计独立的情况;
在无目标的H0假设下,基于矩阵M的PDF的算术平均估计
Figure BDA0002663534820000052
主数据X的第m个分量xm的PDF可表示为
Figure BDA0002663534820000053
在有目标的H1假设下,基于矩阵M的PDF的算术平均估计
Figure BDA0002663534820000054
主数据X的第m个分量xm的PDF可表示为
Figure BDA0002663534820000055
其中,上述两式中,
Figure BDA0002663534820000056
Figure BDA0002663534820000057
分别表示H0和H1两种假设下xm在M已知时的复高斯分布条件PDF;1×P维的复幅度向量a可表示为a=[a1,a2,…,aP],即包含了P个散射点的未知复幅度am(m=1,2,…,P),其中,复标量am表示距离扩展目标第m个散射点的未知复幅度;
根据广义似然比检验准则,距离扩展目标检测统计量可表示为
Figure BDA0002663534820000058
从式(15)可知,由于存在迭代求和运算,难以获得未知量am(m=1,2,…,P)的最大似然估计,因此需对广义似然比检验准则进行简化修正,求解am(m=1,2,…,P)的其他合理估计表达式;若忽略求和运算对am(m=1,2,…,P)的耦合影响,即采用求和前去耦,直接在每次求和运算之前对am(m=1,2,…,P)进行极值求解,可得am(m=1,2,…,P)的K个简化最大似然估计为:
Figure BDA0002663534820000061
将式(16)的
Figure BDA0002663534820000062
代入式(15)中替换各求和项中的未知量am,并对式(15)中Λ2(X)的表达式进行等价数学变换,针对异质杂波环境,可得散射点独立时基于杂波先验信息的距离扩展目标自适应融合检测统计量λ2,即
Figure BDA0002663534820000063
上式中,|·|表示取复数的模值。
针对距离扩展目标不同散射点间统计独立的情况,可令
λ=λ2 (18)
步骤3为保持检测方法的CFAR特性,根据预设的虚警概率设置检测门限T;将主数据对应的检测统计量λ与检测门限T进行比较,若λ≥T,则判定当前待检测距离单元存在距离扩展目标,主数据不作为后续其他待检测距离单元的辅助数据;反之若λ<T,则判定当前待检测距离单元不存在距离扩展目标,主数据作为后续其他待检测距离单元的辅助数据。
与背景技术相比,本发明的有益效果是:1)针对实际异质杂波环境中不同距离单元间杂波协方差矩阵结构非均匀特点,构建了易于数学处理的杂波协方差矩阵结构合理先验分布;2)将杂波先验信息与辅助数据中所含异质杂波信息进行合理的算术平均融合,构建了简洁有效的协方差矩阵结构估计方法,在满足估计准确性的条件下降低了估计算法的计算复杂度,为距离扩展目标检测器CFAR设计奠定扎实基础;3)针对距离扩展目标不同散射点之间独立与否,通过求和前去耦操作,设计简单易行的距离扩展目标散射点复幅度简化最大似然估计方法,分别构建了基于杂波先验信息的距离扩展目标自适应融合检测器,均具有闭合形式的表达式,且在保持CFAR特性的同时降低了算法运算量,进一步提高了宽带雷达对异质杂波环境的适应能力,改善了异质杂波下的目标检测性能,提升复杂电磁环境下宽带雷达对弱小目标的探测能力;4)本发明方法适用于部分非宽带雷达探测情形,例如,使用低/中分辨率雷达检测大目标或检测以相同速度运动的空间邻近点目标群(舰艇编队、飞机编队、车辆编队等情况),具有很好的应用前景。
四、附图说明
图1是本发明所提出的基于杂波先验信息的目标自适应融合检测方法的功能模块图。图1中,1.中间矩阵计算模块,2.杂波协方差矩阵结构概率密度函数的算术平均估计模块,3.基于杂波先验信息的距离扩展目标自适应融合检测器构建模块,4.检测判决模块。
具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明做出的任何修改和改变,都落入本发明的保护范围。
为验证本发明所述方法的有效性,本具体实施方式给出了二个实施例,第一个实施例针对宽带雷达对地探测环境,第二个实施例针对宽带雷达对海探测环境。
实施例1:
参照说明书附图1,实施例1的具体实施方式分为以下几个步骤:
步骤A1利用对地探测宽带雷达,对待检测地域周围的无目标范围进行雷达照射,获得与待检测距离单元邻近的不包含目标的距离单元回波复幅值,构成K个只含纯地杂波的辅助数据yk(k=1,2,...K),将辅助数据送至中间矩阵计算模块(1);在中间矩阵计算模块(1)中,依据式(3)计算矩阵Rk(k=1,2,...K),并将矩阵Rk(k=1,2,…K)送至杂波协方差矩阵结构概率密度函数的算术平均估计模块(2);在杂波协方差矩阵结构概率密度函数的算术平均估计模块(2)中,依据式(6)计算杂波协方差矩阵结构M概率密度函数的算术平均估计
Figure BDA0002663534820000071
并将
Figure BDA0002663534820000072
送至基于杂波先验信息的距离扩展目标自适应融合检测器构建模块(3);
值得注意的是,步骤A1中,利用逆Wishart分布对杂波协方差矩阵结构进行先验分布建模,充分考虑了实际建筑物、林地、草地等不同地杂波环境导致的异质杂波情况,针对异质杂波下不同距离单元间杂波协方差矩阵结构非均匀特点,构建的地杂波协方差矩阵结构先验分布易于数学处理;另外,通过将杂波先验信息与辅助数据中所含异质杂波信息进行简洁有效的算术平均融合,所获得的协方差矩阵结构估计表达式具有闭型形式,兼顾了地杂波协方差矩阵结构估计准确性和算法运算量,便于后续数学处理,为对地探测宽带雷达距离扩展目标检测器CFAR设计奠定扎实基础。
步骤A2在基于杂波先验信息的距离扩展目标自适应融合检测器构建模块(3)中,依据式(11)和式(12)计算异质杂波下基于杂波先验信息的距离扩展目标自适应融合检测的检测统计量λ,并将λ送至检测判决模块(4);
值得注意的是,针对地面目标不同散射点间相关性较强的情况,本实施例中考虑距离扩展目标不同散射点间的非独立性,由于未知复幅度向量a的最大似然估计无法获得闭型解析表达式,进而采用求和前去耦的处理方式,直接在每次求和运算之前进行极值求解,得到a的SML估计表达式(10)进行检测器设计,可充分兼顾不同散射点间相关性的影响,构建了具有闭合形式的距离扩展目标知识辅助自适应融合检测器,提高了特定场景下的检测性能;另外,在步骤A2中,所构建的复杂地面环境异质杂波下基于杂波先验信息的距离扩展目标自适应融合检测器,在保持CFAR特性的同时控制了检测算法运算量,进一步提高了宽带对地雷达在异质地杂波环境下的适应能力,改善了复杂地面环境异质杂波下的距离扩展目标检测性能,提升了复杂电磁环境下宽带雷达对弱小地面目标的探测能力。本实施例方法同样适用于使用低/中分辨率雷达检测以相同速度运动的空间邻近车辆点目标群,具有很好的应用前景。
步骤A3根据预设的虚警概率设置检测门限T:具体为,设定虚警概率为Pfa,根据蒙特卡洛方法,依据前期积累的100/Pfa个实测地杂波数据计算检测阈值T。进一步,检测统计量λ与检测门限T进行比较,若λ≥T,则判定当前待检测距离单元存在距离扩展目标,主数据X不作为后续其他待检测距离单元的辅助数据;反之若λ<T,则判定当前待检测距离单元不存在距离扩展目标,主数据X作为后续其他待检测距离单元的辅助数据。
实施例2:
参照说明书附图1,实施例2的具体实施方式分为以下几个步骤:
步骤B1利用对海探测宽带雷达,对待检测海域周围的无目标范围进行雷达照射,获得与待检测距离单元邻近的不包含目标的距离单元回波复幅值,构成K个只含纯海杂波的辅助数据yk(k=1,2,...K),将辅助数据送至中间矩阵计算模块(1);在中间矩阵计算模块(1)中,依据式(3)计算矩阵Rk(k=1,2,...K),并将矩阵Rk(k=1,2,...K)送至杂波协方差矩阵结构概率密度函数的算术平均估计模块(2);在杂波协方差矩阵结构概率密度函数的算术平均估计模块(2)中,依据式(6)计算杂波协方差矩阵结构M概率密度函数的算术平均估计
Figure BDA0002663534820000081
并将
Figure BDA0002663534820000082
送至基于杂波先验信息的距离扩展目标自适应融合检测器构建模块(3);
值得注意的是,步骤B1中,利用逆Wishart分布对杂波协方差矩阵结构进行先验分布建模,充分考虑了不同海况海洋环境中杂波尖峰等导致的异质杂波情况,针对海洋环境异质杂波下不同距离单元间杂波协方差矩阵结构非均匀特点,构建的海杂波协方差矩阵结构先验分布易于数学处理;另外,通过将杂波先验信息与辅助数据中所含异质杂波信息进行简洁有效的算术平均融合,所获得的协方差矩阵结构估计表达式具有闭型形式,兼顾了海杂波协方差矩阵结构估计准确性和算法运算量,便于后续数学处理,为对海宽带雷达距离扩展目标检测器CFAR设计奠定扎实基础。
步骤B2在基于杂波先验信息的距离扩展目标自适应融合检测器构建模块(3)中,依据式(17)和式(18)计算异质海杂波下基于杂波先验信息的距离扩展目标自适应融合检测的检测统计量λ,并将λ送至检测判决模块(4);
值得注意的是,针对海面目标不同散射点间相互独立的情况,本实施例中考虑距离扩展目标不同散射点间统计独立性,由于存在迭代求和运算,未知量am(m=1,2,…,P)的最大似然估计无法获得闭型解析表达式,进而采用求和前去耦的处理方式,直接在每次求和运算之前进行极值求解,得到未知量am(m=1,2,…,P)的SML估计表达式(16)进行检测器设计,可充分利用不同散射点间的独立性信息,构建了具有闭合形式的距离扩展目标自适应融合检测器,提高了特定场景下的检测性能;另外,在步骤B2中,所构建的海洋环境异质杂波下基于杂波先验信息的距离扩展目标自适应融合检测器,在保持CFAR特性的同时控制了检测算法运算量,进一步提高了宽带对海雷达在异质海杂波环境下的适应能力,改善了复杂海洋环境异质杂波下的距离扩展目标检测性能,提升了复杂电磁环境下宽带雷达对弱小海面目标的探测能力。本发明方法同样适用于使用低/中分辨率雷达检测以相同速度运动的空间邻近舰艇编队点目标群,具有很好的应用前景。
步骤B3根据预设的虚警概率设置检测门限T:具体为,设定虚警概率为Pfa,根据蒙特卡洛方法,依据前期积累的100/Pfa个实测数据计算检测阈值T;考虑到海杂波获取难度大,若实际获得的纯海杂波实测数据量Z少于100/Pfa,则缺少的100/Pfa-Z个杂波数据可利用海杂波仿真模型进行仿真获得,其中的模型参数根据已获得的纯海杂波实测数据进行合理估计设定。进一步将检测统计量λ与检测门限T进行比较,若λ≥T,则判定当前待检测距离单元存在距离扩展目标,主数据X不作为后续其他待检测距离单元的辅助数据;反之若λ<T,则判定当前待检测距离单元不存在距离扩展目标,主数据X作为后续其他待检测距离单元的辅助数据。

Claims (4)

1.基于杂波先验信息的目标自适应融合检测方法,其特征在于,包括以下步骤:
步骤1针对P个待检测距离单元的主数据X和K个辅助数据,利用逆Wishart分布对其中的杂波协方差矩阵结构进行先验分布建模,进而基于K个辅助数据分别构建杂波协方差矩阵结构的K个后验概率密度函数,将杂波协方差矩阵结构的K个后验概率密度函数进行算术平均,构建异质杂波协方差矩阵结构概率密度函数的算术平均估计方法,获得异质杂波协方差矩阵结构概率密度函数的算术平均估计;
步骤2基于异质杂波协方差矩阵结构概率密度函数的算术平均估计,在有目标和无目标两种假设下求解主数据X的概率密度函数,根据简化广义似然比检验准则,忽略求和运算对距离扩展目标的P个散射点未知复幅度的耦合影响,采用求和前去耦,直接在每次求和运算前对P个散射点未知复幅度进行极值求解,融合杂波协方差矩阵结构先验分布信息与单个辅助数据中的杂波信息,逐一求解距离扩展目标散射点未知复幅度的K个简化最大似然估计,针对距离扩展目标不同散射点之间是否统计独立的两种情况,分别构建异质杂波环境下基于杂波先验信息的距离扩展目标自适应融合检测的检测统计量;
步骤3根据预设的虚警概率设置检测门限T;将主数据X对应的检测统计量λ与T进行比较,若λ≥T,则判定当前待检测距离单元存在距离扩展目标,X不作为后续其他待检测距离单元的辅助数据;反之若λ<T,则判定当前待检测距离单元不存在距离扩展目标,X作为后续其他待检测距离单元的辅助数据。
2.根据权利要求1所述的基于杂波先验信息的目标自适应融合检测方法,其特征在于,在所述步骤1中:
基于杂波协方差矩阵结构M的K个后验概率密度函数,采用算术平均方法估计M的概率密度函数,得到异质杂波协方差矩阵结构概率密度函数的算术平均估计
Figure FDA0003498004700000011
Figure FDA0003498004700000012
其中,exp[·]表示指数函数,函数tr(·)和det(·)分别表示对矩阵求迹和求行列式,N表示雷达接收阵元数与相干处理脉冲数的乘积,矩阵Rk,k=1,2,...K,具体表示为
Figure FDA0003498004700000013
上标“H”表示共轭转置,L表示逆Wishart分布的自由度,R0为N×N维的Hermitian复数矩阵,K个辅助数据yk,k=1,2,...K,均为N×1维的复向量;
Figure FDA0003498004700000014
Γ(·)表示Gamma函数。
3.根据权利要求1所述的基于杂波先验信息的目标自适应融合检测方法,其特征在于,在所述步骤2中:
针对距离扩展目标不同散射点之间非统计独立的情况,根据简化广义似然比检验准则,直接在每次求和运算之前对距离扩展目标未知复幅度向量a进行极值求解,融合杂波协方差矩阵结构先验分布信息与单个辅助数据中的杂波信息,得到a的K个简化最大似然估计,进而针对异质杂波环境,得到基于杂波先验信息的目标自适应融合检测的检测统计量λ为
Figure FDA0003498004700000021
其中,IN表示N×N维的单位矩阵;p表示已知的空时导向矢量,是一个N×1维的单位向量,根据雷达系统工作参数确定;P表示待检测距离单元的个数;矩阵Rk,k=1,2,...K,表示为
Figure FDA0003498004700000022
上标“H”表示共轭转置,L表示逆Wishart分布的自由度,R0为N×N维的Hermitian复数矩阵,K个辅助数据yk,k=1,2,...K,均为N×1维的复向量;det(·)表示对矩阵求行列式。
4.根据权利要求1所述的基于杂波先验信息的目标自适应融合检测方法,其特征在于,在所述步骤2中:
针对距离扩展目标不同散射点之间统计独立的情况,根据简化广义似然比检验准则,直接在每次求和运算之前对距离扩展目标所有P个散射点未知复幅度am,m=1,2,…,P,进行极值求解,融合杂波协方差矩阵结构先验分布信息与单个辅助数据中的杂波信息,得到am的K个简化最大似然估计,进而针对异质杂波环境,得到基于杂波先验信息的目标自适应融合检测的检测统计量λ为
Figure FDA0003498004700000023
其中,N×1维复向量xm,m=1,2,…,P,表示第m个主数据分量,|·|表示取复数的模值;矩阵Rk,k=1,2,…K,表示为
Figure FDA0003498004700000024
上标“H”表示共轭转置,L表示逆Wishart分布的自由度,R0为N×N维的Hermitian复数矩阵,K个辅助数据yk,k=1,2,…K,均为N×1维的复向量;p表示已知的空时导向矢量,是一个N×1维的单位向量,根据雷达系统工作参数确定。
CN202010911710.5A 2020-09-02 2020-09-02 基于杂波先验信息的目标自适应融合检测方法 Active CN111999716B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010911710.5A CN111999716B (zh) 2020-09-02 2020-09-02 基于杂波先验信息的目标自适应融合检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010911710.5A CN111999716B (zh) 2020-09-02 2020-09-02 基于杂波先验信息的目标自适应融合检测方法

Publications (2)

Publication Number Publication Date
CN111999716A CN111999716A (zh) 2020-11-27
CN111999716B true CN111999716B (zh) 2022-03-29

Family

ID=73464422

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010911710.5A Active CN111999716B (zh) 2020-09-02 2020-09-02 基于杂波先验信息的目标自适应融合检测方法

Country Status (1)

Country Link
CN (1) CN111999716B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731388B (zh) * 2020-12-14 2023-10-13 北京遥感设备研究所 一种基于有效散射点能量积累的目标检测方法
CN112731307B (zh) * 2020-12-19 2023-01-03 湖北工业大学 基于距离-角度联合估计的ratm-cfar检测器及检测方法
CN115453484B (zh) * 2022-08-29 2023-08-25 电子科技大学 一种自适应散射中心估计的距离扩展目标检测方法
CN115508828B (zh) * 2022-10-20 2024-05-14 中国人民解放军海军航空大学 子空间干扰下雷达目标智能融合检测方法
CN116106829B (zh) * 2022-11-28 2024-06-18 中国人民解放军海军航空大学 部分均匀杂波加干扰下目标智能检测方法
CN116577734B (zh) * 2023-07-13 2023-09-22 中国人民解放军空军预警学院 基于先验知识的机载雷达精细化杂波仿真方法与装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932766A (zh) * 2017-04-27 2017-07-07 中国人民解放军海军航空工程学院 基于变参数广义结构的距离扩展目标自适应检测方法
CN108919224A (zh) * 2018-07-26 2018-11-30 中国人民解放军海军航空大学 基于斜对称结构的宽带雷达目标自适应融合检测方法
CN110940970A (zh) * 2019-11-06 2020-03-31 河海大学 一种针对浮油海面的mimo雷达目标检测方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8724741B2 (en) * 2009-10-02 2014-05-13 Telefonaktiebolaget L M Ericsson (Publ) Signal quality estimation from coupling matrix
US8138963B1 (en) * 2010-09-10 2012-03-20 Mitsubishi Electric Research Laboratories, Inc. Method for detecting targets using space-time adaptive processing and shared knowledge of the environment
US10088555B2 (en) * 2014-12-15 2018-10-02 Airbus Singapore Private Limited Automated method for selecting training areas of sea clutter and detecting ship targets in polarimetric synthetic aperture radar imagery
WO2019047210A1 (zh) * 2017-09-11 2019-03-14 深圳大学 基于知识的稀疏恢复空时自适应处理方法及系统
CN108919225B (zh) * 2018-07-26 2020-06-30 中国人民解放军海军航空大学 部分均匀环境下距离扩展目标多通道融合检测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932766A (zh) * 2017-04-27 2017-07-07 中国人民解放军海军航空工程学院 基于变参数广义结构的距离扩展目标自适应检测方法
CN108919224A (zh) * 2018-07-26 2018-11-30 中国人民解放军海军航空大学 基于斜对称结构的宽带雷达目标自适应融合检测方法
CN110940970A (zh) * 2019-11-06 2020-03-31 河海大学 一种针对浮油海面的mimo雷达目标检测方法及系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
《Adaptive clutter suppression, sea mine detection/classification, and fusion processing string for sonar imagery》;Tom Aridgides, Manuel F. Fernandez, Gerald J. Dobeck;《Detection and Remediation Technologies for Mines and Minelike Targets IV》;19990802;1-8 *
《Adaptive persymmetric detector of generalised likelihood ratio test in homogeneous environment》;Tao Jian, You He, Guisheng Liao, Jian Guan, Yunlong Dong;《IET Signal Processing》;20160328;91-99 *
《空间非均匀海杂波的协方差矩阵估计新算法》;时艳玲;《系统工程与电子技术》;20130806;2239-2243 *
唐波 ; 张玉 ; 李科.《基于先验知识及其定量评估的自适应杂波抑制研究》.《航空学报》.2013, *
简涛 ; 黄晓冬 ; 王捷 ; 何友.《基于修正最大似然估计的距离扩展目标检测器》.《武汉大学学报(信息科学版)》.2016, *

Also Published As

Publication number Publication date
CN111999716A (zh) 2020-11-27

Similar Documents

Publication Publication Date Title
CN111999716B (zh) 基于杂波先验信息的目标自适应融合检测方法
CN111999714B (zh) 基于多散射点估计和杂波知识辅助的自适应融合检测方法
CN111999715B (zh) 异质杂波下目标知识辅助自适应融合检测方法
CN111999717B (zh) 基于协方差矩阵结构统计估计的自适应融合检测方法
CN108919223B (zh) 宽带雷达目标自适应结构化匹配滤波方法
CN108919225B (zh) 部分均匀环境下距离扩展目标多通道融合检测方法
KR102665123B1 (ko) 저피탐지를 위한 통합망 레이더 체류 시간과 방사 전력의 공동 최적화 방법
CN113238211B (zh) 一种干扰条件下参数化自适应阵列信号检测方法与系统
CN110389327A (zh) 接收站位置误差下多站多外辐射源雷达双基距定位方法
Yang et al. Enhanced knowledge-aided space–time adaptive processing exploiting inaccurate prior knowledge of the array manifold
CN115390027A (zh) 非均匀杂波下目标知识辅助智能融合检测方法
Wu et al. NLOS positioning for building layout and target based on association and hypothesis method
CN108562877A (zh) 一种基于信号包络特征的欺骗式干扰抑制方法
CN116299387B (zh) 非均匀杂波下干扰正交抑制的目标智能检测方法
CN116106829B (zh) 部分均匀杂波加干扰下目标智能检测方法
CN111999718B (zh) 基于几何平均估计的知识辅助自适应融合检测方法
CN112014823B (zh) 基于目标幅度综合估计的自适应融合检测方法
CN115524672A (zh) 结构化干扰与杂波下目标鲁棒智能检测方法
Peng et al. Constant false alarm rate detection of slow targets in polarimetric along‐track interferometric synthetic aperture radar imagery
CN118584444A (zh) 基于杂波结构化知识的目标智能融合检测方法
CN115267714B (zh) 一种基于dbf的杂波环境快速感知方法
CN118584445A (zh) 基于全数据两步估计的自适应融合目标智能检测方法
CN118795439A (zh) 基于备择假设全参数估计的目标智能融合检测方法
CN111856426B (zh) 基于中央厄米特结构和非匀质模型的子空间目标检测方法
Robey et al. Array calibration and modeling of steering vectors

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant