CN110687538A - 一种基于近场聚焦的超波束形成方法 - Google Patents
一种基于近场聚焦的超波束形成方法 Download PDFInfo
- Publication number
- CN110687538A CN110687538A CN201911031089.7A CN201911031089A CN110687538A CN 110687538 A CN110687538 A CN 110687538A CN 201911031089 A CN201911031089 A CN 201911031089A CN 110687538 A CN110687538 A CN 110687538A
- Authority
- CN
- China
- Prior art keywords
- focusing
- array
- super
- beam forming
- phase
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000003384 imaging method Methods 0.000 claims abstract description 11
- 230000010363 phase shift Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 4
- 238000003780 insertion Methods 0.000 claims description 3
- 230000037431 insertion Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
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
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
-
- 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/523—Details of pulse systems
- G01S7/526—Receivers
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
本发明公开一种基于近场聚焦的超波束形成方法,属于声纳信号处理技术领域。接收阵元域数据;计算聚焦环系数;将聚焦环系数应用到超波束形成中;对波束形成处理之后的波束域数据进行成像处理。本发明具有超波束窄主瓣、低旁瓣的优点,可以有效地提高声纳***的检测性能,降低虚警概率,提高目标方位的估计精度;采用分段聚焦超波束形成技术,在近场范围内对相位进行聚焦补偿,弥补了相位失配造成的图像模糊的问题。
Description
技术领域
本发明涉及声纳信号处理技术领域,特别涉及一种基于近场聚焦的超波束形成方法。
背景技术
高频多波束声纳是目前海洋目标探测的重要仪器,主要应用于海洋中近距离小目标的探测、定位和分类。但是图像声纳在近场工作区,存在着成像模糊等问题,严重影响设备的工作性能。这是因为传统的波束形成方法是以平面波为假设的,没有对近场范围进行聚焦补偿,造成相位失配,波束主瓣将变宽,旁瓣升高,致使方位分辨率下降,造成声纳图像性能的严重恶化。同时传统的常规波束形成主瓣较宽、旁瓣较高,不利于多目标的分辨识别。
发明内容
本发明的目的在于提供一种基于近场聚焦的超波束形成方法,以解决传统波束形成在近场工作区,成像模糊的问题。
为解决上述技术问题,本发明提供一种基于近场聚焦的超波束形成方法,包括:
接收阵元域数据;
计算聚焦环系数;
将聚焦环系数应用到超波束形成中;
对波束形成处理之后的波束域数据进行成像处理。
可选的,计算聚焦环系数包括:
声聚焦相移波束形成的相移是按球面波规律计算的,对于空间某一点射出的声波进行声聚焦,能够形成左、右阵声聚焦;
计算聚焦点到达接收阵各基元的声程差或相位差,然后在接收处理时补偿声程差或相位差,补偿后形成对接收点传播过来的信号实现同相叠加,使声基阵聚焦于空间该聚焦点;
左半阵相位差φi表达式为:
其中,R0是聚焦距离,λ为波长,N是基元总数,i是基元序号,d是基元间隔,θk是入射角;
在相邻阵元间***相移角ψk,则与各阵元对应的输出为:
t是时间,k是波束号,j是虚数单位,ω0=2*π*f*t,f是中心频率,φih是相移,ψik是相位补偿,其中-ψik=φik,则上式改写为:
主动声纳通过发射换能器向水中发射脉冲信号,经目标反射,利用接收回波与发射脉冲信号的时间差来测量距离,往返的传播时间为:
式中,R为目标距离,c为声速;
回波信号经采样后,采样点N与距离的关系为:
fs是采样率,根据上式计算不同聚焦距离的相位补偿系数的***位置。
可选的,将聚焦环系数应用到超波束形成中包括:
将接收线阵均等分为左、右两个子阵,左、右子阵分别按照基带常规波束形成方法进行波束形成,得到和波束SL和差波束SR,则:
和波束SS=|SL|+|SR|
差波束SD=|SL-SR|
将和波束与差波束相减,加权计算得到超波束SH输出为:
式中n为超波束指数,范围为0.3≤n≤1,用以调整输出波束主瓣宽度和旁瓣高度。
可选的,采用近场模型计算波束,进行超波束形成,采用分段聚焦算法,大幅度降低旁瓣的同时,压缩主瓣宽度。
在本发明中提供了一种基于近场聚焦的超波束形成方法,接收阵元域数据;计算聚焦环系数;将聚焦环系数应用到超波束形成中;对波束形成处理之后的波束域数据进行成像处理。
本发明具有以下有益效果:
(1)具有超波束窄主瓣、低旁瓣的优点,可以有效地提高声纳***的检测性能,降低虚警概率,提高目标方位的估计精度;
(2)采用分段聚焦超波束形成技术,在近场范围内对相位进行聚焦补偿,弥补了相位失配造成的图像模糊的问题;
(3)有效提高声纳对空间中不同方位多个目标的分辨能力,在不损失聚焦超波束形成精度的同时,具有快速计算能力,结构简单、易于实现的优点。
附图说明
图1是本发明提供的基于近场聚焦的超波束形成方法流程示意图;
图2是左、右阵声聚焦示意图;
图3是右波束聚焦A点示意图;
图4是超波束形成处理方法示意图;
图5是进行常规波束形成、未聚焦的效果图;
图6是进行超波束形成、近场聚焦的效果图。
具体实施方式
以下结合附图和具体实施例对本发明提出的一种基于近场聚焦的超波束形成方法作进一步详细说明。根据下面说明和权利要求书,本发明的优点和特征将更清楚。需说明的是,附图均采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施例的目的。
实施例一
本发明提供了一种基于近场聚焦的超波束形成方法,其流程如图1所示,包括:
接收阵元域数据;
计算聚焦环系数;
将聚焦环系数应用到超波束形成中;
对波束形成处理之后的波束域数据进行成像处理。
具体的,步骤一:接收数据:
接收阵元域数据;
步骤二:计算聚焦环系数:
近场条件为:
其中,L为基阵的有效阵长,λ为波长,D为近场范围的距离。
声聚焦相移波束形成的相移是按球面波规律计算的,对于空间某一点射出的声波进行声聚焦,能够形成左、右阵声聚焦,如图2所示;图2中 N是基元数,o是基阵几何中心,d是基元间隔,Fo是法线,右波束焦距 Ao=Ro右,左波束焦距Bo=Ro左,Ao与Fo,Bo与Fo夹角为左、右波束指向角,1<i<j<N。
首先计算聚焦点A点(或B点)声波到达接收阵各基元的声程差或相位差,然后在接收处理时补偿这些声程差或相位差,经过补偿后,形成对接收点传播过来的信号实现同相叠加,使声基阵聚焦于空间A点(或B点)。下面介绍球面波相移补偿角计算公式。为了分析简化,以右波束为例,如图3所示。
在图3中,右波束焦距为Ro,利用三角形勾股定理计算Ai、Aj:聚焦点A到基元i、j的距离Ai、Aj,相应声程Ai-Ro,Aj-Ro和左、右半阵相位差φi、φj。推算结果:φi和φj表达式是完全一样的,相位差表达式为:
其中,R0是聚焦距离,λ为波长,N是基元总数,i是基元序号,d是基元间隔,θk是入射角;
在相邻阵元间***相移角ψk,则与各阵元对应的输出为:
t是时间,k是波束号,j是虚数单位,ω0=2*π*f*t,f是中心频率,φih是相移,ψik是相位补偿,其中-ψik=φik,则上式改写为:
φik是相位差,式(1.4)为球面波相移补偿角的计算公式。
由公式(1.4)可知,近场环境下聚焦波束形成的相位补偿系数是信号入射角θ和声源与阵元中心的距离R的二维函数。对单点聚焦只能保证焦点附件区域的波束性能,其它位置上的目标未得到精确补偿,波束性能将会恶化。因此在距离上,将整个近场区域分成若干段,对每段距离的中点分别进行定焦补偿。分段越多,波束误差越小,聚焦性能越好,而计算量越大,实现越困难。在实际工程中需要考虑计算量和聚焦性能的平衡。
主动声纳通过发射换能器向水中发射脉冲信号,经目标反射,利用接收回波与发射脉冲信号的时间差来测量距离,往返的传播时间为:
式中,R为目标距离,c为声速;
回波信号经采样后,采样点N与距离的关系为:
fs是采样率,根据公式(1.6)计算不同聚焦距离的相位补偿系数的***位置。
步骤三:聚焦超波束形成:
将步骤二中得到的聚焦环系数应用到下面的超波束形成中。
超波束形成处理方法如图4所示,将接收线阵均等分为左、右两个子阵,左、右子阵分别按照基带常规波束形成方法进行波束形成,得到和波束SL和差波束SR,根据超波束算法的定义,
和波束SS=|SL|+|SR| (1.7)
差波束SD=|SL-SR| (1.8)
将和波束与差波束相减,加权计算得到超波束SH输出为:
式中n为超波束指数,范围为0.3≤n≤1,用以调整输出波束主瓣宽度和旁瓣高度。
步骤四:发送数据:
通过上述步骤之后,得到波束形成处理之后的波束域数据,进行成像处理。图5采用远场模型计算波束,进行常规波束形成,未采用聚焦算法,可以看出,较高的旁瓣造成近场成像效果恶化,目标越近,成像效果越差,并且目标距离出现了偏移,定位不准确。图6采用近场模型计算波束,进行超波束形成,采用分段聚焦算法,大幅度降低旁瓣的同时,压缩了主瓣宽度,提高声纳对空间中多个目标的分辨能力,有效地提高***的检测性能,降低虚警概率,并提高目标方位的估计精度,成像效果最优。
上述描述仅是对本发明较佳实施例的描述,并非对本发明范围的任何限定,本发明领域的普通技术人员根据上述揭示内容做的任何变更、修饰,均属于权利要求书的保护范围。
Claims (4)
1.一种基于近场聚焦的超波束形成方法,其特征在于,包括:
接收阵元域数据;
计算聚焦环系数;
将聚焦环系数应用到超波束形成中;
对波束形成处理之后的波束域数据进行成像处理。
2.如权利要求1所述的基于近场聚焦的超波束形成方法,其特征在于,计算聚焦环系数包括:
声聚焦相移波束形成的相移是按球面波规律计算的,对于空间某一点射出的声波进行声聚焦,能够形成左、右阵声聚焦;
计算聚焦点到达接收阵各基元的声程差或相位差,然后在接收处理时补偿声程差或相位差,补偿后形成对接收点传播过来的信号实现同相叠加,使声基阵聚焦于空间该聚焦点;
左半阵相位差φi表达式为:
其中,R0是聚焦距离,λ为波长,N是基元总数,i是基元序号,d是基元间隔,θk是入射角;
在相邻阵元间***相移角ψk,则与各阵元对应的输出为:
t是时间,k是波束号,j是虚数单位,ω0=2*π*f*t,f是中心频率,φih是相移,ψik是相位补偿,其中-ψik=φik,则上式改写为:
φik是相位差,上式为球面波相移补偿角的计算公式;
主动声纳通过发射换能器向水中发射脉冲信号,经目标反射,利用接收回波与发射脉冲信号的时间差来测量距离,往返的传播时间为:
式中,R为目标距离,c为声速;
回波信号经采样后,采样点N与距离的关系为:
fs是采样率,根据上式计算不同聚焦距离的相位补偿系数的***位置。
4.如权利要求3所述的基于近场聚焦的超波束形成方法,其特征在于,采用近场模型计算波束,进行超波束形成,采用分段聚焦算法,大幅度降低旁瓣的同时,压缩主瓣宽度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911031089.7A CN110687538A (zh) | 2019-10-28 | 2019-10-28 | 一种基于近场聚焦的超波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911031089.7A CN110687538A (zh) | 2019-10-28 | 2019-10-28 | 一种基于近场聚焦的超波束形成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110687538A true CN110687538A (zh) | 2020-01-14 |
Family
ID=69114256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911031089.7A Pending CN110687538A (zh) | 2019-10-28 | 2019-10-28 | 一种基于近场聚焦的超波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110687538A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112729742A (zh) * | 2020-12-25 | 2021-04-30 | 西安理工大学 | 一种基于概率成像的冲击定位方法 |
CN113030983A (zh) * | 2021-03-17 | 2021-06-25 | 中国科学院声学研究所 | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 |
CN113391301A (zh) * | 2021-05-26 | 2021-09-14 | 南京理工大学 | 基于时间调制线性阵列的超波束形成方法及*** |
CN115902853A (zh) * | 2023-03-09 | 2023-04-04 | 中国科学院深海科学与工程研究所 | 一种适用高速海底测绘的合成接收孔径聚焦波束形成方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104215965A (zh) * | 2014-09-05 | 2014-12-17 | 中国船舶重工集团公司第七一五研究所 | 高频声纳近场聚焦距离的分档方法 |
CN109116334A (zh) * | 2018-10-19 | 2019-01-01 | 海鹰企业集团有限责任公司 | 基于超波束加权的声纳波束形成方法及*** |
CN109283536A (zh) * | 2018-09-01 | 2019-01-29 | 哈尔滨工程大学 | 一种多波束测深声呐水体成像波束形成算法 |
CN109581388A (zh) * | 2018-12-20 | 2019-04-05 | 华中科技大学 | 一种实时三维成像声纳的近场宽视角波束形成方法 |
CN110235022A (zh) * | 2016-12-15 | 2019-09-13 | 达尔豪斯大学 | 用于使用相干复合菲涅尔聚焦进行超声波束形成的***和方法 |
-
2019
- 2019-10-28 CN CN201911031089.7A patent/CN110687538A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104215965A (zh) * | 2014-09-05 | 2014-12-17 | 中国船舶重工集团公司第七一五研究所 | 高频声纳近场聚焦距离的分档方法 |
CN110235022A (zh) * | 2016-12-15 | 2019-09-13 | 达尔豪斯大学 | 用于使用相干复合菲涅尔聚焦进行超声波束形成的***和方法 |
CN109283536A (zh) * | 2018-09-01 | 2019-01-29 | 哈尔滨工程大学 | 一种多波束测深声呐水体成像波束形成算法 |
CN109116334A (zh) * | 2018-10-19 | 2019-01-01 | 海鹰企业集团有限责任公司 | 基于超波束加权的声纳波束形成方法及*** |
CN109581388A (zh) * | 2018-12-20 | 2019-04-05 | 华中科技大学 | 一种实时三维成像声纳的近场宽视角波束形成方法 |
Non-Patent Citations (2)
Title |
---|
张昌: "图像声纳的聚焦波束形成技术研究及实现", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
杨东方, 《海洋出版社》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112729742A (zh) * | 2020-12-25 | 2021-04-30 | 西安理工大学 | 一种基于概率成像的冲击定位方法 |
CN112729742B (zh) * | 2020-12-25 | 2023-08-22 | 西安理工大学 | 一种基于概率成像的冲击定位方法 |
CN113030983A (zh) * | 2021-03-17 | 2021-06-25 | 中国科学院声学研究所 | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 |
CN113030983B (zh) * | 2021-03-17 | 2021-12-28 | 中国科学院声学研究所 | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 |
CN113391301A (zh) * | 2021-05-26 | 2021-09-14 | 南京理工大学 | 基于时间调制线性阵列的超波束形成方法及*** |
CN115902853A (zh) * | 2023-03-09 | 2023-04-04 | 中国科学院深海科学与工程研究所 | 一种适用高速海底测绘的合成接收孔径聚焦波束形成方法 |
CN115902853B (zh) * | 2023-03-09 | 2023-05-12 | 中国科学院深海科学与工程研究所 | 一种适用高速海底测绘的合成接收孔径聚焦波束形成方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110687538A (zh) | 一种基于近场聚焦的超波束形成方法 | |
CN107831480B (zh) | 弹载雷达和差通道稳健自适应杂波抑制方法 | |
CN109581388B (zh) | 一种实时三维成像声纳的近场宽视角波束形成方法 | |
CN109765521B (zh) | 一种基于子阵划分的波束域成像方法 | |
CN109725285B (zh) | 一种基于mvdr协方差矩阵元素自适应相角转换的doa估计方法 | |
Zhang et al. | Influence of the stop-and-hop assumption on synthetic aperture sonar imagery | |
CN109557526B (zh) | 一种基于压缩感知理论的矢量水听器稀疏布阵方法 | |
Chi et al. | High-resolution real-time underwater 3-D acoustical imaging through designing ultralarge ultrasparse ultra-wideband 2-D arrays | |
CN103513238B (zh) | 一种规整化最小二乘子空间相交的目标方位测向方法 | |
CN102846336B (zh) | 超声成像延时控制方法及装置 | |
JP2533287B2 (ja) | 線形アレイ横方向運動補正方法 | |
CN109491009B (zh) | 一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法 | |
CN105824016B (zh) | 运动平台雷达检测超低空目标的稳健空时自适应处理方法 | |
CN113126030B (zh) | 基于宽带声场干涉结构的深海直达声区目标深度估计方法 | |
CN111983599A (zh) | 一种基于方位-俯仰字典的目标二维doa估计方法 | |
CN109375197B (zh) | 一种小尺寸矢量阵低频散射校正方法 | |
CN116679306A (zh) | 一种基于系数修正的扩展孔径声纳成像方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN113820654B (zh) | 基于波束域降维的s波段雷达目标低仰角doa估计方法 | |
CN111431575B (zh) | 基于常规波束形成的来波方向稀疏重构方法 | |
CN113075633A (zh) | 一种基于距离-俯仰联合字典的目标定位方法 | |
CN109029387B (zh) | 一种波束内拟合多波束测深算法 | |
CN112711018A (zh) | 一种近场空间网状结构聚焦波束形成方法 | |
CN111551943A (zh) | 一种基于damas2的稀疏布阵高分辨三维声成像方法及*** | |
CN110244288A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200114 |