一种CF-mMIMO下的用户分簇及功率分配方法

摘要

一种CF‑mMIMO下的用户分簇及功率分配方法，涉及通信技术领域,是为了解决目前的用户分簇方法的计算量较大导致的不易于实现，以及用户的功率分配方法还需要提供提高服务的公平性的问题。本方法通过最大化所有用户的最小SINR，从而最大化接入网络的用户，最小化用户的掉线率。不需要知道用户的具体位置信息就可以对用户进行分簇，该分簇方法的计算量较低，更易于实现；对于功率控制系数分配，该算法保证了最大化簇内所有用户的最小SINR，保证了用户服务的公平性以及用户的接入。

权利要求

1.一种CF‑mMIMO下的用户分簇及功率分配方法，在Cell‑free Massive MIMO下：所有AP在相同的时频资源下同时向所有用户提供服务；其特征是：Cell‑free Massive MIMO下最大用户信干噪比的用户分簇及功率分配方法的具体方法为：步骤一、初始化：设置初始分簇数I0、最大分簇数Imax；

步骤二、执行第l+1次迭代：对所有用户进行分簇，分成I0+l个用户簇，对每个用户簇均通过SCA算法计算功率控制系数，得到最小的用户信干噪比Γmin；

步骤三、判断步骤二获得的最小的用户信干噪比Γmin用户是否大于或等于预设的停止迭代条件Γmin≥γth，如果判断结果为是，则执行步骤四；如果判断结果为否且I0+l＝Imax，则执行步骤五；

步骤四、令l＝l+1；返回执行步骤二；

步骤五、.用户检测及排除步骤；

令Ei表示用户ui的分簇参考变量，它通过如下公式获得：式中：m和m′的物理含义相同；

ui表示第i个用户；

M表示AP的总数量，m表示第m个AP；

表示第m个AP到用户u i的大尺度衰落系数；

表示第m个AP与用户u i的大尺度衰落系数的信道状态估计值；

L表示每个AP上的天线数量；

ρd表示每个AP的最大归一化发射功率即：由噪声功率N0归一化；

由于已知分簇数为Il＝I0+l，则每个簇的最小用户数Z1为不大于用户数比分簇数的最大取整数，通过解下面的凸优化问题得到每个用户的具体的分簇结果式中：Ei表示用户ui的分簇参考变量；

表示具体的分簇结果， 表示用户ui属于用户簇nq，否则表示为用户集， 表示第i个用户；

表示用户簇集，nq表示隶属于 的第q个簇；

Smax衡量指标，表示信干噪比最大的用户簇的近似信干噪比；

Z1每个簇的最小用户数；

1 1

步骤六、初始化：n＝1，最大迭代次数为N，N为正整数，初始点为(t ,θ)，误差为ε；

各用户的功率分配方法包括以下步骤：步骤七、执行第n次迭代，用GUROBI或者MOSEK求解下式：* *

得到(t ,θ)，

式中：

t表示信干噪比最小的用户的信干噪比；

表示第m个AP和用户uq,k之间的功率控制系数；

为用户分簇结束后用户簇nq中的用户集合，重新定义用户uq,k是隶属于 中的第k个用户；

表示AP的集合，

θ是凸优化过程中引入的松弛变量；

n n n n

F(t,θ；t ,θ)是f(t,θ)在点(t ,θ)处的一阶泰勒展开n n

(t ,θ)表示第n次执行步骤七后得到的最优解；

* n

步骤八、判断|t ‑t |≤ε，或者n＝N，如果判断结果为是，则停止迭代，完成一次Cell‑free Massive MIMO下最大用户信干噪比的用户分簇及功率分配；如果判断结果为否，则执行步骤九；

n

t表示第n次迭代后得到的t的最优解，表示本次求解步骤七得到的最优解；

ε表示允许的最大误差；

n n * *

步骤九、更新：n＝n+1，(t ,θ)＝(t ,θ)，返回执行步骤七。

2.根据权利要求1所述的一种CF‑mMIMO下的用户分簇及功率分配方法，其特征在于用户ui的接收信号信干噪比SINR记为Γi，表示为：式中： 表示为用户集， 表示第i个用户， 表示第j个用户M表示AP的总数量，m表示第m个AP；

Γi表示用户ui的接收信号信干噪比；

表示第m个AP到用户ui的大尺度衰落系数；

表示第m个AP与用户ui的大尺度衰落系数的信道状态估计值；

L表示每个AP上的天线数量表示第m个AP和用户ui之间的功率控制系数表示具体的分簇结果， 表示用户ui属于用户簇nq，否则ρd表示每个AP的最大归一化发射功率即：由噪声功率N0归一化。

3.根据权利要求2所述的一种CF‑mMIMO下的用户分簇及功率分配方法，其特征在于用户ui的接收信号的数据速率表示为，表示用户簇集， 表示实际的用户簇数；

Bmax为系统的总带宽；

Γi表示用户ui的接收信号信干噪比；

Ri表示用户ui的接收信号的数据速率。

说明书

一种CF‑mMIMO下的用户分簇及功率分配方法

技术领域

[0001] 本申请涉及通信技术领域，具体涉及Cell‑free Massive MIMO下的功率控制技术

背景技术

[0002] 随着通信技术的发展，无线接入设备的数量和数据流量的消耗将呈现爆炸式增长。 Massive MIMO是一种很有前途的5G无线接入技术，具有多个天线的基站可以在同一时
频资源中同时为许多用户服务，该技术可以通过简单的信号处理来提供高吞吐量，可靠性
和能效。近年来，Cell‑Free Massive MIMO已作为一种新型的分布式大规模MIMO网络体系
结构被提出，该体系结构配备有大量的分布式接入点(AP)，服务于数量更少的分布式用户。
对于CF‑mMIMO，所有接入点通过前传链路连接到中央处理单元，并且不再存在小区边界的
概念。与现有的集中式Massive MIMO相比，CF‑mMIMO不仅在一定程度上继承了信道硬化的
良好特性，而且具有更高的能量效率和系统部署灵活性。分布式网络拓扑使CF‑mMIMO能够
从宏观多样性中受益，减轻集中式Massive MIMO的边缘影响，并确保更好的系统覆盖率，用
户服务公平性和网络吞吐量。显然，CF‑mMIMO网络有望成为未来的潜在架构之一。

[0003] 在CF‑mMIMO中，上行链路和下行链路以TDD模式工作，并且相干间隔通常分为三个阶段：上行链路训练，下行链路有效载荷数据传输和上行链路有效载荷数据传输。上行链路
训练完成信道的估计，对于TDD，在上行链路和下行链路中的信道增益是相同的，在下行链
路有效载荷数据传输和上行链路有效载荷数据传输的过程中会应用训练过程中得到的信
道估计信息进行数据信号的处理。

[0004] 目前的用户分簇及功率分配方法如：文献[1](Rezaei,A.R.Heidarpour,C.等在2020 年4月在IEEE Communications Letters,公开的《Tellambura and A.Tadaion,
“Underlaid Spectrum Sharing for Cell‑Free Massive MIMO‑NOMA》)的分簇方法存在的
计算量较大导致的不易于实现；功率分配方法还需要提供提高服务的公平性。

发明内容

[0005] 本发明是为了解决目前的用户分簇方法的计算量较大导致的不易于实现，以及用户的功率分配方法还需要提供提高服务的公平性，从而提供一种Cell‑free Massive MIMO
下最大用户信干燥比的用户分簇及功率分配方法。

[0006] Cell‑free Massive MIMO下最大用户信干燥比的用户分簇及功率分配方法，在 Cell‑free Massive MIMO下：所有AP在相同的时频资源下同时向所有用户提供服务；其特
征是：Cell‑free Massive MIMO下最大用户信干燥比的用户分簇及功率分配方法的具体方
法为：

[0007] 步骤一、初始化：设置初始分簇数I0、最大分簇数Imax；

[0008] 步骤二、执行第l+1次迭代：对所有用户进行分簇，分成I0+l个用户簇，对每个用户簇均通过SCA算法计算功率控制系数，得到最小的用户信干噪比Γmin；

[0009] 步骤三、判断步骤二获得的最小的用户信干噪比Γmin用户是否大于或等于预设的停止迭代条件Γmin≥γth，如果判断结果为是，则执行步骤四；如果判断结果为否且 I0+l＝
Imax，则执行步骤五；

[0010] γth表示信干噪比阈值，当用户接受信号的信干噪比大于这个阈值，该信号就可以被识别。

[0011] 步骤四、令l＝l+1；返回执行步骤二；

[0012] 步骤五、用户检测及排除步骤；

[0013] 令Ei表示用户ui的分簇参考变量，它通过如下公式获得：

[0014]

[0015] 式中：m和m′的物理含义相同；

[0016] ui表示第i个用户

[0017] M表示AP的总数量，m表示第m个AP

[0018] βm,ui表示第m个AP到用户ui的大尺度衰落系数

[0019] γm,ui表示第m个AP与用户ui的大尺度衰落系数的信道状态估计值；

[0020] L表示每个AP上的天线数量

[0021] ρd表示每个AP的最大归一化发射功率即：由噪声功率N0归一化；

[0022] 由于已知分簇数为Il＝I0+l，则每个簇的最小用户数Z1为不大于用户数比分簇数的最大取整数，通过解下面的凸优化问题得到每个用户的具体的分簇结果

[0023]

[0024]

[0025]

[0026] 式中：Ei表示用户ui的分簇参考变量

[0027] 表示具体的分簇结果， 表示用户ui属于用户簇nq，否则

[0028] 表示为用户集， 表示第i个用户。

[0029] 表示用户簇集，nq表示隶属于 的第q个簇，N为正整数；

[0030] Smax衡量指标，表示信干噪比最大的用户簇的近似信干噪比

[0031] Z1每个簇的最小用户数；

[0032] 各用户的功率分配方法包括以下步骤：

[0033] 步骤六、初始化：n＝1，最大迭代次数为N，N为正整数，初始点为(t1,θ1)，误差为ε。；

[0034] 步骤七、执行第n次迭代，用GUROBI或者MOSEK求解下式：

[0035]

[0036]

[0037]

[0038]

[0039]

[0040] ，得到(t*,θ*)，

[0041] 式中：

[0042] t表示信干噪比最小的用户的信干噪比

[0043] 表示第m个AP和用户uq,k之间的功率控制系数

[0044] 为用户分簇结束后用户簇nq中的用户集合，重新定义用户uq,k是隶属于 中的第 k个用户

[0045] 表示AP的集合，

[0046] 是凸优化过程中引入的松弛变量

[0047] F(t,θ；tn,θn)是f(t,θ)在点(tn,θn)处的一阶泰勒展开

[0048] (tn,θn)表示第n次执行步骤七后得到的最优解。

[0049] 步骤八、判断|t*‑tn|＜ε，或者n＝N，如果判断结果为是，则停止迭代，完成一次Cell‑free Massive MIMO下最大用户信干燥比的用户分簇及功率分配；如果判断结果为
是，否，则执行步骤九；

[0050] tn表示第n次迭代后得到的t的最优解。

[0051] t*表示本次求解步骤七得到的最优解

[0052] ε表示允许的最大误差

[0053] 步骤九、更新：n＝n+1，(tn,θn)＝(t*,θ*)，返回执行步骤七。

[0054] 通过本方法可以最大化所有用户的最小SINR，从而最大化接入网络的用户，最小化用户的掉线率。本申请的用户分簇的方法，该方法只需要知道用户的大尺度衰落信息，而
不需要知道用户的具体位置信息就可以对用户进行分簇，该分簇方法的计算量较低，更易
于实现；第二个是基于SCA的功率控制系数分配算法，该方法保证了最大化簇内所有用户的
最小SINR，保证了用户服务的公平性以及用户的接入。

附图说明

[0055] 图1是应用本申请所提出的分簇方法与其他方法的仿真对比示意图；其中：位于最上方的曲线是无分簇的仿真曲线图；位于最下方的曲线是本的仿真曲线图；位于中间的曲
线是对比文件1的仿真曲线图；位于最下方的曲线是本的仿真曲线图；

[0056] 图2是应用本申请所提出的分簇方法与其他方法的时间与通过率的仿真对比示意图；

[0057] 图3是应用本申请所提出的分簇方法与其他方法的流量与通过率的仿真对比示意图；

具体实施方式

[0058] 具体实施方式一、1、Cell‑free Massive MIMO下最大用户信干燥比的用户分簇及功率分配方法，在Cell‑free Massive MIMO下：所有AP在相同的时频资源下同时向所有用
户提供服务；其特征是：Cell‑free Massive MIMO下最大用户信干燥比的用户分簇及功率
分配方法的具体方法为：

[0059] 步骤一、初始化：设置初始分簇数I0、最大分簇数Imax；

[0060] 步骤二、执行第l+1次迭代：对所有用户进行分簇，分成I0+l个用户簇，对每个用户簇均通过SCA算法计算功率控制系数，得到最小的用户信干噪比Γmin；

[0061] 步骤三、判断步骤二获得的最小的用户信干噪比Γmin用户是否大于或等于预设的停止迭代条件Γmin≥γth，如果判断结果为是，则执行步骤四；如果判断结果为否且 I0+l＝
Imax，则执行步骤五；

[0062] γth表示信干噪比阈值，当用户接受信号的信干噪比大于这个阈值，该信号就可以被识别。

[0063] 步骤四、令l＝l+1；返回执行步骤二；

[0064] 步骤五、用户检测及排除步骤；

[0065] 令Ei表示用户ui的分簇参考变量，它通过如下公式获得：

[0066]

[0067] 式中：m和m′的物理含义相同；

[0068] ui表示第i个用户

[0069] M表示AP的总数量，m表示第m个AP

[0070] 表示第m个AP到用户ui的大尺度衰落系数

[0071] 表示第m个AP与用户ui的大尺度衰落系数的信道状态估计值；

[0072] L表示每个AP上的天线数量

[0073] ρd表示每个AP的最大归一化发射功率即：由噪声功率N0归一化；

[0074] 由于已知分簇数为Il＝I0+l，则每个簇的最小用户数Z1为不大于用户数比分簇数的最大取整数，通过解下面的凸优化问题得到每个用户的具体的分簇结果

[0075]

[0076]

[0077]

[0078]

[0079] 式中：Ei表示用户ui的分簇参考变量

[0080] 表示具体的分簇结果， 表示用户ui属于用户簇nq，否则

[0081] 表示为用户集， 表示第i个用户。

[0082] 表示用户簇集，nq表示隶属于 的第q个簇，N为正整数；

[0083] Smax衡量指标，表示信干噪比最大的用户簇的近似信干噪比

[0084] Z1每个簇的最小用户数；

[0085] 各用户的功率分配方法包括以下步骤：

[0086] 步骤六、初始化：n＝1，最大迭代次数为N，N为正整数，初始点为(t1,θ1)，误差为ε。；

[0087] 步骤七、执行第n次迭代，用GUROBI或者MOSEK求解下式：

[0088]

[0089]

[0090]

[0091]

[0092]

[0093] ，得到(t*,θ*)，

[0094] 式中：

[0095] t表示信干噪比最小的用户的信干噪比

[0096] 表示第m个AP和用户uq,k之间的功率控制系数

[0097] 为用户分簇结束后用户簇nq中的用户集合，重新定义用户uq,k是隶属于 中的第 k个用户

[0098] 表示AP的集合，

[0099] 是凸优化过程中引入的松弛变量

[0100] F(t,θ；tn,θn)是f(t,θ)在点(tn,θn)处的一阶泰勒展开

[0101] (tn,θn)表示第n次执行步骤七后得到的最优解。

[0102] 步骤八、判断|t*‑tn|＜ε，或者n＝N，如果判断结果为是，则停止迭代，完成一次Cell‑free Massive MIMO下最大用户信干燥比的用户分簇及功率分配；如果判断结果为
是，否，则执行步骤九；

[0103] tn表示第n次迭代后得到的t的最优解。

[0104] t*表示本次求解步骤七得到的最优解

[0105] ε表示允许的最大误差

[0106] 步骤九、更新：n＝n+1，(tn,θn)＝(t*,θ*)，返回执行步骤七。

[0107] 用户ui的接收信号信干噪比SINR记为Γi，表示为：

[0108]

[0109] 表示为用户集， 表示第i个用户， 示第j个用户；

[0110] M表示AP的总数量，m表示第m个AP；

[0111] Γi表示用户ui的接收信号信干噪比；

[0112] 表示第m个AP到用户ui的大尺度衰落系数；

[0113] 表示第m个AP与用户ui的大尺度衰落系数的信道状态估计值；

[0114] L表示每个AP上的天线数量；

[0115] 表示第m个AP和用户ui之间的功率控制系数；

[0116] 表示具体的分簇结果， 表示用户ui属于用户簇nq，否则

[0117] ρd表示每个AP的最大归一化发射功率即：由噪声功率N0归一化。

[0118] 用户ui的接收信号的数据速率表示为，

[0119]

[0120] 表示用户簇集， 表示实际的用户簇数。

[0121] Bmax为系统的总带宽；

[0122] Gi表示用户ui的接收信号信干噪比；

[0123] Ri表示用户ui的接收信号的数据速率。

[0124] 原理：SCA方法

[0125] SCA方法描述了一种解决非凸优化问题的通用方案，其中在每次迭代中，非凸可行集由一个内凸近似值近似。后者是使用非凸约束函数的上限定义的。在适当的条件下，建立
到KKT点的单调收敛。

[0126] 考虑以下通用优化问题：

[0127] minf(x)

[0128] (P)s.t.gi(x)≤0,i＝1,…,m

[0129]

[0130] 其中f，gi(i＝1,…,m)都是 上的连续可微函数。另外，假设函数f和最后m‑p(对于p＜m)约束函数gp+1,…,gm在 上是凸的。因此，问题的“非凸部分”是由于前p个约束函数
g1,…,gp的非凸性引起的。假设gi(i＝1,…,p)存在一个凸的上估计函数，假设存在
以及连续函数Gi： 使得对于所有的 y∈Y都有 gi(x)≤Gi(x,y)，对一个固
定的y，Gi(·,y)是凸函数且连续可微。向量y充当参数向量的角色，并且相应地Y称为允许
参数集。SCA方法的基本思想是，在每次迭代时，将每个非凸函数gi(i＝1,…,p)替换为选择
适当的参数向量y的上凸近似函数 因此，在该方法的步骤k(k≥1)，需要解决
以下形式的凸问题：

[0131] minf(x)

[0132] (Pk)s.t.Gi(x,yk)≤0,i＝1,…,p

[0133] gj(x)≤0,j＝p+1,…,m

[0134]

[0135] 向量yk是一个固定的参数向量，取决于问题的解(Pk‑1)。

[0136] 在CF‑mMIMO环境中，所有AP在相同的时频资源下同时向所有用户提供服务。但是，如果系统中的用户过多，则会发生导频污染，并且用户之间的干扰可能太高而无法满足信
噪比(SINR)阈值要求。因此，本专利将介绍用户分簇的概念，即每个分片中的用户分为多个
用户簇，并且同一簇中的用户共享相同的时频资源。 表示为用户集， 表示第i个用
户。 表示用户簇集，nq表示隶属于 的第q个簇。 表示实际的用户簇数。 是一个二
进制指标， 表示用户ui属于用户簇nq，否则 用户分簇结束后，可以定义
为用户簇nq中的用户集合，重新定义用户uq,k是隶属于 中的第k个用户。

[0137] 表示AP集，每个AP具有最大功率Pmax以及L个天线。Bmax为系统的总带宽，分配给每个簇的带宽为 分配给每个簇的功率为 则用户ui的接收信号信干
噪比(Singal to Interference plus Noise Ratio,SINR)，记为Γi，可表示为

[0138]

[0139] 用户ui的接收信号的数据速率可以表示为，

[0140]

[0141] 对于用户ui而言，接入网络有两个条件，第一个条件是他的信干噪比不能低于信干噪比阈值γth，否则不容易在同一时频资源的用户中将其辨别出来；第二个条件是他的用
户数据速率不能低于他的最小需求速率，在本申请中

[0142] 将用户的最小需求速率同一标记为Rmin。

[0143] 1.用户分簇的方法：

[0144] Ei表示用户ui的分簇参考变量，他可以通过如下公式计算：

[0145]

[0146] 已知分簇数为Il＝I0+l，每个簇的最小用户数Z1为不大于用户数比分簇数的最大取整数，通过解下面的凸优化问题就可以得到具体的分簇方案

[0147]

[0148]

[0149]

[0150]

[0151] 2.SCA算法计算功率控制系数

[0152] 这个部分对所有用户分好簇后，对一个簇内的所有用户执行功率控制系数算法。以簇nq为例说明，其他所有簇同理。该算法的目的是使该用户簇中用户的最小SINR最大化，
从而使超过阈值的用户数量最大化，使得尽可能多的用户进入网络，那么这个问题就可以
被构建为，

[0153]

[0154]

[0155]

[0156] 这个问题明显是非凸的问题，无法通过计算器求解，因此引入松弛变量并且应用SCA 来找到一个近似解，这个问题就会被转换为，

[0157]

[0158]

[0159]

[0160]

[0161]n n

[0162] 其中 F(t,θ；t ,θ)n n
是 f(t,θ)在点(t ,θ)处的一阶泰勒展开。

[0163] 3.初始分簇数和最大分簇数

[0164] 初始分簇时，需要考虑避免簇内用户过多发生导频污染，因此设置 其中为总用户数。设置最大分簇数的时候需要避免分簇过多，导致分配给每个簇的带宽资源
过低，从而导致用户数据速率过低而无法满足最低速率要求，因此最大分簇数需要满足
得到 Imax选择满足这个条件的最大正
整数。

[0165] 五、发明效果：

[0166] 通过本方法可以最大化所有用户的最小SINR，从而最大化接入网络的用户，最小化用户的掉线率。本方法提出了两个子方法，第一个是用户分簇的方法，该方法只需要知道
用户的大尺度衰落信息，而不需要知道用户的具体位置信息就可以对用户进行分簇，该分
簇方法的计算量较低，更易于实现；第二个是基于SCA的功率控制系数分配算法，该算法保
证了最大化簇内所有用户的最小SINR，保证了用户服务的公平性以及用户的接入。

[0167] 本节中将给出算法的数值仿真结果，AP在1km×1km的正方形区域内均匀分布，所有用户随机分布。在这里使用路径损失模型如下：

[0168]

[0169] 其中 表示用户ui和基站m之间的距离，而 表示用户ui和基站m之间的路径损耗。

[0170] d0＝10m，d1＝50m， 表示阴影衰落，并且 具体参数见下表

[0171]符号 参数 参数值
σsh 阴影衰落方差 8dB
N0 噪声功率谱密度 ‑174dBm/Hz
Bmax 系统总带宽 20MHz
Pmax AP最大发射功率 200mW
Tc 相干间隔 200
τ 训练导频长度 40
Rmin 用户最小速率需求 600kbps
γth 信干噪比阈值 8dB

[0172] 图1显示了三种方法在不同访问用户下的掉线率。可以看出，三种方案的掉线率将随访问用户数量的增加而增加，我们提出的方案在不同访问用户数量下的掉线率最低。可
以看出，三种方案的掉线率将随访问用户数量的增加而增加，我们提出的方案在不同访问
用户数量下的掉线率最低。不分簇的方案出现中断的可能性最高，其原因很明显： CF‑
mMIMO体系结构要求天线的数量必须大于用户数量，一旦不满足此要求，掉线率将迅速增
加。而另外两种方案引用了用户分簇，因此掉线率较低。本申请提出的分簇方案具有最低的
掉线率，因为本申请的最小SINR大于其他两个方案中的最小SINR，因此可以方便更多用户，
可以满足SINR阈值。