摘要
提出了一种应用于均匀线阵中近场源定位的降秩Capon算法。该算法能够将经典二维Capon(Two‑dimensional Capon, 2D‑Capon)算法中的二维谱峰搜索转化为一维谱峰搜索,得到自动配对的近场信源角度和距离参数估计。与经典的2D‑Capon算法相比,本文提出的算法无需信源数估计,同时由于避免二维谱峰搜索过程,其计算复杂度大大降低,且该算法参数估计性能与2D‑Capon算法非常接近。仿真结果表明该算法可有效用于近场信源的参数估计。
空间信源定位是阵列信号处理领域中的一项关键技术,它在医学成像、雷达、无线通信、声呐等方面都有着广泛的应用
针对近场信源的定位问题,国内外学者做了大量研究工作,提出了多种应用于近场的信源定位方法,根据原理的不同,可大致分为非谱峰搜索和谱峰搜索两类。非谱峰搜索类的算法一般借助二阶或高阶统计量,通过计算闭式解得到信源的参数估计。近年来学者们提出了多种基于二阶统计量的算法
谱峰搜索类算法的共同特点是估计精度高,但计算量巨大。Swindlehurst等
从上述两类算法的介绍与分析中可知,如何有效地降低计算复杂度,避免谱峰搜索和进行参数配对,同时最大限度地提升参数估计精度,是近场信源定位技术研究的关键点。基于此,本文将矩阵降秩思想与Capon算法结合,对经典的近场源估计方法进行简化,提出了一种均匀线阵中基于降秩(Rank reduce, RARE)思想的近场源参数估计方法。本文的主要贡献如下:(1)提出了基于降秩思想的角度和距离参数联合估计方法;(2)相较于经典二维Capon(Two‑dimensional Capon, 2D‑Capon)算法,本文算法避免了二维谱峰搜索,大大减小了计算复杂度;(3)本文算法的参数估计性能接近经典2D‑Capon算法,具有较高的参数估计精度;(4)本文算法无需信源数估计。
如

图1 近场信源定位均匀线阵模型
Fig.1 Structure of uniform linear array for near‑field sources localization
则第
m个阵元上的接收信号可以表示为
(1) |
式中: ; ; 表示第 k个信源发出的信号被第 m个阵元接收并解调后的基带信号; 表示阵元上的加性噪声; 为第 k个信号入射方向的反方向与 y轴之间的夹角; 为第 k个信号的波长; 为信源到参考阵元之间的距离。将式(1)写成矩阵的形式为
(2) |
式中: 为接收信号矩阵, 为均匀线阵的方向矩阵, 为信源矩阵, 为阵列接收噪声矩阵。 中的列向量 为导向矢量,具有如下形式
(3) |
利用 J个快拍的接收信号,可以计算信号协方差矩阵为
(4) |
为不失一般性,本文做如下假设:
(1) 信源为相互统计独立的窄带随机过程,零均值,具有非零功率,信源的波长归一化为1;
(2) 阵元接收噪声为零均值、白或色高斯噪声,并与信源统计独立;
(3) 对于不同的信源,即 ,相位参数各不相同,即满足 ;
(4) 阵元为全向阵元且响应特性完全相同,无通道不一致、互耦等因素的影响,空间增益为1,阵元间距满足 ;
(5) 阵元个数与信源个数满足 。
在经典的近场2D‑Capon算法中,信源参数 可通过式(5)在空域中进行全局谱峰搜索得到
(5) |
式中 为
(6) |
经典2D‑Capon算法需要全局二维谱峰搜索,复杂度很高。为了降低算法复杂度,本文借鉴降秩思想,提出降秩Capon算法来实现二维参数估计,该算法有效避免了高复杂度的二维谱峰搜索过程。
由于阵列结构的对称性,导向矢量可以分解成如下形式
(7) |
式中
(8) |
式中: 仅包含信源的角度信息, 同时包含角度和距离信息。由式(8)可知 ,故可将式(7)代入式(5),得到
(9) |
式中: ,可见 中只包含信源角度参数信息。又 可知, 为非负定的共轭对称矩阵,因此 成立的充要条件为当且仅当 为奇异矩阵。由假设条件可知,当 时,噪声子空间 的列秩不小于 ,则可知 为满秩矩阵,只有当角度参数信息取到信源的实际位置时,矩阵 会降秩,即 ,此时 就会变成奇异矩阵,正交性成立。因此可以通过式(10)的一维谱峰搜索得到信源的DOA估计为
(10) |
式中: 表示取最大值时对应的角度值; 表示取行列式值; 。
由式(10)得到信源的角度估计参数之后,将 依次逐个代入经典2D‑Capon谱函数中,并构造式(11)中的距离搜索的谱函数,在距离上进行一维谱峰搜索,可得到距离参数的估计 为
(11) |
式中:距离的搜索范围 , ,由于需要将 K个角度估计逐个代入,可知需要进行 次一维搜索。上述搜索过程能使得距离估计 与角度估计 自动配对。
至此,已经完成了均匀线阵中近场信源基于降秩Capon算法的角度和距离参数的估计,该降秩Capon算法的主要步骤总结如下:
步骤1 根据式(4)计算接收信号协方差矩阵 ;
步骤2 根据式(5)构造谱峰搜索函数,并按式(7)将导向矢量 拆分为 ,并构造 ;
步骤3 利用 ,由式(10)构造关于角度信息的一维函数,通过角度搜索得到接收信号的DOA估计;
步骤4 将得到的信源DOA估计结果逐个代入式(11),然后再通过距离的一维谱峰搜索,得到与角度参数配对的距离估计。
降秩Capon算法的复杂度主要包括:计算接收信号的协方差矩阵
需要
,求
需要
,角度搜索需要
,
次距离搜索
,因此总的复杂度为
;传统的经典2D‑Capon算法的复杂度为
。其中,
为角度空间的谱峰搜索次数;
为近场距离区间内的谱峰搜索次数,
和
为搜索步长;
为阵元个数;
;
为快拍数;
为信源个数。

图2 算法复杂度对比图
Fig.2 Complexity comparison of two algorithms
本文所提算法优点总结如下:
(1)该算法能够有效实现近场源角度与距离参数的联合估计,且参数自动配对;
(2)该算法避免二维谱峰搜索,相比较于经典的2D‑Capon 算法,大大降低了计算的复杂度;
(3)该算法的参数估计性能非常接近经典2D‑Capon 算法,具有较高的参数估计精度;
(4)该算法无需信源数估计。
本文采用蒙特卡洛实验仿真,仿真中假设有两个近场信号被阵列所接收,其角度和距离参数分别为 和 。 M、 K、 J分别为阵列阵元数、信源数和接收信号快拍数。为了评估算法的参数估计性能,仿真实验次数为1 000次。角度和距离估计的求根均方误差(Root mean square error, RMSE)分别定义如下
(12) |
(13) |
式中: 和 分别为第 k个信源的角度和距离的实际值, 和 分别为第 i次实验中得到的第 k个信源的角度和距离参数的估计值。
仿真1

图3 角度和距离参数估计的散布图
Fig.3 Angle and range estimation of the proposed algorithm
仿真2

图4 角度的估计性能随快拍数变化情况
Fig.4 Angle estimation performance versus different snapshots

图5 距离的估计性能随快拍数变化情况
Fig.5 Range estimation performance versus different snapshots
仿真3

图6 角度估计性能对比
Fig.6 Angle estimation performance comparison

图7 距离估计性能对比
Fig.7 Range estimation performance comparison
针对均匀线阵中近场信源的角度和距离参数联合估计问题,本文提出了一种降秩Capon算法。该算法无需信源数估计,且由于不需要进行二维谱峰搜素,其计算复杂度远远低于传统的2D‑Capon算法。同时,该算法能够获得自动配对的角度和距离参数估计。仿真表明,其参数估计性能与经典的2D‑Capon算法非常接近,且具有较高的参数估计精度。
参考文献
张小飞, 沈金清, 汪云飞 . 电磁矢量互质阵中基于降维Capon的DOA和极化估计算法[J]. 数据采集与处理, 2018, 33( 6): 953‑ 961. [百度学术]
ZHANG Xiaofei, SHEN Jinqing, WANG Yunfei . DOA and polarization estimation for electromagnetic vector sensor coprime array via reduced‑dimension Capon[J]. Journal of Data Acquisition and Processing, 2018, 33( 6): 953‑ 961. [百度学术]
陈未央, 张小飞, 张立岑 . 声矢量传感器阵中基于Kalman滤波和OPASTd的DOA跟踪算法[J]. 南京航空航天大学学报, 2015, 47( 3): 377‑ 383. [百度学术]
CHEN Weiyang, ZHANG Xiaofei, ZHANG Licen . DOA tracking algorithm for acoustic vector‑sensor array via Kalman filter and OPASTd[J]. Journal of Nanjing University of Aeronautics and Astronautics, 2015, 47( 3): 377‑ 383. [百度学术]
徐乐, 张小飞, 林新平, 等 . 电磁矢量阵中基于PARALIND分解的相干DOA估计算法[J]. 数据采集与处理, 2019, 34( 4): 697‑ 705. [百度学术]
XU Le, ZHANG Xiaofei, LIN Xinping, et al . PARALIND decomposition‑based coherent direction of arrival estimation algorithm for electro‑magnetic vector array[J]. Journal of Data Acquisition and Processing, 2019, 34( 4): 697‑ 705. [百度学术]
CAO Renzheng, ZHANG Xiaofei . A generalized propagator algorithm for localization of non‑circular sources using arbitrary array geometry[J]. Transactions of Nanjing University of Aeronautics & Astronautics, 2019, 36( 2): 341‑ 351. [百度学术]
ZHAI Hui, CHEN Weiyang, ZHANG Xiaofei, et al . Low‑complexity DOA estimation of noncircular signals for coprime sensor arrays[J]. Transactions of Nanjing University of Aeronautics & Astronautics, 2019, 36( 4): 599‑ 608. [百度学术]
YE H, DEGROAT R D . Maximum likelihood DOA estimation and asymptotic Cramer‑Rao bounds for additive unknown colored noise[J]. IEEE Transactions on Signal Processing, 1995, 43( 4): 938‑ 949. [百度学术]
PESAVENTO M, GERSHMAN A B . Maximum‑likelihood direction‑of‑arrival estimation in the presence of unknown nonuniform noise[J]. IEEE Transactions on Signal Processing, 2001, 49( 7): 1310‑ 1324. [百度学术]
ROY R, KAILATH T . ESPRIT‑estimation of signal parameters via rotational invariance techniques[J]. IEEE Transactions on Acoustics Speech & Signal Processing, 1989, 37( 7): 984‑ 995. [百度学术]
SWINDLEHURST A L, OTTERSTEN B, ROY R, et al . Multiple invariance ESPRIT[J]. IEEE Transactions on Signal Processing, 1992, 40( 4): 867‑ 881. [百度学术]
MATHEWS C P, ZOLTOWSKI M D . Eigenstructure techniques for 2‑D angle estimation with uniform circular array[J]. IEEE Transactions on Signal Processing, 1994, 42( 9): 2395‑ 2407. [百度学术]
SCHMIDT R O . Multiple emitter location and signal parameter estimation[J]. IEEE Transactions on Antennas & Propagation, 1986, 34( 3): 276‑ 280. [百度学术]
ZOLTOWSKI M D, KAUTZ G M, SILVERSTEIN S D . Beamspace root‑MUSIC[J]. IEEE Transactions on Signal Processing, 1993, 41( 1): 344‑ 364. [百度学术]
RAO B D, HARI K V S . Weighted subspace methods and spatial smoothing: Analysis and comparison[J]. IEEE Transactions on Signal Processing, 1993, 41( 2): 788‑ 803. [百度学术]
ABED‑MERAIM K, HUA Y, BELOUCHRANI A . Second‑order near‑field source localization: Algorithm and performance analysis[C]//Proceedings of Conference on Signals, Systems & Computers. [S.l.]: [s. n. ], 1996( 1): 723‑ 727. [百度学术]
ABED‑MERAIM K, HUA Y . 3‑D near field source localization using second order statistics[C]//Proceedings of Asilomar Conference on IEEE,[S.l.]: IEEE. 1997( 2): 1307‑ 1311. [百度学术]
CHALLA R N, SHAMSUNDER S . High‑order subspace‑based algorithms for passive localization of near‑field sources[C]// Proceedings of Asilomar Conference on Signals. [S.l.]: [s.n.], 1995(2): 777. [百度学术]
HAARDT M, CHALLA R N, Shamsunder S . Improved bearing and range estimation via high‑order subspace based Unitary ESPRIT[C]//Proceedings of Conference on Signals, Systems & Computers. [S.l.]: [s. n. ], 1996( 1): 380‑ 384. [百度学术]
王波, 王树勋 . 一种基于二阶统计量的近场源三维参数估计方法[J]. 电子与信息学报, 2006, 28( 1): 45‑ 49. [百度学术]
WANG Bo, WANG Shuxun . A three‑dimensional parameter estimation method of near field sources based on second statistics[J]. Journal of Electronics & Information Technology, 2006, 28( 1): 45‑ 49. [百度学术]
梁军利, 杨树元, 王诗俊, 等 . 一种新的基于二阶统计量的近场源定位算法[J]. 电子与信息学报, 2008, 30( 3): 596‑ 599. [百度学术]
LIANG Junli, YANG Shuyuan, WANG Shijun, et al . A new near‑field source localization algorithm using second‑order statistics[J]. Journal of Electronics & Information Technology, 2008, 30( 3): 596‑ 599. [百度学术]
SWINDLEHURST A L, KAILATH T . Passive direction‑of‑arrival and range estimation for near‑field sources[C]//Proceedings of the Workshop on Spectrum Estimation & Modeling.[S.l.]: IEEE, 1988: 123‑ 128. [百度学术]
HUANG Y D, BARKAT M . Near‑field multiple source localization by passive sensor array[J]. IEEE Transactions on Antennas & Propagation, 1991, 39( 7): 968‑ 975. [百度学术]
WEISS A J, FRIEDLANDER B . Range and bearing estimation using polynomial rooting[J]. IEEE Journal of Oceanic Engineering, 1993, 18( 2): 130‑ 137. [百度学术]
STARER D, NEHORAI A . Passive localization of near‑field sources by path following[J]. IEEE Transactions on Signal Process, 2002, 42( 3): 677‑ 680. [百度学术]
GROSICKI E, ABED‑MERAIM K, HUA Y . A weighted linear prediction method for the near‑field source localization[J]. IEEE Transactions on Signal Processing, 2005, 53( 10): 3651‑ 3660. [百度学术]
LEE J H, LEE C M, LEE K K . A modified path‑following algorithm using a known algebraic path[J]. IEEE Transactions on Signal Processing, 1999, 47( 5): 1407‑ 1409. [百度学术]
LIANG J, LIU D . Passive localization of mixed near‑field and far‑field sources using two‑stage MUSIC algorithm[J]. IEEE Transactions on Signal Processing, 2009, 58( 1): 108‑ 120. [百度学术]
XIE J, TAO H, RAO X, et al . Passive localization of noncircular sources in the near‑field[C]//Proceedings of Radar Symposium. [S.l.]: IEEE, 2015: 493‑ 498. [百度学术]