CN107170441B - 圆环阵最优频率不变响应超指向性波束形成方法 - Google Patents
圆环阵最优频率不变响应超指向性波束形成方法 Download PDFInfo
- Publication number
- CN107170441B CN107170441B CN201710479869.2A CN201710479869A CN107170441B CN 107170441 B CN107170441 B CN 107170441B CN 201710479869 A CN201710479869 A CN 201710479869A CN 107170441 B CN107170441 B CN 107170441B
- Authority
- CN
- China
- Prior art keywords
- weight vector
- array
- optimal
- frequency
- circular array
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/18—Methods or devices for transmitting, conducting or directing sound
- G10K11/26—Sound-focusing or directing, e.g. scanning
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明涉及一种圆环阵最优频率不变响应超指向性波束形成方法,针对均匀圆环形阵列提出一种最优频率不变响应超指向性波束形成方法。该方法首先给出了一般化的期望波束表达式,并将不同频率处的权值向量表达成期望波束权值向量的函数,然后构建多约束优化问题,获得满足约束条件的最优期望波束,最后通过函数关系获得其它频率处的权值向量,合成得到宽带频率不变响应超指向性波束。可以灵活地控制波束形状,可以获得更宽频率范围内的频率不变响应超指向性波束,可以获得满足约束条件的最优期望波束。
Description
技术领域
本发明属于声学阵列信号处理和声呐技术等领域,是一种波束形成方法,涉及一种圆环阵最优频率不变响应超指向性波束形成方法,特别涉及一种频率不变响应超指向性波束形成方法,适用于圆环形阵列的宽带目标检测以及宽带目标方位估计。
背景技术
在声呐、雷达、通信以及语音工程等领域,常用阵列接收宽带信号并进行处理,以实现目标检测、定位、分辨等目的。频率不变响应超指向性波束形成方法既能保证宽带信号无畸变,也能获得更好的空间指向性和信噪比增益,可显著提高阵列处理的相关性能,因此备受关注。
已有的频率不变响应超指向性波束形成方法主要有:文献1“Design of robustsuperdirective arrays with a tunable tradeoff between directivity andfrequency-invariance.IEEE Trans.Signal Process.,2011,59(5):2169-2181”公开的加权优化方法,该方法理论上可应用于任意阵型,但其形式较复杂,未能获得简易的闭式解,也难以灵活地控制波束形状。文献2“Design of circular differential microphonearrays.Springer topics in signal processing,ed.J.Benesty andW.Kellermann.vol.12.Berlin:Springer-Verlag,2015.”公开的差分方法,该方法有两种具体形式,一种是在满足相关无畸变和零点约束的条件下获得权值向量,另一种是基于有限阶Jacobi-Anger级数展开拟合期望波束。然而,这两种形式合成的波束仍然不够精确,其频率不变性能可以进一步提高。文献3“一种圆环形阵列恒定束宽波束形成器设计方法.中国:CN103903609 B,2016”公开的期望波束拟合方法,该方法通过最小化合成波束和期望波束之间的均方误差,得到了圆环阵在不同频率上权值向量的精确闭式解,可方便地合成宽带频率不变响应波束。但该方法的期望波束是由相同阵元数的圆环阵合成得到,且只能事先给定,不能保证期望波束是满足需求的最佳选择。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种圆环阵最优频率不变响应超指向性波束形成方法,克服了现有技术频率不变性能有限和难以给出最佳期望波束的不足:不同频率处的权值向量是一般化期望波束权值向量的函数。
技术方案
一种圆环阵最优频率不变响应超指向性波束形成方法,其特征在于适用于均匀圆环形阵列,步骤如下:
步骤1:建立均匀圆环形阵列的一般化期望波束的表达式:
步骤2、计算最优期望波束权值向量:
构建多约束优化问题:
所述ΩSL是波束旁瓣区域,Bd表示期望波束响应向量, 其中,r为圆环形阵列的半径,k=2π/λ,λ表示入射平面波的波长,k1r1是感兴趣的kr范围内的最小值,k2r2是感兴趣的kr范围内的最大值,L为k1r1到k2r2之间离散化点的个数;
当阵元数M为偶数且N≤M/2的情况时最优期望波束权值向量:
Σ2=diag{1-M|b0|2/γ0,1-M|b1|2/γ1,...,1-M|bM/2|2/γM/2,1-M|bM/2-1|2/γM/2-1,...,1-M|b1|2/γ1}
当阵元数M为奇数且N≤(M-1)/2的情况时最优期望波束权值向量:
Σ2=diag{1-M|b0|2/γ0,1-M|b1|2/γ1,...,1-M|b(M-1)/2|2/γ(M-1)/2,1-M|b(M-1)/2-1|2/γ(M-1)/2-1,...,1-M|b1|2/γ1}
其中:
步骤3、计算在所需频率范围内合成波束的权值向量:
考虑当阵元数M为偶数且N≤M/2的情况,期望波束权值向量与合成波束权值向量的函数关系为:
考虑当阵元数M为奇数且N≤(M-1)/2的情况,期望波束权值向量与合成波束权值向量的函数关系为:
步骤4、计算阵元域权值向量:w=Vω
步骤5、合成最终波束:
其中:P(φ)=[p0(φ),p1(φ),...,pM-1(φ)]T为阵列流形向量,ps(φ)为第s个阵元接收到的从方向φ入射的单位幅度平面波信号,表达式为:φs=sβ,(·)*表示求共轭,上标H表示共轭转置。
有益效果
本发明提出的一种圆环阵最优频率不变响应超指向性波束形成方法,针对均匀圆环形阵列提出一种最优频率不变响应超指向性波束形成方法。该方法首先给出了一般化的期望波束表达式,并将不同频率处的权值向量表达成期望波束权值向量的函数,然后构建多约束优化问题,获得满足约束条件的最优期望波束,最后通过函数关系获得其它频率处的权值向量,合成得到宽带频率不变响应超指向性波束。
有益效果体现在:
1.本发明公开的方法将不同频率处的权值向量闭式地表达成一般化期望波束权值向量的函数,也可以灵活地控制波束形状,比文献1公开的方法更优越。
2.本发明公开的方法可以获得更宽频率范围内的频率不变响应超指向性波束,相比于文献2公开的方法,其频率不变性能更好。
本发明公开的方法可以获得满足约束条件的最优期望波束,比文献3公开的方法更优越。
附图说明
图1是圆环形阵列示意图。
图2是最优期望波束图。
图3是不同方法的指向性指数。
图4是不同方法的白噪声增益
图5是合成的频率不变响应超指向性波束图。图5(a)是波束图随频率变化的三维显示,图5(b)是波束图随频率变化的叠加显示。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明实施例可在最小均方误差准则下给出均匀圆环阵不同频率处的权值向量与期望波束权值向量的函数,并获得了满足约束条件的最优期望波束,进而合成得到宽带频率不变响应超指向性波束。其过程为:
1.给出一般化期望波束的表达式。由于圆环形阵列更多关注的是水平面内波束形成的性能,所以仅考虑水平范围内的二维波束。一般化的期望波束可表示为如下形式:
2.参照图1和2。该圆环形阵列包含M个均匀分布的阵元。构建如下多约束优化问题:
所述ΩSL是波束旁瓣区域,Bd表示期望波束响应向量,由式(1)所示元素组成,其中,r为圆环形阵列的半径,k=2π/λ,λ表示入射平面波的波长,k1r1是感兴趣的kr范围内的最小值,k2r2是感兴趣的kr范围内的最大值,L为k1r1到k2r2之间离散化点的个数。
考虑当阵元数M为偶数且N≤M/2的情况,有:
Σ2=diag{1-M|b0|2/γ0,1-M|b1|2/γ1,...,1-M|bM/2|2/γM/2,1-M|bM/2-1|2/γM/2-1,...,1-M|b1|2/γ1}。所述ρs=J0[2kr sin(sβ/2)],β=2π/M,J0(·)是第0阶柱形Bessel函数,diag{·}表示对角阵。
考虑当阵元数M为奇数且N≤(M-1)/2的情况,有:
Σ2=diag{1-M|b0|2/γ0,1-M|b1|2/γ1,...,1-M|b(M-1)/2|2/γ(M-1)/2,1-M|b(M-1)/2-1|2/γ(M-1)/2-1,...,1-M|b1|2/γ1}。
其中:
对于式(2)的使用,有以下说明:μq(q∈{1,2,3,4})中可任选一个作为优化目标,其余三个则为给定的约束值。第一个约束用于控制稳健性。第二个约束用于控制合成的最小均方误差最小。第三个约束用于控制波束旁瓣响应,q2=1,2或者∞,一般取q2=∞,表示约束波束旁瓣,对应的旁瓣级为20lgμ2dB。如果取q2=2,则表示约束均方旁瓣。第四个约束用于控制波束指向性。若要使某一项约束失效,可直接设μq=∞。上式所示的多约束优化问题均可由二阶锥规划进行求解。
以声透明均匀圆环形阵列为例进行仿真,设定相关参数为:M=16,N=3,ka∈[0.1,8],ΩSL={(θ,φ)|θ=90°,φ∈[0°,φ0-Δ]∪[φ0+Δ,360°]},φ0=180°,Δ=80°,q=4,μ1=105,μ2=10-4,μ3=0.1,q2=2,L=80。由式计算得到的最优期望波束权值向量列于表1,最优期望波束如图2所示,其旁瓣级小于-20dB,满足约束条件。
表1 最优期望波束权值向量
a<sub>-3</sub> | a<sub>-2</sub> | a<sub>-1</sub> | a<sub>0</sub> | a<sub>1</sub> | a<sub>2</sub> | a<sub>3</sub> |
0.0184 | 0.1525 | 0.2066 | 0.2451 | 0.2066 | 0.1525 | 0.0184 |
3.计算在所需频率范围内合成波束的权值向量。考虑当阵元数M为偶数且N≤M/2的情况,期望波束权值向量与合成波束权值向量的函数关系由下式给出:
考虑当阵元数M为奇数且N≤(M-1)/2的情况,期望波束权值向量与合成波束权值向量的函数关系由下式给出:
将表1所列权值向量代入式(4)或(5),即可得到所需要的合成波束权值向量ω=[ω0,ω1,...,ωM-1]T。
4.计算阵元域权值向量。将由式(4)或(5)得到的权值向量通过下式转化为阵元域权值向量:
w=Vω (6)
5.合成最终波束。将式(6)得到的权值向量代入式(7),即可得到合成波束:
所述P(φ)=[p0(φ),p1(φ),...,pM-1(φ)]T为阵列流形向量,ps(φ)为第s个阵元接收到的从方向φ入射的单位幅度平面波信号,表达式为:φs=sβ,(·)*表示求共轭,上标H表示共轭转置。
合成得到的频率不变响应超指向性波束的指向性指数和白噪声增益分别如图3和图4所示,波束图如图5所示。当kr<2.6时,合成波束的指向性指数大于常规方法的值,具有较好的超指向性;对应于最小频率kr=0.1的白噪声增益约为-50dB,满足约束条件。在所示频率范围内,合成得到的宽带波束具有很好的频率不变特性,与期望波束的误差小于设定上限10-4。
Claims (1)
1.一种圆环阵最优频率不变响应超指向性波束形成方法,其特征在于适用于均匀圆环形阵列,步骤如下:
步骤1:建立均匀圆环形阵列的一般化期望波束的表达式:
步骤2、计算最优期望波束权值向量:
构建多约束优化问题:
所述ΩSL是波束旁瓣区域,Bd表示期望波束响应向量, 其中,r为圆环形阵列的半径,k=2π/λ,λ表示入射平面波的波长,k1r1是感兴趣的kr范围内的最小值,k2r2是感兴趣的kr范围内的最大值,L为k1r1到k2r2之间离散化点的个数;
当阵元数M为偶数且N≤M/2的情况时最优期望波束权值向量:
当阵元数M为奇数且N≤(M-1)/2的情况时最优期望波束权值向量:
其中:
步骤3、计算在所需频率范围内合成波束的权值向量:
考虑当阵元数M为偶数且N≤M/2的情况,期望波束权值向量与合成波束权值向量的函数关系为:
考虑当阵元数M为奇数且N≤(M-1)/2的情况,期望波束权值向量与合成波束权值向量的函数关系为:
步骤4、计算阵元域权值向量:w=Vω
步骤5、合成最终波束:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710479869.2A CN107170441B (zh) | 2017-06-22 | 2017-06-22 | 圆环阵最优频率不变响应超指向性波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710479869.2A CN107170441B (zh) | 2017-06-22 | 2017-06-22 | 圆环阵最优频率不变响应超指向性波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107170441A CN107170441A (zh) | 2017-09-15 |
CN107170441B true CN107170441B (zh) | 2020-10-27 |
Family
ID=59819096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710479869.2A Active CN107170441B (zh) | 2017-06-22 | 2017-06-22 | 圆环阵最优频率不变响应超指向性波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107170441B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549052B (zh) * | 2018-03-20 | 2021-04-13 | 南京航空航天大学 | 一种时频-空域联合加权的圆谐域伪声强声源定位方法 |
CN108595758B (zh) * | 2018-03-22 | 2021-11-09 | 西北工业大学 | 一种任意形式传感器阵列最优宽带波束图综合方法 |
CN109493844A (zh) * | 2018-10-17 | 2019-03-19 | 南京信息工程大学 | 基于fir滤波器的恒定束宽波束形成方法 |
CN111273302B (zh) * | 2020-03-02 | 2022-08-23 | 西北工业大学 | 一种浅海匀速运动目标初始状态估计方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102998660A (zh) * | 2012-11-26 | 2013-03-27 | 哈尔滨工程大学 | 一种近场范围内鲁棒多波束形成方法 |
CN103069824A (zh) * | 2010-08-20 | 2013-04-24 | 索尼公司 | 用于交互式电视的服务器负荷平衡 |
CN103438987A (zh) * | 2013-08-30 | 2013-12-11 | 西北工业大学 | 基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法 |
CN103630887A (zh) * | 2013-11-11 | 2014-03-12 | 南京航空航天大学 | 一种多波束形成方法及使用该方法的多波束声纳 |
CN103903609A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列恒定束宽波束形成器设计方法 |
CN103902831A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种基于模态分解与综合的超指向性波束形成方法 |
CN103902830A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 |
CN103902832A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列实数加权超指向性波束形成方法 |
CN104768099A (zh) * | 2014-01-02 | 2015-07-08 | 中国科学院声学研究所 | 用于圆环形阵的模态波束形成器及频域宽带实现方法 |
CN105929384A (zh) * | 2016-07-06 | 2016-09-07 | 西北工业大学 | 基于十字阵发射圆环阵接收的联合低旁瓣多波束形成方法 |
KR101718282B1 (ko) * | 2016-04-11 | 2017-03-20 | 한국과학기술원 | 가시선 채널 환경에서 부배열 안테나를 이용한 균일 원형 배열 안테나 시스템 및 그의 빔포밍 방법 |
-
2017
- 2017-06-22 CN CN201710479869.2A patent/CN107170441B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103069824A (zh) * | 2010-08-20 | 2013-04-24 | 索尼公司 | 用于交互式电视的服务器负荷平衡 |
CN102998660A (zh) * | 2012-11-26 | 2013-03-27 | 哈尔滨工程大学 | 一种近场范围内鲁棒多波束形成方法 |
CN103438987A (zh) * | 2013-08-30 | 2013-12-11 | 西北工业大学 | 基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法 |
CN103630887A (zh) * | 2013-11-11 | 2014-03-12 | 南京航空航天大学 | 一种多波束形成方法及使用该方法的多波束声纳 |
CN104768099A (zh) * | 2014-01-02 | 2015-07-08 | 中国科学院声学研究所 | 用于圆环形阵的模态波束形成器及频域宽带实现方法 |
CN103903609A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列恒定束宽波束形成器设计方法 |
CN103902831A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种基于模态分解与综合的超指向性波束形成方法 |
CN103902830A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 |
CN103902832A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列实数加权超指向性波束形成方法 |
KR101718282B1 (ko) * | 2016-04-11 | 2017-03-20 | 한국과학기술원 | 가시선 채널 환경에서 부배열 안테나를 이용한 균일 원형 배열 안테나 시스템 및 그의 빔포밍 방법 |
CN105929384A (zh) * | 2016-07-06 | 2016-09-07 | 西北工业大学 | 基于十字阵发射圆环阵接收的联合低旁瓣多波束形成方法 |
Non-Patent Citations (7)
Title |
---|
"A superdirective beamforming method for linear sensor arrays";Yixin Yang 等;《OCEANS 2014-PAITEI》;IEEE;20141124;共6页 * |
"Adaptive Beamforming Using Frequency Invariant Uniform Concentric Circular Arrays";H. H. Chen 等;《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS》;20071117;第54卷(第9期);第1938-1949页 * |
"Robust Superdirective Frequency-Invariant Beamforming for Circular Sensor Arrays";Yong Wang 等;《IEEE SIGNAL PROCESSING LETTERS》;IEEE;20170605;第24卷(第8期);第1193-1197页 * |
"Uniform Concentric Circular Arrays With Frequency-Invariant Characteristics—Theory,Design, Adaptive Beamforming and DOA Estimation";S. C. Chan 等;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;20061219;第55卷(第1期);第165-177页 * |
"传感器阵列超指向性原理及方法研究";汪勇;《CNKI中国博士学位论文全文数据库信息科技辑》;20170115;I140-19 * |
"刚性球上圆阵波束设计及性能仿真";张成 等;《火力与指挥控制》;20101130;第35卷(第11期);第125-137页 * |
"小尺度传感器阵列超指向性研究及实现";蒋伟;《CNKI中国优秀硕士学位论文全文数据库信息科技辑》;20070115;I140-29 * |
Also Published As
Publication number | Publication date |
---|---|
CN107170441A (zh) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107170441B (zh) | 圆环阵最优频率不变响应超指向性波束形成方法 | |
US9237391B2 (en) | Low noise differential microphone arrays | |
Rasekh et al. | Design of an adaptive wideband beamforming algorithm for conformal arrays | |
CN109581388B (zh) | 一种实时三维成像声纳的近场宽视角波束形成方法 | |
CN103902830B (zh) | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 | |
CN109541526B (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
Xia et al. | Noise reduction method for acoustic sensor arrays in underwater noise | |
CN108595758B (zh) | 一种任意形式传感器阵列最优宽带波束图综合方法 | |
CN104768099B (zh) | 用于圆环形阵的模态波束形成器及频域宽带实现方法 | |
CN111239747B (zh) | 一种基于解卷积的声纳高分辨低旁瓣二维成像方法 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN116679306A (zh) | 一种基于系数修正的扩展孔径声纳成像方法 | |
Chen et al. | SRV constraint based FIB design for wideband linear array | |
Liang et al. | Direction of arrival estimation under near-field interference using matrix filter | |
Zhang et al. | Frequency-invariant beamformer design via ADPM approach | |
CN110489780B (zh) | 一种由指向性声传感器组成的端射直线阵波束形成方法 | |
Xie et al. | Broadband frequency invariant beamformer | |
Zhang et al. | An approach for improving the performance of generalized inverse beamforming | |
CN111610489B (zh) | 一种基于阶数调整的任意阵超指向性波束优化方法 | |
CN112433218B (zh) | 一种舰用共形阵虚拟障板的实现方法 | |
CN113030983B (zh) | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 | |
Zhao et al. | Research on robust time domain wideband null beamforming interference suppression | |
Wang et al. | Optimal design of frequency-invariant beamformers for circular arrays | |
CN118157731A (zh) | 一种主瓣干扰下基于波束域的自适应接收波束保形方法 | |
CN109917358B (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 |