CN110489780A - 一种由指向性声传感器组成的端射直线阵波束形成方法 - Google Patents
一种由指向性声传感器组成的端射直线阵波束形成方法 Download PDFInfo
- Publication number
- CN110489780A CN110489780A CN201910595821.7A CN201910595821A CN110489780A CN 110489780 A CN110489780 A CN 110489780A CN 201910595821 A CN201910595821 A CN 201910595821A CN 110489780 A CN110489780 A CN 110489780A
- Authority
- CN
- China
- Prior art keywords
- array
- line array
- noise
- beam forming
- directive property
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种由指向性声传感器组成的端射直线阵波束形成方法,直接计算波束形成过程中所需要的噪声协方差矩阵,以波束形成的加权向量:w=CTD‑2C·P(θ0)计算直线阵的波束形成响应即波束图:B(θ)=wHP(θ)。具有简单易用、计算速度快的优点。该方法首先使两个传感器同时指向选取的典型方向(端射方向),然后计算传感器间的空间相关性的函数,由此组成噪声协方差矩阵,最后进行波束形成进行输出。
Description
技术领域
本发明属于声学阵列信号处理、声纳信号处理等领域,涉及一种由指向性声传感器组成的端射直线阵波束形成方法,适用于目标信号检测、噪声消除等领域。
背景技术
由具有指向性的声传感器(比如声矢量传感器)组成的阵列(比如直线阵、圆环阵、平面阵和立体阵等)在雷达、声纳、空气中传声器等系统中具有广泛的应用,常 被用来进行定向声信号增强、噪声抑制、微弱目标检测等。早期的由指向性的声传感 器组成阵列的研究有:(1)文献1“Hawkes M,Nehorai A,et al.Acoustic vector-sensorcorrelations in ambient noise[J].IEEE J.Ocean.Eng.,2001,26(3):337-347.”公开的方法, 计算了声矢量传感器各个分量的空间相关性,只给出了质点振速传感器直线阵的波束 形成方法,没有给出与声压分量的联合估计;(2)文献2“D’Spain G L,Luby J C,WilsonG R,et al.Vector sensors and vector sensor line arrays:Comments on optimalarray gain and detection[J].J.Acoust.Soc.Am,2006,120(1):171-185.”在文献1的基础上进一步分 析了矢量传感器各个分量的空间相关性,给出了相关性的解析表达式,但是仅仅停留 在单自由度的矢量传感器(比如质点振速传感器)的分析。(3)文献3“王绪虎,陈建 峰,韩晶等.环境噪声场中二阶组合水听器的空间相关性[J].振动与冲击,2013, 32(21).”给出了矢量水听器的一阶、二阶各个分量之间的相关性,但其表达式涉及到多 个夹角,较为复杂不易计算,并且未涉及关于具有指向性的声传感器的空间相关性和 波束形成。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种由指向性声传感器组成的端射直线阵波束形成方法,解决指向性声传感器阵列在波束形成中的空间噪声相关性计算复 杂、只能通过积分求解且没有解析解的缺陷。本方法具有较为简单直接的关于空间噪 声相关性计算方式,并且波束形成能够有效抑制栅瓣以及非期望方向到达的干扰和噪 声信号。本发明假设多个指向性声传感器位于坐标系的z轴上,同时指向轴线方向, 即可组成端射线列阵,这在麦克风、声纳系统等实际应用具有广泛的使用。
技术方案
一种由指向性声传感器组成的端射直线阵波束形成方法,其特征在于步骤如下:
步骤1:声传感器等间距组成的直线阵在空间均匀噪声场中由噪声相关性组成的噪声协方差矩阵R,其第m行、第m'列的元素为:
其中,m≥m',a、b为实数且a+b=1,0<a<1,0<b<1,表示虚数单 位,x=k·(m-m')·d,k=2πf/c(f为信号频率,c为声速,π为圆周率)为波数,d 为相邻两个传感器间的距离;
步骤2:当m<m'时,噪声协方差矩阵的第m行、第m'列的元素为:
步骤3:计算波束形成的加权向量:
w=CTD-2C·P(θ0)
其中:上标“T”表示转置,θ0表示入射信号的俯仰角,且θ0=0°; P(θ0)=[p1(θ0),p2(θ0),...,pm(θ0),...,pM(θ0)]T为指向θ0的阵列流形向量,M为阵元个数, pm(θ0)=(a+bcosθ0)·exp[-jk(m-1)dcosθ0],D为一个M×M对角矩阵,D-2用来对输 入噪声进行归一化,C是一个正交变换矩阵;
所述C=[c0 c1 … ck … cM-1]T
其中:ck=[ck0 ck1 … ckk 0 … 0]T
其中,上标“*”表示共轭,
步骤4:计算直线阵的波束形成响应即波束图:
B(θ)=wHP(θ)
其中:上标“H”表示共轭转置,P(θ)=[p1(θ),p2(θ),...,pm(θ),...,pM(θ)]T为指向角 度θ的阵列流形向量,θ表示入射信号的俯仰角, pm(θ)=(a+bcosθ)·exp[-jk(m-1)dcosθ]。
有益效果
本发明提出的一种由指向性声传感器组成的端射直线阵波束形成方法,直接计算波束形成过程中所需要的噪声协方差矩阵,具有简单易用、计算速度快的优点。该方 法首先使两个传感器同时指向选取的典型方向(端射方向),然后计算传感器间的空间 相关性的函数,由此组成噪声协方差矩阵,最后进行波束形成进行输出。
有益效果体现在:
1.本发明公开的波束形成设计方法是针对具有指向性传感器之间的相关性,而文献1公开的方法只是处理无指向性的声压传感器。
2.本发明公开的方法设计出的相关性计算方法是关于具有指向性的传感器,而文献2和文献3公开的方法只是针对单自由度的矢量传感器,并没有对更复杂的指向性 传感器的描述。
3.本发明公开的波束形成设计方法可对声传感器间的相关性通过简单的算数求得,具有解析解的形式,避免了繁琐的不精确的积分过程,而文献4的计算涉及到的 关于矢量传感器的一阶、二阶相关性的计算较为复杂、参数过多,不能直接计算。
附图说明
图1是两个具有指向性的声传感器同时指向z轴的轴线方向的示意图。
图2有指向性的声传感器阵元在不同a和b取值时的空间接收响应,图2(a)a=0.2,b=0.8;图2(b)a=0.5,b=0.5;图2(c)a=0.8,b=0.2。从图2中可以看出, 阵元在不同a和b的取值时有不同的空间接收响应。
图3是以第1号传感器为参考得到的两个传感器(第1号和第2号)指向z轴轴 线(端射)方向的相关性的值,作为示例,a=0.2,b=0.8,图3(a)是相关性值的实部, 图3(b)是相关性值的虚部。
图4是以第1号传感器为参考得到的两个传感器(第1号和第2号)指向z轴轴 线(端射)方向的相关性的值,作为示例,a=b=0.5,图4(a)是相关性值的实部,图 4(b)是相关性值的虚部。
图5是以第1号传感器为参考得到的两个传感器(第1号和第2号)指向z轴轴 线(端射)方向的相关性的值,作为示例,a=0.8,b=02,图5(a)是相关性值的实部, 图5(b)是相关性值的虚部。
图6是由M=7个声传感器阵元且x=6时计算得到的波束图,图6(a)中:a=0.2, b=0.8;图6(b)中:a=b=0.5;图6(c)中a=0.8,b=0.2。虚线为传统声压直线阵,实 线为阵元有指向性的直线阵
图7是由M=7个声传感器阵元计算得到的阵增益曲线,虚线为传统声压直线阵,实线为阵元有指向性的直线阵图7(a)中:a=0.2,b=0.8;图7(b)中:a=b=0.5;图7(c)中:a=0.8,b=0.2。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
具有指向性的声传感器组成的端射直线阵的波束形成方法步骤如下:
步骤1、求声传感器等间距组成的直线阵在空间均匀噪声场中由噪声相关性组成的噪声协方差矩阵R,其第m行、第m'列的元素由下式计算::
其中,m≥m',a、b为实数且a+b=1,0<a<1,0<b<1,表示虚数单位, x=k·(m-m')·d,k=2πf/c(f为信号频率,c为声速,π为圆周率)为波数,d为相 邻两个传感器间的距离。
步骤2、当m<m'时,噪声协方差矩阵的第m行、第m'列的元素由下式计算:
图3是以第1号传感器为参考得到的两个传感器(第1号和第2号)指向z轴轴 线(端射)方向的相关性的值,即是噪声协方差矩阵中的第1行、第2列上的元素值, 参数a=0.2,b=0.8,图3(a)是相关性值的实部,图3(b)是相关性值的虚部。
图4是以第1号传感器为参考得到的两个传感器(第1号和第2号)指向z轴轴 线(端射)方向的相关性的值,即是噪声协方差矩阵中的第1行、第2列上的元素值, 参数a=b=0.5,图4(a)是相关性值的实部,图4(b)是相关性值的虚部。
图5是以第1号传感器为参考得到的两个传感器(第1号和第2号)指向z轴轴 线(端射)方向的相关性的值,即是噪声协方差矩阵中的第1行、第2列上的元素值, 参数a=0.8,b=0.2,图5(a)是相关性值的实部,图5(b)是相关性值的虚部。
从图3~图5可以看出,随着a值的增大,b值的减小,相关性值的实部越来越大, 逐渐向传统声压传感器靠近。
步骤3、计算波束形成的加权向量,由下式得到:
w=CTD-2C·P(θ0) (10)
其中,上标“T”表示转置,θ0表示入射信号的俯仰角,且θ0=0°, P(θ0)=[p1(θ0),p2(θ0),...,pm(θ0),...,pM(θ0)]T为指向θ0的阵列流形向量,M为阵元个数, pm(θ0)=(a+bcosθ0)·exp[-jk(m-1)dcosθ0],公式(10)中的D为一个M×M对角矩阵,D-2用来对输入噪声进行归一化,C是一个正交变换矩阵,具有如下形式,
C=[c0 c1 … ck … cM-1]T (11)
其中,ck=[ck0 ck1 … ckk 0 … 0]T,cki可由下式计算,
其中,上标“*”表示共轭,公式(12)中,有
步骤4、计算直线阵的波束形成响应(波束图),由下式得到:
B(θ)=wHP(θ) (14)
其中,上标“H”表示共轭转置,P(θ)=[p1(θ),p2(θ),...,pm(θ),...,pM(θ)]T为指向角度θ的 阵列流形向量,θ表示入射信号的俯仰角,pm(θ)=(a+bcosθ)·exp[-jk(m-1)dcosθ]。
图6是M=7个声传感器阵元且x=6时计算得到的波束图,蓝色虚线为传统声压 直线阵,红色实线为阵元有指向性的直线阵,图6(a)中:a=0.2,b=0.8;图6(b)中:a =b=0.5;图6(c)中a=0.8,b=0.2。两种阵列均指向端射方向,从图6中可以看到, 传统的声压直线阵的波束图已经出现栅瓣,而阵元有指向性的直线阵成功消除了栅瓣, 并且图6(b)很好的抑制了来自90°~270°方向到达的干扰信号。
图7是由M=7个声传感器阵元计算得到的随x变化的阵增益曲线,两种阵列同 时指向端射方向,蓝色线为传统声压直线阵,红色线为阵元有指向性的直线阵,图7(a) 中:a=0.2,b=0.8;图7(b)中:a=b=0.5;图7(c)中:a=0.8,b=0.2。从图7中可以 看到,有指向性阵元组成的直线阵的阵增益在x>3时远远大于传统声压直线阵;在x 值较小时,图7(a)和图7(b)也比声压阵的阵增益要高1dB以上。随着a值的增大,b 值的减小,有指向性阵元组成的直线阵的阵增益逐渐接近于传统声压直线阵,与理论 相符。在x=4左右,两种阵列的阵增益差值最大,而此时阵元间距与波长的比值大于 0.5,这提示我们在采用有指向性的阵元组成阵列时应注意阵元间距与传统声压阵的半 波长是不一样的。
Claims (1)
1.一种由指向性声传感器组成的端射直线阵波束形成方法,其特征在于步骤如下:
步骤1:声传感器等间距组成的直线阵在空间均匀噪声场中由噪声相关性组成的噪声协方差矩阵R,其第m行、第m'列的元素为:
其中,m≥m',a、b为实数且a+b=1,0<a<1,0<b<1,表示虚数单位,x=k·(m-m')·d,k=2πf/c为波数:f为信号频率,c为声速,π为圆周率;d为相邻两个传感器间的距离;
步骤2:当m<m'时,噪声协方差矩阵的第m行、第m'列的元素为:
步骤3:计算波束形成的加权向量:
w=CTD-2C·P(θ0)
其中:上标“T”表示转置,θ0表示入射信号的俯仰角,且θ0=0°;P(θ0)=[p1(θ0),p2(θ0),...,pm(θ0),...,pM(θ0)]T为指向θ0的阵列流形向量,M为阵元个数,pm(θ0)=(a+b cosθ0)·exp[-jk(m-1)d cosθ0],D为一个M×M对角矩阵,D-2用来对输入噪声进行归一化,C是一个正交变换矩阵;
所述C=[c0 c1 … ck … cM-1]T
其中:ck=[ck0 ck1 … ckk 0 … 0]T
其中,上标“*”表示共轭,
步骤4:计算直线阵的波束形成响应即波束图:
B(θ)=wHP(θ)
其中:上标“H”表示共轭转置,P(θ)=[p1(θ),p2(θ),...,pm(θ),...,pM(θ)]T为指向角度θ的阵列流形向量,θ表示入射信号的俯仰角,pm(θ)=(a+b cosθ)·exp[-jk(m-1)d cosθ]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910595821.7A CN110489780B (zh) | 2019-07-03 | 2019-07-03 | 一种由指向性声传感器组成的端射直线阵波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910595821.7A CN110489780B (zh) | 2019-07-03 | 2019-07-03 | 一种由指向性声传感器组成的端射直线阵波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110489780A true CN110489780A (zh) | 2019-11-22 |
CN110489780B CN110489780B (zh) | 2022-05-03 |
Family
ID=68546551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910595821.7A Active CN110489780B (zh) | 2019-07-03 | 2019-07-03 | 一种由指向性声传感器组成的端射直线阵波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110489780B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090022336A1 (en) * | 2007-02-26 | 2009-01-22 | Qualcomm Incorporated | Systems, methods, and apparatus for signal separation |
CN103438987A (zh) * | 2013-08-30 | 2013-12-11 | 西北工业大学 | 基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法 |
CN103902831A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种基于模态分解与综合的超指向性波束形成方法 |
CN105372644A (zh) * | 2015-11-03 | 2016-03-02 | 西北工业大学 | 一种基于动态重修正的自适应波束形成方法及系统 |
CN108663668A (zh) * | 2018-05-18 | 2018-10-16 | 西安电子科技大学 | 基于iaa的干扰加噪声协方差矩阵重构稳健波束形成方法 |
CN109459744A (zh) * | 2018-12-19 | 2019-03-12 | 西北工业大学 | 一种实现多干扰抑制的稳健自适应波束形成方法 |
-
2019
- 2019-07-03 CN CN201910595821.7A patent/CN110489780B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090022336A1 (en) * | 2007-02-26 | 2009-01-22 | Qualcomm Incorporated | Systems, methods, and apparatus for signal separation |
CN103438987A (zh) * | 2013-08-30 | 2013-12-11 | 西北工业大学 | 基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法 |
CN103902831A (zh) * | 2014-04-11 | 2014-07-02 | 西北工业大学 | 一种基于模态分解与综合的超指向性波束形成方法 |
CN105372644A (zh) * | 2015-11-03 | 2016-03-02 | 西北工业大学 | 一种基于动态重修正的自适应波束形成方法及系统 |
CN108663668A (zh) * | 2018-05-18 | 2018-10-16 | 西安电子科技大学 | 基于iaa的干扰加噪声协方差矩阵重构稳健波束形成方法 |
CN109459744A (zh) * | 2018-12-19 | 2019-03-12 | 西北工业大学 | 一种实现多干扰抑制的稳健自适应波束形成方法 |
Non-Patent Citations (5)
Title |
---|
SHEFENG YAN ET AL: "Robust supergain beamforming for circular array via second-order cone programming", 《APPLIED ACOUSTICS》 * |
YIXIN YANG ET AL: "A Superdirective Beamforming Method for Linear Sensor Arrays", 《OCEANS 2014 - TAIPEI》 * |
付留芳: "基于光纤矢量水听器的海洋环境噪声特性研究", 《中国优秀硕士学位论文全文数据库电子期刊 基础科学辑》 * |
朱少豪 等: "直线形质点振速传感器阵列超指向性研究", 《声学技术》 * |
汪勇: "传感器阵列超指向性原理及方法研究", 《中国博士学位论文全文数据库电子期刊 信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110489780B (zh) | 2022-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730491B (zh) | 一种基于l型阵的虚拟阵列doa估计方法 | |
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN106405487B (zh) | 一种通用基于扩展esprit技术的空间谱估计方法 | |
CN107884741A (zh) | 一种多球阵列多宽带声源快速定向方法 | |
CN108225536B (zh) | 基于水听器幅度与相位自校准的稳健自适应波束形成方法 | |
Jo et al. | Parametric direction-of-arrival estimation with three recurrence relations of spherical harmonics | |
CN110346752B (zh) | 基于互质稀疏阵的无模糊测向方法 | |
CN112269164A (zh) | 深海可靠声路径下基于干涉结构匹配处理弱目标定位方法 | |
CN114527427B (zh) | 一种基于球形麦克风阵列的低频波束形成声源定位方法 | |
CN110082712B (zh) | 一种声矢量圆阵相干目标方位估计方法 | |
CN105005038A (zh) | 一种改进的声矢量阵相干源doa估计算法 | |
CN109541548A (zh) | 一种基于匹配场的空气声呐定位方法 | |
CN111323784A (zh) | 一种近场分布式声源定位方法 | |
CN110687538A (zh) | 一种基于近场聚焦的超波束形成方法 | |
CN109375227A (zh) | 一种解卷积波束形成三维声成像方法 | |
CN109932679B (zh) | 一种传感器列系统最大似然角度分辨率估计方法 | |
Jing et al. | Sound source localisation using a single acoustic vector sensor and multichannel microphone phased arrays | |
CN109491009B (zh) | 一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法 | |
Altes | Angle estimation and binaural processing in animal echolocation | |
CN109541526A (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
Ma et al. | Spatiotemporal two-dimensional deconvolution beam imaging technology | |
CN111812607B (zh) | 基于波束空间的米波mimo雷达低仰角估计方法 | |
CN109375197B (zh) | 一种小尺寸矢量阵低频散射校正方法 | |
CN110489780A (zh) | 一种由指向性声传感器组成的端射直线阵波束形成方法 | |
CN109669172A (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 |