CN111610489A - 一种基于阶数调整的任意阵超指向性波束优化方法 - Google Patents
一种基于阶数调整的任意阵超指向性波束优化方法 Download PDFInfo
- Publication number
- CN111610489A CN111610489A CN202010463504.2A CN202010463504A CN111610489A CN 111610489 A CN111610489 A CN 111610489A CN 202010463504 A CN202010463504 A CN 202010463504A CN 111610489 A CN111610489 A CN 111610489A
- Authority
- CN
- China
- Prior art keywords
- directivity
- order
- error sensitivity
- super
- matrix
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic 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/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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radio Transmission System (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明涉及一种基于阶数调整的任意阵超指向性波束优化方法,利用阵元域最优权值向量可转化为不同阶模态域最优权值向量的性质,通过指定指向性因子或误差敏感度函数,求取所需要的阶数,获得优化后的超指向性波束。本发明利用了下述性质克服了没有闭式解或者需要知道误差统计特性或者适用范围有限的不足:阵元域最优权值向量可转化为不同阶模态域最优权值向量,通过调整阶数,可实现给定指向性或稳健性的超指向性波束。有益效果:计算方法更为简单,更少依赖先验信息,适用范围更广。
Description
技术领域
本发明属于波束形成方法,涉及一种基于阶数调整的任意阵超指向性波束优化方法,适用于任意阵列的低信噪比目标检测以及目标方位的高分辨率估计,属于声学、阵列信号处理和声呐技术等领域。
背景技术
超指向性波束形成方法通常需要在指向性和稳健性之间折中,结合优化方法,可以在给定稳健性的条件下获得最大的指向性或者在给定指向性的条件下获得最好的稳健性,相对于最优超指向性方法,操作上具有更大的灵活性,实际应用起来更为方便。文献1“Robust Supergain beamforming for circular array via second-order coneprogramming.Appl.Acoust.,2005,66(9):1018-1032”公开了一种利用二阶锥优化的超指向性波束形成方法,能够在给定稳健性的条件下获得较高指向性,但该方法没有闭式解,只能利用数值方法求得权值向量;文献2“Design of robust superdirective arrays witha tunable tradeoff between directivity and frequency-invariance.IEEETrans.Signal Process.,2011,59(5):2169-2181”公开了一种利用误差统计特性在指向性、稳健性和频率不变性之间折中的超指向性方法,具有闭式解,但该方法需要提前知道误差的统计特性,而这在实际中往往难以获得;文献3“Robust superdirective frequency-invariant beamforming for circular sensor arrays.IEEE Signal Process.Lett.,2017,24(8):1193-1197”公开了一种利用期望波束拟合的超指向性优化方法,可以获得具有频率不变响应的宽带超指向性波束,但该方法需要事先指定期望波束,且仅适用于圆环阵。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于阶数调整的任意阵超指向性波束优化方法,给定期望指向性因子或误差敏感度函数值,可求出对应的阶数,得到符合要求的权值向量以合成超指向性波束。避免现有技术没有闭式解或者需要知道误差统计特性或者适用范围有限的不足。
技术方案
一种基于阶数调整的任意阵超指向性波束优化方法,其特征在于步骤如下:
步骤1、给定期望指向性因子Dd或误差敏感度函数Td:
给定的期望指向性因子Dd应满足0<Dd<Dmax,其中Dmax为最大指向性因子,为:
所述p(θ,φ)为阵列流形向量,其第m号元素为:
所述矩阵ρnn为噪声协方差矩阵,其第i行j列元素为:
其中:dij为第i和j号阵元间的距离,J0为第0阶柱形Bessel函数;
所述矩阵C为下三角矩阵,其元素为:
给定的期望误差敏感度函数Td应满足0<Td<Tmax,其中Tmax为最大误差敏感度函数为:
步骤2、依据给定的期望指向性因子Dd或误差敏感度函数Td计算阶数N′:
1、给定期望指向性因子Dd时,阶数N′为:
所述:N为满足下式成立的最大非负整数:
2、给定期望误差敏感度函数Td时,阶数N′为:
所述:N为满足下式成立的最大非负整数:
所述η由下式计算:
其中:
步骤3、计算最优权值向量:
权值向量w为:
所述ω为:
步骤4、合成超指向性波束,计算指向性指数和误差敏感度指数:
1、阶数N′的超指向性波束合成为:B(θ,φ)=wΗp(θ,φ)
指向性指数为;DI=10log10D;
误差敏感度指数为SI=10log10T。
有益效果
本发明提出的一种基于阶数调整的任意阵超指向性波束优化方法,利用阵元域最优权值向量可转化为不同阶模态域最优权值向量的性质,通过指定指向性因子或误差敏感度函数,求取所需要的阶数,获得优化后的超指向性波束。本发明利用了下述性质克服了没有闭式解或者需要知道误差统计特性或者适用范围有限的不足:阵元域最优权值向量可转化为不同阶模态域最优权值向量,通过调整阶数,可实现给定指向性或稳健性的超指向性波束。
有益效果体现在:
1、本发明公开的方法在给定期望指向性因子或误差敏感度函数的条件下,给出了对应阶数的闭式求解公式,相对于文献1公开的数值计算方法更为简单。
2、本发明公开的方法不需要利用误差的统计特性等信息,比文献2公开的方法更少依赖先验信息。
3、本发明公开的方法适用于任意阵,比文献3公开的方法适用范围更广。
附图说明
图1是任意阵示意图。
图2是给定指向性因子条件下得到的阶数。
图3是给定指向性因子条件下得到的宽带波束图、指向性指数和误差敏感度指数。图3(a)是宽带波束图,图3(b)是指向性指数,图3(c)是误差敏感度指数。
图4是给定误差敏感度函数条件下得到的阶数。
图5是给定误差敏感度函数条件下得到的宽带波束图、指向性指数和误差敏感度指数。图5(a)是宽带波束图,图5(b)是指向性指数,图5(c)是误差敏感度指数。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
(1)参照图1。考虑一任意阵,阵元数M=8,阵元位置坐标如图1所示,平面波入射方向为(θ0,φ0)=(90°,180°),即x轴负方向,频率范围为[50Hz,1000Hz]。
给定期望指向性因子Dd=109/10,在所给频率范围内均满足0<Dd<Dmax,其中Dmax为最大指向性因子,由下式计算:
所述p(θ,φ)为阵列流形向量,其第m号元素为
所述矩阵ρnn为噪声协方差矩阵,其第i行j列元素为
其中dij为第i和j号阵元间的距离,J0为第0阶柱形Bessel函数。后续仿真中仅考虑三维各向同性均匀噪声场。
所述矩阵C为下三角矩阵,其元素为
给定期望误差敏感度函数Td=1015/10,在频率范围[50Hz,400Hz]内满足0<Td<Tmax,在频率范围(400Hz,1000Hz]内Td≥Tmax,其中Tmax为最大误差敏感度函数,由下式计算:
(2)参照图2和图4。依据给定的期望指向性因子Dd或误差敏感度函数Td计算阶数N′。
①给定期望指向性因子Dd=109/10时,阶数N′由下式计算:
所述N为满足下式成立的最大非负整数:
所述η由下式计算:
最终得到的阶数如图2所示,其值在4.3到5.8之间,且随频率升高而变大。
②给定期望误差敏感度函数Td=1015/10时,阶数N′同样为:
但所述N为满足下式成立的最大非负整数:
所述η由下式计算:
其中:
最终得到的阶数如图4所示,其值在0.2到7之间,在频率范围[50Hz,400Hz]内,其值随频率升高而变大,在频率范围(400Hz,1000Hz]内均等于7。
(3)计算最优权值向量。
权值向量w由下式计算:
所述ω为
(4)参照图3和图5。合成超指向性波束,计算指向性指数和误差敏感度指数。
①阶数N′的超指向性波束由下式合成:
B(θ,φ)=wΗp(θ,φ) (13)
②阶数N′的指向性因子D由下式计算:
指向性指数为DI=10log10D。
③阶数N′的误差敏感度函数T由下式计算:
误差敏感度指数为SI=10log10T。
给定期望指向性因子Dd=109/10时,得到的宽带波束图、指向性指数和误差敏感度指数如图3所示。由图3(a)可知,波束图具有一定的不对称性,这是由于阵列的不对称性导致的。图3(b)所示本发明方法的指向性指数在所示频率范围内均等于期望值9dB,且大于常规方法的指向性指数而小于MVDR方法的值,具有较好的超指向性。与之对应的误差敏感度指数如图3(c)所示,本发明方法的值同样大于常规方法且大部分小于MVDR方法的值,意味着其稳健性优于MVDR方法而劣于常规方法,这是指向性指数的提升带来的代价。
给定期望误差敏感度函数Td=1015/10时,得到的宽带波束图、指向性指数和误差敏感度指数如图5所示。由图5(a)可知,波束图同样具有一定的不对称性。图5(b)所示本发明方法的指向性指数在频率范围[50Hz,400Hz]内随频率升高而变大,且大于常规方法的指向性指数而小于MVDR方法的值,具有较好的超指向性;而在频率范围(400Hz,1000Hz]内,其值等于MVDR方法的值。对于图5(c)所示的误差敏感度指数,本发明方法的值在频率范围[50Hz,400Hz]内均等于期望值15dB,且大于常规方法而小于MVDR方法的值,意味着其稳健性优于MVDR方法而劣于常规方法;而在频率范围(400Hz,1000Hz]内,其值等于MVDR方法的值,与指向性指数的规律一致。
由上述仿真结果可知,在给定期望指向性因子或误差敏感度函数后,可直接计算超指向性的对应阶数,且阶数为非负实数,基于得到的阶数,即可得到优化后的波束图、指向性指数和误差敏感度指数。
Claims (1)
1.一种基于阶数调整的任意阵超指向性波束优化方法,其特征在于步骤如下:
步骤1、给定期望指向性因子Dd或误差敏感度函数Td:
给定的期望指向性因子Dd应满足0<Dd<Dmax,其中Dmax为最大指向性因子,为:
所述p(θ,φ)为阵列流形向量,其第m号元素为:
所述矩阵ρnn为噪声协方差矩阵,其第i行j列元素为:
其中:dij为第i和j号阵元间的距离,J0为第0阶柱形Bessel函数;
所述矩阵C为下三角矩阵,其元素为:
给定的期望误差敏感度函数Td应满足0<Td<Tmax,其中Tmax为最大误差敏感度函数为:
步骤2、依据给定的期望指向性因子Dd或误差敏感度函数Td计算阶数N′:
1、给定期望指向性因子Dd时,阶数N′为:
所述:N为满足下式成立的最大非负整数:
2、给定期望误差敏感度函数Td时,阶数N′为:
所述:N为满足下式成立的最大非负整数:
所述η由下式计算:
其中:
步骤3、计算最优权值向量:
权值向量w为:
所述ω为:
步骤4、合成超指向性波束,计算指向性指数和误差敏感度指数:
1、阶数N′的超指向性波束合成为:B(θ,φ)=wΗp(θ,φ)
指向性指数为;DI=10log10D;
误差敏感度指数为SI=10log10T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010463504.2A CN111610489B (zh) | 2020-05-27 | 2020-05-27 | 一种基于阶数调整的任意阵超指向性波束优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010463504.2A CN111610489B (zh) | 2020-05-27 | 2020-05-27 | 一种基于阶数调整的任意阵超指向性波束优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111610489A true CN111610489A (zh) | 2020-09-01 |
CN111610489B CN111610489B (zh) | 2022-08-23 |
Family
ID=72200544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010463504.2A Active CN111610489B (zh) | 2020-05-27 | 2020-05-27 | 一种基于阶数调整的任意阵超指向性波束优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111610489B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101964934A (zh) * | 2010-06-08 | 2011-02-02 | 浙江大学 | 二元麦克风微阵列语音波束形成方法 |
CN103902831A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种基于模态分解与综合的超指向性波束形成方法 |
CN103902830A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 |
CN105759239A (zh) * | 2016-03-09 | 2016-07-13 | 临境声学科技江苏有限公司 | 一种降阶的频不变的鲁棒的超指向波束形成算法 |
CN107238829A (zh) * | 2017-06-22 | 2017-10-10 | 西北工业大学 | 一种双层圆环阵稳健超指向性波束形成方法 |
CN109600152A (zh) * | 2018-12-17 | 2019-04-09 | 西北工业大学 | 一种基于子空间基变换的自适应波束形成方法 |
CN109959917A (zh) * | 2019-03-08 | 2019-07-02 | 南京航空航天大学 | 一种宽带非频变多波束成像声纳的阵列稀疏方法 |
-
2020
- 2020-05-27 CN CN202010463504.2A patent/CN111610489B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101964934A (zh) * | 2010-06-08 | 2011-02-02 | 浙江大学 | 二元麦克风微阵列语音波束形成方法 |
CN103902831A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种基于模态分解与综合的超指向性波束形成方法 |
CN103902830A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 |
CN105759239A (zh) * | 2016-03-09 | 2016-07-13 | 临境声学科技江苏有限公司 | 一种降阶的频不变的鲁棒的超指向波束形成算法 |
CN107238829A (zh) * | 2017-06-22 | 2017-10-10 | 西北工业大学 | 一种双层圆环阵稳健超指向性波束形成方法 |
CN109600152A (zh) * | 2018-12-17 | 2019-04-09 | 西北工业大学 | 一种基于子空间基变换的自适应波束形成方法 |
CN109959917A (zh) * | 2019-03-08 | 2019-07-02 | 南京航空航天大学 | 一种宽带非频变多波束成像声纳的阵列稀疏方法 |
Non-Patent Citations (4)
Title |
---|
WANG YONG ETC.,: "High-Order Superdirectivity of Circular Sensor", 《ACTA ACUSTICA UNITED WITH ACUSTICA》 * |
WANG YONG ETC.,: "Robust Superdirective Frequency-Invariant", 《IEEE SIGNAL PROCESSING LETTERS》 * |
汪勇 等: "圆环阵超指向性方法仿真研究", 《声学技术》 * |
郭俊媛 等: "基于超指向性多极子矢量阵的水下低频声源方位", 《物理学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111610489B (zh) | 2022-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hung | Matrix-construction calibration method for antenna arrays | |
CN111693971B (zh) | 一种用于弱目标检测的宽波束干扰抑制方法 | |
CN103902830A (zh) | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 | |
Xia et al. | Noise reduction method for acoustic sensor arrays in underwater noise | |
CN107170441B (zh) | 圆环阵最优频率不变响应超指向性波束形成方法 | |
CN104768099B (zh) | 用于圆环形阵的模态波束形成器及频域宽带实现方法 | |
CN109541526A (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN108595758A (zh) | 一种任意形式传感器阵列最优宽带波束图综合方法 | |
CN111610489B (zh) | 一种基于阶数调整的任意阵超指向性波束优化方法 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN107241131B (zh) | 一种利用信号非圆特性的波束形成方法 | |
Yang et al. | Robust adaptive beamformer using interpolation technique for conformal antenna array | |
Sun et al. | Robust spherical microphone array beamforming with multi-beam-multi-null steering, and sidelobe control | |
Luo et al. | Constrained maximum directivity beamformers based on uniform linear acoustic vector sensor arrays | |
CN113160843B (zh) | 基于质点振速传感器微阵列的干扰语音抑制方法及装置 | |
CN113820654B (zh) | 基于波束域降维的s波段雷达目标低仰角doa估计方法 | |
CN113624330B (zh) | 一种水下目标辐射噪声测量组合体积阵及测量方法 | |
Zhang et al. | Frequency-invariant beamformer design via ADPM approach | |
CN114384502A (zh) | 一种基于稀疏阵列的相干增益处理方法 | |
CN112433218B (zh) | 一种舰用共形阵虚拟障板的实现方法 | |
CN110489780B (zh) | 一种由指向性声传感器组成的端射直线阵波束形成方法 | |
CN112649787B (zh) | 一种基于低频圆环阵的目标方位估计方法 | |
CN111474533B (zh) | 一种圆环阵任意阶超指向性波束形成方法 | |
Wang et al. | Robust constant-beamwidth beamforming based on focusing approach for acoustic imaging | |
Sun et al. | Robust adaptive beamforming method for active sonar in single snapshot |
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 |