CN102427602B - 稀疏直接定位方法 - Google Patents
稀疏直接定位方法 Download PDFInfo
- Publication number
- CN102427602B CN102427602B CN201110400367.9A CN201110400367A CN102427602B CN 102427602 B CN102427602 B CN 102427602B CN 201110400367 A CN201110400367 A CN 201110400367A CN 102427602 B CN102427602 B CN 102427602B
- Authority
- CN
- China
- Prior art keywords
- lattice point
- base station
- sparse
- position determination
- mobile terminal
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种稀疏直接定位方法。该方法由移动终端向参考基站发出定位请求,各参考基站分别接收移动终端发出的信号,并将信号汇集到定位中心;定位中心首先对定位区域进行格点划分,并将格点坐标组成集合;然后根据已知的基站和格点坐标预先建立过完备原子库,并建立稀疏定位模型;最后根据压缩感知原理对稀疏信号进行重构,所得稀疏向量中非零块所对应格点位置即为所求移动终端位置。本发明的方法由于不需要进行信号时延、角度、强度等特征参数估计,可以比两步定位方法获得更高的定位精度;同时该方法利用定位问题的稀疏性特征,依据接收信号判断稀疏向量中非零值所在位置,从而得出移动终端位置,可以大大降低计算复杂度。
Description
技术领域
本发明属于无线定位技术领域,特别涉及一种利用无线通信设施进行无线定位的方法。
背景技术
无线定位技术在搜索救援、智能交通、物流管理、地质勘探和国土开发、航海/航空导航等诸多领域都有广泛应用。如今在军事行动和日常生活的各个领域,具有空间位置特性的地理信息发挥着越来越重要的作用,高精度的无线定位技术也越来越受到广泛关注。
在众多无线定位***中,最著名的是把无线电发射源设置在各种轨道卫星上的定位***,例如美国的全球定位***(GPS)、欧洲的伽利略(Galileo)***、俄罗斯的GLONASS***以及我国的“北斗”定位***等,凭借着广域覆盖的巨大优势,将无线电定位技术发展到一个新的高度。尽管卫星定位技术已经在国民经济各个方面得到广泛应用,但是在应用领域由于受到各种接收误差的影响,需要通过其它辅助手段(例如建立差分基准站)才能达到所需的定位精度要求;同时在接收信号受到物理遮挡的情况下常常无法完成导航任务。因此,利用现有和即将建设的庞大的民用无线通信设施进行无线定位,不仅可以弥补卫星定位***的不足,而且可以作为无线通信高附加值的服务。尤其是在美国联邦通信委员会颁布了E911(Emergency call 911)强制性定位要求后,加上巨大市场利润的驱动,国内外出现了研究移动通信***终端定位技术的热潮。
受卫星定位***的影响,当前利用民用通信设施的定位方法绝大多数采用两步定位模式(如图2所示),即首先估计接收信号的一个或多个特征参数(如信号强度(RSS)、信号到达时间(TOA)、信号到达时间差(TDOA)和信号到达角度(AOA)等),再根据这些参数估计出移动终端的位置。此类定位模式的定位精度取决于第一步参数估计的精度和第二步位置解算算法的性能。在两步定位模式下,众多研究一直试图从提高参数估计精度和增强定位算法性能两个角度提出改进方案,以达到改善最终定位性能的目的。但不同于卫星定位***,无线通信网络不是专门为定位设计的***,其特征参数测量和估计是在完成通信功能的同时进行的,而且在信号格式上一般也没有为定位目的专门设计的定位信息域,难以保证特征参数测量精度,其误差不仅影响定位参数的估计精度,而且会传递到第二步位置解算阶段,引起更大误差。尽管已有大量关于参数估计的论文发表,但因为实际环境的复杂性和随机性,目前还没有在所有环境下都能进行高精度参数估计的一般方法。直接定位方法(Direct Position Determination,DPD)是无线定位领域近几年发展起来的一种定位新模式(如图3所示),该方法突破了传统两步定位方法的局限,将特征参数检测和定位融为一体,直接利用信号波形进行定位,无需专门的参数估计过程。
目前主流的直接定位方法是Weiss提出的格点搜索方法(Weiss AJ.“Directposition determination of narrowband radio frequency transmitters,”IEEE SignalProcessing Letters,2004,Vol.11,No.5,pp.513-516.)。假设有L个已知自身坐标的参考基站参与定位,该方法将第n个基站接收到待定位移动终端发出的信号表示为:
rn=αnθns(t-τn(p))+vn(t)
其中rn表示接收信号,s(t-τn(p))表示信号波形,p表示待定位移动终端的位置。αn表示移动终端与第n个基站之间的信道衰减。τn(p)表示电磁波从移动终端到基站所经历的时延,而θn=[θn(0),…,θn(M-1)]T是第n个基站的天线阵列响应,M表示天线阵元数目。vn(t)表示噪声矢量。
然后对rn抽样Ns点,并经离散傅里叶变换得:
其中T为采样间隔时间。在此基础上,依据最小二乘原则建立如下的代价函数:
通过求解该函数的最小值即可得到移动终端的位置估计。由于这是一个多维优化问题,计算复杂度很高。为了减少计算量,Weiss等人将定位区域划分成若干个格点,将每个格点作为移动终端可能的潜在位置,分别计算每个格点的代价函数值,其极值所对应位置即为所求位置。
很显然,格点划分得越密,定位精度越高,但计算量也越大。另外由于这种直接定位方法不像两步定位方法可以利用清晰的几何关系,其计算量要远高于传统两步定位方法,尤其不适用于实时性要求较高的场合。这也是当前直接定位方法没有得到广泛应用的重要原因。
发明内容
本发明的目的是针对现有技术中存在的不足,根据定位问题的稀疏性特征,依据压缩感知原理,提出一种基于稀疏恢复的直接定位方法,不仅去除了参数估计环节误差对定位精度的影响,而且从根本上解决现有直接定位方法计算量大的问题。
本发明为实现上述发明目的,所采用的技术方案是:
稀疏直接定位方法,包括如下步骤:
1)定位请求:假设待定位的移动终端位于p=(x,y),L个已知位置的基站位于qn=(xn,yn),n=1,…,L,每个基站均装有含M阵元的阵列天线;移动终端向参考基站发出定位请求,各参考基站分别接收移动终端发出的信号,并将信号汇集到定位中心;
2)格点划分:定位中心将基站所围成的定位区域划分为NΩ个格点,此处可以根据需要进行均匀划分或者非均匀划分,并将格点坐标组成集合由于在某一特定时刻移动终端所在的位置在空间域上是唯一的,也即移动终端只在其中的一个格点处(或其附近),于是,一旦定位区域内的格点位置划定,移动终端在空间上的位置可精确地用一个稀疏向量表示,仅仅移动终端位置所对应格点处的值非零,而其他格点对应索引值为0,如此一来,定位问题可以转变为依据接收信号判断稀疏向量中非零值所在位置的问题;
其中表示第m个阵元上的接收信号,表示信号波形,一般采用训练序列或同步序列等已知信号波形,αn表示移动终端与第n个基站之间的信道衰减,θn(m)是第n个基站上第m个阵元接收到信号的天线阵列响应,表示电磁波从格点到第n个基站所经历的时延;由于基站和格点位置都已知,时延可由两者间距离除以光速精确求得,同样也可以根据所选天线类型由基站位置与格点之间的几何关系求出到达角度值;
此处r(k)=[r(0)T,…,r(M-1)T]T是一个LM×1的矢量,其中第m项 是一个LM×L矩阵,而其中diag(·)表示对角矩阵;α=[α0,α1…,αL-1]T是一个L×1矢量,表示信道衰减系数;v(k)也是一个LM×1矢量,包括噪声和模型误差;
因此对于所有NΩ个格点,过完备原子库可以构造如下:
4)稀疏重构:根据压缩感知原理,本发明的稀疏直接定位方法可以归结求解下述方程:
(5)
本发明的有益效果主要体现在以下几个方面:
(1)本发明的方法由于不需要进行时延、角度、强度等特征参数估计,因而避免了同步、峰值检测等参数估计中的难点问题,同时也去除了参数估计环节误差对定位精度的影响,可以比两步定位方法获得更高的定位精度。
(2)本发明的方法利用定位问题的稀疏性特征,提出稀疏直接定位方法,将直接定位方法从现有的格点搜索方式变为依据接收信号判断稀疏向量中非零值所在位置的方式,一次计算即可得出移动终端位置坐标,可以大大降低计算复杂度。
(3)本发明的方法也不同于指纹(Fingerpring)定位法。本发明方法建立的过完备原子库由确定关系计算得到,无需像指纹数据库那样必须随着环境的改变而相应调整,甚至重建,可以有效地节省建设与维护成本。
附图说明
图1是本发明稀疏直接定位方法的流程图;
图2是现有技术中两步定位方式示意图;
图3是现有技术中直接定位方式示意图;
图4是本发明实施例中格点划分示意图;
图5是本发明实施例中仿真结果图。
具体实施例
为了更好地理解本发明的技术方案,以下将结合附图及具体实施例对本发明的工作流程及有益效果进行详细说明。
在本实施例中,定位区域如图4所示,4个基站的坐标分别为(500m,500m),(500m,-500m),(-500m,-500m),(-500m,500m),用符号n∈{1,2,3,4}表示。选择基站1作为定位中心,并且基站之间可以互相通信。每个基站安装一部均匀线阵天线(UniformLinearArray,ULA),天线阵元数均为10,阵元间隔取为信号波长的一半。信道幅度响应服从均值为0,标准差为0.01的正态分布,信道相位响应服从(-π,π)上的均匀分布,信号载波频率为900MHz。格点划分方式采用均匀划分方式,格点间隔为20m,这样在所围区域内共设定了NΩ=51×51个格点。由于定位区域是四个基站所围区域,因此这些格点位置可以预先得到,记为i∈{1,2,…,NΩ}。
由于基站和格点位置都是已知的,因此从格点到第n个基站所经历的时延可以由公式 直接计算得到,其中c表示光速。同样地,格点和第n个基站的波达角也可以通过计算得到,为于是,当采用均匀线阵时,第n个基站的阵列响应矢量为 其中ζ=2π/λ为波数,λ是信号波长,d为阵元间隔。注意此处的时延和阵列响应是根据几何关系直接计算得到,并非通过参数估计得到,因此可以在定位前预先获得。相应地,根据上述时延和阵列响应建立的过完备原子库也可以预先建立,并且只要基站位置和格点划分方法不改变,过完备原子库也不会发生改变。
在过完备原子库建立后,根据的块稀疏特性,本实施例采用一种块结构的贪婪算法-BCOSAMP算法(Baraniuk,R.G.,Cevher,V.,Duarte M.F.,and Hegde,C.,‘Model-based compressive sensing’,IEEE Transactions on Information Theory,2010,Vol.56,No.4,pp.1982-2001.)进行稀疏重构,所得中非零块所对应格点位置即为所求移动终端位置。一般来说,移动终端位置p不会正好等于格点位置但当格点密度足够大时,可以做到
采用计算机仿真验证本发明的性能(仿真中使用的是奔腾IV(2.4GHz)处理器和1GB内存的计算机),统计运行1000次的平均定位效果,仿真结果如图5所示。如图所示,稀疏直接定位方法的定位性能要优于Weiss方法,尤其在低信噪比时,该方法可以明显地改善定位精度;而在高信噪比时两种定位算法的性能接近。然而,稀疏直接定位方法平均运行一次只需0.021秒,而Weiss方法要7.251秒,可见稀疏直接定位方法的计算时间都要远小于Weiss方法,也即该方法的复杂度要低很多。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (1)
1.稀疏直接定位方法,其特征在于,所述方法包括如下步骤:
1)移动终端向参考基站发出定位请求,各参考基站分别接收移动终端发出的信号,并将信号汇集到定位中心;
3)利用基站和格点坐标,依据确定的几何关系建立过完备原子库,并建立稀疏定位模型;
其中表示第m个阵元上的接收信号,表示信号波形,αn表示移动终端与第n个基站之间的信道衰减,θn(m)是第n个基站上第m个阵元接收到信号的天线阵列响应,表示电磁波从格点到第n个基站所经历的时延L为参考基站的个数,每个基站均装有含M阵元的阵列天线,表示第k个采样时刻第n个基站上第m个阵元受到的噪声干扰以及模型误差;
此处r(k)=[r(0)T,…,r(M-1)T]T是一个LM×1的矢量,其中第m项 是一个LM×L矩阵,表示过完备原子库在格点处的分量,而其中定义diag(·)表示对角矩阵;α=[α0,α1…,αL-1]T是一个L×1矢量,表示信道衰减系数;v(k)是一个LM×1矢量,包括噪声和模型误差;
因此对于所有NΩ个格点,过完备原子库可以建立如下:
从而稀疏直接定位模型可以表示为:
4)进行稀疏重构;
利用压缩感知理论中的算法求解下述方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110400367.9A CN102427602B (zh) | 2011-12-06 | 2011-12-06 | 稀疏直接定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110400367.9A CN102427602B (zh) | 2011-12-06 | 2011-12-06 | 稀疏直接定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102427602A CN102427602A (zh) | 2012-04-25 |
CN102427602B true CN102427602B (zh) | 2014-05-14 |
Family
ID=45961534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110400367.9A Expired - Fee Related CN102427602B (zh) | 2011-12-06 | 2011-12-06 | 稀疏直接定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102427602B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102736056B (zh) * | 2012-06-18 | 2014-06-11 | 华为技术有限公司 | 一种目标定位方法、装置与*** |
CN102752851B (zh) * | 2012-06-29 | 2015-12-09 | 中国科学院深圳先进技术研究院 | 室内定位指纹库的指纹信息收集方法及*** |
CN103096465B (zh) * | 2013-01-11 | 2015-06-24 | 南京师范大学 | 一种环境自适应的多目标直接定位方法 |
CN104902491B (zh) * | 2014-03-06 | 2018-11-13 | 精英电脑(苏州工业园区)有限公司 | 决定无线存取点位置方法、无线存取点及无线存取点组合 |
CN104093202B (zh) * | 2014-07-02 | 2017-11-17 | 南京信息工程大学 | 一种环境自适应的无设备目标定位方法 |
CN106059971A (zh) * | 2016-07-07 | 2016-10-26 | 西北工业大学 | 一种信号相关性衰减情况下基于稀疏重构的相关检测方法 |
CN110095751A (zh) * | 2019-03-22 | 2019-08-06 | 中山大学 | 基于相关向量机实现数据驱动建模的目标定位跟踪*** |
CN113721194B (zh) * | 2021-07-30 | 2024-04-09 | 南京师范大学 | 一种基于mwcs的近场语音信号三维定位方法 |
CN115825863B (zh) * | 2022-12-16 | 2023-12-29 | 南京航空航天大学 | 一种冲击噪声下的非圆信号快速直接定位方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644774A (zh) * | 2009-07-07 | 2010-02-10 | 中国科学院声学研究所 | 一种基于线性规划的近场目标定位方法及其*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7250907B2 (en) * | 2003-06-30 | 2007-07-31 | Microsoft Corporation | System and methods for determining the location dynamics of a portable computing device |
-
2011
- 2011-12-06 CN CN201110400367.9A patent/CN102427602B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644774A (zh) * | 2009-07-07 | 2010-02-10 | 中国科学院声学研究所 | 一种基于线性规划的近场目标定位方法及其*** |
Non-Patent Citations (2)
Title |
---|
《基于非视距误差值直接估计的定位跟踪算法》;柯炜,吴乐南;《信号处理》;20101231;全文 * |
柯炜,吴乐南.《基于非视距误差值直接估计的定位跟踪算法》.《信号处理》.2010, |
Also Published As
Publication number | Publication date |
---|---|
CN102427602A (zh) | 2012-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102427602B (zh) | 稀疏直接定位方法 | |
Woo et al. | Application of WiFi-based indoor positioning system for labor tracking at construction sites: A case study in Guangzhou MTR | |
CN101620271B (zh) | 运用广播调频信号定位地理位置的方法及*** | |
CN102088769B (zh) | 直接估计和消除非视距误差的无线定位方法 | |
CN103501538B (zh) | 基于多径能量指纹的室内定位方法 | |
CN102291817B (zh) | 移动通信网络中的基于位置测量样本的群定位方法 | |
CN103995250B (zh) | 射频标签轨迹追踪方法 | |
US20110287778A1 (en) | Method and system for estimation of mobile station velocity in a cellular system based on geographical data | |
CN102981144A (zh) | 空中运动平台对目标的三维无源定位方法 | |
CN102100058A (zh) | 通过选择最佳wlan-ps方案使用混合卫星和wlan定位***确定定位的方法和*** | |
CN103402258A (zh) | 一种基于Wi-Fi的室内定位***和方法 | |
CN102460202A (zh) | 基于包括多个天线的无线通信装置的定位***和方法 | |
CN102711239A (zh) | 基于rss指纹数据库的二次模糊聚类室内定位方法 | |
CN103796304B (zh) | 一种基于虚拟训练集和马尔科夫链的煤矿井下定位方法 | |
CN109084774A (zh) | 一种基于信道状态信息的无人机集群定位方法 | |
CN102749637A (zh) | 一种车载gps精确定位的实现方法 | |
CN106231552B (zh) | 一种基于混合方式的非合作目标定位方法 | |
CN104038901A (zh) | 一种减少指纹数据采集工作量的室内定位方法 | |
CN103116159A (zh) | 一种多模式自定位组网雷达生命探测方法和装置 | |
CN102288938B (zh) | 一种有效的无线传感器网络节点的三维定位方法 | |
CN110926461A (zh) | 一种基于超宽带室内定位方法和***、导航方法和*** | |
CN102928860A (zh) | 基于局部定位信息提高gps定位精度的方法 | |
CN103096465B (zh) | 一种环境自适应的多目标直接定位方法 | |
CN110954865A (zh) | 一种基于电离层信息的短波时差定位方法 | |
CN103249144A (zh) | 一种基于c型的无线传感器网络节点定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20140514 Termination date: 20161206 |