发明名称 基于分式规划和罚函数方法的高能效资源优化方法
摘要 本发明公开了一种基于分式规划和罚函数方法的高能效资源优化方法,包括建立所有D2D用户能效的数学模型;通过分式规划的思想和罚函数的方法将优化D2D整体能效函数(分式函数)转化为两层优化问题;第一层优化问题由每组D2D对各自单独完成功率控制,为求解函数最小值问题;第二层优化问题中提出了一种启发式资源分配算法,来解决第一层优化问题遗留的资源冲突问题。本发明能够提高移动通信D2D系统能效,满足绿色通信和延长移动终端电池使用时间的要求。
申请公布号 CN103428767A 申请公布日期 2013.12.04
申请号 CN201310411012.9 申请日期 2013.09.11
申请人 东南大学 发明人 蒋雁翔;刘强;尤肖虎
分类号 H04W28/06(2009.01)I;H04W72/04(2009.01)I 主分类号 H04W28/06(2009.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 1.基于分式规划和罚函数方法的高能效资源优化方法,其特征在于:包括如下步骤:(1)建立能效的目标函数,如式(1)、式(2)所示:<maths num="0001"><![CDATA[<math><mrow><msub><mrow><mi>max</mi><mi>U</mi></mrow><mi>EE</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>P</mi><mi>C</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>w</mi><mo>&CenterDot;</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mrow><mi>D</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow><mrow><msub><mi>p</mi><mrow><mi>I</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>+</mo><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>该目标函数包括如下约束条件:①每组D2D对的最低传输速率要求,即最低传输速率不能小于γ<sub>i</sub>,不同组的D2D对的最低传输速率可以不同:<maths num="0003"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>,</mo><mo>&ForAll;</mo><mi>j</mi></mrow></math>]]></maths>②与一组D2D对共享相同资源块的蜂窝用户对该组D2D对中的两个D2D用户的干扰功率必须小于一定值τ:<maths num="0004"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>p</mi><mrow><mi>I</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mrow><mi>k</mi><mo>.</mo><mi>i</mi></mrow></msub><mo>&le;</mo><mi>&tau;</mi></mtd><mtd><mo>&ForAll;</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mn>1</mn></mtd></mtr></mtable></mfenced></math>]]></maths>③x<sub>i,j</sub>取1表示第i组D2D对选择第j个资源块进行复用,x<sub>i,j</sub>取0表示第i组D2D对不选择第j个资源块进行复用:<maths num="0005"><![CDATA[<math><mrow><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&Element;</mo><mo>{</mo><mn>0,1</mn><mo>}</mo><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>,</mo><mo>&ForAll;</mo><mi>j</mi></mrow></math>]]></maths>④每组D2D对能且只能复用一个蜂窝用户的资源块:<maths num="0006"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>i</mi></mrow></math>]]></maths>⑤每个蜂窝用户的资源块最多只能被一组D2D对复用:<maths num="0007"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>j</mi></mrow></math>]]></maths>⑥D2D用户的最大传输功率限定:<maths num="0008"><![CDATA[<math><mrow><mn>0</mn><mo>&le;</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&le;</mo><msub><mi>p</mi><mi>max</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>,</mo><mo>&ForAll;</mo><mi>j</mi></mrow></math>]]></maths>每组D2D对中包含两个D2D用户,其中一个为接收用户,另一个为发送用户;其中:N<sub>d</sub>表示D2D对的组数,M表示可分配资源块的个数,i表示第i组D2D对,j表示第j个资源块,k表示与第i组D2D对复用相同资源块的蜂窝用户的序号;U<sub>EE</sub>表示所有D2D对的能效之和,<img file="FDA0000380334010000021.GIF" wi="345" he="132" />表示所有D2D对传输速率之和,<img file="FDA0000380334010000022.GIF" wi="369" he="132" />表示所有D2D对消耗功率之和,p<sub>i,j</sub>表示第i组D2D对在第j个资源块上的传输功率,P<sub>C</sub>表示移动终端上电路所消耗的功率,w表示资源块的带宽,h<sub>D,i</sub>表示同一组D2D对中发射用户和接收用户之间的信道增益,p<sub>I,k,i</sub>表示与一组D2D对共享资源块的蜂窝用户的发射功率,h<sub>k,i</sub>表示共享同一资源块的蜂窝用户与D2D对中的接收用户间的信道增益,σ<sup>2</sup>表示高斯白噪声的方差;(2)根据分式规划将式(1)等效转化为式(3):<maths num="0009"><![CDATA[<math><mrow><munder><mi>max</mi><mrow><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&Element;</mo><mi>&Omega;</mi></mrow></munder><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>q</mi><mo>*</mo></msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><msub><mi>P</mi><mrow><mi>C</mi><mo>_</mo><mi>ave</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中:Ω表示由约束条件①~⑥共同定义出的可行域,q<sup>*</sup>表示D2D用户整体的最优能效,P<sub>C_ave</sub>=P<sub>C</sub>/N<sub>d</sub>;而式(3)与式(4)为等效问题:<maths num="0010"><![CDATA[<math><mrow><munder><mi>min</mi><mrow><msub><mi>x</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>&Element;</mo><mi>&Omega;</mi></mrow></munder><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>q</mi><mo>*</mo></msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>C</mi><mo>_</mo><mi>ave</mi></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>利用罚函数法去除约束条件①和约束条件⑥,则式(4)可表示为式(5):<img file="FDA0000380334010000025.GIF" wi="1518" he="134" />其中:Ω′表示去掉约束条件①、约束条件②和约束条件⑥所剩下的约束条件所获得的可行域,<img file="FDA0000380334010000026.GIF" wi="40" he="41" />和<img file="FDA0000380334010000027.GIF" wi="45" he="41" />表示很大的正数,<maths num="0011"><![CDATA[<math><mrow><msub><mi>Z</mi><mn>1</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>max</mi><mo>{</mo><mn>0</mn><mo>,</mo><msup><mn>2</mn><mfrac><msub><mi>&gamma;</mi><mrow><mi>i</mi><mo>,</mo><mi>l</mi></mrow></msub><mi>w</mi></mfrac></msup><mo>-</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mrow><mi>D</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow><mrow><msub><mi>p</mi><mrow><mi>I</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mrow><mi>k</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo><mo>}</mo><mo>,</mo></mrow></math>]]></maths><maths num="0012"><![CDATA[<math><mrow><msub><mi>Z</mi><mn>2</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>max</mi><mo>{</mo><mn>0</mn><mo>,</mo><msub><mi>p</mi><mi>max</mi></msub><mo>-</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>}</mo><mo>;</mo></mrow></math>]]></maths>对式(5)进行整理得式(6):<maths num="0013"><![CDATA[<math><mrow><munder><mi>min</mi><mrow><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&Element;</mo><msup><mi>&Omega;</mi><mo>&prime;</mo></msup></mrow></munder><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><mi>A</mi></mrow></math>]]></maths><img file="FDA0000380334010000033.GIF" wi="894" he="605" />将式(5)分解成两层优化问题进行求解,该两层优化问题分别为第一层优化问题和第二层优化问题;其中第一层优化问题为功率控制问题,第二层优化问题为资源分配问题;(3)在解决第一层优化问题时,需要求解N<sub>d</sub>个子问题,即每组D2D对在所有可行资源上分别求f<sub>i</sub>(q)的最小值,q表示能效:<img file="FDA0000380334010000034.GIF" wi="888" he="320" />首先对f<sub>i</sub>(q)的函数求一阶导数并令导数为零,得到的最优功率点有式(7)所示三种可能情况:<img file="FDA00003803340100000411.GIF" wi="1854" he="624" />将三种可能最优功率点分别带入f<sub>i</sub>(q)的函数内,使得f<sub>i</sub>(q)取值最小的点记为最优功率点<img file="FDA0000380334010000042.GIF" wi="91" he="67" />对应的<img file="FDA0000380334010000043.GIF" wi="228" he="71" />即为f<sub>i</sub>(q)的最小值;(4)在解决第二层优化问题时,需要选出一组D2D对,使得该D2D对复用当前的资源块时,式(4)的值是最小的,即解决一个组合优化问题,等效数学模型如式(8):<maths num="0014"><![CDATA[<math><mrow><mi>max</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><msup><mi>M</mi><mo>&prime;</mo></msup><mo></mo></mrow></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mrow><mi>p</mi><mo>=</mo><msubsup><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>*</mo></msubsup></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>该数学模型包括如下约束条件:<maths num="0015"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>i</mi></mrow></math>]]></maths><maths num="0016"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>d</mi></msub></munderover><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>j</mi></mrow></math>]]></maths><maths num="0017"><![CDATA[<math><mrow><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&Element;</mo><mo>{</mo><mn>0,1</mn><mo>}</mo><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>,</mo><mo>&ForAll;</mo><mi>j</mi></mrow></math>]]></maths>采用启发式资源分配算法进行资源分配,确定每组D2D对所复用的资源。
地址 211189 江苏省南京市江宁区东南大学路2号