发明名称 一种无线传感器网络低能耗覆盖优化方法
摘要 本发明公开了一种无线传感器网络低能耗覆盖优化方法,其在保证轮次簇能耗最低的前提下,能够获取满足覆盖能力的骨干节点最优值,同时通过PSO动态调整节点的飞行方向以及飞行速度,逐步迭代获取最优覆盖的骨干节点部署方案,既能有效地采集感知区域的数据信息,又能够充分管理传感器网络资源,有效地解决了无线传感器网络应用中的骨干节点分布不均、骨干节点数目无法达到最优、网络能耗不均衡、网络覆盖能力偏低等问题,为无线传感器网络在保证低能耗的前提下合理有效部署骨干节点以达到较高的网络覆盖均匀性和网络覆盖率提供技术支撑,为在网络低能耗的前提下提高网络的区域覆盖能力提供了一种全新的解决方案。
申请公布号 CN103118373B 申请公布日期 2016.08.10
申请号 CN201310058571.6 申请日期 2013.01.24
申请人 北京理工大学 发明人 何遵文;陈存香;刘阳;匡镜明
分类号 H04W16/18(2009.01)I;H04W52/02(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W16/18(2009.01)I
代理机构 北京中海智圣知识产权代理有限公司 11282 代理人 杨树芬
主权项 一种无线传感器网络低能耗覆盖优化方法,其特征在于,包括以下步骤:1)根据基于保证轮次簇能耗及覆盖能力的骨干节点数目优化算法,获得骨干节点最优值:11)普通节点的能耗来自于感知数据的发送能量,骨干节点的能耗主要来自于接收感知数据、数据融合处理、数据转发至汇聚节点的能耗,因此普通节点的能耗为:<img file="FDA0000962015340000011.GIF" wi="384" he="68" />骨干节点的能耗为:<img file="FDA0000962015340000012.GIF" wi="849" he="62" />则某轮簇内的能耗为:E<sub>C</sub>=E<sub>CH</sub>+(N/k‑1)E<sub>M</sub>,网络总能耗为<img file="FDA0000962015340000013.GIF" wi="826" he="70" />式中,E<sub>elec</sub>为收发电路单位bit数据能耗;ε<sub>fs</sub>、ε<sub>mp</sub>分别为近距离和远距离的功率衰减系数;d<sub>toCH</sub>为各普通节点到骨干节点距离的均值;E<sub>DA</sub>为融合单位bit数据消耗的能量;d<sub>toSink</sub>为各骨干节点到汇聚节点距离的均值;d<sub>0</sub>为参考距离,一般为<img file="FDA0000962015340000014.GIF" wi="267" he="71" />其中,k为能耗最小时的骨干节点;N为节点集;l为监测区域内理想骨干节点部署信息;α为权重系数;ε<sub>fs</sub>、ε<sub>mp</sub>分别为近距离和远距离的功率衰减系数;12)假设监测区域边长为M,簇所占区域为圆形,骨干节点位于簇中心,根据节点覆盖最优模型,即相邻任意三个骨干节点成等边三角形,由此可得到每个簇所占区域的面积近似为<img file="FDA0000962015340000015.GIF" wi="251" he="63" />节点的分布律为<img file="FDA0000962015340000016.GIF" wi="251" he="69" />则普通节点到骨干节点距离平方的数学期望值为:<img file="FDA0000962015340000017.GIF" wi="941" he="79" />13)根据轮次内网络总能耗及普通节点到骨干节点距离平方的数学期望值,可获得能耗最小时骨干节点k的最优值为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>k</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><mroot><mn>27</mn><mn>4</mn></mroot><msqrt><mi>N</mi></msqrt><mi>M</mi></mrow><mrow><mn>2</mn><msubsup><mi>&pi;d</mi><mrow><mi>t</mi><mi>o</mi><mi>S</mi><mi>i</mi><mi>n</mi><mi>k</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><msqrt><mfrac><msub><mi>&epsiv;</mi><mrow><mi>f</mi><mi>s</mi></mrow></msub><msub><mi>&epsiv;</mi><mrow><mi>m</mi><mi>p</mi></mrow></msub></mfrac></msqrt><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>t</mi><mi>o</mi><mi>S</mi><mi>i</mi><mi>n</mi><mi>k</mi></mrow></msub><mo>&GreaterEqual;</mo><msub><mi>d</mi><mn>0</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><mroot><mn>27</mn><mn>4</mn></mroot><msqrt><mi>N</mi></msqrt><mi>M</mi></mrow><mrow><mn>2</mn><msubsup><mi>&pi;d</mi><mrow><mi>t</mi><mi>o</mi><mi>S</mi><mi>i</mi><mi>n</mi><mi>k</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>t</mi><mi>o</mi><mi>S</mi><mi>i</mi><mi>n</mi><mi>k</mi></mrow></msub><mo>&lt;</mo><msub><mi>d</mi><mn>0</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000962015340000018.GIF" wi="638" he="278" /></maths>ε<sub>fs</sub>、ε<sub>mp</sub>分别为近距离和远距离的功率衰减系数;d<sub>0</sub>为参考距离,一般为<img file="FDA0000962015340000019.GIF" wi="267" he="70" />2)根据用于保证网络覆盖均匀性的改进分簇部署算法,获得覆盖均匀性指标:21)覆盖均匀性反映了传感器节点在被监测区域的分布情况,根据覆盖均匀性与节点之间距离的关系,采用改进分簇部署算法来获取覆盖均匀性指标,距离标准差值越小,则覆盖均匀性越高,同时网络中节点的能量消耗越低,具体如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>U</mi><mo>=</mo><mfrac><mn>1</mn><msub><mi>k</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub></munderover><msub><mi>U</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA00009620153400000110.GIF" wi="239" he="125" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>U</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>k</mi><mi>i</mi></msub></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mi>i</mi></msub></munderover><msup><mrow><mo>(</mo><mrow><msub><mi>D</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>M</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mrow>]]></math><img file="FDA00009620153400000111.GIF" wi="446" he="150" /></maths>式中U为覆盖均匀性指标,即距离标准差值指标,k<sub>opt</sub>为骨干节点总数目,k<sub>i</sub>为第i个骨干节点的邻居节点个数,D<sub>i,j</sub>为第i个骨干节点与第j个骨干节点之间的距离,M<sub>i</sub>表示第i个骨干节点与邻居骨干节点距离的平均值;22)邻居节点的选取,根据节点覆盖最优模型以及最优骨干节点数量可得到骨干节点之间的最佳距离为<img file="FDA0000962015340000021.GIF" wi="489" he="84" />为了保证网络覆盖率,选取距离第i个骨干节点为(0,Ropt±2(Rs‑Ropt/2)]的节点为该骨干节点的邻居节点,R<sub>s</sub>为骨干节点的感知距离;3)根据用于解决骨干节点分布不均,平衡网络覆盖率和覆盖均匀性的改进PSO覆盖能效优化算法,获得骨干节点部署结果:31)优化目标是满足网络覆盖率最大化及覆盖均匀性最大化,即网络覆盖率最大化、距离标准差值最小化,由此可得到最小化目标函数为:f(X)=C‑<sup>(1‑α)</sup>U<sup>α</sup>,式中C和U分别为节点集N所对应网络部署状态的网络覆盖率和网络覆盖均匀性指标,α为权重系数,用于调节优化中两项指标的权重,以适应无线传感器网络不同的约束条件;32)初始化微粒群,根据各个微粒的位置信息,在节点集N中搜索距离各个微粒中骨干节点最近的节点,并作为初始微粒群集X;33)根据节点覆盖最优模型,获得监测区域内理想骨干节点部署信息L,根据下式将骨干节点进行迭代进化:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>&omega;v</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>r</mi><mn>1</mn></msub><mrow><mo>(</mo><msubsup><mi>p</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><msub><mi>r</mi><mn>2</mn></msub><mrow><mo>(</mo><msubsup><mi>p</mi><mrow><mi>g</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000962015340000022.GIF" wi="701" he="71" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>,</mo><mrow><mo>(</mo><mi>V</mi><mi> </mi><mi>m</mi><mi>a</mi><mi>x</mi><mo>&Subset;</mo><mo>&lsqb;</mo><mi>R</mi><mi>o</mi><mi>p</mi><mi>t</mi><mo>-</mo><mn>2</mn><mo>(</mo><mi>R</mi><mi>s</mi><mo>-</mo><mi>R</mi><mi>o</mi><mi>p</mi><mi>t</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mi>R</mi><mi>o</mi><mi>p</mi><mi>t</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><mi>R</mi><mi>s</mi><mo>-</mo><mi>R</mi><mi>o</mi><mi>p</mi><mi>t</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>&rsqb;</mo><mo>)</mo><mo>,</mo></mrow>]]></math><img file="FDA0000962015340000023.GIF" wi="1274" he="71" /></maths>每代进化完毕后获得具有最小适应度值的骨干节点集位置信息S,在S±2(Rs‑Ropt/2)范围内搜索是否满足节点集N,若满足则更新相应骨干节点位置信息,否则设置λ(S±L)为该节点下次迭代的飞行方向。
地址 100081 北京市海淀区中关村南大街5号