发明名称 一种提高能量有效性的无线传感器网络二阶一致性方法
摘要 本发明提供了一种提高能量有效性的无线传感器网络二阶一致性方法,包括对WSNs进行邻接权重分布式优化的步骤,二阶状态一致性模型建立及更新步骤,建立二阶状态一致性收敛条件步骤,二阶状态一致性能量消耗优化步骤。基于网络邻接权重的调整来优化二阶一致性算法的能耗,并采用分布式设计优化网络邻接矩阵,在收敛效率和通信能耗两方面均有良好改善;与现有其他方法相比,本发明有效地提高了收敛速度、降低了通信能耗,性能优越。
申请公布号 CN105611611A 申请公布日期 2016.05.25
申请号 CN201511022350.9 申请日期 2015.12.29
申请人 安徽理工大学 发明人 陈珍萍;黄友锐;曲立国;凌六一;唐超礼;徐善永;韩涛
分类号 H04W52/02(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W52/02(2009.01)I
代理机构 南京众联专利代理有限公司 32206 代理人 叶涓涓
主权项 一种提高能量有效性的无线传感器网络二阶一致性方法,其特征在于,包括如下步骤:步骤A,首先对WSNs进行邻接权重分布式优化:步骤A‑1,对矩阵W的第i,j个元素W<sub>ij</sub>按照下式进行初始化:<img file="FDA0000894762500000011.GIF" wi="1374" he="199" />上式中,d<sub>i</sub>、d<sub>j</sub>分别为节点i和节点j的入度;步骤A‑2,对网络G的拉普拉斯矩阵L作线性变换得到矩阵M=I‑L/2,初始化M(0)为I‑L/2,其中,L=D<sup>in</sup>‑W(0),W={W<sub>ij</sub>},<img file="FDA0000894762500000018.GIF" wi="309" he="63" /><img file="FDA0000894762500000012.GIF" wi="236" he="125" />步骤A‑3,采用下式所述的基于梯度的迭代算法进行M<sub>ij</sub>(k)和W<sub>ij</sub>(k)更新:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>M</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>P</mi><mi>X</mi></msub><mo>{</mo><msub><mi>M</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>&gamma;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>E</mi><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mi>X</mi></msub><mo>{</mo><mn>1</mn><mo>-</mo><munder><mo>&Sigma;</mo><mrow><mi>q</mi><mo>&Element;</mo><msub><mi>N</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></munder><msub><mi>M</mi><mrow><mi>i</mi><mi>q</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mi>j</mi></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>W</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mi>j</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>2</mn><msub><mi>M</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000894762500000013.GIF" wi="1494" he="325" /></maths>上式中,γ(k)是迭代步长且满足<img file="FDA0000894762500000014.GIF" wi="455" he="103" />g<sub>ij</sub>(k)是f(M)梯度函数g(k)的第i,j个元素,且g<sub>ij</sub>(k)按照下式计算:<maths num="0002"><math><![CDATA[<mrow><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>&part;</mo><mi>f</mi><mrow><mo>(</mo><mi>M</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><msub><mi>M</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mfrac><mo>=</mo><mn>2</mn><mrow><mo>(</mo><mn>2</mn><msub><mi>M</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mi>M</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mi>M</mi><mrow><mi>j</mi><mi>j</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>+</mo><mi>&eta;</mi><mi>sgn</mi><mrow><mo>(</mo><msub><mi>M</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000894762500000015.GIF" wi="1629" he="127" /></maths>上式中,sgn(x)为一符号函数,当x>0,sgn(x)取值为1;当x<0时,sgn(x)取值为‑1;当x=0时,sgn(x)取值为0;P<sub>X</sub>{y}为y到集合X的投影,X=[‑1,1],投影函数P<sub>X</sub>{y}具体表达为下式:<img file="FDA0000894762500000016.GIF" wi="1486" he="270" />上式中,ε为阈值;<img file="FDA0000894762500000017.GIF" wi="173" he="63" />为节点i的输入邻接集的子集,N<sub>ei</sub>表达为下式:N<sub>ei</sub>={j|j∈N<sub>i</sub>(k),d<sub>i</sub>(k)>2,d<sub>j</sub>(k)>2,j至少是i的第σ<sub>i</sub>(k)远节点}基于式(2)进行多次迭代,迭代过程如下:在进行第k次迭代时,若节点i的某非零M<sub>ij</sub>(k)被置为零,则调整σ<sub>i</sub>(k)=σ<sub>i</sub>(k)‑1、d<sub>i</sub>(k)=d<sub>i</sub>(k)‑1和d<sub>j</sub>(k)=d<sub>j</sub>(k)‑1,且σ<sub>i</sub>(0)=κd<sub>i</sub>(0)、d<sub>i</sub>(0)=d<sub>i</sub>,κ为调整网络G稀疏度的参数,d<sub>i</sub>为网络G中节点i的初始入度,σ<sub>i</sub>(k)为第k次迭代时节点i可以被优化的邻居节点数;在第k次迭代后,节点i调整σ<sub>i</sub>(k+1)=σ<sub>i</sub>(k)、d<sub>i</sub>(k+1)=d<sub>i</sub>(k)和d<sub>j</sub>(k+1)=d<sub>j</sub>(k),且根据剩余邻接节点中与节点i最远的节点来调整发射功率;迭代多次后,得到邻接权重优化矩阵W<sup>o</sup>;步骤B,二阶状态一致性模型建立及更新:网络G中,按照下式更新节点i的位置信息x<sub>1i</sub>(k)和速度信息x<sub>2i</sub>(k):<maths num="0003"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mi>T</mi></mtd><mtd><mrow><mn>0.5</mn><msup><mi>T</mi><mn>2</mn></msup></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>u</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000894762500000021.GIF" wi="1502" he="142" /></maths>上式中,T为网络更新周期,u<sub>1i</sub>(k)和u<sub>2i</sub>(k)为其位置和速度的一致性控制输入,所述u<sub>1i</sub>(k)和u<sub>2i</sub>(k)在第k次状态更新时由节点i根据自身与邻居节点的位置状态不一致构造,具体表示如下式:<maths num="0004"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>u</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr></mtable></mfenced><mo>=</mo><mo>-</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mi>&mu;</mi></mtd></mtr><mtr><mtd><msup><mi>&mu;</mi><mi>2</mi></msup></mtd></mtr></mtable></mfenced><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>W</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>o</mi></msubsup><mrow><mo>(</mo><mo>(</mo><msub><mi>x</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mi>x</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000894762500000022.GIF" wi="1325" he="143" /></maths>上式中,<img file="FDA0000894762500000023.GIF" wi="69" he="71" />即为步骤A得到的邻接权重优化矩阵,μ为控制增益,μ符合以下条件:<maths num="0005"><math><![CDATA[<mrow><mn>0</mn><mo>&lt;</mo><mi>&mu;</mi><mi>T</mi><mo>&lt;</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><mn>2</mn><mo>,</mo><mfrac><mn>2</mn><mrow><msub><mi>&lambda;</mi><mi>N</mi></msub><mrow><mo>(</mo><mi>L</mi><mo>)</mo></mrow></mrow></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000894762500000024.GIF" wi="1238" he="119" /></maths>在第k次状态更新后,节点i将自身位置信息x<sub>1i</sub>(k+1)广播给邻居节点。
地址 232002 安徽省淮南市舜耕中路168号