CN103278799B - 基于均匀线阵改进Toeplitz化的逆波束形成方法 - Google Patents
基于均匀线阵改进Toeplitz化的逆波束形成方法 Download PDFInfo
- Publication number
- CN103278799B CN103278799B CN201310171191.3A CN201310171191A CN103278799B CN 103278799 B CN103278799 B CN 103278799B CN 201310171191 A CN201310171191 A CN 201310171191A CN 103278799 B CN103278799 B CN 103278799B
- Authority
- CN
- China
- Prior art keywords
- array
- mrow
- cross
- matrix
- msub
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 66
- 238000001228 spectrum Methods 0.000 claims abstract description 50
- 230000009466 transformation Effects 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 235000000748 Psidium microphyllum Nutrition 0.000 claims 1
- 244000077974 Psidium microphyllum Species 0.000 claims 1
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000010276 construction Methods 0.000 abstract description 3
- 238000013178 mathematical model Methods 0.000 abstract 1
- 238000003491 array Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明公开一种基于均匀线阵改进Toeplitz化的逆波束形成方法,本发明首先建立阵列信号的数学模型,并由互谱矩阵定义求出接收线阵列输出信号的互谱矩阵,然后设置线性接收阵列的第一个阵元为参考阵元,将其他阵元接收到的信号与其接收信号做相关运算,获得一组阵列输入相关值 。用 为行,为列构造一个Toeplitz 矩阵得到经Toeplitz化的互谱矩阵空间分布。最后对互谱矩阵空间分布的各个矩阵元素进行相位补偿求和实现逆波束形成。本发明保持了线性谱估计优良的稳定性和可靠性,在不损失阵列有效孔径条件下,信号之间的互扰得到有效消除,方位估计精度高,计算量小。
Description
技术领域
本发明涉及声源定位技术领域,特别是涉及一种基于均匀线阵改进Toeplitz化的逆波束形成方法。
背景技术
传统波束形成算法空间角度分辨率主要取阵列的孔径和信噪比,阵列孔径一旦确定其角度分辨率极限(即瑞利限)也随之确定。基于子空间分解的空间谱分析算法引入特征子空间概念,突破了角度分辨率的瑞利限,使得角度分辨率和估计精度有很大改善。但是,在强相关信号、小快拍、低信噪比条件下该类算法的性能就急剧下降,且对多源的DOA进行准确估计时需要提前对空间分布的信源数进行预估计,但在低信噪比条件下提前预估信源数很难办到。基于逆波束形成(IBF)的方位谱估计方法则无需对空间分布的信源进行预估计,能克服高分辨率方法在低信噪比条件下性能急剧下降的不足且保持较高的方位分辨率。等间隔均匀线阵的逆波束形成可分三步实现:(1)对互谱矩阵作Toeplitz平均,即Toeplitz化求互谱阵的空间分布);(2)对空间分布作相移;(3)求和输出。均匀线阵接收信号的互谱矩阵是Toeplitz矩阵。但是,当各信源之间或者信源与噪声之间存在相关性时,互谱矩阵将不再是Toeplitz矩阵,目标之间的互扰得不到有效消除, 在相邻目标间隔较近的情况下, 目标谱估计失真仍然较大,并不能达到理想的高分辨效果。
通常的做法是先用空间平滑技术解相干,这种方法主要的缺点在于:(1) 多子阵空间平滑会减小阵列的有效孔径,减少了可以估计的信源数;(2) 空间平滑技术的运算量较大,从而对实现的***参数要求过大。Toeplitz近似化方法的基本思想是:将协方差矩阵对角线元素进行平均,构造一个Toeplitz矩阵来近似阵列信号相关矩阵,并基于所构造的Toeplitz矩阵进行目标方位估计,但是Toeplitz矩阵构造明显是一种近似的方法,存在较大的估计偏差,是一种非一致性的方位估计,同时矩阵的构造也带来了较大的计算量。
上述传统的Toeplitz化方法实质是将接收信号的协方差矩阵下三角部分各条对角线上的元素求平均,替代相应的对角线元素,存在信号空间能量向噪声空间泄露的缺陷,尤其在相干源的情况下,会产生信号与噪声空间间隔模糊。
发明内容
本发明针对现有技术的不足,提出了一种基于均匀线阵改进Toeplitz化的逆波束形成方法。
基于均匀线阵改进Toeplitz化的逆波束形成方法,具体实现步骤如下:
步骤一:设置接收麦克风阵列,确定接收阵列的信号模型。
将M个各向同性麦克风阵元等间距的分布在直角坐标系中,空间K个目标平面波相对于法线方向的夹角为 。则第m个阵元的输出信号为:
(1)
式中,是阵元指向性或者阵元灵敏度,是入射目标平面波信号,是各阵元接收到的目标辐射信号相对于参考点的时间延迟。是各阵元接收到的高斯白噪声信号,式(1)是阵列信号的一般表达式。
步骤二:由信号模型确定接收线阵列输出信号的互谱矩阵。
阵列信号的互谱矩阵定义为:
(2)
式中
其中是各态历经的,互谱矩阵用信号均值计算,带入阵列信号模型表达式,则接收线阵列的互谱矩阵为:
(3)
步骤三:获得一组阵列输入相关值
设置线性接收阵列的第一个阵元为参考阵元,将其他阵元接收到的信号与其接收信号做相关运算,获得一组阵列输入相关值:
(4)
步骤四:用 为行,为列,构造一个Toeplitz 矩阵得到互谱矩阵的空间分布。
对于均匀阵列,其互谱矩阵中具有相同i-j值的元素,含有相同的目标信息,即
(5)
用代替来进行接收信号互谱矩阵Toeplitz化。用 为行,为列构造一个Toeplitz 矩阵得到互谱矩阵的空间分布
(6)
步骤五:用代替得到经过Toeplitz化的互谱阵空间分布,对互谱矩阵空间分布的各个矩阵元素进行相位补偿求和实现空间域到波数谱域的傅里叶变换,得到阵元数为M时的改进Toeplitz化逆波束形成算法。
令,互谱对所有连续的p都成立,对作傅里叶变换即可将其与频率波束谱联系起来,即:
(7)
其中,,,为接收平面波的频率,,为阵的最高工作频率,为平面波的传播速度。
空间域到波数域的傅里叶变换过程相当于对互谱矩阵元素进行相位补偿求和,即将上式离散化并用代替得:
(8)
表示Hadamard积,∑表示对得到的Hadamard积每个元素值求和,其中是互谱空间分布矩阵元素相应的相位补偿值。
本发明的有益效果:本发明所提出的基于均匀线阵改进Toeplitz化的逆波束形成方法,保持了线性谱估计优良的稳定性和可靠性,在不损失阵列有效孔径条件下,信号之间的互扰得到有效消除,方位估计精度高,计算量小。
附图说明
图1是均匀线阵接收平面波示意图;
图2是CBF与IBF波束模式图。
具体实施方式
本发明的接收阵列位于目标信号源的远场,各阵元的几何尺寸远小于入射平面波波长λ,各阵元间距Δ远大于阵元尺寸且大于或者等于输入平面波半波长,即,各阵元接收的噪声是时、空不相关的高斯白噪声,其均值为零,方差为。均匀线阵接收平面波示意图1所示。
基于均匀线阵改进Toeplitz化的逆波束形成方法,具体实施方式如下:
步骤一:设置接收麦克风阵列,确定接收阵列的信号模型。
将M个各向同性麦克风阵元等间距的分布在直角坐标系中,空间K个目标平面波相对于参考方向(也称法线方向)的夹角为。则第m个阵元的输出信号为:
(1)
式中,是阵元指向性或者阵元灵敏度,是入射目标平面波信号,是各阵元接收到的目标辐射信号相对于参考点的时间延迟。是各阵元接收到的高斯白噪声信号,式(1)是阵列信号的一般表达式,适用于窄带及宽带信号。
步骤二:由信号模型确定接收线阵列输出信号的互谱矩阵。
阵列信号的互谱矩阵定义为:
(2)
式中
其中是各态历经的,互谱矩阵用信号均值计算,带入阵列信号模型表达式,则接收线阵列的互谱矩阵为:
(3)
步骤三:获得一组阵列输入相关值
设置线性接收阵列的第一个阵元为参考阵元,将其他阵元接收到的信号与其接收信号做相关运算,获得一组阵列输入相关值:
(4)
步骤四:用 为行,为列,构造一个Toeplitz 矩阵得到互谱矩阵的空间分布
对于均匀阵列,其互谱矩阵中具有相同i-j值的元素,含有相同的目标信息,即
(5)
所以可用代替来进行接收信号互谱矩阵Toeplitz化。用 为行,为列构造一个Toeplitz 矩阵得到互谱矩阵的空间分布
(6)
步骤五:用代替得到经过Toeplitz化的互谱阵空间分布,对互谱矩阵空间分布的各个矩阵元素进行相位补偿求和实现空间域到波数谱域的傅里叶变换,得到阵元数为M时的改进Toeplitz化逆波束形成算法。
设,假定互谱对所有连续的p都成立,对作傅里叶变换即可将其与频率波束谱联系起来,即:
(7)
其中,,,为接收平面波的频率,,为阵的最高工作频率,为平面波的传播速度。
把互谱阵的空间分布通过傅里叶变换与声场的频率波数谱联系起来,间接实现了波束形成,空间域到波数域的傅里叶变换过程相当于对互谱矩阵元素进行相位补偿求和,上式离散化并用代替得:
(8)
表示Hadamard积,∑表示对得到的Hadamard积每个元素值求和,其中是互谱空间分布矩阵元素相应的相位补偿值。相对于原来的互谱矩阵,经过改进的Toeplitz矩阵只用到了一组阵元输入的相关值,使计算量大大减少,减少的计算量为1/M。
常规波束形成(CBF)与本文所改进的逆波束形成(IBF)相应的波束模式对比如图2所示,由图可以看出,IBF在互谱域上实现波数形成,主瓣约窄1/3,阵增益提高约3dB,相当于阵长增加1倍,相应的噪声抑制能力提高,具有更高的方位分辨力。采用改进的Toeplitz方法,避免了传统Toeplitz方法将相关矩阵下三角部分各条对角线上的元素求平均,替代相应的对角线元素,所造成的信号空间能量向噪声空间泄露的缺陷。
Claims (1)
1.基于均匀线阵改进Toeplitz化的逆波束形成方法,其特征在于,该方法具体包括以下步骤:
步骤一:设置接收麦克风阵列,确定接收阵列的信号模型;
将M个各向同性麦克风阵元等间距的分布在直角坐标系中,空间K个目标平面波相对于法线方向的夹角为θi,i=1,2,…,K;则第m个阵元的输出信号为:
式中,gm(θi)是阵元指向性或者阵元灵敏度,si(t)是入射目标平面波信号,τm(θi)是各阵元接收到的目标辐射信号相对于参考点的时间延迟;nm(t)是各阵元接收到的高斯白噪声信号,式(1)是阵列信号的一般表达式;
步骤二:由信号模型确定接收线阵列输出信号的互谱矩阵;
阵列信号的互谱矩阵定义为:
Rxx=E{x(t)x(t)H} (2)
式中x(t)=[x1(t),…,xm(t),…,xM(t)]T
其中x(t)是各态历经的,互谱矩阵用信号均值计算,带入阵列信号模型表达式,则接收线阵列的互谱矩阵为:
步骤三:获得一组阵列输入相关值
设置线性接收阵列的第一个阵元为参考阵元,将其他阵元接收到的信号与其接收信号做相关运算,获得一组阵列输入相关值:
rxx=[rxx[1,1],rxx[1,2],…,rxx[1,M]] (4)
步骤四:用r* xx为行,rxx为列,构造一个Toeplitz矩阵得到互谱矩阵的空间分布RT;
对于均匀阵列,其互谱矩阵Rxx中具有相同i-j值的元素,含有相同的目标信息,即
rxx[i,j]≌rxx[1,m]m-1=|i-j|,1≤m,i,j≤M (5)
用rxx[1,m]代替rxx[i,j]来进行接收信号互谱矩阵Toeplitz化;用r* xx为行,rxx为列构造一个Toeplitz矩阵得到互谱矩阵的空间分布
步骤五:用RT代替得到经过Toeplitz化的互谱阵空间分布,对互谱矩阵空间分布RT的各个矩阵元素进行相位补偿求和实现空间域到波数谱域的傅里叶变换,得到阵元数为M时的改进Toeplitz化逆波束形成算法;
令p=i-j,互谱rxx(p)对所有连续的p都成立,对rxx(p)作傅里叶变换即可将其与频率波束谱联系起来,即:
其中,u=sinθ,α=π(f/f0),f为接收平面波的频率,f0=c/2λ,为阵的最高工作频率,c为平面波的传播速度,λ为波长;
空间域到波数域的傅里叶变换过程相当于对互谱矩阵元素进行相位补偿求和,即将上式离散化并用RT代替得:
B(u)=ΣRTоA(u),-1≤u≤1 (8)
其中о表示Hadamard积,∑表示对得到的Hadamard积每个元素值求和,其中是互谱空间分布矩阵RT元素相应的相位补偿值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310171191.3A CN103278799B (zh) | 2013-05-10 | 2013-05-10 | 基于均匀线阵改进Toeplitz化的逆波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310171191.3A CN103278799B (zh) | 2013-05-10 | 2013-05-10 | 基于均匀线阵改进Toeplitz化的逆波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103278799A CN103278799A (zh) | 2013-09-04 |
CN103278799B true CN103278799B (zh) | 2015-04-22 |
Family
ID=49061374
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310171191.3A Expired - Fee Related CN103278799B (zh) | 2013-05-10 | 2013-05-10 | 基于均匀线阵改进Toeplitz化的逆波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103278799B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107153172B (zh) * | 2017-05-08 | 2020-04-21 | 重庆大学 | 一种基于互谱优化的互谱广义逆波束形成方法 |
CN109493844A (zh) * | 2018-10-17 | 2019-03-19 | 南京信息工程大学 | 基于fir滤波器的恒定束宽波束形成方法 |
CN113281727B (zh) * | 2021-06-02 | 2021-12-07 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其*** |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100813998B1 (ko) * | 2006-10-17 | 2008-03-14 | (주)펜앤프리 | 3차원 위치 추적 방법 및 장치 |
CN101192869B (zh) * | 2006-11-24 | 2011-06-22 | 中兴通讯股份有限公司 | 无线通信***的多业务波束赋形方法 |
CN102237922B (zh) * | 2011-08-04 | 2014-06-04 | 北京北方烽火科技有限公司 | 一种波束赋形的方法和装置 |
CN102946288B (zh) * | 2012-11-23 | 2014-08-06 | 西安电子科技大学 | 基于自相关矩阵重构的压缩频谱感知方法 |
-
2013
- 2013-05-10 CN CN201310171191.3A patent/CN103278799B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103278799A (zh) | 2013-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730491A (zh) | 一种基于l型阵的虚拟阵列doa估计方法 | |
CN106443587B (zh) | 一种高分辨率的快速反卷积声源成像算法 | |
CN106093921B (zh) | 基于稀疏分解理论的声矢量阵宽带测向方法 | |
CN103323832B (zh) | 一种相控阵三维摄像声纳***换能器阵列的幅相误差校正方法 | |
CN103353588B (zh) | 基于天线均匀平面阵的二维波达方向角估计方法 | |
CN104730513A (zh) | 一种分级子阵聚焦mvdr波束形成方法 | |
CN111736118B (zh) | 一种线列阵阵列扩展方法 | |
CN110837074A (zh) | 一种基于数字波束形成的多同频信源相位干涉仪测向方法 | |
CN103278799B (zh) | 基于均匀线阵改进Toeplitz化的逆波束形成方法 | |
CN105158751A (zh) | 一种声矢量阵快速方位估计方法 | |
CN106908754B (zh) | L型声矢量传感器阵列esprit解相干参数估计方法 | |
CN109491009B (zh) | 一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法 | |
CN111175727B (zh) | 一种基于条件波数谱密度的宽带信号方位估计的方法 | |
CN103323810B (zh) | 一种l阵方位角和俯仰角配对的信号处理方法 | |
CN109375197B (zh) | 一种小尺寸矢量阵低频散射校正方法 | |
CN107064862B (zh) | 电磁矢量传感器线阵解相干music参数估计方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN107241131B (zh) | 一种利用信号非圆特性的波束形成方法 | |
CN114371441A (zh) | 虚拟阵列波达方向估计方法、装置、产品及存储介质 | |
CN111830458B (zh) | 一种平行线阵单快拍二维测向方法 | |
CN106950530B (zh) | 基于互协方差的声矢量传感器阵列解相干参数估计方法 | |
CN110632579B (zh) | 一种利用子阵波束域特征的迭代波束形成方法 | |
Zhang et al. | Research on linear sparse sonar array for targets azimuth bearing estimation | |
CN105785315B (zh) | 一种直线阵协方差矩阵主对角线重构的测角方法 | |
CN109061564B (zh) | 基于高阶累积量的简化近场定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150422 Termination date: 20180510 |