发明名称 一种用于OFDMA系统中的动态资源分配方法
摘要 本发明公开了一种用于OFDMA系统中的动态资源分配方法,包括:(1)参数初始化步骤;(2)子载波调度步骤;(3)功率分配步骤。本发明的方法无需大量逐步逼近搜索即可求解出最优功率,有效地降低了计算复杂度。当频谱效率达到系统要求后,此时能效作为主要考虑因素,通过本发明的方法可以求解出传输单位比特所耗能量最小,并能很好的取得频谱效率与能量效率的平衡折中,进一步提高能源利用率,使得通信更为节能环保。
申请公布号 CN102724761B 申请公布日期 2014.08.13
申请号 CN201210212902.2 申请日期 2012.06.23
申请人 电子科技大学 发明人 郑正广;冀秋月;但黎琳;武刚;李少谦
分类号 H04W72/04(2009.01)I;H04L27/26(2006.01)I 主分类号 H04W72/04(2009.01)I
代理机构 电子科技大学专利中心 51203 代理人 周永宏
主权项 一种用于OFDMA系统中的动态资源分配方法,包括:(1)参数初始化步骤;(2)子载波调度步骤;(3)功率分配步骤;其中,所述参数初始化步骤具体包括初始化如下参数:系统带宽B,子载波数目N,用户数目K,各用户要求的最小速率<img file="FDA0000516161050000011.GIF" wi="91" he="95" />各用户在各子信道的信道增益h<sub>k,n</sub>,其中,k=1,2,…,K;n=1,2,…,N,噪声功率σ<sup>2</sup>,目标误比特率P<sub>k,e</sub>,各用户允许的最大传输功率P<sub>k,max</sub>,各用户电路功率P<sub>k,ckt</sub>,信噪比差额Γ,单个子载波带宽W,各用户在各信道的信道增益噪声比γ<sub>k,n</sub>,功率放大器效率倒数α,及各个用户电路功率的静态功率P<sub>k,s</sub>,k=1,2,…,K;所述子载波调度步骤对所有子载波进行分配,分配给不同的用户,具体包括如下分步骤:(21).把每个用户已分配的子载波集合N<sub>k</sub>设为空集<img file="FDA0000516161050000012.GIF" wi="86" he="60" />未分配子载波集合为全集S={1,2,…,N},各用户已分配的子载波数目m<sub>k</sub>=0,未分配子载波的用户集合记为K<sub>E</sub>,临时用户集合为D={1,2,…,K};(22).选取每个用户最差的子信道,计算出各用户采用最大功率发射时所能达到的最大传输速率,记为<img file="FDA0000516161050000019.GIF" wi="96" he="63" />并将该传输速率作为参考速率;(23).选取临时用户集合D中所有用户在子载波集合S中信道质量最好的子载波序号,记为n<sub>k</sub>,并计算各用户最大传输功率发射所能达到的最大传输速率<img file="FDA0000516161050000013.GIF" wi="90" he="96" />(24).计算所有用户的速率增量<img file="FDA0000516161050000014.GIF" wi="310" he="94" />从临时用户集合D中选取速率增量最大的用户序号<img file="FDA0000516161050000015.GIF" wi="76" he="85" />选取步骤(23)中对应的子载波<img file="FDA0000516161050000016.GIF" wi="72" he="74" />分给用户<img file="FDA0000516161050000017.GIF" wi="73" he="85" />更新该用户已分配的子载波集合、各用户已分配的子载波数目、参考速率、未分配子载波集合、未分配子载波的用户集合;当剩余未分配子载波数目等于集合K<sub>E</sub>中用户数目时,D=K<sub>E</sub>;(25).重复步骤(23)至(24),直至所有子载波都分配完毕;所述功率分配步骤根据已分配的子载波集合,进行各用户功率分配,具体包括如下步骤:(31).计算各个用户在已分配子载波集合发送最小速率<img file="FDA0000516161050000018.GIF" wi="61" he="81" />时所需的最小功率P<sub>k,min</sub>;(32).计算出各个用户在已分配子载波集合上无功率限制要求及无速率约束条件下的功率<img file="FDA0000516161050000021.GIF" wi="86" he="78" /><img file="FDA0000516161050000022.GIF" wi="65" he="78" />的具体计算过程如下:求解单个用户的能效表达式u<sub>k,EE</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>u</mi><mrow><mi>k</mi><mo>,</mo><mi>EE</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>un</mi></mrow></msub><mrow><msub><mi>&alpha;P</mi><mrow><mi>k</mi><mo>,</mo><mi>un</mi></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>,</mo><mi>ckt</mi></mrow></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000516161050000023.GIF" wi="444" he="157" /></maths>其中,R<sub>k,un</sub>、P<sub>k,un</sub>分别表示用户k无最小速率要求及功率限制时所能达到的速率、功率;然后对能效进行求导,令导数为零,<img file="FDA0000516161050000024.GIF" wi="287" he="175" />μ<sub>k,un</sub>为用户k所对应的功率注水线;即可求解出用户k在已分配子载波集合上无功率限制且无速率约束的功率<img file="FDA0000516161050000025.GIF" wi="220" he="94" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>k</mi><mo>,</mo><mi>opt</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msup><mi>n</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><msub><mi>W</mi><mn>0</mn></msub><mo>[</mo><mrow><mo>(</mo><mfrac><msub><mi>P</mi><mrow><mi>k</mi><mo>,</mo><mi>s</mi></mrow></msub><mrow><msup><mi>n</mi><mo>&prime;</mo></msup><mi>&alpha;</mi></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><msup><mi>n</mi><mo>&prime;</mo></msup></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>n</mi><mo>&prime;</mo></msup></munderover><mfrac><mn>1</mn><msub><mover><mi>g</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub></mfrac><mo>)</mo></mrow><msup><mrow><mo>(</mo><munderover><mo>&prod;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>n</mi><mo>&prime;</mo></msup></munderover><msub><mover><mi>g</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><msup><mi>n</mi><mo>&prime;</mo></msup></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>]</mo><mo>+</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><msup><mi>n</mi><mo>&prime;</mo></msup></mfrac><mi>ln</mi><mrow><mo>(</mo><munderover><mo>&prod;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>n</mi><mo>&prime;</mo></msup></munderover><msub><mover><mi>g</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000516161050000026.GIF" wi="1520" he="205" /></maths>其中,W<sub>0</sub>(·)为Lambert函数,e为自然常数,n'为用户k已分配子载波集合中分配功率大于零的子载波数目,<img file="FDA0000516161050000027.GIF" wi="82" he="76" />为按照降序排列后的信道增益噪声比,n'满足如下条件<img file="FDA00005161610500000213.GIF" wi="1235" he="190" />其中,<img file="FDA00005161610500000210.GIF" wi="73" he="60" />表示正实数集,|N<sub>k</sub>|表示用户k已分配的子载波数目,求解满足条件的n',此时无最大传输功率限制及速率约束条件下的总传输功率为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>k</mi><mo>*</mo></msubsup><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>n</mi><mo>&Element;</mo><msub><mi>N</mi><mi>k</mi></msub></mrow></munder><msub><mi>p</mi><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><msup><mi>n</mi><mo>&prime;</mo></msup><msubsup><mi>&mu;</mi><mrow><mi>k</mi><mo>,</mo><mi>opt</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msup><mi>n</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><munder><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>&Element;</mo><msub><mi>N</mi><mi>k</mi></msub></mrow></munder><msup><mi>n</mi><mo>&prime;</mo></msup></munderover><mfrac><mn>1</mn><msub><mover><mi>g</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA00005161610500000211.GIF" wi="962" he="239" /></maths>(33).计算出各个用户满足频谱效率要求下的最优能效时的功率P<sub>k,opt</sub>(N<sub>k</sub>):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>k</mi><mo>,</mo><mi>opt</mi></mrow></msub><mrow><mo>(</mo><msub><mi>N</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>min</mi><mrow><mo>(</mo><mi>max</mi><mrow><mo>(</mo><msubsup><mi>P</mi><mi>k</mi><mo>*</mo></msubsup><mo>,</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA00005161610500000212.GIF" wi="1025" he="112" /></maths>
地址 611731 四川省成都市高新区(西区)西源大道2006号