CN105425206A - 一种非同步无线网络中的稳健最小二乘定位方法 - Google Patents
一种非同步无线网络中的稳健最小二乘定位方法 Download PDFInfo
- Publication number
- CN105425206A CN105425206A CN201510740797.3A CN201510740797A CN105425206A CN 105425206 A CN105425206 A CN 105425206A CN 201510740797 A CN201510740797 A CN 201510740797A CN 105425206 A CN105425206 A CN 105425206A
- Authority
- CN
- China
- Prior art keywords
- sigma
- rho
- eta
- sensor
- delta
- 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.)
- Granted
Links
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种非同步无线网络中的稳健最小二乘定位方法,其先获取未知目标源发射的测量信号经传播到达传感器网络中各个传感器再经中转处理后转发返回到未知目标源时的基于往返到达时间的测量信号等效传输距离测量值;然后根据每个传感器对应的测量信号等效传输距离测量值,获取每个传感器相对应的等效距离测量模型;接着根据重新描述后的距离测量模型建立稳健最小二乘问题;之后通过引入优化变量及利用二阶锥松弛技术,将稳健最小二乘问题松弛为二阶锥规划问题;最后利用内点法技术对二阶锥规划问题求解,得到未知目标源的坐标估计值;优点是能够有效地抑制时钟漂移与中转时间对定位精度的影响,定位精度高,并且有较高的高精度定位稳定性。
Description
技术领域
本发明涉及一种目标定位方法,尤其是涉及一种非同步无线网络中的稳健最小二乘定位方法。
背景技术
目标定位在军事领域诸如精确军事打击中有着不可或缺的作用,现代社会移动互联时代随着基于位置的服务等商业化市场应用的巨大发展,也使高效精确的目标定位研究获得越来越多的关注;同时,目标定位技术在军事侦察、交通监视、家庭自动化、工农业控制、生物医疗、抢险救灾等领域都有广阔的应用前景,因此,研究目标定位方法具有十分重要的意义。而作为GPS等定位***的良好补充,无线网络中的目标定位是一个经典的研究课题。
在目标定位中,基于到达时间的测量值的定位方法占了很大一部分,然而,这种定位方法实现目标定位的前提是假定整个定位网络在时间上完全同步,没有考虑定位网络的非同步性对定位效果的影响,而实际上,由于硬件条件等各种因素,实际网络通常是非同步的,或者说,不可能完全同步,因此这种定位方法难以应用于实际网络。现有的其他常用基于到达时间定位方法也要求精确已知网络的初始传输时间,但这种要求比较难实现或者说代价比较大。
为解决基于到达时间的测量值的定位方法存在的技术问题,非同步无线网络中的定位方法应运而生。在非同步无线网络中的定位方法中,由于传感器时钟存在时钟偏差和时钟漂移,因此在完全未知时钟偏差与时钟漂移的情况下将它们与目标位置联合估计的问题是难以求解的。为了克服这个问题,人们提出了一些方案,比较主流的有:一种是基于到达时间差(TDOA)的定位;另一种是基于往返到达时间(TW-TOA)的定位,图1给出了典型的基于往返到达时间(TW-TOA)的定位环境的示意图。在基于到达时间差的定位方法中,时钟偏差被移除,最终只需对时钟漂移和目标位置进行联合估计,但需要注意的是,这种方法仍然需要定位网络中的传感器在时间上同步;在基于往返到达时间的定位方法中,它不需要定位网络中的任何节点间同步,但需要考虑的是测量中的中转时间对定位结果的影响,如有学者提出的总可行域方法与最小二乘方法,这两种方法都首先对中转时间进行估计,最终通过对时钟漂移和目标源位置的联合估计求得目标位置估计值,但其中对中转时间估计的误差会对定位的性能产生较大的影响。
发明内容
本发明所要解决的技术问题是提供一种非同步无线网络中的稳健最小二乘定位方法,其定位基于往返到达时间,且定位精度高。
本发明解决上述技术问题所采用的技术方案为:一种非同步无线网络中的稳健最小二乘定位方法,其特征在于包括以下步骤:
①在非同步无线网络环境中建立一个二维坐标系或三维坐标系作为参考坐标系,并假设在非同步无线网络环境中存在一个未知目标源和N个位置已知的传感器,且未知目标源在参考坐标系中的坐标为x,N个传感器在参考坐标系中的坐标对应为s1,s2,...,sN,其中,N≥n+1,n表示参考坐标系的维数,s1表示第1个传感器在参考坐标系中的坐标,s2表示第2个传感器在参考坐标系中的坐标,sN表示第N个传感器在参考坐标系中的坐标;
②在非同步无线网络环境中,由未知目标源发射测量信号,测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源,首先确定未知目标源发射的测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源时的时间点与未知目标源发射测量信号时的时间点的时间差,未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源时的时间点与未知目标源发射测量信号时的时间点的时间差为2ti,单位为秒,则其中,1≤i≤N,w表示未知目标源的时钟漂移,si表示第i个传感器在参考坐标系中的坐标,c为光速,Ti表示第i个传感器中转处理未知目标源发射的测量信号所需的中转时间,表示未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源的整条传输路径上的方差为的高斯分布噪声,符号“‖‖”为欧几里德2范数;然后计算未知目标源发射的测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源时的基于往返到达时间的测量信号等效传输距离测量值,未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源时的基于往返到达时间的测量信号等效传输距离测量值为2di,单位为米,则 其中, 表示Ti对di的影响,ni表示di中的噪声,ni服从高斯分布,且ni的方差为
③获取每个传感器相对应的距离测量模型,对于第i个传感器,其相对应的距离测量模型的获取过程为:令w=1+δ,且要求δ满足条件|δ|≤δmax<<1,并确定的取值范围为然后联合和w=1+δ,得到 再联合 和w=1+δ,得到接着根据 和|δ|≤δmax<<1,得到 再假设 则根据 和,得到 之后联合 和 得到 再对 的约等号两边减去的中值得到 最后令 并令 将 简化为并将作为第i个传感器相对应的距离测量模型;其中,δ表示未知目标源相对标准时钟的时钟漂移量,符号“||”为取绝对值符号,符号“<<”为远小于符号,δmax表示未知目标源相对标准时钟的时钟漂移量的最大值,ai和bi对应表示取值的上界和下界,
④对每个传感器相对应的距离测量模型进行重新描述,对于第i个传感器相对应的距离测量模型对其进行重新描述的具体过程为:将 转变为 然后对 的约等号两边进行平方,并假设则省略ni的二次方项得到 再将 转变为: 即 重新描述为
⑤根据重新描述后的距离测量模型,建立一个稳健最小二乘问题,描述为: 然后令 根据 和 将 转变为 再根据 将稳健最小二乘问题描述为: 其中, 表示取使得 的值最小的x, 表示取使得的值最大的{ei},{ei}是指由e1,e2,…,eN组成的集合,f(ei)表示取使得f(ei)的值最大的ei;
⑥确定f(ei)的最大值,如果则f(ei)的最大值为max(f(-ρi),f(ρi));如果则f(ei)的最大值为然后根据和f(ei)的最大值,得到的上镜图形式,描述
为: 其中,符号“||”为取绝对值符号,max()为取最大值函数,其中 表示取使得的值最小的x,{ηi},ηi为 中引入的第i个优化变量,{ηi}为引入的N个优化变量的集合,“s.t.”表示“服从于条件为”;
⑦联合 及 和
⑧在 中引入优化变量y,y=||x||2,然后利用二阶锥松弛技术将y=||x||2松弛为||x||2≤y,得到二阶锥规划问题,描述为: 其中,表示取使得的值最小的x,y,{ηi},符号“[]”为向量表示符号,为si的转置向量;
⑨利用内点法技术对 进行求解,得到x,y,{ηi}对应的估计值,对应记为
所述的步骤③中的ai的值和bi的值的确定过程为:
③-1、假设测试的非同步无线网络中存在位置信息已知的N个传感器;
③-2、由第j'个传感器向第i'个传感器发送测量信号,计算第j'个传感器发射的测量信号经传播到达第i'个传感器再经中转处理后转发返回到第j'个传感器时的时间点与第j'个传感器初始发射测量信号时的时间点的时间差2tj',i'与光速c的乘积,记为 其中,1≤i'≤N,1≤j'≤N,i'≠j';
③-3、根据N个传感器的位置信息,计算第j'个传感器与第i'个传感器之间的真实距离,记为dj',i';然后计算第j'个传感器发射的测量信号经传播到达第i'个传感器后在第i'个传感器的中转时间内的等效传输距离,记为dT,j',i',
③-4、从N个传感器中任意选定一个传感器,假设选定的传感器为第i'个传感器,则令ai=min(dT,j',i'|1≤j'≤N,j'≠i'),并令bi=max(dT,j',i'|1≤j'≤N,j'≠i'),其中,min()为取最小值函数,max()为取最大值函数。
与现有技术相比,本发明的优点在于:相比现有的总可行域方法和最小二乘方法,本发明利用稳健最小二乘方法将时钟漂移与中转时间作为一个无关变量进行稳健处理,只去估计未知目标源在参考坐标系中的坐标值,从而能够有效地抑制时钟漂移与中转时间对定位精度的影响,定位精度高;同时,利用二阶锥松弛技术将稳健最小二乘问题的描述松弛为二阶锥规划问题,这样可以确保得到全局最优解而不受局部收敛的影响,从而有较高的高精度定位稳定性。
附图说明
图1为典型的基于往返到达时间(TW-TOA)的定位环境的示意图;
图2为本发明方法的总体流程框图;
图3为本发明方法及现有的总可行域方法和最小二乘法在测量中均方根误差随噪声大小的变化示意图;
图4为本发明方法及现有的总可行域方法和最小二乘法在测量中均方根误差随传感器(锚节点)数目的变化示意图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明提出的一种非同步无线网络中的稳健最小二乘定位方法,其总体流程框图如图2所示,其包括以下步骤:
①在非同步无线网络环境中建立一个二维坐标系或三维坐标系作为参考坐标系,并假设在非同步无线网络环境中存在一个未知目标源和N个位置已知的传感器,且未知目标源在参考坐标系中的坐标为x,N个传感器在参考坐标系中的坐标对应为s1,s2,...,sN,其中,N≥n+1,n表示参考坐标系的维数,n=2或n=3,即参考坐标系为二维坐标系时n=2,参考坐标系为三维坐标系时n=3,s1表示第1个传感器在参考坐标系中的坐标,s2表示第2个传感器在参考坐标系中的坐标,sN表示第N个传感器在参考坐标系中的坐标。
②在非同步无线网络环境中,如图1所示,由未知目标源发射测量信号,测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源,首先确定未知目标源发射的测量信号经传播到达每个传感器再经中转处理后转发返回至未知目标源时的时间点与未知目标源发射测量信号时的时间点的时间差,未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回至未知目标源时的时间点与未知目标源初始发射测量信号时的时间点的时间差为2ti,单位为秒,则其中,1≤i≤N,w表示未知目标源的时钟漂移,在此w的值未知,si表示第i个传感器在参考坐标系中的坐标,c为光速,Ti表示第i个传感器中转处理未知目标源发射的测量信号所需的中转时间,在此Ti的值未知,表示未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源的整条传输路径上的方差为的高斯噪声,符号“‖‖”为欧几里德2范数;然后计算未知目标源发射的测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源时的基于往返到达时间的测量信号等效传输距离测量值,未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源时的基于往返到达时间的测量信号等效传输距离测量值为2di,单位为米,则 其中, 表示Ti对di的影响,ni表示di中的噪声,ni服从高斯分布,且ni的方差为(即2di中噪声ni的功率),
③获取每个传感器相对应的距离测量模型,对于第i个传感器,其相对应的距离测量模型的获取过程为:令w=1+δ,且要求δ满足条件|δ|≤δmax<<1,并确定的取值范围为然后联合和w=1+δ,得到 再联合 和w=1+δ,得到 接着根据 和|δ|≤δmax<<1,得到 再假设则根据 和得到 之后联合 和 得到 再对 的约等号两边减去的中值得到 最后令 并令 将 简化为并将作为第i个传感器相对应的距离测量模型;其中,δ表示未知目标源相对标准时钟的时钟漂移量,符号“||”为取绝对值符号,符号“”为远小于符号,δmax表示未知目标源相对标准时钟的时钟漂移量的最大值,在此δmax的值已知,a和b对应表示取值的上界和下界,ai的值和bi的值已知,|ei|≤ρi,
在此具体实施例中,步骤③中的ai的值和bi的值的确定过程为:
③-1、假设测试的非同步无线网络存在位置信息已知的N个传感器,传感器的位置信息可由GPS定位得到。
③-2、由第j'个传感器向第i'个传感器发送测量信号,计算第j'个传感器发射的测量信号经传播到达第i'个传感器再经中转处理后转发返回到第j'个传感器时的时间点与第j'个传感器初始发射测量信号时的时间点的时间差2tj',i'与光速c的乘积,记为 其中,1≤i'≤N,1≤j'≤N,i'≠j'。
③-3、根据N个传感器的位置信息,计算第j'个传感器与第i'个传感器之间的真实距离,记为dj',i';然后计算第j'个传感器发射的测量信号经传播到达第i'个传感器后在第i'个传感器的中转时间内的等效传输距离,记为dT,j',i',
③-4、从N个传感器中任意选定一个传感器,假设选定的传感器为第i'个传感器,则令ai=min(dT,j',i'|1≤j'≤N,j'≠i'),并令bi=max(dT,j',i'|1≤j'≤N,j'≠i'),其中,min()为取最小值函数,max()为取最大值函数。
④对每个传感器相对应的距离测量模型进行重新描述,对于第i个传感器相对应的距离测量模型对其进行重新描述的具体过程为:将 转变为 然后对 的约等号两边进行平方,并假设则可以省略ni的二次方项得到 再将 转变为: 即 重新描述为
⑤根据重新描述后的距离测量模型,建立一个稳健最小二乘问题,描述为: 然后令 根据 和 将 转变为 再根据 将稳健最小二乘问题描述为: 其中, 表示取使得 的值最小的x, 表示取使得的值最大的{ei},{ei}是指由e1,e2,…,eN组成的集合,符号“||”为取绝对值符号,f(ei)表示取使得f(ei)的值最大的ei。
⑥确定f(ei)的最大值,如果则f(ei)的最大值为max(f(-ρi),f(ρi));如果则f(ei)的最大值为然后根据和f(ei)的最大值,得到的上镜图形式,描述为: 其中,符号“||”为取绝对值符号,max()为取最大值函数, 表示取使得的值最小的x,{ηi},ηi为 中引入的第i个优化变量,{ηi}为引入的N个优化变量的集合,“s.t.”表示“服从于条件为”。
⑦联合 及 和 得到
⑧在 中引入优化变量y,y=||x||2,然后利用二阶锥松弛技术将y=||x||2松弛为||x||2≤y,得到二阶锥规划问题,描述为: 其中,表示取使得的值最小的x,y,{ηi},符号“[]”为向量表示符号,为si的转置向量。
⑨利用内点法技术对 进行求解,得到x,y,{ηi}对应的估计值,对应记为
为验证本发明方法的可行性和有效性,对本发明方法进行仿真试验。
1)测试本发明方法的性能随测量噪声大小的变化情况。假设使用8个传感器来进行测量,测量的方法为:首先建立一个平面直角坐标系,8个传感器分别在(40,40),(40,-40),(-40,40),(-40,-40),(40,0),(0,40),(40,0),(0,-40)处(单位:m),未知目标源则随机分布在(-40,40)×(-40,40)m2的坐标区域内。在仿真中,未知目标源的时钟漂移w则随机分布在[0.99,1.01]的范围内,也即未知目标源相对标准时钟的时钟漂移量的最大值δmax=0.01,而则假设随机分布的范围为[24,36]m,1≤i≤8,另外,假设所有传感器各自对应的测量信号等效传输距离测量值中的噪声的功率相同,即为其中,σ指图3中的横坐标代表的噪声的标准差差。
图3给出了本发明方法及现有的总可行域方法和线性最小二乘方法在测量中均方根误差随噪声大小的变化示意图。从图3中可以看出在噪声由小至大的变化过程中,本发明方法的定位性能要优于现有的总可行域方法(GTR)和线性最小二乘(LLS)算法。具体来说,在噪声的标准差为1m和9m时,均方根误差能够降低0.5m和2.1m。
2)测试本发明方法的定位精度随着传感器个数增加的变化情况。测量的方法为:在一个平面直角坐标系中假设有8个传感器且8个传感器分别在(40,40),(40,-40),(-40,40),(-40,-40),(40,0),(0,40),(40,0),(0,-40)处(单位:m),未知目标源则随机分布在(-40,40)×(-40,40)m2的坐标区域内,分别取前5至前8个传感器进行定位测试。另外,假设所有传感器各自对应的测量信号等效传输距离测量值中的噪声的标准差相同,即为σ1=σ2=...=σ8=4m。
图4给出了本发明方法及现有的总可行域方法(GTR)和线性最小二乘方法(LLS)在测量中均方根误差随传感器(锚节点)数目的变化示意图。在传感器数目由5个至8个逐渐增加的变化过程中,本发明方法在定位精度方面要优于现有的总可行域方法和线性最小二乘算法。具体来说,在传感器数目少于8个时,总可行域方法和线性最小二乘的性能急剧变差,以致无法完成定位,而本发明方法仍能完成较准确的定位。
由上仿真结果可以看出,本发明方法具有良好的性能。与现有的总可行域方法和线性最小二乘方法相比,本发明方法能够有效的减小均方根误差,提高定位的精度,并且噪声功率的增大并不会显著的减弱定位的性能,体现了定位的稳健性;此外,在网络中传感器比较少的情况下仍能相对准确定位,进一步说明了本发明方法的可行性及有效性。
Claims (2)
1.一种非同步无线网络中的稳健最小二乘定位方法,其特征在于包括以下步骤:
①在非同步无线网络环境中建立一个二维坐标系或三维坐标系作为参考坐标系,并假设在非同步无线网络环境中存在一个未知目标源和N个位置已知的传感器,且未知目标源在参考坐标系中的坐标为x,N个传感器在参考坐标系中的坐标对应为s1,s2,...,sN,其中,N≥n+1,n表示参考坐标系的维数,s1表示第1个传感器在参考坐标系中的坐标,s2表示第2个传感器在参考坐标系中的坐标,sN表示第N个传感器在参考坐标系中的坐标;
②在非同步无线网络环境中,由未知目标源发射测量信号,测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源,首先确定未知目标源发射的测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源时的时间点与未知目标源发射测量信号时的时间点的时间差,未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源时的时间点与未知目标源发射测量信号时的时间点的时间差为2ti,单位为秒,则其中,1≤i≤N,w表示未知目标源的时钟漂移,si表示第i个传感器在参考坐标系中的坐标,c为光速,Ti表示第i个传感器中转处理未知目标源发射的测量信号所需的中转时间,表示未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源的整条传输路径上的方差为的高斯分布噪声,符号“||||”为欧几里德2范数;然后计算未知目标源发射的测量信号经传播到达每个传感器再经中转处理后转发返回到未知目标源时的基于往返到达时间的测量信号等效传输距离测量值,未知目标源发射的测量信号经传播到达第i个传感器再经中转处理后转发返回到未知目标源时的基于往返到达时间的测量信号等效传输距离测量值为2di,单位为米,则 其中, 表示Ti对di的影响,ni表示di中的噪声,ni服从高斯分布,且ni的方差为
③获取每个传感器相对应的距离测量模型,对于第i个传感器,其相对应的距离测量模型的获取过程为:令w=1+δ,且要求δ满足条件|δ|≤δmax<<1,并确定的取值范围为然后联合 和w=1+δ,得到 再联合 和w=1+δ,得到 接着根据 和|δ|≤δmax<<1,得到 再假设则根据 和得到 之后联合 和 得到 再对 的约等号两边减去的中值得到 最后令 并令 将 简化为并将作为第i个传感器相对应的距离测量模型;其中,δ表示未知目标源相对标准时钟的时钟漂移量,符号“||”为取绝对值符号,符号“<<”为远小于符号,δmax表示未知目标源相对标准时钟的时钟漂移量的最大值,ai和bi对应表示取值的上界和下界,
④对每个传感器相对应的距离测量模型进行重新描述,对于第i个传感器相对应的距离测量模型对其进行重新描述的具体过程为:将 转变为 然后对 的约等号两边进行平方,并假设则省略ni的二次方项得到 再将 转变为: 即 重新描述为
⑤根据重新描述后的距离测量模型,建立一个稳健最小二乘问题,描述为: 然后令 根据 和 将 转变为 再根据将稳健最小二乘问题描述为:其中, 表示取使得 的值最小的x, 表示取使得 的值最大的{ei},{ei}是指由e1,e2,…,eN组成的集合,表示取使得f(ei)的值最大的ei;
⑥确定f(ei)的最大值,如果则f(ei)的最大值为max(f(-ρi),f(ρi));如果则f(ei)的最大值为然后根据和f(ei)的最大值,得到的上镜图形式,描述为: 其中,符号“|\”为取绝对值符号,max()为取最大值函数,其中 表示取使得的值最小的x,{ηi},ηi为 中引入的第i个优化变量,{ηi}为引入的N个优化变量的集合,“s.t.”表示“服从于条件为”;
⑦联合 及 和 得到
⑧在 中引入优化变量y,y=||x||2,然后利用二阶锥松弛技术将y=||x||2松弛为||x||2≤y,得到二阶锥规划问题,描述为: 其中,表示取使得的值最小的x,y,{ηi},符号“[]”为向量表示符号,为si的转置向量;
⑨利用内点法技术对 进行求解,得到x,y,{ηi}对应的估计值,对应记为
2.根据权利要求1所述的一种非同步无线网络中的稳健最小二乘定位方法,其特征在于所述的步骤③中的ai的值和bi的值的确定过程为:
③-1、假设测试的非同步无线网络中存在位置信息已知的N个传感器;
③-2、由第j'个传感器向第i'个传感器发送测量信号,计算第j'个传感器发射的测量信号经传播到达第i'个传感器再经中转处理后转发返回到第j'个传感器时的时间点与第j'个传感器初始发射测量信号时的时间点的时间差2tj',i'与光速c的乘积,记为 其中,1≤i'≤N,1≤j'≤N,i'≠j';
③-3、根据N个传感器的位置信息,计算第j'个传感器与第i'个传感器之间的真实距离,记为dj',i';然后计算第j'个传感器发射的测量信号经传播到达第i'个传感器后在第i'个传感器的中转时间内的等效传输距离,记为dT,j',i',
③-4、从N个传感器中任意选定一个传感器,假设选定的传感器为第i'个传感器,则令ai=min(dT,j',i'|1≤j'≤N,j'≠i'),并令bi=max(dT,j',i'|1≤j'≤N,j'≠i'),其中,min()为取最小值函数,max()为取最大值函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510740797.3A CN105425206B (zh) | 2015-11-04 | 2015-11-04 | 一种非同步无线网络中的稳健最小二乘定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510740797.3A CN105425206B (zh) | 2015-11-04 | 2015-11-04 | 一种非同步无线网络中的稳健最小二乘定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105425206A true CN105425206A (zh) | 2016-03-23 |
CN105425206B CN105425206B (zh) | 2017-09-05 |
Family
ID=55503534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510740797.3A Active CN105425206B (zh) | 2015-11-04 | 2015-11-04 | 一种非同步无线网络中的稳健最小二乘定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105425206B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106842121A (zh) * | 2016-11-07 | 2017-06-13 | 宁波大学 | 视距与非视距混合环境中基于到达时间差的稳健定位方法 |
CN106950534A (zh) * | 2017-02-27 | 2017-07-14 | 广东小天才科技有限公司 | 一种空间位置检测方法、***及vr穿戴设备 |
CN107271956A (zh) * | 2017-04-24 | 2017-10-20 | 宁波大学 | 非视距环境中未知起始时间的基于到达时间的定位方法 |
CN107367709A (zh) * | 2017-06-05 | 2017-11-21 | 宁波大学 | 混合环境中基于到达时间鲁棒加权最小二乘定位方法 |
CN107765216A (zh) * | 2017-08-29 | 2018-03-06 | 宁波大学 | 非同步无线网络中的目标位置和时钟参数联合估计方法 |
CN110850366A (zh) * | 2019-10-31 | 2020-02-28 | 宁波大学 | 混合视距和非视距环境下基于接收信号强度的定位方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03122582A (ja) * | 1989-10-05 | 1991-05-24 | Mitsubishi Electric Corp | 追尾装置 |
CN1215474A (zh) * | 1997-01-31 | 1999-04-28 | 摩托罗拉公司 | 在加权最小二乘法位置解决方案中协方差矩阵估计的方法和设备 |
WO2000075681A1 (en) * | 1999-06-09 | 2000-12-14 | Honeywell Inc. | Narrowband based navigation scheme |
CN101680940A (zh) * | 2007-04-13 | 2010-03-24 | 长联通讯技术公司 | 基于加权最小二乘估算来确定地球上的发射器的地理位置解 |
CN103245941A (zh) * | 2013-04-16 | 2013-08-14 | 哈尔滨工程大学 | 一种基于稳健最小二乘的鲁棒波束形成方法 |
CN104820204A (zh) * | 2015-04-13 | 2015-08-05 | 宁波大学 | 一种减小偏差的加权最小二乘定位方法 |
-
2015
- 2015-11-04 CN CN201510740797.3A patent/CN105425206B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03122582A (ja) * | 1989-10-05 | 1991-05-24 | Mitsubishi Electric Corp | 追尾装置 |
CN1215474A (zh) * | 1997-01-31 | 1999-04-28 | 摩托罗拉公司 | 在加权最小二乘法位置解决方案中协方差矩阵估计的方法和设备 |
WO2000075681A1 (en) * | 1999-06-09 | 2000-12-14 | Honeywell Inc. | Narrowband based navigation scheme |
CN101680940A (zh) * | 2007-04-13 | 2010-03-24 | 长联通讯技术公司 | 基于加权最小二乘估算来确定地球上的发射器的地理位置解 |
CN103245941A (zh) * | 2013-04-16 | 2013-08-14 | 哈尔滨工程大学 | 一种基于稳健最小二乘的鲁棒波束形成方法 |
CN104820204A (zh) * | 2015-04-13 | 2015-08-05 | 宁波大学 | 一种减小偏差的加权最小二乘定位方法 |
Non-Patent Citations (3)
Title |
---|
G WANG等: ""Second-order cone relaxation for TOA-Based source localization with unknown start transmission time"", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》 * |
SHENGJIN ZHANG等: ""Robust NLOS error mitigation method for TOA-Based localization via second-order cone relaxation"", 《IEEE COMMUNICATIONS LETTERS》 * |
王刚: ""无线传感器网络中定位与跟踪算法的研究"", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106842121A (zh) * | 2016-11-07 | 2017-06-13 | 宁波大学 | 视距与非视距混合环境中基于到达时间差的稳健定位方法 |
CN106950534A (zh) * | 2017-02-27 | 2017-07-14 | 广东小天才科技有限公司 | 一种空间位置检测方法、***及vr穿戴设备 |
CN106950534B (zh) * | 2017-02-27 | 2020-07-03 | 广东小天才科技有限公司 | 一种空间位置检测方法、***及vr穿戴设备 |
CN107271956A (zh) * | 2017-04-24 | 2017-10-20 | 宁波大学 | 非视距环境中未知起始时间的基于到达时间的定位方法 |
CN107367709A (zh) * | 2017-06-05 | 2017-11-21 | 宁波大学 | 混合环境中基于到达时间鲁棒加权最小二乘定位方法 |
CN107765216A (zh) * | 2017-08-29 | 2018-03-06 | 宁波大学 | 非同步无线网络中的目标位置和时钟参数联合估计方法 |
CN107765216B (zh) * | 2017-08-29 | 2019-10-25 | 宁波大学 | 非同步无线网络中的目标位置和时钟参数联合估计方法 |
CN110850366A (zh) * | 2019-10-31 | 2020-02-28 | 宁波大学 | 混合视距和非视距环境下基于接收信号强度的定位方法 |
CN110850366B (zh) * | 2019-10-31 | 2021-07-16 | 宁波大学 | 混合视距和非视距环境下基于接收信号强度的定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105425206B (zh) | 2017-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105425206A (zh) | 一种非同步无线网络中的稳健最小二乘定位方法 | |
CN108668358B (zh) | 一种应用于无线传感网络的基于到达时间的协作定位方法 | |
CN105334495B (zh) | 一种无线网络中基于信号到达时间的非视距稳健定位方法 | |
CN107690184A (zh) | 联合tdoa‑aoa无线传感器网络半定规划定位方法 | |
EP2817652B1 (en) | Method and system for simultaneous receiver calibration and object localisation for multilateration | |
CN110673089B (zh) | 未知视距和非视距分布情况下基于到达时间的定位方法 | |
CN105607039A (zh) | 非视距环境下基于到达时间差的稳健最小二乘定位方法 | |
CN106446422A (zh) | 一种基于对数似然估计的无源定位跟踪新方法 | |
CN103376447A (zh) | 一种实现非合作双多基地雷达目标三维定位的方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
CN102196559A (zh) | 基于tdoa定位的通道时延误差消除方法 | |
CN107576936A (zh) | 一种去除宽带噪声干扰信号交叉定位虚假点的方法 | |
CN103152826A (zh) | 一种基于nlos状态检测补偿的移动目标追踪方法 | |
CN109581281A (zh) | 基于到达时间差和到达频率差的移动目标定位方法 | |
CN107367709A (zh) | 混合环境中基于到达时间鲁棒加权最小二乘定位方法 | |
CN104820204A (zh) | 一种减小偏差的加权最小二乘定位方法 | |
CN106851821A (zh) | 一种基于无线通信基站的室内三维定位方法 | |
CN103558602A (zh) | 一种用于多基地声纳配置方式的模拟退火定位方法 | |
CN110662163A (zh) | 基于rss和aoa的三维无线传感网络协作定位方法 | |
CN107765216B (zh) | 非同步无线网络中的目标位置和时钟参数联合估计方法 | |
CN111157943A (zh) | 在异步网络中基于toa的传感器位置误差抑制方法 | |
Huang et al. | Analysis of TOA localization with heteroscedastic noises | |
CN104363649A (zh) | 带有约束条件的ukf的wsn节点定位方法 | |
CN114585082B (zh) | 一种电力物联网设备的无线定位方法、装置及存储介质 | |
CN105866811A (zh) | 一种基于地面合作信号的双星座定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |