[发明专利]一种基于迭代的最小二乘定位方法有效
申请号: | 201210229997.9 | 申请日: | 2012-07-04 |
公开(公告)号: | CN102811419A | 公开(公告)日: | 2012-12-05 |
发明(设计)人: | 江成能;刘珩;卢继华;卜祥元 | 申请(专利权)人: | 北京理工大学 |
主分类号: | H04W4/02 | 分类号: | H04W4/02;H04W64/00 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 100081 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及一种基于迭代的最小二乘定位方法,具体涉及一种在蜂窝通信网络中实现移动站可靠定位的方法,属于无线电定位技术领域。本方法首先对蜂窝基站数据进行分组,再利用球面相交(SSI)技术分别进行单次SSI-LS估计,得到中间估计值,然后计算残差,得到相应的权值,并归一化加权,得到移动站的位置初始估计值,代入残差Taylor法,通过迭代得到最终位置估计值。本发明方法经残差Taylor法的迭代求解,能提升抗NLOS能力;能够实现在蜂窝通信系统中可靠的移动站定位功能以及较强的抗NLOS能力。 | ||
搜索关键词: | 一种 基于 最小 定位 方法 | ||
【主权项】:
1.一种基于迭代的最小二乘定位方法,其特征在于:包括如下步骤:步骤一、对蜂窝基站数据进行分组;参与移动站定位且位置已知的蜂窝基站数量为M个,M≥3,其中服务基站为BS1,其它各蜂窝基站以BSm表示,2≤m≤M;以BS1作为定位参考点,设置BS1位置坐标为二维列向量x1=(x1,y1)T=(0,0)T;第m个基站的位置坐标为已知向量xm=(xm,ym)T;利用服务基站BS1接收移动站信号的到达时延为τ1,第m个基站接收到的来自移动站的信号的到达时延为τm,将所有τm分别与τ1作差,获得其余M-1个蜂窝基站相对于服务基站BS1的M-1个TDOA数据
从M-1个TDOA数据中,任取其中3个TDOA数据作为一个组合,将该组合中TDOA数据的下标首位作为集合Sk,共有N个满足该条件的集合Sk,1≤k≤N,
步骤二、对步骤一得到的分组数据进行中间估计;对N个下标首位集合Sk所对应的TDOA数据利用球面相交技术分别进行单次SSI-LS估计,得到N个中间估计值:
步骤三、利用步骤二得到的中间估计值
计算残差,得到相应的权值;对每个Sk的中间估计值
求全集残差:J ( x ^ k ) = Σ m = 2 M J m ( x ^ k ) ]]>= Σ m = 2 M { [ r ^ m 1 - ( | | x ^ k - x m | | - | | x ^ k - x 1 | | ) 2 ] } ]]> 得到N个全集残差,对其求倒,得到相应的权值wk:w k = 1 J ( x ^ k ) ∀ k ; ]]> 步骤四、对中间估计值
做归一化加权,得到该移动站的3BS2估计值;对步骤二获得的中间估计值
以步骤三得到的相应的权值wk加权,并作归一化处理,得到移动站的位置(x,y)的估计值![]()
x ^ 3 BS 2 = Σ k = 1 N x ^ k w k Σ k = 1 N w k ; ]]> 步骤五、将位置估计值
作为初始值代入残差Taylor法,并通过迭代得到最终位置估计值;所述残差Taylor法为:步骤5.1,将位置估计值
作为初始值,计算定位误差:Δx Δy = ( G t T W - 1 G t ) - 1 G t T W - 1 h t ]]> W=Q·F其中,W为M-1维加权矩阵,符号·为矩阵乘积;Q为M-1维噪声协方差矩阵;F为M-1维残差向量构成的对角方阵:F = diag { J 2 ( x ^ 3 BS 2 ) , J 3 ( x ^ 3 BS 2 ) , . . . , J m ( x ^ 3 BS 2 ) , . . . , J M ( x ^ 3 BS 2 ) } ]]> 数据向量ht和系数矩阵Gt分别为:h t = r ^ 21 - ( r ^ 2 - r ^ 1 ) r ^ 31 - ( r ^ 3 - r ^ 1 ) . . . r ^ M 1 - ( r ^ M - r ^ 1 ) ]]>G t = ( x 1 - x ) / r ^ 1 - ( x 2 - x ) / r ^ 2 ( y 1 - y ) / r ^ 1 - ( y 2 - y ) / r ^ 2 ( x 1 - x ) / r ^ 1 - ( x 3 - x ) / r ^ 3 ( y 1 - y ) / r ^ 1 - ( y 3 - y ) / r ^ 3 . . . . . . ( x 1 - x ) / r ^ 1 - ( x M - x ) / r ^ M ( y 1 - y ) / r ^ 1 - ( y M - y ) / r ^ M . ]]> 距离值
n=1,2,...,M,x=x3BS2,y=y3BS2;步骤5.2,令x=x3BS2+Δx,y=y3BS2+Δy,重复步骤5.1,一直迭代到Δx和Δy满足系统要求的精度;此时所得的x3BS2和y3BS2即为横坐标值和纵坐标值最终的估计值x ^ = ( x 3 BS 2 final , y 3 BS 2 final ) ; ]]> 至此,所述的TR3BS2法处理完毕。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京理工大学,未经北京理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201210229997.9/,转载请声明来源钻瓜专利网。