CN112327305A - 一种快速频域宽带mvdr声纳波束形成方法 - Google Patents
一种快速频域宽带mvdr声纳波束形成方法 Download PDFInfo
- Publication number
- CN112327305A CN112327305A CN202011227777.3A CN202011227777A CN112327305A CN 112327305 A CN112327305 A CN 112327305A CN 202011227777 A CN202011227777 A CN 202011227777A CN 112327305 A CN112327305 A CN 112327305A
- Authority
- CN
- China
- Prior art keywords
- mvdr
- frequency domain
- sonar
- data
- fast
- 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
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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
-
- 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/537—Counter-measures or counter-counter-measures, e.g. jamming, anti-jamming
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种快速频域宽带MVDR声呐波束形成方法,通过采集设备获取声纳阵元域数据,做FFT将其转换到频域,取对应滤波截止频率范围内的数据,所得数据对应的频点乘以对应相位差,构建协方差矩阵,计算输出功率。经与实际装备比对测试,本发明方法的波束宽度为CBF的1/2左右,对相干噪声的抑制能力比CBF算法提升约6‑9dB,信号检测和方位分辨性能较优,和普通宽带MVDR相比,可在单次快拍下获得稳定优化解,运算量较小,有助于改善声纳的探测性能。
Description
技术领域
本发明属于水声工程领域,涉及一种快速频域宽带MVDR声纳波束形成方法。
背景技术
水声阵列信号处理主要是利用声信号的空域特性来增强信号及有效提取信号空域信息。与传统的单个定向传感器相比,水声阵列信号处理具有灵活的波束控制、高信号增益、极强的抗干扰能力及高空间分辨能力等优点。
波束形成技术是声纳信号处理的核心技术,是获取目标方位信息和进行信号检测的主要手段。波束形成技术的核心思想是将一定几何形状(如直线阵或圆阵)排列的多元基阵各阵元输出经过处理,如对多个阵元接收信号进行加权、时延或相位补偿、求和处理等,使其对预定来波方向的入射信号形成同相相加。因此,一个波束形成器可以看成是一个空间滤波器,它可以滤去空间某些方位的信号,只让指定方位的信号通过。
目前各类声纳主要采用常规波束形成技术(conventional beamforming,简称CBF),即通过对传感器阵列在时域上的延时累加或者在频域上的移相相加获得空间增益从而获得目标方位信息。然而常规波束形成技术形成的波束分辨率、检测能力受到阵列尺寸和工作频率所决定的瑞利极限限制而无法提升,主要体现在以下三个方面的问题:
(1)波束宽
根据CBF原理,在基阵孔径和工作频率一定的情况下,波束宽度是无法进一步改进的。在目标数量较少时该问题不突出,然而在多目标时问题显得突出。波束宽度较宽将极大程度影响多目标下的噪声听测,容易导致混合目标听测,较大程度影响听测效果,难以做出准确识别。
(2)旁瓣高
CBF的旁瓣高度约为主瓣的10%,根据CBF原理,旁瓣是无法避免的。旁瓣的存在容易导致假目标的产生,难以区分是旁瓣还是弱目标,因此对于强目标干扰下的弱目标检测尤其不利。
(3)相干噪声抑制能力弱
CBF的基本原理是在假定各路水听器获取的目标信号是相关的,而他们的背景噪声是不相关的,该假定与大多数的实际情况是基本吻合的,因此CBF总体效果不错。然而,海上实际情况表明,在某些情况下,各路水听器的背景噪声具有一定的相关性,此时CBF不再是最佳的检测手段。
提高波束形成干扰抑制能力的一种有效措施是采用自适应波束形成方法,其中最为典型的是Capton提出的最小方差无畸变响应(Minimum Variance DistortionlessResponse Beamforming,MVDR)波束形成技术(Capon J.High-resolution frequency-wavenumber spectrum analysis.Proceedings of the IEEE 1969;57(8):1408-18.)。它在保持波束指向方向信号无失真的条件下,通过使基阵输出功率最小来实现对干扰的有效抑制,具有较好的方位分辨能力和较强的干扰抑制能力,尤其在低信噪比条件下,具有良好的强相关噪声抑制能力。
宽带信号的MVDR频域实现方法(简称频域MVDR)为首先将数据进行分块(每块数据记为一次快拍),然后将每次快拍内的数据通过傅里叶变换分解为多个窄带,在每个窄带上利用多次快拍估计阵列协方差矩阵,利用该协方差矩阵估计优化权向量得到窄带上的优化输出,最后综合各个窄带的输出结果。对于固定短阵列而言,频域MVDR在实际应用中效果并不理想,主要表现在①由于频域MVDR方法为获得满秩和稳定的协方差矩阵估计,需要对接收数据进行多次分块处理,因而要求信号平稳时间较长。对于方位快变目标而言,在该时间段内,方位变化可能达到几个波束宽度,而自适应波束形成方法的权向量不能足够快地适应非平稳信号时,其波束形成性能将急剧降低(Vorobyov SA,Gershman AB,Luo ZQ,MaN.Adaptive beamforming with joint robustness against mismatched signalsteering vector and interference nonstationarity.IEEE Signal ProcessingLetters 2004;11(2Part 1):108-11.)。而当快拍数较少时,协方差矩阵估计不稳定,难以获得稳定的优化权向量(Reed IS,Mallet J,Brennan L.Rapid convergence rate inadaptive antennas.IEEE Trans.Aerosp.Electron.Syst 1974;10.)。②在存在干扰源扰动时,MVDR不能很好的将零点对准干扰方向,将导致弱信号检测能力下降(Song H,Kuperman WA,Hodgkiss WS,Gerstoft P,Kim JS.Null broadening with snapshot-deficient covariance matrices in passive sonar.IEEE Journal of OceanicEngineering 2003;28(2):250-61.)。
总结近年来有关MVDR方法的研究成果,基本都是围绕在减小频域MVDR方法对快拍数的需求上,即如何在短数据条件下获得稳定的优化权向量估计。虽然以上方法都在一定程度上作出了改善,但由于频域MVDR方法本身难以获得足够多的快拍次数,因此在改善的同时也付出了代价。然而,由于MVDR方法具有良好的干扰和噪声抑制能力,在低信噪比条件下也能获得良好的方位估计性能,这对于被动声纳探测而言,无疑具有重要的军事意义,因此,如何有效改善目前的频域MVDR方法在短数据条件下的方位估计性能,具有重要的现实意义。
发明内容
针对上述问题,本发明的目的是提供一种形成波束宽度窄、对相干噪声的抑制能力强、信号检测和方位分辨性能高的快速频域宽带MVDR声纳波束形成方法,以提升声纳对水下目标的探测能力。
本发明的技术方案是:
步骤一:通过采集设备获取声纳M路N点阵元域数据,其中,N≥32M且为2的整次幂;
步骤二:将M路信号xm(n),m=1,2,...M,n=1,2,...N做N点的快速傅里叶变换(Fast Fourier Transform,简称FFT,)将xm(n)转换到频域,取对应滤波截止频率范围[NL,NH]内的数据:X(fn)=fft[x(n)];
设与[fL,fH]对应的fn分别为[NL,NH],则有:NL=fL/fs×N,NH=fH/fs×N;
只取滤波频段内数据并进行处理得到:X'(fn)=2·X(fn)NL<n<NH;
步骤三:将步骤二所得X'(fn)对应的频点乘以对应相位差exp(j·2π·fn·τm(θi)),即Y(fn)=X'(fn)exp[j2πfnτm(θi)];
式中,τm(θi)是第i个信号源到达第m号阵元相对于参考阵元的时间延迟,它只和目标信号所处的方位有关;
本发明的优点:
经与实际装备比对测试,本发明方法的波束宽度为CBF的1/2左右,对相干噪声的抑制能力比CBF算法提升约6-9dB,信号检测和方位分辨性能较优;与普通宽带MVDR相比,可在单次快拍下获得稳定优化解,运算量较小,有助于改善声纳的探测性能。
附图说明
图1是基于快速频域MVDR算法的声纳波束形成方法的流程图。
图2是不同信噪比下本发明方法处理结果与CBF处理结果对比图。
图3是利用本发明方法处理7000s海上实际声纳数据的结果与CBF结果对比图。
具体实施方式
为了便于本领域技术人员理解,下面结合附图和具体实施例对本发明的技术方案进行清楚、完整地描述。显然,以下所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
实施例1:
步骤一:仿真生成16路4096点阵元域数据。均匀间隔直线线列阵,阵元数目N为16,阵元间距为0.07米,存在2个目标;目标1信号为宽频色噪声,入射角为70度,目标强度10dB保持不变;目标2信号为宽频色噪声,入射角为40-100度变化,目标强度10dB保持不变;采样频率为20kHz,处理带宽为3kHz-8kHz,背景噪声为非相干高斯白噪声。
步骤二:将16路信号xm(n),m=1,2,...16,n=1,2,...N做4096点的FFT,将xm(n)转换到频域X(fn)=fft[x(n)],取对应滤波截止频率范围[614,1638]内的数据进行处理得到:X'(fn)=2·X(fn)614<n<1638;
步骤三:将步骤二所得X'(fn)对应的频点乘以对应相位差exp(j·2π·fn·τm(θi)),即Y(fn)=X'(fn)exp[j2πfnτm(θi)];
分别对比分析不同信噪比时CBF方法与本发明方法的处理结果,如图2所示。从图中可以看出:
(1)在信噪比较低时,CBF方法与本发明方法的性能基本相当;在信噪比较高时,快本发明方法的性能明显变好,波束明显变窄,背景抑制能力明显优于CBF算法。
(2)几乎所有条件下,本发明方法的结果输出曲线低于CBF算法,体现了MVDR算法在最小功率输出约束条件下的自适应加权效果,其性能不低于CBF算法。
实施例2:
步骤一:通过采集设备获取实际圆阵声纳32路4096点阵元域数据,采样率40kHz。
步骤二:将32路信号xm(n),m=1,2,...32,n=1,2,...N做4096点的FFT,将xm(n)转换到频域X(fn)=fft[x(n)],取对应滤波截止频率范围[1228,3276]内的数据进行处理得到:X'(fn)=2·X(fn)1228<n<3276
步骤三:将步骤二所得X'(fn)对应的频点乘以对应相位差exp(j·2π·fn·τm(θi)),即Y(fn)=X'(fn)exp[j2πfnτm(θi)]
分别利用本发明方法和CBF处理7000s海上实际声纳数据,其结果如图3所示。从图中可以看出:本发明方法形成的波束宽度是CBF的1/2左右,方位分辨性能较优,对于多目标分辨以及多目标下的弱目标探测具有帮助,并且表现出较好的检测能力。
对于本领域的普通技术人员而言,根据本发明的教导,在不脱离本发明的原理与精神的情况下,对实施方式所进行的改变、修改、替换和变型仍落入本发明的保护范围之内。
Claims (4)
1.一种快速频域宽带MVDR声纳波束形成方法,其特征在于包括以下步骤:
步骤一:通过采集设备获取声纳M路N点阵元域数据;
步骤二:将M路信号xm(n),m=1,2,...M,n=1,2,...N做N点的快速傅里叶变换,将xm(n)转换到频域,取对应滤波截止频率范围[NL,NH]内的数据:X(fn)=fft[x(n)];
设与[fL,fH]对应的fn分别为[NL,NH],则有:NL=fL/fs×N,NH=fH/fs×N;
步骤三:将步骤二所得X'(fn)对应的频点乘以对应相位差exp(j·2π·fn·τm(θi)),即
Y(fn)=X'(fn)exp[j2πfnτm(θi)];
式中,τm(θi)是第i个信号源到达第m号阵元相对于参考阵元的时间延迟,它只和目标信号所处的方位有关;
2.如权利要求1所述的快速频域宽带MVDR声纳波束形成方法,其特征在于:步骤一中N≥32M且为2的整次幂。
3.如权利要求1所述的快速频域宽带MVDR声纳波束形成方法,其特征在于:步骤二中频域带通滤波后只取滤波频段内数据并进行处理得到:X'(fn)=2·X(fn)NL<n<NH。
4.如权利要求1所述的快速频域宽带MVDR声纳波束形成方法,其特征在于:步骤五中的I=[1,1,…,1],上标T表示转置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011227777.3A CN112327305B (zh) | 2020-11-06 | 2020-11-06 | 一种快速频域宽带mvdr声纳波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011227777.3A CN112327305B (zh) | 2020-11-06 | 2020-11-06 | 一种快速频域宽带mvdr声纳波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112327305A true CN112327305A (zh) | 2021-02-05 |
CN112327305B CN112327305B (zh) | 2022-10-04 |
Family
ID=74315646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011227777.3A Active CN112327305B (zh) | 2020-11-06 | 2020-11-06 | 一种快速频域宽带mvdr声纳波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112327305B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114740424A (zh) * | 2022-04-08 | 2022-07-12 | 中国人民解放军海军潜艇学院 | 一种基于自适应波束合成的频域非失真听音方法 |
CN114879170A (zh) * | 2022-04-26 | 2022-08-09 | 中国人民解放军海军潜艇学院 | 一种用于mvdr波束形成技术的数据坏道消除方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060133211A1 (en) * | 2004-12-17 | 2006-06-22 | Yang Tsih C | Method and apparatus for acoustic source tracking using a horizontal line array |
EP1919251A1 (en) * | 2006-10-30 | 2008-05-07 | Mitel Networks Corporation | Beamforming weights conditioning for efficient implementations of broadband beamformers |
US20080240463A1 (en) * | 2007-03-29 | 2008-10-02 | Microsoft Corporation | Enhanced Beamforming for Arrays of Directional Microphones |
WO2013180269A1 (ja) * | 2012-05-31 | 2013-12-05 | 日立アロカメディカル株式会社 | 超音波撮像装置 |
CN104101871A (zh) * | 2013-04-15 | 2014-10-15 | 中国科学院声学研究所 | 一种用于被动合成孔径的抑制窄带干扰方法及系统 |
CN104730513A (zh) * | 2013-12-19 | 2015-06-24 | 中国科学院声学研究所 | 一种分级子阵聚焦mvdr波束形成方法 |
CN204559608U (zh) * | 2015-02-15 | 2015-08-12 | 重庆电子工程职业学院 | 一种自启动风扇散热家用路由器 |
CN105445723A (zh) * | 2015-12-22 | 2016-03-30 | 中国船舶重工集团公司第七一五研究所 | 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法 |
CA2941477A1 (en) * | 2015-09-17 | 2017-03-17 | Navico Holding As | Adaptive beamformer for sonar imaging |
CN107783081A (zh) * | 2017-10-19 | 2018-03-09 | 杭州睿达汽车科技有限公司 | 一种针对宽带雷达的sstmv最小方差波束形成方法 |
US20190122686A1 (en) * | 2017-10-19 | 2019-04-25 | Kardome Technology Ltd. | Speech enhancement using clustering of cues |
CN110868356A (zh) * | 2019-11-26 | 2020-03-06 | 东北林业大学 | 一种智能安全路由的控制方法及其智能安全路由器 |
CN210769463U (zh) * | 2019-08-30 | 2020-06-16 | 广州金升阳科技有限公司 | 一种风扇的供电控制电路 |
-
2020
- 2020-11-06 CN CN202011227777.3A patent/CN112327305B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060133211A1 (en) * | 2004-12-17 | 2006-06-22 | Yang Tsih C | Method and apparatus for acoustic source tracking using a horizontal line array |
EP1919251A1 (en) * | 2006-10-30 | 2008-05-07 | Mitel Networks Corporation | Beamforming weights conditioning for efficient implementations of broadband beamformers |
US20080240463A1 (en) * | 2007-03-29 | 2008-10-02 | Microsoft Corporation | Enhanced Beamforming for Arrays of Directional Microphones |
WO2013180269A1 (ja) * | 2012-05-31 | 2013-12-05 | 日立アロカメディカル株式会社 | 超音波撮像装置 |
CN104101871A (zh) * | 2013-04-15 | 2014-10-15 | 中国科学院声学研究所 | 一种用于被动合成孔径的抑制窄带干扰方法及系统 |
CN104730513A (zh) * | 2013-12-19 | 2015-06-24 | 中国科学院声学研究所 | 一种分级子阵聚焦mvdr波束形成方法 |
CN204559608U (zh) * | 2015-02-15 | 2015-08-12 | 重庆电子工程职业学院 | 一种自启动风扇散热家用路由器 |
CA2941477A1 (en) * | 2015-09-17 | 2017-03-17 | Navico Holding As | Adaptive beamformer for sonar imaging |
CN105445723A (zh) * | 2015-12-22 | 2016-03-30 | 中国船舶重工集团公司第七一五研究所 | 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法 |
CN107783081A (zh) * | 2017-10-19 | 2018-03-09 | 杭州睿达汽车科技有限公司 | 一种针对宽带雷达的sstmv最小方差波束形成方法 |
US20190122686A1 (en) * | 2017-10-19 | 2019-04-25 | Kardome Technology Ltd. | Speech enhancement using clustering of cues |
CN210769463U (zh) * | 2019-08-30 | 2020-06-16 | 广州金升阳科技有限公司 | 一种风扇的供电控制电路 |
CN110868356A (zh) * | 2019-11-26 | 2020-03-06 | 东北林业大学 | 一种智能安全路由的控制方法及其智能安全路由器 |
Non-Patent Citations (2)
Title |
---|
VAN VEEN B D ET.AL: "《Minimum variance beam forming with soft response constraints》", 《IEEE TRANS SIGNAL PROCESSING》 * |
黄清: "《自适应宽带稳健波束形成及神经网络算法》", 《东南大学学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114740424A (zh) * | 2022-04-08 | 2022-07-12 | 中国人民解放军海军潜艇学院 | 一种基于自适应波束合成的频域非失真听音方法 |
CN114740424B (zh) * | 2022-04-08 | 2024-10-25 | 中国人民解放军海军潜艇学院 | 一种基于自适应波束合成的频域非失真听音方法 |
CN114879170A (zh) * | 2022-04-26 | 2022-08-09 | 中国人民解放军海军潜艇学院 | 一种用于mvdr波束形成技术的数据坏道消除方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112327305B (zh) | 2022-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106788653A (zh) | 一种基于协方差矩阵重构的自适应波束形成方法 | |
Su et al. | Adaptive beamforming for nonstationary HF interference cancellation in skywave over-the-horizon radar | |
KR20060044285A (ko) | 위상 배열 방식 시스템용 잡음 억제 시스템 및 방법 | |
KR20040102039A (ko) | 레이더의 스펙트럼 생성 시스템 및 방법 | |
CN106842140B (zh) | 一种基于和差波束降维的主瓣干扰抑制方法 | |
CN109407055A (zh) | 基于多径利用的波束形成方法 | |
CN113835068B (zh) | 一种基于独立成分分析的盲源分离实时抗主瓣干扰方法 | |
CN112327305B (zh) | 一种快速频域宽带mvdr声纳波束形成方法 | |
CN110261826A (zh) | 一种零陷展宽的相干干扰抑制方法 | |
CN111693971B (zh) | 一种用于弱目标检测的宽波束干扰抑制方法 | |
CN110687528A (zh) | 自适应波束形成器生成方法及系统 | |
CN106842135B (zh) | 基于干扰加噪声协方差矩阵重构的自适应波束形成方法 | |
CN115085784B (zh) | 一种分级加权稀疏子阵天线阵列稳健波束合成方法 | |
CN110208757B (zh) | 一种抑制主瓣干扰的稳健自适应波束形成方法及装置 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN114487985B (zh) | 一种基于差-和信号的波束锐化方法及系统 | |
CN107728132B (zh) | 一种改善干扰阻塞算法输出信噪比的方法 | |
CN110261814A (zh) | 基于空间谱重构和导向矢量直接估计的波束形成方法 | |
CN113156392B (zh) | 一种基于俯仰域自适应处理的杂波抑制方法 | |
CN116226611A (zh) | 基于分数域反卷积波束形成的啁啾信号波达方向估计方法 | |
CN114371441A (zh) | 虚拟阵列波达方向估计方法、装置、产品及存储介质 | |
CN112444790A (zh) | 一种在强干扰情况下目标检测的方法 | |
CN114152918A (zh) | 基于压缩感知的抗间歇式主瓣干扰方法 | |
Xu et al. | Coherent source direction estimation for three-row bathymetric sidescan sonars | |
Wang et al. | A long line array broadband adaptive beamforming with post-filtering and incoherent noise reduction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |