CN101335574A - 多途信道中的声屏蔽方法 - Google Patents

多途信道中的声屏蔽方法 Download PDF

Info

Publication number
CN101335574A
CN101335574A CNA2008101368231A CN200810136823A CN101335574A CN 101335574 A CN101335574 A CN 101335574A CN A2008101368231 A CNA2008101368231 A CN A2008101368231A CN 200810136823 A CN200810136823 A CN 200810136823A CN 101335574 A CN101335574 A CN 101335574A
Authority
CN
China
Prior art keywords
sound
shielding
focusing
power
output
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.)
Granted
Application number
CNA2008101368231A
Other languages
English (en)
Other versions
CN101335574B (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN200810136823A priority Critical patent/CN101335574B/zh
Publication of CN101335574A publication Critical patent/CN101335574A/zh
Application granted granted Critical
Publication of CN101335574B publication Critical patent/CN101335574B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Abstract

本发明提供的是一种多途信道中的声屏蔽方法。由声速仪测量水的声速;水平接收阵接收声场中的多个点源的辐射信号,放大滤波,然后进行数据采集;实现球面波声场中阵的声聚焦,得到的权值即为最佳聚焦权;在最佳聚焦权的基础上,增加屏蔽因子,得到最佳聚焦屏蔽权;进行反傅立叶变换可得到输出波形,也可进行频域积分得到输出功率;后续相干处理或输出经声屏蔽、聚焦后的功率。本发明以声屏蔽原理、方法为基础,适用于水声工程中抗多点源相干干扰,并且工程上可采用水平阵,更适合装载在运动平台上或海底固定使适用。

Description

多途信道中的声屏蔽方法
(一)技术领域
本发明涉及利用信号处理技术在水声相干多途信道中分辨多目标,特别是涉及同一方向上多目标的分辨,即探测同一方向上不同距离处的任意声源的方法。
(二)背景技术
基阵视野中若有多个声源,它们将相互干扰而降低对其中任何一个目标的探测能力,称之“多目标相互干扰”。信号处理的理想决策是对有兴趣探测的点源(目标)实施聚焦,对其余声源实施屏蔽。在平面波条件下,既不能实现聚焦,也不能实现屏蔽,所能实现的最优处理是阵处理技术,即将波束对准目标,而将指向性零点对准其余的(干扰)点源。波束零点(或极小值)抑制了多目标干扰,但零点方向也是探测盲区。也就是说,传统的阵处理技术无法分辨同一方向上的两个点源。
在平面波声场中,既不能实现声屏蔽,也不能实现声聚焦,首选的信号处理技术是(传统的)阵处理技术,该技术无法分辨同一方向上的两个点源。
在球面波条件下,可采用聚焦波束形成技术(参见:聚焦波束形成声图测量原理研究,《声学学报》,2007;32(4):356-361),在中、近程实现声聚焦。各阵元输出按球面波传播时延作时延差补偿并加权求和,就可得到聚焦输出。聚焦和屏蔽均是基于声场相干性的利用。若按球面波时延并使各通道反相输出进行补偿,如此,各通道反相相加,阵必在屏蔽点输出极小值,从而实现“声屏蔽”。
时间反转镜(Time Reversal Mirror)(参见:Phase conjugation in the ocean:Experimental demonstration of an acoustic time-reversal mirror.《J.Acoust.soc.Am》,1998;103(1):25-40;用时间反转法在水下波导介质中实现自适应聚焦的研究,《声学学报》,2002,27(6):541-548;矢量反转镜时空滤波技术研究,《声学学报》,2005,30(3):271-278。)信号处理能在多途信道(或称波导)中实现声聚焦和声屏蔽。M.A.Kuperman用垂直阵时间反转镜抑制混响并增强回波(参见:Environmentally adaptive reverberation nulling using a time reversal mirror.《Acoust.soc.Am》,August 2004;116(2):762-768),称之为“Reverberation Nulling”。海试证明方法是科学的,但要求在目标处设置一PS(Probe Source)源或采用物理时反镜,这不实用,无法在位置未知的非合作目标处加一PS源。但他们的研究有重要的科学价值。
(三)发明内容
本发明的目的在于提供一种不仅用于设置零点抗混响,而且用于设置屏蔽区抗相干噪声源干扰,并且可应用于水平阵(或单水听器),可适合装载在运动平台或海底固定式使用的多途信道中的声屏蔽方法。
本发明的目的是这样实现的:
(1)由声速仪测量水的声速,通过建模预报目标和干扰的信道冲激响应函数,或者通过接收探测声源的发射信号通过拷贝相关估计信道;
(2)水平接收阵接收声场中的多个点源的辐射信号,放大滤波,然后进行数据采集,采集包括需要探测的目标即被聚焦目标的辐射信号和干扰信号即被屏蔽目标的辐射信号;
(3)对接收信号作快速傅立叶变换,在频域采用相移补偿的球面波声场阵聚焦波束形成,或基于步骤(1)估计的信道得到被动时间反转镜(TRM)相位共轭谱,实现球面波声场中阵的声聚焦,得到的权值即为最佳聚焦权;
(4)在最佳聚焦权的基础上,增加屏蔽因子,得到最佳聚焦屏蔽权,寻找该最优权的准则为:在干扰处输出为0,在目标处输出最大;
(5)根据一系列专门的最佳屏蔽权频域矩阵算法,进行反傅立叶变换可得到输出波形,也可进行频域积分得到输出功率;
(6)后续相干处理或输出经声屏蔽、聚焦后的功率。
本发明以声屏蔽原理、方法为基础,适用于水声工程中抗多点源相干干扰,并且工程上可采用水平阵,更适合装载在运动平台上或海底固定使适用。
本发明研究局部空间(而不是扇面内的)的声屏蔽技术,抑制屏蔽区内声源产生的到达信号,改善对该区域外目标的探测性能,甚至可以屏蔽同一方向上不同距离处的任一声源,而有效检测另一声源。
本发明所谓的声屏蔽是:利用信号处理技术可在选定的局部区域内抑制其中的声源产生的相干干扰,从而改善对屏蔽区外目标的探测性能。
举例来说,拖船干扰是拖曳线列阵声纳的主要干扰,在拖船方向常造成约±20°的强干扰扇面,采用自适应波束形成技术可抑制该干扰,但在该扇面内仍是探测盲区。声屏蔽技术可在拖船区域设置屏蔽区,从而抑制拖船干扰,但在该扇面内屏蔽区外仍能有效探测目标。此例说明了本发明的声屏蔽技术与传统的阵处理技术的本质差别,声屏蔽技术可有效抗多目标干扰。
本发明有广泛的应用前景。例如,多种主动式和被动式水声设备在同频段同时工作时,会相互干扰,以致不能同时有效工作,可用于屏蔽拖船干扰,使拖曳线列阵声纳能有效消除拖船方向的探测盲区。
(四)附图说明
图1(a)和图1(b)为球面波声场中四元阵声屏蔽及声聚焦几何图。
图2为多途信道中水平阵声屏蔽及聚焦系统框图。
图3为抗拖船干扰声屏蔽几何态势图。
(五)具体实施方式
下面结合附图举例对本发明做更详细地描述:
本发明的目的是不仅用于设置零点抗混响,而且用于设置屏蔽区抗相干噪声源干扰。并且本发明可应用于水平阵(或单水听器),可适合装载在运动平台或海底固定式使用,从而应用更方便。
为实现本发明的目的,本发明采取的多途信道中的声屏蔽及声聚焦实现包括如下步骤:
(1)由声速仪测量水的声速,通过建模预报目标和干扰的信道冲激响应函数,或者通过接收探测声源的发射信号通过拷贝相关估计信道。
(2)水平接收阵接收声场中的多个点源的辐射信号,放大滤波,然后进行数据采集,包括需要探测的目标(被聚焦目标)的辐射信号,和干扰信号(被屏蔽目标)的辐射信号。
(3)对接收信号作快速傅立叶变换,在频域采用相移补偿的球面波声场阵聚焦波束形成,或基于步骤(1)估计的信道得到被动时间反转镜(TRM)相位共轭谱,两种方法是异曲同工的,均实现了球面波声场中阵的声聚焦,此时得到的权值即为最佳聚焦权。
(4)在最佳聚焦权的基础上,增加屏蔽因子,即得到最佳聚焦屏蔽权,寻找该最优权的准则为:在干扰处输出为0,在目标处输出最大,从而实现抗干扰的同时,对目标聚焦。
(5)根据一系列专门的最佳屏蔽权频域矩阵算法,进行反傅立叶变换可得到输出波形,也可进行频域积分得到输出功率。
(6)允许后续相干处理,从而提高水声系统的探测能力。也可输出经声屏蔽、聚焦后的功率。
为简明说明发明方法的物理机理,以球面波声场中四元水平阵的声聚焦及声屏蔽为例述如下。
图1(a)和图1(b)为球面波声场中四元阵声屏蔽及声聚焦几何图。声场中若有两个点源S1和S2,它们的辐射信号分别为zs(t)和zn(t)。希望探测zs(t),而屏蔽zn(t),即屏蔽点设置在S2处,而聚焦点在S1处。
各阵元输出信号的相位差是由球面波传播时延T决定的。rsi为S1到第i阵元的距离,rni为S2到第i阵元的距离,c为声速。Tsi和Tni为传播时延。
对各阵元输出进行相位差(或者时延差)补偿后加权相加即得到聚焦输出。可表述为:
Hs=[Hs1,Hs2,Hs3,Hs4]T    Hn=[Hn1,Hn2,Hn3,Hn4]T    (1)
可求得相应的聚焦峰的功率。从而得到聚焦输出的频谱,即时间反转镜(TRM)相位共轭谱表示法。换而言之,TRM实现了球面波声场中阵的声聚焦。
在多途信道中,(由Experimental demonstration of adaptive reverberationnulling using time reversal.《J.Acoust.soc.Am》,september 2005;118(3):1381-1387可知,)对S2屏蔽,且对S1聚焦的最佳权Wop为:
W op = [ I - H n * ( H n T H n * ) - 1 H n T ] H s * - - - ( 2 )
时空滤波器Wop的输出为:
y = y s + y n = Z s H s T [ I - H n * ( H n T H n * ) - 1 H n T ] H s * + Z s H n T [ I - H n * ( H n T H n * ) - 1 H n T ] H s * - - - ( 3 )
其中被屏蔽分量yn为:
y n = Z n [ H n T - H n T H n * ( H n T H n * ) - 1 H n T ] H s * = Z n [ H n T - H n T ] H s * = 0 - - - ( 4 )
上式说明S2辐射的Zn(t),它确实被屏蔽了。
聚焦输出ys,由(3)式可知:
y s = Z s [ H s T H s * - H s T H n * ( H n T H n * ) - 1 H n T H s * ] - - - ( 5 )
上式右边第一项即TRM共轭相位谱输出,它是聚焦输出;第二项为小量。只要与S1和S2相联系的信道频率响应Hs和Hn相差较大,即二者的互谱较小,甚至二者是正交的,则(5)式可达到良好的聚焦效果。
HsHn=0    (6)
当(6)式成立时,(5)式达到最佳效果。
上述说明了声屏蔽及声聚焦的原理。在图1(b)的情况下,两个点源均在阵的端向连线上,在球面波条件下无法区分S1和S2,但在多途信道中Hs和Hn的差别很大,二者的相关性很小,可以分别实现对S1聚焦,对S2屏蔽;或反之。
根据上述发明研究成果,本发明采用下述技术路线:
结合图2。考察多途信道中有两个点源S1和S2,分别辐射Zs(t)和Zn(t),二者互相独立。用一个水平基阵屏蔽Zn(t),而对Zs(t)聚焦,阵元数为N。
Hs1,Hs2,......,HsN为S1到各阵元间的信道频率响应。
Hn1,Hn2,......,HnN为S2到各阵元间的信道频率响应。
n1(t),n2(t),......,nN(t)为各阵元接收的各向同性干扰噪声,它们互相独立。
图2中的最佳权Wop为:
W op = [ I - H n * ( H n T H n * ) - 1 H n T ] H 0 * - - - ( 7 )
上式中的H0为搜索信道矢量,当搜索点与S1位置重合时有:
H0=HS    (8)
此时得到的声屏蔽及声聚焦输出在形式上同(4)式和(5)式。
这里以拖曳阵抗拖船干扰为例进行的说明本发明的研究应用,研究表明声屏蔽抗拖船干扰效果良好,可以消除拖船方向的探测盲区。
结合图3。跟前面所述一致,S1为需要探测的目标(聚焦目标),而S2为拖船产生的干扰(屏蔽目标),将其视为一点源,目标和干扰都在水平阵的端射方向。
把目标S1辐射的信号与干扰S2辐射的信号的功率比取分贝形式定义为“相干信噪比”,记为SNR;把目标S1、干扰S2辐射的信号与本地噪声的功率比取分贝形式定义为“非相干信噪比”,简称“信噪比”,分别记为snr1和snr2。
由声速仪测量水的声速。接收阵对接收水声信号放大滤波,然后进行数据采集。用拷贝相关器估计信道的冲击响应,它的Fourier变换即为Hs和Hn,在PC机中按(7)式做声屏蔽及聚焦分析。

Claims (1)

1、一种多途信道中的声屏蔽方法,其特征是:
(1)由声速仪测量水的声速,通过建模预报目标和干扰的信道冲激响应函数,或者通过接收探测声源的发射信号通过拷贝相关估计信道;
(2)水平接收阵接收声场中的多个点源的辐射信号,放大滤波,然后进行数据采集,采集包括需要探测的目标即被聚焦目标的辐射信号和干扰信号即被屏蔽目标的辐射信号;
(3)对接收信号作快速傅立叶变换,在频域采用相移补偿的球面波声场阵聚焦波束形成,或基于步骤(1)估计的信道得到被动时间反转镜(TRM)相位共轭谱,实现球面波声场中阵的声聚焦,得到的权值即为最佳聚焦权;
(4)在最佳聚焦权的基础上,增加屏蔽因子,得到最佳聚焦屏蔽权,寻找该最优权的准则为:在干扰处输出为0,在目标处输出最大;
(5)根据一系列专门的最佳屏蔽权频域矩阵算法,进行反傅立叶变换可得到输出波形,也可进行频域积分得到输出功率;
(6)后续相干处理或输出经声屏蔽、聚焦后的功率。
CN200810136823A 2008-07-30 2008-07-30 多途信道中的声屏蔽方法 Expired - Fee Related CN101335574B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810136823A CN101335574B (zh) 2008-07-30 2008-07-30 多途信道中的声屏蔽方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810136823A CN101335574B (zh) 2008-07-30 2008-07-30 多途信道中的声屏蔽方法

Publications (2)

Publication Number Publication Date
CN101335574A true CN101335574A (zh) 2008-12-31
CN101335574B CN101335574B (zh) 2012-09-26

Family

ID=40197921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810136823A Expired - Fee Related CN101335574B (zh) 2008-07-30 2008-07-30 多途信道中的声屏蔽方法

Country Status (1)

Country Link
CN (1) CN101335574B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101645264B (zh) * 2009-07-17 2012-03-14 中国科学院声学研究所 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法
CN103915102A (zh) * 2014-03-12 2014-07-09 哈尔滨工程大学 一种lfm水声多途信号的噪声抑制方法
CN103983946A (zh) * 2014-05-23 2014-08-13 北京神州普惠科技股份有限公司 一种声源定位过程中的多测量通道信号处理方法
CN104538019A (zh) * 2014-12-25 2015-04-22 哈尔滨工程大学 一种广义声子晶体隔声罩
CN109870695A (zh) * 2019-02-21 2019-06-11 哈尔滨工程大学 基于深海海底反射声的非合作目标多节点水声定位方法
CN110492950A (zh) * 2019-07-16 2019-11-22 西北工业大学 面向链路间干扰抑制的时间反转水声网络多址接入方法
CN115242583A (zh) * 2022-07-27 2022-10-25 中国科学院声学研究所 一种基于水平线列阵的信道脉冲响应被动估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101166065B (zh) * 2007-07-24 2010-10-06 哈尔滨工程大学 深海远程水声通信方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101645264B (zh) * 2009-07-17 2012-03-14 中国科学院声学研究所 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法
CN103915102A (zh) * 2014-03-12 2014-07-09 哈尔滨工程大学 一种lfm水声多途信号的噪声抑制方法
CN103915102B (zh) * 2014-03-12 2017-01-18 哈尔滨工程大学 一种lfm水声多途信号的噪声抑制方法
CN103983946A (zh) * 2014-05-23 2014-08-13 北京神州普惠科技股份有限公司 一种声源定位过程中的多测量通道信号处理方法
CN104538019A (zh) * 2014-12-25 2015-04-22 哈尔滨工程大学 一种广义声子晶体隔声罩
CN109870695A (zh) * 2019-02-21 2019-06-11 哈尔滨工程大学 基于深海海底反射声的非合作目标多节点水声定位方法
CN109870695B (zh) * 2019-02-21 2023-02-03 哈尔滨工程大学 基于深海海底反射声的非合作目标多节点水声定位方法
CN110492950A (zh) * 2019-07-16 2019-11-22 西北工业大学 面向链路间干扰抑制的时间反转水声网络多址接入方法
CN115242583A (zh) * 2022-07-27 2022-10-25 中国科学院声学研究所 一种基于水平线列阵的信道脉冲响应被动估计方法

Also Published As

Publication number Publication date
CN101335574B (zh) 2012-09-26

Similar Documents

Publication Publication Date Title
CN101335574B (zh) 多途信道中的声屏蔽方法
Blomberg et al. Improving sonar performance in shallow water using adaptive beamforming
Leigsnering et al. Multipath exploitation in through-the-wall radar imaging using sparse reconstruction
Kuperman et al. Shallow-water acoustics
Ahmad et al. Synthetic aperture beamformer for imaging through a dielectric wall
CN107561486B (zh) 一种基于主动时反的浅海目标波达方向估计方法
CN108845309A (zh) 一种基于主动虚拟时反方法的水下目标doa估计方法
CN106814360B (zh) 一种基于线性调频信号的多波束测深系统
Stergiopoulos Implementation of adaptive and synthetic-aperture processing schemes in integrated active-passive sonar systems
Liu et al. Low sidelobe range profile synthesis for sonar imaging using stepped-frequency pulses
Chavali et al. Multiplicative and min processing of experimental passive sonar data from thinned arrays
JP3367462B2 (ja) アクティブソーナー及びその目標検出方法
Godin et al. Application of time reversal to passive acoustic remote sensing of the ocean
Chavali et al. Coprime processing for the Elba Island sonar data set
Borcea et al. Coherent interferometric imaging, time gating and beamforming
Song et al. High-resolution modal wavenumber estimation in range-dependent shallow water waveguides using vertical line arrays
Bradaric et al. Waveform diversity for different multistatic radar configurations
CN102353960A (zh) 前视避碰声纳浅海稳健高分辨方位估计方法
Heaney Shallow water narrowband coherence measurements in the Florida Strait
Xionghou et al. Improved capon beamformer using Fourier integral method (FIM)
Foroozan et al. Time reversal MIMO radar for angle-Doppler estimation
Pan et al. Robust time-reversal is combined with distributed multiple-input multiple-output sonar for detection of small targets in shallow water environments
Zhou et al. Fluctuating modal interference characteristics in shallow water with a seasonal thermocline
Sabra et al. Effect of ocean currents on the performance of a time-reversing array in shallow water
Chen et al. Range sidelobe suppression for FD-MIMO sonar imaging with FDCW pulses

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20081231

RA01 Restoration of patent right

Former decision: deemed withdrawal of patent application after publication

Former decision publication date: 20111109

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120926

Termination date: 20180730

CF01 Termination of patent right due to non-payment of annual fee