CN103792528A - 一种基于对角减载的水声阵列Bartlett波束形成的方法 - Google Patents

一种基于对角减载的水声阵列Bartlett波束形成的方法 Download PDF

Info

Publication number
CN103792528A
CN103792528A CN201410047825.9A CN201410047825A CN103792528A CN 103792528 A CN103792528 A CN 103792528A CN 201410047825 A CN201410047825 A CN 201410047825A CN 103792528 A CN103792528 A CN 103792528A
Authority
CN
China
Prior art keywords
array
signal
underwater sound
output power
load
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
CN201410047825.9A
Other languages
English (en)
Other versions
CN103792528B (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 CN201410047825.9A priority Critical patent/CN103792528B/zh
Publication of CN103792528A publication Critical patent/CN103792528A/zh
Application granted granted Critical
Publication of CN103792528B publication Critical patent/CN103792528B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046—Techniques for image enhancement involving transmitter or receiver
    • G01S7/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Signal Processing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明属于水声信号处理领域，具体涉及一种基于对角减载的水声阵列Bartlett波束形成的方法。本发明通过数据采集得到M次独立快拍的阵列信号，得到自相关矩阵；分别采集并计算各个通道有信号和无信号时的输出功率，对所有通道计算结果取平均，得到有信号时的平均输出功率和无信号时的平均输出功率；计算得到最佳的减载系数λ；得到对角减载后的自相关矩阵；进行常规波束形成，得到对角减载后的波束形成结果。本发明提供了一种基于对角减载的水声阵列Bartlett波束形成的方法，提高了阵列的输出信噪比，计算简单，适用于复杂的水声信道环境中，工程上可广泛应用于舷侧声呐等声呐基阵阵元位置固定已知的水声阵列。

Description

一种基于对角减载的水声阵列Bartlett波束形成的方法
技术领域
本发明属于水声信号处理领域，具体涉及一种基于对角减载的水声阵列Bartlett波束形成的方法。
背景技术
波束形成是指将一定几何形状排列的多元基阵各阵元输出经过处理形成空间指向性的方法，波束形成使得声呐基阵能量集中，获得阵列增益，得到抗噪声及混响等的作用。由于波束形成技术在声呐信号处理过程中的重要作用，多年来国内外学者一直致力于对其技术进行改进和完善，其中提高阵列增益和波束的主旁瓣比是要解决的核心问题。
在水声信道中，由于信道的复杂性和低信噪比的影响，多重信号分类算法（MUSIC）（参见：Adaptivity to Background Noise Spatial Coherence for High Resolution Passive Methods.《Proc.ICASSP》,1980,1:307-310）和旋转子空间不变算法（ESPRIT）（参见：ESPRIT-Estimation of Signal Parameters via Rotational Invariance Techniques.《IEEE Trans.ASSP》,1989,37(7):984-995）等要求有较高信噪比的信号子空间类波束形成技术便不能发挥其作用，Bartlett波束形成技术目前仍然是水声阵列信号处理的常用手段，并且在高斯白噪声背景下，Bartlett是最佳波束形成器（参见：Bartlett波束形成的波束零限权设计,《哈尔滨工程大学学报》，2008，29(12)：36-38；On the Performance of Energy Detection Using Bartlett’sEstimate for Spectrum Sensing in Cognitive Radio Systems.《IEEE TRANSACTIONS ONSIGNAL PROCESSING》,2012,60(7):3394-3404）,对角加载波束形成的方法（参见：Robustadaptive beamforming using sequential quadratic programming:An iterative solution to themismatch problem.《IEEE Signal Processing Letters》，2008，15：733-736；On robust Caponbeamforming and diagonal loading.《IEEE Trans.on Signal Processing》,2003,51(7):1702-1715）虽然能够很好的抑制噪声，然而当快拍数少和阵元位置存在误差时，该方法的性能显著下降，同时还存在加载量难以确定的问题。因而，研究一种既适用于水声信道的复杂性和低信噪比的情况，又能改善以往研究方法不足之处的方法是有必要的。
发明内容
本发明的目的在于针对水声阵列信号处理中受复杂的水声信道和低信噪比影响的问题，提出了一种提高了阵列的输出信噪比，减载量的选取简单灵活的基于对角减载的水声阵列Bartlett波束形成的方法。
本发明的目的是这样实现的：
（1）通过数据采集得到M次独立快拍的阵列信号X，为N×M维向量，由Rx＝E[XXH]得到自相关矩阵Rx，为N×N维向量，N为水听器阵列的阵元个数；
（2）分别采集并计算各个通道有信号和无信号时的输出功率，对所有通道计算结果取平均，得到有信号时的平均输出功率
Figure BDA0000465000590000021
和无信号时的平均输出功率
Figure BDA0000465000590000022
（3）基于步骤（2）计算得到的有信号和无信号时的平均输出功率，根据对阵列输出增益或波束主旁瓣比的需求，计算得到最佳的减载系数λ；
（4）在步骤（1）中得到的自相关矩阵Rx的主对角线上乘以步骤（3）中求得的最佳减载系数，得到对角减载后的自相关矩阵
Figure BDA0000465000590000023
（5）进行常规波束形成，
Figure BDA0000465000590000024
得到对角减载后的波束形成结果。
有信号时的平均输出功率为
Figure BDA0000465000590000025
无信号时的平均输出功率为
Figure BDA0000465000590000026
本发明的有益效果在于：本发明提供了一种基于对角减载的水声阵列Bartlett波束形成的方法，对角减载技术和Bartlett波束形成技术的结合，能够起到很好的抗噪声的作用，很大程度上提高了阵列的输出信噪比，并且本发明中对角减载采用的减载量可根据实际工作中对于阵列输出增益或者波束主旁瓣比的需求灵活选取，计算简单，适用于复杂的水声信道环境中，工程上可广泛应用于舷侧声呐等声呐基阵阵元位置固定已知的水声阵列。
附图说明
图1为基于对角减载的Bartlett波束形成流程图。
具体实施方式
为了更加详细的对本发明进行阐述，下面结合附图对本发明实施过程进行举例说明。
声纳基阵的阵元位置是固定并且已知的，基阵接收信号。参照图1，实施过程如下：
（1）首先由数据采集器采集得到声纳基阵中各个阵元M次独立快拍的阵列信号X，为N×M维向量，对阵列信号作自相关Rx＝E[XXH]，得到阵列信号的自相关矩阵Rx，为N×N维向量，N为水听器阵列的阵元个数；
（2）分别采集并计算各个通道有信号和无信号时的输出功率，记为和
Figure BDA0000465000590000028
对所有通道计算结果取平均，得到有信号时的平均输出功率为
Figure BDA0000465000590000029
无信号时的平均输出功率为
Figure BDA00004650005900000210
（3）下面进行对角线项上的减载系数的选择，以实际工作中需满足一定的阵列输出增益为例，波束形成后目标方位上的阵列输出增益为
G = 10 lg ( | ω H ρ s ω | | ω H ρ n ω | )
其中ρs、ρn分别为阵列输入平面波信号和噪声的归一化自相关矩阵，ω为Bartlett波束形成权矢量，
Figure BDA0000465000590000032
υk为Bartlett导向矢量。
若进行减载系数为λ的对角减载，则减载后得到 ρ s λ = ( υ s υ s H ) + ( λ - 1 ) × I , ρ n λ = λ × I , υs为平面波入射信号的方向矢量。当导向矢量对准目标时，即υk＝υs，得到基于对角减载的Bartlett波束形成的阵列增益为
G λ = 10 lg ( N - 1 + λ λ )
因而，得到了阵列输出增益和减载系数的关系式，选取一定的阵列输出增益后，由该式计算出最佳的减载系数。
（4）计算得到最佳减载系数后，在步骤（1）中得到的自相关矩阵Rx上乘以该减载系数，得到对角减载后的自相关矩阵
Figure BDA0000465000590000035
（5）进行Bartlett波束形成，
Figure BDA0000465000590000036
得到基于对角减载的Bartlett波束形成后输出平均功率。
理想状态下，自相关矩阵可写为
Figure BDA0000465000590000037
则不同通道间噪声信号互相关为零，仅剩下对角线上的噪声自相关分量，
Figure BDA0000465000590000038
为信号功率，
Figure BDA0000465000590000039
Figure BDA00004650005900000310
为噪声功率，因为噪声分量主要集中在对角线上，因而按照上述方式进行对角减载，很大程度上减小了噪声影响，提高阵列的输出信噪比，并且减载系数可根据设定的阵列输出增益计算得到。
本发明的优点是：对阵列信号的自相关矩阵进行适当的对角减载，减载系数由实际工作中对阵列输出增益或对主旁瓣比的要求进行优化选取，提高了Bartlett波束形成器的阵列增益，改善了主旁瓣比，克服了以往波束形成技术的不足，提高了声纳基阵在水声信道中的检测性能，方法简单易行，可广泛应用于阵元位置固定已知的所有声纳基阵中，本发明具有一定的应用前景。

Claims (2)

1.一种基于对角减载的水声阵列Bartlett波束形成的方法，其特征在于：
（1）通过数据采集得到M次独立快拍的阵列信号X，为N×M维向量，由Rx＝E[XXH]得到自相关矩阵Rx，为N×N维向量，N为水听器阵列的阵元个数；
（2）分别采集并计算各个通道有信号和无信号时的输出功率，对所有通道计算结果取平均，得到有信号时的平均输出功率
Figure FDA0000465000580000011
和无信号时的平均输出功率
Figure FDA0000465000580000012
（3）基于步骤（2）计算得到的有信号和无信号时的平均输出功率，根据对阵列输出增益或波束主旁瓣比的需求，计算得到最佳的减载系数λ；
（4）在步骤（1）中得到的自相关矩阵Rx的主对角线上乘以步骤（3）中求得的最佳减载系数，得到对角减载后的自相关矩阵
Figure FDA0000465000580000013
（5）进行常规波束形成，得到对角减载后的波束形成结果。
2.根据权利要求1所述的一种基于对角减载的水声阵列Bartlett波束形成的方法，其特征在于：所述的有信号时的平均输出功率为
Figure FDA0000465000580000015
无信号时的平均输出功率为 σ n 2 = 1 N Σ i = 1 N σ ni 2 .
CN201410047825.9A 2014-02-11 2014-02-11 一种基于对角减载的水声阵列Bartlett波束形成的方法 Active CN103792528B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410047825.9A CN103792528B (zh) 2014-02-11 2014-02-11 一种基于对角减载的水声阵列Bartlett波束形成的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410047825.9A CN103792528B (zh) 2014-02-11 2014-02-11 一种基于对角减载的水声阵列Bartlett波束形成的方法

Publications (2)

Publication Number Publication Date
CN103792528A true CN103792528A (zh) 2014-05-14
CN103792528B CN103792528B (zh) 2016-05-04

Family

ID=50668388

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410047825.9A Active CN103792528B (zh) 2014-02-11 2014-02-11 一种基于对角减载的水声阵列Bartlett波束形成的方法

Country Status (1)

Country Link
CN (1) CN103792528B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105445723A (zh) * 2015-12-22 2016-03-30 中国船舶重工集团公司第七一五研究所 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1361679A2 (en) * 2002-05-07 2003-11-12 Matsushita Electric Industrial Co., Ltd. Radio communication device and arrival direction estimation method
JP2005062058A (ja) * 2003-08-18 2005-03-10 Toshiba Corp 捜索レーダ装置
CN101609150A (zh) * 2009-07-07 2009-12-23 哈尔滨工程大学 一种提高基阵分辨力和增益的快速波束形成方法
US20110007606A1 (en) * 2008-02-18 2011-01-13 Thomas Edgar Curtis Underwater Surveillance
CN103308889A (zh) * 2013-05-13 2013-09-18 辽宁工业大学 复杂环境下被动声源二维doa估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1361679A2 (en) * 2002-05-07 2003-11-12 Matsushita Electric Industrial Co., Ltd. Radio communication device and arrival direction estimation method
JP2005062058A (ja) * 2003-08-18 2005-03-10 Toshiba Corp 捜索レーダ装置
US20110007606A1 (en) * 2008-02-18 2011-01-13 Thomas Edgar Curtis Underwater Surveillance
CN101609150A (zh) * 2009-07-07 2009-12-23 哈尔滨工程大学 一种提高基阵分辨力和增益的快速波束形成方法
CN103308889A (zh) * 2013-05-13 2013-09-18 辽宁工业大学 复杂环境下被动声源二维doa估计方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105445723A (zh) * 2015-12-22 2016-03-30 中国船舶重工集团公司第七一五研究所 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法

Also Published As

Publication number Publication date
CN103792528B (zh) 2016-05-04

Similar Documents

Publication Publication Date Title
CN103901395B (zh) 一种冲击噪声环境下相干信号波达方向动态跟踪方法
CN101369014B (zh) 应用于多输入多输出雷达的双边约束自适应波束形成方法
CN111123192B (zh) 一种基于圆形阵列和虚拟扩展的二维doa定位方法
CN106788653A (zh) 一种基于协方差矩阵重构的自适应波束形成方法
CN102393525B (zh) 子空间投影的导航干扰抑制与信号增强方法
CN103116162B (zh) 基于目标空间稀疏性的高分辨声呐定位方法
CN103176168B (zh) 一种机载非正侧视阵雷达近程杂波对消方法
CN103984676A (zh) 一种基于协方差矩阵重构的正交投影自适应波束形成方法
CN105301580A (zh) 一种基于分裂阵互谱相位差方差加权的被动探测方法
CN104360310A (zh) 一种多目标近场源定位方法和装置
CN103308904A (zh) 一种兼顾近场和远场性能的三维摄像声纳系统换能器阵列的稀疏优化方法
CN105022050A (zh) 一种多传感器阵列的水声信道离散噪声源抑制方法
CN105005038A (zh) 一种改进的声矢量阵相干源doa估计算法
CN108931766A (zh) 一种基于稀疏重构的非均匀stap干扰目标滤除方法
CN104502904A (zh) 一种鱼雷自导波束锐化方法
CN103902830A (zh) 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法
CN107703478A (zh) 基于互相关矩阵的扩展孔径二维doa估计方法
CN102064892A (zh) 子阵级自适应单脉冲的两级干扰抑制方法
CN103558584A (zh) 一种波达方向的检测方法及装置
CN105353340A (zh) 一种双层圆柱阵水下被动目标检测方法
CN102087357A (zh) 传感器阵列回波方向估计及多波束回波测深底检测方法
CN104156553A (zh) 无需信源数估计的相干信号波达方向估计方法及系统
CN103792528A (zh) 一种基于对角减载的水声阵列Bartlett波束形成的方法
CN114563760B (zh) 一种基于sca阵型的二阶超波束形成方法、设备及介质
CN102183755A (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