CN107367709B - 混合环境中基于到达时间鲁棒加权最小二乘定位方法 - Google Patents
混合环境中基于到达时间鲁棒加权最小二乘定位方法 Download PDFInfo
- Publication number
- CN107367709B CN107367709B CN201710411344.5A CN201710411344A CN107367709B CN 107367709 B CN107367709 B CN 107367709B CN 201710411344 A CN201710411344 A CN 201710411344A CN 107367709 B CN107367709 B CN 107367709B
- Authority
- CN
- China
- Prior art keywords
- sight
- line
- target source
- sensor
- coordinate system
- 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 82
- 230000004807 localization Effects 0.000 title abstract description 3
- 230000008054 signal transmission Effects 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims description 27
- 238000006467 substitution reaction Methods 0.000 claims description 21
- 239000000463 material Substances 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 6
- YVPYQUNUQOZFHG-UHFFFAOYSA-N amidotrizoic acid Chemical compound CC(=O)NC1=C(I)C(NC(C)=O)=C(I)C(C(O)=O)=C1I YVPYQUNUQOZFHG-UHFFFAOYSA-N 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 101100391182 Dictyostelium discoideum forI gene Proteins 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 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/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- 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/0257—Hybrid positioning
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种混合环境中基于到达时间鲁棒加权最小二乘定位方法,其在无线传感器网络中建立参考坐标系;然后将目标源与每个传感器之间的信号传输距离以视距模型方式和以非视距模型方式分别进行描述;接着根据视距模型和非视距模型,并采用最坏情况下的鲁棒加权最小二乘方法,得到混合环境中基于到达时间的定位问题,进而得到非凸的定位问题;之后在非凸的定位问题中引入辅助变量,并采用二阶锥松弛方法,得到二阶锥规划问题;最后利用内点法对二阶锥规划问题进行求解,得到目标源在参考坐标系中的坐标位置的最终估计值;优点是其全面的考虑了视距和非视距两种状态,充分利用了路径状态信息并更有针对性的处理了非视距误差,且定位精度高。
Description
技术领域
本发明涉及一种目标定位方法,尤其是涉及一种混合环境中基于到达时间鲁棒加权最小二乘定位方法。
背景技术
随着科学技术的快速发展,目标定位技术在各个领域扮演着越来越重要的角色。例如,诸多的手机应用、目标跟踪与导航、紧急救援、军事作战等。但是,目前的目标定位方法仍然有很多缺陷,如精度不高、易受干扰等。因此,对于精确的目标定位方法的研究仍然十分有价值。
在众多基础的目标定位方法中,使用最为广泛的是基于到达时间(TimeofArrival)的定位方法,其基本网络结构如图1所示。相对于其他方法,基于到达时间的定位方法可以使用较低的成本实现较高的定位精度,同时,该定位方法的计算复杂度相对较低。因此,基于到达时间的定位方法已经得到了普遍的应用。
然而,在现实生活中,对定位精度的影响因素有很多。其中,最主要的因素是由于信号在传输路径上被高大物体阻挡,导致信号的传输发生散射和折射,使信号的传输时间延长,到达时间的测量不准确。本领域内将这种路径叫做非视距路径,导致信号测量产生的距离误差叫做非视距误差,非视距误差的产生原因如图1所示。这种现象在高楼林立的都市,封闭的室内环境中是普遍存在的。因此,要实现非视距环境中精确定位,必须找到合理的高效的定位方法来减少非视距误差带来的不利影响。
目前,基于到达时间的非视距误差的处理方法有很多,主要有以下两种思想:(1)已知路径状态且在视距路径数目大于或等于3条的条件下,将非视距测量全部抛弃,单纯利用视距测量信息进行目标定位,但是此种方法在视距路径数目小于3条的情况下就完全不适用。(2)利用非视距的某些信息,对所有路径进行有关于非视距误差的鲁棒处理,虽然此种方法的适用性与路径状态无关,但是此种方法忽略了部分的路径状态的先验信息,不能有效的提高定位性能。针对上述两种方法的优缺点,应该寻求一种既可以有效的利用路径状态的先验信息,又可以稳健处理非视距误差的方法。然而,目前将视距与非视距联合起来的应用较少,因此,这种研究是十分有必要的。
发明内容
本发明所要解决的技术问题是提供一种混合环境中基于到达时间鲁棒加权最小二乘定位方法,其全面的考虑了视距和非视距两种状态,充分利用了路径状态信息并更有针对性的处理了非视距误差,且定位精度高。
本发明解决上述技术问题所采用的技术方案为:一种混合环境中基于到达时间鲁棒加权最小二乘定位方法,其特征在于包括以下步骤:
①构建一个无线传感器网络,该无线传感器网络中存在一个用于发射测量信号的目标源和N个用于接收测量信号的传感器;然后在该无线传感器网络中,建立一个平面坐标系或空间坐标系作为参考坐标系;接着将N个传感器在参考坐标系中的坐标位置对应记为s1,…,sN,将目标源在参考坐标系中的坐标位置记为x;其中,N≥3,s1表示第1个传感器在参考坐标系中的坐标位置,sN表示第N个传感器在参考坐标系中的坐标位置,s1,…,sN,x分别为2维的列向量;
②计算目标源与每个传感器之间的信号传输距离,将目标源与第i个传感器之间的信号传输距离记为di,di=c×ti,其中,1≤i≤N,c表示光速,ti表示测量信号从目标源发出到第i个传感器接收所经历的时间;
③将目标源与每个传感器之间的信号传输距离以视距模型方式进行描述,将di的视距模型表示为:di=||x-si||+ni,i∈Φa;然后将di=||x-si||+ni,i∈Φa转变为di-ni=||x-si||,i∈Φa;接着对di-ni=||x-si||,i∈Φa的等式两边进行平方并展开,得到再忽略中的得到最后将整理成
将目标源与每个传感器之间的信号传输距离以非视距模型方式进行描述,将di的非视距模型表示为:di=||x-si||+ei+ni,i∈Φb;然后将di=||x-si||+ei+ni,i∈Φb转变为di-ei=||x-si||+ni,i∈Φb;接着对di-ei=||x-si||+ni,i∈Φb的等式两边进行平方并展开,得到再忽略中的得到最后将整理成
基于di=||x-si||+ei+ni,i∈Φb,并且根据已知条件|ni|<<ei,获得一个非视距的约束条件:di≥||x-si||,i∈Φb;然后对di≥||x-si||,i∈Φb的不等式两边进行平方并展开,得到接着将转变为
其中,符号“|| ||”为求欧几里德范数符号,si表示第i个传感器在参考坐标系中的坐标位置,ni表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的测量噪声,ni服从零均值的高斯分布 表示ni的功率,ei表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的非视距误差,符号“| |”为取绝对值符号,ρi表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的非视距误差的上限,Φa表示视距集合,Φa中的元素为测量信号从目标源发出到所有传感器接收所经历的路径中的所有视距路径对应的传感器的序号,Φb表示非视距集合,Φb中的元素为测量信号从目标源发出到所有传感器接收所经历的路径中的所有非视距路径对应的传感器的序号,为si的转置;
④根据和并采用最坏情况下的鲁棒加权最小二乘方法,得到视距与非视距混合环境中基于到达时间的定位问题,描述为:
然后令将转变为
接着在区间ei∈[0,ρi]内确定f(ei)的最大值,当ρi≤di时有当ρi>di时有最后将如果ρi≤di和如果ρi>di代入中,得到非凸的定位问题,描述为:
如果ρi≤di
如果ρi>di
其中,中的是针对测量信号从目标源发出到第i个传感器接收所经历的路径为视距路径时添加的权值,中的是针对测量信号从目标源发出到第i个传感器接收所经历的路径为非视距路径时添加的权值,pi表示针对鲁棒处理非视距误差所增加的权值,min()为取最小值函数,max()为取最大值函数,“s.t.”表示“受约束于……”,f(0)是将ei=0代入得到的,f(ρi)是将ei=ρi代入得到的,f(di)是将ei=di代入得到的,
⑤在非凸的定位问题中引入辅助变量τ1,…,τi,…,τN、η1,…,ηi,…,ηN和辅助变量y,将等价描述为:
然后将代入中,得到
接着采用二阶锥松弛方法将
中的||x||2=y松弛为||x||2≤y,得到二阶锥规划问题,描述为:
⑥利用内点法对二阶锥规划问题进行求解,得到全局最优解,记为x*,x*即为目标源在参考坐标系中的坐标位置的最终估计值。
所述的步骤④中的的获取过程为:
④_1、使其中,max()为取最大值函数,sup()为取上确界函数,符号“| |”为取绝对值符号;
④_2、确定各自的最大值,具体过程为:
④_2a、在ei>>|ni|的前提下,近似得到di≈||x-si||+ei;然后令ri=||x-si||,将di≈||x-si||+ei转变为di≈ri+ei;
④_2b、对于有然后根据di≈ri+ei,近似得到接着根据0≤ei≤ρi,得到最终得到的最大值为
对于有然后根据di≈ri+ei,近似得到接着将求解的最大值分为以下两种情况:当ρi≤di时,始终成立,当0≤ei≤ρi时,得到当ρi>di时,得到而只可能在ei=0或ei=di处取得最大值,进一步有最终得到的最大值为其中,
对于有且ei为正数、ei≤di恒成立,进一步有最终得到的最大值为
④_3、根据ei>>|ni|的情况,从di=ri+ei+ni中得到ri≤di;同时,根据ni服从高斯分布,并根据3σ准则,认为|ni|≤3σi是成立的;然后将ri≤di和|ni|≤3σi代入和中,得到和
④_4、将和代入中,得到
④_5、当ρi≤di时,是始终成立的,这样在ρi≤di情况下,根据将如果ρi≤di简化成如果ρi≤di;
当ρi>di时,得到也是始终成立的;且将转换成等价形式:然后根据和ρi>di的条件,得到接着根据 和将如果ρi>di简化成如果ρi>di;
④_6、综合ρi≤di和ρi>di两种情况,始终有
与现有技术相比,本发明的优点在于:
1)本发明方法在建立目标源与每个传感器之间的信号传输距离的模型的过程中,综合的考虑了视距和非视距两种路径状态,有针对型的处理非视距误差的存在问题,更加具体的运用了先验信息,从而有效的提升了本发明方法的目标定位性能。
2)本发明方法针对视距环境下的距离测量模型处理方式与非视距环境下的距离测量模型处理方式的不平衡性问题,为非视距处理方式选择合适的权值,在保证充分利用视距测量的同时,充分利用非视距信息,从而提高了本发明方法的目标定位性能,这样即使在视距路径非常少的情况下,也具有充足的性能优势,因此本发明方法更加贴近于实际应用。
附图说明
图1为视距与非视距混合环境中基于到达时间(TOA)的目标定位环境的示意图;
图2为本发明方法的总体流程框图;
图3为在视距与非视距混合环境中,已知路径状态的条件下,非视距路径的数目为4条时,利用本发明方法及现有的未区分路径状态的最小二乘方法和现有的区分路径状态的最小二乘方法,获得的目标源的坐标位置估计值与坐标位置真实值的均方根误差随测量噪声增加的变化图;
图4为在视距与非视距混合环境中,已知路径状态的条件下,非视距路径的数目为4条时,利用本发明方法及现有的未区分路径状态的最小二乘方法和现有的区分路径状态的最小二乘方法,获得的目标源的坐标位置估计值与坐标位置真实值的均方根误差随非视距误差的上限增加的变化图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明提出的一种混合环境中基于到达时间鲁棒加权最小二乘定位方法,其总体流程框图如图2所示,其包括以下步骤:
①构建一个无线传感器网络,该无线传感器网络中存在一个用于发射测量信号的目标源和N个用于接收测量信号的传感器;然后在该无线传感器网络中,建立一个平面坐标系或空间坐标系作为参考坐标系;接着将N个传感器在参考坐标系中的坐标位置对应记为s1,…,sN,将目标源在参考坐标系中的坐标位置记为x;其中,N≥3,在本实施例中取N=7,s1表示第1个传感器在参考坐标系中的坐标位置,sN表示第N个传感器在参考坐标系中的坐标位置,s1,…,sN,x分别为2维的列向量。
②计算目标源与每个传感器之间的信号传输距离,将目标源与第i个传感器之间的信号传输距离记为di,di=c×ti,其中,1≤i≤N,c表示光速,ti表示测量信号从目标源发出到第i个传感器接收所经历的时间。
③将目标源与每个传感器之间的信号传输距离以视距模型方式进行描述,将di的视距模型表示为:di=||x-si||+ni,i∈Φa;然后将di=||x-si||+ni,i∈Φa转变为di-ni=||x-si||,i∈Φa;接着对di-ni=||x-si||,i∈Φa的等式两边进行平方并展开,得到由于在现实中一般情况下di>>|ni|是成立的,因此此处再忽略中的得到最后将整理成
将目标源与每个传感器之间的信号传输距离以非视距模型方式进行描述,将di的非视距模型表示为:di=||x-si||+ei+ni,i∈Φb;然后将di=||x-si||+ei+ni,i∈Φb转变为di-ei=||x-si||+ni,i∈Φb;接着对di-ei=||x-si||+ni,i∈Φb的等式两边进行平方并展开,得到由于在现实中一般情况下di>>|ni|是成立的,因此此处再忽略中的得到最后将整理成
基于di=||x-si||+ei+ni,i∈Φb,并且根据已知条件|ni|<<ei,获得一个非视距的约束条件:di≥||x-si||,i∈Φb;然后对di≥||x-si||,i∈Φb的不等式两边进行平方并展开,得到接着将转变为
其中,符号“|| ||”为求欧几里德范数符号,si表示第i个传感器在参考坐标系中的坐标位置,ni表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的测量噪声,ni服从零均值的高斯分布 表示ni的功率,一般情况下0≤|ni|<<ei≤ρi,ei表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的非视距误差,符号“| |”为取绝对值符号,ρi表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的非视距误差的上限,ρi的值在进行目标定位之前已经测量完成,ρi为常量,Φa表示视距集合,Φa中的元素为测量信号从目标源发出到所有传感器接收所经历的路径中的所有视距路径对应的传感器的序号,Φa已知,Φb表示非视距集合,Φb中的元素为测量信号从目标源发出到所有传感器接收所经历的路径中的所有非视距路径对应的传感器的序号,Φb已知,Φa∪Φb=Φ,符号“∪”为集合并运算符号,Φ表示由N个传感器的序号构成的集合,为si的转置。
④根据和并采用最坏情况下的鲁棒加权最小二乘方法,得到视距与非视距混合环境中基于到达时间的定位问题,描述为:
然后令将转变为
接着在区间ei∈[0,ρi]内确定f(ei)的最大值,当ρi≤di时有当ρi>di时有
最后将如果ρi≤di和如果ρi>di代入中,得到非凸的定位问题,描述为:
其中,中的是针对测量信号从目标源发出到第i个传感器接收所经历的路径为视距路径时添加的权值,中的是针对测量信号从目标源发出到第i个传感器接收所经历的路径为非视距路径时添加的权值,pi表示针对鲁棒处理非视距误差所增加的权值,min()为取最小值函数,max()为取最大值函数,“s.t.”表示“受约束于……”,f(0)是将ei=0代入得到的,f(ρi)是将ei=ρi代入得到的,f(di)是将ei=di代入得到的,
在此具体实施例中,步骤④中的的获取过程为:
④_1、为了保证在缓解非视距测量作用的条件下,充分发挥视距测量对目标定位性能的贡献,进行适当的理论推导,使其中,max()为取最大值函数,sup()为取上确界函数,符号“| |”为取绝对值符号。
④_2、确定各自的最大值,具体过程为:
④_2a、在ei>>|ni|的前提下,近似得到di≈||x-si||+ei;然后令ri=||x-si||,将di≈||x-si||+ei转变为di≈ri+ei。
④_2b、对于有然后根据di≈ri+ei,近似得到接着根据0≤ei≤ρi,得到最终得到的最大值为
对于有然后根据di≈ri+ei,近似得到接着将求解的最大值分为以下两种情况:当ρi≤di时,始终成立,因此当0≤ei≤ρi时,得到当ρi>di时,得到而只可能在ei=0或ei=di处取得最大值,进一步有最终得到的最大值为其中,
对于有且ei为正数、ei≤di恒成立,进一步有最终得到的最大值为
④_3、根据ei>>|ni|的情况,从di=ri+ei+ni中得到ri≤di;同时,根据ni服从高斯分布,并根据3σ准则,将有99.73%的概率认为|ni|≤3σi是成立的;然后将ri≤di和|ni|≤3σi代入和中,得到 和
④_4、将和代入中,得到
④_5、接下来通过对分子各项进行比较分析,对不同情况下的pi进行简化,当ρi≤di时,是始终成立的,这样在ρi≤di情况下,根据将如果ρi≤di简化成如果ρi≤di。
当ρi>di时,得到也是始终成立的;且将转换成等价形式:然后根据和ρi>di的条件,得到接着根据 和将如果ρi>di简化成如果ρi>di。
④_6、综合ρi≤di和ρi>di两种情况,可以确定无论在ρi≤di还是ρi>di情况下,始终有
⑤在非凸的定位问题中引入辅助变量τ1,…,τi,…,τN、η1,…,ηi,…,ηN和辅助变量y,将等价描述为:
然后将代入中,得到
接着采用二阶锥松弛方法将
中的||x||2=y松弛为||x||2≤y,得到二阶锥规划问题,描述为:
⑥利用内点法对二阶锥规划问题进行求解,得到全局最优解,记为x*,x*即为目标源在参考坐标系中的坐标位置的最终估计值。
为验证本发明方法的可行性和有效性,对本发明方法进行仿真试验。
假设使用N=7个传感器对本发明方法的性能进行测量,各个传感器均匀分布在以原点(0,0)为中心且半径为5米的圆上,目标源的位置是在30×30m2的方形区域内随机选择。进一步假设目标源到所有传感器之间的7条路径中有4条为非视距路径,其余的路径为视距路径。假设所有传感器的测量噪声的功率(方差)相同,即为非视距路径的非视距误差的上限为ρ1=ρ2=…=ρ7=ρ,即所有非视距路径的非视距误差的上限相同,都为ρ。
图3给出了在视距与非视距混合环境中,已知路径状态的条件下,非视距路径的数目为4条时,利用本发明方法及现有的未区分路径状态的最小二乘方法和现有的区分路径状态的最小二乘方法,获得的目标源的坐标位置估计值与坐标位置真实值的均方根误差(RMSE)随测量噪声增加的变化图;图4给出了在视距与非视距混合环境中,已知路径状态的条件下,非视距路径的数目为4条时,利用本发明方法及现有的未区分路径状态的最小二乘方法和现有的区分路径状态的最小二乘方法,获得的目标源的坐标位置估计值与坐标位置真实值的均方根误差(RMSE)随非视距误差的上限增加的变化图。在仿真图3和图4中,“未区分路径状态的最小二乘方法”表示忽略已知路径状态的前提下,将所有路径视为非视距路径处理的鲁棒最小二乘方法;“区分路径状态的最小二乘方法”表示已知视距与非视距混合环境的前提下,将视距与非视距路径状态分开处理,但未添加合理权值的鲁棒最小二乘算法;本发明方法即表示已知视距与非视距混合环境的前提下,将视距与非视距状态分开处理,同时添加合理权值的鲁棒加权最小二乘算法。从图3中可以看出本发明方法的性能优于其他两种方法;从图4中可以看出,随着非视距误差的上限的增大,本发明方法的优良性能体现的越来越明显,此外本发明方法始终优于未区分路径状态的最小二乘方法和区分路径状态的最小二乘方法。
Claims (2)
1.一种混合环境中基于到达时间鲁棒加权最小二乘定位方法,包括以下步骤①和步骤②:
①构建一个无线传感器网络,该无线传感器网络中存在一个用于发射测量信号的目标源和N个用于接收测量信号的传感器;然后在该无线传感器网络中,建立一个平面坐标系或空间坐标系作为参考坐标系;接着将N个传感器在参考坐标系中的坐标位置对应记为s1,...,sN,将目标源在参考坐标系中的坐标位置记为x;其中,N≥3,s1表示第1个传感器在参考坐标系中的坐标位置,sN表示第N个传感器在参考坐标系中的坐标位置,s1,...,sN,x分别为2维的列向量;
②计算目标源与每个传感器之间的信号传输距离,将目标源与第i个传感器之间的信号传输距离记为di,di=c×ti,其中,1≤i≤N,c表示光速,ti表示测量信号从目标源发出到第i个传感器接收所经历的时间;
其特征在于还包括以下步骤③、步骤④、步骤⑤和步骤⑥:
③将目标源与每个传感器之间的信号传输距离以视距模型方式进行描述,将di的视距模型表示为:di=||x-si||+ni,i∈Φa;然后将di=||x-si||+ni,i∈Φa转变为di-ni=||x-si||,i∈Φa;接着对di-ni=||x-si||,i∈Φa的等式两边进行平方并展开,得到再忽略中的得到最后将整理成
将目标源与每个传感器之间的信号传输距离以非视距模型方式进行描述,将di的非视距模型表示为:di=||x-si||+ei+ni,i∈Φb;然后将di=||x-si||+ei+ni,i∈Φb转变为di-ei=||x-si||+ni,i∈Φb;接着对di-ei=||x-si||+ni,i∈Φb的等式两边进行平方并展开,得到再忽略中的得到最后将整理成
基于di=||x-si||+ei+ni,i∈Φb,并且根据已知条件|ni|<<ei,获得一个非视距的约束条件:di≥||x-si||,i∈Φb;然后对di≥||x-si||,i∈Φb的不等式两边进行平方并展开,得到接着将转变为
其中,符号“|| ||”为求欧几里德范数符号,si表示第i个传感器在参考坐标系中的坐标位置,ni表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的测量噪声,ni服从零均值的高斯分布 表示ni的功率,ei表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的非视距误差,符号“| |”为取绝对值符号,ρi表示测量信号从目标源发出到第i个传感器接收所经历的路径上存在的非视距误差的上限,Φa表示视距集合,Φa中的元素为测量信号从目标源发出到所有传感器接收所经历的路径中的所有视距路径对应的传感器的序号,Φb表示非视距集合,Φb中的元素为测量信号从目标源发出到所有传感器接收所经历的路径中的所有非视距路径对应的传感器的序号,为si的转置;
④根据和并采用最坏情况下的鲁棒加权最小二乘方法,得到视距与非视距混合环境中基于到达时间的定位问题,描述为:
然后令将转变为
接着在区间ei∈[0,ρi]内确定f(ei)的最大值,当ρi≤di时有当ρi>di时有
最后将如果ρi≤di和如果ρi>di代入中,得到非凸的定位问题,描述为:
其中,中的是针对测量信号从目标源发出到第i个传感器接收所经历的路径为视距路径时添加的权值,中的是针对测量信号从目标源发出到第i个传感器接收所经历的路径为非视距路径时添加的权值,pi表示针对鲁棒处理非视距误差所增加的权值,min()为取最小值函数,max()为取最大值函数,“s.t.”表示“受约束于……”,f(0)是将ei=0代入得到的,f(ρi)是将ei=ρi代入得到的,f(di)是将ei=di代入得到的,
⑤在非凸的定位问题中引入辅助变量τ1,...,τi,...,τN、η1,...,ηi,...,ηN和辅助变量y,将等价描述为:
然后将代入中,得到
接着采用二阶锥松弛方法将
中的||x||2=y松弛为||x||2≤y,得到二阶锥规划问题,描述为:
⑥利用内点法对二阶锥规划问题进行求解,得到全局最优解,记为x*,x*即为目标源在参考坐标系中的坐标位置的最终估计值。
2.根据权利要求1所述的混合环境中基于到达时间鲁棒加权最小二乘定位方法,其特征在于所述的步骤④中的的获取过程为:
④_1、使其中,max()为取最大值函数,sup()为取上确界函数,符号“||”为取绝对值符号;
④_2、确定各自的最大值,具体过程为:
④_2a、在ei>>|ni|的前提下,近似得到di≈||x-si||+ei;然后令ri=||x-si||,将di≈||x-si||+ei转变为di≈ri+ei;
④_2b、对于有然后根据di≈ri+ei,近似得到接着根据0≤ei≤ρi,得到最终得到的最大值为
对于有然后根据di≈ri+ei,近似得到接着将求解的最大值分为以下两种情况:当ρi≤di时,始终成立,当0≤ei≤ρi时,得到当ρi>di时,得到而只可能在ei=0或ei=di处取得最大值,进一步有最终得到的最大值为其中,
对于有且ei为正数、ei≤di恒成立,进一步有最终得到的最大值为
④_3、根据ei>>|ni|的情况,从di=ri+ei+ni中得到ri≤di;同时,根据ni服从高斯分布,并根据3σ准则,认为|ni|≤3σi是成立的;然后将ri≤di和|ni|≤3σi代入和中,得到和
④_4、将和代入中,得到
④_5、当ρi≤di时,是始终成立的,这样在ρi≤di情况下,根据将如果ρi≤di简化成如果ρi≤di;
当ρi>di时,得到 也是始终成立的;且将转换成等价形式:然后根据和ρi>di的条件,得到接着根据 将如果ρi>di简化成如果ρi>di;
④_6、综合ρi≤di和ρi>di两种情况,始终有
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710411344.5A CN107367709B (zh) | 2017-06-05 | 2017-06-05 | 混合环境中基于到达时间鲁棒加权最小二乘定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710411344.5A CN107367709B (zh) | 2017-06-05 | 2017-06-05 | 混合环境中基于到达时间鲁棒加权最小二乘定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107367709A CN107367709A (zh) | 2017-11-21 |
CN107367709B true CN107367709B (zh) | 2019-07-26 |
Family
ID=60305297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710411344.5A Active CN107367709B (zh) | 2017-06-05 | 2017-06-05 | 混合环境中基于到达时间鲁棒加权最小二乘定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107367709B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109471061B (zh) * | 2018-09-18 | 2020-07-14 | 宁波大学 | 一种鲁棒处理模型参数误差的接收信号强度差定位方法 |
CN109597028B (zh) * | 2018-11-05 | 2022-11-29 | 宁波大学 | 一种混合视距非视距条件下基于到达时间的鲁棒定位方法 |
CN109975745B (zh) * | 2019-02-28 | 2023-11-14 | 徐州臣赐网络科技有限公司 | 一种基于到达时间差的近远场统一定位方法 |
CN110321902B (zh) * | 2019-05-09 | 2021-07-13 | 哈尔滨工业大学 | 一种基于socp的室内自动视觉指纹采集方法 |
CN110221244B (zh) * | 2019-05-24 | 2022-04-19 | 宁波大学 | 非视距条件下基于到达时间差的鲁棒定位方法 |
CN110221245B (zh) * | 2019-05-28 | 2022-04-19 | 宁波大学 | 联合估计目标位置和非视距误差的鲁棒tdoa定位方法 |
CN110673088B (zh) * | 2019-08-23 | 2021-06-15 | 宁波大学 | 混合视距和非视距环境中基于到达时间的目标定位方法 |
CN110888111B (zh) * | 2019-11-08 | 2021-08-03 | 宁波大学 | 非同步网络中基于toa的鲁棒非视距目标自定位方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101282577B (zh) * | 2008-04-09 | 2011-04-27 | 山东大学 | 基于Ad-hoc网络技术的NLOS定位误差消除方法 |
US9070236B1 (en) * | 2013-01-28 | 2015-06-30 | The United States Of America As Represented By The Secretary Of The Navy | Method and articles of manufacture for time-space-position-information (TSPI) |
CN105425206B (zh) * | 2015-11-04 | 2017-09-05 | 宁波大学 | 一种非同步无线网络中的稳健最小二乘定位方法 |
CN105607039B (zh) * | 2016-01-27 | 2017-11-17 | 宁波大学 | 非视距环境下基于到达时间差的稳健最小二乘定位方法 |
-
2017
- 2017-06-05 CN CN201710411344.5A patent/CN107367709B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107367709A (zh) | 2017-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107367709B (zh) | 混合环境中基于到达时间鲁棒加权最小二乘定位方法 | |
CN105334495B (zh) | 一种无线网络中基于信号到达时间的非视距稳健定位方法 | |
US10145933B2 (en) | Angle determining system and method | |
CN105607039B (zh) | 非视距环境下基于到达时间差的稳健最小二乘定位方法 | |
CN106842121B (zh) | 视距与非视距混合环境中基于到达时间差的稳健定位方法 | |
CN109597028B (zh) | 一种混合视距非视距条件下基于到达时间的鲁棒定位方法 | |
CN108668358B (zh) | 一种应用于无线传感网络的基于到达时间的协作定位方法 | |
CN102209382A (zh) | 一种基于rssi的无线传感器网络节点定位方法 | |
CN107690184A (zh) | 联合tdoa‑aoa无线传感器网络半定规划定位方法 | |
CN110673089B (zh) | 未知视距和非视距分布情况下基于到达时间的定位方法 | |
CN103618997B (zh) | 一种基于信号强度概率的室内定位方法和装置 | |
CN107360552B (zh) | 一种多分类器全局动态融合的室内定位方法 | |
CN103415071A (zh) | 一种水下无线传感器网络定位方法 | |
CN108737952A (zh) | 基于rssi测距改进的多边形加权质心定位方法 | |
CN109471061A (zh) | 一种鲁棒处理模型参数误差的接收信号强度差定位方法 | |
CN111194000B (zh) | 基于蓝牙融合混合滤波与神经网络的测距方法与系统 | |
CN108650629A (zh) | 一种基于无线通信基站的室内三维定位算法 | |
CN106792516B (zh) | 基于无线通信基站的三维定位方法 | |
Bandiera et al. | TDOA localization in asynchronous WSNs | |
CN110673088B (zh) | 混合视距和非视距环境中基于到达时间的目标定位方法 | |
CN110221245B (zh) | 联合估计目标位置和非视距误差的鲁棒tdoa定位方法 | |
CN107994940A (zh) | 一种基于禁忌搜索的可见光定位方法 | |
CN110850366B (zh) | 混合视距和非视距环境下基于接收信号强度的定位方法 | |
He et al. | Hybrid RSS-DOA fingerprint indoor localization algorithm with circular reference points layout | |
CN113573231B (zh) | 一种基于改进三边测量法的室内rssi多点定位方法 |
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 |