CN103487784B - 一种基于信号到达时间的定位方法 - Google Patents
一种基于信号到达时间的定位方法 Download PDFInfo
- Publication number
- CN103487784B CN103487784B CN201310398801.3A CN201310398801A CN103487784B CN 103487784 B CN103487784 B CN 103487784B CN 201310398801 A CN201310398801 A CN 201310398801A CN 103487784 B CN103487784 B CN 103487784B
- Authority
- CN
- China
- Prior art keywords
- base station
- terminal
- coordinate system
- distance
- positioning
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000004807 localization Effects 0.000 title abstract 2
- 238000005315 distribution function Methods 0.000 claims abstract description 15
- 230000001413 cellular effect Effects 0.000 claims abstract description 8
- 238000005259 measurement Methods 0.000 abstract description 13
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/14—Determining absolute distances from a plurality of spaced points of known location
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0284—Relative positioning
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mobile Radio Communication Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于信号到达时间的定位方法,其将非线性的定位问题转换为多重积分问题,在非线性的定位问题转换为多重积分问题的过程中,根据蜂窝移动网络环境中的每个基站的位置、每个基站的信号的测量噪声方差以及待定位终端到每个基站的距离,构造一个高斯分布函数,然后从高斯分布函数中随机抽取若干个样本,再根据抽取到的样本获取待定位终端的位置,优点在于构造的高斯分布函数能提高计算多重积分的效率及积分值的精度,且在求解多重积分问题的过程中只需要抽取少数样本就能得到待定位终端的精确位置,计算复杂度低。
Description
技术领域
本发明涉及一种定位方法,尤其是涉及一种基于信号到达时间的定位方法。
背景技术
在雷达和声纳的研究中,终端定位是一个经典的研究课题。最近30年来,尤其在911事件之后,人们对终端定位的要求越来越高,从而使终端定位受到了越来越多的关注。终端定位技术在军事侦察、交通监视、工农业控制、生物医疗、环境监测、抢险救灾及危险区域远程控制等领域都有广阔的应用前景,因此研究定位方法具有十分重要的意义。图1给出了一个典型的定位问题的示意图,待定位终端通过测量不同基站的下行导频信号,得到不同基站的下行导频信号的TOA(TimeofArrival,到达时间)或TDOA(TimeDifferenceofArrival,到达时间差),根据该测量结果同时结合每个基站的坐标,就能够计算出待定位终端的位置。与GPS定位相比,蜂窝移动***中的定位以其定位速度快、成本低(不需要移动终端上添加额外的硬件)、耗电少、室内可用等优势,而被广泛应用到终端定位中。
尽管已经有学者提出了很多种方法解决定位问题,但是由于定位问题的非凸性,即定位问题具有多个稳定的平衡态,无法同时满足定位高精度化以及算法简单化的要求。在早期的研究中,通常采用基于泰勒展开的线性化方法来解决定位问题,其根据给定的初始估计将原定位问题进行泰勒展开并进行线性化近似,求解到线性解后,再根据此线性解进行泰勒展开,如此迭代直到达到最后的解,该方法要求有一个精确的初始估计,否则将无法收敛到全局最优解甚至得到的解发散。针对基于泰勒展开的线性化方法所存在的问题,人们提出了球形插值(SphericalInterpolation,SI)法,其通过引入一个中间变量来线性化非线性方程组,设定的中间变量不是想要的,但它与待定位终端的位置的关系已知,该方法基于线性最小二乘(LinearLeastSquares,LLS)准则求解线性方程组,但完全忽略了上面提到的中间变量与终端位置之间已知的关系,导致获得的待定位终端的位置不精确。为提高SI方法的性能,人们提出了一种两步加权最小二乘(Two-StageWeightedLeastSquares,Two-StageWLS)方法,该方法明确考虑了中间变量与待定位终端的位置之间的关系,因而能够取得比SI方法更好的性能,但是该方法有门限效应,即当测量时的噪声达到一定程度时,定位精度将急剧下降。近年来人们尝试采用一些非线性规划的理论来解决定位问题,由于定位问题的非凸性,人们尝试采用松弛方法将非凸问题松弛为凸问题(如半正定规划),求解该凸问题后,再采用一些后续处理方法来得到原定位问题的近似解,这种方法的特点是得到的待定位终端的位置的精度高,但其算法的复杂度远高于线性化方法。
发明内容
本发明所要解决的技术问题是提供一种定位精度高、计算复杂度低的基于信号到达时间的定位方法。
本发明解决上述技术问题所采用的技术方案为:一种基于信号到达时间的定位方法,其特征在于:包括以下步骤:
①首先在蜂窝移动网络环境中建立一个平面直角坐标系或空间直角坐标系作为参考坐标系,然后假设该蜂窝移动网络环境中存在一个待定位终端和N个基站,则该待定位终端通过接收N个基站的信号获取每个基站到待定位终端的到达时间,再根据每个基站到待定位终端的到达时间获取待定位终端到每个基站的距离,将待定位终端到第i个基站的距离记为zi,zi=c1ti,其中,c1为光速,ti为第i个基站到待定位终端的到达时间,i=1,2,…,N,N≥3;
②根据每个基站在参考坐标系中的位置以及待定位终端到每个基站的距离,构造一个高斯分布函数,记为其中,R=(BTQ-1B)-1,BT表示B的转置, Q=DQ1DT, DT表示D的转置, “||||”为欧几里得范数,为第1个基站的信号的测量噪声方差,为第2个基站的信号的测量噪声方差,为第N个基站的信号的测量噪声方差,s1为第1个基站在参考坐标系中的位置,s2为第2个基站在参考坐标系中的位置,sN为第N个基站在参考坐标系中的位置,z1为待定位终端到第1个基站的距离,z2为待定位终端到第2个基站的距离,zN为待定位终端到第N个基站的距离;
③在高斯分布函数中随机抽取M个样本,将随机抽取到的M个样本中的第m个样本记为xm,其中,m=1,2,…,M,M≥50;
④根据随机抽取到的M个样本,获取待定位终端在参考坐标系中的位置,记为 其中, n表示xm的维数,“||”表示行列式符号,exp表示自然基数,表示的转置,∫…∫表示n重积分符号,λ为一个常数,x1,…,xm表示积分变量,f(xm)=(z-dm)-1Q1 -1(z-dm),z=[z1,z2,…,zN]T,dm=[||xm-s1||,||xm-s2||,…,||xm-sN||]T,“||||”为欧几里得范数。
与现有技术相比,本发明的优点在于:将非线性的定位问题转换为多重积分问题,在非线性的定位问题转换为多重积分问题的过程中,首先构造一个高斯分布函数,然后在该高斯分布函数中随机抽取少量样本来解决该多重积分问题,从而获取到待定位终端的精确位置;高斯分布函数的构建,能提高计算多重积分的效率及积分值的精度,从而提高了获取到的待定位终端的位置的精度,且降低了计算的复杂度。
附图说明
图1为典型的定位环境示意图;
图2为本发明的总体流程图;
图3为本发明的待定位终端位于(1000,1000)时定位精度随噪声大小的变化图,其中,dBm2=10lg(m2),m2(平方米)是σ2的单位,σ表示基站的信号的测量噪声的标准差,σ2表示基站的信号的测量噪声的方差,MSE表示待定位终端的定位均方误差,用于衡量定位精度;
图4为本发明的待定位终端位于(5000,6000)时定位精度随噪声大小的变化图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明提出了一种基于信号到达时间的定位方法,图2给出了其总体流程示意图,其主要包括以下步骤:
①首先在蜂窝移动网络环境中建立一个平面直角坐标系或空间直角坐标系作为参考坐标系,假设该蜂窝移动网络环境中存在一个待定位终端和N个基站,则该待定位终端通过接收N个基站的信号获取每个基站到待定位终端的到达时间,再根据每个基站到待定位终端的到达时间获取待定位终端到每个基站的距离,将待定位终端到第i个基站的距离记为zi,zi=c1ti,其中,c1为光速,ti为第i个基站到待定位终端的到达时间,i=1,2,…,N,N≥3。
在此过程中,获取第i个基站信号到达待定位终端所需要的时间的过程如下:在某一时刻t,第i个基站在其时钟的控制下发出某一结构的测距码(该测距码为伪随机码,在待定位终端可复制),与此同时待定位终端则在待定位终端时钟的控制下复制出结构完全相同的测距码(以下简称复制码)。由第i个基站所产生的测距码经ti时间的传播后到达待定位终端并被待定位终端所接收。在待定位终端,由第i个基站所产生的复制码则经过一个时间延迟器延迟时间τ后与接收到的基站信号做相关运算,直到检测到相关运算输出的峰值为止,此时测量获得的延迟时间τ即为到达时间(TOA)ti。要得到准确的待定位终端到达基站所需要的时间,就需要基站与待定位终端的时钟保持高精度同步。但在实际中,由于基站与待定位终端时钟相对于***时钟的偏移,导致测得的到达时间有误差,这也是TOA测量噪声的主要来源。
②根据每个基站在参考坐标系中的位置以及待定位终端到每个基站的距离,构造一个高斯分布函数,记为其中,R=(BTQ-1B)-1,BT表示B的转置, Q=DQ1DT, DT表示D的转置, “||||”为欧几里得范数,为第1个基站的信号的测量噪声方差,为第2个基站的信号的测量噪声方差,为第N个基站的信号的测量噪声方差,s1为第1个基站在参考坐标系中的位置,s2为第2个基站在参考坐标系中的位置,sN为第N个基站在参考坐标系中的位置,z1为待定位终端到第1个基站的距离,z2为待定位终端到第2个基站的距离,zN为待定位终端到第N个基站的距离。本发明中,将定位问题转换成多重积分问题,而高斯分布函数的构造就是为了解决该多重积分问题。
③在高斯分布函数中随机抽取M个样本,将随机抽取到的M个样本中的第m个样本记为xm,其中,m=1,2,…,M,而考虑到定位的精度以及运算的复杂度,取M≥50。
④根据随机抽取到的M个样本,获取待定位终端在参考坐标系中的位置,记为 其中, n表示xm的维数,“||”表示行列式符号,exp表示自然基数,表示的转置,∫…∫表示n重积分符号,λ为一个常数,如取λ=1000,x1,…,xn表示积分变量,f(xm)=(z-dm)-1Q1 -1(z-dm),z=[z1,z2,...,zN]T,dm=[||xm-s1||,||xm-s2||,…,||xm-sN||]T,“||||”为欧几里得范数。
以下通过仿真来验证本发明的定位方法的可行性、有效性以及定位性能。假设安置5个基站来进行测量,测量的方法为:首先建立一个平面直角坐标系,5个基站在建立的平面直角坐标系上的位置依次为(0,0)、(0,6000)、 这5个基站形成一个测量区域。另外,在测量过程中,取M=100,且取λ=1000。克拉美-罗界(CRB)是目前最常用于评价参数估计性能的一个性能下界,故定位方法在衡量定位精度的时候,引入了克拉美-罗界,即任何无偏定位方法的定位性能均不可能低于克拉美-罗界。
图3给出了当待定位终端位于(1000,1000)时的定位精度,从中可以看出,当待定位终端位于测量区域中心时,在基站的信号的测量噪声方差由小至大的变化过程中,本发明提出的定位方法的定位性能都逼近于克拉美-罗界,因此具有很高的测量精度。
图4给出了当待定位终端位于(5000,6000)时的定位精度,从中可以看出,当待定位终端在测量区域边缘时,在整个基站的信号的测量噪声方差范围内,本发明公开的方法的定位性能仅在基站的信号的测量噪声方差较高的情况下才会略微高于克拉美-罗界,因此即使对于远场终端,在较大的噪声环境下本发明提出的定位方法仍然可以有较高的定位精度。
由仿真结果可以看出,本发明提出的定位方法具有良好的定位性能,能够很好地满足定位高精度的需求,且本发明提出的定位方法只需要通过构造一个高斯分布函数并在该高斯分布函数中取出少量样本就能得到待定位终端的精确位置,操作简单且计算复杂度低。
Claims (1)
1.一种基于信号到达时间的定位方法,其特征在于:包括以下步骤:
①首先在蜂窝移动网络环境中建立一个平面直角坐标系或空间直角坐标系作为参考坐标系,然后假设该蜂窝移动网络环境中存在一个待定位终端和N个基站,则该待定位终端通过接收N个基站的信号获取每个基站到待定位终端的到达时间,再根据每个基站到待定位终端的到达时间获取待定位终端到每个基站的距离,将待定位终端到第i个基站的距离记为zi,zi=c1ti,其中,c1为光速,ti为第i个基站到待定位终端的到达时间,i=1,2,…,N,N≥3;
②根据每个基站在参考坐标系中的位置以及待定位终端到每个基站的距离,构造一个高斯分布函数,记为其中,R=(BTQ-1B)-1,BT表示B的转置, Q=DQ1DT, DT表示D的转置, “||||”为欧几里得范数,为第1个基站的信号的测量噪声方差,为第2个基站的信号的测量噪声方差,为第N个基站的信号的测量噪声方差,s1为第1个基站在参考坐标系中的位置,s2为第2个基站在参考坐标系中的位置,sN为第N个基站在参考坐标系中的位置,z1为待定位终端到第1个基站的距离,z2为待定位终端到第2个基站的距离,zN为待定位终端到第N个基站的距离;
③在高斯分布函数中随机抽取M个样本,将随机抽取到的M个样本中的第m个样本记为xm,其中,m=1,2,…,M,M≥50;
④根据随机抽取到的M个样本,获取待定位终端在参考坐标系中的位置,记为 其中, n表示xm的维数,“||”表示行列式符号,exp表示自然基数,表示的转置, ∫…∫表示n重积分符号,λ为一个常数,x1,…,xn表示积分变量,f(xm)=(z-dm)-1Q1 -1(z-dm),z=[z1,z2,…,zN]T,dm=[||xm-s1||,||xm-s2||,…,||xm-sN||]T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310398801.3A CN103487784B (zh) | 2013-09-04 | 2013-09-04 | 一种基于信号到达时间的定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310398801.3A CN103487784B (zh) | 2013-09-04 | 2013-09-04 | 一种基于信号到达时间的定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103487784A CN103487784A (zh) | 2014-01-01 |
CN103487784B true CN103487784B (zh) | 2016-01-20 |
Family
ID=49828154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310398801.3A Expired - Fee Related CN103487784B (zh) | 2013-09-04 | 2013-09-04 | 一种基于信号到达时间的定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103487784B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106507472A (zh) * | 2015-09-07 | 2017-03-15 | 北京信威通信技术股份有限公司 | 一种无线通信***中的定位方法 |
CN110596688B (zh) * | 2019-08-29 | 2021-06-08 | 电子科技大学 | 一种分布式mimo雷达***中的高精度定位方法及装置 |
CN110855840B (zh) * | 2019-11-29 | 2022-02-11 | 合肥开元埃尔软件有限公司 | 一种基于网络推送和主动扫描相结合的手机安全切换*** |
CN111263321A (zh) * | 2019-12-16 | 2020-06-09 | 重庆邮电大学 | 一种提高toa室内超宽带定位精度的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW201042279A (en) * | 2009-05-19 | 2010-12-01 | Ralink Technology Corp | Method and apparatus of using soft information for enhancing accuracy of position location estimation for a wireless communication system |
CN102262234A (zh) * | 2010-05-27 | 2011-11-30 | 精工爱普生株式会社 | 位置计算方法及位置计算装置 |
CN102788967A (zh) * | 2011-05-16 | 2012-11-21 | Ls产电株式会社 | 使用到达时间差的实时定位*** |
CN103152824A (zh) * | 2013-03-06 | 2013-06-12 | 中国电子科技集团公司第二十八研究所 | 一种无线传感器网络中节点定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100434259B1 (ko) * | 1998-08-21 | 2004-07-16 | 엘지전자 주식회사 | 이동통신시스템에서의이동단말기위치탐지방법 |
US8711038B2 (en) * | 2006-10-05 | 2014-04-29 | Her Majesty The Queen In Right Of Canada As Represented By The Minister Of Industry, Through The Communications Research Centre Canada | High-resolution ranging and location finding using multicarrier signals |
-
2013
- 2013-09-04 CN CN201310398801.3A patent/CN103487784B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW201042279A (en) * | 2009-05-19 | 2010-12-01 | Ralink Technology Corp | Method and apparatus of using soft information for enhancing accuracy of position location estimation for a wireless communication system |
CN102262234A (zh) * | 2010-05-27 | 2011-11-30 | 精工爱普生株式会社 | 位置计算方法及位置计算装置 |
CN102788967A (zh) * | 2011-05-16 | 2012-11-21 | Ls产电株式会社 | 使用到达时间差的实时定位*** |
CN103152824A (zh) * | 2013-03-06 | 2013-06-12 | 中国电子科技集团公司第二十八研究所 | 一种无线传感器网络中节点定位方法 |
Non-Patent Citations (2)
Title |
---|
New Semidefinite Relaxation Method for Acoustic Energy-Based Source Localization;Gang Wang等;《IEEE SENSORS JOURNAL》;20130531;第13卷(第5期);1514-1521 * |
无线传感器网络中定位与跟踪算法的研究;王刚;《中国博士学位论文全文数据库(电子期刊)信息科技辑》;20130515(第5期);I136-46 * |
Also Published As
Publication number | Publication date |
---|---|
CN103487784A (zh) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lin et al. | A new constrained weighted least squares algorithm for TDOA-based localization | |
Wang et al. | TDOA source localization in the presence of synchronization clock bias and sensor position errors | |
CN108668358B (zh) | 一种应用于无线传感网络的基于到达时间的协作定位方法 | |
Wang et al. | A bias-reduced nonlinear WLS method for TDOA/FDOA-based source localization | |
Montillet et al. | Extracting white noise statistics in GPS coordinate time series | |
CN103487784B (zh) | 一种基于信号到达时间的定位方法 | |
CN101588328B (zh) | 一种高精度无线信道参数化模型的联合估计方法 | |
CN110673089B (zh) | 未知视距和非视距分布情况下基于到达时间的定位方法 | |
Yang | On adding and removing sensors in a solar irradiance monitoring network for areal forecasting and PV system performance evaluation | |
Qiao et al. | An improved method of moments estimator for TOA based localization | |
Arias-de-Reyna et al. | Crowd-based learning of spatial fields for the internet of things: From harvesting of data to inference | |
Forsythe et al. | Ionospheric horizontal correlation distances: Estimation, analysis, and implications for ionospheric data assimilation | |
Hua et al. | Geometry-based non-line-of-sight error mitigation and localization in wireless communications | |
Mahfouz et al. | Kernel-based localization using fingerprinting in wireless sensor networks | |
Rife et al. | Symmetric overbounding of correlated errors | |
CN111007457B (zh) | 一种基于块稀疏贝叶斯模型的辐射源直接定位方法 | |
Cao et al. | Multidimensional scaling‐based passive emitter localisation from time difference of arrival measurements with sensor position uncertainties | |
Ma et al. | Super-resolution time delay estimation using exponential kernel correlation in impulsive noise and multipath environments | |
Sun et al. | Computationally attractive and statistically efficient estimator for noise resilient TOA localization | |
CN106383333A (zh) | 一种基于互相关的改进时延估计方法 | |
CN110536410B (zh) | 非视距环境下基于rss和tdoa测量的定位方法 | |
Huang et al. | Research on regional ionospheric TEC modeling using RBF neural network | |
CN113923590B (zh) | 一种锚节点位置不确定情况下的toa定位方法 | |
Wang et al. | Efficient semidefinite relaxation for energy-based source localization in sensor networks | |
CN113347705B (zh) | 基于ica和cs的超宽带aoa室内定位方法 |
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: 20160120 Termination date: 20180904 |