CN111965629B - 一种主动声呐非均匀背景抑制恒虚警检测方法及装置 - Google Patents

一种主动声呐非均匀背景抑制恒虚警检测方法及装置 Download PDF

Info

Publication number
CN111965629B
CN111965629B CN202011044952.5A CN202011044952A CN111965629B CN 111965629 B CN111965629 B CN 111965629B CN 202011044952 A CN202011044952 A CN 202011044952A CN 111965629 B CN111965629 B CN 111965629B
Authority
CN
China
Prior art keywords
amplitude signal
determining
unit
small
signal
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
CN202011044952.5A
Other languages
English (en)
Other versions
CN111965629A (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.)
Beijing Zhongke Haixun Digital Technology Co ltd
Original Assignee
Beijing Zhongke Haixun Digital Technology Co ltd
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 Beijing Zhongke Haixun Digital Technology Co ltd filed Critical Beijing Zhongke Haixun Digital Technology Co ltd
Priority to CN202011044952.5A priority Critical patent/CN111965629B/zh
Publication of CN111965629A publication Critical patent/CN111965629A/zh
Application granted granted Critical
Publication of CN111965629B publication Critical patent/CN111965629B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本申请实施例公开了一种主动声呐非均匀背景抑制恒虚警检测方法及装置,方法包括:接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,并估计背景功率;根据估计的背景功率,确定所述待检测单元中是否含有目标信号。本申请实施例能自适应的抑制非均匀背景,在多目标干扰背景环境中的目标检测有更强的鲁棒性,在混响边缘背景环境中有更好的虚警控制能力。

Description

一种主动声呐非均匀背景抑制恒虚警检测方法及装置
技术领域
本发明涉及主动声呐混响抑制领域。尤其涉及一种主动声呐非均匀背景抑制恒虚警检测方法及装置。
背景技术
在主动声呐目标检测中,混响是影响检测性能的主要干扰之一。混响是声波经水下散射体散射叠加而成的,在浅海环境下,海底的礁石、山脊及水中的鱼群等会造成强混响干扰;水体混响和界面混响等不同混响之间的混响背景功率突变会造成混响边缘干扰;高分辨率声呐当以低掠射角探测时,如果水中存在较大物体,由于其遮挡而在阴影部分没有回波从而形成阴影干扰。以上三种非均混响会对有源声呐的目标检测造成很大的影响,因此如何消除非均匀背景的影响、实现主动声纳在浅海环境下对水下目标的稳健检测是声呐工作者亟待解决的一个问题。
目前,混响抑制方法主要包括:声呐波形设计,空域滤波,空时联合处理,抗混响检测器设计等。由于干扰目标的多普勒频率和浅海慢速目标多普勒频率接近,且由于多普勒频率分辨率和空间角度分辨率的限制,声呐波形设计,空域滤波和空时联合处理很难对混响抑制做到普适应和彻底性。现有抗混响检测器主要包括有序统计恒虚警(OrderStatistic-Constant False Alarm Rate,OS-CFAR)检测器和自动删除平均恒虚警(Automatic Censored Cell Averaging-Constant False Alarm Rate,ACCA-CFAR)检测器,OS-CFAR需要干扰目标数的先验知识,这限制了该检测器在实际工作中的使用,ACCA-CFAR在混响边缘中虚警率控制能力不理想。
发明内容
由于现有方法存在上述问题,本申请实施例提出一种主动声呐非均匀背景抑制恒虚警检测方法及装置。
第一方面,本申请实施例提出一种主动声呐非均匀背景抑制恒虚警检测方法及装置,包括:
接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;
确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;
根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,并估计背景功率;
根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
在另一个可能的实现中,所述确定所述参考单元中的小幅值信号的的参考单元,包括:
将所有参考单元从小到大排序;
根据前n1 *个参考单元,确定小幅值信号的功率σ1,n1 *≤N;所述小幅值信号包括均匀混响或阴影干扰;
确定第n1 *+1个参考单元的功率与所述小幅值信号的功率σ1的比值R1
根据所述比值R1和功率比阈值ξ(ρ),确定第n1 *+1个参考单元是否为大幅值信号,所述大幅值信号包括干扰信号或均匀混响;
当第n1 *+1个参考单元是大幅值信号时,确定第n1 *+1个参考单元以后的所有参考单元都是大幅值信号,并且确定前n1 *个参考单元为小幅值信号参考单元;
当第n1 *+1个参考单元不是大幅值信号时,使n1 *加1,重复上述步骤。
在另一个可能的实现中,所述确定所述参考单元中的大幅值信号的参考单元,包括:
将所有参考单元从小到大排序;
根据最后n2 *个参考单元,确定大幅值信号的功率σ2,n2 *≤N;所述大幅值信号包括均匀混响或干扰目标数量多于n2 *时的多目标干扰信号;
确定第N-n2 *-1个参考单元的功率与所述大幅值信号的功率σ2的比值R2
根据所述比值R2和功率比阈值ξ(ρ),确定第N-n2 *-1个参考单元是否为小幅值信号,当第N-n2 *-1个参考单元是小幅值信号时,确定幅值小于第N-n2 *-1个参考单元的所有参考单元都是小幅值信号,并且确定最后n2 *个参考单元为大幅值信号的参考单元;
当第N-n2 *-1个参考单元不是小幅值信号时,使n2 *加1,重复上述步骤。
在另一个可能的实现中,所述根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,包括:
将所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *分别与N进行比较;
当N1 *和N2 *分别满足第一预设条件和第二预设条件时,确定小幅值信号的参考单元作为背景功率估计的数据;
当N1 *和N2 *分别满足第三预设条件和第四预设条件时,确定大幅值信号的参考单元作为背景功率估计的数据;
当N1 *和N2 *分别满足第五预设条件和第六预设条件时,确定大于第一预设阈值且小于第二预设阈值的参考单元作为背景功率估计的数据。
在另一个可能的实现中,所述估计背景功率,包括:
当N1 *和N2 *分别满足第一预设条件和第二预设条件时,将所述小幅值信号的参考单元取均值作为背景功率;
当N1 *和N2 *分别满足第三预设条件和第四预设条件时,将所述大幅值信号的参考单元取均值作为背景功率;
当N1 *和N2 *分别满足第五预设条件和第六预设条件时,将大于第一预设阈值且小于第二预设阈值的参考单元取均值作为背景功率。
在另一个可能的实现中,所述根据估计的背景功率,确定所述待检测单元中是否含有目标信号,包括:
根据估计的背景功率,确定自适应检测阈值;
比较所述待检测单元与所述自适应检测阈值的大小;若所述待检测单元大于或等于所述自适应检测阈值,则确定所述待检测单元中含有目标信号;否则,确定所述待检测单元中不含有目标信号。
在另一个可能的实现中,所述根据估计的背景功率,确定自适应检测阈值,包括:
将估计的背景功率与阈值因子相乘的结果作为所述自适应检测阈值。
第二方面,本申请实施例提出一种主动声呐非均匀背景抑制恒虚警检测装置,包括:
接收模块和确定模块;
所述接收模块,用于接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;
所述确定模块,用于确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;
所述确定模块,还用于根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,并估计背景功率;
所述确定模块,还用于根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
第三方面,本申请实施例还提出一种主动声呐非均匀背景抑制恒虚警检测装置,包括至少一个处理器,所述处理器用于执行存储器中存储的程序,当所述程序被执行时,使得所述装置执行:
接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;
确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;
根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定背景功率估计的数据,并估计背景功率;
根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
第四方面,本申请实施例还提出一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如下步骤:
接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;
确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;
根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定背景功率估计的数据,并估计背景功率;
根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
由上述技术方案可知,本申请实施例能自适应的识别出更多的干扰,在多目标干扰背景环境中的目标检测有更强的鲁棒性,在混响边缘背景环境中有更好的虚警控制能力。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些图获得其他的附图。
图1(a)为本申请实施例提供的小幅值信号的参考单元数量大于总参考单元数量的一半的示意图;
图1(b)为本申请实施例提供的大幅值信号的参考单元数量大于总参考单元数量的一半的示意图;
图1(c)为本申请实施例提供的小幅值和大幅值信号的参考单元数量同时小于总参考单元数量的一半的示意图;
图2为本申请实施例提供的非均匀背景抑制恒虚警检测器的模型框图;
图3为本申请实施例提供的均匀混响背景下本申请和单元平均恒虚警(CA-CFAR)检测器,有序统计恒虚警(OS-CFAR)检测器,自动删除平均恒虚警(ACCA-CFAR)检测器的仿真结果示意图;
图4(a)-图4(d)为本申请实施例提供的分别加6个干扰,10个干扰,14个干扰和18个干扰目标时本申请和CA-CFAR检测器,OS-CFAR检测器,ACCA-CFAR检测器的仿真结果示意图;
图5为本申请实施例提供的本申请和CA-CFAR检测器,OS-CFAR检测器,ACCA-CFAR检测器的混响边缘虚警率性能仿真结果示意图;
图6为本申请实施例提供的一种主动声呐非均匀背景抑制恒虚警检测方法的流程示意图;
图7为本申请实施例提供的一种主动声呐非均匀背景抑制恒虚警检测装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
需要说明的是,本申请实施例的说明书和权利要求书中的术语“第一”和“第二”等是用于区别不同的对象,而不是用于描述对象的特定顺序。例如,第一预设条件和第二预设条件等是用于区别不同的预设条件,而不是用于描述目标对象的特定顺序。在本申请实施例中,“示例性的”、“举例来说”或者“例如”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性的”、“举例来说”或者“例如”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言,使用“示例性的”、“举例来说”或者“例如”等词旨在以具体方式呈现相关概念。
在主动声呐目标检测中,主动声呐会发出探测信号,相应地也会接收到声呐回波信号。若恒虚警检测器接收的声呐回波信号中含有非均匀背景,则无法准确判断出声呐回波信号中是否含有目标信号。因此,为提高非均匀背景下恒虚警检测器检测目标的性能,本申请实施例提出非均匀背景抑制恒虚警检测器,下面对该检测器进行详细介绍。
首先,非均匀背景抑制恒虚警检测器接收N+1个声呐回波信号xi,N为正整数且为偶数。中间的一个声呐回波信号D作为待检测单元,用于检测中间的声呐回波信号D中是否含有目标信号。左边的声呐回波信号x1及右边的声呐回波信号/>到xN作为参考单元,用于估计背景功率。
其次,非均匀背景抑制恒虚警检测器会确定出小幅值信号的参考单元和大幅值信号的参考单元。首先介绍小幅值信号的参考单元的确定过程。具体地,将所有参考单元按幅值从小到大排序,即为y(1)<y(2)...<y(N)。使用前n1 *个参考单元初步估计小幅值信号的功率σ1,n1 *≤N;其中,小幅值信号包括均匀混响或阴影干扰。具体地,k1表示第k1个参考单元。计算第n1 *+1个参考单元的功率,并计算第n1 *+1个参考单元的功率与小幅值信号的功率σ1的比值为R1。当R1大于功率比阈值ξ(ρ)时,则确定第n1 *+1个参考单元是大幅值信号,该大幅值信号包括干扰信号或均匀混响,此时,确定第n1 *+1个参考单元以后的所有参考单元都是大幅值信号,停止判断,并且确定前n1 *个参考单元为小幅值信号的参考单元;当R1小于或等于功率比阈值ξ(ρ)时,则确定第n1 *+1个参考单元不是大幅值信号,使n1 *加1,并重复上述步骤。其中,/>在此,ρ为设定的参考单元功率与小幅值信号的功率的最低功率比。需要说明的是,参考单元的功率、小幅值信号的功率和功率比阈值满足如下公式:
然后介绍大幅值信号的参考单元的确定过程。具体地,将所有参考单元按幅值从小到大排序,即为y(1)<y(2)...<y(N)。使用最后n2 *个参考单元初步估计大幅值信号的功率σ2,n2 *≤N;其中,大幅值信号包括均匀混响或干扰目标数量多于n2 *时的多目标干扰信号。具体地,k2表示第k2个参考单元。计算第N-n2 *-1个参考单元的功率,并计算第N-n2 *-1个参考单元的功率与大幅值信号的功率σ2的比值为R2。当R2小于功率比阈值ξ(ρ)时,则确定第N-n2 *-1个参考单元是小幅值信号,该小幅值信号包括阴影干扰或均匀混响,在此,ρ为大幅值信号的功率与参考单元功率的最高功率比。此时,确定第N-n2 *-1个参考单元以后的所有参考单元都是小幅值信号,停止判断,并且确定最后n2 *个参考单元为大幅值信号的参考单元;当R2大于或等于功率比阈值ξ(ρ)时,则确定第N-n2 *-1个参考单元不是小幅值信号,使n2 *加1,并重复上述步骤。需要说明的是,参考单元功率、大幅值信号的功率和功率比阈值满足如下公式:
还需要说明的是,其中,N为参考单元数量,Q在正向判断中为干扰或强混响存在的概率,在反向判断中为小混响或阴影干扰存在的概率。
然后,非均匀背景抑制恒虚警检测器根据确定出的小幅值信号的参考单元可以确定其数目为N1 *并根据确定出的大幅值信号的参考单元可以确定其数目为N2 *。根据N1 *和N2 *,确定用于背景功率估计的数据,并使用确定出的背景功率估计的数据估计背景功率。具体地,将小幅值信号的参考单元的数目N1 *和大幅值信号的参考单元的数目N2 *分别与N进行比较;当N1 *和N2 *分别满足第一预设条件和第二预设条件时,确定小幅值信号的参考单元的均值作为背景功率估计的数据;举例来说,当且/>时,表明小幅值信号为均匀混响,确定选择小幅值信号的参考单元的均值作为背景功率估计的数据,以抑制多目标干扰和混响边缘中的大功率混响;此时,背景功率/>i表示第i个小幅值信号的参考单元;当N1 *和N2 *分别满足第三预设条件和第四预设条件时,确定大幅值信号的参考单元的均值作为背景功率估计的数据;举例来说,当/>且/>时,表明大幅值信号为均匀混响,确定选择大幅值信号的参考单元的均值作为背景功率估计的数据,以抑制阴影干扰(包括热噪声和环境噪声)和混响边缘背景中的小功率混响;此时,背景功率/>i表示第i个大幅值混响的参考单元;当N1 *和N2 *分别满足第五预设条件和第六预设条件时,确定大于第一预设阈值且小于第二预设阈值的参考单元的均值作为背景功率估计的数据;举例来说,当/>且/>时,表明中间的信号为均匀混响,确定选择第N1 *+1个参考单元至第N-N2 *个参考单元作为背景功率估计的数据,以抑制多目标干扰和混响边缘;此时,背景功率/>i表示第i个小幅值信号的参考单元。/>表示向下取整,/>表示向上取整。
最后,非均匀背景抑制恒虚警检测器将上述估计出的背景功率Z与阈值因子T相乘的结果作为自适应检测阈值。其中,阈值因子Pfa为恒虚警率值,R为去掉强混响后均匀背景参考单元的数量。需要说明的是,当/>且/>时,R=N1 *;当/>时,R=N2 *;当/>且/>时,R=N-N1 *-N2 *。比较待检测单元D与自适应检测阈值的大小。若待检测单元大于或等于自适应检测阈值,则确定待检测单元中含有目标信号;否则,确定待检测单元中不含有目标信号。
从声呐在实际的使用中可知非均匀背景的参考单元排序后可以分为图1所示的三种情况。其中N1 *表示小幅值信号的参考单元的数量,N2 *表示大幅值信号的参考单元的数量。图1(a)表示小幅值信号的参考单元的数量大于参考单元总数的一半,对应的非均匀背景包括多目标干扰或者强混响参考单元数量少于总参考单元数量一半的混响边缘,此时应该选择小幅值信号的参考单元的均值估计背景功率;图1(b)表示大幅值信号的参考单元的数量大于参考单元总数的一半,对应的非均匀背景包括阴影干扰或者强混响数量大于总参考单元一半的混响边缘,此时应该选择大幅值信号的参考单元的均值估计背景功率;图1(c)表示小幅值和大幅值信号的参考单元数量同时小于参考单元总数的一半,当N1 *=0,N2 *=0时对应的为均匀背景,此时应该使用全部的参考单元数据估计背景功率。当N1 *和N2 *不等于0时对应的为同时含有干扰(又称强混响)和阴影干扰(包括热噪声和环境噪声)的非均匀背景,此时应该把大幅值和小幅值数据去掉选择中间的数据,并使用选择的中间的数据的均值估计背景功率。
但在仿真中发现,当N1 *和N2 *以N/2作为分界时,混响边缘位置在N/2+1时虚警率会上升。这主要是因为混响边缘在这个位置时大幅值和小幅值数据各占一半,此时应该选择大幅值参考单元作为背景。但对于N1 *和N2 *的估计难度较大,检测器估计出的N1 *有可能大于N/2导致选择小幅值混响的参考单元作为背景功率估计的数据,导致虚警率上升。为了减小上述的影响,选择以N1 *>3/5N和N2 *>2/5N作为背景选择的判断标准,以使混响边缘位置在N/2+1时有更大的概率选择大幅值信号的参考单元作为背景功率估计的数据。
因此该检测器可以通过估计小幅值和大幅值信号的参考单元数量的估计值N1 *和N2 *,根据N1 *和N2 *的大小动态选择估计背景功率的数据以消除多目标干扰和混响边缘等非均匀背景对主动声呐目标恒虚警检测的影响。非均匀背景抑制恒虚警检测器的模型框图如图2所示。
图3显示了均匀混响背景下本发明和单元平均恒虚警(Cell Averaging-ConstantFalse Alarm Rate,CA-CFAR)检测器,有序统计恒虚警(OS-CFAR)检测器,自动删除平均恒虚警(ACCA-CFAR)检测器的仿真结果;其中横坐标SRR为信号和混响功率比,纵坐标pd为目标发现概率。图4(a)-图4(d)显示了分别加6个干扰(又称强混响),10个干扰,14个干扰和18个干扰目标时本申请和CA-CFAR检测器,OS-CFAR检测器,ACCA-CFAR检测器的仿真结果;其中横坐标SRR为信号和混响功率比,纵坐标pd为目标发现概率。图5显示了本申请和CA-CFAR检测器,OS-CFAR检测器,ACCA-CFAR检测器的混响边缘虚警率性能仿真结果;其中横坐标n代表混响边缘在参考单元中的位置,纵坐标pfa代表虚警概率。可以发现本申请所提的非均匀背景抑制恒虚警检测器不但在均匀背景和多目标干扰背景下具有较好的目标检测性能,在混响边缘情况下也有优异的虚警控制能力。
图6为本申请提供的一种主动声呐非均匀背景抑制恒虚警检测方法的流程示意图,该流程示意图包括:S601—S604,具体包括:
S601,接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数。
在本申请实施例中,非均匀背景抑制恒虚警检测器接收N+1个声呐回波信号xi,N为正整数且为偶数。将中间的一个声呐回波信号D作为待检测单元,其余的N个声呐回波信号,即声呐回波信号x1及声呐回波信号/>到xN作为参考单元。
S602,确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元。
在本申请实施例中,非均匀背景抑制恒虚警检测器确定出N个参考单元中的小幅值信号的参考单元和大幅值信号的参考单元。具体地,将所有参考单元从小到大排序;根据前n1 *个参考单元,确定小幅值信号的功率σ1,n1 *≤N;确定第n1 *+1个参考单元的功率与所述小幅值信号的功率σ1的比值R1;根据所述比值R1和功率比阈值ξ(ρ),确定第n1 *+1个参考单元是否为大幅值信号,所述大幅值信号包括干扰信号或均匀混响;当第n1 *+1个参考单元是大幅值信号时,确定第n1 *+1个参考单元以后的所有参考单元都是大幅值信号,并且确定前n1 *个参考单元为小幅值信号的参考单元;当第n1 *+1个参考单元不是大幅值信号时,使n1 *加1,重复上述步骤。将所有参考单元从小到大排序;根据最后n2 *个参考单元,确定大幅值信号的功率σ2,n2 *≤N;确定第N-n2 *-1个参考单元的功率与所述大幅值信号的功率σ2的比值R2;根据所述比值R2和功率比阈值ξ(ρ),确定第N-n2 *-1个参考单元是否为小幅值信号,所述小幅值信号包括阴影干扰或均匀混响;当第N-n2 *-1个参考单元是小幅值信号时,确定第N-n2 *-1个参考单元以前的所有参考单元都是小幅值信号,并且确定最后n2 *个参考单元为大幅值信号的参考单元;当第N-n2 *-1个参考单元不是小幅值信号时,使n2 *加1,重复上述步骤。
S603,根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,并估计背景功率。
在本申请实施例中,非均匀背景抑制恒虚警检测器根据小幅值信号的参考单元的数目N1 *和大幅值信号的参考单元的数目N2 *,确定背景功率估计的数据,并使用背景功率估计的数据估计背景功率。具体地,将所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *分别与N进行比较;当N1 *和N2 *分别满足第一预设条件和第二预设条件时,确定小幅值信号的参考单元的均值作为背景功率估计的数据;当N1 *和N2 *分别满足第三预设条件和第四预设条件时,确定大幅值信号的参考单元的均值作为背景功率估计的数据;当N1 *和N2 *分别满足第五预设条件和第六预设条件时,确定大于第一预设阈值且小于第二预设阈值的参考单元的均值作为背景功率估计的数据。在使用背景功率估计的数据估计背景功率时,具体做法为:当N1 *和N2 *分别满足第一预设条件和第二预设条件时,将小幅值信号的参考单元取均值作为背景功率;当N1 *和N2 *分别满足第三预设条件和第四预设条件时,将大幅值信号的参考单元取均值作为背景功率;当N1 *和N2 *分别满足第五预设条件和第六预设条件时,将大于第一预设阈值且小于第二预设阈值的参考单元取均值作为背景功率。
S604,根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
在本申请实施例中,将S603中估计的背景功率与阈值因子相乘的结果作为自适应检测阈值。比较待检测单元与自适应检测阈值的大小。若待检测单元大于或等于自适应检测阈值,则确定待检测单元中含有目标信号;否则,确定待检测单元中不含有目标信号。
图7为本申请提供的一种主动声呐非均匀背景抑制恒虚警检测装置的结构示意图,该结构示意图包括:接收模块701和确定模块702;
接收模块701,用于用于接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;
确定模块702,用于确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;
确定模块702,还用于根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,并估计背景功率;
确定模块702,还用于根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
本申请实施例提供一种主动声呐非均匀背景抑制恒虚警检测装置,其特征在于,包括至少一个处理器,所述处理器用于执行存储器中存储的程序,当所述程序被执行时,使得所述装置执行:
接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,并估计背景功率;根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
本申请实施例提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如下步骤:
接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;根据所述小幅值信号的参考单元的数目N1 *和所述大幅值信号的参考单元的数目N2 *,确定用于背景功率估计的数据,并估计背景功率;根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
应说明的是:以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。

Claims (9)

1.一种主动声呐非均匀背景抑制恒虚警检测方法,其特征在于,包括:
接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;
确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;
根据所述小幅值信号的参考单元的数目和所述大幅值信号的参考单元的数目/>确定用于背景功率估计的数据,并估计背景功率,其中将所述小幅值信号的参考单元的数目和所述大幅值信号的参考单元的数目/>分别与N进行比较;
和/>分别满足第一预设条件和第二预设条件时,确定小幅值信号的参考单元作为背景功率估计的数据;
和/>分别满足第三预设条件和第四预设条件时,确定大幅值信号的参考单元作为背景功率估计的数据;
和/>分别满足第五预设条件和第六预设条件时,确定大于第一预设阈值且小于第二预设阈值的参考单元作为背景功率估计的数据;
根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
2.根据权利要求1所述的方法,其特征在于,所述确定所述参考单元中的小幅值信号的参考单元,包括:
将所有参考单元从小到大排序;
根据前个参考单元,确定小幅值信号的功率σ1,/>所述小幅值信号包括均匀混响或阴影干扰;
确定第个参考单元的功率与所述小幅值信号的功率σ1的比值R1
根据所述比值R1和功率比阈值ξ(ρ),确定第个参考单元是否为大幅值信号,所述大幅值信号包括干扰信号或均匀混响;
当第个参考单元是大幅值信号时,确定第/>个参考单元以后的所有参考单元都是大幅值信号,并且确定前/>个参考单元为小幅值信号的参考单元;
当第个参考单元不是大幅值信号时,使/>加1,重复上述步骤。
3.根据权利要求1所述的方法,其特征在于,所述确定所述参考单元中的大幅值信号的参考单元,包括:
将所有参考单元从小到大排序;
根据最后个参考单元,确定大幅值信号的功率σ2,/>所述大幅值信号包括均匀混响或干扰目标数量多于/>时的多目标干扰信号;
确定第个参考单元的功率与所述大幅值信号的功率σ2的比值R2
根据所述比值R2和功率比阈值ξ(ρ),确定第个参考单元是否为小幅值信号,当第个参考单元是小幅值信号时,确定幅值小于第/>个参考单元的所有参考单元都是小幅值信号,并且确定最后/>个参考单元为大幅值信号的参考单元;
当第个参考单元不是小幅值信号时,使/>加1,重复上述步骤。
4.根据权利要求1所述的方法,其特征在于,所述估计背景功率,包括:
和/>分别满足第一预设条件和第二预设条件时,将所述小幅值信号的参考单元取均值作为背景功率;
和/>分别满足第三预设条件和第四预设条件时,将所述大幅值信号的参考单元取均值作为背景功率;
和/>分别满足第五预设条件和第六预设条件时,将大于第一预设阈值且小于第二预设阈值的参考单元取均值作为背景功率。
5.根据权利要求1所述的方法,其特征在于,所述根据估计的背景功率,确定所述待检测单元中是否含有目标信号,包括:
根据估计的背景功率,确定自适应检测阈值;
比较所述待检测单元与所述自适应检测阈值的大小;若所述待检测单元大于或等于所述自适应检测阈值,则确定所述待检测单元中含有目标信号;否则,确定所述待检测单元中不含有目标信号。
6.根据权利要求5所述的方法,其特征在于,所述根据估计的背景功率,确定自适应检测阈值,包括:
将估计的背景功率与阈值因子相乘的结果作为所述自适应检测阈值。
7.一种主动声呐非均匀背景抑制恒虚警检测装置,其特征在于,包括:
接收模块和确定模块;
所述接收模块,用于接收N+1个声呐回波信号,中间的一个声呐回波信号作为待检测单元,其余的N个声呐回波信号作为参考单元,N为正整数且为偶数;
所述确定模块,用于确定所述参考单元中的小幅值信号的参考单元和大幅值信号的参考单元;
所述确定模块,还用于根据所述小幅值信号的参考单元的数目和所述大幅值信号的参考单元的数目/>确定用于背景功率估计的数据,并估计背景功率,其中将所述小幅值信号的参考单元的数目/>和所述大幅值信号的参考单元的数目/>分别与N进行比较;当/>分别满足第一预设条件和第二预设条件时,确定小幅值信号的参考单元作为背景功率估计的数据;当/>和/>分别满足第三预设条件和第四预设条件时,确定大幅值信号的参考单元作为背景功率估计的数据;当/>和/>分别满足第五预设条件和第六预设条件时,确定大于第一预设阈值且小于第二预设阈值的参考单元作为背景功率估计的数据;
所述确定模块,还用于根据估计的背景功率,确定所述待检测单元中是否含有目标信号。
8.一种主动声呐非均匀背景抑制恒虚警检测装置,其特征在于,包括至少一个处理器,所述处理器用于执行存储器中存储的程序,当所述程序被执行时,使得所述装置执行:
如权利要求1-6任一项所述的方法。
9.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现如权利要求1至6任一所述的主动声呐非均匀背景抑制恒虚警检测方法。
CN202011044952.5A 2020-09-28 2020-09-28 一种主动声呐非均匀背景抑制恒虚警检测方法及装置 Active CN111965629B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011044952.5A CN111965629B (zh) 2020-09-28 2020-09-28 一种主动声呐非均匀背景抑制恒虚警检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011044952.5A CN111965629B (zh) 2020-09-28 2020-09-28 一种主动声呐非均匀背景抑制恒虚警检测方法及装置

Publications (2)

Publication Number Publication Date
CN111965629A CN111965629A (zh) 2020-11-20
CN111965629B true CN111965629B (zh) 2023-10-03

Family

ID=73387538

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011044952.5A Active CN111965629B (zh) 2020-09-28 2020-09-28 一种主动声呐非均匀背景抑制恒虚警检测方法及装置

Country Status (1)

Country Link
CN (1) CN111965629B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114280612B (zh) * 2021-12-13 2023-03-31 珠海安自达科技有限公司 一种绝缘体目标的毫米波雷达恒虚警检测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104391290A (zh) * 2014-11-17 2015-03-04 电子科技大学 一种适用于复杂非均匀杂波下的cfar检测器
KR20160039383A (ko) * 2014-10-01 2016-04-11 국방과학연구소 적응 배경 선택에 의한 일정 오경보율 검파 방법 및 장치
CN109188388A (zh) * 2018-09-03 2019-01-11 中国科学院声学研究所 一种对抗多目标干扰的恒虚警检测方法
CN109239677A (zh) * 2017-07-10 2019-01-18 中国科学院声学研究所 一种环境自适应恒虚警检测门限确定方法
CN109655834A (zh) * 2018-12-21 2019-04-19 中国海洋石油集团有限公司 基于恒虚警检测的多波束声呐测深方法及系统
CN111650591A (zh) * 2020-04-08 2020-09-11 中国船舶重工集团公司第七一五研究所 一种主动声纳宽带空间谱回波亮点增强及自动提取方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20160039383A (ko) * 2014-10-01 2016-04-11 국방과학연구소 적응 배경 선택에 의한 일정 오경보율 검파 방법 및 장치
CN104391290A (zh) * 2014-11-17 2015-03-04 电子科技大学 一种适用于复杂非均匀杂波下的cfar检测器
CN109239677A (zh) * 2017-07-10 2019-01-18 中国科学院声学研究所 一种环境自适应恒虚警检测门限确定方法
CN109188388A (zh) * 2018-09-03 2019-01-11 中国科学院声学研究所 一种对抗多目标干扰的恒虚警检测方法
CN109655834A (zh) * 2018-12-21 2019-04-19 中国海洋石油集团有限公司 基于恒虚警检测的多波束声呐测深方法及系统
CN111650591A (zh) * 2020-04-08 2020-09-11 中国船舶重工集团公司第七一五研究所 一种主动声纳宽带空间谱回波亮点增强及自动提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Fuzzy Statistical Normalization CFAR Detector for Non-Rayleigh Data;Yanwei Xu 等;《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》;第51卷(第1期);第383-396页 *
均匀混响背景下抗多目标干扰恒虚警检测器设计;殷超然 等;《水下无人系统学报》;第27卷(第4期);第434-441页 *
李清亮 等.《雷达地海杂波测量与建模》.国防工业出版社,2017,第430-435页. *

Also Published As

Publication number Publication date
CN111965629A (zh) 2020-11-20

Similar Documents

Publication Publication Date Title
CN104569948B (zh) 海杂波背景下子带自适应glrt‑ltd检测方法
US5900835A (en) Coherent hidden markov model
CN105785330B (zh) 一种认知型副瓣干扰抑制方法
Abraham et al. Reverberation envelope statistics and their dependence on sonar bandwidth and scattering patch size
CN103117064A (zh) 处理信号
CN109521430B (zh) 一种抑制窄带干扰的距离扩展目标检测方法
CN109188388B (zh) 一种对抗多目标干扰的恒虚警检测方法
RU2711406C1 (ru) Способ классификации гидроакустических сигналов шумоизлучения морских объектов
CN111965629B (zh) 一种主动声呐非均匀背景抑制恒虚警检测方法及装置
CN112180341B (zh) 一种背景自适应cfar算法选择的实现方法
KR101817011B1 (ko) 군집화 특성 기반의 능동소나 클러터 제거 방법 및 장치
RU2305297C2 (ru) Устройство с корреляционным формирователем характеристики направленности для обнаружения сигналов и определения направления на их источник
JP3785715B2 (ja) 目標検出装置
CN110441761B (zh) 基于分布式浮标探测的多源信息融合方法
CN111796253A (zh) 基于稀疏信号处理的雷达目标恒虚警检测方法
CN110673118A (zh) 一种主动声纳单频脉冲串波形设计及检测算法
CN115616514A (zh) 一种基于背景熵的密集目标快速检测方法及系统
CN111123274B (zh) 一种水下声呐成像系统目标检测方法
RU2736188C1 (ru) Способ отображения гидроакустической информации
Park et al. Feasibility of range estimation using sonar LPI
JP2002168937A (ja) 水中目標物位置検出装置及び方法
Wage et al. Performance analysis of dominant mode rejection beamforming
CN108680904A (zh) 一种船用导航雷达同频干扰的抑制方法
Garvanov et al. Sensitivity of API CFAR detectors towards change of input parameters of pulse jamming
CN111580056A (zh) 适用于高分辨雷达的迭代形式的恒虚警检测方法

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20220129

Address after: 100000 room s6306, floor 3, building 1, yard 33, Dijin Road, Haidian District, Beijing

Applicant after: BEIJING ZHONGKE HAIXUN DIGITAL TECHNOLOGY Co.,Ltd.

Address before: 100190, No. 21 West Fourth Ring Road, Beijing, Haidian District

Applicant before: INSTITUTE OF ACOUSTICS, CHINESE ACADEMY OF SCIENCES

GR01 Patent grant
GR01 Patent grant