发明名称 人体行走姿态下无线体域网生存周期的汇聚节点放置方法
摘要 本发明涉及一种人体行走姿态下最大化无线体域网生存周期的汇聚节点放置方法,属于无线通信领域,包括构建无线体域网网络生存期最大化模型;利用无线体域网信道模型、功率控制模型以及能量消耗模型,推导最大化无线体域网生存周期的汇聚节点放置方法;提出汇聚节点放置方法的近似求解算法。本发明针对无线体域网中汇聚节点的放置位置问题,克服了传统的将汇聚节点固定在人体腰部位置引起能量消耗增加的缺点,同时也首次将该方法应用到拓扑动态变化的无线体域网中。在满足比特误码率以及功率约束条件下,该方法最大化了网络的生存期,更加适用于人体生理参数的全天候监测。
申请公布号 CN105722102A 申请公布日期 2016.06.29
申请号 CN201610070941.1 申请日期 2016.02.02
申请人 吉林大学 发明人 胡封晔;王丰;刘小兰;熊威;王珊珊;邵美琪;隋丹;王丽恒
分类号 H04W16/18(2009.01)I;H04W16/22(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W16/18(2009.01)I
代理机构 吉林长春新纪元专利代理有限责任公司 22100 代理人 魏征骥
主权项 一种人体行走姿态下无线体域网生存周期的汇聚节点放置方法,其特征在于包括下列步骤:(1)无线体域网网络生存期最大化模型,(1.1)初始化初始化的过程包括:采集人体生理参数的传感器节点按照其自身功能放置在人体的特定位置;(1.2)传感器节点选择需要采集人体的某项生理参数信息时,一个传感器节点醒来进行信息传递,其它传感器节点仍旧保持睡眠状态,节点的选择遵照基于姿态感知的生存期最大化协议(PA‑DPLM,Posture‑Aware Dynamic Protocol for Lifetime Maximization),代表该协议的高能效函数f(e<sub>i</sub>,c<sub>i</sub>,p<sub>i</sub>)如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>EI</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow><mrow><msub><mi>e</mi><mi>c</mi></msub><mo>+</mo><mn>1</mn><mo>/</mo><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow></mfrac><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>N</mi><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>L</mi></mrow>]]></math><img file="FDA0000920279900000011.GIF" wi="1205" he="155" /></maths>其中EI<sub>i,j</sub>为传感器节点i在j次传输时的高能效系数,e<sub>i,j</sub>为节点i当前的能量,c<sub>i,j</sub>为节点i的信道状态信息,p<sub>i,j</sub>为节点i的姿态信息,e<sub>c</sub>为传输链路所消耗的能量,N为无线体域网中传感器节点的数量,L为无线体域网的网络生存期;(1.3)网络死亡判断当采集某类人体生理参数的全部传感器节点都不能进行数据采集时,则无线体域网进入死亡状态;(2)汇聚节点放置方法汇聚节点根据下文阐述的最优放置方法,放置到最优位置使得无线体域网生存期最大化,步骤如下:(2.1)路径损耗模型L<sub>path</sub>(d<sub>i</sub>)=a·log<sub>10</sub>d<sub>i</sub>+b+N<sub>0</sub>  [dB]其中L<sub>path</sub>(d<sub>i</sub>)是节点i当前到汇聚节点的路径损耗,a,b,N<sub>0</sub>是给定的参数a=‑8.6,b=‑20.3,N<sub>0</sub>=10<sup>‑7</sup>,行走姿态下d<sub>i</sub>的表达式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>a</mi><mi>p</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mrow><mi>a</mi><mi>p</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>z</mi><mrow><mi>a</mi><mi>p</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FDA0000920279900000012.GIF" wi="927" he="107" /></maths>其中(x<sub>ap</sub>,y<sub>ap</sub>,z<sub>ap</sub>)和(x<sub>i</sub>,y<sub>i</sub>,z<sub>i</sub>)分别是汇聚节点和节点i的坐标;(2.2)功率控制模型比特误码率(BER,Bit Error Rate)约束条件下的功率控制模型如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mrow><mo>{</mo><mrow><mfrac><mrow><msub><mi>AM</mi><mn>0</mn></msub></mrow><mrow><msubsup><mi>q&xi;&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><msubsup><mi>log</mi><mn>2</mn><mi>M</mi></msubsup></mrow></mfrac><mo>,</mo><msub><mi>P</mi><mi>max</mi></msub></mrow><mo>}</mo></mrow><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub></mrow>]]></math><img file="FDA0000920279900000021.GIF" wi="608" he="187" /></maths>其中,P<sub>i,j</sub>为节点i在j次传输时的功率;A=0.25,q=1,M<sub>0</sub>=10<sup>‑7</sup>;ξ为比特误码率约束条件;由于调制方式是二进制相移键控,所以M=2;P<sub>max</sub>是最大传输功率。<img file="FDA0000920279900000029.GIF" wi="77" he="83" />为节点i在j次传输时的信道方差,表达式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>&delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>&eta;d</mi><mi>i</mi><mrow><mo>-</mo><mi>&alpha;</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000920279900000022.GIF" wi="247" he="91" /></maths>其中η=1,α是路径损耗因子,可以通过路径损耗模型计算得到;(2.3)能量消耗模型在数据包的发送时间为T<sub>d</sub>时,节点i在j次传输时的能量消耗模型如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup><mo>=</mo><msub><mi>e</mi><mi>c</mi></msub><mo>+</mo><mn>1</mn><mo>/</mo><msub><mi>c</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>T</mi><mi>d</mi></msub></mrow>]]></math><img file="FDA0000920279900000023.GIF" wi="538" he="91" /></maths>其中<img file="FDA0000920279900000024.GIF" wi="67" he="82" />是节点i在j次传输时消耗的能量,c<sub>i,j</sub>为此时的信道状态,将上述的路径损耗模型以及功率控制模型代入上式,可以得到能量消耗模型的完全表达式:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup><mo>=</mo><msub><mi>e</mi><mi>c</mi></msub><mo>+</mo><msub><mrow><mo>{</mo><mrow><mfrac><mrow><msub><mi>AM</mi><mn>0</mn></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mrow><msubsup><mi>q&xi;&eta;d</mi><mi>i</mi><mrow><mo>-</mo><mi>&alpha;</mi></mrow></msubsup><msubsup><mi>log</mi><mn>2</mn><mi>M</mi></msubsup></mrow></mfrac><mo>,</mo><msub><mi>P</mi><mi>max</mi></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mo>}</mo></mrow><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub></mrow>]]></math><img file="FDA0000920279900000025.GIF" wi="803" he="175" /></maths>(2.4)优化目标优化目标是最大化无线体域网的生存期,优化目标和约束条件表示如下:maximize Lsubject to<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msub><mi>e</mi><mn>0</mn></msub><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&le;</mo><msub><mi>e</mi><mi>c</mi></msub></mrow>]]></math><img file="FDA0000920279900000026.GIF" wi="496" he="148" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>/</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup><mo>=</mo><mi>m</mi><mi>a</mi><mi>x</mi><mo>{</mo><msub><mi>EI</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>}</mo><mo>=</mo><mi>m</mi><mi>a</mi><mi>x</mi><mo>{</mo><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>/</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup><mo>}</mo></mrow>]]></math><img file="FDA0000920279900000027.GIF" wi="934" he="97" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&gt;</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000920279900000028.GIF" wi="611" he="91" /></maths>0≤i≤N0≤j≤L其中e<sub>0</sub>为传感器节点的初始能量;X<sub>i,j</sub>为选择函数,若节点i在j次传输时被选中传输数据,其值为1,否则为0;max{EI<sub>i,j</sub>}表示EI<sub>i,j</sub>的最大值;(2.5)最优放置方法将(2.3)中的能量消耗模型的完全表达式代入(2.4)中的优化目标和约束条件,可以得到:maximize Lsubject to<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msub><mi>e</mi><mn>0</mn></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mrow><mo>(</mo><mrow><msub><mi>e</mi><mi>c</mi></msub><mo>+</mo><mrow><mo>{</mo><mrow><mfrac><mrow><msub><mi>AM</mi><mn>0</mn></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mrow><msubsup><mi>q&xi;&eta;d</mi><mi>i</mi><mrow><mo>-</mo><mi>&alpha;</mi></mrow></msubsup><msubsup><mi>log</mi><mn>2</mn><mi>M</mi></msubsup></mrow></mfrac><mo>,</mo><msub><mi>P</mi><mi>max</mi></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mo>}</mo></mrow></mrow><mo>)</mo></mrow><mi>min</mi></msub><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&le;</mo><msub><mi>e</mi><mi>c</mi></msub></mrow>]]></math><img file="FDA0000920279900000031.GIF" wi="1184" he="185" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mfrac><mrow><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup></mfrac><mo>=</mo><mi>max</mi><mrow><mo>{</mo><mrow><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>/</mo><mrow><mo>(</mo><mrow><msub><mi>e</mi><mi>c</mi></msub><mo>+</mo><msub><mrow><mo>{</mo><mrow><mfrac><mrow><msub><mi>AM</mi><mn>0</mn></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mrow><msubsup><mi>q&xi;&eta;d</mi><mi>i</mi><mrow><mo>-</mo><mi>&alpha;</mi></mrow></msubsup><msubsup><mi>log</mi><mn>2</mn><mi>M</mi></msubsup></mrow></mfrac><mo>,</mo><msub><mi>P</mi><mi>max</mi></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mo>}</mo></mrow><mi>min</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow></mrow>]]></math><img file="FDA0000920279900000032.GIF" wi="1198" he="219" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><mrow><mo>(</mo><mrow><msub><mi>e</mi><mi>c</mi></msub><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mrow><mo>{</mo><mrow><mfrac><mrow><msub><mi>AM</mi><mn>0</mn></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mrow><msubsup><mi>q&xi;&eta;d</mi><mi>i</mi><mrow><mo>-</mo><mi>&alpha;</mi></mrow></msubsup><msubsup><mi>log</mi><mn>2</mn><mi>M</mi></msubsup></mrow></mfrac><mo>,</mo><msub><mi>P</mi><mi>max</mi></msub><msub><mi>T</mi><mi>d</mi></msub></mrow><mo>}</mo></mrow><mi>min</mi></msub><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mo>)</mo></mrow><mo>&gt;</mo><msubsup><mi>e</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>t</mi><mi>x</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000920279900000033.GIF" wi="1256" he="199" /></maths>0≤i≤N0≤j≤L通过求解上述优化问题,可以得到汇聚节点的最优放置位置;(3)求解汇聚节点放置方法采用一种近似求解算法,通过枚举的方式来求解汇聚节点的坐标,具体算法如下:首先,假定传感器节点分布在1m×1m的人体躯干上,鉴于目前绝大多数生理参数传感器的面积大于1cm×1cm,将该1m×1m平面等分成100份,这样,该平面类似于围棋盘,在其它传感器节点位置固定的前提下,需要找到一个最优的位置放置汇聚节点这一棋子,也就是说,给定了汇聚节点的坐标的取值集合;然后,将集合内的坐标值逐一代入(2.5)中的优化问题,进而得到与坐标集合一一对应的无线体域网生存期的取值集合,这种对应关系类似于三维坐标系中x‑y面与z轴的关系,那么利用这种对应关系绘制一个三维曲面,使用MATLAB软件中的surf(x,y,z)函数,该函数的输入参数x、y、z分别代表三维坐标系中的三个坐标轴上的取值,给定这些取值后surf(x,y,z)函数就能绘制相应的三维曲面,所以,利用已有生存期的取值集合以及汇聚节点的坐标的取值集合,就可以绘制以汇聚节点的坐标为自变量的生存期三维曲面;最后,在这一曲面上找到生存期的最大值以及该值对应的坐标,该坐标即为汇聚节点的最优放置位置。
地址 130012 吉林省长春市前进大街2699号