一种适用于混响边缘环境的水下多目标检测方法

1.本技术涉及声纳探测技术领域，尤其涉及一种适用于混响边缘环境的水下多目标检测方法。


背景技术：

2.海洋混响作为主动声纳的永恒背景干扰，对海洋探测装备性能发挥起着主要限制作用。对于浅海主动声纳探测，尤其当声纳载体具有一定运动速度时，平台接收混响谱在角度-多普勒域呈现明显的非线性映射关系，称为空时耦合性。水下空时自适应检测(stad)作为融合空时二维滤波和恒虚警处理的水声探测技术，是实现浅海混响抑制和目标检测的关键手段。
3.然而，受海底地形起伏、海水内波、界面起伏等非均匀介质和不规则散射面调制，混响能量随参考单元渐变或突变呈现显著非均匀性，又称为混响边缘环境。该环境特征导致运动声纳接收样本不再满足统计一致的要求，促使现有水下stad模型的干扰估计精度和感兴趣信号急剧衰减，大幅降低声纳探测距离并引发目标漏检。此外，实际的水下环境极为复杂，水下目标多以集群形式出现。如何实现水下stad技术在混响边缘环境下的多目标精确探测，在声纳探测领域具有重要的理论意义和应用价值。
4.传统水下stad模型对混响空时建模偏于理想化，忽略了浅海地形变化、界面起伏等因素导致的统计非均匀性，缺乏对混响边缘环境影响机理、统计特征和场景类别的充分认知，对混响边缘位置、参考单元数据能量分布等信息挖掘不足，无法满足实际应用需求。针对该问题，国外学者pia addabbo等利用期望最大化(em)算法提出了混响边缘背景下的回波样本分类方案，实现了对不同参考单元样本所属能量区域分类以及边缘位置查找，达到提升非均匀混响抑制的目的。
5.然而，混响抑制仅是声纳探测的中间过程，而非最终目标。现有非均匀环境下的协方差矩阵分类方案仅考虑了对混响边缘位置、能量区划分的实现，未考虑由多目标后向散射结构化回波造成的辅助数据污染，不具备不同参考单元上的多个点目标探测能力，从而不能满足实际声纳探测需求。


技术实现要素：

6.为解决浅海混响边缘环境中的运动声纳多目标检测困难问题，本发明目的在于提供一种适用于混响边缘环境的水下多目标检测方法，切实提高非均匀混响背景下的水下stad性能。
7.所述一种适用于混响边缘环境的水下多目标检测方法，包括：
8.建立存在多点目标的混响边缘统计分布及离散时间信号模型；基于所述模型，给出水下多目标stad的二元假设检验问题；
9.利用em算法求解所述模型，完成对背景未知参数的强化估计；基于所述背景未知参量的强化估计实现对不同参考单元非均匀混响样本所属能量区域的分类，边缘位置查找
以及所述多目标所在单元的位置估计；
10.将所述背景未知参量的强化估计结果应用于所述水下多目标stad的似然比检验统计量求解，完成对混响边缘环境下的水下多目标存在性的恒虚警检测。
11.优选的，所述建立存在多点目标的混响边缘统计分布及离散时间信号模型，在主动声纳回波经信号预处理后，得到k个参考单元样本矢量表示复数域，n为系统空时维度。
12.将声纳回波划分为l个非均匀混响能量区，单个能量区内混响区间协方差矩阵表示为其中，l＝1，...，l为能量区类别标志，k
l
为单个能量区内参考单元样本数，满足
13.考虑非均匀混响样本服从n维0均值、协方差矩阵为m
l
的多元复高斯分布，将所述水下多目标stad的二元假设检验问题描述如下：
[0014][0015]
式中h0和h1分别表示无目标和有目标假设。
[0016]
优选的，所述h1的假设，首先利用利用声纳参考单元样本矩阵z＝[z1，...，zk]的联合对数似然函数求解jensen不等式得到em算法e步计算结果：
[0017][0018]
式中(h1)表示em算法第h-1次迭代的计算结果，l0 l和l1 l分别表示第l个混响能量区中无目标和有目标的类别标识，为h1假设下第h-1次迭代的参数估计集合，为h0假设下第h-1
次迭代的参数估计集合，e-h(s)
为基于模型阶数选择准则的惩罚函数，h(s)是与未知参量数量有关的惩罚因子，为第h-1次迭代中参考单元样本zk属于ls l类别的概率，为参考单元样本zk的条件概率密度函数。
[0019]
其次利用最大似然估计方法求解em算法m步计算结果，得到第h次迭代中未知参量和估计：
[0020][0021]
em迭代结束后，取各参考单元样本的最大值为最终分类结果。
[0022]
优选的，所述h0的假设，首先利用参考单元样本矩阵z＝[z1，...，zk]的联合对数似然函数求解jensen不等式得到em算法的e步计算结果：
[0023][0024]
其次利用最大似然估计方法求解em算法m步计算结果，得到第h次迭代中未知参量和估计：
[0025][0026]
em迭代结束后，取各参考单元样本的最大值为最终分类结果。
[0027]
优选的，将所述背景未知参量的强化估计结果应用于所述水下多目标stad的似然比检验统计量求解，得到一种混响边缘环境下的水下多目标stad方法：
[0028][0029]
本发明提供一种适用于混响边缘环境的水下多目标检测方法，将不同混响能量区、点目标存在时各参考单元统计特征变化规律用于em算法，具备准确的样本分类性能。由于样本分类及目标检测均是在各参考单元上进行，实现目标检测的同时还可估算出各目标在参考单元的位置。相比于传统技术，本发明利用不同参考单元混响协方差矩阵的潜在统计信息实现混响所属能量区分类，同时基于分类结果指导完成未知参量强化估计以及水下stad似然比检验求解，实现混响边缘位置查找以及混响能量区划分，达到提升多个点目标检测性能的目的。
附图说明
[0030]
为了更简单说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0031]
图1为本技术实施例中适用于混响边缘环境的水下多目标检测方法的流程示意图；
[0032]
图2为本技术实施例中采用蒙特-卡罗仿真方法得到的样本分类性能图；
[0033]
图3为本技术实施例中采用蒙特-卡罗仿真方法得到的目标检测概率随sinr变化曲线图。
具体实施方式
[0034]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例，本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0035]
本发明将混响边缘环境下的多目标检测问题转换为水下stad二元假设检验求解，
包括有目标假设和无目标假设两部分。为提高检测精度，首先利用混响边缘统计分布建立非均匀混响空时模型，基于不同参考单元混响数据的统计特征变化规律，引入代表混响能量区域类别的隐藏变量。采用em算法求解模型，完成对背景中混响协方差矩阵、目标幅值等未知参量的强化估计，以及对不同参考单元非均匀混响样本所属能量区域的分类、边缘位置查找以及所述多目标所在单元的位置估计，并将背景未知参量的强化估计结果用于水下stad似然比检测统计量求解，完成对混响边缘环境下的水下多目标存在性的恒虚警检测。
[0036]
图1中，本发明提出的一种适用于混响边缘环境的水下多目标检测方法的步骤如下：
[0037]
步骤s110：建立存在多点目标的混响边缘统计分布及离散时间信号模型，给出水下多目标stad的二元假设检验问题。
[0038]
首先介绍多个点目标存在时混响边缘统计分布及离散时间信号模型，在此基础上，给出水下多目标stad的二元假设检验问题。
[0039]
假设主动声纳回波经信号预处理后，得到k个参考单元样本矢量表示复数域，n为系统空时维度。受海底地形变化等因素影响，声纳回波划分为l个非均匀混响能量区，不同混响区间协方差矩阵特性不一致，单个能量区内混响强度变化相对平稳，表示为为能量区类别标志，k
l
为单个能量区内参考单元样本数，满足
[0040]
考虑非均匀混响样本服从n维0均值、协方差矩阵为m
l
的多元复高斯分布：当有目标存在时，声纳回波样本表示为t表示有目标存在的参考单元位置下标，t是目标总数，αi是第i个目标的未知幅值因子，v是目标归一化空时导向矢量。基于上述模型，水下多目标stad的二元假设检验问题描述如下：
[0041][0042]
式中h0和h1分别表示无目标和有目标假设。
[0043]
步骤s120：利用em算法求解模型，完成对背景未知参数的强化估计，实现对不同参考单元非均匀混响样本所属能量区域的分类、边缘位置查找以及多目标所在单元的位置估计。
[0044]
基于上式中模型描述，引入代表混响能量区域类别的隐藏随机变量ck，k＝1，...，k，其中每个ck均与zk相对应，p(ck＝l)＝p
l
表示参考单元样本zk属于能量区l的概率，又称ck的概率质量函数。下面利用em算法求解上式中模型，联合实现非均匀混响样本分类、背景未知参量估计和多目标似然比检测。
[0045]
一方面，基于h1的假设，首先利用声纳参考单元样本矩阵z＝[z1，...，zk]的联合对数似然函数求解jensen不等式得到em算法e步计算结果：
[0046][0047]
式中(h1)表示em算法第h-1次迭代的计算结果，l0 l和l1 l分别表示第l个混响能
量区中无目标和有目标的类别标识，为h1假设下第h-1次迭代的参数估计集合，为h0假设下第h-1次迭代的参数估计集合，e-h(s)
为基于模型阶数选择准则的惩罚函数，h(s)是与未知参量数量有关的惩罚因子，为第h-1次迭代中参考单元样本zk属于ls l类别的概率，为参考单元样本zk的条件概率密度函数，表示为：
[0048][0049][0050]
其中，|
·
|表示矩阵行列式，exp{
·
}表示指数运算，tr[
·
]表示矩阵求迹，(
·
)h表示共轭转置。
[0051]
其次利用最大似然估计方法求解em算法m步计算结果，得到第h次迭代中未知参量和估计：
[0052][0053]
需说明的是，由于和估计过程中存在参数嵌套问题，因此采用循环迭代方式得到它们的次优解，步骤如下：
[0054]
(1)已知，估计
[0055][0056]
(2)已知，估计
[0057][0058]
(3)反复循环步骤(1)～(2)，直至满足迭代停止条件。
[0059]
em迭代结束后，取各参考单元样本的最大值为最终分类结果。
[0060]
基于h1的假设，可实现对不同参考单元非均匀混响样本的能量区域划分，边缘位置查找以及各目标所在单元位置估计，其中对目标估计的实现也可以作为步骤s130中多目标检测的辅助验证
[0061]
另一方面，基于h0的假设，同样利用参考单元样本矩阵z＝[z1，...，zk]的联合对数似然函数求解jensen不等式得到em算法的e步计算结果：
[0062]
[0063]
其中
[0064]
其次利用最大似然估计方法求解em算法m步计算结果，得到第h次迭代中未知参量和估计：
[0065][0066][0067]
em迭代结束后，取各参考单元样本的最大值为最终分类结果。
[0068]
基于h0的假设，可实现对不同参考单元非均匀混响样本的能量区域划分，边缘位置查找。
[0069]
步骤s130：将背景未知参量的强化估计结果应用于水下多目标stad的似然比检验统计量求解，完成对混响边缘环境下的水下多目标存在性的恒虚警检测，得到一种混响边缘环境下的水下多目标stad方法：
[0070][0071]
在一个实施例中，对本发明提出的一种适用于混响边缘环境的水下多目标检测方法做仿真分析，假设声纳阵元数n＝8，接收参考单元样本数k＝64，混响回波被划分为l＝2个非均匀混响能量区，各能量区样本数k1＝k2＝32，混响协方差矩阵其中mc＝ρ
|i-j|
，i，j＝1，...，n，ρ＝0.9为一步滞后相关系数，分别为第1和第2个混响能量区功率，两个目标分别位于第15和38个参考单元处，目标波达角度均匀0
°
，信干噪比下面采用蒙特-卡罗仿真方法对所提出方法的非均匀样本分类、目标检测性能加以分析，设定对检测概率pd的独立实验次数为103次，虚警概率p
fa
和蒙特-卡罗试验次数分别为10-2
和100/p
fa
，em循环最大迭代次数为15，总分类类别数为4，具体包括：
[0072]
类别1：参考单元位于混响能量区1，无目标存在；
[0073]
类别2：参考单元位于混响能量区2，无目标存在；
[0074]
类别3：参考单元位于混响能量区1，有一个目标存在；
[0075]
类别4：参考单元位于混响能量区2，有一个目标存在。
[0076]
图2给出各参考单元分别在sinr＝25db和35db时的样本分类性能，可以看出所提出方法的分类结果与场景设置的真实结果基本一致，可实现准确的混响边缘位置查找、混响能量区分类以及目标所在参考单元位置估计。图3给出pd随sinr的变化曲线，可看出在sinr＝20db时，所提出方法的目标检测概率达到1，具有优越的目标检测性能。
[0077]
针对浅海混响边缘环境下的运动声纳多目标检测问题，现有方法未将多目标回波纳入混响样本分类方案，不具备不同参考单元上的多个点目标探测能力，无法满足实际声
纳探测需求。本发明提出一种适用于混响边缘环境的水下多目标检测方法，利用目标存在及不同混响能量区的参考单元统计特征变化规律实现样本分类，同时基于分类结果指导完成未知参量强化估计以及水下stad似然比检验统计量求解，实现混响边缘位置查找、混响能量区划分的同时，具备优越的多个点目标检测能力。
[0078]
需要说明的是，在此提供的方法不与任何特定计算机、虚拟装置或者其它设备固有相关。各种通用装置也可以与基于在此的示教一起使用。根据上面的描述，构造这类装置所要求的结构是显而易见的。此外，本发明也不针对任何特定的编程语言。应当明白，可以利用各种编程语言实现在此描述的本发明内容，并且上面对特定语言、系统功能模块的调用所做的描述仅仅是为了披露发明的最佳实施方式。
[0079]
在此处所提供的说明书中，说明了大量的具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下完成实现。在一些示例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
[0080]
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若对本发明的这些修改和变型属于本发明权利要去及其同等技术的范围之内，则本发明也意图包含这些改动和变型在内。