发明名称 无线传感网中速率和可靠性协同的跨层优化方法
摘要 本发明公开了一种无线传感网中速率和可靠性协同的跨层优化方法,本发明通过利用变量替换以及引入中间变量的方法将一个非凸的无线传感网速率和可靠性协同的跨层优化问题转化成一个凸问题,然后再利用对偶分解和次梯度方法,设计了分布式优化算法,可分布式求解转化后的凸问题。该方法兼顾了无线传感网中速率和可靠性这两个重要性能指标,提出的分布式优化算法便于转化成无线传感网实际实施的协议。
申请公布号 CN105792258A 申请公布日期 2016.07.20
申请号 CN201610227725.3 申请日期 2016.04.13
申请人 浙江理工大学 发明人 徐伟强;魏良晓;史清江;王成群;吕文涛
分类号 H04W24/02(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W24/02(2009.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 邱启旺
主权项 一种无线传感网络中速率和可靠性协同的跨层优化方法,其特征在于,包括以下步骤:(1)建立无线传感网速率和可靠性协同的跨层优化问题P<sub>1</sub>:<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><munder><mi>max</mi><mrow><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>,</mo><msub><mi>x</mi><mi>f</mi></msub><mo>,</mo><msub><mi>p</mi><mi>l</mi></msub></mrow></munder></mtd><mtd><mrow><msub><mi>&omega;</mi><mn>1</mn></msub><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><msub><mi>U</mi><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&omega;</mi><mn>1</mn></msub><mo>)</mo></mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&omega;</mi><mi>f</mi></msub><mi>E</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000011.GIF" wi="1142" he="127" /></maths><maths num="0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><mfrac><msub><mi>x</mi><mi>f</mi></msub><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub></mfrac><mo>&le;</mo><mi>B</mi><mi>log</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msub><mi>&sigma;</mi><mi>l</mi></msub></mfrac><msub><mi>p</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000012.GIF" wi="957" he="165" /></maths><maths num="0003"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>f</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>x</mi><mi>f</mi></msub><mo>&le;</mo><msubsup><mi>x</mi><mi>f</mi><mi>max</mi></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow>]]></math><img file="FDA0000964538510000013.GIF" wi="677" he="103" /></maths><maths num="0004"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><msub><mi>p</mi><mi>l</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>n</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>n</mi><mo>&Element;</mo><mi>N</mi></mrow>]]></math><img file="FDA0000964538510000014.GIF" wi="590" he="150" /></maths><maths num="0005"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>,</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000964538510000015.GIF" wi="606" he="84" /></maths>其中,x<sub>f</sub>表示会话流的传输速率,U<sub>f</sub>(x<sub>f</sub>)表示速率的效用函数,r<sub>l,f</sub>表示码率,B表示链路l的带宽,σ<sub>l</sub>表示噪声干扰,p<sub>l</sub>表示在链路l上消耗的功率,p<sub>n</sub>表示进行传输数据时节点n消耗的功率,E(r<sub>l,f</sub>)表示会话f使用链路l的出错概率,定义为一个关于码率r<sub>l,f</sub>的函数,为了权衡会话流f的传输速率和错误概率,设定权重值w<sub>1</sub>,w<sub>1</sub>可取[0,1]之间的任何一个值;为了保持前后数量级一致,添加w<sub>f</sub>。(2)将非凸的问题P<sub>1</sub>转化为凸问题P<sub>2</sub>:<maths num="0006"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><munder><mi>max</mi><mrow><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup><mo>,</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>,</mo><msub><mi>p</mi><mi>l</mi></msub></mrow></munder></mtd><mtd><mrow><msub><mi>&omega;</mi><mn>1</mn></msub><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><msub><mi>U</mi><mi>f</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup></msup><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&omega;</mi><mi>f</mi></msub><mi>E</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000016.GIF" wi="1198" he="135" /></maths><maths num="0007"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><msubsup><mi>x</mi><mi>f</mi><mrow><mo>&prime;</mo><mi>min</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>f</mi><mrow><mo>&prime;</mo><mi>max</mi></mrow></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000017.GIF" wi="822" he="103" /></maths><maths num="0008"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>logr</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><msub><mi>logc</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi><mo>,</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000964538510000018.GIF" wi="1029" he="95" /></maths><maths num="0009"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>,</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000964538510000019.GIF" wi="605" he="85" /></maths><maths num="0010"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mi>log</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msub><mi>&sigma;</mi><mi>l</mi></msub></mfrac><msub><mi>p</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow>]]></math><img file="FDA00009645385100000110.GIF" wi="844" he="167" /></maths><maths num="0011"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><msub><mi>p</mi><mi>l</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>n</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>n</mi><mo>&Element;</mo><mi>N</mi></mrow>]]></math><img file="FDA0000964538510000021.GIF" wi="589" he="143" /></maths>其中x′<sub>f</sub>=log(x<sub>f</sub>),<img file="FDA0000964538510000022.GIF" wi="866" he="87" />c<sub>l,f</sub>表示l链路上f会话的中间变量;(3)采用对偶分解方法和次梯度方法对问题P<sub>2</sub>进行分布式求解。问题P<sub>2</sub>的对偶问题P<sub>3</sub>为:max D(β)s.t. β≥0其中,β是对偶因子。使用次梯度算法求解该问题。对偶问题P<sub>3</sub>的目标函数如下所示:<maths num="0012"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>x</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>c</mi><mo>,</mo><mi>p</mi><mo>&Element;</mo><mi>&Gamma;</mi></mrow></munder><mi>L</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>c</mi><mo>,</mo><mi>p</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000964538510000023.GIF" wi="732" he="103" /></maths><maths num="0013"><math><![CDATA[<mrow><mi>&Gamma;</mi><mo>=</mo><mfenced open = "{" close = "}"><mtable><mtr><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msub><mi>&sigma;</mi><mi>l</mi></msub></mfrac><msub><mi>p</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow></mtd></mtr><mtr><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msub></mrow></munder><msub><mi>p</mi><mi>l</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>n</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>n</mi><mo>&Element;</mo><mi>N</mi></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mi>f</mi><mrow><mo>&prime;</mo><mi>min</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>f</mi><mrow><mo>&prime;</mo><mi>max</mi></mrow></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>&le;</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi><mo>,</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000964538510000024.GIF" wi="939" he="542" /></maths>因此,问题P<sub>2</sub>的拉格朗日函数如下:<maths num="0014"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>L</mi><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>c</mi><mo>,</mo><mi>p</mi><mo>,</mo><mi>&beta;</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>&omega;</mi><mn>1</mn></msub><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><msub><mi>U</mi><mi>f</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup></msup><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&omega;</mi><mi>f</mi></msub><mi>E</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow></munder><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>logr</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>+</mo><msub><mi>logc</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>-</mo><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000025.GIF" wi="1454" he="294" /></maths>(4)将D(β)分解为可分布式求解的三类子问题:子问题一:<maths num="0015"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><munder><mi>max</mi><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup></munder></mtd><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><mo>{</mo><msub><mi>&omega;</mi><mn>1</mn></msub><msub><mi>U</mi><mi>f</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup></msup><mo>)</mo></mrow><mo>-</mo><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup><mo>}</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000026.GIF" wi="862" he="183" /></maths><maths num="0016"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><msubsup><mi>x</mi><mi>f</mi><mrow><mo>&prime;</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>f</mi><mrow><mo>&prime;</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000027.GIF" wi="829" he="110" /></maths>子问题二:<maths num="0017"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><munder><mi>max</mi><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub></munder><mi></mi><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&omega;</mi><mn>1</mn></msub><mo>)</mo></mrow><munder><mi>&Sigma;</mi><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&omega;</mi><mi>f</mi></msub><mi>E</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mi>log</mi><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&DoubleLeftRightArrow;</mo><munder><mi>max</mi><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub></munder><munder><mi>&Sigma;</mi><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><mo>{</mo><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mi>log</mi><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&omega;</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>f</mi></msub><mi>E</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>)</mo></mrow><mo>}</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000031.GIF" wi="1341" he="286" /></maths><maths num="0018"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mn>0</mn><mo>&le;</mo><msub><mi>r</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi><mo>,</mo><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000032.GIF" wi="846" he="86" /></maths>子问题三:<maths num="0019"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><munder><mi>max</mi><mrow><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>,</mo><msub><mi>p</mi><mi>l</mi></msub></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><mrow><mo>(</mo><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mi>log</mi><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>&DoubleLeftRightArrow;</mo><munder><mi>max</mi><mrow><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>,</mo><msub><mi>p</mi><mi>l</mi></msub></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>n</mi><mo>&Element;</mo><mi>N</mi></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>L</mi><mi>out</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mi>log</mi><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000033.GIF" wi="942" he="311" /></maths><maths num="0020"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mi>B</mi><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msub><mi>&sigma;</mi><mi>l</mi></msub></mfrac><msub><mi>p</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000964538510000034.GIF" wi="998" he="167" /></maths><maths num="0021"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><msub><mi>p</mi><mi>l</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>n</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>n</mi><mo>&Element;</mo><mi>N</mi></mrow>]]></math><img file="FDA0000964538510000035.GIF" wi="590" he="143" /></maths>其中,β<sub>l,f</sub>是对偶因子。(5)使用分布式方法将P<sub>3</sub>进行分布式求解,具体包括以下子步骤:(5.1)初始化β<sub>l,f</sub>(0),迭代次数t=1;将β<sub>l,f</sub>(0)带入三个对偶子问题中,得到c<sub>l,f</sub>(0)、p<sub>l</sub>(0)、r<sub>l,f</sub>(0)和x′<sub>f</sub>(0);(5.2)使用次梯度算法求解P<sub>3</sub>,即通过下式求出第t次迭代的β<sub>l,f</sub>即β<sub>l,f</sub>(t):β<sub>l,f</sub>(t+1)=[β<sub>l,f</sub>(t)‑k<sub>β</sub>(t)(log c<sub>l,f</sub>(t)+log r<sub>l,f</sub>(t)‑x′<sub>f</sub>(t))]<sup>+</sup>其中k<sub>β</sub>(t)表示的是步长,[z]<sup>+</sup>=max{0,z};(5.3)在会话流上更新会话f的速率x′<sub>f</sub>,最大化:<maths num="0022"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mn>1</mn></msub><msub><mi>U</mi><mi>f</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup></msup><mo>)</mo></mrow><mo>-</mo><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><mi>L</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><msubsup><mi>x</mi><mi>f</mi><mo>&prime;</mo></msubsup></mrow>]]></math><img file="FDA0000964538510000036.GIF" wi="501" he="127" /></maths>且x′<sub>f</sub>在<img file="FDA0000964538510000037.GIF" wi="266" he="83" />范围之内,其中<img file="FDA0000964538510000038.GIF" wi="155" he="79" />为速率的效用函数,L(f)表示的是会话f流过链路的链路集合,设定一个权重值w<sub>1</sub>;(5.4)在链路和链路上的会话更新码率r<sub>l,f</sub>,最大化:β<sub>l,f</sub> log r<sub>l,f</sub>‑(1‑ω<sub>1</sub>)ω<sub>f</sub>E(r<sub>l,f</sub>)且r<sub>l,f</sub>在[0,1]的范围之内,其中E(r<sub>l,f</sub>)是错误概率的函数,为了保持前后数量级一致,从而添加一个w<sub>f</sub>;(5<sub>.</sub>5)在节点上更新中间变量c<sub>l,f</sub>和功率p<sub>l</sub>,最大化:<maths num="0023"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>&beta;</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><msub><mi>logc</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub></mrow>]]></math><img file="FDA0000964538510000041.GIF" wi="629" he="143" /></maths>且c<sub>l,f</sub>在两个约束之下,分别为:<maths num="0024"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>f</mi><mo>&Element;</mo><mi>F</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>c</mi><mrow><mi>l</mi><mo>,</mo><mi>f</mi></mrow></msub><mo>&le;</mo><mi>B</mi><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msub><mi>&sigma;</mi><mi>l</mi></msub></mfrac><msub><mi>p</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&Element;</mo><mi>L</mi></mrow>]]></math><img file="FDA0000964538510000042.GIF" wi="865" he="166" /></maths><maths num="0025"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>l</mi><mo>&Element;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><msub><mi>p</mi><mi>l</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>n</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>n</mi><mo>&Element;</mo><mi>N</mi></mrow>]]></math><img file="FDA0000964538510000043.GIF" wi="582" he="143" /></maths>(5.6)重复步骤(5.1)到步骤(5.5),直到目标函数收敛,得到无线传感网络的最优码率r<sup>*</sup>、会话速率x<sup>*</sup>以及链路l上消耗的功率p<sup>*</sup>,从而实现无线传感网络速率和可靠性的协同跨层优化。
地址 310018 浙江省杭州市江干经济开发区白杨街道2号大街928号