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

1.本发明属于领域，具体涉及一种基于单水听器接收信号估计运动声源径向速度与深度的方法。

背景技术：

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

技术实现要素：

5.本发明提供一种基于单水听器接收信号估计运动声源径向速度与深度的方法，用以实现声源速度与深度的估计。
6.本发明通过以下技术方案实现：
7.一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述方法具体包括以下步骤：
8.步骤一、根据单水听器接收信号的声强谱图，提取其干涉结构在频率-模态多普勒频移域上的特征分布；
9.步骤二、根据步骤一的频率-模态多普勒频移域特征分布，建立相应的匹配目标函数；
10.步骤三、根据步骤二的目标函数，采用模拟退火方法提取声源径向速度-深度域上目标函数最大值的坐标，以此作为声源速度与深度的估计值。
11.一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤一具体为，浅海水平无关波导远场条件下的声场可写为一组简正波累加的形式：
[0012][0013]
其中，f为角频率，k
rn
(f)是第n阶简正波的本征值，ψn(z)是第n阶简正波的本征函数，zs和zr分别是声源深度和接收水听器深度，r是声源与接收水听器之间的距离，ρ(zs)是声源深度处的海水密度；
[0014]
分析声场的干涉结构时，式(1)简写为：
[0015][0016]
基于式子(2)，远距离声场强度写为如式子所示的形式：
[0017][0018]
其中本征值之差k
rm-k
rn
是与频率有关的物理量，因此exp[
·
]部分形成了距离r 与频率f上的干涉结构。
[0019]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，将式子(3) 写为各阶简正波声场强度累加的形式，则各阶简正波干涉声强的表达式为：
[0020][0021]
假定声源处于径向运动，不同时间的接收距离写为：
[0022]
r＝r0+vrt
ꢀꢀ
(5)
[0023]
则给出各阶简正波声强在接收场时域上的表达式：
[0024][0025]
基于式子的干涉结构，采用傅里叶谱估计的方法，得到时频干涉结构在频率
ꢀ‑
模态多普勒频移域上的特征分布：
[0026][0027]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤二具体为，根据已知的海洋环境参数，利用kraken简正波程序，计算宽频带条件下不同深度声源在接收深度上激发产生的声场结果p
rplc
(zs,zr,r,f)，由此得到匹配场的声场强度i
rplc
(zs,zr,r,f)，其中接收深度zr与频率f是已知的；
[0028]
根据f(f,ν)对径向速度vr、声源深度zs的敏感特性以及对初始距离r0的不敏感特性，最终将匹配场声强写为如式子(8)所示的形式：
[0029]irplc
(zs,zr,r,f)＝i
rplc
(zr,t,f；zs,vr)
ꢀꢀ
(8)
[0030]
其中zs与vr是待匹配的物理量。
[0031]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，基于式子 (7)的方法，将式子(8)中获得的时间-频率域上的匹配场声强干涉结构进行时域上的傅里叶分析，得到匹配场的频率-模态多普勒频移域输出结果：
[0032][0033]
设定观测场信号处理后的频率-模态多普勒频移域输出结果为f
obs
(ν,f)；在匹配处理中使用的是f
obs
(ν,f)与f
rplc
(ν,f；zs,vr)的模值，根据最小二乘原理，设计目标函数如式子所示：
[0034][0035]
其中nf与n
ν
分别是f轴与ν轴上的数据点数；最终将式子(11)中的目标函数取分贝处理，如式子所示：
[0036][0037]
其中，j
db
(zs,vr)为目标函数取分贝处理结果，j(zs,vr)为目标函数值。
[0038]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤三具体为，目标函数构造在(zs,vr)上的模糊平面，设定其中的最大值点坐标其表达式为：
[0039][0040]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述模拟退火方法为自适应单纯形模拟退火算法。
[0041]
本发明的有益效果是：
[0042]
本发明根据声强的时频干涉谱图，给出其时域干涉结构的傅里叶谱估计，得到接收信号在频率-模态多普勒频移(f.ν)域上的能量分布特征。在此基础上分析不同声源径向速度与声源深度对频率-模态多普勒频移域上能量分布的影响，证明以此为匹配场特征实现声源速度-深度估计的可行性，并由此设计目标函数。随后采用了模拟退火算法在待估计的声源径向速度-深度域上快速寻目标函数的最大值坐标，以此作为声源速度与深度的估计值。
附图说明
[0043]
图1是本发明的目标函数的敏感性分析图，其中，(a)径向速度的敏感性分析图，(b)声源深度的敏感性分析图，(c)初始距离的敏感性分析图。
[0044]
图2是本发明的接收信号声强lofar谱图，其中，(a)声源接近水听器算例接收信号声强lofar谱图，(b)声源远离水听器算例接收信号声强lofar 谱图。
[0045]
图3是本发明时频干涉结构在频率-模态多普勒频移域上的特征分布，其中， (a)声源接近水听器算例接收信号声强在频率-模态多普勒频移域上的特征分布， (b)声源远离水听器算例接收信号声强在频率-模态多普勒频移域上的特征分布。
[0046]
图4是本发明的目标函数在声源径向速度-声源深度域上的模糊函数，其中， (a)声源接近水听器算例目标函数在声源径向速度-声源深度域上的模糊函数， (b)声源远离水听器算例目标函数在声源径向速度-声源深度域上的模糊函数。
[0047]
图5是本发明的模拟退火算法迭代后目标函数变化示意图，其中，(a)声源接近水听器算例模拟退火算法迭代后目标函数变化，(b)声源远离水听器算例模拟退火算法迭代后目标函数变化。图6是本发明的方法流程图。
具体实施方式
[0048]
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0049]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述方法具体包括以下步骤：
[0050]
步骤一、根据单水听器接收信号的声强谱图，提取其干涉结构在频率-模态多普勒频移域上的特征分布；
[0051]
步骤二、根据步骤一的频率-模态多普勒频移域特征分布，建立相应的匹配目标函数；
[0052]
步骤三、根据步骤二的目标函数，采用模拟退火方法提取声源径向速度-深度域上目标函数最大值的坐标，以此作为声源速度与深度的估计值。
[0053]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤一具体为，浅海水平无关波导远场条件下的声场可写为一组简正波累加的形式：
[0054][0055]
其中，f为角频率，k
rn
(f)是第n阶简正波的本征值，ψn(z)是第n阶简正波的本征函数，zs和zr分别是声源深度和接收水听器深度，r是声源与接收水听器之间的距离，ρ(zs)是声源深度处的海水密度；
[0056]
分析声场的干涉结构时，式(1)简写为：
[0057][0058]
基于式子(2)，远距离声场强度写为如式子所示的形式：
[0059][0060]
其中本征值之差k
rm-k
rn
是与频率有关的物理量，因此exp[
·
]部分形成了距离r 与频率f上的干涉结构。
[0061]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，将式子(3) 写为各阶简正波声场强度累加的形式，则各阶简正波干涉声强的表达式为：
[0062][0063]
假定声源处于径向运动，不同时间的接收距离写为：
[0064]
r＝r0+vrt
ꢀꢀ
(5)
[0065]
则给出各阶简正波声强在接收场时域上的表达式：
[0066][0067]
基于式子的干涉结构，采用傅里叶谱估计的方法，得到时频干涉结构在频率
ꢀ‑
模态多普勒频移域上的特征分布：
[0068][0069]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤二具体为，根据已知的海洋环境参数，利用kraken简正波程序，计算宽频带条件下不同深度声源在接收深度上激发产生的声场结果p
rplc
(zs,zr,r,f)，由此得到匹配场的声场强度i
rplc
(zs,zr,r,f)，其中接收深度zr与频率f是已知的；
[0070]
根据f(f,ν)对径向速度vr、声源深度zs的敏感特性以及对初始距离r0的不敏感特性，最终将匹配场声强写为如式子(8)所示的形式：
[0071]irplc
(zs,zr,r,f)＝i
rplc
(zr,t,f；zs,vr)
ꢀꢀ
(8)
[0072]
其中zs与vr是待匹配的物理量。
[0073]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，基于式子 (7)的方法，将式子(8)中获得的时间-频率域上的匹配场声强干涉结构进行时域上的傅里叶分析，得到匹配场的频率-模态多普勒频移域输出结果：
[0074][0075]
设定观测场信号处理后的频率-模态多普勒频移域输出结果为f
obs
(ν,f)；在匹配处理中使用的是f
obs
(ν,f)与f
rplc
(ν,f；zs,vr)的模值，根据最小二乘原理，设计目标函数如式子所示：
[0076][0077]
其中nf与n
ν
分别是f轴与ν轴上的数据点数；最终将式子(11)中的目标函数取分贝处理，如式子所示：
[0078][0079]
其中，j
db
(zs,vr)为目标函数取分贝处理结果，j(zs,vr)为目标函数。
[0080]
目标函数的敏感性分析如图1所示。
[0081]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述步骤三具体为，目标函数构造在(zs,vr)上的模糊平面，设定其中的最大值点坐标其表达式为：
[0082]
[0083]
一种基于单水听器接收信号估计运动声源径向速度与深度的方法，所述模拟退火方法为自适应单纯形模拟退火算法。
[0084]
设定海水层声速保持在1500m/s，海水层深度50m；海底半无限空间中纵波波速c
p
为1800m/s，，密度ρ2为1.5g/cm3，纵波传播损失α
p
为0.4db/λ。设定声源深度为10m，径向速度为4m/s。设定两个算例：1.初始距离为10km，声源接近水听器；2.初始距离为6.8km，声源远离水听器。设定水听器接收深度为50m 于海底。由此得到的接收信号声强谱图如图2所示。
[0085]
根据式子的傅里叶谱估计，可以得到如图2的时频干涉结构在频率-模态多普勒频移域上的特征分布，如图3所示：
[0086]
式子计算的目标函数在声源径向速度-声源深度域上的模糊函数如图4所示：
[0087]
由于遍历得到的模糊函数在计算较大范围时速度较慢，且估计精度较低，因此采用模拟退火方法进行估计，其迭代后目标函数变化如图5所示：
[0088]
根据模拟退火算法估计得到的声源径向速度与深度结果如表1所示：
[0089]
表1声源径向速度与深度估计结果
[0090][0091][0092]
以上对本发明提出的基于干涉结构在频率-模态多普勒频移域上特征的声源径向速度-深度估计方法，进行了详尽的介绍，本文中采用了不同运动状态的声源算例对本发明的原理及实施方式进行了阐述。

技术特征：

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
,r,f)，由此得到匹配场的声场强度i
rplc
(z
s
,z
r
,r,f)，其中接收深度z
r
与频率f是已知的；根据f(f,ν)对径向速度v
r
、声源深度z
s
的敏感特性以及对初始距离r0的不敏感特性，最终将匹配场声强写为如式子(8)所示的形式：i
rplc
(z
s
,z
r
,r,f)＝i
rplc
(z
r
,t,f；z
s
,v
r
)
ꢀꢀꢀꢀꢀ
(8)其中z
s
与v
r
是待匹配的物理量。5.根据权利要求4所述一种基于单水听器接收信号估计运动声源径向速度与深度的方法，其特征在于，基于式子(7)的方法，将式子(8)中获得的时间-频率域上的匹配场声强干涉结构进行时域上的傅里叶分析，得到匹配场的频率-模态多普勒频移域输出结果：设定观测场信号处理后的频率-模态多普勒频移域输出结果为f
obs
(ν,f)；在匹配处理中使用的是f
obs
(ν,f)与f
rplc
(ν,f；z
s
,v
r
)的模值，根据最小二乘原理，设计目标函数如式子所示：其中n
f
与n
ν
分别是f轴与ν轴上的数据点数；最终将式子(10)中的目标函数取分贝处理，如式子所示：其中，j
db
()为目标函数取分贝处理结果，j为目标函数值。6.根据权利要求5所述一种基于单水听器接收信号估计运动声源径向速度与深度的方法，其特征在于，所述步骤三具体为，目标函数构造在(z
s
,v
r
)上的模糊平面，设定其中的最大值点坐标其表达式为：7.根据权利要求1所述一种基于单水听器接收信号估计运动声源径向速度与深度的方法，其特征在于，所述模拟退火方法为自适应单纯形模拟退火算法。

技术总结

本发明提供一种基于单水听器接收信号估计运动声源径向速度与深度的方法。步骤一、根据单水听器接收信号的声强谱图，提取其干涉结构在频率-模态多普勒频移域上的特征分布；步骤二、根据步骤一的频率-模态多普勒频移域特征分布，建立相应的匹配目标函数；步骤三、根据步骤二的目标函数，采用模拟退火方法提取声源径向速度-深度域上目标函数最大值的坐标，以此作为声源速度与深度的估计值。本发明用以实现声源速度与深度的估计。现声源速度与深度的估计。现声源速度与深度的估计。