一种基于单水听器接收信号估计运动声源径向速度与深度的方法技术

本发明专利技术提供一种基于单水听器接收信号估计运动声源径向速度与深度的方法。步骤一、根据单水听器接收信号的声强谱图，提取其干涉结构在频率

【技术实现步骤摘要】
一种基于单水听器接收信号估计运动声源径向速度与深度的方法


[0001]本专利技术属于领域，具体涉及一种基于单水听器接收信号估计运动声源径向速 度与深度的方法。

技术介绍

[0002]辐射噪声的分析是侦测水下声源的有效手段之一，减振降噪技术的主要进步 方向是抑制辐射噪声中的线谱特征，而低频连续谱特征尚未受到有效的抑制。近 年来，基于宽频带信号的处理方法得到了学者们的重视。
[0003]俄罗斯学者Kuznetsov利用单水听器在时间
‑
频率域上的声强干涉结构，实 现了运动声源径向速度的估计。在进一步的研究中，他们基于波导不变量理论与 二维傅里叶变换方法，提出了一种利用时频干涉结构的声源运动参数估计算法。
[0004]然而，Kuznetsov等人的方法需要具体确定时间
‑
频率域声强干涉结构中起 到主要贡献的简正波模态，这一特征的先验信息包括海洋环境参数与声源深度。 若声源的深度不是先验信息，则会严重影响了这一算法在处理时频干涉结构上的 应用。但是，Nicolas等人利用水平阵提取了接收信号在频率
‑
波数域上的频散 特征，这一特征包含了声源的深度信息。T.C.Yang利用垂直阵获得了不同接收 深度信号的模态特征，通过模态匹配也获得了声源深度的准确估计。由此可见， 利用接收信号的简正波特征是能够实现声源深度的估计的。

技术实现思路

[0005]本专利技术提供一种基于单水听器接收信号估计运动声源径向速度与深度的方 法，用以实现声源速度与深度的估计。
[0006]本专利技术通过以下技术方案实现：
[0007]一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述方法 具体包括以下步骤：
[0008]步骤一、根据单水听器接收信号的声强谱图，提取其干涉结构在频率
‑
模态 多普勒频移域上的特征分布；
[0009]步骤二、根据步骤一的频率
‑
模态多普勒频移域特征分布，建立相应的匹配 目标函数；
[0010]步骤三、根据步骤二的目标函数，采用模拟退火方法提取声源径向速度
‑
深 度域上目标函数最大值的坐标，以此作为声源速度与深度的估计值。
[0011]一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤 一具体为，浅海水平无关波导远场条件下的声场可写为一组简正波累加的形式：
[0012][0013]其中，f为角频率，k
rn
(f)是第n阶简正波的本征值，ψ
n
(z)是第n阶简正 波的本征函数，z
s
和z
r
分别是声源深度和接收水听器深度，r是声源与接收水听 器之间的距离，ρ(z
s
)是声源深度处的海水密度；
[0014]分析声场的干涉结构时，式(1)简写为：
[0015][0016]基于式子(2)，远距离声场强度写为如式子所示的形式：
[0017][0018]其中本征值之差k
rm
‑
k
rn
是与频率有关的物理量，因此exp[
·
]部分形成了距离r 与频率f上的干涉结构。
[0019]一种基于单水听器接收信号估计运动声源径向速度与深度的方法，将式子(3) 写为各阶简正波声场强度累加的形式，则各阶简正波干涉声强的表达式为：
[0020][0021]假定声源处于径向运动，不同时间的接收距离写为：
[0022]r＝r0+v
r
t
ꢀꢀ
(5)
[0023]则给出各阶简正波声强在接收场时域上的表达式：
[0024][0025]基于式子的干涉结构，采用傅里叶谱估计的方法，得到时频干涉结构在频率
ꢀ‑
模态多普勒频移域上的特征分布：
[0026][0027]一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤 二具体为，根据已知的海洋环境参数，利用Kraken简正波程序，计算宽频带条 件下不同深度声源在接收深度上激发产生的声场结果p
rplc
(z
s
,z
r
,r,f)，由此得到 匹配场的声场强度I
rplc
(z
s
,z
r
,r,f)，其中接收深度z
r
与频率f是已知的；
[0028]根据F(f,ν)对径向速度v
r
、声源深度z
s
的敏感特性以及对初始距离r0的不 敏感特性，最终将匹配场声强写为如式子(8)所示的形式：
[0029]I
rplc
(z
s
,z
r
,r,f)＝I
rplc
(z
r
,t,f；z
s
,v
r
)
ꢀꢀ
(8)
[0030]其中z
s
与v
r
是待匹配的物理量。
[0031]一种基于单水听器接收信号估计运动声源径向速度与深度的方法，基于式子 (7)的方法，将式子(8)中获得的时间
‑
频率域上的匹配场声强干涉结构进行 时域上的傅里叶分析，得到匹配场的频率
‑
模态多普勒频移域输出结果：
[0032][0033]设定观测场信号处理后的频率
‑
模态多普勒频移域输出结果为F
obs
(ν,f)；在 匹配处理中使用的是F
obs
(ν,f)与F
rplc
(ν,f；z
s
,v
r
)的模值，根据最小二乘原理， 设计目标函数如式子所示：
[0034][0035]其中N
f
与N
ν
分别是f轴与ν轴上的数据点数；最终将式子(11)中的目标函数 取分贝处理，如式子所示：
[0036][0037]其中，J
dB
(z
s
,v
r
)为目标函数取分贝处理结果，J(z
s
,v
r
)为目标函数值。
[0038]一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤 三具体为，目标函数构造在(z
s
,v
r
)上的模糊平面，设定其中的最大值点坐标 其表达式为：
[0039][0040]一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述模拟 退火方法为自适应单纯形模拟退火算法。
[0041]本专利技术的有益效果是：
[0042]本专利技术根据声强的时频干涉谱图，给出其时域干涉结构的傅里叶谱估计，得 到接收信号在频率
‑
模态多普勒频移(f.ν)域上的能量分布特征。在此基础上分析 本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于单水听器接收信号估计运动声源径向速度与深度的方法，其特征在于，所述方法具体包括以下步骤：步骤一、根据单水听器接收信号的声强谱图，提取其干涉结构在频率
‑
模态多普勒频移域上的特征分布；步骤二、根据步骤一的频率
‑
模态多普勒频移域特征分布，建立相应的匹配目标函数；步骤三、根据步骤二的目标函数，采用模拟退火方法提取声源径向速度
‑
深度域上目标函数最大值的坐标，以此作为声源速度与深度的估计值。2.根据权利要求1所述一种基于单水听器接收信号估计运动声源径向速度与深度的方法，其特征在于，所述步骤一具体为，浅海水平无关波导远场条件下的声场可写为一组简正波累加的形式：其中，f为角频率，k
rn
(f)是第n阶简正波的本征值，ψ
n
(z)是第n阶简正波的本征函数，z
s
和z
r
分别是声源深度和接收水听器深度，r是声源与接收水听器之间的距离，ρ(z
s
)是声源深度处的海水密度；分析声场的干涉结构时，式(1)简写为：基于式子(2)，远距离声场强度写为如式子所示的形式：其中本征值之差k
rm
‑
k
rn
是与频率有关的物理量，因此exp[
·
]部分形成了距离r与频率f上的干涉结构。3.根据权利要求2所述一种基于单水听器接收信号估计运动声源径向速度与深度的方法，其特征在于，将式子(3)写为各阶简正波声场强度累加的形式，则各阶简正波干涉声强的表达式为：假定声源处于径向运动，不同时间的接收距离写为：r＝r0+v
r
t
ꢀꢀꢀꢀꢀ
(5)则给出各阶简正波声强在接收场时域上的表达式：基于式子的干涉结构，采用傅里叶谱估计的方法，得到时频干涉结构在频率
‑
模态多普勒频移域上的特征分布：
4.根据权利要求1所述一种基于单水听器接收信号估计运动声源径向速度与深度的方法，其特征在于，所述步骤二具体为，根据已知的海洋环境参数，利用Kraken简正波程序，计算宽频带条件下不同深度声源在接收深度上激发产生的声场结果p
rplc
(z
s
,z
r
...

【专利技术属性】
技术研发人员：张明辉，叶扬，朴胜春，张海刚，付金山，宋扬，郭俊媛，龚李佳，李楠松，王笑寒，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：