一种移动目标RSSI定位的粗大误差识别方法技术

本发明专利技术公开了一种移动目标RSSI定位的粗大误差识别方法，针对接收信号强度识别(RSSI)的最小二乘定位方法对粗大误差敏感性高的问题，提出了一种在只知道锚节点的与未知节点之间的直接测量距离的情况下对该锚节点是否存在粗大误差进行快速识别的方法。本发明专利技术通过计算待识别锚节点与未知节点之间的间接测量距离，进而对直接测量距离分布函数中的参数进行估计，在此情况下给出粗大误差上下限的计算方法，根据粗大误差上下限得到置信区间，然后对粗大误差进行识别。该发明专利技术无需对锚节点的位置进行多次测量，简单方便，特别是在移动目标定位的条件下具有一定的实用价值。

【技术实现步骤摘要】

本专利技术涉及RSSI定位领域，具体涉及一种移动目标RSSI定位的粗大误差识别方法。
技术介绍
当前，基于RSSI最小二乘定位的方法中，很少考虑距离测量粗大误差的识别与去除方法。在常见的误差处理技术中，往往使未知节点在同一位置多次测量与锚节点的位置，利用多次测量值用滤波的方法去除粗大误差。然而，当未知节点是移动的时候，一个时刻可以测量未知节点到各锚节点的位置，但在下一时刻测量时，未知节点已经移动了，无法取得与前一时刻相同的测量结果，也就无法做到重复测量。
技术实现思路
鉴于此，本专利技术的目的是提供一种移动目标RSSI定位的粗大误差识别方法。本专利技术的目的是通过以下技术方案实现的，一种移动目标RSSI定位的粗大误差识别方法，包括以下步骤：S1.选定待评估的锚节点，计算间接测量距离S2.求出间接测量距离的均值和方差；S3.求锚节点与未知节点之间有干扰时的直接测量距离di'的概率分布函数；S4.利用分布函数计算直接测量距离di'的均值和方差；S5.估计未知节点到锚节点之间直接测量距离di'分布函数参数di和标准差σ；S6.根据给定的显著性水平，利用直接测量距离di'的分布函数确定粗大误差上下限，并根据上下限对粗大误差进行识别。进一步，所述步骤S1中，间接测量距离通过以下方法获取：设有四个锚节点i,j,k,m,前三个节点到未知节点的距离已知，最后一个锚节点m(xm,ym)到未知节点的距离未知，利用i,j,k点的信息，确定出未知节点的位置(x,y)，其中未知节点(x,y)到3个锚节点(xi,yi)、(xj,yj)、(xk,yk)的距离di、dj、dk，因此m锚节点到未知节点的距离为：当环境存在干扰的时候，任意三个锚节点i、j、k(i≠j≠k≠m)根据计算出来的距离会各不相同，相当于在真实的距离上增加了一个干扰信号，是一个随机变量。假设有n个锚节点，那么每三个锚节点计算出的间接测量距离，就是该随机变量的一个观察值，该随机变量的分布函数与直接测量距离的分布函数是一致的。进一步，在步骤S2中，间接测量距离的均值和方差分别为：其中，表示从n个锚节点中每次取三个锚节点的组数。进一步，在步骤S3中，直接测量距离di'的概率分布函数为：其中，di表示锚节点和第i个未知节点之间的实际距离，n'为信道衰减指数。进一步，在步骤S4中，直接测量距离di'的均值Edi'和方差D(d′i)分别为：σ表示干扰εi的标准差。进一步，在步骤S5中，未知节点到锚节点之间直接测量距离di'分布函数参数di和标准差σ的估计方法为：因为直接测量距离和间接测量距离的分布函数是一致的，所以其均值和方差也是相同的，所以有：由此得出参数di标准差σ的表达式，进一步，在步骤S6中，粗大误差上下限的确定方法如下：令di'＝y，λ为显著性水平，则粗大误差下限dL用确定，粗大误差的上限dH用确定，由此得到粗大误差上下限分别为，进一步，所述置信区间为：由于采用了上述技术方案，本专利技术具有如下的优点：1、只需对待评估锚节点的位置测量一次就能快速识别粗大误差；2、在基于RSSI的最小二乘定位中运用此方法剔除粗大误差后能显著提高定位精度。附图说明为了使本专利技术的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术作进一步的详细描述，其中：图1为一种移动目标RSSI定位的粗大误差识别方法算法流程图；图2为三边定位示意图；图3为测量的距离无误差时多个锚节点定位示意图；图4为测量的距离有误差时多个锚节点定位示意图；图5为各节点位置分布图；图6为无大干扰锚节点时定位出的未知节点分布图；图7为有一个大干扰锚节点时定位出的未知节点分布图；图8为未知节点到锚节点位置的间接测量方法示意图；图9为粗大误差识别流程图。具体实施方式以下将结合附图，对本专利技术的优选实施例进行详细的描述；应当理解，优选实施例仅为了说明本专利技术，而不是为了限制本专利技术的保护范围。在RSSI定位中，通过测量未知节点(x,y)到3个锚节点(xi,yi)、(xj,yj)、(xk,yk)的距离di、dj、dk，就可以确定出未知节点的坐标。三边定位示意图如图2所示。未知节点坐标根据公式(1)计算。当有多个锚节点时，若测量距离没有误差，以锚节点为圆心、距离为半径的所有圆必定相 交于一点上，任意选择3个锚节点就可以确定出未知节点的位置。如图3所示。当测量的距离有误差的时候，分别以各锚节点为圆心、距离为半径的圆的相交点会有很多个，这时候通常用最小二乘的方法估计未知节点的位置，如图4所示。i个锚节点下的未知节点位置计算方法如下：令则锚节点到未知节点的距离可以根据信号强度测量出来，环境无干扰时，信号强度与实际距离的关系如下：P(di)＝P0-10n'lg(di) (3)当环境有干扰后，测量得到的功率为：P(d′i)＝P0-10n'lg(di)+εi (4)其中εi为干扰，一般假设它是白噪声随机信号，因此测量到的距离也是随机信号，由文献《WSN定位中的RSSI概率质心计算方法》得分布函数如式(5)所示。当其中任意一个锚节点的干扰较大时，采用最小二乘估计未知节点，如果利用这个干扰大的锚节点，反而会降低定位的精度。设有一个未知节点位于(240,270)，各锚节点的位置如图5所示。根据文献《无线网络衰落和损耗的建模与仿真研究》，设锚节点的初始功率为P0＝-44dB，信道衰减系数n'＝2，背景噪声服从的均值为0，标准差为0.8的正态分布。表1列出了9个锚节点在不同的噪声情况下的实际距离和干扰后的距离值。表1根据干扰后的距离利用最小二乘法确定出的未知节点位置为：(248,272)。若其中一个节点的干扰过大，超出了一般的干扰范围，如第I点，干扰标准差为3.5，由此确定出的距离为262，根据这个有粗大误差计算出的位置为：(272,268)。当锚节点I的干扰不大时利用最小二乘估计定位，仿真运行100次，得出的未知节点的位置分布如图6所示。当锚节点I的干扰较大时利用最小二乘估计定位，仿真运行100次，得出的未知节点的位置分布如图7所示。图6和图7中的叉点表示未知节点的真实位置，圆点表示定位出的未知节点的位置。由此可见，将某一个干扰大的锚节点纳入最小二乘计算会对定位精度造成很大影响，应该在计算之前予以识别并剔除。当前，基于RSSI最小二乘定位的方法中，很少考虑距离测量粗大误差的识别与去除方法。在常见的误差处理技术中，往往使未知节点在同一位置多次测量与锚节点的位置，利用多次测量值用滤波的方法去除粗大误差。然而，当未知节点是移动的时候，一个时刻可以测量未知节点到各锚节点的位置，但在下一时刻测量时，未知节点已经移动了，无法取得与前一时刻相同的测量结果，也就无法做到重复测量。为此，本专利技术提出了一种粗大误差识别方法来对移动目标RSSI定位的粗大误差进行识别，该方法包括以下步骤：S1.选定待评估的锚节点，计算间接测量距离设锚节点到未知节点位置的间接测量距离为锚节点到未知节点的实际距离为di，锚节点到未知节点的直接测量距离为d′i。如图8，有四个锚节点i,j,k,m,前三个节点到未知节点的距离已知，最后一个锚节点m到未知节点的距离可以通过前三个节点的信息计算出来。根据公式(1)，利用i,j,k点的信息，可以确定出未知节点的位置(x,y),因此m锚节点到未知节点的距离如式(6本文档来自技高网...

【技术保护点】
一种移动目标RSSI定位的粗大误差识别方法，其特征在于：包括以下步骤：S1.选定待评估的锚节点，计算间接测量距离S2.求出间接测量距离的均值和方差；S3.求锚节点与未知节点之间有干扰时的直接测量距离di'的概率分布函数；S4.利用分布函数计算直接测量距离di'的均值和方差；S5.估计未知节点到锚节点之间直接测量距离di'分布函数参数di和标准差σ；S6.根据给定的显著性水平，利用直接测量距离di'的分布函数确定粗大误差上下限，并根据上下限对粗大误差进行识别。

【技术特征摘要】
1.一种移动目标RSSI定位的粗大误差识别方法，其特征在于：包括以下步骤：S1.选定待评估的锚节点，计算间接测量距离S2.求出间接测量距离的均值和方差；S3.求锚节点与未知节点之间有干扰时的直接测量距离di'的概率分布函数；S4.利用分布函数计算直接测量距离di'的均值和方差；S5.估计未知节点到锚节点之间直接测量距离di'分布函数参数di和标准差σ；S6.根据给定的显著性水平，利用直接测量距离di'的分布函数确定粗大误差上下限，并根据上下限对粗大误差进行识别。2.根据权利要求1所述的移动目标RSSI定位的粗大误差识别方法，其特征在于：所述步骤S1中，间接测量距离通过以下方法获取：设有四个锚节点i,j,k,m,前三个节点到未知节点的距离已知，最后一个锚节点m(xm,ym)到未知节点的距离未知，利用i,j,k点的信息，确定出未知节点的位置(x,y)， x = 1 2 ( d i 2 - d j 2 + x j 2 - x i 2 + y j 2 - y i 2 ) ( y k - y i ) - 1 2 ( d i 2 - d k 2 + x k 2 - x i 2 + y k 2 - y i 2 ) ( y j - y i ) ( x j - x i ) ( y k - y j ) - ( x k - x i ) ( y j - y i ) y = 1 2 ( d i 2 - d j 2 + x j 2 - x i 2 + y j 2 - y i 2 ) ( x k - x i ) - 1 2 ( d i 2 - d k 2 + x k 2 - x i 2 + y k 2 - y i 2 ) ...

【专利技术属性】
技术研发人员：程森林，靳双，孙棣华，王川，韩雨松，徐洋，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆;50