饶中初1，2，郭贤生2
（1.同方电子科技有限公司，江西九江332007；2.电子科技大学电子工程学院，四川成都610054）
摘要： 为了解决复杂室内环境下的近场源定位问题，利用室内小型天线阵列，提出一种近场多径模型下的协方差矩阵匹配及信号子空间匹配定位方法，较好地解决了传统近场源定位算法中无法解决的多径传播问题。一方面，该方法改善了室内定位环境中由于模型失配引起的定位性能恶化问题；另一方面，它也为近场多径环境下的目标定位提供了一种新的解决思路。仿真结果表明，该算法不仅在定位精度上有了较大改善，而且在定位成功概率方面具有较大优势。
教育期刊网 http://www.jyqkw.com
关键词 ： 室内定位；空间谱；子空间匹配；二阶统计量
中图分类号：TN911.72?34 文献标识码：A 文章编号：1004?373X（2015）15?0018?04
收稿日期：2015?01?28
0 引言
辐射源的波达方向（DOA）估计在电子侦察、智能天线、雷达和声纳等领域有着广泛的应用，已受到人们的高度重视。大量的高分辨阵列测向技术都是运用于室外环境，例如MUSIC算法[1]和ESPRIT算法[2]。然而远场定位技术在近场源场景并不适用，尤其是存在大量多径的室内环境。
指纹库定位技术是一种常见的算法，文献[3?4]将该方法运用于定位远场窄带信号源。由于RSS值在静止环境中都会发生变化，这导致RSS算法的定位精度非常低。文献[5?6]提出了建立信号子空间指纹库提高室内定位精度，但是该算法仅仅考虑了平面波（远场窄带环境），然而在室内环境中，到达天线阵列的波形是球面波而不是平面波，在这种情况下，近场窄带模型比远场窄带模型更有效。
在复杂室内环境中，信道特征表征为复杂的多径传播，信号到达接收阵列为单个信号的不同延时复本，表现为相干多径情况。同时，对有些频段的信号源，其所表现出来的模型特征不再适合于远场窄带假设，如声信号等。
近场信号模型是一种较为有效地描述这种复杂环境信号传播特征的模型，同时考虑到多径及信道的非平稳性等特征[7?8]，其环境背景噪声可建模为非高斯随机过程，如满足冲击噪声等，这时候，传统的基于二阶统计量分析的经典信号处理方法将会带来性能的严重下降，因此，有必要发展基于分数低阶统计量的室内近场源定位理论及高精度定位方法。
本文将对近场源模型下的子空间匹配定位方法进行研究，为室内复杂环境下的近场源定位提供一种可行性解决方案。室内近场多径信号传播模型如图1所示。
1 均匀圆阵下近场源多径信号模型
考虑N 根天线的均匀圆阵，圆阵半径为R 。一个窄带信号包含q 条多径到达阵列天线。第l 个阵元的接收信号为：
式中：xkl (t)代表接收信号的第k 次反射到达第l个阵元在第t 次快拍时的数据；γk (t)是复系数，代表了第k次反射的衰减和相移；随机过程nl (t)是第l 个阵元接收到的高斯噪声，该噪声是零均值、未知方差σ2；L 代表快拍总数。
对于均匀圆阵，近场源的方向矢量a(r ) k ,θk ,?k 可表示为：
式中：ω = 2π λ，λ 是信号波长；rk ,θk∈[0, π 2],?k∈ [0,2π]分别表示信号源离基站的距离，到达阵列的俯仰角和到达阵列的方位角；(?)T 代表共轭转置。
第k 个信号和第l 个阵元的距离是dl (rk ,θk ,?k) ：
式中：ρl (θk ,?k) = sin θk cos(?k - (l - 1)θ0 )，l = 1,2,…,N,θ0 =2π/N 。使用泰勒级数展开，dl (rk ,θk ,?k) 可表示为：
如果信号源处于菲涅尔区域[9]，即信号源到接收阵列之间的距离满足：
则方向矢量可表示为[10]：
接收信号可表示为：
2 最大似然估计
假设接收信号x(t) 采样L 次，噪声是独立同分布的。高斯随机变量具有零均值、未知方差σ2，因此概率密度函数可以表达为[11]：
最大似然估计表达式为：
经过求导和消除常数因子后，上面求最大值的问题等价于如下求解最小值的问题：
上式的最小值问题，可以通过求解最小二乘解得到s? 的估计：
式中: (?)H 代表赫米特共轭转置。将s? 代入[A ?,s?]，可以得到：
式中：PA = A(AH A)-1 AH 是矩阵A 的投影矩阵。这样上式可以写成：
式中：tr{?} 是求秩运算；R ? 是协方差矩阵估计：
为了定位信号源，搜索指纹库找到与估计R ? 最相似的向量对应的坐标位置：
式中：R ? 代表实测阶段获得的协方差矩阵；Pi 是建库阶段第i 个格点得到的投影矩阵。
3 基于二阶统计量和子空间匹配的定位方法
3.1 信号子空间匹配定位方法
信号子空间定位算法的主要思想是：从协方差矩阵的特征值分解中得到信号子空间，然后建立信号子空间的指纹库。建库阶段，按格点储存最大特征值对应的特征向量S1；实测阶段，按照相同的策略得到实测信号子空间S2。下面将仔细说明信号子空间匹配技术。
式中：αk 代表S1 和S2 之间第k 个主角，如果S1 = S2，则αk = 0。对比指纹库中的每一个集合，挑出所得值最大的那个集合对应的坐标为发射源的位置。
3.2 协方差矩阵匹配定位方法
协方差匹配算法的主要思想如下：
（1）在建库阶段，计算k 个格点的协方差矩阵并用其建立指纹库：
（2）在实测阶段，接收数据为xT，则协方差矩阵为：
4 仿真分析
仿真在一个20 m×20 m×3 m的区域进行，参考位置的索引如图2所示，发射终端放置在地板上，每个格点之间的距离为1 m。
均匀圆阵半径R = 2λ, 该圆阵配备了12根天线。考虑加性高斯白噪声环境，仿真中包含了1个直达信号和3 条多径信号，3 条多径的复衰减系数分别是(0.4 + 0.8j), (-0.3 - 0.7j), (0.5 - 0.6j) 。500 次蒙特卡洛实验运用于不同信噪比下测试算法的性能。信噪比的变化为-10~8 dB，实验快拍数为500，测试点的位置在区域中随机选取。图3 给出了传统基于接收信号强度（RSS）的定位方法、基于协方差矩阵匹配的定位方法、最大似然定位方法以及基于信号子空间匹配的定位方法的定位误差曲线。
从图3可以看出，基于协方差矩阵匹配和子空间匹配的定位方法比传统的RSS 方法具有更高的定位性能。这两种方法定位的均方根误差接近于最大似然估计的定位结果。同时，细节性分析可以发现，子空间匹配方法的定位性能要优于协方差矩阵匹配定位方法，这是因为和协方差矩阵比较而言，协方差矩阵进行特征分解获取的信号子空间存储了信号的主要信息，能更有效地对抗多径传播及环境变化引起的信道时变特性。另外，由于传统RSS方法只能提供接收信号强度信息，而接收信号强度只是协方差矩阵的对角线元素的值，对协方差矩阵的非对角线元素，基于RSS的定位直接舍去这些携带大量信道环境信息的数据，其定位精度较完全利用其他自相关值的自相关矩阵匹配方法大大降低。从图3还可以看出，随着信噪比的增加，本文提出的协方差矩阵方法和子空间匹配定位方法的定位性能几乎和最大似然估计方法的定位结果重合，这也说明了本文提出方法的有效性。
为了进一步验证本文所提出方法的定位性能，图4给出了不同信噪比下几种方法的定位成功概率比较结果。值得说明的是，图中成功概率的定义为：
从图4可以看出，本文所提方法在SNR较大时，其成功概率较高，当信噪比为0 dB时，定位成功概率大于98%。
5 结论
本文主要介绍了多径模型下的近场源定位问题，利用子空间匹配和协方差匹配方法解决了传统近场源无法解决多径传播的情况。仿真结果表明，本文所提算法不仅在定位精度上有了较大改善，还在定位成功概率方面具有较大优势。
教育期刊网 http://www.jyqkw.com
参考文献
[1] SCHMIDT R O. Multiple emitter location and signal parameter estimation [J]. IEEE Transactions on Antennas and Propaga?tion，1986，34（3）：276?280.
[2] ROY R，KAILATH T. Esprit?estimation of signal parameters via rotational invariance techniques [J]. IEEE Transactions on ASSP，1989，37（7）：984?995.
[3] HATAMI A，PAHLAVAN K. In?building intruder detection for WLAN access [C]// 2004 Position Location and Navigation Sym?posium. [S.l.]：IEEE，2004：592?597.
[4] WEISS A J. On the accuracy of a cellular location system based on RSS measurements [J]. IEEE Transactions on Vehicu?lar Technology，2003，52（6）：1508?1518.
[5] TSUJI H，KIKUCHI S，KAVEH M. Indoor localization using subspace matching：an experimental evaluation [C]// Procee ?dings of 2006 IEEE Workshop on Sensor Array and Multichan?nel. [S.l.]：IEEE，2006：541?545.
[6] NEZAFAT M，KAVEH M，TSUJI H. Indoor localization using a spatial channel signature database [J]. IEEE Antennas and Wireless Propagation Letters，2006，5（1）：406?409.
[7] ZHI Wanjun，CHIA M Y. Near ? field source localization via symmetric subarrays [J]. IEEE Signal Processing Letters，2007，14（6）：409?412.
[8] HE Jin，SWAMY M N S，OMAIR AHMAD M. Efficient appli?cation of MUSIC algorithm under the coexistence of far?field and near?field sources [J]. IEEE Transactions on Signal Proces?sing，2012，60（4）：2066?2070.
[9] JOHNSON R C，JASIK H. Antenna engineering handbook [M].3rd ed. New York：McGraw?Hill，1993.
[10] HUANG Y D，BARKAT M. Near?field multiple source loca?lization by passive sensor array [J]. IEEE Transactions on An?tennas Propagation，1991，39（7）：968?975.
[11] ZISKIND I，WAX M. Maximum likelihood localization of mul?tiple sources by alternating projection [J]. IEEE Transactions on Acoustics，Speech and Signal Processing，1988，36（10）：1553?1560.
作者简介：饶中初（1969—），男，硕士，高级工程师，2013—2014年在电子科技大学从事国内访问学者。研究方向为无线通信、技术侦察和电子对抗技术等。
郭贤生（1979—），男，博士，副教授。研究方向为无线定位、阵列信号处理技术等。