发明名称 一种多无线电系统中基于公平性和精细化带宽分配的资源分配方法
摘要 本发明针对基于OFDMA(Orthogonal Frequency Division Multiple Access)的多载波系统中用户比例公平性和系统效率问题进行了研究,提出了一种联合资源分配方法,不仅保证了用户比例公平性下系统的吞吐量,还充分考虑了分配的带宽是子信道带宽整倍数的特点,对分配给每个终端的带宽进行子信道整数倍调整。最后通过仿真对比,从吞吐量和公平性两方面说明了本发明的优越性。
申请公布号 CN103888234A 申请公布日期 2014.06.25
申请号 CN201410079975.8 申请日期 2014.03.05
申请人 南京邮电大学 发明人 潘甦;曹跑跑;张磊
分类号 H04L5/00(2006.01)I;H04L27/26(2006.01)I;H04W72/04(2009.01)I 主分类号 H04L5/00(2006.01)I
代理机构 南京知识律师事务所 32207 代理人 胡玲
主权项 1.一种OFDMA多系统中基于公平性和精细化带宽分配的联合资源分配方法,其特征在于,所述方法包括以下步骤:步骤1、将用户端功率和网络端带宽建模为基于公平性的联合最优化模型:数据传输过程中,每一种NT都会为各个UE分配不同的带宽资源和发射功率,假设信道衰落是慢变的,即信道在一个资源分配周期内保持固定,再假设网络t分配给用户s的子信道带宽是连续的,则异构系统中UE<sub>s</sub>(s=1,2,3…S)的吞吐量可表示为<sup>[3]</sup>:<maths num="0001"><![CDATA[<math><mrow><msub><mi>R</mi><mi>s</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>r</mi><mi>st</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&eta;</mi><mi>st</mi></msub><mo>)</mo></mrow><msub><mi>&beta;</mi><mi>t</mi></msub><msub><mi>x</mi><mi>st</mi></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msup><mrow><mo>|</mo><mi>H</mi><mo>|</mo></mrow><mn>2</mn></msup><msub><mi>p</mi><mi>st</mi></msub></mrow><mrow><msub><mi>N</mi><mi>st</mi></msub><msub><mi>x</mi><mi>st</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中r<sub>st</sub>为UE<sub>s</sub>在NT<sub>t</sub>内的数据速率,η<sub>st</sub>(0≤ηst≤1)为UE<sub>s</sub>到NT<sub>t</sub>的平均误比特率<sup>[8]</sup>,其值为远远小于1的正数,β<sub>t</sub>(0<β<sub>t</sub><1)为NT<sub>t</sub>的系统效率<sup>[9]</sup>,其值和信道的编码方式有关,x<sub>st</sub>为网络NT<sub>t</sub>分配给UE<sub>s</sub>的带宽,p<sub>st</sub>为UE<sub>s</sub>在NT<sub>t</sub>的发射功率,|H<sub>st</sub>|<sup>2</sup>表示信道传输增益,N<sub>st</sub>表示此时的噪声功率谱密度。我们取β<sub>t</sub>=1,ηst=0,令<img file="FDA0000473077160000012.GIF" wi="232" he="156" />因此式(1)可简化为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>R</mi><mi>s</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>r</mi><mi>st</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn><mtext></mtext></mrow><mi>T</mi></munderover><msub><mi>x</mi><mi>st</mi></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mi>st</mi></msub><msub><mi>p</mi><mi>st</mi></msub></mrow><msub><mi>x</mi><mi>st</mi></msub></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>设置一种比例公平性方法,使得速率R<sub>s</sub>在不同的用户间在某种意义上是公平的,同时使得系统总速率之和在一定意义上最大,也就是说此方法同时兼顾了公平性和系统效率:将此问题写成最大化用户的传输数据速率的对数和形式,则可将最优化问题建模为如式(3)所示,其近似最优解易通过最优化方法中的梯度法得到,且其局部最优解即为全局最优解;<maths num="0003"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><mi>max</mi><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><mi>ln</mi><msub><mi>R</mi><mi>s</mi></msub><mo>=</mo><mi>max</mi><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><mi>ln</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>r</mi><mi>st</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><mi>max</mi><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><mi>ln</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>x</mi><mi>st</mi></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mi>st</mi></msub><msub><mi>p</mi><mi>st</mi></msub></mrow><msub><mi>x</mi><mi>st</mi></msub></mfrac><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths><img file="FDA0000473077160000015.GIF" wi="1475" he="467" />步骤2、忽略条件d的前提下,通过梯度法求得近似最优解,所述方法如下:最优化问题的拉格朗日函数为:<maths num="0004"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>L</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>p</mi><mo>;</mo><mi>&lambda;</mi><mo>,</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>x</mi><mi>st</mi></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mi>st</mi></msub><msub><mi>p</mi><mi>st</mi></msub></mrow><msub><mi>x</mi><mi>st</mi></msub></mfrac><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>&lambda;</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>t</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><msub><mi>x</mi><mi>st</mi></msub><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><msub><mi>&mu;</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>s</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>p</mi><mi>st</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><mo>[</mo><mi>ln</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>x</mi><mi>st</mi></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mi>st</mi></msub><msub><mi>p</mi><mi>st</mi></msub></mrow><msub><mi>x</mi><mi>st</mi></msub></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>&lambda;</mi><mi>t</mi></msub><msub><mi>x</mi><mi>st</mi></msub><mo>+</mo><msub><mi>&mu;</mi><mi>s</mi></msub><msub><mi>P</mi><mi>s</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>s</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>p</mi><mi>st</mi></msub><mo>]</mo></mtd></mtr><mtr><mtd><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>&lambda;</mi><mi>t</mi></msub><msub><mi>X</mi><mi>t</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中λ=[λ<sub>1</sub>,λ<sub>2</sub>,...λ<sub>m</sub>]是与接入点带宽分配约束相关的乘子向量,而μ=[μ<sub>1</sub>,μ<sub>2</sub>,...μ<sub>n</sub>]是与用户端发送功率约束相关的乘子向量;利用KKT条件<sup>[10]</sup>求得带宽分配与功率分配之间的关系:<maths num="0005"><![CDATA[<math><mrow><msub><mi>p</mi><mi>st</mi></msub><mo>=</mo><msub><mi>x</mi><mi>st</mi></msub><mo>&CenterDot;</mo><msup><mrow><mo>[</mo><mfrac><mn>1</mn><mrow><mi>ln</mi><mn>2</mn><mo>&CenterDot;</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>&mu;</mi><mi>S</mi></msub></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><msub><mi>g</mi><mi>st</mi></msub></mfrac><mo>]</mo></mrow><mo>+</mo></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中[a]<sup>+</sup>=max{a,0};本方法采用梯度法更新带宽<sup>[10]</sup>,选取合适的步长a则x<sub>st</sub>,<img file="FDA0000473077160000027.GIF" wi="83" he="48" />可收敛:<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>x</mi><mi>st</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mi>x</mi><mi>st</mi><mi>k</mi></msubsup><mo>-</mo><mi>a</mi><mfrac><mrow><mo>&PartialD;</mo><mi>L</mi></mrow><mrow><mo>&PartialD;</mo><msub><mi>x</mi><mi>st</mi></msub></mrow></mfrac><mo>]</mo></mrow><mo>+</mo></msup><mo>,</mo><mo>&ForAll;</mo><mi>s</mi><mo>,</mo><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>在已知x<sub>st</sub>,<img file="FDA0000473077160000028.GIF" wi="82" he="46" />下,根据式(5),则可求得p<sub>st</sub>,<img file="FDA0000473077160000029.GIF" wi="106" he="47" />同样采用梯度法更新拉格朗日乘子λ和μ,如下:<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>&lambda;</mi><mi>t</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mi>&lambda;</mi><mi>t</mi><mi>k</mi></msubsup><mo>-</mo><mi>b</mi><mrow><mo>(</mo><msub><mi>X</mi><mi>t</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><msub><mi>x</mi><mi>st</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mo>+</mo></msup><mo>,</mo><mo>&ForAll;</mo><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>&mu;</mi><mi>s</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><msup><mrow><mo>[</mo><msubsup><mi>&mu;</mi><mi>s</mi><mi>k</mi></msubsup><mo>-</mo><mi>c</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>s</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>p</mi><mi>st</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mo>+</mo></msup><mo>,</mo><mo>&ForAll;</mo><mi>s</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中k为迭代轮次,b和c为迭代步长。通过选取合适的迭代步长理论上可以保证收敛;由式(5)我们可以看到,功率分配p<sub>st</sub>、带宽分配x<sub>st</sub>和用户当前速率R<sub>s</sub>之间相互关联,因此需要采用如下的迭代方法求解;R<sub>s</sub>更新如下:<maths num="0009"><![CDATA[<math><mrow><msubsup><mi>R</mi><mi>s</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>R</mi><mi>s</mi><mi>m</mi></msubsup><mo>-</mo><msup><mi>d</mi><mi>m</mi></msup><mrow><mo>(</mo><msubsup><mi>R</mi><mi>s</mi><mi>m</mi></msubsup><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><mi>ln</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>x</mi><mi>st</mi></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mi>st</mi></msub><msub><mi>p</mi><mi>st</mi></msub></mrow><msub><mi>x</mi><mi>st</mi></msub></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中m为速率迭代次数,d是UE<sub>m</sub>速率更新迭代步长。通过选取合适的迭代步长可以保证速率的收敛性;步骤3、根据步骤2求出的近似最优解,对所有求得的带宽χ<sub>st</sub>取向下的最大子信道带宽整数倍:<img file="FDA0000473077160000031.GIF" wi="323" he="153" />其中d<sub>t</sub>代表网络NT<sub>t</sub>的子信道带宽,然后对于每一个NT,统计剩余带宽资源<img file="FDA0000473077160000032.GIF" wi="394" he="122" />为了保证公平性,对于每一个NT,将剩余带宽资源按照子信道带宽为单元分配给带宽最小的UE<sub>s</sub>:For NT t=1:T<maths num="0010"><![CDATA[<math><mrow><mi>While</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mi>t</mi><mo>*</mo></msubsup><mo>></mo><msub><mi>d</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>选取带宽最小的UE<sub>s</sub>,分配一个dt给它,更新UE<sub>s</sub>带宽以及剩余带宽资源;End WhileEnd For步骤4、为了衡量用户比例公平性,引入Jain公平性因子(Fairness index,FI)<sup>[11]</sup>来表示方法的公平性,FI定义为:<maths num="0011"><![CDATA[<math><mrow><mi>FI</mi><mo>=</mo><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><msub><mi>R</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><mrow><mo>(</mo><mi>S</mi><mo>&CenterDot;</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><msubsup><mi>R</mi><mi>s</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>
地址 210046 江苏省南京市亚东新城区文苑路9号