CN112099017A - 一种基于过扫描的圆环阵低旁瓣波束优化方法 - Google Patents
一种基于过扫描的圆环阵低旁瓣波束优化方法 Download PDFInfo
- Publication number
- CN112099017A CN112099017A CN202010890410.3A CN202010890410A CN112099017A CN 112099017 A CN112099017 A CN 112099017A CN 202010890410 A CN202010890410 A CN 202010890410A CN 112099017 A CN112099017 A CN 112099017A
- Authority
- CN
- China
- Prior art keywords
- overscan
- amount
- sidelobe
- side lobe
- level
- 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/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
-
- 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/006—Theoretical aspects
-
- 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
-
- 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/534—Details of non-pulse systems
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)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明涉及一种基于过扫描的圆环阵低旁瓣波束优化方法,通过给定过扫描量,在波束成形操作中引入额外的延迟,使得波束响应旁瓣级降低,既可以实现给定旁瓣级情况时的波束优化,也给出了最优过扫描量选取的指标函数。本发明较二阶锥优化方法计算简单方便,相对于圆环阵MVDR方法旁瓣级有明显降低,且仍具有超指向性。获取的低旁瓣波束与相应旁瓣约束下二阶锥优化得到的波束也较为接近。存在过扫描量时,波束图较不存在过扫描量时的旁瓣级明显下降。通过选取合适的过扫描量,可以实现给定的旁瓣级约束,显著降低旁瓣级,达到预期的性能。选取最优过扫描量时,可以实现在最小水平波束和垂直波束的主瓣展宽下使得水平波束旁瓣级最低。
Description
技术领域
本发明属于声学阵列信号处理、语音信号处理和声呐技术等领域,涉及一种基于过扫描的圆环阵低旁瓣波束优化方法,特别是涉及圆环形传感器阵列的波束优化方法。
背景技术
在信号处理问题中,波束图优化设计包括两方面研究内容,一方面是控制波束旁瓣,另一方面是设计波束主瓣响应。常用的时延求和方法是对各个通道的信号进行时延补偿,然后相加,使得各个阵元信号同相相加,从而提高信噪比,但是该方法提供的空间指向性十分有限。而对于MVDR波束形成方法得到的波束响应,虽然获得了最大指向性因子,但不能控制旁瓣。
现有的旁瓣控制波束形成方法多需要复杂的算法对旁瓣进行约束求解。公布号为CN 106682405A的专利公开了一种基于凸优化的低旁瓣波束图综合设计方法,其利用交替方向乘子法和凸优化约束来对旁瓣进行控制,计算相对复杂。文献2“任意几何形状和阵元指向性的传感器阵列优化波束形成方法,声学学报,2005,vol.30(3),p.264-270”公开了一种基于二阶锥规划(SOCP:Second-Order Cone Programming)的白噪声增益约束方法,能得到较稳健的低旁瓣波束,但同样使用了凸优化的方法。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于过扫描的圆环阵低旁瓣波束优化方法,解决现有的低旁瓣波束优化方法需要二阶锥优化等复杂算法的问题,通过给定过扫描量,在波束成形操作中引入额外的延迟,使得波束响应旁瓣级降低,既可以实现给定旁瓣级情况时的波束优化,也给出了最优过扫描量选取的指标函数,过程简单,计算方便。
技术方案
一种基于过扫描的圆环阵低旁瓣波束优化方法,其特征在于步骤如下:
步骤1、计算基于过扫描的波束权值向量:
wσ(θ0,φ0)=Tωσ(θ0,φ0)
T=[v0,v1,…,vM-1]
ωσ(θ0,φ0)=[ω0,σ(θ0,φ0),ω1,σ(θ0,φ0),…,ωM-1,σ(θ0,φ0)]T
vm=M-1/2[1,eimβ,…,ei(M-1)mβ]T
其中:σ是过扫描量,M为阵元数,β是相邻阵元间夹角,(θ0,φ0)为波束指向角,上标T表示转置;
步骤2:阵列流形向量
P(θ,φ)=[p0(θ,φ),p1(θ,φ),…,pM-1(θ,φ)]T
合成波束由下式计算:
步骤3、过扫描量的选取:
1、当旁瓣级给定时,根据下述旁瓣级和过扫描量的函数关系,选取此时对应的过扫描量:
SLh,σ=max(Bσ(Ω)|Ω∈ΩSL)
ΩSL={(θ,φ)|θ=θ0,φ∈[0,φ0-ΔSL]∪[φ0+ΔSL,360]}
其中:Ω=(θ,φ),ΔSL为旁瓣级波束宽度即波束主瓣功率下降到与旁瓣级相等时的两方向间夹角的一半;
2、最优过扫描量的选取:根据函数Q衡量不同过扫描量取值时的合成波束Bσ性能,以函数Q取最大值时,对应的过扫描量σop即为最优过扫描量;
其中:SLh,σ=0是过扫描量σ=0时,水平方向上合成波束的最高旁瓣峰值;
ΔSLh,σ=|SLh,σ-SLh,σ=0|是水平方向上,不同过扫描量取值下与扫描量为零时合成波束的最高旁瓣峰值变化量;
MLh,σ=0和MLe,σ=0分别是水平方向和垂直方向上,过扫描量σ=0时合成波束的-3dB主瓣宽度;
MLh,σ和MLe,σ分别是水平方向和垂直方向上,不同过扫描量σ取值时合成波束的-3dB主瓣宽度;
步骤4:合成波束,并计算对应的指向性因数和白噪声增益,将步骤3所得过扫描量σ代入加权向量wσ和合成波束Bσ;
其中:ρn是声透明圆环阵三维归一化噪声互谱矩阵
第(m,m')号元素为:ρmm'=ρs=sinc(k·Δrs)
Δrs=2asin(sβ/2),s=|m-m'|。
有益效果
本发明提出的一种基于过扫描的圆环阵低旁瓣波束优化方法,通过给定过扫描量,在波束成形操作中引入额外的延迟,使得波束响应旁瓣级降低,既可以实现给定旁瓣级情况时的波束优化,也给出了最优过扫描量选取的指标函数,过程简单,计算方便。本发明较二阶锥优化方法计算简单方便,相对于圆环阵MVDR方法旁瓣级有明显降低,且仍具有超指向性。本文获取的低旁瓣波束与相应旁瓣约束下二阶锥优化得到的波束也较为接近。存在过扫描量时,波束图较不存在过扫描量时的旁瓣级明显下降。通过选取合适的过扫描量,可以实现给定的旁瓣级约束,显著降低旁瓣级,达到预期的性能。选取最优过扫描量时,可以实现在最小水平波束和垂直波束的主瓣展宽下使得水平波束旁瓣级最低。
附图说明
图1为不同过扫描量取值对应的水平方向上合成波束旁瓣级;
图2为不同过扫描量取值对应的指标函数Q;
图3为过扫描量σ=0时的三维波束图;
图4为-20dB旁瓣约束下,过扫描量取σ=0.26和此时二阶锥优化的三维波束图;
图5为-20dB旁瓣约束下,过扫描量σ=0、过扫描量σ=0.26和此时二阶锥优化得到的水平波束和垂直波束;
图6为最优过扫描量σ=0.31和此时二阶锥优化的三维波束图;
图7为过扫描量σ=0、最优过扫描量σ=0.31和此时二阶锥优化得到的水平波束和垂直波束;
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明包括以下步骤:
步骤1:计算基于过扫描的波束权值向量:
wσ(θ0,φ0)=Tωσ(θ0,φ0)
其中:wσ(θ0,φ0)=[w0,σ(θ0,φ0),w1,σ(θ0,φ0),…,wM-1,σ(θ0,φ0)]T为阵元域权值向量,(θ0,φ0)为波束指向角,ωσ(θ0,φ0)=[ω0,σ(θ0,φ0),ω1,σ(θ0,φ0),…,ωM-1,σ(θ0,φ0)]T,ωσ(θ0,φ0)的第m个元素是:
α是使指向(θ0,φ0)方向的波束响应为1的参数,由下式计算:
所述:
ρs=sinc(k·Δrs)
T=[v0,v1,…,vM-1]
vm=M-1/2[1,eimβ,…,ei(M-1)mβ]T
其中σ是过扫描量,a为圆环阵半径,k=2π/λ为波数,λ是信号波长。φs是s号阵元方位角,M为阵元数,β是相邻阵元间夹角,φs=sβ,β=2π/M,Δrs是m和m'号阵元的间距,Δrs=2asin(sβ/2),s=|m-m'|。|·|表示求模值,上标T表示转置,上标H表示共轭转置。
步骤2:当加权向量为wσ,阵列流形向量P(θ,φ)=[p0(θ,φ),p1(θ,φ),…,pM-1(θ,φ)]T,合成波束由下式计算:
其中P(θ,φ)的第m个元素为:
步骤3:过扫描量的选取:
1、给定旁瓣级情况时的波束优化
旁瓣级与过扫描量的关系可由下式计算:
SLh,σ=max(Bσ(Ω)|Ω∈ΩSL)
ΩSL={(θ,φ)|θ=θ0,φ∈[0,φ0-ΔSL]∪[φ0+ΔSL,360]}
其中Ω=(θ,φ),ΔSL为旁瓣级波束宽度(即波束主瓣功率下降到与旁瓣级相等时的两方向间夹角)的一半,对于不同的过扫描量,ΔSL可取不同值。
当旁瓣级给定时,通过旁瓣级和过扫描量的函数关系,可以选取此时对应的过扫描量。
旁瓣级约束为-20dB,参照图1,此时对应的过扫描量σ=0.26。
2、最优过扫描量的选取
构建函数Q,来衡量不同过扫描量取值时的合成波束Bσ性能。
其中SLh,σ=0是过扫描量σ=0时,水平方向上合成波束的最高旁瓣峰值。
ΔSLh,σ=|SLh,σ-SLh,σ=0|是水平方向上,不同过扫描量取值下与扫描量为零时合成波束的最高旁瓣峰值变化量。
MLh,σ=0和MLe,σ=0分别是水平方向和垂直方向上,过扫描量σ=0时合成波束的-3dB主瓣宽度。
MLh,σ和MLe,σ分别是水平方向和垂直方向上,不同过扫描量σ取值时合成波束的-3dB主瓣宽度。
过扫描量σ取不同值时的-3dB主瓣宽度可由下式计算:
当函数Q取最大值时,对应的过扫描量σop即为最优过扫描量。
参照图2,当函数Q取最大值时,对应的过扫描量σop=0.31即为最优过扫描量。σ∈[0,0.4],均匀取值41个离散点。
步骤4:合成波束,并计算对应的指向性因数和白噪声增益。将步骤3所得过扫描量σ代入步骤1和2可直接得到加权向量wσ和合成波束Bσ。
计算过扫描量σ=0.26,σ=0.31时的加权向量wσ和合成波束Bσ:
此时白噪声增益由下式计算:
指向性因数由下式计算:
其中ρn是声透明圆环阵三维归一化噪声互谱矩阵,有:
其中第(m,m')号元素为:
ρmm'=ρs=sinc(k·Δrs)
式中s=|m-m'|。
当σ=0,即不存在过扫描量时为最优超指向性波束形成,此时得到的合成波束指向性指数为14.42dB,白噪声增益为-63.31dB,三维波束图参照图3。
在-20dB旁瓣约束下,过扫描量取值σ=0.26。此时过扫描波束优化和二阶锥优化得到的合成波束指向性指数分别为13.09dB、13.26dB,白噪声增益分别为-56.03dB、-59.16dB。由图4和图5可得,旁瓣级均为-20dB,满足约束条件。且过扫描波束优化得到的垂直波束旁瓣级低于二阶锥优化得到的垂直波束的旁瓣级。
最优过扫描量σ=0.31和二阶锥优化得到的合成波束指向性指数分别为12.50dB,12.73dB,白噪声增益分别为-53.98dB,-56.49dB。由图6和图7可以看出,本文优化方法得到的合成波束与二阶锥优化得到的合成波束在波束图上相接近,且过扫描波束优化得到的垂直波束旁瓣级低于二阶锥优化得到的垂直波束的旁瓣级,且此时指向性指数和白噪声增益大小近似。
本文提出的基于过扫描低旁瓣波束优化方法过程简单,计算方便,且在一定范围内可以得到与二阶锥优化方法近似的合成波束。
Claims (1)
1.一种基于过扫描的圆环阵低旁瓣波束优化方法,其特征在于步骤如下:
步骤1、计算基于过扫描的波束权值向量:
wσ(θ0,φ0)=Tωσ(θ0,φ0)
T=[v0,v1,…,vM-1]
ωσ(θ0,φ0)=[ω0,σ(θ0,φ0),ω1,σ(θ0,φ0),…,ωM-1,σ(θ0,φ0)]T
vm=M-1/2[1,eimβ,…,ei(M-1)mβ]T
其中:σ是过扫描量,M为阵元数,β是相邻阵元间夹角,(θ0,φ0)为波束指向角,上标T表示转置;
步骤2:阵列流形向量
P(θ,φ)=[p0(θ,φ),p1(θ,φ),...,pM-1(θ,φ)]T
合成波束由下式计算:
步骤3、过扫描量的选取:
1、当旁瓣级给定时,根据下述旁瓣级和过扫描量的函数关系,选取此时对应的过扫描量:
SLh,σ=max(Bσ(Ω)|Ω∈ΩSL)
ΩSL={(θ,φ)|θ=θ0,φ∈[0°,φ0-ΔSL]∪[φ0+ΔSL,360°]}
其中:Ω=(θ,φ),ΔSL为旁瓣级波束宽度即波束主瓣功率下降到与旁瓣级相等时的两方向间夹角的一半;
2、最优过扫描量的选取:根据函数Q衡量不同过扫描量取值时的合成波束Bσ性能,以函数Q取最大值时,对应的过扫描量σop即为最优过扫描量;
其中:SLh,σ=0是过扫描量σ=0时,水平方向上合成波束的最高旁瓣峰值;
ΔSLh,σ=|SLh,σ-SLh,σ=0|是水平方向上,不同过扫描量取值下与扫描量为零时合成波束的最高旁瓣峰值变化量;
MLh,σ=0和MLe,σ=0分别是水平方向和垂直方向上,过扫描量σ=0时合成波束的-3dB主瓣宽度;
MLh,σ和MLe,σ分别是水平方向和垂直方向上,不同过扫描量σ取值时合成波束的-3dB主瓣宽度;
步骤4:合成波束,并计算对应的指向性因数和白噪声增益,将步骤3所得过扫描量σ代入加权向量wσ和合成波束Bσ;
其中:ρn是声透明圆环阵三维归一化噪声互谱矩阵
第(m,m')号元素为:ρmm'=ρs=sinc(k·Δrs)
Δrs=2asin(sβ/2),s=|m-m'|。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010890410.3A CN112099017B (zh) | 2020-08-29 | 2020-08-29 | 一种基于过扫描的圆环阵低旁瓣波束优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010890410.3A CN112099017B (zh) | 2020-08-29 | 2020-08-29 | 一种基于过扫描的圆环阵低旁瓣波束优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112099017A true CN112099017A (zh) | 2020-12-18 |
CN112099017B CN112099017B (zh) | 2022-10-14 |
Family
ID=73758277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010890410.3A Active CN112099017B (zh) | 2020-08-29 | 2020-08-29 | 一种基于过扫描的圆环阵低旁瓣波束优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112099017B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070109179A1 (en) * | 2005-11-15 | 2007-05-17 | Werntz Paul C | Monostatic radar beam optimization |
CN104035064A (zh) * | 2014-06-05 | 2014-09-10 | 哈尔滨工程大学 | 适用于任意阵型的稳健宽带导向最小方差波束形成方法 |
CN106019236A (zh) * | 2016-05-24 | 2016-10-12 | 南京理工大学 | 一种基于数据重构的稀布阵数字波束形成方法 |
CN107271989A (zh) * | 2017-05-27 | 2017-10-20 | 西北工业大学 | 一种具有最优旋转角的低波束旁瓣mimo声纳阵型 |
CN111430915A (zh) * | 2020-04-29 | 2020-07-17 | 电子科技大学 | 一种基于方向图重构单元的阵列波束赋形方法 |
CN111551923A (zh) * | 2020-05-27 | 2020-08-18 | 电子科技大学 | 一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 |
-
2020
- 2020-08-29 CN CN202010890410.3A patent/CN112099017B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070109179A1 (en) * | 2005-11-15 | 2007-05-17 | Werntz Paul C | Monostatic radar beam optimization |
CN104035064A (zh) * | 2014-06-05 | 2014-09-10 | 哈尔滨工程大学 | 适用于任意阵型的稳健宽带导向最小方差波束形成方法 |
CN106019236A (zh) * | 2016-05-24 | 2016-10-12 | 南京理工大学 | 一种基于数据重构的稀布阵数字波束形成方法 |
CN107271989A (zh) * | 2017-05-27 | 2017-10-20 | 西北工业大学 | 一种具有最优旋转角的低波束旁瓣mimo声纳阵型 |
CN111430915A (zh) * | 2020-04-29 | 2020-07-17 | 电子科技大学 | 一种基于方向图重构单元的阵列波束赋形方法 |
CN111551923A (zh) * | 2020-05-27 | 2020-08-18 | 电子科技大学 | 一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 |
Non-Patent Citations (2)
Title |
---|
王英民等: "一种基于混合优化算法的阵列波束形成方法", 《哈尔滨工程大学学报》 * |
鄢社锋等: "任意几何形状和阵元指向性的传感器阵列优化波束形成方法", 《声学学报(中文版)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112099017B (zh) | 2022-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Robust high-order superdirectivity of circular sensor arrays | |
CN112162266B (zh) | 一种基于凸优化理论的共形阵二维波束优化方法 | |
JP4491081B2 (ja) | アンテナ | |
CN109541526B (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN112099017B (zh) | 一种基于过扫描的圆环阵低旁瓣波束优化方法 | |
CN107170441A (zh) | 圆环阵最优频率不变响应超指向性波束形成方法 | |
CN101149435A (zh) | 一种u型阵波束形成加权方法 | |
Wang et al. | Arbitrary-order superdirectivity of circular sensor arrays | |
Wang et al. | High-order superdirectivity of circular sensor arrays mounted on baffles | |
CN107238829B (zh) | 一种双层圆环阵稳健超指向性波束形成方法 | |
Ferguson | Sharpness applied to the adaptive beamforming of acoustic data from a towed array of unknown shape | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN116564265A (zh) | 基于参量次级声源的水下结构噪声控制方法 | |
JP4248294B2 (ja) | 不定項を用いたマイクロホンによるビームフォーミング | |
CN116679306A (zh) | 一种基于系数修正的扩展孔径声纳成像方法 | |
CN115085784A (zh) | 一种分级加权稀疏子阵天线阵列稳健波束合成方法 | |
CN112073873B (zh) | 一种无冗余阵元的一阶可调差分阵列的优化设计方法 | |
Yu et al. | A robust wavenumber-domain superdirective beamforming for endfire arrays | |
CN111830458B (zh) | 一种平行线阵单快拍二维测向方法 | |
Patel et al. | On the design of optimal linear microphone array geometries | |
CN109839611B (zh) | 一种适用平面阵的加权傅里叶积分方法 | |
CN113624330A (zh) | 一种水下目标辐射噪声测量组合体积阵及测量方法 | |
CN110489780B (zh) | 一种由指向性声传感器组成的端射直线阵波束形成方法 | |
CN111610489B (zh) | 一种基于阶数调整的任意阵超指向性波束优化方法 | |
Kuznetsov et al. | Equations for Calculating the Amplitude–Frequency and Phase–Frequency Responses of a Tripole-Type Vector–Scalar Receiver with a Time Delay of a Monopole Signal |
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 |