发明名称 一种多小区中继OFDMA系统中分布式资源分配方法
摘要 本发明是提供一种针对带有中继协作的OFDMA多小区系统中的低复杂度分布式资源分配方法。该方法主要包括自适应载波分配和注水功率分配,载波分配和功率分配交替迭代进行,直到算法收敛。在功率分配阶段,给出了一种迭代多平面注水法,通过调节不同的注水平面实现容量的最大化。在系统的公平性方面,可以实现对任意等级用户的灵活覆盖,并可以根据用户的优先等级提供差异化服务。本发明给出的方法收敛性快,并且系统容量得到提高。
申请公布号 CN103220114A 申请公布日期 2013.07.24
申请号 CN201310146746.9 申请日期 2013.04.24
申请人 南京邮电大学 发明人 朱琦;吕林输
分类号 H04L5/00(2006.01)I;H04W16/14(2009.01)I 主分类号 H04L5/00(2006.01)I
代理机构 江苏爱信律师事务所 32241 代理人 唐小红
主权项 1.一种多小区中继OFDMA系统中分布式资源分配方法,其特征在于包括以下步骤:a.初始化最大外层循环次数I<sub>max</sub>;为第l小区的第k用户设置优先权值<img file="FDA00003098280100011.GIF" wi="72" he="91" /><img file="FDA00003098280100012.GIF" wi="195" he="86" />初始化外层循环变量i=0;初始化所有基站功率变量<img file="FDA00003098280100013.GIF" wi="84" he="98" />和中继功率变量<img file="FDA00003098280100014.GIF" wi="394" he="95" />为0,其中:m是小区内中继的编号,n是小区内子载波的编号;b.初始化小区编号l=1;c.对任意子载波n计算链路(m,k)<img file="FDA00003098280100015.GIF" wi="200" he="81" />的信道容量,并从中选出具有最大加权数据速率的链路<img file="FDA00003098280100016.GIF" wi="144" he="111" />如下式:<maths num="0001"><![CDATA[<math><mrow><msub><mrow><mo>(</mo><mover><mi>m</mi><mo>^</mo></mover><mo>,</mo><mover><mrow><mi>k</mi><mo>,</mo></mrow><mo>^</mo></mover><mo>)</mo></mrow><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><munder><mi>arg</mi><mrow><mi>m</mi><mo>&Element;</mo><msubsup><mi>&Phi;</mi><mi>l</mi><mi>r</mi></msubsup><mo>,</mo></mrow></munder><munder><mi>max</mi><mrow><mi>k</mi><mo>&Element;</mo><msubsup><mi>&Phi;</mi><mi>l</mi><mi>u</mi></msubsup></mrow></munder><mo>{</mo><msubsup><mi>w</mi><mi>l</mi><mi>k</mi></msubsup><mo>&CenterDot;</mo><mi>min</mi><mo>[</mo><msubsup><mi>R</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><mo>,</mo><msubsup><mi>R</mi><mrow><mi>m</mi><mo>,</mo><mi>k</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow></msubsup><mo>]</mo><mo>}</mo><mo>,</mo><mo>&ForAll;</mo><mi>n</mi></mrow></math>]]></maths>当有多于1条链路的加权数据速率最大时,则选取第1条最大链路;其中:<img file="FDA00003098280100018.GIF" wi="80" he="88" />为链路基站l到中继m在子载波n上的数据速率,即:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>R</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msubsup><mi>p</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><msup><mrow><mo>|</mo><msubsup><mi>H</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>l</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>I</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>其中:<img file="FDA000030982801000110.GIF" wi="206" he="114" />为小区l的基站到该小区中继m在第n子载波上的信道增益,<img file="FDA000030982801000111.GIF" wi="103" he="91" />为小区l基站在子载波n上为中继m分配的功率,σ<sup>2</sup>为子载波n上的噪声功率,<img file="FDA000030982801000112.GIF" wi="70" he="85" />为小区l第m中继测量得到的来自其他小区在子载波n上的干扰功率;<img file="FDA000030982801000113.GIF" wi="103" he="91" />为链路中继m到用户k在子载波n上的数据速率:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>R</mi><mrow><mi>m</mi><mo>,</mo><mi>k</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow></msubsup><mo>=</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msubsup><mi>p</mi><mrow><mi>m</mi><mo>,</mo><mi>k</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow></msubsup><msup><mrow><mo>|</mo><msubsup><mi>H</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>I</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow><mi>n</mi></msubsup></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>其中:<img file="FDA000030982801000115.GIF" wi="206" he="111" />为l小区的中继m到该小区用户k在第n子载波上的信道增益,<img file="FDA000030982801000116.GIF" wi="90" he="99" />为l小区中继m在子载波n上为用户k分配的功率,σ<sup>2</sup>为子载波n上的噪声功率,<img file="FDA00003098280100021.GIF" wi="80" he="91" />为l小区第k用户测量得到的来自其他小区在子载波n上的干扰功率;d.采用多平面迭代注水法完成功率分配,步骤如下:d1.初始化最大内层循环次数S<sub>max</sub>;设置循环初始变量s=0;随机初始化基站功率对偶变量λ<sub>l</sub>和中继功率对偶变量<img file="FDA00003098280100022.GIF" wi="264" he="82" />d2.测量当前链路在子载波n上的噪声功率σ<sup>2</sup>;测量中继m在子载波n上受到的同频干扰功率<img file="FDA00003098280100023.GIF" wi="104" he="89" />测量用户k在子载波n上受到的干扰功率<img file="FDA00003098280100024.GIF" wi="96" he="106" />测量第一跳链路基站l到中继m在子载波n上的信道增益<img file="FDA00003098280100025.GIF" wi="238" he="114" />测量第二跳链路中继m到用户k在子载波n上的信道增益<img file="FDA00003098280100026.GIF" wi="232" he="122" />并计算广义信道质量比:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&eta;</mi><mi>l</mi><mi>n</mi></msubsup><mo>=</mo><mfrac><mrow><msup><mrow><mo>|</mo><msubsup><mi>H</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>l</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mrow><mo>(</mo><msubsup><mi>I</mi><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow><mi>k</mi></msubsup><mo>+</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mrow><msup><mrow><mo>|</mo><msubsup><mi>H</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mrow><mo>(</mo><msubsup><mi>I</mi><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow><mi>m</mi></msubsup><mo>+</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>d3.计算基站l在子载波n上为中继m分配的最优功率:<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>p</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msup><mrow><mo>[</mo><mfrac><msubsup><mi>w</mi><mi>l</mi><mi>k</mi></msubsup><mrow><mi>ln</mi><mn>2</mn><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>l</mi></msub><mo>+</mo><msubsup><mi>&mu;</mi><mi>l</mi><mi>m</mi></msubsup><msubsup><mi>&eta;</mi><mi>l</mi><mi>n</mi></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mrow><msubsup><mi>I</mi><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow><mi>m</mi></msubsup><mo>+</mo><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow><msup><mrow><mo>|</mo><msubsup><mi>H</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><mo>|</mo></mrow><mn>2</mn></msup></mfrac><mo>]</mo></mrow><mo>+</mo></msup></mrow></math>]]></maths>上式是一种多平面注水(multilevel water-filling)形式,其中函数x<sup>+</sup>=max(0,x),进一步求出中继m在子载波n上为用户k分配的功率<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>p</mi><mrow><mi>m</mi><mo>,</mo><mi>k</mi></mrow><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow></msubsup><mo>=</mo><msubsup><mi>p</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><msubsup><mi>&eta;</mi><mi>l</mi><mi>n</mi></msubsup><mo>;</mo></mrow></math>]]></maths>d4.选取适当的步长ε(s),为了迭代能够收敛,需要使步长随迭代次数s减小,如可以选择步长为迭代次数s的倒数,ε(s)=1s;然后沿次梯度方向更新对偶变量:<maths num="0007"><![CDATA[<math><mrow><msub><mi>&lambda;</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>&lambda;</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mi>&epsiv;</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>p</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msub><mi>P</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mo>+</mo></msup></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>&mu;</mi><mi>l</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>s</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mi>&mu;</mi><mi>l</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mi>&epsiv;</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>&Sigma;</mi><mrow><mi>n</mi><mo>&Element;</mo><msubsup><mi>&Psi;</mi><mi>l</mi><mi>m</mi></msubsup></mrow></msub><msubsup><mi>p</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow><mi>n</mi></msubsup><msubsup><mi>&eta;</mi><mi>l</mi><mi>n</mi></msubsup><mo>-</mo><msub><mi>P</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mo>+</mo></msup></mrow></math>]]></maths>式中:P<sub>B</sub>为基站的最大发射功率,P<sub>R</sub>为中继的最大发射功率;d5.更新循环变量s=s+1;d6.判断是否终止迭代,如达到功率收敛条件或已达到最大迭代次数S<sub>max</sub>,则终止迭代,完成功率分配过程,转步骤e;否则,返回步骤d3;所述达到功率收敛条件的判断标准为对偶变量的相对增加量不大于收敛指标δ<sub>1</sub>,max{△λ<sub>l</sub>λ<sub>l</sub>,△μ<sub>l</sub>μ<sub>l</sub>}≤δ<sub>1</sub>,所述达到最大迭代次数的判断标准为s≥S<sub>max</sub>;e.当前小区l资源分配已完成,继续进行下一个小区的分配,更新小区编号l=l+1,转步骤f;f.判断当前小区编号l是否大于小区总数L,若是,则已完成一轮所有小区的资源分配,应结束本轮分配,转步骤g;否则,转步骤c,继续对当前小区l进行子载波和功率分配;g.更新外层循环变量i=i+1;h.判断外层循环是否终止:如达到容量收敛条件或已达到最大迭代次数I<sub>max</sub>,则终止外层循环,完成本方法的分配过程;否则,返回步骤b;所述达到容量收敛条件的判断标准为所有小区的总容量的增量不大于收敛指标δ<sub>2</sub>,所述达到最大迭代次数的判断标准为i≥I<sub>max</sub>。
地址 210003 江苏省南京市新模范马路66号