发明名称 一种无线传感器网络节点定位方法及装置
摘要 本发明公开了一种无线传感器网络定位方法,所述方法包括:A:在预设区域内随机产生一个坐标,将该坐标作为该区域内所有未知节点的初始估计坐标,并将初始估计坐标形成一个矩阵作为初始解;B:将初始解作为旧解,根据目标函数生成旧解对应的旧目标函数值;C:根据旧解生成一个新解,并生成新解对应的新目标函数值;D:比较新目标函数值是否比旧目标函数值小,如果是,则将新解作为旧解,并将该新目标函数值作为旧目标函数值;E:判断本次循环是否符合第一预设条件,如果是,则返回C;如果否,则结束循环,并将D中目标函数值较小的解作为当前局部最优解,根据当前局部最优解确定未知节点的坐标,以实现无线传感器网络定位。
申请公布号 CN103517411B 申请公布日期 2017.02.15
申请号 CN201310375549.4 申请日期 2013.08.22
申请人 东北大学 发明人 闻英友;陈锡民;赵大哲;孟颖辉;陈剑
分类号 H04W64/00(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W64/00(2009.01)I
代理机构 北京集佳知识产权代理有限公司 11227 代理人 黄熊;王宝筠
主权项 一种无线传感器网络定位方法,其特征在于,所述方法包括:步骤A:在预设区域内随机产生一个坐标,将该坐标作为该区域内所有未知节点的初始估计坐标,将所述所有未知节点的初始估计坐标形成一个矩阵作为初始解;步骤B:将所述初始解作为旧解,根据目标函数生成旧解对应的旧目标函数值;步骤C:根据所述旧解生成一个新解,并根据目标函数生成所述新解对应的新目标函数值;步骤D:比较新目标函数值是否比旧目标函数值小,如果是,则将所述新解作为旧解,并将该新目标函数值作为旧目标函数值;步骤E:判断本次循环是否符合第一预设条件,如果是,则返回步骤C;如果否,则结束循环,并将步骤D中目标函数值较小的解作为当前局部最优解,根据当前局部最优解确定所述区域内未知节点的坐标,以实现无线传感器网络定位;所述区域共有n个节点,其中m个锚节点,n‑m个未知节点,所述目标函数由CF表示,<maths num="0001"><math><![CDATA[<mrow><mi>C</mi><mi>F</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mi>m</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><msub><mi>M</mi><mi>i</mi></msub></mrow></munder><msup><mrow><mo>(</mo><msub><mover><mi>d</mi><mo>~</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA0001136930370000011.GIF" wi="501" he="149" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0001136930370000012.GIF" wi="630" he="87" /></maths>其中,M<sub>i</sub>为未知节点i的邻居节点组合,<img file="FDA0001136930370000013.GIF" wi="53" he="83" />为邻居节点i和j之间的实际距离,所述邻居节点之间的实际距离<img file="FDA0001136930370000014.GIF" wi="52" he="79" />具体由如下步骤获得:获取邻居节点i和j之间的接近度,所述接近度是指邻居节点i和j之间的接近程度;根据邻居节点i和j之间的接近度与标准距离获得邻居节点i和j之间的估计距离,将所述估计距离作为邻居节点i和j之间的实际距离<img file="FDA0001136930370000015.GIF" wi="79" he="82" />所述标准距离是指接近度为1的实际距离;所述<img file="FDA0001136930370000016.GIF" wi="54" he="77" />为根据邻居节点i和j的估计坐标计算的邻居节点i和j之间的估计距离,<img file="FDA0001136930370000017.GIF" wi="141" he="55" />和<img file="FDA0001136930370000018.GIF" wi="154" he="63" />分别表示邻居节点i和j的估计坐标,所述邻居节点i和j的估计坐标为旧解中邻居节点i和j的坐标,如果邻居节点是锚节点,则其估计坐标为实际坐标;所述根据所述旧解生成一个新解具体为:<maths num="0003"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msup><mover><mi>U</mi><mo>&OverBar;</mo></mover><mrow><mi>n</mi><mi>e</mi><mi>w</mi></mrow></msup><mo>=</mo><mi>i</mi><mi>n</mi><mi>v</mi><mrow><mo>(</mo><mi>N</mi><mn>1</mn><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>-</mo><mi>m</mi></mrow><mo>)</mo><mo>+</mo><mi>N</mi><mn>2</mn><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>-</mo><mi>m</mi></mrow><mo>)</mo><mo>)</mo></mrow><mo>.</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>&lsqb;</mo><mo>(</mo><mi>D</mi><mn>1</mn><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>-</mo><mi>m</mi></mrow><mo>)</mo></mrow><mo>+</mo><mi>D</mi><mn>2</mn><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>-</mo><mi>m</mi></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mover><mi>U</mi><mo>&OverBar;</mo></mover><mrow><mi>o</mi><mi>l</mi><mi>d</mi></mrow></msup><mo>+</mo><mrow><mo>(</mo><mrow><mi>N</mi><mi>A</mi><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>m</mi></mrow><mo>)</mo></mrow><mo>-</mo><mi>D</mi><mn>3</mn><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>m</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>A</mi><mo>&rsqb;</mo><mo>,</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001136930370000021.GIF" wi="1558" he="199" /></maths>其中,<img file="FDA0001136930370000022.GIF" wi="110" he="55" />和<img file="FDA0001136930370000023.GIF" wi="90" he="55" />分别表示未知节点在新、旧解中的估计坐标矩阵;矩阵A表示所述区域内所有锚节点的坐标;NA(n‑m,m)矩阵表示未知节点和锚节点之间的真实邻居关系,如果未知节点i和锚节点s是邻居节点,则对应的na<sub>is</sub>值为1,否则为0;其余参数矩阵N1(n‑m,n‑m)、N2(n‑m,n‑m)、D1(n‑m,n‑m)、D2(n‑m,n‑m)和D3(n‑m,m)的计算方法如下:<maths num="0004"><math><![CDATA[<mrow><mi>N</mi><mn>1</mn><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mi>m</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>nu</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mi>j</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>nu</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mtd><mtd><mrow><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>N</mi><mn>2</mn><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mo>{</mo><mrow><mtable><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>na</mi><mrow><mi>i</mi><mi>s</mi></mrow></msub></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mi>j</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></mtd></mtr></mtable><mo>,</mo></mrow></mrow>]]></math><img file="FDA0001136930370000024.GIF" wi="1542" he="215" /></maths><img file="FDA0001136930370000025.GIF" wi="1141" he="472" /><img file="FDA0001136930370000026.GIF" wi="1063" he="223" /><img file="FDA0001136930370000027.GIF" wi="779" he="213" />其中,nu<sub>ij</sub>是指未知节点i与未知节点j之间的邻居关系;<img file="FDA0001136930370000028.GIF" wi="46" he="62" />是指邻居节点i和j之间的实际距离,<img file="FDA0001136930370000029.GIF" wi="69" he="67" />是指根据当前未知节点i的估计坐标和未知节点j估计坐标计算的估计距离;<img file="FDA00011369303700000210.GIF" wi="53" he="71" />是指未知节点i和锚节点s之间的实际距离,<img file="FDA00011369303700000211.GIF" wi="81" he="70" />是指根据当前未知节点i的估计坐标和锚节点s估计坐标计算的估计距离。
地址 110003 辽宁省沈阳市和平区文化路3号巷11号