CN109541526A - 一种利用矩阵变换的圆环阵方位估计方法 - Google Patents
一种利用矩阵变换的圆环阵方位估计方法 Download PDFInfo
- Publication number
- CN109541526A CN109541526A CN201811637625.3A CN201811637625A CN109541526A CN 109541526 A CN109541526 A CN 109541526A CN 201811637625 A CN201811637625 A CN 201811637625A CN 109541526 A CN109541526 A CN 109541526A
- Authority
- CN
- China
- Prior art keywords
- array
- covariance matrix
- vector
- array manifold
- manifold vector
- 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
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)
其中,M是阵元个数,上标“T”表示转置,θ为俯仰角,φ为水平方位角,pm(θ,φ)=exp(-j2πfr sinθcos(φ-φm)/c),r为圆环阵半径,f为信号频率,c为声速;
步骤2:设计新的第m阶阵列流形向量:
其中,上标“H”表示复共轭转置,vm=M-1/2[1,ejmβ,…,ej(M-1)mβ]T特征向量,β=2π/M是相邻阵元之间的夹角;
所述E=[E0,E1,...,EM-1]T (2)
步骤3:计算第m阶新的噪声协方差矩阵:
其中,Rx为圆环阵实际接收数据计算得到的噪声协方差矩阵,由下式得到:
步骤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)下的方位谱。
图3是在信噪比为10dB(噪声为高斯白噪声),频率为f=1000Hz,引入接收信号的随机幅度误差(σg)和相位误差(σψ)的方差分别为10-3,10-2和10-1时,不同阶数m(m=0,1,…,6)下的方位谱,为避免随机误差的影响,这是200次随机实验进行平均后的结果。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明涉及一种适用于圆环阵目标方位估计的设计方法。该方法结合空间均匀噪声场中的噪声协方差矩阵的特征向量与原始阵列流形向量,设计出一组新的阵列流形向量,并定义了不同阶数的阵列流形向量的表达式,然后通过矩阵变换计算出新的噪声协方差矩阵,再利用新得到的阵列流形向量和新的噪声协方差矩阵进行目标方位估计,得到方位谱。不同阶数的方位谱的稳健性不同,阶数越低,稳健性越高,在实际应用中根据需要可灵活选择某一阶数的方位谱。其过程为:
1.计算得到圆环阵的原始阵列流形向量,由下式得到:
P(θ,φ)=[p0(θ,φ),...,pm(θ,φ),...,pM-1(θ,φ)]T (7)
其中,上标“T”表示转置,参照图1,M是阵元个数,θ为俯仰角,φ为水平方位角,pm(θ,φ)=exp(-j2πfrsinθcos(φ-φm)/c),r为圆环阵半径,f为信号频率,c为声速。
2.设计新的阵列流形向量,见下式:
E=[E0,E1,...,EM-1]T (8)
其中,上标“H”表示复共轭转置,vm=M-1/2[1,ejmβ,…,ej(M-1)mβ]T特征向量,β=2π/M是相邻阵元之间的夹角。定义第m阶阵列流形向量为:
3.计算第m阶新的噪声协方差矩阵,由下式的矩阵变换得到:
其中,Rx为圆环阵实际接收数据计算得到的噪声协方差矩阵,由下式得到:
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),r为圆环阵半径,f为信号频率,c为声速;
步骤2:设计新的第m阶阵列流形向量:
其中,上标“H”表示复共轭转置,vm=M-1/2[1,ejmβ,…,ej(M-1)mβ]T特征向量,β=2π/M是相邻阵元之间的夹角;
所述E=[E0,E1,...,EM-1]T
步骤3:计算第m阶新的噪声协方差矩阵:
其中,Rx为圆环阵实际接收数据计算得到的噪声协方差矩阵,由下式得到:
步骤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 true CN109541526A (zh) | 2019-03-29 |
CN109541526B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111562584A (zh) * | 2020-04-30 | 2020-08-21 | 北京中安智能信息科技有限公司 | 被动声纳方位历程图处理方法和装置及设备 |
CN112649787A (zh) * | 2020-11-30 | 2021-04-13 | 海鹰企业集团有限责任公司 | 一种基于低频圆环阵的目标方位估计方法 |
CN113608192A (zh) * | 2021-08-09 | 2021-11-05 | 广东工业大学 | 一种探地雷达远场定位方法、装置及计算机可读存储介质 |
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》, vol. 27, no. 6, 31 December 2016 (2016-12-31), pages 1142 - 1150 * |
毛维平等: "均匀圆阵相干信源二维波达方向估计", 《系统工程与电子技术》, vol. 35, no. 8, 31 August 2013 (2013-08-31), pages 1596 - 1601 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111562584A (zh) * | 2020-04-30 | 2020-08-21 | 北京中安智能信息科技有限公司 | 被动声纳方位历程图处理方法和装置及设备 |
CN111562584B (zh) * | 2020-04-30 | 2022-08-16 | 北京中安智能信息科技有限公司 | 被动声纳方位历程图处理方法和装置及设备 |
CN112649787A (zh) * | 2020-11-30 | 2021-04-13 | 海鹰企业集团有限责任公司 | 一种基于低频圆环阵的目标方位估计方法 |
CN112649787B (zh) * | 2020-11-30 | 2022-04-12 | 海鹰企业集团有限责任公司 | 一种基于低频圆环阵的目标方位估计方法 |
CN113608192A (zh) * | 2021-08-09 | 2021-11-05 | 广东工业大学 | 一种探地雷达远场定位方法、装置及计算机可读存储介质 |
CN113608192B (zh) * | 2021-08-09 | 2022-02-18 | 广东工业大学 | 一种探地雷达远场定位方法、装置及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109541526B (zh) | 2022-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN107315162B (zh) | 基于内插变换和波束形成的远场相干信号doa估计方法 | |
CN108845325B (zh) | 拖曳线列阵声纳子阵误差失配估计方法 | |
CN104730491A (zh) | 一种基于l型阵的虚拟阵列doa估计方法 | |
CN105005038A (zh) | 一种改进的声矢量阵相干源doa估计算法 | |
CN109541526A (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN111025273A (zh) | 一种畸变拖曳阵线谱特征增强方法及系统 | |
CN109696657A (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN107170441B (zh) | 圆环阵最优频率不变响应超指向性波束形成方法 | |
CN109541573B (zh) | 一种弯曲水听器阵列的阵元位置校准方法 | |
CN109541572B (zh) | 一种基于线性环境噪声模型的子空间方位估计方法 | |
CN108957389A (zh) | 一种实数域多通道信号目标方位估计方法 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN109814065B (zh) | 基于相位因子加权的波束形成方法 | |
CN109061594B (zh) | 一种强杂波环境下基于子阵空时自适应处理的测角方法 | |
CN115061089B (zh) | 一种声源定位方法、系统、介质、设备及装置 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN103278799B (zh) | 基于均匀线阵改进Toeplitz化的逆波束形成方法 | |
CN102183755A (zh) | 一种基于柯西-高斯模型的新型高分辨方位估计方法 | |
CN111157951B (zh) | 一种基于差分麦克风阵列的三维声源定位方法 | |
CN110632579B (zh) | 一种利用子阵波束域特征的迭代波束形成方法 | |
Liang et al. | Direction of arrival estimation under near-field interference using matrix filter | |
CN110244288A (zh) | 一种基于聚焦归位原理的水声阵列信号处理方法 | |
CN110824484A (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 |