CN102721943B - 基于l型传感器阵列的非相干信号二维波达方向跟踪方法及装置 - Google Patents
基于l型传感器阵列的非相干信号二维波达方向跟踪方法及装置 Download PDFInfo
- Publication number
- CN102721943B CN102721943B CN201210190774.6A CN201210190774A CN102721943B CN 102721943 B CN102721943 B CN 102721943B CN 201210190774 A CN201210190774 A CN 201210190774A CN 102721943 B CN102721943 B CN 102721943B
- Authority
- CN
- China
- Prior art keywords
- overbar
- theta
- phi
- zeta
- delta
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于L型传感器阵列的非相干信号二维波达方向跟踪方法及装置,该方法利用RLS算法自适应估计噪声子空间从而建立瞬态目标函数,再运用近似牛顿迭代算法求解所得目标函数,从而自适应更新角度。通过充分探究信号的运动模型,利用伦伯格状态观测器解决信号角度交叉的问题,同时避免了前后时刻间参数的关联运算。另一方面,该方法将二维波达方向估计问题转化为一维的估计问题,实现参数降维从而减小了运算量,提供自动配对的信号方位角和仰角估计值。
Description
技术领域:
本发明属于通信领域,涉及二维波达方向的估计与跟踪,具体涉及一种基于L型传感器阵列的非相干信号二维波达方向跟踪方法及装置。
背景技术:
阵列信号处理的实际应用包括雷达、声纳、通信等,多入射信号的二维(2-D)波达方向(DOA)(例如,方位角和仰角)必须从众多的测量中估计得到。随着维数的增加,二维估计问题比一维更加复杂,并且估计的方位角必须与同一入射信号的仰角对应起来。最大似然法(ML)为二维问题提供了一种最优解,但是在参数空间计算繁重的多维搜索常常需要获得似然函数的全局最大值。由于他们比最大似然法的实现简单,许多基于特征结构和平面阵列(如,均匀圆形阵列和矩形阵列)的次优的方法已经由解决一维问题发展为解决二维波达角度估计问题。
二维波达方向估计的计算复杂度通常会受到阵列几何结构的影响,这些特殊的构造可以经过研究发展为计算高效的DOA估计方法。近年来,一种L型阵列受到了很大的关注,这种阵列由两列均匀线性阵列(ULA)直角连接组成,与传统的平面阵列相比它具有一些几何结构和实现方面的优势,且比其他具有两列或更多ULA的简单结构的平面阵列具有更高的估计精确度。很多需要特征分解的二维DOA估计方法都是使用的L阵,其中阵列放置在x-y平面或者x-z平面。传统的二维估计问题可分解为两个独立的一维估计问题,并且经过深入研究的一维基于子空间的方法可以直接应用于避免计算复杂的多维搜索与优化。ESPRIT这类方法提供了通过矩阵的特征值对方位角和仰角自配对的封闭型的估计,其它方法都是通过配对过程。然而,现存的配对技术的计算代价很高,并且经常不能提供正确的配对结果,所以,配对失败会造成估计失败。对于上述估计失败的措施在上面提到的二维估计方法中还没有考虑。另外,这些方法都存在计算密集和费时的特征分解过程,因此不能应用到实时过程中。
另一方面,众多的一维波达方向(如,方位角)跟踪的算法都已经载入文献,但是多运动目标的二维跟踪问题还没有经过广泛深入的研究,且难点主要在方位角和仰角的配对以及在连续两个时间间隔下这些估计的结合。通过扩展一维递归算法或基于渐变的迭代搜索算法到二维场景下,目前提出的两种不同的二维DOA跟踪算法利用的是矩形平面阵列。即使方位角和仰角的连接关系就在估计本身中,它们需要特征值分解的MUSIC方法来确定噪声方差、信号协方差矩阵、仰角方位角的初步估计。因此它们的跟踪效果受参数精度的影响,当信噪比(SNR)较低时估计性能会降低。另外,联合二维DOA跟踪算法在估计DOA及在每个时间间隔最小化关联仰角和方位角的估计过程中仍旧存在计算计算量很大的特征值分解过程,即使子空间的更新是通过使用低秩自适应滤波器完成的。
发明内容:
本发明提供了基于L型传感器阵列的非相干信号二维波达方向跟踪方法及装置,在两个连续的时间间隔内的方位角和仰角的关联可由动态模型和伦贝格状态观测器解决。。该方法无需特征值分解,计算简单有效,同时能够解决入射信号轨迹有交叉情况下的信号二维波达方向跟踪。
具体如下:
一种基于L型传感器阵列的非相干信号二维波达方向跟踪方法,包括以下步骤:
2)初始化伦伯格观测器,由已估计的状态向量计算定义的状态向量以及定义的方向角;
3)在方向向量更新的间隔内,估计瞬时互协方差矩阵;
4)计算估计的正交投影;
5)预测方向角之间的增量;
一种方向估计与跟踪装置,用于利用阵列天线来估计和跟踪入射信号的到达方向,该阵列天线具有两列成直角的阵列传感器组成,每列传感器以相同的间隔设置在不同空间位置上的多个天线阵子,该装置包括:
初始值估计装置,估计空间缓慢移动的信号的波达方向初始值;
定时装置,记录状态更新的索引间隔;
瞬时互协方差矩阵计算装置,在更新的时间间隔内计算瞬时互协方差矩阵;
增量计算装置,利用迭代的方法计算方位角和仰角的增量;
状态向量更新装置,利用伦伯格观测器对状态向量进行更新并估计出新的仰角和方位。
本发明的有益效果是:
相比于已有的跟踪算法,本算法无需特征值分解,计算简单有效,同时能够解决入射信号轨迹有交叉情况下的信号二维波达方向跟踪。因此这种在线算法对于波达方向随时间变化的信号源有很强的跟踪能力,并且所估计的仰角和方位角和真实值十分接近。
附图说明:
图1是表示根据本发明实施例的二维方向估计方法的L型阵列的几何结构图。
图2(a)是表示根据本发明实施例的仰角θk和方位角φk的可行域。
具体实施方式:
下面结合附图对本发明做详细描述。
数据模型以及问题描述:
如图1所示,L型传感器阵列放置在x-z平面,并由两列均匀线阵组成,每列均匀线阵包含M个全向的传感器,间距为d,在坐标系原点z0处的传感器作为每列均匀线阵的参考,z0和x1传感器间距也是d。假设p个不相干的远场窄带信号{sk(n)},波长λ,从不同的仰角和方位角{θk(n),φk(n)}入射到阵列上。如图1所示,仰角θk(n)和方位角φk(n)定义为相对于z轴或者x轴的顺时针测量,投影方位角为在x-y平面内相对于x轴逆时针测量。这里两列均匀线阵接收到的信号可表示为
z(n)=A(θ(n))s(n)+wz(n) (1)
x(n)=A(φ(n))s(n)+wx(n) (2)
这里A(θ(n))和A(φ(n))式沿z轴和x放置的线阵的响应矩阵,wz(n)和wx(n)是对应的噪声矢量。另外定义和时刻的速度为和其相应的加速度为和动态模型的状态向量可以表示为 及我们可将沿z轴放置的均匀线阵分割为两个不相重叠的子阵列,分别包括p个传感器和M-p个传感器,式(1)中的z(n)可以表示为其中z1(n)和z2(n)为子阵列的接收数据。同时假设θk(n)和φk(n)随时间缓慢变化(相对于采样率1/Ts),即满足 其中在DOA更新的时间间隔T上的快拍数Ns已知,例如T=NsTs。因此,二维DOA跟踪问题是由Ns个快拍数据{z(n),x(n)}估计仰角-方位角对并使得同一个入射信号现有估计和先前估计保持正确的数据关联,其中
下面概要描述根据本发明实施例的联合方位角和仰角的波达方向估计。
其中B(θ)=diag(a2(θ),IM),a2(θ)=[ej2πpdcosθ/λ,…,ej2π(M-1)dcosθ/λ],k=1,2,…,p。定义 构造多项式 求解多项式 的根,计算这些根中最接近单位圆的根的相位,而由z=ej2πdcosφ/λ计算相应的仰角估计值
这里gθk,gφk为观测增益,并满足矩阵F-gθkcT和F-gφkcT的特征值在单位圆内。
Claims (6)
1.一种基于L型传感器阵列的非相干信号二维波达方向跟踪方法,其特征在于:包括以下步骤:
其中B(θ)=diag(a2(θ),IM),;a2(θ)=[ej2πpdcosθ/λ,…,ej2π(M-1)dcosθ/λ],k=1,2,…,p;定义构造多项式求解多项式pφk(z)=0的根,计算这些根中最接近单位圆的根的相位,而由z=ej2πdcosφ/λ计算相应的仰角估计值
2)初始化伦伯格观测器,由已估计的状态向量计算定义的状态向量以及定义的方向角;
3)在方向向量更新的间隔内,估计瞬时互协方差矩阵;
4)计算估计的正交投影;
5)预测方向角之间的增量;
2.根据权利要求1所述的方法,其特征在于:针对单个或多个运动或静止的信号进行自适应二维波达方向的跟踪估计。
6.一种方向估计与跟踪装置,用于利用阵列天线来估计和跟踪入射信号的到达方向,包括阵列天线,该阵列天线具有两列成直角的阵列传感器组成,每列传感器以相同的间隔设置在不同空间位置上的多个天线阵子,其特征在于,该装置还包括:
初始值估计装置,估计空间缓慢移动的信号的波达方向初始值;
定时装置,记录状态更新的索引间隔;
瞬时互协方差矩阵计算装置,在更新的时间间隔内计算瞬时互协方差矩阵;
增量计算装置,利用迭代的方法计算方位角和仰角的增量;
状态向量更新装置,利用伦伯格观测器对状态向量进行更新并估计出新的仰角和方位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210190774.6A CN102721943B (zh) | 2012-06-11 | 2012-06-11 | 基于l型传感器阵列的非相干信号二维波达方向跟踪方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210190774.6A CN102721943B (zh) | 2012-06-11 | 2012-06-11 | 基于l型传感器阵列的非相干信号二维波达方向跟踪方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102721943A CN102721943A (zh) | 2012-10-10 |
CN102721943B true CN102721943B (zh) | 2014-03-12 |
Family
ID=46947759
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210190774.6A Active CN102721943B (zh) | 2012-06-11 | 2012-06-11 | 基于l型传感器阵列的非相干信号二维波达方向跟踪方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102721943B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049003B (zh) * | 2012-10-19 | 2016-03-30 | 西安交通大学 | 一种基于平行均匀线阵的相干信号二维波达角度跟踪方法及装置 |
CN103607217B (zh) * | 2013-10-11 | 2015-07-01 | 中国人民解放军重庆通信学院 | 基于子空间跟踪的脉冲uwb通信***干扰抑制方法 |
US9759807B2 (en) * | 2013-10-25 | 2017-09-12 | Texas Instruments Incorporated | Techniques for angle resolution in radar |
KR102115950B1 (ko) | 2013-11-07 | 2020-06-05 | 삼성전자주식회사 | 보행 보조 로봇 및 보행 보조 로봇의 제어 방법 |
CN103954931B (zh) * | 2014-04-28 | 2016-06-08 | 西安交通大学 | 一种远场和近场混合信号源的定位方法 |
CN104898085B (zh) * | 2015-05-14 | 2017-11-17 | 电子科技大学 | 一种极化敏感阵列参数估计的降维music算法 |
GB2539722B (en) * | 2015-06-25 | 2021-10-13 | Airspan Ip Holdco Llc | Bearing calculation |
CN109655783B (zh) * | 2018-12-26 | 2023-07-21 | 西安云脉智能技术有限公司 | 一种感应器列阵来波方向的估计方法 |
CN110018438B (zh) * | 2019-04-23 | 2020-09-25 | 北京邮电大学 | 一种波达方向估计方法和装置 |
CN112327244B (zh) * | 2020-10-22 | 2022-06-24 | 中国电子科技集团公司第五十四研究所 | 一种基于l型阵列的二维非相干分布式目标参数估计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5459668A (en) * | 1993-06-04 | 1995-10-17 | University Of Southern California | Method and apparatus for signal analysis employing a virtual cross-correlation computer |
ATE463749T1 (de) * | 2001-04-27 | 2010-04-15 | Mitsubishi Elec R&D Ct Europe | Verfahren zur einfallsrichtungsschätzung |
KR100703322B1 (ko) * | 2004-05-07 | 2007-04-03 | 삼성전자주식회사 | 배열 안테나 시스템을 위한 빔 형성 장치 및 방법 |
-
2012
- 2012-06-11 CN CN201210190774.6A patent/CN102721943B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102721943A (zh) | 2012-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102721943B (zh) | 基于l型传感器阵列的非相干信号二维波达方向跟踪方法及装置 | |
CN103353596B (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN102707258B (zh) | 基于l型传感器阵列的信号方位角和仰角联合估计方法 | |
CN103049003B (zh) | 一种基于平行均匀线阵的相干信号二维波达角度跟踪方法及装置 | |
CN107015191A (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
CN106405487B (zh) | 一种通用基于扩展esprit技术的空间谱估计方法 | |
CN109239646B (zh) | 一种冲击噪声环境下连续量子水蒸发的二维动态测向方法 | |
CN102253363A (zh) | 基于l型阵列的相干信号二维波达方向估计装置及其方法 | |
CN101644760B (zh) | 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法 | |
CN107656239B (zh) | 一种基于极化敏感阵列的相干信源测向方法 | |
CN103116162B (zh) | 基于目标空间稀疏性的高分辨声呐定位方法 | |
CN104391183A (zh) | 近场测量天线远区场特性快速计算方法 | |
CN104950297A (zh) | 基于矩阵1范数拟合的阵元误差估计方法 | |
CN106970348B (zh) | 电磁矢量传感器阵列解相干二维music参数估计方法 | |
Liu et al. | Joint angle and range estimation for bistatic FDA-MIMO radar via real-valued subspace decomposition | |
CN110376547B (zh) | 基于二阶统计量的近场源定位方法 | |
CN109212466B (zh) | 一种基于量子蜻蜓演化机制的宽带测向方法 | |
Liang et al. | Two-dimensional DOA estimation method of acoustic vector sensor array based on sparse recovery | |
CN108614234B (zh) | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 | |
CN109633521A (zh) | 基于子空间重构的面阵二维波达方向估计方法 | |
Liu et al. | 2D DOA estimation for coherent signals with acoustic vector-sensor array | |
CN104678350A (zh) | 大规模mimo***中基于tls-esprtt算法的2d doa估计 | |
Cao et al. | Two-dimensional direction of arrival estimation using generalized ESPRIT algorithm with non-uniform L-shaped array | |
CN113821907B (zh) | 针对大型平面天线阵列***的幅度和相位自动校准方法 | |
Yang et al. | Coprime L‐shaped array connected by a triangular spatially‐spread electromagnetic‐vector‐sensor for two‐dimensional direction of arrival estimation |
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 | ||
GR01 | Patent grant |