CN108668358B - 一种应用于无线传感网络的基于到达时间的协作定位方法 - Google Patents
一种应用于无线传感网络的基于到达时间的协作定位方法 Download PDFInfo
- Publication number
- CN108668358B CN108668358B CN201810436213.7A CN201810436213A CN108668358B CN 108668358 B CN108668358 B CN 108668358B CN 201810436213 A CN201810436213 A CN 201810436213A CN 108668358 B CN108668358 B CN 108668358B
- Authority
- CN
- China
- Prior art keywords
- unknown target
- target source
- coordinate system
- reference coordinate
- ith
- 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
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000005259 measurement Methods 0.000 claims abstract description 54
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 11
- 230000005540 biological transmission Effects 0.000 claims description 35
- 238000005457 optimization Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000005516 engineering process Methods 0.000 abstract description 14
- 238000004891 communication Methods 0.000 description 6
- 230000001186 cumulative effect Effects 0.000 description 4
- 238000005315 distribution function Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 208000032370 Secondary transmission Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/006—Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
-
- 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/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W84/00—Network topologies
- H04W84/18—Self-organising networks, e.g. ad-hoc networks or sensor networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种应用于无线传感网络的基于到达时间的协作定位方法,其利用二阶锥松弛技术和半正定松弛技术对通过构造最大似然函数得到的最小化问题进行松弛,得到混合半正定和二阶锥规划的问题的描述,这样可以确保得到全局最优解而不受局部收敛的影响,定位精度高;并且能够有效地抑制测量噪声误差的影响;同时可利用现有技术进行求解出未知目标源的坐标的估计值,降低了锚节点的分布密度,减少了花费成本。
Description
技术领域
本发明涉及一种目标定位方法,尤其是涉及一种应用于无线传感网络的基于到达时间的协作定位方法。
背景技术
无线传感网络技术已广泛应用于军事、工业和商业活动中,其中确定无线传感网络中的每个传感器的具体位置是该技术要解决的核心问题之一。GPS定位仅在室外环境下有效,在室内等特殊场景,必须寻求新的定位方法。通常的方法是,先通过人工或其它手段部署位置已知的节点,这些位置已知的节点称为锚节点(Anchor-Node),而预先不知道自身位置的节点则通过锚节点来定位,这些位置未知的节点称为未知节点(Unknown-Node),即为未知目标源,再通常利用带有噪声的测量值来确定大量未知节点的位置,该方法也称为目标定位方法。
目前,目标定位方法有基于到达时间差(TDoA)的目标定位方法、基于接收信号强度(RSS)的目标定位方法、基于到达角度(AoA)的目标定位方法、基于到达时间(ToA)的目标定位方法等。与前三种方法相比,基于到达时间的目标定位方法只需要未知节点和锚节点的传输时间就可以得到未知节点与锚节点之间的距离,定位精度较高,从而实现精确的目标定位。
目前,基于到达时间的目标定位方法主要分为非协作定位和协作定位两种。由于非协作定位技术只允许锚节点和未知节点之间通信,因此当通信距离较大时,这种定位技术容易出现随时间断通信的问题,不能保证通信的连续性;另外,在非协作定位技术中,每个未知节点需要连接三个及以上锚节点才能完成未知节点定位,这就要求锚节点的分布密度高,导致成本花费大。为了避免非协作定位技术的弊端,并提高定位精度,一个重要的解决方法就是协作定位技术,与非协作定位技术相比,协作定位技术不但有锚节点和未知节点之间的通信连接,而且还有未知节点之间相互通信连接,因此,协作定位技术可以降低对锚节点的分布密度高的需求,并且能大大提高估计精度和鲁棒性。基于到达时间的协作定位方法最常用的估计方式是最大似然(ML)估计,其可以渐近地提供最优解。然而,求解最大似然估计是一个非常困难的工作,这是因为最大似然估计函数是非线性和非凸性的,并没有闭式解,可是当可用的锚节点数量有限或未知节点位于锚节点的凸包外时,虽然凸优化的复杂度有点增加,但是可以提供相当高的精度。因此,非常有必要研究一种能够求得全局最优解、保证定位精度的基于到达时间的协作定位方法。
发明内容
本发明所要解决的技术问题是提供一种应用于无线传感网络的基于到达时间的协作定位方法,其能够求得全局最优解、保证定位的精度。
本发明解决上述技术问题所采用的技术方案为:一种应用于无线传感网络的基于到达时间的协作定位方法,其特征在于包括以下步骤:
步骤一:在无线传感网络环境中,建立一个平面直角坐标系或空间坐标系作为参考坐标系,并设定存在M个锚节点和N个未知目标源;然后将M个锚节点在参考坐标系中的坐标对应记为s1,s2,…,sM,将N个未知目标源在参考坐标系中的坐标对应记为x1,x2,…,xN;其中,M≥2,N≥2,s1为第1个锚节点在参考坐标系中的坐标,s2为第2个锚节点在参考坐标系中的坐标,sM为第M个锚节点在参考坐标系中的坐标,x1为第1个未知目标源在参考坐标系中的坐标,x2为第2个未知目标源在参考坐标系中的坐标,xN为第N个未知目标源在参考坐标系中的坐标;
步骤二:在无线传感网络环境中,由每个未知目标源发射测量信号,每个未知目标源发射的测量信号经过无线传感网络环境传播后由所有锚节点和其它未知目标源接收;先确定每个未知目标源发射的测量信号从发射到各个锚节点接收所经历的时间,并确定每个未知目标源发射的测量信号从发射到其它各个未知目标源接收所经历的时间,将第i个未知目标源发射的测量信号从发射到第j个锚节点接收所经历的时间记为ti,j,将第i个未知目标源发射的测量信号从发射到第k个未知目标源接收所经历的时间记为ti,k;然后计算每个未知目标源发射的测量信号从发射到各个锚节点接收的传输距离测量值,并计算每个未知目标源发射的测量信号从发射到其它各个未知目标源接收的传输距离测量值,将第i个未知目标源发射的测量信号从发射到第j个锚节点接收的传输距离测量值记为di,j,di,j=c×ti,j,将第i个未知目标源发射的测量信号从发射到第k个未知目标源接收的传输距离测量值记为di,k,di,k=c×ti,k;其中,i、j、k均为正整数,1≤i≤N,1≤j≤M,1≤k≤N,k≠i,ti,j和ti,k的单位均为秒,di,j和di,k的单位均为米,c表示光速;
步骤三:将协作定位模型描述为:然后将协作定位模型中的di,j=||xi-sj||+ni,j中的ni,j移到等式左边后对等式两边平方,同样将di,k=||xi-xk||+ni,k中的ni,k移到等式左边后对等式两边平方,得到:接着忽略中的(ni,j)2和(ni,k)2,得到:再将转换为:其中,符号“|| ||”为求欧氏距离符号,||xi-sj||表示第i个未知目标源与第j个锚节点之间的真实距离,||xi-xk||表示第i个未知目标源与第k个未知目标源之间的真实距离,xi为第i个未知目标源在参考坐标系中的坐标,xk为第k个未知目标源在参考坐标系中的坐标,sj为第j个锚节点在参考坐标系中的坐标,ni,j表示第i个未知目标源与第j个锚节点之间的测量噪声,ni,k表示第i个未知目标源与第k个未知目标源之间的测量噪声;
步骤四:根据构造一个最大似然函数,得到最小化问题,描述为:然后令X=[x1,x2,…,xN],其中,min()为取最小值函数,σi,j表示ni,j的标准差,σi,k表示ni,k的标准差,yi,j为引入的变量,为xi的转置,gi,k为引入的变量,为xk的转置,X为引入的变量,符号“[]”为向量或矩阵表示符号,F为引入的变量,I2为二阶的单位矩阵,[X I2]T为[X I2]的转置,XT为X的转置;
步骤五:将最小化问题转化成一个优化估计问题,描述为:约束条件为: gi,k=Fi,i-Fi,k-Fk,i+Fk,k、然后利用二阶锥松弛技术对约束条件和进行松弛,将松弛为||[2((di,j)2-(yi,j)2)4(di,j)2×(σi,j)2-ui,j]||≤4(di,j)2×(σi,j)2+ui,j,将松弛为||[2((di,k)2-(gi,k)2)4(di,k)2(σi,k)2-hi,k]||≤4(di,k)2(σi,k)2+hi,k;接着利用半正定松弛技术对约束条件进行松弛,将改写为线性矩阵不等式形式:其中,ui,j为引入的优化变量,hi,k为引入的优化变量,Fi,i表示F中的第i行第i列元素,为sj的转置,FN+1:N+2,i表示由F中的第N+1行至第N+2行第i列的元素组成的维数为2×1的向量,Fi,k表示F中的第i行第k列元素,Fk,i表示F中的第k行第i列元素,Fk,k表示F中的第k行第k列元素,rank(F)表示求F的秩;
步骤六:确定最终的凸优化估计问题,描述为:约束条件为:||[2((di,j)2-(yi,j)2)4(di,j)2×(σi,j)2-ui,j]||≤4(di,j)2×(σi,j)2+ui,j、||[2((di,k)2-(gi,k)2)4(di,k)2(σi,k)2-hi,k]||≤4(di,k)2(σi,k)2+hi,k、gi,k=Fi,i-Fi,k-Fk,i+Fk,k、 再求解最终的凸优化估计问题得到X中的各个元素的估计值,即x1,x2,…,xN的估计值。
与现有技术相比,本发明的优点在于:本发明方法利用二阶锥松弛技术和半正定松弛技术对通过构造最大似然函数得到的最小化问题进行松弛,得到混合半正定和二阶锥规划的问题的描述,这样可以确保得到全局最优解而不受局部收敛的影响,定位精度高;并且能够有效地抑制测量噪声误差的影响;同时可利用现有技术进行求解出未知目标源的坐标的估计值,降低了锚节点的分布密度,减少了花费成本。
附图说明
图1为本发明方法的总体实现框图;
图2为利用本发明方法和两种现有的基于到达时间的协作定位方法的均方根误差随测量噪声的标准差的变化曲线图;
图3为利用本发明方法和两种现有的基于到达时间的协作定位方法的累计分布函数随估计误差的变化曲线图;
图4为当未知目标源的数目为5个时,利用本发明方法和两种现有的基于到达时间的协作定位方法的均方根误差随锚节点数目变化的曲线图;
图5为当锚节点数目为4个时,利用本发明方法和两种现有的基于到达时间的协作定位方法的均方根误差随未知目标源数目变化的曲线图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明提出的一种应用于无线传感网络的基于到达时间的协作定位方法,其总体实现框图如图1所示,其包括以下步骤:
步骤一:在无线传感网络环境中,建立一个平面直角坐标系或空间坐标系作为参考坐标系,并设定存在M个锚节点和N个未知目标源;然后将M个锚节点在参考坐标系中的坐标对应记为s1,s2,…,sM,将N个未知目标源在参考坐标系中的坐标对应记为x1,x2,…,xN;其中,M≥2,在本实施例中取M=4,N≥2,在本实施例中取N=5,s1为第1个锚节点在参考坐标系中的坐标,s2为第2个锚节点在参考坐标系中的坐标,sM为第M个锚节点在参考坐标系中的坐标,x1为第1个未知目标源在参考坐标系中的坐标,x2为第2个未知目标源在参考坐标系中的坐标,xN为第N个未知目标源在参考坐标系中的坐标。
步骤二:在无线传感网络环境中,由每个未知目标源发射测量信号,每个未知目标源发射的测量信号经过无线传感网络环境传播后由所有锚节点和其它未知目标源接收;先确定每个未知目标源发射的测量信号从发射到各个锚节点接收所经历的时间,并确定每个未知目标源发射的测量信号从发射到其它各个未知目标源接收所经历的时间,将第i个未知目标源发射的测量信号从发射到第j个锚节点接收所经历的时间记为ti,j,将第i个未知目标源发射的测量信号从发射到第k个未知目标源接收所经历的时间记为ti,k;然后计算每个未知目标源发射的测量信号从发射到各个锚节点接收的传输距离测量值,并计算每个未知目标源发射的测量信号从发射到其它各个未知目标源接收的传输距离测量值,将第i个未知目标源发射的测量信号从发射到第j个锚节点接收的传输距离测量值记为di,j,di,j=c×ti,j,将第i个未知目标源发射的测量信号从发射到第k个未知目标源接收的传输距离测量值记为di,k,di,k=c×ti,k;其中,i、j、k均为正整数,1≤i≤N,1≤j≤M,1≤k≤N,k≠i,ti,j和ti,k的单位均为秒,di,j和di,k的单位均为米,c表示光速。
步骤三:将协作定位模型描述为:然后将协作定位模型中的di,j=||xi-sj||+ni,j中的ni,j移到等式左边后对等式两边平方,同样将di,k=||xi-xk||+ni,k中的ni,k移到等式左边后对等式两边平方,得到:由于(ni,j)2远远小于2di,j×ni,j,(ni,k)2远远小于2di,k×ni,k,因此可忽略测量噪声的二次项,即接着忽略中的(ni,j)2和(ni,k)2,得到:再将转换为:其中,符号“|| ||”为求欧氏距离符号,||xi-sj||表示第i个未知目标源与第j个锚节点之间的真实距离,||xi-xk||表示第i个未知目标源与第k个未知目标源之间的真实距离,xi为第i个未知目标源在参考坐标系中的坐标,xk为第k个未知目标源在参考坐标系中的坐标,sj为第j个锚节点在参考坐标系中的坐标,ni,j表示第i个未知目标源与第j个锚节点之间的测量噪声,ni,k表示第i个未知目标源与第k个未知目标源之间的测量噪声。
步骤四:根据构造一个最大似然(ML)函数,得到最小化问题,描述为:很明显最小化问题的描述是一个严重非线性和非凸的函数,无法找到一个闭式解,因此为了将其转化为凸的形式,在此引入一些变量,即然后令 X=[x1,x2,…,xN],其中,min()为取最小值函数,σi,j表示ni,j的标准差,σi,k表示ni,k的标准差,yi,j为引入的变量,为xi的转置,gi,k为引入的变量,为xk的转置,X为引入的变量,符号“[]”为向量或矩阵表示符号,F为引入的变量,I2为二阶的单位矩阵,[X I2]T为[X I2]的转置,XT为X的转置。
步骤五:为了避免最大似然收敛问题,将最小化问题转化成一个优化估计问题,描述为:约束条件为: gi,k=Fi,i-Fi,k-Fk,i+Fk,k、然后利用二阶锥松弛技术对约束条件和进行松弛,将松弛为||[2((di,j)2-(yi,j)2)4(di,j)2×(σi,j)2-ui,j]||≤4(di,j)2×(σi,j)2+ui,j,将松弛为||[2((di,k)2-(gi,k)2)4(di,k)2(σi,k)2-hi,k]||≤4(di,k)2(σi,k)2+hi,k;然而在F上仍然有一个非凸约束,导致该优化估计问题实际仍是非凸的,因此接着利用半正定松弛技术对约束条件进行松弛,将改写为线性矩阵不等式(LMI)形式:其中,ui,j为引入的优化变量,hi,k为引入的优化变量,Fi,i表示F中的第i行第i列元素,为sj的转置,FN+1:N+2,i表示由F中的第N+1行至第N+2行第i列的元素组成的维数为2×1的向量,Fi,k表示F中的第i行第k列元素,Fk,i表示F中的第k行第i列元素,Fk,k表示F中的第k行第k列元素,rank(F)表示求F的秩。
步骤六:确定最终的凸优化估计问题,描述为:约束条件为:||[2((di,j)2-(yi,j)2)4(di,j)2×(σi,j)2-ui,j]||≤4(di,j)2×(σi,j)2+ui,j、||[2((di,k)2-(gi,k)2)4(di,k)2(σi,k)2-hi,k]||≤4(di,k)2(σi,k)2+hi,k、gi,k=Fi,i-Fi,k-Fk,i+Fk,k、 由于最终的凸优化估计问题是混合半正定和二阶锥规划的问题,是一个凸问题,因此再求解最终的凸优化估计问题得到X中的各个元素的估计值,即x1,x2,…,xN的估计值,在MATLAB中可以使用CVX工具箱来求解得到X。
为了验证本发明方法的可行性和有效性,对本发明方法进行仿真试验。
1)测试本发明方法的性能随测量噪声大小的变化情况。设定在无线传感网络环境中存在M=4个锚节点、N=5个未知目标源,测量的方法为:建立一个平面直角坐标系作为参考坐标系,锚节点的坐标随机分布在以(0,0)为圆心,以5米为半径的圆内,未知目标源的位置随机分布在20×20米的正方形内,仿真中假设未知目标源与锚节点之间的测量噪声的功率及未知目标源与未知目标源之间的测量噪声的功率是相同的。
图2给出了利用本发明方法和两种现有的基于到达时间的协作定位方法的均方根误差(RMSE)随测量噪声的标准差的变化曲线图。从图2中可以看出,随着测量噪声的标准差的变大,三种方法以及克拉美-罗界(CRB)的均方根误差均成上升趋势,并且当测量噪声的标准差变大时,三种方法以及克拉美-罗界的均方根误差逐渐增大,本发明方法与克拉美-罗界的均方根误差最接近。
图3给出了利用本发明方法和两种现有的基于到达时间的协作定位方法的累计分布函数(cumulative distribution function,CDF)随估计误差的变化曲线图。从图3中可以看出,在所有估计误差变化范围内,利用本发明方法都具有良好的性能,具体来说,当估计误差为3米时,利用本发明方法的累计分布函数可以达到88.50%,而其它两种方法均达到79.60%;当估计误差为10米时,利用本发明方法的累积分布函数可达98.88%,而现有方法一达到97.24%,现有方法二达到96.44%。因此,足以说明在估计误差同等假设条件下,利用本发明方法比两种现有的基于到达时间的协作定位方法的定位性能更好。
2)分别测试本发明方法的性能随着锚节点或未知目标源个数增加的变化情况。先设定在无线传感网络环境中锚节点的数目从4个增加到10个,再设定未知目标源的数目从4个增加到10个。测量的方法为:建立一个平面直角坐标系作为参考坐标系,锚节点的坐标随机分布在以(0,0)为圆心,以5米为半径的圆内,未知目标源的位置随机分布在20×20米的正方形内。
图4给出了当未知目标源的数目为5个时,利用本发明方法和两种现有的基于到达时间的协作定位方法的均方根误差(RMSE)随锚节点数目变化的曲线图。从图4中可以看出,利用本发明方法在不同锚节点数目下均具有较低的均方根误差,并且随着锚节点数目的增加,利用本发明方法在均方根误差上的降低要优于两种现有的基于到达时间的协作定位方法,并且随着锚节点数目的增加,其均方根误差要小于其它方法更接近克拉美-罗界(CRB)。
图5给出了当锚节点数目为4个时,利用本发明方法和两种现有的基于到达时间的协作定位方法的均方根误差(RMSE)随未知目标源数目变化的曲线图。从图5中可以看出,随着未知目标源数目的增加,利用本发明方法具有更小的均方根误差,特别的,当未知目标源数目大于9时,利用本发明方法对于估计均方根误差这一性能有了更大的提升。
图2至图5中第一种现有的基于到达时间的协作定位方法为Biswas P,Lian T C,Wang T C,et al.Semidefinite programming based algorithms for sensor networklocalization[J].ACM Transactions on Sensor Networks,2006,2(2):188-220.(基于半正定规划的传感器网络定位算法),简称为现有方法一;第二种基于到达时间的协作定位方法为Wang N,Yang L.Further results on cooperative localization viasemidefinite programming[C]//Information Sciences and Systems.IEEE,2011:1-6.(基于半正定规划的协作定位的进一步研究),简称为现有方法二。
由上述仿真结果可以看出,本发明方法具有良好的性能,能够很好地满足定位高精度的需求,并且能够有效地抑制噪声误差的影响,最终求解的是带有较少未知参数的混合半正定和二阶锥规划的问题,因此其计算复杂度低,这充分说明了本发明方法可行且有效。
Claims (1)
1.一种应用于无线传感网络的基于到达时间的协作定位方法,其特征在于包括以下步骤:
步骤一:在无线传感网络环境中,建立一个平面直角坐标系或空间坐标系作为参考坐标系,并设定存在M个锚节点和N个未知目标源;然后将M个锚节点在参考坐标系中的坐标对应记为s1,s2,…,sM,将N个未知目标源在参考坐标系中的坐标对应记为x1,x2,…,xN;其中,M≥2,N≥2,s1为第1个锚节点在参考坐标系中的坐标,s2为第2个锚节点在参考坐标系中的坐标,sM为第M个锚节点在参考坐标系中的坐标,x1为第1个未知目标源在参考坐标系中的坐标,x2为第2个未知目标源在参考坐标系中的坐标,xN为第N个未知目标源在参考坐标系中的坐标;
步骤二:在无线传感网络环境中,由每个未知目标源发射测量信号,每个未知目标源发射的测量信号经过无线传感网络环境传播后由所有锚节点和其它未知目标源接收;先确定每个未知目标源发射的测量信号从发射到各个锚节点接收所经历的时间,并确定每个未知目标源发射的测量信号从发射到其它各个未知目标源接收所经历的时间,将第i个未知目标源发射的测量信号从发射到第j个锚节点接收所经历的时间记为ti,j,将第i个未知目标源发射的测量信号从发射到第k个未知目标源接收所经历的时间记为ti,k;然后计算每个未知目标源发射的测量信号从发射到各个锚节点接收的传输距离测量值,并计算每个未知目标源发射的测量信号从发射到其它各个未知目标源接收的传输距离测量值,将第i个未知目标源发射的测量信号从发射到第j个锚节点接收的传输距离测量值记为di,j,di,j=c×ti,j,将第i个未知目标源发射的测量信号从发射到第k个未知目标源接收的传输距离测量值记为di,k,di,k=c×ti,k;其中,i、j、k均为正整数,1≤i≤N,1≤j≤M,1≤k≤N,k≠i,ti,j和ti,k的单位均为秒,di,j和di,k的单位均为米,c表示光速;
步骤三:将协作定位模型描述为:然后将协作定位模型中的di,j=||xi-sj||+ni,j中的ni,j移到等式左边后对等式两边平方,同样将di,k=||xi-xk||+ni,k中的ni,k移到等式左边后对等式两边平方,得到:接着忽略中的(ni,j)2和(ni,k)2,得到:再将转换为:其中,符号“|| ||”为求欧氏距离符号,||xi-sj||表示第i个未知目标源与第j个锚节点之间的真实距离,||xi-xk||表示第i个未知目标源与第k个未知目标源之间的真实距离,xi为第i个未知目标源在参考坐标系中的坐标,xk为第k个未知目标源在参考坐标系中的坐标,sj为第j个锚节点在参考坐标系中的坐标,ni,j表示第i个未知目标源与第j个锚节点之间的测量噪声,ni,k表示第i个未知目标源与第k个未知目标源之间的测量噪声;
步骤四:根据构造一个最大似然函数,得到最小化问题,描述为:然后令X=[x1,x2,…,xN],其中,min( )为取最小值函数,σi,j表示ni,j的标准差,σi,k表示ni,k的标准差,yi,j为引入的变量,yi,j用于表示第i个未知目标源在参考坐标系中的坐标到第j个锚节点在参考坐标系中的坐标的范数距离,为xi的转置,gi,k为引入的变量,gi,k用于表示第i个未知目标源在参考坐标系中的坐标到第k个未知目标源在参考坐标系中的坐标的范数距离,为xk的转置,X为引入的变量,X用于表示由N个未知目标源在参考坐标系中的坐标构成的矩阵,符号“[]”为向量或矩阵表示符号,F为引入的变量,F用于表示由N个未知目标源在参考坐标系中的坐标构成的矩阵与二阶的单位矩阵构成的特殊矩阵,I2为二阶的单位矩阵,[X I2]T为[X I2]的转置,XT为X的转置;
步骤五:将最小化问题转化成一个优化估计问题,描述为:约束条件为: gi,k=Fi,i-Fi,k-Fk,i+Fk,k、然后利用二阶锥松弛技术对约束条件和进行松弛,将松弛为||[2((di,j)2-(yi,j)2) 4(di,j)2×(σi,j)2-ui,j]||≤4(di,j)2×(σi,j)2+ui,j,将松弛为||[2((di,k)2-(gi,k)2) 4(di,k)2(σi,k)2-hi,k]||≤4(di,k)2(σi,k)2+hi,k;接着利用半正定松弛技术对约束条件进行松弛,将改写为线性矩阵不等式形式:其中,ui,j为引入的优化变量,ui,j用于表示满足第i个未知目标源和第j个锚节点约束条件的待求解的未知整数变量,hi,k为引入的优化变量,hi,k用于表示满足第i个未知目标源和第k个未知目标源约束条件的待求解的未知整数变量,Fi,i表示F中的第i行第i列元素,为sj的转置,FN+1:N+2,i表示由F中的第N+1行至第N+2行第i列的元素组成的维数为2×1的向量,Fi,k表示F中的第i行第k列元素,Fk,i表示F中的第k行第i列元素,Fk,k表示F中的第k行第k列元素,rank(F)表示求F的秩;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810436213.7A CN108668358B (zh) | 2018-05-09 | 2018-05-09 | 一种应用于无线传感网络的基于到达时间的协作定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810436213.7A CN108668358B (zh) | 2018-05-09 | 2018-05-09 | 一种应用于无线传感网络的基于到达时间的协作定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108668358A CN108668358A (zh) | 2018-10-16 |
CN108668358B true CN108668358B (zh) | 2020-08-18 |
Family
ID=63778916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810436213.7A Active CN108668358B (zh) | 2018-05-09 | 2018-05-09 | 一种应用于无线传感网络的基于到达时间的协作定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108668358B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109889971B (zh) * | 2019-01-30 | 2020-09-15 | 浙江大学 | 应用于大型室内环境的基站三维协作定位方法 |
CN110658490B (zh) * | 2019-08-23 | 2021-11-12 | 宁波大学 | 基于rss和aoa的三维无线传感网络非协作定位方法 |
CN110673088B (zh) * | 2019-08-23 | 2021-06-15 | 宁波大学 | 混合视距和非视距环境中基于到达时间的目标定位方法 |
CN110519771A (zh) * | 2019-09-02 | 2019-11-29 | 西北工业大学 | 无线移动传感器网络中基于通信信标的k近邻覆盖方法 |
CN111157943B (zh) * | 2019-12-31 | 2021-10-08 | 宁波大学 | 在异步网络中基于toa的传感器位置误差抑制方法 |
CN111505576B (zh) * | 2020-03-23 | 2022-01-18 | 宁波大学 | 一种针对tdoa定位的传感器选择方法 |
CN111432368B (zh) * | 2020-03-31 | 2022-10-21 | 中国人民解放军国防科技大学 | 一种适用于稀疏锚节点wsn的测距定位方法 |
CN112596026B (zh) * | 2020-12-08 | 2024-01-19 | 江苏智海智能科技有限公司 | 一种未知发射机位置的椭圆目标定位方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105334495A (zh) * | 2015-11-04 | 2016-02-17 | 宁波大学 | 一种无线网络中基于信号到达时间的非视距稳健定位方法 |
CN106019217A (zh) * | 2016-05-12 | 2016-10-12 | 天津大学 | 基于aoa的二维无线传感器网络半定规划定位方法 |
-
2018
- 2018-05-09 CN CN201810436213.7A patent/CN108668358B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105334495A (zh) * | 2015-11-04 | 2016-02-17 | 宁波大学 | 一种无线网络中基于信号到达时间的非视距稳健定位方法 |
CN106019217A (zh) * | 2016-05-12 | 2016-10-12 | 天津大学 | 基于aoa的二维无线传感器网络半定规划定位方法 |
Non-Patent Citations (4)
Title |
---|
"RSS测距模式下无线传感器网络凸松驰定位算法";周雄,陈国彬;《计算机工程》;20150815;全文 * |
Gang Wang;H. Chen."NLOS Error Mitigation for TOA-Based Localization via Convex Relaxation".《IEEE Transactions on Wireless Communications 》.2014, * |
Gang Wang;Kehu Yang."Efficient semidefinite relaxation for energy-based source localization in sensor networks".《2009 IEEE International Conference on Acoustics, Speech and Signal Processing 》.2009, * |
Shengjin Zhang;Shangchao Gao."Robust NLOS Error Mitigation Method for TOA-Based Localization via Second-Order Cone Relaxation".《IEEE Communications Letters》.2015, * |
Also Published As
Publication number | Publication date |
---|---|
CN108668358A (zh) | 2018-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108668358B (zh) | 一种应用于无线传感网络的基于到达时间的协作定位方法 | |
CN110673089B (zh) | 未知视距和非视距分布情况下基于到达时间的定位方法 | |
CN105607039B (zh) | 非视距环境下基于到达时间差的稳健最小二乘定位方法 | |
CN109342993B (zh) | 基于RSS-AoA混合测量的无线传感器网络目标定位方法 | |
CN109917333B (zh) | 融合aoa观测量与tdoa观测量的无源定位方法 | |
CN105334495B (zh) | 一种无线网络中基于信号到达时间的非视距稳健定位方法 | |
CN105425206B (zh) | 一种非同步无线网络中的稳健最小二乘定位方法 | |
CN109597028B (zh) | 一种混合视距非视距条件下基于到达时间的鲁棒定位方法 | |
CN109581281B (zh) | 基于到达时间差和到达频率差的移动目标定位方法 | |
CN110658490B (zh) | 基于rss和aoa的三维无线传感网络非协作定位方法 | |
CN107690184A (zh) | 联合tdoa‑aoa无线传感器网络半定规划定位方法 | |
CN111157943B (zh) | 在异步网络中基于toa的传感器位置误差抑制方法 | |
CN109471061B (zh) | 一种鲁棒处理模型参数误差的接收信号强度差定位方法 | |
CN111929640B (zh) | 一种发送功率未知条件下的传感器网络定位方法 | |
Ampeliotis et al. | Low complexity multiple acoustic source localization in sensor networks based on energy measurements | |
CN110286353B (zh) | 基于非视距环境下RSS-ToA的无线传感器网络目标定位方法 | |
CN110662163A (zh) | 基于rss和aoa的三维无线传感网络协作定位方法 | |
CN112887901A (zh) | 一种基于量化toa量测的凸优化目标定位方法 | |
Bandiera et al. | TDOA localization in asynchronous WSNs | |
CN111505575B (zh) | 基于转化toa模型的针对tdoa定位的传感器选择方法 | |
CN110221245B (zh) | 联合估计目标位置和非视距误差的鲁棒tdoa定位方法 | |
CN103487784B (zh) | 一种基于信号到达时间的定位方法 | |
CN110536410B (zh) | 非视距环境下基于rss和tdoa测量的定位方法 | |
CN110996248A (zh) | 一种基于rss的凸优化无线传感器网络定位方法 | |
CN113923590B (zh) | 一种锚节点位置不确定情况下的toa定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240710 Address after: 2405 Caiwuwei Development Building, 3083 Bao'an South Road, Songyuan Community, Guiyuan Street, Luohu District, Shenzhen City, Guangdong Province 518022 Patentee after: Shenzhen zhenzhiyun Technology Co.,Ltd. Country or region after: China Address before: 315211, Fenghua Road, Jiangbei District, Zhejiang, Ningbo 818 Patentee before: Ningbo University Country or region before: China |
|
TR01 | Transfer of patent right |