CN103885028B - 一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法 - Google Patents

Abstract

一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法，属于无线传感器定位技术领域。本发明针对现有的分布式传感器定位方法获得位置坐标信息精度不高，计算量大、节点的功率消耗大、计算复杂度高，可定位范围小，锚节点较少的区域中未知节点的定位精度差，甚至有些边缘节点无法定位的问题，提出了一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法：通过无线传感器网络中的已有的锚节点首先完成部分位于锚节点内部的未知节点的定位，其次通过质心定位算法，进行二次定位，并利用误差概率分布校正位置信息，实现了锚节点贫瘠区未知节点的高精度定位。本发明适用于水下传感器定位，森林防火等具体的应用场景。

Description

一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法

技术领域

本发明涉及一种无线传感器网络节点的定位方法，属于无线传感器定位技术领域。

背景技术

无线传感器网络中定位问题是一个非常活跃的研究领域，目前国内外在该方向上已经有很多的综述和论文发表。这些文献的作者讨论了在各种情况下无线传感器网络中的定位方法，分析和评价了这些方法的优劣和适用的条件。

在诸多的应用中，无线传感器网络一个最为重要的用途就是对事件位置信息的获取，譬如水下传感器定位，森林防火等具体的应用场景，因为对于这些事件就是需要知道较为准确的位置信息，否则传感器获得的参数是没有意义的。当利用无线传感器网络对外界进行感知和测量时，传感器网络节点必须先知道自己的位置。传统的传感器节点位置的获取一般包括人工给传感器节点输入位置信息和利用传感器上的GPS设备来获得位置信息，考虑到无线传感器网络的规模和部署的环境，人工注入位置信息的方式并不可取，而为每个传感器配备GPS设备则大大增加了无线传感器网络的成本和功耗。因此，研究传感器网络中仅有部分节点已知自身位置信息的情况下，如何让其余节点准确的获得自身位置信息显得至关重要。

无线传感器网络中的定位方法按照网络定位是否需要获取各节点之间的物理距离，定位方法又可以分为基于测距的定位方法和无需测距的定位方法。

相比于其他形式的网络，无线传感器网络没有基础设施支撑，节点能量供应和成本的限制更为严格。对于存在中心节点的集中式传感器网络，由于其他节点和中心节点之间不断进行控制信息和数据信息的交换，中心节点及其周围节点无疑是整个网络通信和能量的瓶颈，同时中心的损毁必然导致整个网络的不能使用。与之相比，分布式传感器网络在网络健壮性和能量利用方面有着更大的优势。因此，研究在分布式无线传感器网络中的定位问题更有实际意义和应用价值。

很多文献和专利中提出的方法都是考虑如何利用迭代的手段来获取较为精准的位置信息，但是没有考虑利用误差概率分布来校正获得的位置信息，从而得到精度更高的位置信息，本发明提出的基于误差校正的联合质心定位方法则充分考虑这一点，实现了高精度定位。

发明内容

本发明的目的是提出一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法，以解决针对现有的分布式传感器定位方法获得位置坐标信息精度不高，计算量大、节点的功率消耗大、计算复杂度高，可定位范围小，锚节点（位置已知的传感器节点）较少的区域中未知节点的定位精度差，甚至有些边缘节点无法定位的问题。

本发明为解决上述技术问题所采用的技术方案是：

本发明利用先定位出的传感器节点，将第一部分定位出的未知节点视为准锚节点，参与二次迭代定位。由于二次迭代定位等效锚节点增多，因此可以利用质心算法定位锚节点贫瘠的区域，这样可以在较少的分布不均匀的无线传感器网络中扩大定位区域。之后，将利用第一部分的误差概率分布方法校正定位坐标。因此本发明提出的方法对于锚节点较少，且无线传感器节点分布不均匀的无线传感器网络具有十分重要的意义。尤其对于相对偏远的节点，经典的质心定位算法会因为信道条件的恶化，存在很大的迭代累积误差，而基于误差校正的联合质心方法则可以提高定位精度。

本发明所述的一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法，是按照以下步骤实现的：

步骤一、利用误差概率分布实现部分未知节点的高精度定位，具体过程为：

步骤一（一）：通过接收信号强度测算出未知节点与□个锚节点之间的距离,其中n表示在未知节点通信半径内的锚节点数目；

步骤一（二）：解算未知节点的坐标，推导误差函数；

步骤一（三）：通过测量得到误差概率密度函数；

步骤一（四）：利用步骤一（三）所得函数，计算联合概率分布：

每个锚节点到未知节点的距离是独立测量的，联合概率分布函数f为

$f=\underset{i=1}{\overset{2}{\mathrm{\Π}}}f\left({\mathrm{\ρ}}_i\right)---\left(10\right);$

步骤一（五）：计算测量误差的期望值：

$\left\{\begin{array}{c}E(x\_\mathrm{error})=\∫\∫x\_\mathrm{error}\×\mathrm{fd}{\mathrm{\ρ}}_{1}d{\mathrm{\ρ}}_2\\ E(y\_\mathrm{error})=\∫\∫y\_\mathrm{error}\×\mathrm{fd}{\mathrm{\ρ}}_{1}d{\mathrm{\ρ}}_2\end{array}\right.---\left(11\right)$

E(x_error)，E(y_error)分别表示横纵坐标的误差期望值；

步骤一（六）：利用其他组合的锚节点，重复步骤一（四）和步骤一（五），得到M组误差均值，之后求平均值，即

$\left\{\begin{array}{c}\stackrel{\‾}{x\_\mathrm{error}}=\frac{1}{M}\underset{j=1}{\overset{M}{\mathrm{\Σ}}}{E}_j(x\_\mathrm{error})\\ \stackrel{\‾}{y\_\mathrm{error}}=\frac{1}{M}\underset{j=1}{\overset{M}{\mathrm{\Σ}}}{E}_j(y\_\mathrm{error})\end{array}\right.---\left(12\right)$

其中 $M=\left(\begin{array}{c}n\\ 2\end{array}\right)=\frac{n!}{2!\×(n-2)!},$ 分别表示横纵坐标测量误差的平均值；

步骤一（七）：利用未知节点通信半径内三个等效锚节点坐标(x_i,y_i;x_j,y_j;x_k,y_k)估算未知节点的粗略坐标(x_estimate,y_estimate)的均值；

步骤一（八）：求最终坐标(x_end,y_end)：

$\left\{\begin{array}{c}x\_\mathrm{end}=\stackrel{\‾}{x\_\mathrm{estimate}}-\stackrel{\‾}{x\_\mathrm{error}}\\ y\_\mathrm{end}=\stackrel{\‾}{y\_\mathrm{estimate}}-\stackrel{\‾}{y\_\mathrm{error}}\end{array}\right.---\left(15\right);$

步骤一（九）精度评估：利用指标均方根距离来衡量定位精度，其表达式为

rmsd=(||x_end-x_true||²+||y_end-y_true||²)^1/2(16)

其中(x_true,y_true)表示未知节点的真实位置坐标，rmsd为指标均方根距离；

步骤二：利用基于误差校正的联合质心算法完成剩余节点的定位，具体过程为：

步骤二（一）：选择部分已定位出的未知节点，作为准锚节点，和已有的锚节点共同作为等效锚节点，参与其他未知节点的二次定位；

步骤二（二）：通过接收信号强度测算出未知节点与等效锚节点之间的距离；

步骤二（三）：利用未知节点通信半径内的三个等效锚节点坐标，算出未知节点的粗

略坐标： $\left\{\begin{array}{c}x\_\mathrm{estimate}=\frac{a_0{a}_1-a_1{a}_2}{a_3{b}_1-a_2{b}_2}\\ y\_\mathrm{estimate}=\frac{a_0{b}_2-a_1{b}_1}{a_2{b}_2-a_3{b}_1}\end{array}\right.$

$\left(13\right)$

其中

$\left\{\begin{array}{c}{a}_0=0.5[(x_i^2-x_j^2)+(y_i^2-y_j^2)]-0.5(m_i^2-m_j^2)\\ a_1=0.5[(x_j^2-x_k^2)+(y_j^2-y_k^2)]-0.5(m_j^2-m_k^2)\\ a_2=y_i-y_j\\ a_3=y_j-y_k\\ b_1=x_i-x_j\\ b_2=x_j-x_k\end{array}\right.;$

步骤二（四）：利用其它等效锚节点的重新组合，重复步骤二（三）并获得N组坐标粗略估算值，记作(x_eatimate,y_eatimate)_i,i=1,2,...N，其中 $N=\begin{array}{c}\left(\begin{array}{c}q\\ 3\end{array}\right)\end{array}=\frac{q!}{3!\×(q-3)!},$ q表示未知节点通信半径内的锚节点和第一次定位选出的准锚节点的总和；

步骤二（五）：从N组估算出的位置点中选取G个点，其中3≤G≤N，选出的G构成一个点数最多的凸区域；

步骤二（六）：计算凸区域的质心作为未知节点的最终的坐标：

$\left\{\begin{array}{c}x\_\mathrm{end}=\frac{1}{G}\underset{i=1}{\overset{G}{\mathrm{\Σ}}}{(x\_\mathrm{eatimate})}_i\\ y\_\mathrm{end}=\frac{1}{G}\underset{i=1}{\overset{G}{\mathrm{\Σ}}}{(y\_\mathrm{eatimate})}_i\end{array}\right.---\left(17\right);$

步骤二（七）：重复步骤二（二）至步骤二（六），直到完成网络中剩余未知节点的定位。

本发明的有益效果是：

本发明与现有的分布式传感器定位方法相比优点在于：

一、在基于RSSI（接收信号强度）获得的带有误差的位置坐标的基础上，利用误差概率分布特性，校正了获得的坐标信息，使获得的位置坐标信息精度提高了28%～32%。

二、与现有的方法利用迭代的手段相比降低了计算量，节点的功率消耗降低，计算复杂度低。

三、扩大了定位范围，并解决了锚节点（位置已知的传感器节点）较少的区域中未知节点的定位精度差，甚至有些边缘节点无法定位的问题。

附图说明

图1为本发明的差错概率分布定位原理示意图，其中A₁，A₂，A₃，A₄为锚节点，S为未知节点，d₁(1+ρ₁)，d₂(1+ρ₂)，d₃(1+ρ₃)，d₄(1+ρ₄)分别表示未知节点S到锚节点A₁，A₂，A₃，A₄的测量距离；

图2为本发明的基于误差校正的联合质心方法定位原理图；

图3为本发明的定位节点分布图；

图4为本发明的定位误差分布图；

图5(a)为服从均匀分布时，现有的最小二乘估计算法估计误差分布直方图；

图5(b)为服从均匀分布时，本发明的基于误差校正的联合质心定位方法的估计误差分布直方图；

图6为服从均匀分布时，现有的最小二乘估计算法和本发明的基于误差校正的联合质心定位方法的累积分布函数，CDF表示累积分布函数；

图7(a)为服从正态分布时，现有的最小二乘估计算法估计误差分布直方图；

图7(b)为服从正态分布时，本发明的基于误差校正的联合质心定位方法的估计误差分布直方图；

图8为服从正态分布时，现有的最小二乘估计算法和本发明的基于误差校正的联合质心定位方法的累积分布函数；

图9(a)为现有的质心定位算法估计误差分布直方图；

图9(b)为本发明的基于误差校正的联合质心定位方法估计误差分布直方图；

图10为现有的质心算法与本发明的基于误差校正的联合质心定位方法的累积概率密度函数。

具体实施方式

具体实施方式一：本实施方式所述的一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法，是按照以下步骤实现的：

步骤一、利用误差概率分布实现部分未知节点的高精度定位，具体过程为：

步骤一（一）：通过接收信号强度测算出未知节点与□个锚节点之间的距离,其中n表示在未知节点通信半径内的锚节点数目；

步骤一（二）：解算未知节点的坐标，推导误差函数；

步骤一（三）：通过测量得到误差概率密度函数；

步骤一（四）：利用步骤一（三）所得函数，计算联合概率分布：

每个锚节点到未知节点的距离是独立测量的，联合概率分布函数f为

$f=\underset{i=1}{\overset{2}{\mathrm{\Π}}}f\left({\mathrm{\ρ}}_i\right)---\left(10\right);$

步骤一（五）：计算测量误差的期望值：

$\left\{\begin{array}{c}E(x\_\mathrm{error})=\∫\∫x\_\mathrm{error}\×\mathrm{fd}{\mathrm{\ρ}}_{1}d{\mathrm{\ρ}}_2\\ E(y\_\mathrm{error})=\∫\∫y\_\mathrm{error}\×\mathrm{fd}{\mathrm{\ρ}}_{1}d{\mathrm{\ρ}}_2\end{array}\right.---\left(11\right)$

E(x_error)，E(y_error)分别表示横纵坐标的误差期望值；

步骤一（六）：利用其他组合的锚节点，重复步骤一（四）和步骤一（五），得到M组误差均值，之后求平均值，即

$\left\{\begin{array}{c}\stackrel{\‾}{x\_\mathrm{error}}=\frac{1}{M}\underset{j=1}{\overset{M}{\mathrm{\Σ}}}{E}_j(x\_\mathrm{error})\\ \stackrel{\‾}{y\_\mathrm{error}}=\frac{1}{M}\underset{j=1}{\overset{M}{\mathrm{\Σ}}}{E}_j(y\_\mathrm{error})\end{array}\right.---\left(12\right)$

其中 $M=\left(\begin{array}{c}n\\ 2\end{array}\right)=\frac{n!}{2!\×(n-2)!},$ 分别表示横纵坐标测量误差的平均值；

步骤一（七）：利用未知节点通信半径内三个等效锚节点坐标(x_i,y_i;x_j,y_j;x_k,y_k)估算未知节点的粗略坐标(x_estimate,y_estimate)的均值；

步骤一（八）：求最终坐标(x_end,y_end)：

$\left\{\begin{array}{c}x\_\mathrm{end}=\stackrel{\‾}{x\_\mathrm{estimate}}-\stackrel{\‾}{x\_\mathrm{error}}\\ y\_\mathrm{end}=\stackrel{\‾}{y\_\mathrm{estimate}}-\stackrel{\‾}{y\_\mathrm{error}}\end{array}\right.---\left(15\right);$

步骤一（九）精度评估：利用指标均方根距离来衡量定位精度，其表达式为

rmsd=(||x_end-x_true||²+||y_end-y_true||²)^1/2(16)

其中(x_true,y_true)表示未知节点的真实位置坐标，rmsd为指标均方根距离；

步骤二：利用基于误差校正的联合质心算法完成剩余节点的定位，具体过程为：

步骤二（一）：选择部分已定位出的未知节点，作为准锚节点，和已有的锚节点共同作为等效锚节点，参与其他未知节点的二次定位；

步骤二（二）：通过接收信号强度测算出未知节点与等效锚节点之间的距离；

步骤二（三）：利用未知节点通信半径内的三个等效锚节点坐标，算出未知节点的粗

略坐标： $\left\{\begin{array}{c}x\_\mathrm{estimate}=\frac{a_0{a}_1-a_1{a}_2}{a_3{b}_1-a_2{b}_2}\\ y\_\mathrm{estimate}=\frac{a_0{b}_2-a_1{b}_1}{a_2{b}_2-a_3{b}_1}\end{array}\right.$

$\left(13\right)$

其中

$\left\{\begin{array}{c}{a}_0=0.5[(x_i^2-x_j^2)+(y_i^2-y_j^2)]-0.5(m_i^2-m_j^2)\\ a_1=0.5[(x_j^2-x_k^2)+(y_j^2-y_k^2)]-0.5(m_j^2-m_k^2)\\ a_2=y_i-y_j\\ a_3=y_j-y_k\\ b_1=x_i-x_j\\ b_2=x_j-x_k\end{array}\right.;$

步骤二（四）：利用其它等效锚节点的重新组合，重复步骤二（三）并获得N组坐标粗略估算值，记作(x_eatimate,y_eatimate)_i,i=1,2,...N，其中 $N=\begin{array}{c}\left(\begin{array}{c}q\\ 3\end{array}\right)\end{array}=\frac{q!}{3!\×(q-3)!},$ q表示未知节点通信半径内的锚节点和第一次定位选出的准锚节点的总和；

步骤二（五）：从N组估算出的位置点中选取G个点，其中3≤G≤N，选出的G构成一个点数最多的凸区域；

步骤二（六）：计算凸区域的质心作为未知节点的最终的坐标：

$\left\{\begin{array}{c}x\_\mathrm{end}=\frac{1}{G}\underset{i=1}{\overset{G}{\mathrm{\Σ}}}{(x\_\mathrm{eatimate})}_i\\ y\_\mathrm{end}=\frac{1}{G}\underset{i=1}{\overset{G}{\mathrm{\Σ}}}{(y\_\mathrm{eatimate})}_i\end{array}\right.---\left(17\right);$

步骤二（七）：重复步骤二（二）至步骤二（六），直到完成网络中剩余未知节点的定位。

具体实施方式二：本实施方式与具体实施方式一不同的是：步骤一（二）所述的解算未知节点的坐标过程为：

利用未知节点通信半径内的任意两个锚节点的距离解算如下：

$\left\{\begin{array}{c}{m}_i=\sqrt{{(x_i-x)}^2+{(y_i-y)}^2}\\ m_j=\sqrt{{(x_j-x)}^2+{(y_j-y)}^2}\end{array}\right.---\left(1\right)$

(x_i,y_i;x_j,y_j)表示锚节点的坐标，m_i、m_j表示未知节点到锚节点i,j的测量距离,x和y表示未知节点的坐标，其中i,j=1,2,...,n，且

m_k=d_k+λ_k=d_k(1+ρ_k)(2)

其中，k=1,2,...,n，m_k表示未知节点到锚节点k的测量距离，d_k表示未知节点到锚节点k的真实距离，λ_k表示未知节点到锚节点k测量误差，ρ_k表示未知节点到锚节点k测量的不确定度；

结合(1)式和(2)式，得出未知节点的坐标为：

$x=\frac{1}{4}[(d_i^2+2d_i{\mathrm{\λ}}_i+{\mathrm{\λ}}_i^2)-(d_j^2+2d_j{\mathrm{\λ}}_j+{\mathrm{\λ}}_j^2)]---\left(18\right)$

$\begin{array}{c}y=1\±\frac{1}{4}(-16+8d_i^2+8d_j^2-d_i^4-d_j^4+2d_j^2{d}_i^2\\ +16d_i{\mathrm{\λ}}_i+8{\mathrm{\λ}}_i^2+8d_i{\mathrm{\λ}}_i{d}_j{\mathrm{\λ}}_j-4d_j^3{\mathrm{\λ}}_j\\ -16d_j^2{\mathrm{\λ}}_j^2+2d_j^2{\mathrm{\λ}}_i^2-4d_j{\mathrm{\λ}}_j^3+2d_i^2{\mathrm{\λ}}_j^2\\ +2{\mathrm{\λ}}_i^2{\mathrm{\λ}}_j^2-4d_i^3{\mathrm{\λ}}_i-6d_i^2{\mathrm{\λ}}_i^2-4d_i{\mathrm{\λ}}_i^3\\ -{\mathrm{\λ}}_i^4-{\mathrm{\λ}}_j^4+16d_j{\mathrm{\λ}}_j+8{\mathrm{\λ}}_j^2+4d_j^2{d}_i{\mathrm{\λ}}_i\\ +4d_i^2{d}_j{\mathrm{\λ}}_j+4{\mathrm{\λ}}_i^2{d}_j{\mathrm{\λ}}_j+4{\mathrm{\λ}}_j^2{d}_i{\mathrm{\λ}}_i{)}^{1/2}\end{array}---\left(19\right).$

其它步骤及参数与具体实施方式一相同。

具体实施方式三：本实施方式与具体实施方式一或二不同的是：步骤一（二）所述的推导误差函数过程为：

由于误差相对于距离而言是很小的，因此误差的平方项作差更小，可以忽略，由(18)和（19）式可得

$\begin{array}{c}x=\frac{1}{4}[(d_i^2+2d_i{\mathrm{\λ}}_i+{\mathrm{\λ}}_i^2)-(d_j^2+2d_j{\mathrm{\λ}}_j+{\mathrm{\λ}}_j^2)]\\ =\frac{1}{4}(d_i^2-d_j^2)+\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)+\frac{1}{4}({\mathrm{\λ}}_i^2-{\mathrm{\λ}}_j^2)\\ \≅\frac{1}{4}(d_i^2-d_j^2)+\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)\\ =\mathrm{const}+\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)\end{array}---\left(3\right)$

其中 $\mathrm{const}=\frac{1}{4}(d_i^2-d_j^2)$ 表示常数，

则坐标的误差为：

$\begin{array}{c}x\_\mathrm{error}=\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)\\ =\frac{1}{2}(d_i^2{\mathrm{\ρ}}_i-d_j^2{\mathrm{\ρ}}_j)\end{array}---\left(4\right)$

同理可得

$\begin{array}{c}y\_\mathrm{error}=\frac{1}{4}(\frac{{\mathrm{\β}}_1}{2\sqrt{\mathrm{\α}}}{\mathrm{\λ}}_i+\frac{{\mathrm{\β}}_2}{2\sqrt{\mathrm{\α}}}{\mathrm{\λ}}_j)\\ =\frac{1}{4}(\frac{d_1{\mathrm{\β}}_1}{2\sqrt{\mathrm{\α}}}{\mathrm{\ρ}}_i+\frac{d_2{\mathrm{\β}}_2}{2\sqrt{\mathrm{\α}}}{\mathrm{\ρ}}_j)\end{array}---\left(5\right)$

其中x_error,y_error分别表示坐标x,y的误差，且令

$\left\{\begin{array}{c}\mathrm{\α}=-16+8d_j^2+8d_i^2-d_j^4-d_i^4+2d_i^2{d}_j^2\\ {\mathrm{\β}}_1=16d_i-4d_i^3+4d_j^3{d}_i\\ {\mathrm{\β}}_2=16d_j-4d_j^3+4d_i^3{d}_j\end{array}\right.,$

将d_i=m_i/(1+ρ_i)和带入(4)和(5)式得

$\left\{\begin{array}{c}x\_\mathrm{error}=\frac{1}{2}(\frac{m_i^2{\mathrm{\ρ}}_i}{(1+{\mathrm{\ρ}}_i)}-\frac{m_j^2{\mathrm{\ρ}}_j}{(1+{\mathrm{\ρ}}_j)})\\ y\_\mathrm{error}=\frac{1}{4}(g\left({\mathrm{\ρ}}_i\right)+g\left({\mathrm{\ρ}}_j\right))\end{array}\right.---\left(6\right).$

其它步骤及参数与具体实施方式一或二相同。

具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：步骤一（三）所述的误差概率密度函数如下：

(1)无遮挡的环境下服从正态分布： $f\left({\mathrm{\ρ}}_i\right)=\frac{1}{\mathrm{\σ}\sqrt{2\mathrm{\π}}}{e}^{-\frac{{\mathrm{\ρ}}_i^2}{2{\mathrm{\σ}}^2}},{\mathrm{\ρ}}_i\∈(-\mathrm{\δ},\mathrm{\δ})---\left(7\right)$

(2)有遮挡的环境下服从瑞利分布： $f\left({\mathrm{\ρ}}_i\right)=\frac{{\mathrm{\ρ}}_i}{{\mathrm{\σ}}^2}\exp (-\frac{{\mathrm{\ρ}}_i^2}{2{\mathrm{\σ}}^2}),{\mathrm{\ρ}}_i\≥0---\left(8\right)$

(3)空旷环境下服从均匀分布： $f\left({\mathrm{\ρ}}_i\right)=\frac{1}{2\mathrm{\δ}},{\mathrm{\ρ}}_i\∈(-\mathrm{\δ},\mathrm{\δ})---\left(9\right)$

其中f(ρ_i)表示第i个测量不确定度ρ_i的概率密度函数,i=1,2，σ表示标准差，δ表示误差范围。

其它步骤及参数与具体实施方式一至三之一相同。

具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：步骤一（七）所述的估算未知节点的粗略坐标(x_estimate,y_estimate)的均值过程为：根据式（13）进行 $M^{\′}=\left(\begin{array}{c}n\\ 3\end{array}\right)$ 次估算，可得估算均值为

$\left\{\begin{array}{c}\stackrel{\‾}{x\_\mathrm{estimate}}=\frac{1}{M^{\′}}\underset{j=1}{\overset{M^{\′}}{\mathrm{\Σ}}}{(x\_\mathrm{estimate})}_j\\ \stackrel{\‾}{y\_\mathrm{estimate}}=\frac{1}{M^{\′}}\underset{j=1}{\overset{M^{\′}}{\mathrm{\Σ}}}{(y\_\mathrm{estimate})}_j\end{array}\right.---\left(14\right).$

其它步骤及参数与具体实施方式一至四之一相同。

本发明提出的一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法，在每一个未知节点定位时必须在其通信半径范围内有至少两个位置信息已知的传感器节点，即锚节点；且通过测量获得误差概率分布情况的前提下适用于节点分布不均匀的无线传感器网路。基于误差校正的联合质心定位方法主要是利用误差概率分布来校正通过RSSI（接收信号强度）估算出的位置信息，提高了定位精度。通过无线传感器网络中的已有的锚节点首先完成部分位于锚节点内部的未知节点的定位。其次通过质心定位算法，进行二次定位，并利用误差概率分布校正位置信息，实现了锚节点贫瘠区未知节点的高精度定位。

实施例一、定位效果：图3中“□”表示锚节点，“☆”表示未知节点，共30个节点，其中24个未知节点，分布在100×100m的区域内。图4表示定位结束后误差分布图，当“圆圈”和“五角星”重合时表示定位成功，没有圆圈表示未定位成功。很明显20个未知节点定位成功。

实施例二、当误差服从均匀分布时，基于100×100m区域内随机部署30个点，仍然是24个未知节点，采用现有的最小二乘估计算法和本发明方法分别进行定位，重复此过程10000次，统计得到了图5与图6的仿真结果。

图5给出了现有的经典的最小二乘估计算法的误差分布和本发明提出的基于误差校正的联合质心定位方法误差分布，很明显本发明方法的定位精度要好于经典的方法。

图6中的两种算法的累积分布函数可以更加清晰地看出本发明的优势。

比较衡量指标RMSD（均方根距离），最小二乘法与基于误差校正的联合质心定位方法的比值为：0.017：0.0078。

实施例三、当误差服从正态分布时，基于100×100m区域内随机部署30个点，其中24个未知节点，采用现有的最小二乘估计算法和本发明方法分别进行定位，重复此过程10000次，统计得到了图7与图8的仿真结果。

从图7和图8可以看出，当传感器网络处于的环境服从正态分布时，本发明提出的定位方法优势更加明显。

比较衡量指标RMSD（均方根距离），现有的最小二乘法与本发明提出的基于误差校正的联合质心定位方法比值为：0.00133：0.00097。

实施例四、基于100×100m区域内随机部署30个点，其中24个未知节点，采用本发明提出的基于误差校正的联合质心定位方法与现有的质心定位方法分别进行定位，重复此过程10000次，统计得到了图9与图10的仿真结果。

本发明提出的基于误差校正的联合质心定位方法与现有的质心定位方法的仿真结果比较分析：由图9和图10可以看出本发明提出的基于误差校正的联合质心定位方法的定位精度更高，并且扩大了定位的范围。

Claims (2)

1.一种适用于无线传感器网络节点定位的基于误差校正的联合质心定位方法，其特征在于所述方法是按照以下步骤实现的：

步骤一、利用误差概率分布实现部分未知节点的高精度定位，具体过程为：

步骤一(一)：通过接收信号强度测算出未知节点与n个锚节点之间的距离,其中n表示在未知节点通信半径内的锚节点数目；

步骤一(二)：解算未知节点的坐标，推导误差函数；

利用未知节点通信半径内的任意两个锚节点的距离解算如下：

$\left\{\begin{array}{c}{m}_i=\sqrt{{(x_i-x)}^2+{(y_i-y)}^2}\\ m_j=\sqrt{{(x_j-x)}^2+{(y_j-y)}^2}\end{array}\right.---\left(1\right)$

(x_i,y_i；x_j,y_j)表示锚节点的坐标，m_i、m_j表示未知节点到锚节点i,j的测量距离,其中i,j＝1,2,…,n；

m_k＝d_k+λ_k＝d_k(1+ρ_k)(2)

其中，k＝1,2,…,n，m_k表示未知节点到锚节点k的测量距离，d_k表示未知节点到锚节点k的真实距离，λ_k表示未知节点到锚节点k测量误差，ρ_k表示未知节点到锚节点k测量的不确定度；

x和y表示未知节点的坐标，结合(1)式和(2)式，得出未知节点的坐标为：

$x=\frac{1}{4}\[(d_i^2+2d_i{\mathrm{\λ}}_i+{\mathrm{\λ}}_i^2)-(d_j^2+2d_j{\mathrm{\λ}}_j+{\mathrm{\λ}}_j^2)\]---\left(18\right)$

$\begin{array}{c}y=1\±\frac{1}{4}(-16+8d_i^2+8d_j^2-d_i^4-d_j^4+2d_j^2{d}_i^2\\ +16d_i{\mathrm{\λ}}_i+8{\mathrm{\λ}}_i^2+8d_i{\mathrm{\λ}}_i{d}_j{\mathrm{\λ}}_j-4d_j^3{\mathrm{\λ}}_j\\ -16d_j^2{\mathrm{\λ}}_j^2+2d_j^2{\mathrm{\λ}}_i^2-4d_j{\mathrm{\λ}}_j^3+2d_i^2{\mathrm{\λ}}_j^2\\ +2{\mathrm{\λ}}_i^2{\mathrm{\λ}}_j^2-4d_i^3{\mathrm{\λ}}_i-6d_i^2{\mathrm{\λ}}_i^2-4d_i{\mathrm{\λ}}_i^3\\ -{\mathrm{\λ}}_i^4-{\mathrm{\λ}}_j^4+16d_j{\mathrm{\λ}}_j+8{\mathrm{\λ}}_j^2+4d_j^2{d}_i{\mathrm{\λ}}_i\\ +4d_i^2{d}_j{\mathrm{\λ}}_j+4{\mathrm{\λ}}_i^2{d}_j{\mathrm{\λ}}_j+4{\mathrm{\λ}}_j^2{d}_i{\mathrm{\λ}}_i{)}^{1/2}\end{array}---\left(19\right)$

由于误差相对于距离而言是很小的，因此误差的平方项作差更小，可以忽略，由(18)和(19)式可得

$\begin{array}{c}x=\frac{1}{4}\[(d_i^2+2d_i{\mathrm{\λ}}_i+{\mathrm{\λ}}_i^2)-(d_j^2+2d_j{\mathrm{\λ}}_j+{\mathrm{\λ}}_j^2)\]\\ =\frac{1}{4}(d_i^2-d_j^2)+\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)+\frac{1}{4}({\mathrm{\λ}}_i^2-{\mathrm{\λ}}_j^2)\\ \≅\frac{1}{4}(d_i^2-d_j^2)+\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)\\ =const+\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)\end{array}---\left(3\right)$

其中 $const=\frac{1}{4}(d_i^2-d_j^2)$ 表示常数，

则坐标的误差为：

$\begin{array}{c}x\_error=\frac{1}{2}(d_i{\mathrm{\λ}}_i-d_j{\mathrm{\λ}}_j)\\ =\frac{1}{2}(d_i^2{\mathrm{\ρ}}_i-d_j^2{\mathrm{\ρ}}_j)\end{array}---\left(4\right)$

同理可得

$\begin{array}{c}y\_error=\frac{1}{4}\left(\frac{{\mathrm{\β}}_1}{2\sqrt{\mathrm{\α}}}{\mathrm{\λ}}_i+\frac{{\mathrm{\β}}_2}{2\sqrt{\mathrm{\α}}}{\mathrm{\λ}}_j\right)\\ =\frac{1}{4}\left(\frac{d_i{\mathrm{\β}}_1}{2\sqrt{\mathrm{\α}}}{\mathrm{\ρ}}_i+\frac{d_i{\mathrm{\β}}_2}{2\sqrt{\mathrm{\α}}}{\mathrm{\ρ}}_j\right)\end{array}---\left(5\right)$

其中x_error,y_error分别表示坐标x,y的误差，且令

$\left\{\begin{array}{c}\mathrm{\α}=-16+8d_j^2+8d_i^2-d_j^4-d_i^4+2d_i^2{d}_j^2\\ {\mathrm{\β}}_1=16d_i-4d_i^3+4d_j^3{d}_i\\ {\mathrm{\β}}_2=16d_j-4d_j^3+4d_i^3{d}_j\end{array}\right.,$

步骤一(三)：通过测量得到误差概率密度函数；

(1)无遮挡的环境下服从正态分布： $f\left({\mathrm{\ρ}}_i\right)=\frac{1}{\mathrm{\σ}\sqrt{2\mathrm{\π}}}{e}^{-\frac{{\mathrm{\ρ}}_i^2}{2{\mathrm{\σ}}^2}},{\mathrm{\ρ}}_i\∈(-\mathrm{\δ},\mathrm{\δ})---\left(7\right)$

(2)有遮挡的环境下服从瑞利分布： $f\left({\mathrm{\ρ}}_i\right)=\frac{{\mathrm{\ρ}}_i}{{\mathrm{\σ}}^2}\exp (-\frac{{\mathrm{\ρ}}_i^2}{2{\mathrm{\σ}}^2}),{\mathrm{\ρ}}_i\≥0---\left(8\right)$

(3)空旷环境下服从均匀分布： $f\left({\mathrm{\ρ}}_i\right)=\frac{1}{2\mathrm{\δ}},{\mathrm{\ρ}}_i\∈(-\mathrm{\δ},\mathrm{\δ})---\left(9\right)$

其中f(ρ_i)表示第i个测量不确定度ρ_i的概率密度函数,i＝1,2，σ表示标准差，δ表示误差范围；

步骤一(四)：利用步骤一(三)所得函数，计算联合概率分布：

每个锚节点到未知节点的距离是独立测量的，联合概率分布函数f为

$f=\underset{i=1}{\overset{2}{\Π}}f\left({\mathrm{\ρ}}_i\right)---\left(10\right)$

步骤一(五)：计算测量误差的期望值：

$\left\{\begin{array}{c}E(x\_error)=\∫\∫x\_error\×{\mathrm{fd\ρ}}_1{\mathrm{d\ρ}}_2\\ E(y\_error)=\∫\∫y\_error\×{\mathrm{fd\ρ}}_1{\mathrm{d\ρ}}_2\end{array}\right.---\left(11\right)$

x_error,y_error分别表示坐标x,y的误差，E(x_error)，E(y_error)分别表示横纵坐标的误差期望值；

步骤一(六)：利用其他组合的锚节点，重复步骤一(四)和步骤一(五)，得到M组误差期望值，之后求平均值，即

$\left\{\begin{array}{c}\stackrel{\‾}{x\_error}=\frac{1}{M}\underset{j=1}{\overset{M}{\Σ}}{E}_j(x\_error)\\ \stackrel{\‾}{y\_error}=\frac{1}{M}\underset{j=1}{\overset{M}{\Σ}}{E}_j(y\_error)\end{array}\right.---\left(12\right)$

其中 分别表示横纵坐标测量误差的平均值；

步骤一(七)：利用未知节点通信半径内三个等效锚节点坐标(x_i,y_i；x_j,y_j；x_k,y_k)估算未知节点的粗略坐标(x_estimate,y_estimate)的均值；

步骤一(八)：求最终坐标(x_end,y_end)：

$\left\{\begin{array}{c}x\_end=\stackrel{\‾}{x\_estimate}-\stackrel{\‾}{x\_error}\\ y\_end=\stackrel{\‾}{y\_estimate}-\stackrel{\‾}{y\_error}\end{array}\right.---\left(15\right)$

和为未知节点的粗略坐标(x_estimate,y_estimate)的均值；

步骤一(九)精度评估：利用指标均方根距离来衡量定位精度，其表达式为

rmsd＝(||x_end-x_true||²+||y_end-y_true||²)^1/2(16)

其中(x_true,y_true)表示未知节点的真实位置坐标，rmsd为指标均方根距离；

步骤二：利用基于误差校正的联合质心算法完成剩余节点的定位，具体过程为：

步骤二(一)：选择部分已定位出的未知节点，作为准锚节点，和已有的锚节点共同作为等效锚节点，参与其他未知节点的二次定位；

步骤二(二)：通过接收信号强度测算出未知节点与等效锚节点之间的距离；

步骤二(三)：利用未知节点通信半径内的三个等效锚节点坐标，算出未知节点的粗略坐标： $\left\{\begin{array}{c}x\_estimate=\frac{a_0{a}_1-a_1{a}_2}{a_3{b}_1-a_2{b}_2}\\ y\_estimate=\frac{a_0{b}_2-a_1{b}_1}{a_2{b}_2-a_3{b}_1}\end{array}\right.---\left(13\right)$

其中

$\left\{\begin{array}{c}{a}_0=0.5\[(x_i^2-x_j^2)+(y_i^2-y_j^2)\]-0.5(m_i^2-m_j^2)\\ a_1=0.5\[(x_j^2-x_k^2)+(y_j^2-y_k^2)\]-0.5(m_j^2-m_k^2)\\ a_2=y_i-y_j\\ a_3=y_j-y_k\\ b_1=x_i-x_j\\ b_2=x_j-x_k\end{array}\right.$

步骤二(四)：利用其它等效锚节点的重新组合，重复步骤二(三)并获得N组坐标粗略估算值，记作(x_eatimate,y_eatimate)_i,i＝1,2,…N，其中q表示未知节点通信半径内的锚节点和第一次定位选出的准锚节点的总和；

步骤二(五)：从N组估算出的位置点中选取G个点，其中3≤G≤N，选出的G个点构成一个点数最多的凸区域；

步骤二(六)：计算凸区域的质心作为未知节点的最终的坐标：

$\left\{\begin{array}{c}x\_end=\frac{1}{G}\underset{i=1}{\overset{G}{\Σ}}{(x\_eatimate)}_i\\ y\_end=\frac{1}{G}\underset{i=1}{\overset{G}{\Σ}}{(y\_eatimate)}_i\end{array}\right.---\left(17\right);$

步骤二(七)：重复步骤二(二)至步骤二(六)，直到完成网络中剩余未知节点的定位。

2.根据权利要求1所述的一种适用于传感器网络节点定位的基于误差校正的联合质心定位方法，其特征在于步骤一(七)所述的估算未知节点的粗略坐标(x_estimate,y_estimate)的均值过程为：根据式(13)进行 $M^{\′}=\left(\begin{array}{c}n\\ 3\end{array}\right)$ 次估算，可得估算均值为

$\left\{\begin{array}{c}\stackrel{\‾}{x\_estimate}=\frac{1}{M^{\′}}\underset{j=1}{\overset{M^{\′}}{\Σ}}{(x\_estimate)}_j\\ \stackrel{\‾}{y\_estimate}=\frac{1}{M^{\′}}\underset{j=1}{\overset{M^{\′}}{\Σ}}{(y\_estimate)}_j\end{array}\right.---\left(14\right).$