发明名称 有干扰环境下移动通信网络中能量效率最优化的方法
摘要 本发明提出一种有干扰环境下移动通信网络中能量效率(energy efficiency,EE)最优化的方法。首先将EE最优化问题建模为一个非凸优化问题。通过使用分式规划,可将此非凸优化问题从分式形式转化为减法形式。再使用凸函数之差规划即可将减法形式的EE最优化问题从非凸优化问题转化为凸优化问题。最后使用内点法求解此问题,即可得到EE最优值。此方法简单且结果准确。
申请公布号 CN103945402A 申请公布日期 2014.07.23
申请号 CN201410206508.7 申请日期 2014.05.15
申请人 东南大学 发明人 尤肖虎;李知航;潘志文;童恩
分类号 H04W16/18(2009.01)I;H04W24/02(2009.01)I;H04W72/04(2009.01)I 主分类号 H04W16/18(2009.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种有干扰环境下移动通信网络中能量效率最优化的方法,包括以下步骤:步骤1:建立通信网络,其中所述通信网络中包括M个基站,每个基站覆盖范围内均只有一个用户设备,设定迭代次数t,初次迭代时,令t=1;容忍度δ=0.01和通信网络能量效率的初始最优值为q<sub>1</sub>;步骤2:定义通信网络中所有用户设备获得的速率和为R<sub>T</sub>(P),根据公式<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>R</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>lo</mi><msub><mi>g</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>p</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>ii</mi></msub></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow><mi>M</mi></msubsup><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>ji</mi></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>z</mi><mn>2</mn></msubsup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000505520110000011.GIF" wi="861" he="167" /></maths>计算得到所有用户设备获得的速率和为R<sub>T</sub>(P),其中P=[p<sub>1</sub>,p<sub>2</sub>,...,p<sub>M</sub>]为每个基站分配给自身用户设备的发射功率向量,对于<img file="FDA0000505520110000016.GIF" wi="179" he="57" />p<sub>i</sub>表示基站i分配给自身用户设备的发射功率,g<sub>ji</sub>为基站j与用户设备i之间的信道冲击响应,<img file="FDA0000505520110000017.GIF" wi="62" he="64" />为复可加白高斯噪声功率,通过上行链路或周期报告,用户设备i可将g<sub>ji</sub>反馈给自身服务基站i;步骤3:根据公式<img file="FDA0000505520110000012.GIF" wi="422" he="94" />获得通信网络消耗总功率P<sub>T</sub>(P),其中P<sub>C</sub>为通信网络中电路消耗的功率;根据公式<img file="FDA0000505520110000013.GIF" wi="358" he="140" />获得通信网络能量效率η<sub>EE</sub>(P),定义通信网络能量效率的最优化函数为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>max</mi><mo>{</mo><msub><mi>&eta;</mi><mi>EE</mi></msub><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>}</mo></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mn>0</mn><mo>&le;</mo><msub><mi>p</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>max</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>&Element;</mo><mi>M</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000505520110000014.GIF" wi="550" he="147" /></maths>其中p<sub>max</sub>为任意基站分配给自身用户设备的最大发射功率,p<sub>max</sub>是移动通信网络中给定的系统参数;步骤4:通过使用分式规划法,可将通信网络能量效率的最优化函数转化为如下问题:(P1)  <maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>max</mi><mo>{</mo><msub><mi>R</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>q</mi><mi>t</mi></msub><mo>&CenterDot;</mo><msub><mi>P</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>}</mo></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mn>0</mn><mo>&le;</mo><msub><mi>P</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>max</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>&Element;</mo><mi>M</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000505520110000015.GIF" wi="551" he="147" /></maths>其中q<sub>t</sub>为网络能量效率的最大值,在初次计算时,令q<sub>t</sub>=q<sub>1</sub>;步骤5:将问题(P1)的目标函数重写为:R<sub>T</sub>(P)‑q<sub>t</sub>·P<sub>T</sub>(P)=f(P)‑h(P)其中<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mtext>=</mtext><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>lo</mi><msub><mi>g</mi><mn>2</mn></msub><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>ji</mi></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>q</mi><mi>t</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>p</mi><mi>i</mi></msub><mo>+</mo><msub><mi>P</mi><mi>C</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000505520110000021.GIF" wi="1021" he="147" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>h</mi><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>lo</mi><msub><mi>g</mi><mn>2</mn></msub><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow><mi>M</mi></munderover><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>ji</mi></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000505520110000022.GIF" wi="704" he="149" /></maths>因此,问题(P1)等价于如下问题:(P2)  <maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>max</mi><mo>{</mo><mi>f</mi><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>-</mo><mi>h</mi><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>}</mo></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mn>0</mn><mo>&le;</mo><msub><mi>p</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>max</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>&Element;</mo><mi>M</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000505520110000023.GIF" wi="551" he="150" /></maths>问题(P2)是一个标准的凸函数之差规划;步骤6:令k=1,随机设置P<sup>k</sup>,只需使P<sup>k</sup>中的每个元素均位于区间[0,p<sub>max</sub>]内即可,计算I<sup>k</sup>=f(P<sup>k</sup>)‑h(P<sup>k</sup>);步骤7:计算h(P)的梯度,其可记为<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mo>&dtri;</mo><mi>h</mi><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><msub><mi>e</mi><mi>i</mi></msub><mrow><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow><mi>M</mi></munderover><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>ji</mi></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>ln</mi><mn>2</mn></mrow></mfrac></mrow>]]></math><img file="FDA0000505520110000024.GIF" wi="757" he="219" /></maths>其中e<sub>i</sub>是一个维数为M的列向量,具体而言:e<sub>i</sub>(i)=0,当j≠i时,e<sub>i</sub>(j)=g<sub>ji</sub>;步骤8:将h(P)用它的一阶泰勒公式展开为:h(P)=h(P<sup>k</sup>)+▽h<sup>T</sup>(P<sup>k</sup>)·(P‑P<sup>k</sup>)此时可将问题(P2)转换为下述问题:(P3)  <maths num="0008" id="cmaths0008"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>max</mi><mo>{</mo><mi>f</mi><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>-</mo><mo>[</mo><mi>h</mi><mrow><mo>(</mo><msup><mi>P</mi><mi>k</mi></msup><mo>)</mo></mrow><mo>+</mo><mo>&dtri;</mo><msup><mi>h</mi><mi>T</mi></msup><mrow><mo>(</mo><msup><mi>P</mi><mi>k</mi></msup><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>P</mi><mo>-</mo><msup><mi>P</mi><mi>k</mi></msup><mo>)</mo></mrow><mo>]</mo><mo>}</mo></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mn>0</mn><mo>&le;</mo><msub><mi>p</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>p</mi><mi>max</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>&Element;</mo><mi>M</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000505520110000025.GIF" wi="926" he="153" /></maths>步骤9:问题(P3)是一个标准的凸优化问题,通过使用内点法可以进行求解,具体为:令F(P)=f(P)‑[h(P<sup>k</sup>)+▽h<sup>T</sup>(P<sup>k</sup>)·(P‑P<sup>k</sup>)],初始化P<sup>*</sup>=P<sup>k</sup>,s为一个大于0的正数,u为一个小于1的正数;步骤10:以P<sup>*</sup>为起始点,计算<img file="FDA0000505520110000026.GIF" wi="610" he="96" />的最大值,记为<img file="FDA0000505520110000027.GIF" wi="93" he="86" />步骤11:更新P<sup>*</sup>,令<img file="FDA0000505520110000031.GIF" wi="182" he="65" />如果s/M<δ,则执行步骤12。否则令s=s·u,并返回执行步骤10;步骤12:令k=k+1,且P<sup>k</sup>=P<sup>*</sup>,计算I<sup>k</sup>=f(P<sup>k</sup>)‑h(P<sup>k</sup>);步骤13:如果|I<sup>k</sup>‑I<sup>k‑1</sup>|≤δ,则执行步骤14;否则返回执行步骤8;步骤14:令P<sup>t</sup>=P<sup>k</sup>,如果R<sub>T</sub>(P<sub>t</sub>)‑q<sub>t</sub>·P<sub>T</sub>(P<sub>t</sub>)>δ,则执行步骤15;否则执行步骤16;步骤15:令q<sub>t+1</sub>=R<sub>T</sub>(P<sub>t</sub>)/P<sub>T</sub>(P<sub>t</sub>),t=t+1,其中t为当前迭代次数,返回执行步骤4;步骤16:结束,获得通信网络能量效率的最大值为q<sub>t</sub>,相应的用户设备获得的最优发射功率向量为P<sub>t</sub>。
地址 210096 江苏省南京市四牌楼2号