一种基于曲线拟合的无线传感器网络节点定位方法技术

本发明专利技术公开了一种基于曲线拟合的无线传感器网络节点定位方法，在双移动锚节点的静态路径规划基础上实现。根据待定位节点和移动锚节点运动轨迹距离曲线拟合，再拟合后曲线极点，计算待定位节点位置。在同等路径规划条件下，仿真实验中比较了三边定位法和本发明专利技术所述方法的节点定位精度，结果表明本方法的优越性。

A node location method of wireless sensor network based on curve fitting

【技术实现步骤摘要】
一种基于曲线拟合的无线传感器网络节点定位方法
本专利技术涉及无线传感器应用领域，尤其涉及一种双移动锚节点路径规划下的基于曲线拟合的无线传感器网络节点定位方法。
技术介绍
无线传感器网络是由许多微小的、低成本的、低功耗的、能够互相通信的传感器节点组织而成，具有自组织、自适应和高冗余等特性。它在目标定位、跟踪领域、信息采集、信息管理，环境监测等领域都有广泛的应用前景[1]。而这些应用的前提是传感器网络自身节点的定位准确，因此定位技术对无线传感器网络后续工作的开展起着关键的支撑作用。为节点安装GPS接收器是节点知道自己位置最便捷的方式之一，但是为每个传感器节点配备GPS接收器成本高昂，尤其是大型无线传感器网络。无线传感器网络通常设置一个或多个配备GPS接收器的特殊节点作为信标节点，其他未知节点根据和信标节点的跳数、距离或其他信息推算自身的位置。目前基于移动锚节点节点定位算法可以分为两类，和距离相关以及和距离无关。和距离相关的定位方法比较普遍的是三边定位法。和距离相关的定位算法大多对测距精度要求比较高。本专利技术所述的基于曲线拟合的无线传感器网络节点定位方法也属于和距离相关的定位算法，但是对测距精度要求不高，对测量误差有较好的包容度。仿真结果验证了本算法的有效性。
技术实现思路
有鉴于此，确有必要提供一种双移动锚节点路径规划下的基于曲线拟合的无线传感器网络节点定位方法。本专利技术提出的定位方法是基于2个锚节点按预定路径规划下实施的。两个锚节点按照指定路径遍历无传感器网络，按预定时间间隔发送自身位置和其他信息。待定位节点接收到2段来自移动锚节点的运动轨迹相关信息后，用曲线拟合方法计算自身和锚节点轨迹的最小距离，结合锚节点的位置信息，计算自身的位置坐标。但在某些受限情况，还需要采用三边定位法辅助定位。本专利技术所述的基于曲线拟合的无线传感器网络节点定位方法对测量误差有较好的鲁棒性，定位精度对比三边测量方法提升很多。仿真结果验证了本算法的有效性。为了克服现有技术的缺陷，本专利技术的技术方案如下：一种无线传感器网络下双移动锚节点路径规划和节点定位方法，包括以下步骤：步骤S1：在无线传感器网络中设置两个移动锚节点并预先设定移动锚节点的移动路线，两个移动锚节点按照既定路线移动并按既定时间间隔发送广播信息，所述广播信息至少包括移动锚节点编号和自身位置信息；步骤S2：待定位节点等待接收锚节点数据并记录。如果收到2条或以上来自不同编号的锚节点轨迹，执行步骤3。其中，锚节点轨迹表示按既定时间间隔连续收到锚节点信息所构成锚节点信息集合。如果超过2个既定时间间隔没有再收到锚节点信息，那么认为该锚节点轨迹接收完毕；步骤S3：分别用曲线拟合方法分别求待定位节点和检测到的来自s1和s2两段轨迹的最小距离d1和d2。如果和某一段轨迹的最小距离小于阈值，那么节点采用三边定位法定位，之后加权融合。否则进入下一步；其中，所述步骤S3进一步包括：步骤S31：用曲线拟合方法求待定位节点p和检测到的来自锚节点s1轨迹的最小距离d1。其中，所述步骤S31进一步包括：步骤S311：假设待定位节点p检测到N1个来自锚节点s1的信息，把不同时刻来自移动锚节点s1信息看成不同虚拟锚节点的信息。那么该节点p和由锚节点s1信息构成的虚拟锚节点距离可由RSSI计算得到，但是它们之间的距离是带有噪声的。记录这些距离信息；步骤S312：用多项式最小二乘法对改组距离拟合成2阶多项式曲线。f(x)＝a0+a1x+a2x2(1)其中x表示节点p和虚拟锚节点距离样本，f(x)是对应的2阶多项式函数。步骤S313：定义一个损失函数Q＝∑j(a0+a1x+a2x2-dj)2(2)当损失函数Q取得最小值时，曲线拟合的误差也最小。Q分别对a0，a1，a2的求偏导，函数Q的极值点为偏导为0的点。求解方程，即可得到a0，a1，a2的值。A＝(XTX)-1XTY(3)其中，A＝[a0a1a2]T；X为m×3维的矩阵，N1表示样本的数量；Y是m×1维样本的输出向量。步骤S314：待定位节点p和锚节点s1的距离曲线拟合成抛物线函数f(x)。根据抛物线性质可求出f1(x)的极值点坐标，在xs1处取最小值f1(xs1)。也就是说，待定位节点和锚节点s1轨迹的最短距离d1＝f1(x0)。步骤S32：用同样的方法求待定位节点p和检测到的来自锚节点s2轨迹的最小距离d2＝f2(x0)。步骤S33：判断是否d1＜D，或者d2＜D，这里D表示预设距离阈值。如果满足d1＜D，或者d2＜D，执行步骤4，否则执行步骤5；步骤4：若d1＜D，或者d2＜D执行该步骤。假设待定位节点p检测到M组不共线虚拟锚节点组合，每组包含3个不共线虚拟锚节点。采用三边定位法定位算法，根据每组虚拟锚节点位置和相关信息计算待定位节点坐标，之后将M组的定位结果求平均，得到待定位节点的坐标。其中，所述步骤S4进一步包括：步骤S41：设第j组三边定位信息中虚拟锚节点坐标分别为(xji，yji)，相应的距离待定位节点距离为dji，i＝1，2，3根据Pythagoras定理，可以得到：步骤S42：由式(4)构成的方程组可以推导得到：令A＝2a，根据第j组虚拟锚节点可计算得到该组数据所得的待定位节点位置坐标(xj，yj)为：：步骤S43：将重复步骤S41和S42，直到求出M组虚拟锚节点的所对应的待定位所有估计坐标(xj，yj)，j＝1，2，...，M。再将这M组虚拟锚节点定位结果求平均，最后得到待定位节点p估计坐标如下：步骤5：若d1＜D和d2＜D都不满足，执行步骤5。根据步骤3得到d1和d2可用于推测自身位置。根据公式(9)(10)得到待定位节点p的4个可能定位结果，分别表示为(xp1，yp1a)，(xp1，yp1b)，(xp2，yp2a)，(xp1，yp2b)，具体如下：其中，其中分别表示锚节点s1，s2的坐标，即锚节点si移动过程中，分别对它们的距离曲线拟合后，求出的距离节点p最近的锚节点坐标值。步骤6：计算待定位节点p的横坐标：xp＝(xp1+xp2)/2(11)步骤7：确定待定位节点p的纵坐标。其中，所述步骤S7进一步包括：步骤S71：计算之前求得的最短距离d1，d2的平均误差ε。ε＝(d-d1-d2)/2(12)其中d表示s1和s2运动轨迹之间的距离。步骤S72：对距离d1，d2做误差修正，修正后的最短距离分别表示为dx1，dx2：dx1＝d1-ε(13)dx2＝d2-ε(14)步骤S73：比较yp1a，yp1b和yp2a，yp2b的值，找出数值最相近的2个值。假设yp1b和yp2a取值最相近，那么节点p纵坐标计算如下：步骤8：如果待定位节点又检测到一条来自移动锚节点s1本文档来自技高网...

【技术保护点】
1.一种基于曲线拟合的无线传感器网络节点定位方法，在一定区域中随机布置无线传感器网络，其特征在于，该方法至少包括以下步骤：/n步骤S1：在无线传感器网络中设置两个移动锚节点并预先设定移动锚节点的移动路线，两个移动锚节点按照既定路线移动并按既定时间间隔发送广播信息，所述广播信息至少包括移动锚节点编号和自身位置信息；/n步骤S2：待定位节点等待接收锚节点数据并记录；如果收到2条或以上来自不同编号的锚节点轨迹，执行步骤S3；/n其中，锚节点轨迹表示按既定时间间隔连续收到锚节点信息所构成锚节点信息集合；如果超过2个既定时间间隔没有再收到锚节点信息，那么认为该锚节点轨迹接收完毕；/n步骤S3：分别用曲线拟合方法分别求待定位节点和检测到的来自s

【技术特征摘要】
1.一种基于曲线拟合的无线传感器网络节点定位方法，在一定区域中随机布置无线传感器网络，其特征在于，该方法至少包括以下步骤：
步骤S1：在无线传感器网络中设置两个移动锚节点并预先设定移动锚节点的移动路线，两个移动锚节点按照既定路线移动并按既定时间间隔发送广播信息，所述广播信息至少包括移动锚节点编号和自身位置信息；
步骤S2：待定位节点等待接收锚节点数据并记录；如果收到2条或以上来自不同编号的锚节点轨迹，执行步骤S3；
其中，锚节点轨迹表示按既定时间间隔连续收到锚节点信息所构成锚节点信息集合；如果超过2个既定时间间隔没有再收到锚节点信息，那么认为该锚节点轨迹接收完毕；
步骤S3：分别用曲线拟合方法分别求待定位节点和检测到的来自s1和s2两段轨迹的最小距离d1和d2；如果和某一段轨迹的最小距离小于阈值，那么节点采用三边定位法定位，之后加权融合；否则进入步骤S4；
其中，所述步骤S3进一步包括：
步骤S31：用曲线拟合方法求待定位节点p和检测到的来自锚节点s1轨迹的最小距离d1；
其中，所述步骤S31进一步包括：
步骤S311：假设待定位节点p检测到N1个来自锚节点s1的信息，把不同时刻来自移动锚节点s1信息看成不同虚拟锚节点的信息；
那么该节点p和由锚节点s1信息构成的虚拟锚节点距离可由RSSI计算得到，记录这些距离信息；
步骤S312：用多项式最小二乘法对改组距离拟合成2阶多项式曲线：
f(x)＝a0+a1x+a2x2(1)
其中x表示节点p和虚拟锚节点距离样本，f(x)是对应的2阶多项式函数；
步骤S313：定义一个如下损失函数：
Q＝∑j(a0+a1x+a2x2-dj)2(2)
当损失函数Q取得最小值时，曲线拟合的误差也最小；Q分别对a0，a1，a2的求偏导，函数Q的极值点为偏导为0的点；求解方程，即可得到a0，a1，a2的值；
A＝(XTX)-1XTY(3)
其中，A＝[a0a1a2]T；X为N1×3维的矩阵，N1表示样本的数量；Y是N1×1维样本的输出向量；
步骤S314：待定位节点p和锚节点s1的距离曲线拟合成抛物线函数f1(x)；根据抛物线性质可求出f1(x)的极值点坐标，在xs1处取最小值即，待定位节点和锚节点s1轨迹的最短距离
步骤S32：求待定位节点p和检测到的来自锚节点s2轨迹的最小距离
步骤S33：判断是否d1＜D，或者d2＜D，这里D表示预设距离阈值；如果满足d1＜D，或者d2＜D，执行步骤4，否则执行步骤S5；
步骤S4：若d1＜D，或者d2＜D执行该步骤；假设待定位节点p检测到M组不共线虚拟锚节点组合，每组包含3个不共线虚拟锚节点；采用三边定位法定...

【专利技术属性】
技术研发人员：胡炜薇，陈军相，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33