发明名称 一种无线传感器网络网关优化部署方法
摘要 一种无线传感器网络网关优化部署方法,涉及一种网关优化部署方法。为了解决几何K中心下的无线传感器网络网关部署问题,以缩小覆盖半径,提高网络服务质量。主要步骤:网关位置向量初始化、网关位置向量变异操作、交叉操作、选择操作、重复上述步骤直到迭代次数到达P=500,在第500代种群中分别计算各个目标向量<img file="DDA0000724768830000011.GIF" wi="369" he="74" />对应的适应值,适应值最小的一个目标向量即为无线传感器网络中网关的最优部署位置坐标。实验结果表明,通过微分进化算法求解无线传感器网络中网关的部署位置,能够比现有基于粒子群的算法收敛速度提高50%左右,覆盖半径缩小20%,因此该方法能够显著提高网络服务质量。
申请公布号 CN104883702A 申请公布日期 2015.09.02
申请号 CN201510274852.4 申请日期 2015.05.26
申请人 哈尔滨工业大学 发明人 杨京礼;许永辉;姜守达
分类号 H04W24/02(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W24/02(2009.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种无线传感器网络网关优化部署方法,在传感器节点数量为n的无线传感器网络G中设置K个网关节点,G的邻接矩阵为A=(e<sub>ij</sub>)<sub>n×n</sub>,最短距离矩阵为D=(d<sub>ij</sub>)<sub>n×n</sub>,d<sub>ij</sub>表示从节点v<sub>i</sub>路由到v<sub>j</sub>所需要的最小跳数,最短距离矩阵由Floyd算法求得;节点v<sub>i</sub>选择网关u<sub>k</sub>作为其服务网关,则满足式(1)的要求:d(v<sub>i</sub>,u<sub>k</sub>)≤d(v<sub>i</sub>,u<sub>l</sub>),k,l≤K,j≠l          (1)式中:v<sub>i</sub>包含在网关u<sub>k</sub>的服务集U<sub>k</sub>中,即v<sub>i</sub>∈U<sub>k</sub>,u<sub>k</sub>与服务集U<sub>k</sub>中节点之间的最大距离为<img file="FDA0000724768800000011.GIF" wi="355" he="107" />称为网关u<sub>k</sub>的覆盖半径;所有网关节点中的最大覆盖半径<img file="FDA0000724768800000012.GIF" wi="413" he="107" />称为网关集{u<sub>k</sub>}<sub>K</sub>的覆盖半径;所述无线传感器网络网关优化部署方法使得网关集的覆盖半径最小,如式(2)所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><mi>min</mi></mtd><mtd><mi></mi><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>K</mi></mrow></munder><munder><mi>max</mi><mrow><msub><mi>v</mi><mi>i</mi></msub><mo>&Element;</mo><msub><mi>U</mi><mi>k</mi></msub></mrow></munder><mrow><mo>(</mo><mi>d</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>u</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mtd><mtd><msub><mrow><mo>{</mo><msub><mi>u</mi><mi>k</mi></msub><mo>}</mo></mrow><mi>K</mi></msub><mo>&Subset;</mo><msup><mi>R</mi><mn>2</mn></msup></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000724768800000013.GIF" wi="1317" he="192" /></maths>使用微分进化算法进行上述式(2)的寻优求解,定义对于K个网关节点,其在二维平面的坐标为:(a<sub>k</sub>,b<sub>k</sub>),k=1,2,...,K,用网关坐标组成的目标向量为:X=(x<sub>1</sub>,x<sub>2</sub>,...,x<sub>m</sub>),其中m=2K,x<sub>2k‑1</sub>=a<sub>k</sub>,x<sub>2k</sub>=b<sub>k</sub>,第t组目标向量为X<sub>t</sub>=(x<sub>t,1</sub>,x<sub>t,2</sub>,...,x<sub>t,m</sub>);寻优求解的过程如下:步骤一、网关位置向量初始化在网络有效区域内,随机产生T组目标向量X<sub>1</sub>,X<sub>2</sub>,...,X<sub>T</sub>组成第一代种群,种群中的每个目标向量表示一组可能的网关位置坐标;设置交叉因子F=0.8,交叉概率为CR=0.4,最大迭代次数P=500,按照网络区域大小设置目标向量中各维数据的上下限范围[down_limit,up_limit];步骤二、网关位置向量变异操作对于由网关位置坐标组成的第p代种群任意一个目标向量<img file="FDA0000724768800000014.GIF" wi="481" he="93" />其中:p=1,2,...,P,t=1,2,...,T;根据微分进化算法按式(3)产生下一代网关位置坐标组成的变异向量<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><msub><mi>V</mi><mi>t</mi></msub><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>=</mo><mrow><mo>(</mo><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mn>1</mn></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>,</mo><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mn>2</mn></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>m</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000724768800000015.GIF" wi="535" he="100" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>x</mi><mrow><msub><mi>r</mi><mn>1</mn></msub><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>+</mo><mi>F</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mrow><msub><mi>r</mi><mn>2</mn></msub><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><msub><mi>r</mi><mn>3</mn></msub><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>)</mo></mrow><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000724768800000016.GIF" wi="1446" he="96" /></maths>其中,<img file="FDA0000724768800000021.GIF" wi="230" he="78" />和<img file="FDA0000724768800000022.GIF" wi="88" he="78" />为第p代种群中随机选择的3个个体目标向量第j位的元素,并且为3个不同的个体;交叉因子F是一个实数,用于控制差值的放大倍数,F=0.8;按照式(3)进行元素变异操作之后,再按照式(4)以初始化过程设定的目标向量中各维数据上下限依据,将进行元素变异操作之后超过有效界限的元素拉回到边界处:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>down</mi><mo>_</mo><mi>limit</mi></mtd><mtd><mi>if</mi></mtd><mtd><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&lt;</mo><mi>down</mi><mo>_</mo><mi>limit</mi></mtd></mtr><mtr><mtd><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd><mtd><mi>if</mi></mtd><mtd><mi>down</mi><mo>_</mo><mi>limit</mi><mo>&le;</mo><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><mi>up</mi><mo>_</mo><mi>limit</mi></mtd></mtr><mtr><mtd><mi>up</mi><mo>_</mo><mi>limit</mi></mtd><mtd><mi>if</mi></mtd><mtd><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>></mo><mi>up</mi><mo>_</mo><mi>limit</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000724768800000023.GIF" wi="1581" he="244" /></maths>步骤三、交叉操作在完成网关位置向量的变异操作后,产生下一代由网关位置坐标组成的交叉向量<img file="FDA0000724768800000024.GIF" wi="586" he="93" />其中向量的每位元素按照式(5)进行计算<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>&psi;</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>v</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd><mtd><mi>rand</mi><mo>&le;</mo><mi>CR</mi><mo>|</mo><mo>|</mo><mi>j</mi><mo>=</mo><msub><mi>j</mi><mi>rand</mi></msub></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd><mtd><mi>rand</mi><mo>></mo><mi>CR</mi><mo>|</mo><mo>|</mo><mi>j</mi><mo>&NotEqual;</mo><msub><mi>j</mi><mi>rand</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000724768800000025.GIF" wi="1533" he="179" /></maths>其中,rand为0~1之间的随机数,j<sub>rand</sub>为1~m中的随机整数;CR是交叉概率,CR=0.4;步骤四、选择操作d(v<sub>i</sub>,u<sub>k</sub>)为传感器节点v<sub>i</sub>到网关节点u<sub>k</sub>的跳数;对于第k个网关u<sub>k</sub>,其节点位置为(a<sub>k</sub>,b<sub>k</sub>),k=1,2,...,K;距离该网关节点距离小于通信半径的传感器节点组成的集合为Θ<sub>k</sub>,该网关节点到Θ<sub>k</sub>中任意节点的跳数均为1;节点v<sub>i</sub>到Θ<sub>k</sub>的距离为<img file="FDA0000724768800000026.GIF" wi="567" he="109" />则节点v<sub>i</sub>到网关节点u<sub>k</sub>的跳数可按式(6)计算:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>u</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>v</mi><mi>i</mi></msub><mo>&Element;</mo><msub><mi>&Theta;</mi><mi>k</mi></msub></mtd></mtr><mtr><mtd><munder><mi>min</mi><mrow><msub><mi>v</mi><mi>l</mi></msub><mo>&Element;</mo><msub><mi>&Theta;</mi><mi>k</mi></msub></mrow></munder><mrow><mo>(</mo><mi>d</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>v</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mtd><mtd><msub><mi>v</mi><mi>i</mi></msub><mo>&NotElement;</mo><msub><mi>&Theta;</mi><mi>k</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000724768800000027.GIF" wi="1472" he="181" /></maths>无线传感器网络网关优化部署的适应值计算函数为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><msub><mrow><mo>{</mo><msub><mi>u</mi><mi>k</mi></msub><mo>}</mo></mrow><mi>K</mi></msub><mo>)</mo></mrow><mo>=</mo><munder><mi>max</mi><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>n</mi></mrow></munder><munder><mi>min</mi><mrow><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>K</mi></mrow></munder><mi>d</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&Theta;</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000724768800000028.GIF" wi="1389" he="109" /></maths>在交叉操作完成后,按照式(8)进行选择操作:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>X</mi><mi>y</mi><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&Psi;</mi><mi>t</mi><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd><mtd><mi>f</mi><mrow><mo>(</mo><msubsup><mi>&Psi;</mi><mi>t</mi><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>&lt;</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mi>t</mi><mi>p</mi></msubsup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>X</mi><mi>t</mi><mi>p</mi></msubsup></mtd><mtd><mi>f</mi><mrow><mo>(</mo><msubsup><mi>&Psi;</mi><mi>t</mi><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mi>t</mi><mi>p</mi></msubsup><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000724768800000029.GIF" wi="1401" he="166" /></maths>上述选择过程是指如果新产生的个体向量由于父代中相应位置的个体向量,则将其取代父代中的个体向量,进入到新一代的种群中,使得无线传感器网络的网关位置坐标种群能够得到持续的优化;步骤五、重复步骤二到步骤四直到迭代次数到达P=500,在第500代种群中分别按照式(7)计算各个目标向量<img file="FDA0000724768800000031.GIF" wi="376" he="84" />对应的适应值,适应值最小的一个目标向量即为无线传感器网络中网关的最优部署位置坐标。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号
您可能感兴趣的专利