CN103278799A - 基于均匀线阵改进Toeplitz化的逆波束形成方法 - Google Patents
基于均匀线阵改进Toeplitz化的逆波束形成方法 Download PDFInfo
- Publication number
- CN103278799A CN103278799A CN2013101711913A CN201310171191A CN103278799A CN 103278799 A CN103278799 A CN 103278799A CN 2013101711913 A CN2013101711913 A CN 2013101711913A CN 201310171191 A CN201310171191 A CN 201310171191A CN 103278799 A CN103278799 A CN 103278799A
- Authority
- CN
- China
- Prior art keywords
- array
- cross
- signal
- matrix
- spectrum
- 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
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
Description
技术领域
本发明涉及声源定位技术领域,特别是涉及一种基于均匀线阵改进Toeplitz化的逆波束形成方法。
背景技术
传统波束形成算法空间角度分辨率主要取阵列的孔径和信噪比,阵列孔径一旦确定其角度分辨率极限(即瑞利限)也随之确定。基于子空间分解的空间谱分析算法引入特征子空间概念,突破了角度分辨率的瑞利限,使得角度分辨率和估计精度有很大改善。但是,在强相关信号、小快拍、低信噪比条件下该类算法的性能就急剧下降,且对多源的DOA进行准确估计时需要提前对空间分布的信源数进行预估计,但在低信噪比条件下提前预估信源数很难办到。基于逆波束形成(IBF)的方位谱估计方法则无需对空间分布的信源进行预估计,能克服高分辨率方法在低信噪比条件下性能急剧下降的不足且保持较高的方位分辨率。等间隔均匀线阵的逆波束形成可分三步实现:(1)对互谱矩阵作Toeplitz平均,即Toeplitz化求互谱阵的空间分布);(2)对空间分布作相移;(3)求和输出。均匀线阵接收信号的互谱矩阵是Toeplitz矩阵。但是,当各信源之间或者信源与噪声之间存在相关性时,互谱矩阵将不再是Toeplitz矩阵,目标之间的互扰得不到有效消除, 在相邻目标间隔较近的情况下, 目标谱估计失真仍然较大,并不能达到理想的高分辨效果。
通常的做法是先用空间平滑技术解相干,这种方法主要的缺点在于:(1) 多子阵空间平滑会减小阵列的有效孔径,减少了可以估计的信源数;(2) 空间平滑技术的运算量较大,从而对实现的系统参数要求过大。Toeplitz近似化方法的基本思想是:将协方差矩阵对角线元素进行平均,构造一个Toeplitz矩阵来近似阵列信号相关矩阵,并基于所构造的Toeplitz矩阵进行目标方位估计,但是Toeplitz矩阵构造明显是一种近似的方法,存在较大的估计偏差,是一种非一致性的方位估计,同时矩阵的构造也带来了较大的计算量。
上述传统的Toeplitz化方法实质是将接收信号的协方差矩阵下三角部分各条对角线上的元素求平均,替代相应的对角线元素,存在信号空间能量向噪声空间泄露的缺陷,尤其在相干源的情况下,会产生信号与噪声空间间隔模糊。
发明内容
本发明针对现有技术的不足,提出了一种基于均匀线阵改进Toeplitz化的逆波束形成方法。
基于均匀线阵改进Toeplitz化的逆波束形成方法,具体实现步骤如下:
步骤一:设置接收麦克风阵列,确定接收阵列的信号模型。
步骤二:由信号模型确定接收线阵列输出信号的互谱矩阵。
阵列信号的互谱矩阵定义为:
步骤三:获得一组阵列输入相关值
设置线性接收阵列的第一个阵元为参考阵元,将其他阵元接收到的信号与其接收信号做相关运算,获得一组阵列输入相关值:
本发明的有益效果:本发明所提出的基于均匀线阵改进Toeplitz化的逆波束形成方法,保持了线性谱估计优良的稳定性和可靠性,在不损失阵列有效孔径条件下,信号之间的互扰得到有效消除,方位估计精度高,计算量小。
附图说明
图1是均匀线阵接收平面波示意图;
图2是CBF与IBF波束模式图。
具体实施方式
本发明的接收阵列位于目标信号源的远场,各阵元的几何尺寸远小于入射平面波波长λ,各阵元间距Δ远大于阵元尺寸且大于或者等于输入平面波半波长,即,各阵元接收的噪声是时、空不相关的高斯白噪声,其均值为零,方差为。均匀线阵接收平面波示意图1所示。
基于均匀线阵改进Toeplitz化的逆波束形成方法,具体实施方式如下:
步骤一:设置接收麦克风阵列,确定接收阵列的信号模型。
将M个各向同性麦克风阵元等间距的分布在直角坐标系中,空间K个目标平面波相对于参考方向(也称法线方向)的夹角为。则第m个阵元的输出信号为:
步骤二:由信号模型确定接收线阵列输出信号的互谱矩阵。
阵列信号的互谱矩阵定义为:
(2)
步骤三:获得一组阵列输入相关值
设置线性接收阵列的第一个阵元为参考阵元,将其他阵元接收到的信号与其接收信号做相关运算,获得一组阵列输入相关值:
(4)
(5)
表示Hadamard积,∑表示对得到的Hadamard积每个元素值求和,其中是互谱空间分布矩阵元素相应的相位补偿值。相对于原来的互谱矩阵,经过改进的Toeplitz矩阵只用到了一组阵元输入的相关值,使计算量大大减少,减少的计算量为1/M。
Claims (1)
1.基于均匀线阵改进Toeplitz化的逆波束形成方法,其特征在于,该方法具体包括以下步骤:
步骤一:设置接收麦克风阵列,确定接收阵列的信号模型;
将M个各向同性麦克风阵元等间距的分布在直角坐标系中,空间K个目标平面波相对于法线方向的夹角为 ;则第m个阵元的输出信号为:
步骤二:由信号模型确定接收线阵列输出信号的互谱矩阵;
阵列信号的互谱矩阵定义为:
步骤三:获得一组阵列输入相关值
设置线性接收阵列的第一个阵元为参考阵元,将其他阵元接收到的信号与其接收信号做相关运算,获得一组阵列输入相关值:
(4)
对于均匀阵列,其互谱矩阵中具有相同i-j值的元素,含有相同的目标信息,即
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 true CN103278799A (zh) | 2013-09-04 |
CN103278799B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107153172A (zh) * | 2017-05-08 | 2017-09-12 | 重庆大学 | 一种基于互谱优化的互谱广义逆波束形成方法 |
CN109493844A (zh) * | 2018-10-17 | 2019-03-19 | 南京信息工程大学 | 基于fir滤波器的恒定束宽波束形成方法 |
CN113281727A (zh) * | 2021-06-02 | 2021-08-20 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100813998B1 (ko) * | 2006-10-17 | 2008-03-14 | (주)펜앤프리 | 3차원 위치 추적 방법 및 장치 |
CN101192869A (zh) * | 2006-11-24 | 2008-06-04 | 中兴通讯股份有限公司 | 无线通信系统的多业务波束赋形方法 |
CN102237922A (zh) * | 2011-08-04 | 2011-11-09 | 北京北方烽火科技有限公司 | 一种波束赋形的方法和装置 |
CN102946288A (zh) * | 2012-11-23 | 2013-02-27 | 西安电子科技大学 | 基于自相关矩阵重构的压缩频谱感知方法 |
-
2013
- 2013-05-10 CN CN201310171191.3A patent/CN103278799B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100813998B1 (ko) * | 2006-10-17 | 2008-03-14 | (주)펜앤프리 | 3차원 위치 추적 방법 및 장치 |
CN101192869A (zh) * | 2006-11-24 | 2008-06-04 | 中兴通讯股份有限公司 | 无线通信系统的多业务波束赋形方法 |
CN102237922A (zh) * | 2011-08-04 | 2011-11-09 | 北京北方烽火科技有限公司 | 一种波束赋形的方法和装置 |
CN102946288A (zh) * | 2012-11-23 | 2013-02-27 | 西安电子科技大学 | 基于自相关矩阵重构的压缩频谱感知方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107153172A (zh) * | 2017-05-08 | 2017-09-12 | 重庆大学 | 一种基于互谱优化的互谱广义逆波束形成方法 |
CN107153172B (zh) * | 2017-05-08 | 2020-04-21 | 重庆大学 | 一种基于互谱优化的互谱广义逆波束形成方法 |
CN109493844A (zh) * | 2018-10-17 | 2019-03-19 | 南京信息工程大学 | 基于fir滤波器的恒定束宽波束形成方法 |
CN113281727A (zh) * | 2021-06-02 | 2021-08-20 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
CN113281727B (zh) * | 2021-06-02 | 2021-12-07 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103278799B (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730491B (zh) | 一种基于l型阵的虚拟阵列doa估计方法 | |
CN110007266B (zh) | 一种冲击噪声下的任意阵列相干源测向方法 | |
CN107315162B (zh) | 基于内插变换和波束形成的远场相干信号doa估计方法 | |
US6594201B2 (en) | System and method for localizing targets using multiple arrays | |
KR101274554B1 (ko) | 도래각 추정 방법 및 이를 이용한 배열 안테나 시스템 | |
CN106526530A (zh) | 基于传播算子的2‑l型阵列二维doa估计算法 | |
CN103018730A (zh) | 分布式子阵波达方向估计方法 | |
CN103383452A (zh) | 分布式阵列目标到达角估计方法 | |
CN111123192A (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN105005038A (zh) | 一种改进的声矢量阵相干源doa估计算法 | |
CN105158741A (zh) | 基于矩阵重构的自适应抗干扰多径多波束形成方法 | |
CN103885049A (zh) | 基于最小冗余线性稀疏子阵的米波雷达低仰角估计方法 | |
CN109932679B (zh) | 一种传感器列系统最大似然角度分辨率估计方法 | |
CN103278799A (zh) | 基于均匀线阵改进Toeplitz化的逆波束形成方法 | |
CN109491009B (zh) | 一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法 | |
CN108490428B (zh) | 一种用于抗主瓣干扰的降维子阵比相跟踪测角方法 | |
CN104156553A (zh) | 无需信源数估计的相干信号波达方向估计方法及系统 | |
Yang | On conventional beamforming and deconvolution | |
WO2009140998A1 (en) | A method for obtaining an array antenna | |
CN117471397A (zh) | 一种基于图信号处理的圆形阵列二维doa估计方法 | |
CN108872930B (zh) | 扩展孔径二维联合对角化doa估计方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN114167346B (zh) | 基于协方差矩阵拟合阵元扩展的doa估计方法及系统 | |
CN114371441A (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 |
Granted publication date: 20150422 Termination date: 20180510 |
|
CF01 | Termination of patent right due to non-payment of annual fee |