发明名称 无线传感器网络静态节点最优覆盖控制集求解方法
摘要 本发明公布了一种无线传感器网络静态节点最优覆盖控制集求解方法,本发明在保持一定覆盖率的前提下,令一部分节点进入休眠状态,从而减少工作节点数,降低网络能耗,延长网络生存时间。进一步地,不考虑节点间通信和计算能量消耗,仅考虑节点传感模块的能量消耗,利用节点具有可调节传感半径的属性,寻找一组工作节点数少,同时传感半径组合合适,能耗低,并且满足一定覆盖要求的节点集合。建立了网络节点传感模块能量消耗、网络覆盖率,以及网络节点休眠率三个子函数,由这3个子函数共同构成遗传算法的适应度函数,采用遗传算法求解该模型。获得了一个覆盖范围大、工作节点数少和能量消耗低的最优工作节点组合和工作节点的最优传感半径组合。
申请公布号 CN102761881A 申请公布日期 2012.10.31
申请号 CN201210203418.3 申请日期 2012.06.19
申请人 江苏科技大学 发明人 朱志宇;张冰;李阳;伍雪冬;王建华;冯友兵;王敏;杨官校;戴晓强
分类号 H04W16/18(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W16/18(2009.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 楼高潮
主权项 1.一种无线传感器网络静态节点最优覆盖控制集求解方法,其特征在于所述方法如下:无线传感器网络中的节点n<sub>i</sub>的坐标为(x<sub>i</sub>,y<sub>i</sub>),二维平面监测区域中任意点p的坐标为(x<sub>p</sub>,y<sub>p</sub>),则节点n<sub>i</sub>对目标点p的检测概率为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>C</mi><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mfenced open='' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><mi>ifd</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>R</mi><mi>s</mi></msub><mo>-</mo><msub><mi>R</mi><mi>e</mi></msub></mtd></mtr><mtr><mtd><msup><mi>e</mi><mrow><mrow><mo>(</mo><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><msubsup><mi>&alpha;</mi><mn>1</mn><msub><mi>&beta;</mi><mn>1</mn></msub></msubsup><mo>/</mo><mrow><mo>(</mo><msubsup><mi>&alpha;</mi><mn>2</mn><msub><mi>&beta;</mi><mn>2</mn></msub></msubsup><mo>+</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow></msup></mtd><mtd><mi>if</mi><msub><mi>R</mi><mi>s</mi></msub><mo>-</mo><msub><mi>R</mi><mi>e</mi></msub><mo>&lt;</mo><mi>d</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>R</mi><mi>s</mi></msub><mo>+</mo><msub><mi>R</mi><mi>e</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><mi>otherwise</mi></mtd></mtr></mtable></mfenced></mfenced></mrow></math>]]></maths>其中,d(n<sub>i</sub>,p)为传感器节点n<sub>i</sub>与目标点p的欧式距离;R<sub>e</sub>(0&lt;R<sub>e</sub>&lt;R<sub>s</sub>)是传感节点测量可靠性参数;α<sub>1</sub>=R<sub>e</sub>-R<sub>s</sub>+d(n<sub>i</sub>,p),α<sub>2</sub>=R<sub>e</sub>+R<sub>s</sub>-d(n<sub>i</sub>,p);λ<sub>1</sub>、λ<sub>2</sub>、β<sub>1</sub>、β<sub>2</sub>是与传感节点特性有关的测量参数;因此,整个二维平面监测区域内的n个传感器节点对目标点p的联合检测概率为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>C</mi><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mi>all</mi></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><munder><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mi>n</mi></mrow></munder><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>C</mi><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,n<sub>all</sub>为测量目标点的传感器节点集合,联合检测概率不低于根据网络需求所设定的阈值c<sub>th</sub>,则目标点能被有效检测到的条件为:<maths num="0003"><![CDATA[<math><mrow><munder><mi>min</mi><mrow><msub><mi>x</mi><mi>p</mi></msub><mo>,</mo><msub><mi>y</mi><mi>p</mi></msub></mrow></munder><mo>{</mo><msub><mi>C</mi><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mi>all</mi></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>}</mo><mo>&GreaterEqual;</mo><msub><mi>C</mi><mi>th</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>将待测区域划分成m×n个网格,再将单元格简化为像素点,网络覆盖率定义为满足式(2)要求的单元格数量占总的单元格数量的比例,即:<maths num="0004"><![CDATA[<math><mrow><msub><mi>C</mi><mi>r</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><msub><mi>x</mi><mi>p</mi></msub><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><msub><mi>y</mi><mi>p</mi></msub><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>C</mi><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mi>all</mi></msub><mo>,</mo><mi>p</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>根据式(2)和式(3),可以得到网络覆盖率函数f<sub>1</sub>:<maths num="0005"><![CDATA[<math><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><msub><mi>x</mi><mi>p</mi></msub><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><msub><mi>y</mi><mi>p</mi></msub><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>C</mi><mi>p</mi></msub><mrow><mo>(</mo><msup><mi>S</mi><mo>&prime;</mo></msup><mo>,</mo><mi>p</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,监测区域划分成m×n个网格单元,监测区域内点p联合检测概率须大于概率阈值c<sub>th</sub>才能被计算在网络覆盖率内。节点休眠率函数f<sub>2</sub>:<maths num="0006"><![CDATA[<math><mrow><msub><mi>f</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mo>|</mo><mi>S</mi><mo>|</mo><mo>-</mo><mo>|</mo><msup><mi>S</mi><mo>&prime;</mo></msup><mo>|</mo></mrow><mrow><mo>|</mo><mi>S</mi><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中|S′|表示工作的节点数,|S|为布设的总的节点数;在监测区域内布设的节点中,工作节点集合为N={n<sub>1</sub>,n<sub>2</sub>,...n<sub>n</sub>},节点传感半径为<img file="FDA00001783791700021.GIF" wi="412" he="74" />其中<img file="FDA00001783791700022.GIF" wi="75" he="66" />表示节点n<sub>i</sub>的传感半径,且<img file="FDA00001783791700023.GIF" wi="347" he="61" /><img file="FDA00001783791700024.GIF" wi="58" he="67" />表示节点n<sub>i</sub>的覆盖面积,<img file="FDA00001783791700025.GIF" wi="64" he="66" />表示节点n<sub>i</sub>的能量消耗;工作节点的覆盖面积为所有工作节点覆盖面积的并集为:<maths num="0007"><![CDATA[<math><mrow><msub><mi>S</mi><msub><mi>n</mi><mn>1</mn></msub></msub><mo>&cup;</mo><msub><mi>S</mi><msub><mi>n</mi><mn>2</mn></msub></msub><mo>&cup;</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>&cup;</mo><msub><mi>S</mi><msub><mi>n</mi><mi>n</mi></msub></msub><mo>=</mo><munderover><mrow><mi></mi><mo>&cup;</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>S</mi><msub><mi>n</mi><mi>i</mi></msub></msub></mrow></math>]]></maths>覆盖区域的平均能耗为:<maths num="0008"><![CDATA[<math><mrow><mover><mi>E</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msub><mi>E</mi><msub><mi>n</mi><mn>1</mn></msub></msub><mo>+</mo><msub><mi>E</mi><msub><mi>n</mi><mn>2</mn></msub></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>E</mi><msub><mi>n</mi><mi>n</mi></msub></msub></mrow><mrow><msub><mi>S</mi><msub><mi>n</mi><mn>1</mn></msub></msub><mo>&cup;</mo><msub><mi>S</mi><msub><mi>n</mi><mn>2</mn></msub></msub><mo>&cup;</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>&cup;</mo><msub><mi>S</mi><msub><mi>n</mi><mi>n</mi></msub></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>R</mi><msub><mi>sn</mi><mi>i</mi></msub></msub></mrow><mrow><munderover><mo>&cup;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>S</mi><msub><mi>n</mi><mi>i</mi></msub></msub></mrow></mfrac><mi>u</mi></mrow></math>]]></maths>描述整个无线传感器网络的覆盖消耗的覆盖区域平均能耗子函数f<sub>3</sub>为:<maths num="0009"><![CDATA[<math><mrow><msub><mi>f</mi><mn>3</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>R</mi><msub><mi>sn</mi><mi>i</mi></msub><mn>2</mn></msubsup><mo>/</mo><mi>S</mi><mo>_</mo><mi>work</mi></mrow></math>]]></maths>其中,<img file="FDA00001783791700029.GIF" wi="71" he="63" />表示工作节点n<sub>i</sub>的传感半径,S_work表示工作节点集合所覆盖的区域面积;总体目标函数为:f=w<sub>1</sub>f<sub>1</sub>+w<sub>2</sub>f<sub>2</sub>-w<sub>3</sub>f<sub>3</sub>                                            (7)其中,w<sub>1</sub>+w<sub>2</sub>+w<sub>3</sub>=1,总体目标函数值f介于0~1之间,值越大,表明覆盖优化效果越优。
地址 212003 江苏省镇江市梦溪路2号