CN112711014A - 一种非均匀布阵舷侧阵声纳波束形成快速方法 - Google Patents
一种非均匀布阵舷侧阵声纳波束形成快速方法 Download PDFInfo
- Publication number
- CN112711014A CN112711014A CN202011464992.5A CN202011464992A CN112711014A CN 112711014 A CN112711014 A CN 112711014A CN 202011464992 A CN202011464992 A CN 202011464992A CN 112711014 A CN112711014 A CN 112711014A
- Authority
- CN
- China
- Prior art keywords
- array
- degrees
- elements
- frequency
- beams
- 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
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Discrete Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种非均匀布阵舷侧阵声纳波束形成快速方法,主要过程为首先将频域波束形成每个频点的补相位求和的计算过程转化为两个DFT卷积,然后进行筛选做一个近似计算。本发明进行工程可接受的近似计算,计算量小;对阵元位置没有严格要求,克服了现有快速算法通常要均匀布阵的限制;形成的波束个数与阵元数无关,波束个数及波束指向可任意设置。
Description
技术领域
本发明涉及声纳阵列信号处理中非均匀布阵舷侧阵的频域波束形成的领域,具体涉及一种非均匀布阵舷侧阵声纳波束形成快速方法。
背景技术
频域常规波束形成作为其信号处理流程必不可少的一个环节,是最消耗计算资源的处理模块之一。在工程实现过程中通常采用多块DSP板或多核CPU并行处理技术保证其实时性,这大大增加了经济成本和程序编写复杂度。舷侧阵声纳是潜艇最常用声纳之一,为了提高舷侧阵声纳的远程警戒和被动定位能力,其面积越来越大,阵元数也越来越多,同时在工程应用中,严格按照等间距布阵存在一定困难。这更加给常规频域波束形成的实现提出了巨大的挑战。
当前已有的快速算法主要有以下几种:
一是将导向矢量(以均匀线列阵为例)换为(其中f表示频率,d表示阵元间距,c表示声速,m表示阵元序号,θ表示波束角度,k表示频点序号,N表示波束数)(这里以端射方向为0度),将计算过程转换为FFT运算。这需要根据频率和阵元间距对每个频点的结果进行修正,对于宽带波束形成实现起来比较受限制,一般只适用于窄带。
二是基于FRFT的思想,通过引入布鲁斯坦公式将波束形成结果转化为分数傅里叶变换(FRFT)的形式,然后采用FFT来实现。这类算法实现过程没有具体描述,在使用FFT计算FRFT的过程中,如果直接采用与数据长度等长度点数的两次FFT和一次FFT-1来实现,由于没有考虑线性卷积与圆周卷积的差别,得出的结果与直接计算有一定的差别,对后面的波束收听和目标跟踪识别等处理产生不利影响。
三是首先将线阵频域波束形成每个频点的补相求和的计算过程转化为线性卷积运算,然后用两次FFT和一次IFFT(逆FFT)计算M(阵元数)点以后的结果,采用数据平移构造新的序列,用一次FFT和一次IFFT计算前M点的结果,将两次结果进行拼接得到全部方位的结果,最后将各频点的结果模的平方求和得到波束形成结果。但是该方法要求阵列严格按照等间距排列,并且波束数必须大于二倍阵元数。
发明内容
本发明的目的在于克服现有技术存在的不足,而提供一种非均匀布阵舷侧阵声纳波束形成快速方法。
本发明的目的是通过如下技术方案来完成的:这种非均匀布阵舷侧阵声纳波束形成快速方法,首先将频域波束形成每个频点的补相位求和的计算过程先转化为两个DFT卷积,即用两次DFT卷积计算频域波束形成每个频点的补相位求和过程,然后进行筛选做一个近似计算;包括以下步骤:
1)设舷侧阵为平面阵,水平维阵元数为n1,垂直维阵元数为n2,总阵元数N=n1×n2,水平维和垂直维均为非均匀布阵,每个阵元的坐标分别为(xn,yn,zn),n=0,1,…,N-1,其中n=0,1,…,n1-1为第一排阵元,n=n1,n1+1,,…,2n1-1为第二排阵元,以此类推;
3)第l个波束的第p个频点输出为:
式中k表示筛选点数,具有傅立叶变换的形式,进一步写成:
得到:
本发明的有益效果为:本发明首先将频域波束形成每个频点的补相位求和的计算过程转化为两个DFT卷积,然后进行筛选做一个近似计算,进行工程可接受的近似计算,计算量小;对阵元位置没有严格要求,克服了现有快速算法通常要均匀布阵的限制;形成的波束个数与阵元数无关,波束个数及波束指向可任意设置。
附图说明
图1为本发明的舷侧阵示意图。
图2为本发明和常规算法波束形成结果在方位角0°、频率1.5kHz时比较示意图。
图3为本发明的常规算法波束形成结果在方位角-45°、频率1.5kHz时比较示意图。
图4为本发明和常规算法波束形成结果在方位角0°、频率0.9kHz时比较示意图。
图5为本发明的常规算法波束形成结果在方位角-45°、频率0.9kHz时比较示意图。
图6为本发明和常规算法波束形成结果在方位角0°、频率0.3kHz时比较示意图。
图7为本发明的常规算法波束形成结果在方位角-45°、频率0.3kHz时比较示意图。
具体实施方式
下面将结合附图对本发明做详细的介绍:
如附图1所示,这种非均匀布阵舷侧阵声纳波束形成快速方法,主要包括以下步骤:
1)设舷侧阵为平面阵,水平维阵元数为n1,垂直维阵元数为n2,总阵元数N=n1×n2,水平维和垂直维均为非均匀布阵,每个阵元的坐标分别为(xn,yn,zn),n=0,1,…,N-1,其中n=0,1,…,n1-1为第一排阵元,n=n1,n1+1,,…,2n1-1为第二排阵元,以此类推;
3)第l个波束的第p个频点输出为:
式中k表示筛选点数,具有傅立叶变换的形式,进一步写成:
得到:
式中k2-k1+1<N,因此采用上式可大量节省计算量。筛选点数k1,k2的具体值与波束号、工作频率、阵元数量等有关,但这些值可提前计算好。
本发明计算过程如下:
3)计算出所有频点所有预成波束的补偿系数,并存储在数据库中;
常规算法与本发明快速算法性能比较分析:
这里常规算法指的是对每个频点的阵元数据补相位然后将各阵元补相位的结果求和。在波束形成计算的过程中复数的乘法是影响算法速度的关键,这里主要考虑复数乘法的影响。假设阵元数为N,波束数为L,频点数为p的频域常规波束形成复数乘法为p*N*L。快速算法增加了一次n点的FFT,故复数乘法次数为p*n/2*log2(n)+p*L*(k2-k1+1),可以看出快速算法的乘法次数与阵元数无关。对常规算法和快速算法进行比较,可得:
从上式可以看出,算法性能提高与频点数无关,在实际中还涉及复数加法和一些其他的运算,实际运算速度会有些差别。上式改写为:
从上式可以看出,该快速算法适合阵元数多或者波束数多基阵。
实施例:空间舷侧阵由四条线阵组成,每条线阵由140个阵元非均匀排列组成,总阵元数为560个,形状如附图1所示。
预成181个波束,快速算法中筛选点数取61,常规算法和快速算法不同角度不同频率下的波束形成结果如附图2所示,频率分别为1.5kHz、0.9kHz、0.3kHz,角度分别为0°、-45°。由附图2中可以看出,在不同角度快速算法和常规算法波束结果存在微小的误差,这是由于筛选点数进行近似计算造成的。两种方法的误差非常微小,而且它们的目标方位完全一致,主瓣宽度和旁瓣高度完全一致,不影响工程使用。对两种算法不同阵元数和不同波束数运算时间进行比较,其中运算次数按式进行计算。利用matlab 7.11.0进行仿真计算,得到运算时间。
表1不同阵元数,预成181个波束,频点数为1,两种方法运行次数及时间比较
表2不同波束数,阵元数为280(1/2阵),频点数为1,两种方法运行次数及时间比较
从表1、表2中可以看出,快速算法的运算次数和运算时间均少于常规算法。当阵元数成倍增加时,快速算法的运算速度明显优于常规算法。
可以理解的是,对本领域技术人员来说,对本发明的技术方案及发明构思加以等同替换或改变都应属于本发明所附的权利要求的保护范围。
Claims (1)
1.一种非均匀布阵舷侧阵声纳波束形成快速方法,其特征在于:首先将频域波束形成每个频点的补相位求和的计算过程先转化为两个DFT卷积,即用两次DFT卷积计算频域波束形成每个频点的补相位求和过程,然后进行筛选做一个近似计算;具体包括以下步骤:
1)设舷侧阵为平面阵,水平维阵元数为n1,垂直维阵元数为n2,总阵元数N=n1×n2,水平维和垂直维均为非均匀布阵,每个阵元的坐标分别为(xn,yn,zn),n=0,1,…,N-1,其中n=0,1,…,n1-1为第一排阵元,n=n1,n1+1,,…,2n1-1为第二排阵元,以此类推;
3)第l个波束的第p个频点输出为:
式中k为筛选点数,具有傅立叶变换的形式,进一步写成:
得到:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011464992.5A CN112711014B (zh) | 2020-12-14 | 2020-12-14 | 一种非均匀布阵舷侧阵声纳波束形成快速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011464992.5A CN112711014B (zh) | 2020-12-14 | 2020-12-14 | 一种非均匀布阵舷侧阵声纳波束形成快速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112711014A true CN112711014A (zh) | 2021-04-27 |
CN112711014B CN112711014B (zh) | 2022-11-01 |
Family
ID=75541869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011464992.5A Active CN112711014B (zh) | 2020-12-14 | 2020-12-14 | 一种非均匀布阵舷侧阵声纳波束形成快速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112711014B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101430380A (zh) * | 2008-12-19 | 2009-05-13 | 北京航空航天大学 | 基于非均匀采样的大斜视角机载sar聚束模式成像方法 |
WO2009112798A1 (en) * | 2008-02-18 | 2009-09-17 | Thomas Edgar Curtis | Underwater surveillance |
CN101813772A (zh) * | 2009-12-31 | 2010-08-25 | 中国科学院声学研究所 | 一种快速宽带频域扩展拖曳阵波束形成方法 |
CN102043145A (zh) * | 2010-11-03 | 2011-05-04 | 中国科学院声学研究所 | 基于声矢量传感器均匀直线阵的快速宽带频域波束形成方法 |
WO2012155990A1 (en) * | 2011-05-18 | 2012-11-22 | Lambda:4 Entwicklungen Gmbh | Method to determine the location of a receiver |
CN104678384A (zh) * | 2013-11-28 | 2015-06-03 | 中国科学院声学研究所 | 一种波束域的声压差互相关谱分析水下目标速度估计方法 |
CN104678376A (zh) * | 2013-11-28 | 2015-06-03 | 中国科学院声学研究所 | 一种扩展的拖曳阵宽带频域逆波束形成方法 |
CN108931766A (zh) * | 2018-04-28 | 2018-12-04 | 河海大学 | 一种基于稀疏重构的非均匀stap干扰目标滤除方法 |
CN109521401A (zh) * | 2018-09-27 | 2019-03-26 | 北京大学 | 一种合成孔径成像的快速波束形成方法 |
CN110515062A (zh) * | 2019-07-04 | 2019-11-29 | 中国船舶重工集团公司第七一五研究所 | 一种全艇多螺旋共形稀疏布阵方法 |
CN110632605A (zh) * | 2019-08-01 | 2019-12-31 | 中国船舶重工集团公司第七一五研究所 | 一种宽容的大孔径拖线阵时域单波束处理方法 |
-
2020
- 2020-12-14 CN CN202011464992.5A patent/CN112711014B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009112798A1 (en) * | 2008-02-18 | 2009-09-17 | Thomas Edgar Curtis | Underwater surveillance |
CN101430380A (zh) * | 2008-12-19 | 2009-05-13 | 北京航空航天大学 | 基于非均匀采样的大斜视角机载sar聚束模式成像方法 |
CN101813772A (zh) * | 2009-12-31 | 2010-08-25 | 中国科学院声学研究所 | 一种快速宽带频域扩展拖曳阵波束形成方法 |
CN102043145A (zh) * | 2010-11-03 | 2011-05-04 | 中国科学院声学研究所 | 基于声矢量传感器均匀直线阵的快速宽带频域波束形成方法 |
WO2012155990A1 (en) * | 2011-05-18 | 2012-11-22 | Lambda:4 Entwicklungen Gmbh | Method to determine the location of a receiver |
CN104678384A (zh) * | 2013-11-28 | 2015-06-03 | 中国科学院声学研究所 | 一种波束域的声压差互相关谱分析水下目标速度估计方法 |
CN104678376A (zh) * | 2013-11-28 | 2015-06-03 | 中国科学院声学研究所 | 一种扩展的拖曳阵宽带频域逆波束形成方法 |
CN108931766A (zh) * | 2018-04-28 | 2018-12-04 | 河海大学 | 一种基于稀疏重构的非均匀stap干扰目标滤除方法 |
CN109521401A (zh) * | 2018-09-27 | 2019-03-26 | 北京大学 | 一种合成孔径成像的快速波束形成方法 |
CN110515062A (zh) * | 2019-07-04 | 2019-11-29 | 中国船舶重工集团公司第七一五研究所 | 一种全艇多螺旋共形稀疏布阵方法 |
CN110632605A (zh) * | 2019-08-01 | 2019-12-31 | 中国船舶重工集团公司第七一五研究所 | 一种宽容的大孔径拖线阵时域单波束处理方法 |
Non-Patent Citations (3)
Title |
---|
P.WEBB ET AL: ""Analysis of fast accurate low ambiguity beam forming for non λ/2 ultrasonic arrays"", 《ULTRASONIC》 * |
刘辉: ""MIMO稀疏阵列SAR天线设计与三维成像技术研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
陈晓: ""基于潜标的目标探测信号处理算法研究及软件实现"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112711014B (zh) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108508423B (zh) | 基于异型阵的子阵数字和差单脉冲测角方法 | |
CN108344967B (zh) | 基于互质面阵的二维波达方向快速估计方法 | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN110197112B (zh) | 一种基于协方差修正的波束域Root-MUSIC方法 | |
CN109635240A (zh) | 大规模数字阵列零陷展宽自适应波束形成方法 | |
CN110346752B (zh) | 基于互质稀疏阵的无模糊测向方法 | |
CN109507643B (zh) | 旁瓣相消中展宽零陷的方法 | |
CN113311397B (zh) | 基于卷积神经网络的大型阵列快速自适应抗干扰方法 | |
CN108120953A (zh) | 一种基于波达方向估计的无线电定位方法 | |
CN110736976B (zh) | 一种任意阵形的声纳波束形成器性能估计方法 | |
CN108490428B (zh) | 一种用于抗主瓣干扰的降维子阵比相跟踪测角方法 | |
CN105572629B (zh) | 一种适用于任意阵列结构的低运算复杂度的二维测向方法 | |
CN109541526A (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN112711014B (zh) | 一种非均匀布阵舷侧阵声纳波束形成快速方法 | |
CN111368256A (zh) | 一种基于均匀圆阵的单快拍测向方法 | |
CN111722178A (zh) | 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN112799008B (zh) | 一种声速无关的快速二维波达方向估计方法 | |
CN115470446A (zh) | 一种适用于多波束测深系统的任意角度快速波束形成方法 | |
CN111366891B (zh) | 一种基于伪协方差矩阵的均匀圆阵单快拍测向方法 | |
CN111830458A (zh) | 一种平行线阵单快拍二维测向方法 | |
CN110907888A (zh) | 一种基于波束形成的快速极化-doa估计方法 | |
Chitre et al. | Optimisation and beamforming of a two dimensional sparse array | |
Weiß et al. | Single snapshot DoA estimation from a Rotman lens using machine learning techniques | |
CN112068071B (zh) | 一种用于任意子阵排布的波束域波达方向矩阵法 |
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 |