发明名称 一种基于RSS和测距无偏估计的WSN节点定位方法
摘要 一种基于RSS和测距无偏估计的WSN节点定位方法,本发明涉及节点定位方法;本发明是要解决基于RSS测距的WSN节点定位算法受噪声的影响而引起测距误差的分析不充分,不能从本质上消除由于噪声所引起的测距偏差的问题而提出的一种基于RSS和测距无偏估计的WSN节点定位方法。该方法是通过1、设未知节点U=(x,y),B<sub>i</sub>信号的坐标为x<sub>i</sub>,y<sub>i</sub>;2、计算未知节点与信标节点B<sub>i</sub>的距离v<sub>i</sub>;3、得到无偏估计量<img file="DDA0000657512330000011.GIF" wi="80" he="88" />4、计算未知节点U到信标节点B<sub>i</sub>的距离<img file="DDA0000657512330000012.GIF" wi="72" he="73" />5、利用牛顿迭代法求解未知节点U的位置x,y;6、获得△X7、得到未知节点定位位置x<sub>k</sub>,y<sub>k</sub>等步骤实现的。本发明应用于节点定位领域。
申请公布号 CN104507164A 申请公布日期 2015.04.08
申请号 CN201510023979.9 申请日期 2015.01.16
申请人 哈尔滨工业大学 发明人 孟维晓;岳晋;韩帅;巩紫君;张德坤
分类号 H04W64/00(2009.01)I 主分类号 H04W64/00(2009.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种基于RSS和测距无偏估计的WSN节点定位方法,其特征在于一种基于RSS和测距无偏估计的WSN节点定位方法具体是按照以下步骤进行的:步骤一、在传感器网络工作环境内,预先部署M个信标节点;设传感器网络工作环境内任意的未知节点U=(x,y),设U接收到来自信标节点B<sub>i</sub>信号的坐标为x<sub>i</sub>,y<sub>i</sub>,i=1,2,...,N,N≤M;其中,N表示对于未知节点U可见的信标节点的数目;步骤二、根据一般的信号传播模型,计算出在受高斯噪声n的影响下,未知节点与信标节点B<sub>i</sub>的距离<img file="FDA0000657512300000011.GIF" wi="512" he="150" />其中,在距离信标节点d<sub>0</sub>处设置参考节点,参考节点接收信标节点B<sub>i</sub>的信号功率为P<sub>i</sub>(d<sub>0</sub>);未知节点U与信标节点B<sub>i</sub>的距离为d<sub>i</sub>;在不受噪声的情况下,未知节点U接收到信标节点B<sub>i</sub>的信号功率为P<sub>i</sub>(d<sub>i</sub>);n为均值为0,方差为<img file="FDA0000657512300000012.GIF" wi="68" he="86" />的高斯噪声;α为路径损耗指数;步骤三、对步骤二的v<sub>i</sub>进行期望值分析得到<img file="FDA0000657512300000013.GIF" wi="491" he="151" />为了使得<img file="FDA0000657512300000014.GIF" wi="242" he="106" />引入未知节点U到信标节点B<sub>i</sub>真实距离的无偏估计量<img file="FDA0000657512300000015.GIF" wi="451" he="155" />其中,<img file="FDA0000657512300000016.GIF" wi="562" he="124" />i=1,2,...,N,n为均值为0,方差为<img file="FDA0000657512300000017.GIF" wi="68" he="78" />的高斯噪声;步骤四、根据二维平面内未知节点U=(x,y)和信标节点B<sub>i</sub>=(x<sub>i</sub>,y<sub>i</sub>)的坐标,计算未知节点U到信标节点B<sub>i</sub>的距离<img file="FDA0000657512300000018.GIF" wi="54" he="82" />满足下式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>d</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000657512300000019.GIF" wi="673" he="97" /></maths>其中,i=1,2,...,N,N≥3;步骤五、利用牛顿迭代法求解未知节点U的位置x,y;将式(1)中在(x<sub>k‑1</sub>,y<sub>k‑1</sub>)<sup>T</sup>处进行泰勒级数展开,并略去二次项,得:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>d</mi><mo>^</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mi>r</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub></mrow><msub><mi>r</mi><mi>i</mi></msub></mfrac><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub></mrow><msub><mi>r</mi><mi>i</mi></msub></mfrac><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006575123000000110.GIF" wi="982" he="137" /></maths>其中,<img file="FDA00006575123000000111.GIF" wi="634" he="112" />i=1,2,...,N;x<sub>k</sub>,y<sub>k</sub>为第k次牛顿迭代法求解未知节点U位置坐标;步骤六、结合最小二乘法对式(2)进行求解,得到矩阵形式:G·ΔX=b  (3)其中,ΔX=[(x‑x<sub>k‑1</sub>) (y‑y<sub>k‑1</sub>)]<sup>T</sup>,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>G</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mfrac><mrow><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac></mtd><mtd><mfrac><mrow><msub><mi>y</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><mi>y</mi></mrow><msub><mi>r</mi><mn>1</mn></msub></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac></mtd><mtd><mfrac><mrow><msub><mi>y</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub></mrow><msub><mi>r</mi><mn>2</mn></msub></mfrac></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mi>N</mi></msub></mrow><msub><mi>r</mi><mi>N</mi></msub></mfrac></mtd><mtd><mfrac><mrow><msub><mi>y</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mi>N</mi></msub></mrow><msub><mi>r</mi><mi>N</mi></msub></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000657512300000021.GIF" wi="741" he="528" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>d</mi><mo>^</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mi>r</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>d</mi><mo>^</mo></mover><mn>2</mn></msub><mo>-</mo><msub><mi>r</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mover><mi>d</mi><mo>^</mo></mover><mi>N</mi></msub><mo>-</mo><msub><mi>r</mi><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000657512300000022.GIF" wi="518" he="344" /></maths>结合最小二乘理论,获得ΔX为:ΔX=(G<sup>T</sup>G)<sup>‑1</sup>G<sup>T</sup>b (6)步骤七、取ΔX=[(x<sub>k</sub>‑x<sub>k‑1</sub>) (y<sub>k</sub>‑y<sub>k‑1</sub>)]<sup>T</sup>,根据公式(4)、(5)和(6)进行迭代求解得到最终的未知节点定位位置x<sub>k</sub>,y<sub>k</sub>;即完成了一种基于RSS和测距无偏估计的WSN节点定位方法。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号