下行NOMA系统中公平的功率分配方案

田心记，张丹青，李晓静

(1.河南理工大学物理与电子信息学院 河南 焦作 454150；2.黄河科技学院信息工程学院 郑州 450000)

随着移动通信的迅猛发展，传统的多址接入技术已难以满足无线数据业务量的爆炸式增长[1-2]。因此，第五代移动通信采用具有更高系统吞吐量和更高频谱效率的NOMA技术[3-5]。相较于传统多址接入技术在时域、频域和码域的研究，NOMA技术引入了一个新的维度-功率域[6]，在基站端为多个用户分配不同的功率，然后将这些用户的信号叠加在相同的时频资源上，用户接收到信号后采用串行干扰消除技术检测期望接收的信号[7-8]。功率分配不仅关系到各用户信号的检测次序，还影响到系统的可靠性和有效性，NOMA系统中的功率分配是近年的研究热点之一。

很多学者对单小区下行NOMA系统中的功率分配方案进行了研究，其中功率分配的目标有3种：最大化和速率、最大化能量效率以及最大化公平性。文献[9]研究了包含两用户的单个NOMA簇在用户最低速率需求约束下的功率分配方案，目标是最大化两个用户的和速率，给出了功率分配方案的闭式解。文献[10]将文献[9]中的两用户扩展到任意用户场景，在总功率约束和用户最低速率需求约束下，提出了最大化单输入单输出(single input single output,SISO)NOMA 系统和速率的功率分配方案。文献[11]研究了包含任意用户的单个NOMA簇中最大化能量效率的功率分配方案，文献[12]扩展了文献[11]的场景，提出了包含多个簇且每个簇包含任意用户的NOMA系统中最大化能量效率的功率分配方案。对于包含多个用户的单个NOMA簇，文献[13]以最大最小公平性为准则，在总功率和强用户最低速率需求约束下，提出了一种既能满足强用户最低速率需求又能最大化弱用户的最低速率的功率分配方案。文献[14]以最大公平性(maximin fairness,MMF)为准则，提出了一种能最大化用户最低速率的功率分配方案。该方案中所有用户的速率都相同，实现了用户在速率上的公平性。然而，该方案未考虑单个用户的速率需求，因此，有可能导致部分用户的速率高于该用户的速率需求，一部分用户的速率低于该用户的速率需求。

针对文献[14]的不足，提出了NOMA系统中公平地提高用户速率的功率分配方案。根据信道条件以及每个用户的最低速率需求计算每个簇所需的最低功率，在满足用户最低速率需求的基础上，推导了单个簇内两个用户提高的最低速率的最大值与该簇的总功率之间的关系，基于此建立满足用户最低速率需求的情况下公平地提高每个用户速率的功率分配优化问题，通过多次调整部分簇的功率得到在满足用户最低速率需求的基础上能最大化每个用户提高的最低速率的功率分配方案。仿真结果显示，当用户的最低速率需求不同时，所提方案中用户提高的最低速率和中断概率均优于MMF方案。

1 系统模型

包含1个基站和2K个用户的下行NOMA系统，如图1所示，基站和用户都配置单根天线。用户被分为K个簇，每个簇包含两个用户，分别用uk1和uk2表示第k个簇中的两个用户，k=1,2,···,K。假定uk1是近距离用户，uk2是远距离用户。基站到uk1和uk2的信道分别为hk1和hk2，|hk1|2≥|hk2|2。基站为第k个簇分配的功率为pk，其中uk1和uk2的功率分别为pk1和pk2，pk=pk1+pk2，pk1＜pk2。基站为每个簇分配一个子频段，簇间子频段正交。

分别用yk1和yk2表示uk1和uk2的接收信号，yk1和yk2的表达形式分别为：

式中，xk1和xk2分别是uk1和uk2的期望接收信号；nk1和nk2分别是uk1和uk2接收到的高斯白噪声，均值为零，方差为 σ2。

近距离用户uk1首先检测出uk2的期望接收信号xk2，并消除xk2对yk1的干扰。uk1译码xk2时的信干噪比(signal to interference and noise ratio,SINR)为：

uk2是远距离用户，可直接译码自身的期望接收信号xk2。uk2译码xk2时的SINR为：

uk1和uk2的单位带宽速率Rk1和Rk2分别为：

分别用rk1和rk2表示uk1和uk2所需的最低单位带宽速率。用P0表示满足所有用户最低速率需求所需的最低总功率，用Pmax表 示基站的总功率。Pmax≥P0时，基站的总功率能够满足所有用户的最低单位带宽速率需求。所提方案的目标是：在满足每个用户最低单位带宽速率需求的情况下，公平地提高每个用户的速率。

根据以上所述，Pmax≥P0时，功率分配的目标用公式表示为：

式中，C1表示基站的总功率为Pmax；C2表示uk1的单位带宽速率不低于rk1；C3表示uk2的单位带宽速率不低于rk2；C4表示uk1译码xk2时对SINR的要求。

2 功率分配方案

本节首先推导了满足所有用户最低速率需求时所需的最低总功率P0，然后再给出Pmax≥P0时公平地提高每个用户速率的功率分配方案。

2.1 最低总功率 P0的推导

首先推导每个簇所需的最低功率。由式(8)中的C2和C3可得：

由于a0是uk1正确译码xk2时对SINR的最低要求，因此a0=ak2即可。是x的单调递增函数且 |hk1|2≥|hk2|2，式(10)成立时，式(8)中的C4必定成立。因此，pk1和pk2满足式(9)和式(10)时，式(8)中C2、C3和C4成立。

用pk0表示第k个簇所需的最低功率，由式(9)和式(10)知pk0的取值为：

式(11)是单个簇满足用户最低速率需求时所需的最低功率，考虑到每个簇都有最低功率要求，因此满足所有用户的最低速率需求时所需的最低总功率P0的取值为：

2.2 功率分配

求解式(8)就能得到Pmax≥P0时的功率分配。若由式(8)直接求解功率，复杂度极高。

所提功率分配方案的优化目标是公平地提高每个用户的速率，文献[14]中MMF方案的优化目标是最大化系统的公平性，文献[15]的优化目标是最大化系统的能量效率。虽然所提功率分配方案与文献[14-15]的优化目标不同，然而三者的场景一致，故式(8)的求解可借鉴文献[14-15]的方法，即先求解第k个簇的总功率为pk时公平地提高该簇内两用户速率的功率分配，然后简化式(8)并求解簇间功率分配，最后根据簇间功率分配的结果为单个用户分配功率。

第k个簇的总功率为pk时，该簇内功率分配的目标用公式表式为：

式中，vki表示用户实际的单位带宽速率与所需的最低单位带宽速率的差值(即提高的速率)，vki=Rki−rki，i=1,2；C1表示该簇的总功率约束；C2和C3表示uki的 速率不低于rki时pki需要满足的条件。

在第k个簇的总功率pk保持不变的情况下，由于pk=pk1+pk2及vki=Rki−rki，i=1,2，增 大pk1时，vk1增大且vk2减小，增大pk2时，vk2增大且vk1减小。所以只有当vk1=vk2时，才能最大化min{vk1,vk2}，此时在满足用户最低单位带宽速率需求的基础上，两个用户提高的速率相等，这样就公平地提高了两个用户的速率。vk1=vk2等价于：

式(14)成立时，pk1与该簇总功率pk之间的关系为：

假定第k个簇的总功率为pk，MMF方案中该簇内两个用户的速率相等即Rk1=Rk2时，uk1的功率与该簇总功率pk之间的关系为：

根据式(15)和式(17)，可得：

当rk1＞rk2时，d=2rk1−rk2＞1，pk1＞；

当rk1=rk2时，d=2rk1−rk2=1，pk1=；

当rk1＜rk2时，d=2rk1−rk2＜ 1，pk1＜。

由文献[14]知，在单个簇的总功率保持不变时，为簇内强用户分配的功率越大，簇内两用户的和速率越高。因此，rk1＞rk2时，所提方案单簇和速率高于MMF方案；rk1=rk2时，所提方案单簇和速率等于MMF方案；rk1＜rk2时，所提方案单簇和速率低于MMF方案。

令vk=max{min{vk1,vk2}}，式(16)给出了单个簇内的vk与该簇总功率pk之间的关系。接下来简化式(8)并求解簇间功率分配。此时，式(8)可转化为：

式中，C1表示基站的总功率为Pmax；C2表示满足单个簇中用户最低速率需求时该簇的总功率需要满足的条件。当式(18)中的C2成立且按式(15)分配功率时，式(8)中C2、C3和C4必定成立。

式(8)要求解2K个用户的功率分配，而式(18)仅求解K个簇的功率分配，因此，式(18)是式(8)的一种简化表达形式。但是，此时仍无法直接给出式(18)的闭合解，为此接下来给出一种迭代的簇间功率分配方案。该方案的思路如下：首先为第k个簇分配满足用户最低速率需求的最低功率pk0，则是用于提高所有用户速率的功率，再将这些功率平均分配给每个簇作为该簇功率的初始值，然后开始迭代；在每次迭代过程中，计算第k个簇中用户提高的速率vk并组成向量 m v={vk}，找出mv中的最大元素对应的簇和最小元素对应的簇，分别用簇m和簇j表示，在满足该簇最低速率需求的条件下，减少第m个簇的功率同时增加第j个簇的功率，继续下次迭代并并计算得到mv；若本次迭代中的min{mv}大于等于前一次迭代中的min{mv}，则继续迭代，否则停止迭代，前一次迭代时的功率即为所提方案为每个簇分配的功率。迭代的簇间功率分配方案的具体步骤如下。

1) 根据信道条件以及用户所需的最低单位带宽速率计算第k个簇所需的最低功率pk0，令且pk=pk0+β，令minrateup=0且i=1，令矩阵是簇的总数，i表示迭代次数，U的每行用来存放每次迭代过程中的pk。

2)根据pk计算vk，若 min{mv}≥minrateup，令i=i+1、index=0且 minrateup=min{mv}，执行步骤3)及其后面的步骤，若 min{mv}＜minrateup，则U(i−1,k)就是所提方案为第k个簇分配的功率，即令pk=U(i−1,k)，无需执行步骤3)，停止迭代。

3)找出 min{mv}对应的簇，用簇j表示，令pj=pj+θ，0 ＜θ＜β。

4)找出 m ax{mv}对应的簇，用簇m表示，若pm−θ≥pm0，将pm−θ赋值给pm且令index=1，若pm−θ＜pm0，从mv中删除vm，再次执行该步骤，直到index[=1；]

5)将[p1,p2,···,pK]赋值给U的第i行，再次执行步骤2)。

minrateup表示前一次迭代过程中mv的最小元素值，min{mv}表示本次迭代过程中mv的最小元素值。若min{mv}＜minrateup，则表示前一次迭代时的功率能更公平地提高用户的速率，停止迭代；若min{mv}≥minrateup，则表示此次迭代的功率能更公平地提高用户的速率，继续进行功率调整以更公平地提高每个用户的速率，即执行步骤3)～5)。在步骤3)中找到min{mv}对应的簇，增加该簇的功率；在步骤4)中，找出 m ax{mv}对应的簇，用簇m表示，若将第m个簇的功率减少 θ后仍不低于该簇所需的最低功率，则将pm−θ赋值给pm，否则不对簇m的功率进行调整，并且从mv中删除vm，采用同样的方法找出可以减少功率的簇。

θ表示每次迭代过程中功率的调整量，index表示每次迭代过程中是否调整过功率，index=1表示已经调整过功率，可以进行下一次迭代。

采用上述方法得到第k个簇的功率pk后，为每个用户分配功率，其中，uk1的功率pk1与pk的关系如式(15)所示，uk2的功率为pk2=pk−pk1。

3 性能比较及仿真

本节仿真了所提功率分配方案的性能，并与文献[14]中的MMF方案进行了对比。假定单个基站服务了10个用户，每个簇中有2个用户，即簇数目K=5，信道服从独立的瑞利分布，高斯白噪声的均值为0、方差为1。由2.2节中的分析可知：单个簇的总功率相等的情况下，当rk1≥rk2时，所提方案的单簇和速率不低于MMF方案，当rk1＜rk2时，所提方案的单簇和速率低于MMF方案，故本节分别仿真了rk1＞rk2和rk1=rk2时，两种方案下系统提高的最低速率及中断概率。

图2仿真了rk1＞rk2时两种功率分配方案下用户提高的最低速率，即min{vk}。可以看出，参数相同时，所提方案用户提高的最低单位带宽速率高于 MMF 方案。信噪比(signal to noise ratio,SNR)范围为[10,20]dB时，两种参数下所提方案用户提高的最低速率比MMF方案高出了0.5 bps/Hz。原因在于所提方案在进行功率分配时考虑了每个用户的最低单位带宽速率需求，在满足用户最低单位带宽速率需求的基础上公平地提高了每个用户的速率，而MMF方案没有考虑用户的最低单位带宽速率需求。从图中还能看出，rk1=2，rk2=1时用户提高的最低速率高于rk1=3，rk2=2时用户提高的最低速率。原因在于当基站总功率Pmax保持不变时，一部分功率用来满足用户的最低速率需求，另一部分功率用于提高用户的速率，当用户的最低速率需求越低，用于提高用户速率的功率越高。

图3仿真了rk1和rk2相等时两种功率分配方案下用户提高的最低速率。可以看出，参数相同时，所提方案提高的最低速率与MMF方案相同。原因在于MMF方案功率分配的目标是最大化用户的最低速率，所提方案的目标是最大化用户提高的最低速率，当rk1=rk2时，两种功率分配方案的目标一致。从图3中还能看出，rk1=rk2=1时用户提高的最低速率高于rk1=rk2=2时用户提高的最低速率。原因在于当基站总功率Pmax保持不变时，用户的最低速率需求越高，P0越高，从而用于提高用户速率的功率Pmax−P0越低。因此，用户的最低速率需求越高，提高的速率越低。

图4仿真了rk1＞rk2时两种功率分配方案的中断概率。中断概率定义为：用户的速率Rki＜rki的概率，即P(Rki＜rki)=1−P(Rki≥rki)，i=1,2。仿真参数与图2中的相同，rk1=3且rk2=2时，所提方案用户uk1和用户uk2的门限值分别为3和2，MMF方案中两用户门限值均为2.5；rk1=2且rk2=1时，所提方案用户uk1和用户uk2的门限值分别为2和1，MMF方案两用户门限值均为1.5。MMF方案两用户的门限值之和与所提方案两用户的门限值之和相等。从图4能看出，参数相同时，所提方案的中断概率小于MMF方案，原因在于所提方案在满足用户最低速率需求的基础上公平地提高了每个用户的速率，MMF方案没有考虑用户的最低速率需求，当所有用户的速率都相同时，有可能部分用户的速率超出了该用户所需的速率，而另一部分用户的速率低于该用户所需的速率。从图4中还能看出，rk1=2且rk2=1时两种方案的中断概率均低于rk1=3且rk2=2时的中断概率。原因在于在其他条件相同的情况下，门限值越低，中断概率越低。

图5仿真了rk1和rk2相等时两种功率分配方案的中断概率。仿真参数与图3中的相同，rk1=rk2=2时，所提方案和MMF方案用户的门限值均为2；rk1=rk2=1时，所提方案和MMF方案的门限值均为1。从图5中能看出，参数相同时，所提方案的中断概率与MMF方案相同，原因在于当rk1和rk2相等时，两种功率分配方案的目标一致。从图中还能看出，rk1=rk2=1时两种功率分配方案的中断概率均低于rk1=rk2=2时的中断概率，原因如前所述。

4 结 束 语

本文研究了单小区下行NOMA系统中的功率分配方案。基于用户的最低速率需求，建立公平的提高用户速率的功率分配优化问题，简化该问题并用迭代算法进行求解。所提功率分配方案既满足了用户的最低速率需求，又公平地提高了用户的速率。由于所提功率分配方案仅考虑了每个NOMA簇包含两用户的场景，如何将该方案扩展到多簇且每个簇包含任意用户的场景有待于进一步研究。