单天线下行NOMA系统中的功率分配方法与流程

本发明属于通信领域，尤其涉及一种单天线下行noma系统中的功率分配方法。

背景技术：

面对海量数据业务的需求，下一代移动通信系统要求传输速率更高、频谱效率更高、接入设备数量更多、时延更低。面对未来更高的挑战，业界纷纷提出了多种新型多址接入技术，非正交多址接入(non-orthogonalmultipleaccess,noma)技术就是其中的一种。noma在发送端引入了功率域复用，将多个用户的信号复用到同一时频资源上。接收端利用串行干扰消除(successiveinterferencecancellation,sic)方法消除部分多用户干扰进而检测自身的期望接收信号。根据noma的实现原理，信道越差的用户的功率越高，按照用户信道条件由差到好的顺序依次解码。因此，在发送端需要为用户分配功率，以获得最大的性能增益，同时达到接收端采用sic消除部分多址干扰的目的。因此，noma中的功率分配是近年的研究热点之一。

传统的功率分配方案分为：全空间搜索功率分配方案、固定功率分配方案、分数阶功率分配方案、平均功率分配方案以及注水功率分配方案。后来，有学者提出了按照某种原则进行功率分配，如最大化速率、最大化能量效率及最大公平。文献“onoptimalpowerallocationfordownlinknon-orthogonalmultipleaccesssystems”以系统总功率作为约束条件，提出了最大化系统能量效率的功率分配方法。然而，该方法仅限于每个簇包含两用户的场景，不适用于每个簇包含多个用户的场景。

技术实现要素：

本发明提出了单天线下行noma系统中的功率分配方法，适用于包括1个基站和mk个用户的单天线下行noma系统且用户和基站都配置单天线。

实现本发明的技术思路是：基站将用户分簇并为用户簇分配正交的子频段，基站根据信道条件计算满足串行干扰消除时每个簇所需的最低总功率以及系统所需的最低总功率，以每个簇所需的最低总功率作为约束条件，构建最大化单个簇内能量效率的功率分配优化问题，求得单个簇的能量效率，然后以系统所需的总功率作为约束条件，构建最大化系统能量效率的簇间功率分配优化问题，采用迭代的方法求解簇间功率分配，最后根据簇间功率分配的结果为簇内的单个用户分配功率。

综上所述，一种单天线下行noma中的功率分配方法，适用于包括1个基站和mk个用户的单天线下行noma系统且用户和基站都配置单天线，包括如下步骤：

a，基站将mk个用户分为k个簇，每个簇包含m个用户，基站为每个簇分配一个子频段，簇间子频段正交；

b，用ukm表示第k个簇中的第m个用户，k＝1,2,…,k，m＝1,2,…,m，基站到ukm的信道为hkm，|hk1|²≥|hk2|²≥…≥|hkm|²，用r0表示正确检测信号时对信干噪比(signaltointerferenceandnoiseratio,sinr)的最低要求，基站计算满足串行干扰消除(successiveinterferencecancellation,sic)时ukm所需的最低功率，用pkm0表示，m＝1时，m＝2,3,4,…,m时，σ²是用户接收到的噪声方差，k是簇的总数，m是每个簇中包含的用户数；

c，基站计算满足sic时第k个簇所需的最低总功率pk0以及系统所需的最低总功率pmin，pk0的表达形式为

pmin的表达形式为

其中，k＝1,2,…,k，m＝1,2,…,m，k是簇的总数，m是每个簇中包含的用户数；

d，基站计算第k个簇的总功率为pk时为ukm分配的功率pkm以及此时的能量效率，k＝1,2,…,k，m＝1,2,…,m；

e，用pmax表示基站的总功率，令pmax≥pmin，构建最大化系统能量效率的簇间功率分配优化问题，

其中，约束条件c1表示系统的总功率为pmax，约束条件c2表示单个簇的总功率不能低于该簇所需的最低总功率；

f，求解步骤e中的优化问题，得到簇间功率分配pk，k＝1,2,…,k；

g，将步骤f得到的pk带入步骤d中的式(1)得到为第k个簇的第m个用户分配的功率pkm，k＝1,2,…,k，m＝1,2,…,m，k是簇的总数，m是每个簇中包含的用户数。

进一步，所述步骤d具体包括：

d1，构建第k个簇的总功率为pk时最大化该簇内用户速率之和的优化问题，

其中，约束条件c1表示该簇的总功率不能低于该簇所需的最低总功率，约束条件c2表示满足sic时对单个用户的功率要求；

d2，求解步骤d1中的优化问题，得到一组次优解，

其中，

d3，基站计算得到该簇内的所有用户的速率之和基站计算该簇的能量效率其中，

进一步，所述步骤f具体包括：

f1，引入函数q(pk,k＝1,2,…,k,α)，表达式如下，

其中，α是一个正数；

f2，构造拉格朗日函数f(pk,k＝1,2,…,k,λ)，

其中，λ是拉格朗日乘子；

f3，求f(pk,k＝1,2,…,k,λ)关于pk和λ的导数，得到如下方程组，

f4，求解步骤f3中的方程组，得到

其中，λ的取值满足式(3)中，如果则否则pk＝pk0；

f5，设置α的初始值为0，设置q(pk,k＝1,2,…,k,α)的初始值为∞且参数δ>0；

f6，由式(3)计算得到pk，k＝1,2,…,k；

f7，将pk和α代入式(2)计算得到q(pk,k＝1,2,…,k,α)；

f8，更新α，令α＝η(pk,k＝1,2,…,k)；

f9，如果|q(pk,k＝1,2,…,k,α)|>δ，则执行步骤f6、步骤f7和步骤f8，直到|q(pk,k＝1,2,…,k,α)|≤δ。

有益效果：

本发明公开的方法将最大化能量效率的功率分配方案扩展到每个簇中包含任意用户的场景，推导了在满足用户速率需求时每个簇所需的最低总功率以及系统所需的最低总功率，在系统所需的最低总功率约束下给出了最大化系统能量效率的功率分配方案。该方法不仅能最大化系统能量效率，还满足了每个用户的速率需求。

附图说明

图1是本发明实施例的系统模型；

图2是本发明的流程图。

具体实施方式

下面给出本发明的一种实施例，对本发明做进一步详细的说明。

如图1所示，考虑包含1个基站和mk个用户的下行noma系统，基站和用户都配置单根天线。用户被分为k个簇，每个簇包含m个用户，用ukm表示第k个簇中的第m个用户，k＝1,2,…,k，m＝1,2,…,m。基站到ukm的信道为hkm，|hk1|²≥|hk2|²≥…≥|hkm|²。基站为第k个簇分配的总功率为pk，其中ukm的功率为pkm，pk1≤pk2≤…≤pkm，基站为每个簇分配一个子频段，簇间子频段正交。

用ykm表示ukm的接收信号，ykm的表达形式为

其中，xkm是ukm的期望接收信号，nkm是ukm接收到的高斯白噪声，均值为零方差为σ²。

uk1进行串行干扰消除(successiveinterferencecancellation,sic)，即首先检测出xkm，并消除该信号对yk1造成的干扰，然后再检测xk(m-1)，并消除该信号对yk1造成的干扰，依次检测其他信号并消除这些信号对yk1造成的干扰，直至检测出xk1。uk1检测xkm时的信干噪比(signaltointerferenceandnoiseratio,sinr)为

同理，ukj检测xkm时的sinr为

其中，j≤m，m＝1,2,…,m，j＝1,2,…,m。

假定r0是正确检测信号时对sinr的最低要求，为了执行sic，ukj检测xkm时的sinr必须不低于r0，因此，要求下式成立

由此可推出，pkm的取值满足

令l(|hkj|²)是|hkj|²的单调递减函数。由于|hk1|²≥|hk2|²≥…≥|hkm|²，当j＝m时，l(|hkj|²)达到最大值，即式(5)可化为

m＝1时，pk1的取值范围为

m＝2时，pk2的取值范围为

m＝3时，pk3的取值范围为

m＝4时，pk4的取值范围为

采用归纳法可得，m＝2,3,…,m时，pkm的取值满足下式

用pkm0表示ukm进行sic以及正确检测期望信号所需的最低功率。m＝1时，m＝2,3,…,m时，pkm0的取值为

由式(12)可得出，当i<j时，pki0<pkj0，即在同一个簇内，用户的信道增益越低，所需的最低功率越高。用pk0表示第k个簇内的所有用户进行sic以及正确检测期望信号所需的最低总功率，则pk0的取值为

若第k个簇的总功率低于pk0，则不能保证该簇内sic的顺利执行，从而无法正确检测所有的期望信号。用pmin表示满足sic以及正确检测期望信号时系统所需的最低总功率，pmin的取值为

用rkj表示ukj的单位带宽速率，rkj表示为

第k个簇内所有用户的单位带宽速率之和为

系统中mk个用户的单位带宽速率之和为

用pmax表示基站的总功率，假设pmax≥pmin，否则无法实现sic。不考虑每个用户的速率约束时，最大化系统能量效率的功率分配的目标函数用公式表示为

其中，约束条件c1表示系统的总功率为pmax，约束条件c2表示单个簇的总功率不能低于该簇所需的最低总功率，约束条件c3用于保证sic的顺利执行。不考虑用户的速率约束并不意味着对每个用户的功率和每个簇的总功率没有要求，否则无法执行sic，更无法检测期望接收的信号。约束条件c3是对单个用户功率的要求，约束条件c2是对单个簇的总功率的要求。

式(18)中的优化问题是非凸优化问题，直接求解的复杂度极高。为此，先求解第k个簇的总功率为pk时，最大化该簇内能量效率的功率分配以及此时的能量效率，然后简化式(18)中的优化问题，求解系统的总功率为pmax时最大化系统能量效率的簇间功率分配，最后基于簇间的功率分配为每个用户分配功率。

先求解第k个簇的总功率为pk时，最大化该簇内能量效率的功率分配，用公式表示为

其中，约束条件c1表示簇的总功率为pk，约束条件c2表示该簇的总功率不能低于该簇所需的最低总功率，约束条件c3用于保证sic的顺利执行。

在第k个簇的总功率为pk时，只有最大化该簇内用户的速率之和，才能最大化该簇的能量效率。接下来求最大化该簇内用户速率之和的功率分配，如式(20)所示，

可采用拉格朗日方法求解式(20)，与文献“dynamicuserclusteringandpowerallocationforuplinkanddownlinknon-orthogonalmultipleaccess(noma)systems”中第v节的方法类似，然而，该方法很复杂，m＝4时，功率分配分8种情况，没有固定的表达式，不利于接下来求解簇间的功率分配。此处给出一种次优的求解方法，该方法使得每个簇内第2个用户至第m个用户刚好满足sic以及正确检测期望信号的需求，同时尽可能地提高第一个用户的速率。

由前面的分析知，式(21)成立时，

第k个簇内的第2个用户至第m个用户刚好满足sic以及正确检测期望信号的需求。将式(21)带入可得

由式(22)可得

其中，

第k个簇的总功率为pk且按照式(21)和式(23)为单个用户分配功率时，该簇内的所有用户的速率之和为该簇的能量效率为其中，这样就求解出了单个簇内最大化能量效率的功率分配。基于该结果，将式(18)简化为

(24)

其中，式(24)中的约束条件c1等价于式(18)中的约束条件c1，式(22)中的约束条件c2等价于式(18)中的约束条件c2，当pk≥pk0且按照式(21)和式(23)为单个用户分配功率时，式(18)中的约束条件c3必定成立，因此式(24)中没有约束条件c3。式(18)考虑mk个用户间的功率分配，式(22)只需要考虑k个簇之间的功率分配，这样就将系统内用户间的功率分配简化为簇间的功率分配，简化了式(18)中的优化问题。

式(24)中的最大化是非凸优化问题。为解决该问题，引入如下目标函数，

其中，α是一个正数。在α确知的情况下，考虑如下凸优化问题，

(26)

有文献“onoptimalpowerallocationfordownlinknon-orthogonalmultipleaccesssystems”知，如果pk^*(α)是式(26)的最优解，k＝1,2,…,k，且q^*(pk,k＝1,2,…,k,α)是式(26)中目标函数的最大取值，则只有当q^*(pk,k＝1,2,…,k,α)＝0时，pk^*(α)是式(24)中优化问题的最优解。因此，求解式(26)中的优化问题，并反复更新α使得q^*(pk,k＝1,2,…,k,α)＝0就能求得式(24)中优化问题。

假定式(26)中的α是一个固定的值。根据式(26)，构造拉格朗日函数f(pk,k＝1,2,…,k,λ)，

其中，λ是拉格朗日乘子。求f(pk,k＝1,2,…,k,λ)关于pk和λ的导数，得到如下方程组，

由式(28)可得

其中，λ的取值满足上式中，如果则否则pk＝pk0。

式(29)中的pk就是式(26)的最优解，找到使得q(pk,k＝1,2,…,k,α)＝0的α就等价于找到了(24)中优化问题的最优解，采用如下迭代算法找到α的取值，迭代算法的具体步骤为：

初始化：设置α的初始值为0，设置q(pk,k＝1,2,…,k,α)的初始值为∞且参数δ>0。while|q(pk,k＝1,2,…,k,α)|>δ

由式(29)计算得到pk，k＝1,2,…,k；

将pk和α代入式(25)计算得到q(pk,k＝1,2,…,k,α)；

更新α，令α＝η(pk,k＝1,2,…,k)；

end

上述迭代算法得到的α就是最大的能量效率。按照上述迭代算法得到pk后，再根据式(21)和式(23)得到pkm，pkm就是为ukm分配的功率。

结合本发明的流程图即图2，单天线noma系统中最大化能量效率的功率分配方法的具体步骤如下：

a，基站将mk个用户分为k个簇，每个簇包含m个用户，基站为每个簇分配一个子频段，簇间子频段正交；

b，用ukm表示第k个簇中的第m个用户，k＝1,2,…,k，m＝1,2,…,m，基站到ukm的信道为hkm，|hk1|²≥|hk2|²≥…≥|hkm|²，用r0表示正确检测信号时对信干噪比(signaltointerferenceandnoiseratio,sinr)的最低要求，基站计算满足串行干扰消除(successiveinterferencecancellation,sic)时ukm所需的最低功率，用pkm0表示，m＝1时，m＝2,3,4,…,m时，σ²是用户接收到的噪声方差，k是簇的总数，m是每个簇中包含的用户数；

c，基站计算满足sic时第k个簇所需的最低总功率pk0以及系统所需的最低总功率pmin，pk0的表达形式为

pmin的表达形式为

其中，k＝1,2,…,k，m＝1,2,…,m，k是簇的总数，m是每个簇中包含的用户数；

d，基站计算第k个簇的总功率为pk时为ukm分配的功率pkm以及此时的能量效率，k＝1,2,…,k，m＝1,2,…,m；

e，用pmax表示基站的总功率，令pmax≥pmin，构建最大化系统能量效率的簇间功率分配优化问题，

其中，约束条件c1表示系统的总功率为pmax，约束条件c2表示单个簇的总功率不能低于该簇所需的最低总功率；

f，求解步骤e中的优化问题，得到簇间功率分配pk，k＝1,2,…,k；

g，将步骤f得到的pk带入步骤d中的式(1)得到为第k个簇的第m个用户分配的功率pkm，k＝1,2,…,k，m＝1,2,…,m，k是簇的总数，m是每个簇中包含的用户数。

以上实施例仅仅是对本发明的举例说明，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。