CN109541526B - 一种利用矩阵变换的圆环阵方位估计方法 - Google Patents
一种利用矩阵变换的圆环阵方位估计方法 Download PDFInfo
- Publication number
- CN109541526B CN109541526B CN201811637625.3A CN201811637625A CN109541526B CN 109541526 B CN109541526 B CN 109541526B CN 201811637625 A CN201811637625 A CN 201811637625A CN 109541526 B CN109541526 B CN 109541526B
- Authority
- CN
- China
- Prior art keywords
- array
- azimuth
- covariance matrix
- vector
- array manifold
- 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
- 239000011159 matrix material Substances 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000009466 transformation Effects 0.000 title claims abstract description 10
- 239000013598 vector Substances 0.000 claims abstract description 41
- 238000001228 spectrum Methods 0.000 claims abstract description 23
- 230000017105 transposition Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000001629 suppression 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
- 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
- G01S3/802—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种利用矩阵变换的圆环阵方位估计方法,首先计算圆环阵的原始阵列流形向量,接着计算得到空间均匀噪声场中噪声协方差矩阵的特征向量,然后利用特征向量和原始阵列流形向量设计出一个新的阵列流形向量,同时利用特征向量对原始数据协方差矩阵进行变换,得到降维的协方差矩阵,最后用设计出的新的阵列流形向量和降维的数据协方差矩阵进行目标方位估计。稳健性高,并且灵活性很高,可以根据不同阶数的结果挑选稳健性较好、估计较优的方位谱。实现了矩阵的降维,计算量大大减小。设计出的阵列流形向量在理论上是一个精确的解析公式。
Description
技术领域
本发明属于声学阵列信号处理、声纳信号处理等领域,涉及一种利用矩阵变换的圆环阵方位估计方法,适用于目标方位估计等领域。
背景技术
由无指向性声压传感器组成的圆环阵列在空气中传声器、水下声纳等***中具有广泛的应用,常被用来进行目标定位、噪声抑制等,AI智能音箱上的多个传声器大多也都均匀分布在圆弧上。早期由无指向性传感器组成的阵列实现目标方位估计的方法主要有:(1)文献1“Frequency-Wavenumber Spectrum Analysis.Proceedings of the IEEE,1969,57(8):1408-1418.”公开的方法,适用于任意形状的阵列,但是其算法稳健性较差,常常无法正确估计方位;(2)文献2“一种改进的稀疏近似最小方差DOA估计算法研究.声学学报,2016(4):465-476.”利用一个折衷参数进行最大似然估计值和稀疏性能的折衷处理,在迭代过程中改变稀疏近似最小方差算法,但是此方法计算量过大,并不实用;(3)文献3“Optimaldesign of modal beamformers for circular arrays.Journal of the AcousticalSociety of America,2015,138(4):2140-2151.”针对圆环形采用模态域波束形成方法,但是对模态阶数进行了截断,在理论上并不是一个闭式解析解,并且过程较为复杂。
现有方法稳健性差,当误差存在时,性能衰退严重的的缺陷。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种利用矩阵变换的圆环阵方位估计方法。
技术方案
一种利用矩阵变换的圆环阵方位估计方法,其特征在于步骤如下:
步骤1:计算得到圆环阵的原始阵列流形向量:
P(θ,φ)=[p0(θ,φ),...,pm(θ,φ),...,pM-1(θ,φ)]T (1)
步骤2:设计新的第m阶阵列流形向量:
所述E=[E0,E1,...,EM-1]T (2)
步骤3:计算第m阶新的噪声协方差矩阵:
步骤4:利用新的阵列流形向量和噪声协方差矩阵,计算第m阶方位谱以进行目标方位估计:
有益效果
本发明提出的一种利用矩阵变换的圆环阵方位估计方法,首先计算圆环阵的原始阵列流形向量,接着计算得到空间均匀噪声场中噪声协方差矩阵的特征向量,然后利用特征向量和原始阵列流形向量设计出一个新的阵列流形向量,同时利用特征向量对原始数据协方差矩阵进行变换,得到降维的协方差矩阵,最后用设计出的新的阵列流形向量和降维的数据协方差矩阵进行目标方位估计。
本发明结合空间均匀噪声场中的噪声协方差矩阵的特征向量与原始阵列流形向量,设计出一组新的阵列流形向量,并定义了不同阶数的阵列流形向量的表达式,然后计算出新的噪声协方差矩阵,再利用新得到的阵列流形向量和新的噪声协方差矩阵进行目标方位估计,得到方位谱。
有益效果体现在:
本发明提出一种适用于圆环阵目标方位估计的设计方法,有益效果体现在:
1.本发明公开的方法设计出的圆环阵目标方位估计的方法的稳健性高于文献1公开的方法,并且灵活性很高,可以根据不同阶数的结果挑选稳健性较好、估计较优的方位谱。
2.本发明公开的方法将空间均匀噪声场中噪声协方差矩阵的特征向量用于目标方位估计,实现了矩阵的降维,计算量大大减小,而文献2的方法迭代次数多,计算量大,不适合实际应用。
3.本发明公开的方法更加简单灵活且设计出的阵列流形向量在理论上是一个精确的解析公式,文献3中方法模态波束形成器计算复杂,并且不是精确的解析解。
附图说明
图1是均匀分布有M个阵元的圆环阵的坐标示意图。仿真中假设水下声速c=1500m/s,M=12,半径r=0.5m。
图2给出了在信噪比为10dB(噪声为高斯白噪声),频率分别为f=200Hz,1000Hz,2000Hz和4000Hz时的不同阶数m(m=0,1,…,6)下的方位谱。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明涉及一种适用于圆环阵目标方位估计的设计方法。该方法结合空间均匀噪声场中的噪声协方差矩阵的特征向量与原始阵列流形向量,设计出一组新的阵列流形向量,并定义了不同阶数的阵列流形向量的表达式,然后通过矩阵变换计算出新的噪声协方差矩阵,再利用新得到的阵列流形向量和新的噪声协方差矩阵进行目标方位估计,得到方位谱。不同阶数的方位谱的稳健性不同,阶数越低,稳健性越高,在实际应用中根据需要可灵活选择某一阶数的方位谱。其过程为:
1.计算得到圆环阵的原始阵列流形向量,由下式得到:
P(θ,φ)=[p0(θ,φ),...,pm(θ,φ),...,pM-1(θ,φ)]T (7)
2.设计新的阵列流形向量,见下式:
E=[E0,E1,...,EM-1]T (8)
3.计算第m阶新的噪声协方差矩阵,由下式的矩阵变换得到:
4.利用新的阵列流形向量和噪声协方差矩阵进行目标方位估计,第m阶方位谱的由下式计算:
参照图2。图2给出了在信噪比为10dB(噪声为高斯白噪声),信号频率分别为f=200Hz,1000Hz,2000Hz和4000Hz时的不同阶数m(m=0,1,…,6)下的方位谱。图2(a)中信号频率为f=200Hz,除了m=0时的方位谱是一条直线外,其他阶数m的方位谱几乎重合;图2(b)中信号频率为f=1000Hz,m=1阶的方位谱已经独立出来;图2(c)中信号频率为f=2000Hz,m=1,2,3,4阶的方位谱都已经独立出来,m=5,6的方位谱仍然重合;图2(d)中信号频率为f=4000Hz,所有阶的方位谱都已经独立出来。这说明在没有误差的情况下,随着频率的升高,高阶方位谱的估计效果更好。与之相比,文献(1)仅能给出一个方位估计的结果;文献(2)计算过程过于复杂,不实用;文献(3)的截取的各阶模态所计算的阵列流形向量仍然是近似的结果,在理论上并不精确,也不是一个高分辨方位估计算法。
参照图3。图3是在信噪比为10dB(噪声为高斯白噪声),信号频率为f=1000Hz,引入接收信号的随机幅度误差(σg)和相位误差(σψ)的方差分别为10-3,10-2和10-1(分别对应图3(a),(b),(c)和(d))时,不同阶数m(m=0,1,…,6)下的方位谱,为避免随机误差的影响,这是200次随机实验进行平均后的结果。参照图3(a),当随机误差较小时,这对方位谱估计的结果影响不大;随着随机误差的增大,参照图3(b),(c)和(d),方位估计的效果越来越差,但可以看到,较低的阶数估计效果要优于高阶,并且阶数为m=2时的结果为最优。与之相比,文献(1)的稳健性较差,不能选择一个恰当且稳健的方位谱结果;文献(2)并没有给出有误差时方位估计方法稳健性的讨论;文献(3)只给出了1%失配误差存在下的波束图结果,并没有讨论误差存在下的方位谱的情况。
Claims (1)
1.一种利用矩阵变换的圆环阵方位估计方法,其特征在于步骤如下:
步骤1:计算得到圆环阵的原始阵列流形向量:
P(θ,φ)=[p0(θ,φ),...,pm(θ,φ),...,pM-1(θ,φ)]T
其中,M是阵元个数,上标“T”表示转置,θ为俯仰角,φ为水平方位角,pm(θ,φ)=exp(-j2πfrsinθcos(φ-φm)/c),φm是第m号阵元对应的方位角,即r为圆环阵半径,f为信号频率,c为声速;
步骤2:设计新的第m阶阵列流形向量:
所述E=[E0,E1,...,EM-1]T
步骤3:计算第m阶新的噪声协方差矩阵:
步骤4:利用新的阵列流形向量和噪声协方差矩阵,计算第m阶方位谱以进行目标方位估计:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811637625.3A CN109541526B (zh) | 2018-12-29 | 2018-12-29 | 一种利用矩阵变换的圆环阵方位估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811637625.3A CN109541526B (zh) | 2018-12-29 | 2018-12-29 | 一种利用矩阵变换的圆环阵方位估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109541526A CN109541526A (zh) | 2019-03-29 |
CN109541526B true CN109541526B (zh) | 2022-09-02 |
Family
ID=65831229
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811637625.3A Active CN109541526B (zh) | 2018-12-29 | 2018-12-29 | 一种利用矩阵变换的圆环阵方位估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109541526B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111562584B (zh) * | 2020-04-30 | 2022-08-16 | 北京中安智能信息科技有限公司 | 被动声纳方位历程图处理方法和装置及设备 |
CN112649787B (zh) * | 2020-11-30 | 2022-04-12 | 海鹰企业集团有限责任公司 | 一种基于低频圆环阵的目标方位估计方法 |
CN113608192B (zh) * | 2021-08-09 | 2022-02-18 | 广东工业大学 | 一种探地雷达远场定位方法、装置及计算机可读存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902830A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 |
CN104515969A (zh) * | 2013-09-29 | 2015-04-15 | 电子科技大学 | 一种基于六角形阵列的相干信号二维doa估计方法 |
CN106199505A (zh) * | 2016-06-28 | 2016-12-07 | 哈尔滨工程大学 | 一种声矢量圆阵模态域稳健方位估计方法 |
CN107132503A (zh) * | 2017-03-23 | 2017-09-05 | 哈尔滨工程大学 | 基于矢量奇异值分解的声矢量圆阵宽带相干源方位估计方法 |
CN107238829A (zh) * | 2017-06-22 | 2017-10-10 | 西北工业大学 | 一种双层圆环阵稳健超指向性波束形成方法 |
CN107907853A (zh) * | 2017-11-07 | 2018-04-13 | 中国人民解放军信息工程大学 | 一种基于均匀圆阵差分相位的单分布源doa估计方法 |
-
2018
- 2018-12-29 CN CN201811637625.3A patent/CN109541526B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104515969A (zh) * | 2013-09-29 | 2015-04-15 | 电子科技大学 | 一种基于六角形阵列的相干信号二维doa估计方法 |
CN103902830A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 |
CN106199505A (zh) * | 2016-06-28 | 2016-12-07 | 哈尔滨工程大学 | 一种声矢量圆阵模态域稳健方位估计方法 |
CN107132503A (zh) * | 2017-03-23 | 2017-09-05 | 哈尔滨工程大学 | 基于矢量奇异值分解的声矢量圆阵宽带相干源方位估计方法 |
CN107238829A (zh) * | 2017-06-22 | 2017-10-10 | 西北工业大学 | 一种双层圆环阵稳健超指向性波束形成方法 |
CN107907853A (zh) * | 2017-11-07 | 2018-04-13 | 中国人民解放军信息工程大学 | 一种基于均匀圆阵差分相位的单分布源doa估计方法 |
Non-Patent Citations (2)
Title |
---|
Guojun Jiang et al..Direction-of-arrival estimation for uniform circular arrays under small sample size.《Journal of Systems Engineering and Electronics》.2016,第27卷(第6期), * |
毛维平等.均匀圆阵相干信源二维波达方向估计.《***工程与电子技术》.2013,第35卷(第8期), * |
Also Published As
Publication number | Publication date |
---|---|
CN109541526A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9237391B2 (en) | Low noise differential microphone arrays | |
CN109541526B (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN112180329B (zh) | 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法 | |
CN108845325B (zh) | 拖曳线列阵声纳子阵误差失配估计方法 | |
CN107966677B (zh) | 一种基于空间稀疏约束的圆阵模态域方位估计方法 | |
CN109459744B (zh) | 一种实现多干扰抑制的稳健自适应波束形成方法 | |
CN111025273B (zh) | 一种畸变拖曳阵线谱特征增强方法及*** | |
CN111693971B (zh) | 一种用于弱目标检测的宽波束干扰抑制方法 | |
Xu et al. | A modified differential beamforming and its application for DOA estimation of low frequency underwater signal | |
CN107170441B (zh) | 圆环阵最优频率不变响应超指向性波束形成方法 | |
CN108761394A (zh) | 一种基于时空处理的高分辨低旁瓣去卷积谱估计方法 | |
Li et al. | Broadband constant beamwidth beamforming for suppressing mainlobe and sidelobe interferences | |
CN110457796B (zh) | 无规平面差分传声器阵列阵元位置的优化方法 | |
CN112285647A (zh) | 一种基于稀疏表示与重构的信号方位高分辨估计方法 | |
CN111175727B (zh) | 一种基于条件波数谱密度的宽带信号方位估计的方法 | |
CN108595758B (zh) | 一种任意形式传感器阵列最优宽带波束图综合方法 | |
CN109541572B (zh) | 一种基于线性环境噪声模型的子空间方位估计方法 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
Ferguson | Sharpness applied to the adaptive beamforming of acoustic data from a towed array of unknown shape | |
CN109814065B (zh) | 基于相位因子加权的波束形成方法 | |
CN109061597B (zh) | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 | |
CN110489780B (zh) | 一种由指向性声传感器组成的端射直线阵波束形成方法 | |
Chan et al. | Adaptive beamforming using uniform concentric circular arrays with frequency invariant characteristics | |
CN111610489B (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 |