CN103033815B - 基于混响协方差矩阵的距离扩展目标的检测方法和装置 - Google Patents

基于混响协方差矩阵的距离扩展目标的检测方法和装置 Download PDF

Info

Publication number
CN103033815B
CN103033815B CN201210556458.6A CN201210556458A CN103033815B CN 103033815 B CN103033815 B CN 103033815B CN 201210556458 A CN201210556458 A CN 201210556458A CN 103033815 B CN103033815 B CN 103033815B
Authority
CN
China
Prior art keywords
covariance matrix
pseudo
data
formulation information
reverberation
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
CN201210556458.6A
Other languages
English (en)
Other versions
CN103033815A (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.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
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 Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201210556458.6A priority Critical patent/CN103033815B/zh
Publication of CN103033815A publication Critical patent/CN103033815A/zh
Application granted granted Critical
Publication of CN103033815B publication Critical patent/CN103033815B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种基于混响协方差矩阵的距离扩展目标的检测方法,包括:通过由声纳阵获取一组回波数据并从中选取据有相同的混响协方差矩阵的主数据和辅数据;根据主数据、辅数据和混响协方差矩阵的伪对称特性计算采样协方差矩阵的伪对称表示信息,主数据的伪对称表示信息和混响协方差矩阵的最大似然估计的伪对称表示信息,混响协方差矩阵的结构具有预先确定的伪对称特性,然后辅之以距离扩展目标的标称导向向量计算检测统计量;根据检测统计量与预定阈值的比较结果判断距离扩展目标是否存在。本发明利用了混响协方差矩阵的先验特性,在小样本辅助数据情况下具有非常好的检测性能,有效地提升了高分辨主动声纳在非均匀环境下的性能。

Description

基于混响协方差矩阵的距离扩展目标的检测方法和装置
技术领域
本发明涉及声纳技术,尤其涉及一种基于混响协方差矩阵的距离扩展目标的检测的方法和装置
背景技术
低分辨率主动声纳的距离分辨率远大于军事目标尺寸,实际目标具有点目标的特征,其回波与发射信号具有相同的形式。因此,混响背景下的目标回波检测问题可以归结为具有随机参量的已知波形的检测问题,相应的检测器以发射信号的匹配滤波器为核心。
相对于深海环境,工作于浅海的主动声纳受到混响的影响更大,尤其是海底混响。为了补偿检测性能上的损失,高分辨率主动声纳应运而生。由于距离分辨率的提高,当距离分辨单元尺寸明显小于目标尺寸时,目标回波将占据多个距离单元,这种情况下的检测目标被称为距离扩展目标。此时,对发射信号的匹配滤波只是对各径向散射中心回波的最佳检测,而没有充分利用目标回波的全部能量,并不是对整个目标回波的最佳检测方法。要充分体现高分辨率声纳的性能优势,就需要研究适用于距离扩展目标的检测方法。因此,国内外学者针对不同的混响背景,在点目标检测的基础上对距离扩展目标检测进行了深入而又广泛的研究。其中较为经典的是Rao检验方法,该方法假设可以获得均匀的辅助数据来估计待测距离单元(主数据)的混响协方差矩阵,从而构造自适应检测统计量。辅助数据一般从主数据临近的距离单元获得,如果它和主数据具有相同的混响协方差矩阵,则称其为均匀的。如果可以获得充足的辅助数据,Rao检验方法将具有非常好的检测性能。
然而,声纳在实际应用中通常工作在非均匀环境。造成非均匀现象的原因有很多种,如多变的海底地表类型、鱼群、孤立干扰等。在这种情况下,均匀的辅助数据很难获得,即使可以获得,其样本数也非常有限,从而导致现有的距离扩展目标的检测性能大大下降。
发明内容
本发明的目的是提供一种能够克服上述缺陷的基于先验知识的距离扩展目标检测方法和装置。
在本发明的第一方面,提供了一种基于混响协方差矩阵的距离扩展目标的检测方法,包括:通过由声纳阵接收的回波获取一组回波数据,选取所述一组回波数据中的多个回波数据作为主数据,选取所述一组回波数据中除所述主数据外的多个回波数据作为辅数据,所述主数据和所述辅数据具有相同的混响协方差矩阵;根据所述主数据、所述辅数据和混响协方差矩阵的伪对称特性来计算采样协方差矩阵的伪对称表示信息,所述主数据的伪对称表示信息和所述混响协方差矩阵的最大似然估计的伪对称表示信息,其中,所述混响协方差矩阵的结构具有预先确定的伪对称特性;以及根据所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量,计算检测统计量;以及将所述检测统计量与预定阈值进行比较,并且根据比较结果判断所述距离扩展目标是否存在。
在本发明的第二方面,提供了一种基于混响协方差矩阵的距离扩展目标的检测装置,包括:回波数据获取单元,用于通过由声纳阵接收的回波获取一组回波数据,选取所述一组回波数据中的多个回波数据作为主数据,选取所述一组回波数据中除所述主数据外的多个回波数据作为辅数据,所述主数据和所述辅数据具有相同的混响协方差矩阵;第一计算单元,用于根据所述辅数据计算采样协方差矩阵;第二计算单元,用于分别计算所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息和混响协方差矩阵的最大似然估计的伪对称表示信息,其中,所述混响协方差矩阵的结构具有预先确定的伪对称特性;以及第三计算单元,用于根据所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量,计算检测统计量;以及比较判断单元,用于将所述检测统计量与预定阈值进行比较,并且根据比较结果判断所述距离扩展目标是否存在。
在第三方面,本发明提供了一种声纳系统,包括:回波数据获取单元,用于通过由声纳阵接收的回波获取一组回波数据,选取所述一组回波数据中的多个回波数据作为主数据,选取所述一组回波数据中除所述主数据外的多个回波数据作为辅数据,所述主数据和所述辅数据具有相同的混响协方差矩阵;第一计算单元,用于根据所述辅数据计算采样协方差矩阵;第二计算单元,用于分别计算所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息和混响协方差矩阵的最大似然估计的伪对称表示信息,其中,所述混响协方差矩阵的结构具有预先确定的伪对称特性;以及第三计算单元,用于根据所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量,计算检测统计量;比较判断单元,用于将所述检测统计量与预定阈值进行比较,并且根据比较结果判断所述距离扩展目标是否存在;以及声纳阵,所述声纳阵连接到所述回波数据获取单元,所述声纳阵用于发射探测声波并接收回波。
本发明充分利用了混响协方差矩阵的先验信息,在小样本辅助数据情况下具有非常好的检测性能,有效地提升了高分辨主动声纳在非均匀环境下的性能。
附图说明
图1是根据本发明实施例的距离扩展目标的检测流程图;
图2是根据本发明实施例的距离扩展目标的检测装置的示意图
图3是根据本发明实施例的第一种情况下检测性能比较的曲线图;以及
图4是根据本发明实施例的第二种情况下检验性能比较的曲线图。
具体实施方式
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
本发明是利用外部环境提供的先验知识,来克服背景的非均匀性,提高声纳系统性能的。众所周知,对于我们常用的等间距线阵,其混响协方差矩阵满足伪对称特性。为此,本发明充分利用这一先验信息,基于Rao检验准则,提出一种知识基距离扩展目标的检测方法。
1.问题描述
假设目标回波是由N个阵元组成的等间距线阵接收,则距离扩展目标检测可归结为如下二元假设:
H 0 : z t = n t , t = 1 , . . . , L x t = n t , t = 1 , . . . , K H 1 : z t = α t v + n t , t = 1 , . . . , L x t = n t , t = 1 , . . . , K - - - ( 1 )
其中H0和H1分别代表目标不存在假设和目标存在假设;zt,t=1,...,L表示待检测单元的采样数据,又称为主数据,长度为L;xt,t=1,...,K表示K个均匀的辅助数据,它们与主数据具有相同的混响协方差矩阵;nt是独立的、零均值复合高斯混响,其协方差矩阵为t=1,...,L+K,其中H表示矩阵或向量的共轭转置;αt,t=1,...,L是目标的幅度;v是目标的标称导向向量。
为了便于检测方法的设计,可以定义一些简化表示,具体有:
·ZL=[z1,…,zL]和ZK=[x1,…,xK]分别表示主数据矩阵和辅数据矩阵;
·α=[α1,…,αL]T是一个L维的列向量,其中(·)T表示矩阵或向量的转置;
·θr=[α1,R,I2,R,α2,I,…,αL,RL,I]T是一个2L维的信号参数向量,其中αt,R和αt,I分别是αt,t=1,…,L的实部和虚部;·为N2+2L+1维的列矢量,其中θs为N2维的多余参数列向量。
为了利用混响协方差矩阵的先验信息,我们假设M是一个正定的伪对称矩阵,即它隶属于由式(2)定义的集合P
M∈P如果M=JM*J    (2)
其中*代表矩阵或向量的复共轭,矩阵J是一个N维的反对角阵,即
J = 0 0 . . . 0 1 0 0 . . . 1 0 . . . . . . . . . . . . . . . 1 0 . . . 0 0 - - - ( 3 )
同时,我们假设目标标称导向向量v也具有伪对称结构,即v=Jv*。对于等间距线阵,上述假设是非常合理的。
2.检测器设计
在上诉假设条件下,Rao检验可以表示为:
其中
·η是由虚警概率确定的检测门限值;
·f(ZL,ZK|θ,M)是H1下主、辅数据的联合概率密度函数,即
f(ZL,ZK|θ,M)=π-N(L+K)||M||-(K+H)
×exp{-tr[M-1((ZL-vαH)(ZL-vαH)H+S)]},
其中S是采样协方差矩阵,即||·||和tr(·)分别代表矩阵的行列式和矩阵迹;·
∂ f ( Z L , Z K | θ , M ) ∂ θ r = [ ∂ f ( Z L , Z K | θ , M ) ∂ α 1 , R , ∂ f ( Z L , Z K | θ , M ) ∂ α 1 , I
. . . , ∂ f ( Z L , Z K | θ , M ) ∂ α L , R , ∂ f ( Z L , Z K | θ , M ) ∂ α L , I ] T ,
其中分别代表相对于αt,R和αt,I的偏导数;
·J(θ)=J(θrs)是Fisher信息矩阵(Fisher Information Matrix),满足
[ J - 1 ( θ ) ] θ r , θ r = 1 2 ( v H M - 1 v ) - 1 I 2 H ,
其中I2H表示2H×2H维的单位矩阵;·是H0下相关参数的最大似然估计。
对于我们的问题,经推导可得
∂ f ( Z L , Z K | θ , M ) ∂ α t , R = 2 Re { v H M - 1 ( z t - α t v ) }
∂ f ( Z L , Z K | θ , M ) ∂ α t , I = 2 Im { v H M - 1 ( z t - α t v ) } - - - ( 5 )
其中Re(·)和Im(·)分别表示复数的实部和虚部。混响协方差矩阵M的伪对称最大似然估计为 M ^ p = 1 K + H [ S p + X p X p H ] , - - - ( 6 )
其中Sp和Xp分别代表采样协方差矩阵的伪对称表示和主数据的伪对称表示,即Sp=(S+JS*J)/2,Xp=[ze1,…,zeL,zo1,…,zoL], t=1,…,L。
将(5)、(6)代入(2),经推导,知识基Rao方法的检测统计量为:
T = v H ( S p + X p X p H ) - 1 X p X p H ( S p + X p + X p H ) - 1 v v H ( S p + X p X p H ) - 1 v - - - ( 8 )
可以看出T的计算是在复数域完成的。实际上,T的计算公式还可以进一步简化。为此,可以引入酉矩阵U,其中U的实部和虚部分别为(I+J)/2和(I-J)/2,I为N×N维的单位矩阵,J如(3)所示。容易得到
vr=Uv=Re(v)-Im(v),
Sr=USpUH=Re(Sp)-JIm(Sp),
Xr=Re(UXp)+Im(UXp)=[zer1,…,zerL,zor1,…,zorL]    (9)
可以看出,vr是实数向量,而Sr和Xr均为实数矩阵。
将(9)带入(8),T可以写为
T = v r H M ^ r - 1 X r X r H M ^ r - 1 v r v r H M ^ r - 1 v r - - - ( 10 )
其中,显而易见,T的计算全部由实数操作完成。
3.性能分析
我们借助蒙特-卡罗仿真方法来分析本发明方法的性能,并与现有的Rao检验方法进行比较。仿真中考虑的是一个由N个阵元组成的等间距线阵,具体的参数设置为N=9,L=3,虚警概率Pfa=10-4,目标标称导向向量v=[1,...,1]T/N,信混比混响为常见的指数相关复合高斯向量,其协方差矩阵M=0.9|i-j|,其中(i,j)为矩阵元素的坐标。
下面,通过具体实施例并结合附图对本发明的检测方法和装置进行详细描述。
图1是根据本发明实施例的距离扩展目标的检测流程图。
首先,声纳阵接收回波数据。距离扩展目标的回波可能占据多个距离单元。回波数据中包含主数据和多个辅数据,其中,主数据为待检测单元数据,主数据是不一定包含目标数据的待检测数据,例如对于目标不存在假设,主数据就不包含目标能量,而辅数据可以为主数据临近的不包含目标能量的数据。例如,声纳阵可以有N个阵元,每个阵元接收L个待检测单元数据(主数据)rt,t=1,…,L,以及K个不包含目标能量的辅助数据rt,t=L+1,…,L+K。
接下来,根据辅助数据计算采样协方差矩阵S,分别计算主数据的伪对称表示、采样协方差矩阵的伪对称表示和混响协方差矩阵最大似然估计的伪对称表示,分别为Xp,Sp混响协方差矩阵最大似然估计的伪对称表示可以根据主数据的伪对称表示和采样协方差矩阵的伪对称表示来计算,
M ^ p = 1 L + K [ S p + X p + X p H ] .
优选地,可以采用酉矩阵U对复数域的主数据的伪对称表示、采样协方差矩阵的伪对称表示、混响协方差矩阵的最大似然估计的伪对称表示以及待检测目标的标称导向向量进行变换,其中,目标的标称导向向量由已知的目标方向计算得到。酉矩阵U的实部和虚部分别为(I+J)/2和(I-J)/2,I为N×N维的单位矩阵,J是一个N×N维的反对角阵,即
J = 0 0 . . . 0 1 0 0 . . . 1 0 . . . . . . . . . . . . . . . 1 0 . . . 0 0
通过变换,将Xp,Sp及目标的标称导向向量v变换为实数矩阵Xr,Sr及实数向量vr
然后,根据上面得到的实数域的主数据矩阵、采样协方差矩阵、混响协方差矩阵以及待检测目标的标称导向向量,计算检测统计量T。如上所述,将公式(9)带入公式(8),可以得出T的计算表达式。
最后,将检测统计量T与预设定的检测门限η进行比较,以便确定距离扩展目标是否存在。如果检测统计量T大于或等于该检测门限η,则距离扩展目标存在;反之,如果检测统计量T小于该检测门限η,则距离扩展目标不存在。检测门限可以采用蒙特-卡罗仿真方法基于实验数据预先获得。
图2是根据本发明实施例的距离扩展目标的检测装置的示意图。
首先,通过声纳阵中的多个阵元接收回波,距离扩展目标的回波可能占据多个距离单元。回波数据中包含主数据和多个辅数据,其中,主数据为待检测单元数据,辅数据可以为主数据临近的不包含目标能量的数据。例如,声纳阵可以有N个阵元,每个阵元接收L个待检测单元数据(主数据)rt,t=1,…,L,以及K个不包含目标能量的辅助数据rt,t=L+1,…,L+K。
接下来,第一计算单元根据辅助数据计算采样协方差矩阵S。第二计算单元分别计算主数据的伪对称表示、采样协方差矩阵的伪对称表示和混响协方差矩阵最大似然估计的伪对称表示,分别为Xp,Sp混响协方差矩阵最大似然估计的伪对称表示可以根据主数据的伪对称表示和采样协方差矩阵的的伪对称表示来计算, M ^ p = 1 L + K [ S p + X p + X p H ] .
优选地,酉变换单元采用酉矩阵U对复数域的主数据、采样协方差矩阵、混响协方差矩阵的最大似然估计的伪对称表示以及待检测目标的标称导向向量进行变换。酉矩阵U的实部和虚部分别为(I+J)/2和(I-J)/2,I为N×N维的单位矩阵,J是一个N×N维的反对角阵,即
J = 0 0 . . . 0 1 0 0 . . . 1 0 . . . . . . . . . . . . . . . 1 0 . . . 0 0
通过变换,将Xp,Sp及目标的标称导向向量v变换为实数矩阵Xr,Sr及实数向量vr
然后,第三计算单元根据上面得到的实数域的主数据矩阵、采样协方差矩阵、混响协方差矩阵以及待检测目标的标称导向向量,计算检测统计量T。如上所述,将公式(9)带入公式(8),可以得出T的计算表达式。
最后,比较判断单元将检测统计量T与预设定的检测门限η进行比较,以便确定距离扩展目标是否存在。如果检测统计量T大于或等于该检测门限η,则距离扩展目标存在;反之,如果检测统计量T小于该检测门限η,则距离扩展目标不存在。检测门限可以采用蒙特-卡罗仿真方法基于实验数据预先获得。
图3和图4示出在两种不同辅助数据长度情况下,采用本发明方法与采用Rao检验方法时检测概率Pd和SRR的关系曲线,本发明方法在图3和图4中也称作知识基方法,意指本本发明的方法是基于混响协方差矩阵的先验的伪对称特性的一种方法。在图3中,设定K=15,而在图4中设定K=25,K为辅助数据的数量;在图3和图4中,设定N都为9,L都为3,N为声纳阵的阵元数量,L为待检测单元的数量。在图3中的情况下,N=9,K=15,L=3,为满足检测概率Pd=0.5,本发明方法所需的信混比为13.1dB,而Rao检验方法所需的信混比为17.6dB,即Rao检验方法相对于本发明方法的检测损失为4.5dB。同理,在图4中,Rao检验方法相对于本发明方法的检测损失要低于图3。由此可以看出,可获得的辅助数据(K)越少,本发明的方法的性能优势越明显。可见,本发明方法合理地利用了外部环境提供的先验知识,有效地提高了高分辨主动声纳系统在非均匀环境下的性能。
专业人员应该还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于混响协方差矩阵的距离扩展目标的检测方法,包括:通过由声纳阵接收的回波获取一组回波数据,选取所述一组回波数据中的多个回波数据作为主数据,选取所述一组回波数据中除所述主数据外的多个回波数据作为辅数据,所述主数据和所述辅数据具有相同的混响协方差矩阵,其特征在于,
根据所述主数据、所述辅数据和混响协方差矩阵的伪对称特性来计算采样协方差矩阵的伪对称表示信息,所述主数据的伪对称表示信息和所述混响协方差矩阵的最大似然估计的伪对称表示信息,其中,所述混响协方差矩阵的结构具有预先确定的伪对称特性;以及
根据所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量,计算检测统计量;
将所述检测统计量与预定阈值进行比较,并且根据比较结果判断所述距离扩展目标是否存在。
2.根据权利要求1所述的检测方法,所述采样协方差矩阵是根据所述辅数据计算得到的。
3.根据权利要求1所述的检测方法,其中,所述混响协方差矩阵的最大似然估计的伪对称表示信息是基于所述主数据的伪对称表示信息和所述采样协方差矩阵的伪对称表示信息计算得到的。
4.根据权利要求1所述的检测方法,其中,所述根据所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量,计算检测统计量的步骤包括:
利用酉矩阵将所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量从复数域变换至实数域;以及
根据所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量变换结果,计算所述检测统计量。
5.根据权利要求3所述的检测方法,其中,所述酉矩阵U
U=(I+J)/2,其中I为N×N维的单位矩阵,J是一个N×N维的反对角阵。
J = 0 0 . . . 0 1 0 0 . . . 1 0 . . . . . . . . . . . . . . . 1 0 . . . 0 0
6.根据权利要求1所述的检测方法,其中,所述预定阈值是采用蒙特卡罗仿真得到的。
7.根据权利要求1所述的检测方法,其中,所述距离扩展目标的标称导向向量是通过所述距离扩展目标的目标方向计算得到的。
8.一种基于混响协方差矩阵的距离扩展目标的检测装置,包括:回波数据获取单元,用于通过由声纳阵接收的回波获取一组回波数据,选取所述一组回波数据中的多个回波数据作为主数据,选取所述一组回波数据中除所述主数据外的多个回波数据作为辅数据,所述主数据和所述辅数据具有相同的混响协方差矩阵,其特征在于,
第一计算单元,用于根据所述辅数据计算采样协方差矩阵;
第二计算单元,用于分别计算所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息和混响协方差矩阵的最大似然估计的伪对称表示信息,其中,所述混响协方差矩阵的结构具有预先确定的伪对称特性;以及
第三计算单元,用于根据所述采样协方差矩阵的伪对称表示信息、所述主数据的伪对称表示信息、所述混响协方差矩阵的最大似然估计的伪对称表示信息和所述距离扩展目标的标称导向向量,计算检测统计量;
比较判断单元,用于将所述检测统计量与预定阈值进行比较,并且根据比较结果判断所述距离扩展目标是否存在。
9.一种声纳系统,包括如权利要求8所述的基于混响协方差矩阵的距离扩展目标的检测装置和声纳阵,所述声纳阵连接到所述回波数据获取单元,所述声纳阵用于发射探测声波并接收回波。
CN201210556458.6A 2012-12-19 2012-12-19 基于混响协方差矩阵的距离扩展目标的检测方法和装置 Active CN103033815B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210556458.6A CN103033815B (zh) 2012-12-19 2012-12-19 基于混响协方差矩阵的距离扩展目标的检测方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210556458.6A CN103033815B (zh) 2012-12-19 2012-12-19 基于混响协方差矩阵的距离扩展目标的检测方法和装置

Publications (2)

Publication Number Publication Date
CN103033815A CN103033815A (zh) 2013-04-10
CN103033815B true CN103033815B (zh) 2014-11-05

Family

ID=48020887

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210556458.6A Active CN103033815B (zh) 2012-12-19 2012-12-19 基于混响协方差矩阵的距离扩展目标的检测方法和装置

Country Status (1)

Country Link
CN (1) CN103033815B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237864B (zh) * 2014-09-29 2017-11-03 西安电子科技大学 基于匹配模糊函数的距离扩展目标检测方法
CN104977585B (zh) * 2015-06-11 2017-07-28 中国科学院声学研究所 一种鲁棒的运动声纳目标检测方法
CN104977579B (zh) * 2015-07-03 2017-07-28 中国科学院声学研究所 一种基于随机协方差矩阵的多亮点目标空时检测方法
CN105931648B (zh) * 2016-06-24 2019-05-03 百度在线网络技术(北京)有限公司 音频信号解混响方法和装置
CN106932766B (zh) * 2017-04-27 2019-07-09 中国人民解放军海军航空大学 基于变参数广义结构的距离扩展目标自适应检测方法
CN107561540B (zh) * 2017-08-13 2020-02-14 中国科学院声学研究所 基于混响对称谱特性的声纳泄露目标的检测方法及装置
CN107479050B (zh) * 2017-08-13 2020-04-14 中国科学院声学研究所 基于对称谱特性和次对称特性的目标检测方法和装置
CN107703495B (zh) * 2017-09-01 2019-10-01 中国科学院声学研究所 一种目标信号检测方法及系统
CN108646249B (zh) * 2018-05-11 2020-06-09 中国科学院声学研究所 一种适用于部分均匀混响背景的参数化泄露目标检测方法
CN108896967B (zh) * 2018-05-11 2020-11-17 清华大学 基于杂波协方差矩阵估计的距离扩展目标检测方法及装置
CN108919224B (zh) * 2018-07-26 2020-06-30 中国人民解放军海军航空大学 基于斜对称结构的宽带雷达目标自适应融合检测方法
RU2736097C1 (ru) * 2019-12-30 2020-11-11 Акционерное Общество "Концерн "Океанприбор" Способ обнаружения локального объекта на фоне распределенной помехи
RU2736567C1 (ru) * 2019-12-31 2020-11-18 Акционерное Общество "Концерн "Океанприбор" Способ обнаружения локального объекта на фоне распределенной помехи при бистатической гидролокации
CN112505665B (zh) * 2020-11-10 2023-06-23 中国科学院声学研究所 适用于部分均匀混响环境的空时自适应检测方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009110578A1 (ja) * 2008-03-03 2009-09-11 日本電信電話株式会社 残響除去装置、残響除去方法、残響除去プログラム、および記録媒体
JP2009212599A (ja) * 2008-02-29 2009-09-17 Nippon Telegr & Teleph Corp <Ntt> 残響除去装置とその方法と、そのプログラムと記録媒体

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5129794B2 (ja) * 2009-08-11 2013-01-30 日本電信電話株式会社 目的信号強調装置とその方法と、プログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009212599A (ja) * 2008-02-29 2009-09-17 Nippon Telegr & Teleph Corp <Ntt> 残響除去装置とその方法と、そのプログラムと記録媒体
WO2009110578A1 (ja) * 2008-03-03 2009-09-11 日本電信電話株式会社 残響除去装置、残響除去方法、残響除去プログラム、および記録媒体

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
ACOUSTIC SOURCE LOCATION IN A THREE-DIMENSIONAL SPACE USING CROSSPOWER SPECTRUM PHASE;Piergiorgio Svaizer et al.;《1997 IEEE》;20071231;231-234 *
Experimental demonstration of adaptive reverberation nulling using time reversal;H. C. Song et al.;《2005 Acoustical Society of America》;20050930;1381-1387 *
主动声呐混响抑制的空时相三维处理方法;艾锐峰等;《电子与信息学报》;20120531;第34卷(第5期);1263-1267 *
水中高速小目标被动检测模型及其应用;胡桥等;《鱼雷技术》;20120831;第20卷(第4期);261-267 *
混响空时协方差矩阵的两种计算方法比较与分析;吕维等;《鱼雷技术》;20120831;第20卷(第4期);251-256 *

Also Published As

Publication number Publication date
CN103033815A (zh) 2013-04-10

Similar Documents

Publication Publication Date Title
CN103033815B (zh) 基于混响协方差矩阵的距离扩展目标的检测方法和装置
Soukissian et al. Effects of different wind data sources in offshore wind power assessment
CN104898103B (zh) 基于多通道杂波图的低速目标检测方法
CN107085205B (zh) 基于杂波协方差矩阵结构信息的自适应检测方法
CN105425223B (zh) 广义帕累托杂波下稀疏距离扩展雷达目标的检测方法
Ganske et al. Identification of extreme storm surges with high-impact potential along the German North Sea coastline
CN105425230A (zh) 逆高斯纹理复合高斯杂波下距离扩展目标自适应检测方法
CN110515052B (zh) 一种基于时间反演的超宽带频域非等间隔采样目标检测方法
Tagle et al. A non‐Gaussian spatio‐temporal model for daily wind speeds based on a multi‐variate skew‐t distribution
CN105259541A (zh) 一种多站雷达抗有源欺骗式干扰的方法
CN109597046A (zh) 基于一维卷积神经网络的米波雷达doa估计方法
CN103777189A (zh) 一种基于信息几何复自回归模型的雷达弱目标检测方法
CN103729846A (zh) 基于不规则三角网的LiDAR点云数据边缘检测方法
Bernier et al. Mapping the return periods of extreme sea levels: allowing for short sea level records, seasonality, and climate change
CN116105697A (zh) 一种基于深度学习的星载gnss-r全球海面有效波高反演方法
CN110850420B (zh) 基于弹球损失的FisherSVM声纳信号判别方法
CN113536948A (zh) 一种渔船捕捞行为识别方法及装置
Cosoli et al. A real-time and offline quality control methodology for SeaSonde high-frequency radar currents
CN107748361A (zh) 基于截断杂波统计的sar图像双参数cfar检测方法
CN106093925A (zh) 一种多雷达联合检测方法
CN104407336B (zh) 一种方位敏感的目标电磁回波模拟方法
CN102034233B (zh) 一种基于contourlet变换的SAR图像波群参数检测方法
Li et al. Quantifying uncertainty in extreme values of design parameters with resampling techniques
CN104614597A (zh) 一种雷暴预警方法
CN107909595A (zh) 基于幅度与能量聚集度乘积的海陆杂波场景分割方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant