发明名称 无线传感网中能量采集节点的联合路径选择和功率分配方法
摘要 本发明公开了一种无线传感网中能量采集节点的联合路径选择和功率分配方法,属于协作通信技术领域。包括步骤:系统场景分析,问题描述;系统数学模型建立;然后利用优化方法求出最优解。本发明针对特殊的应用场景,来源实际应用,区别与以往的单独的传感器节点或者网关的资源分配,综合考虑联合传感器节点和网关的联合功率分配和路径选择,利用网关充当中继站的作用,最大化通信节点之间的吞吐量性能。本发明针对最优化问题的求解,采用凸优化处理,转化优化问题的目标函数,不经过近似计算,不影响问题的精度的同时极大的降低的计算复杂度,减少系统开销产生的时延,寻优过程采用拉格朗日乘子方法,寻优速度快,算法迭代过程中采用次梯度方法,并选用渐进步长,寻优更加精确。
申请公布号 CN106131918A 申请公布日期 2016.11.16
申请号 CN201610662163.5 申请日期 2016.08.12
申请人 梁广俊 发明人 梁广俊
分类号 H04W40/22(2009.01)I;H04W52/24(2009.01)I;H04W52/26(2009.01)I;H04W52/46(2009.01)I 主分类号 H04W40/22(2009.01)I
代理机构 代理人
主权项 一种无线传感网中能量采集节点的联合路径选择和功率分配方法,其特征在于,包括:步骤1:系统场景分析,问题描述;场景中有一个稳定供电的传感器节点A,一个能量采集的网关节点R和一个目标通信传感器节点B,稳定供电的传感器节点A和目标通信传感器节点B之间存在直达路径,能量采集的网关节点R采用解码转发工作方式,一个传输过程包括N个数据块,稳定供电的传感器节点A和能量采集的网关节点R占用相等的带宽W,联合考虑该场景下的稳定供电的传感器节点A和能量采集的网关节点R的路径选择和功率分配问题;假设能量采集模型采用伯努利过程,用于能量采集的电池容量足够大,除了用于传输消耗的能量忽略不计,一个数据块的传输时间是T<sub>P</sub>,<img file="FDA0001077467420000011.GIF" wi="62" he="71" />表示第i数据包传输时传感器节点A与中继节点R之间的信道系数,<img file="FDA0001077467420000012.GIF" wi="62" he="70" />表示第i数据包传输时传感器节点A与目标通信传感器节点B之间的信道系数,<img file="FDA0001077467420000013.GIF" wi="61" he="69" />表示第i数据包传输时中继节点R与目标通信传感器节点B之间的信道系数,采用半双工中继协作通信方式,每个数据包的传输时间内包含两个时隙,在第i个数据包传输时间的第一个时隙,如果传感器节点A传输符号x采用固定功率p<sub>a</sub>,同时x满足x~CN(0,1),那么中继R和目的通信传感器节点B的接收信号可以分别表示为:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>y</mi><mrow><mi>a</mi><mo>,</mo><mi>r</mi></mrow><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>A</mi><mn>0</mn><mfrac><mn>1</mn><mn>2</mn></mfrac></msubsup><msubsup><mi>d</mi><mn>0</mn><mfrac><mi>&alpha;</mi><mn>2</mn></mfrac></msubsup><msubsup><mi>d</mi><mo>-</mo><mrow><mo>-</mo><mfrac><mi>&alpha;</mi><mn>2</mn></mfrac></mrow></msubsup><msubsup><mi>h</mi><mrow><mi>a</mi><mo>,</mo><mi>r</mi></mrow><mi>i</mi></msubsup><msqrt><msub><mi>p</mi><mi>a</mi></msub></msqrt><mi>x</mi><mo>+</mo><msubsup><mi>n</mi><mi>r</mi><mi>i</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mrow><mi>a</mi><mo>,</mo><mi>b</mi></mrow><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>A</mi><mn>0</mn><mfrac><mn>1</mn><mn>2</mn></mfrac></msubsup><msubsup><mi>d</mi><mn>0</mn><mfrac><mi>&alpha;</mi><mn>2</mn></mfrac></msubsup><msup><mi>d</mi><mrow><mo>-</mo><mfrac><mi>&alpha;</mi><mn>2</mn></mfrac></mrow></msup><msubsup><mi>h</mi><mrow><mi>a</mi><mo>,</mo><mi>b</mi></mrow><mi>i</mi></msubsup><msqrt><msub><mi>p</mi><mi>a</mi></msub></msqrt><mi>x</mi><mo>+</mo><msubsup><mi>n</mi><mi>b</mi><mi>i</mi></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>.</mo></mrow>]]></math><img file="FDA0001077467420000014.GIF" wi="1550" he="95" /></maths>其中:d<sub>0</sub>和A<sub>0</sub>分别表示大尺度衰落的参考距离和参考功率,d<sub>1</sub>和d分别表示传感器节点A和中继节点R之间的距离以及传感器节点A和目标通信传感器节点B之间的距离,<img file="FDA0001077467420000015.GIF" wi="45" he="61" />和<img file="FDA0001077467420000016.GIF" wi="38" he="63" />分别表示第i个数据包传输时间的中继节点R和目标通信传感器节点B的接收噪声,α表示大尺度衰落因子,<img file="FDA0001077467420000017.GIF" wi="62" he="70" />和<img file="FDA0001077467420000018.GIF" wi="64" he="69" />分别表示第i个数据包传输时间的中继节点R和目标通信传感器节点B的信道增益;第一时隙中继接收的信噪比可以表示为<img file="FDA0001077467420000019.GIF" wi="643" he="94" />其中:N<sub>0</sub>表示归一化的噪声功率,W表示分配的带宽,在第二个时隙,目标通信传感器节点B采用最大比接受方式,目的通信节点可以分别接收到传感器节点A的直传信号,以及中继节点R的转发信号,第二时隙目的目标通信传感器节点B收到的来自中继节点R的信号表示为<img file="FDA00010774674200000110.GIF" wi="786" he="95" />其中:d<sub>2</sub>和<img file="FDA00010774674200000111.GIF" wi="59" he="69" />分别表示中继节点R和目标通信传感器节点B之间的距离和信道增益,<img file="FDA00010774674200000112.GIF" wi="38" he="63" />表示第i个数据包传输时间的第二时隙时目标通信传感器节点B的接收噪声,p<sub>r</sub>表示相应的中继站R的发射功率;在接收端,采用最大比接受准则的目标通信传感器节点B收到的来自传感器节点A和中继节点R的的信噪比表示<img file="FDA0001077467420000021.GIF" wi="962" he="118" />若选用直达链路,目标通信传感器节点B收到的传感器节点A的信号表示为<img file="FDA0001077467420000022.GIF" wi="789" he="102" />相应的,目标通信传感器节点B收到的传感器节点A的信噪比表示为<img file="FDA0001077467420000023.GIF" wi="684" he="93" />步骤2:平均中断概率推导;所述平均中断概率为<img file="FDA0001077467420000024.GIF" wi="266" he="126" />其中Ο<sup>i</sup>表示第i次传输的中断概率;步骤3:最优化问题归结;在上述假设前提和约束条件下,归结出最优化问题如下:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>P</mi><mn>1</mn><mo>:</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><msup><mi>r</mi><mi>i</mi></msup><mo>,</mo><msubsup><mi>p</mi><mi>r</mi><mi>i</mi></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>i</mi></mrow></munder><mover><mi>O</mi><mo>&OverBar;</mo></mover></mrow></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>i</mi></munderover><msubsup><mi>p</mi><mi>r</mi><mi>k</mi></msubsup><msub><mi>T</mi><mi>p</mi></msub><mo>&le;</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>i</mi></munderover><msubsup><mi>E</mi><mi>H</mi><mi>n</mi></msubsup><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>p</mi><mi>r</mi><mi>i</mi></msubsup><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>r</mi><mi>i</mi></msup><mo>&Element;</mo><mo>{</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>r</mi><mi>i</mi></msup><mo>&le;</mo><msup><mi>D</mi><mi>i</mi></msup><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>}</mo></mrow></mtd></mtr></mtable><mo>,</mo></mrow>]]></math><img file="FDA0001077467420000025.GIF" wi="573" he="406" /></maths>其中<img file="FDA0001077467420000026.GIF" wi="321" he="119" />表示能量因果约束,优化问题的目标函数是最小化中断概率<img file="FDA00010774674200000213.GIF" wi="63" he="47" />优化变量是每次传输的模式选择因子r<sup>i</sup>以及相应的功率分配方案<img file="FDA0001077467420000027.GIF" wi="75" he="61" />约束条件是<img file="FDA0001077467420000028.GIF" wi="1118" he="119" />步骤4:最优化问题求解;最优化问题P1的优化变量是每次传输的模式选择因子r<sup>i</sup>以及相应的功率分配方案<img file="FDA0001077467420000029.GIF" wi="67" he="64" />r<sup>i</sup>是离散变量取值范围为{1,0},<img file="FDA00010774674200000210.GIF" wi="46" he="70" />连续变量,取值范围大于等于0,采用代价函数的方法,定义S<sup>i</sup>=&lt;D<sup>i</sup>,ε<sup>i</sup>&gt;表示状态集,<img file="FDA00010774674200000211.GIF" wi="222" he="86" />表示行为集,同时定义代价函数<img file="FDA00010774674200000212.GIF" wi="830" he="204" />利用凸优化理论关于代价函数的迭代方法,可以得知:<maths num="0003"><math><![CDATA[<mrow><msup><mi>J</mi><mi>i</mi></msup><mrow><mo>(</mo><msup><mi>S</mi><mi>i</mi></msup><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>min</mi><mrow><msup><mi>a</mi><mi>N</mi></msup><mo>&Element;</mo><mi>A</mi><mrow><mo>(</mo><msup><mi>S</mi><mi>N</mi></msup><mo>)</mo></mrow></mrow></msub><msup><mi>O</mi><mi>N</mi></msup><mo>,</mo><mi>i</mi><mo>=</mo><mi>N</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>min</mi><mrow><msup><mi>a</mi><mi>N</mi></msup><mo>&Element;</mo><mi>A</mi><mrow><mo>(</mo><msup><mi>S</mi><mi>i</mi></msup><mo>)</mo></mrow></mrow></msub><msup><mi>O</mi><mi>i</mi></msup><mo>+</mo><msup><mi>J</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>&lt;</mo><mi>N</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0001077467420000031.GIF" wi="806" he="199" /></maths><maths num="0004"><math><![CDATA[<mrow><mi>A</mi><mrow><mo>(</mo><msup><mi>S</mi><mi>i</mi></msup><mo>=</mo><mo>&lt;</mo><msup><mi>D</mi><mi>i</mi></msup><mo>,</mo><msup><mi>&epsiv;</mi><mi>i</mi></msup><mo>&gt;</mo><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mo>{</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>}</mo><mo>&times;</mo><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><msup><mi>&epsiv;</mi><mi>i</mi></msup><msubsup><mi>T</mi><mi>p</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>&rsqb;</mo><mo>,</mo><msup><mi>D</mi><mi>i</mi></msup><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><mo>{</mo><mn>0</mn><mo>}</mo><mo>&times;</mo><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><msup><mi>&epsiv;</mi><mi>i</mi></msup><msubsup><mi>T</mi><mi>p</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>&rsqb;</mo><mo>,</mo><msup><mi>D</mi><mi>i</mi></msup><mo>=</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0001077467420000032.GIF" wi="807" he="167" /></maths>最终,最优解可以表示成<img file="FDA0001077467420000033.GIF" wi="170" he="137" />
地址 241000 安徽省芜湖市镜湖区绿地镜湖世纪城伊顿公馆3-2-1602