一种在线鲁棒的分布式多小区大规模MIMO预编码方法

一种在线鲁棒的分布式多小区大规模mimo预编码方法
技术领域
1.本发明涉及一种分布式多小区大规模mimo预编码方法。


背景技术：

2.在大规模mimo系统中，基站端布置大规模天线阵列以同时服务多个用户。采用大规模mimo技术可以有效降低用户间干扰，大幅提高无线通信系统的能效和谱效。波束域传输是指基站侧通过统一的酉变换将发射信号转换到波束域，充分利用大规模天线阵列的空间角度分辨率和用户信道在波束域中的局部特性。
3.在多小区协同传输中，几个邻近小区的基站交互信息，共同协作进行预编码设计。相比于单小区传输，能够显著的减少共信道干扰，提升系统的传输速率。若在协作过程中，基站之间交互所有的信道状态信息和传输信号，会造成较大的反馈开销。相反的，若基站之间只进行有限的信息交换，每个小区相对独立的进行本小区的预编码设计，则可以节省反馈开销并降低计算复杂度。这种传输方法称为分布式多小区协同传输。
4.在实际传输过程中，由于用户的移动性，信道状态信息会随着时间而变化。当信道状态信息发生较大的改变时，发射端需要对传输策略进行相应的调整。在线优化方法是一种顺序决策过程，可用于这种动态传输设计。具体而言，在每个信道状态信息变化的阶段，发射端根据过去所有的信道状态信息选择这个阶段的传输策略。此后，发射端收到相应的反馈，并继续在下一阶段决定新的传输策略。为了权衡通信各小区的系统容量和能源消耗，提出以能效(各小区用户速率与总功率消耗的比值)为优化目标的问题。传统的能效优化问题常基于瞬时信道状态信息，而在大规模mimo系统中，瞬时信道状态信息的获取有难度。相较而言，统计信道状态信息更容易获取。在对统计信道状态信息进行估计时，由于导频污染等原因，估计值会有误差。因此，传输策略对于信道估计误差的鲁棒性是重要的研究方向。


技术实现要素：

5.发明目的：针对上述现有技术，提出一种低复杂度的利用统计信道状态信息的在线鲁棒的分布式多小区大规模mimo预编码方法。
6.技术方案：一种在线鲁棒的分布式多小区大规模mimo预编码方法，包括如下步骤：
7.1)构建小区簇并配置天线阵列，通过波束赋形生成波束集合：
8.将若干个小区构建为一个小区簇，且各小区的基站配置大规模天线阵列，通过波束赋形生成能够覆盖全部小区的大规模波束集合；
9.2)各小区基站动态获取本小区内用户的统计信道状态信息：
10.所述统计信道状态信息的获取途径包括用户反馈、上行探测信号，或基站直接估计；当用户在通信过程中发生移动时，基站与用户之间的统计信道状态信息发生变化，各基站动态获取统计信道状态信息；
11.3)各小区基站根据动态获取的统计信道状态信息进行在线鲁棒的分布式多小区大规模mimo预编码设计：
12.在每一个统计信道状态信息变化的阶段，各基站利用确定性等同方法、charnes
‑
cooper变换和在线梯度上升算法，根据这一阶段及之前获取的所有信道状态信息预测下一阶段的功率分配矩阵。
13.进一步的，所述基站通过相同的酉变换生成能够覆盖整个小区的大规模波束，实现对空间资源的波束域划分，且基站在生成的波束域与该小区的用户进行通信。
14.进一步的，所述统计信道状态信息为波束域特征模式能量耦合矩阵，小区的用户在上行信道探测阶段发送上行探测信号至基站，基站根据接收到的探测信号估计波束域统计信道状态信息，第ν个小区的基站到第u个小区的用户的统计信道状态信息为：
[0015][0016]
其中，g
u,v
为第ν个小区的基站到第u个小区的用户的波束域信道矩阵，表示期望运算，
⊙
表示矩阵的阿达玛积，为矩阵g
u,v
的共轭矩阵。
[0017]
进一步的，所述在线鲁棒的分布式多小区大规模mimo预编码方法为分布式算法，各基站进行功率分配时，暂时将来自其他小区的干扰看作常数；功率分配完成，基站之间交互信息时，更新干扰信号；干扰加噪声协方差矩阵为将其与波束域信道矩阵合并，生成等效信道矩阵，表示为且等效统计信道状态信息表示为其中σ2为噪声方差，i
n
为n
×
n的单位矩阵，λ
j
为第j个小区的基站的功率分配矩阵，表示对矩阵g
u,j
进行共轭转置。
[0018]
进一步的，当用户在通信过程中发生移动时，所述统计信道状态信息发生变化，基站根据变化的统计信道状态信息实行在线鲁棒的分布式多小区大规模mimo预编码方法，其优化目标为各小区的能效值，优化变量为各个基站发送信号的协方差矩阵，约束条件为各基站发送信号的协方差矩阵满足功率约束；设第s个阶段的等效信道为则第u个小区的能效最大化在线鲁棒大规模mimo预编码问题表示为：
[0019][0020]
s.t.tr{λ
u
}≤p
max,u
,
[0021]
λ
u
≥0,λ
u diagonal,
[0022]
且，
[0023][0024]
p
u
(λ
u
)＝ξ
u tr{λ
u
}+mp
c,u
+p
s,u
[0025]
其中，m为小区的基站天线数，n为每个用户的天线数，为小区u的能效，为用户u的遍历速率，p
u
(λ
u
)为第u个小区的功耗，p
max,u
为基站的发送功率约束，ξ
u
为基站放大系数，p
c,u
为每根天线上耗散的动态功率，p
s,u
为在硬件中耗散的静态电路功率，log{.}表示对数运算，tr{
·
}表示取矩阵迹的运算，det{.}表示取矩阵的行列式运算。
[0026]
进一步的，所述步骤3)具体包括：
[0027]
3a)利用确定性等同方法求期望运算：
[0028]
基于大维随机矩阵理论，通过用户的等效波束域统计信道状态信息，迭代计算目标函数中用户速率的确定性等同辅助变量直至收敛；
[0029]
基于迭代计算得到的确定性等同辅助变量计算目标函数中用户速率的确定性等同表达；
[0030]
将用户速率的确定性等同表达代入能效最大化的在线鲁棒大规模mimo预编码问题中；
[0031]
3b)charnes
‑
cooper变换：
[0032]
对能效最大化的在线鲁棒大规模mimo预编码问题进行变量代换，变换后的问题是与原问题等效的凸优化问题；
[0033]
利用在线梯度上升算法求解变换后的等效问题，得到此问题的解，将此解代入到变量代换式中，计算出对应的功率分配矩阵；
[0034]
3c)在线梯度上升算法：
[0035]
经charnes
‑
cooper变换后的等效问题是凸问题，其目标函数为凸函数；求出目标函数对变量的梯度函数；
[0036]
将当前阶段变量的值加上设定步长和梯度函数的乘积，并将所得结果投影到问题可行域中，得到下一个阶段的变量的解，作为所述等效问题的解。
[0037]
有益效果：(1)基站与小区内的用户在波束域通信，能够与大规模mimo无线信道的空间特性相匹配，从而获取使用大规模天线阵列所带来的能效和谱效的提升；
[0038]
(2)各基站通过稀疏的探测信号获得各个用户的波束域统计信道状态信息，该方法适用于时分双工和频分双工系统，充分利用了本小区用户和其他小区用户的波束域统计信道状态信息；
[0039]
(3)利用确定性等同方法、charnes
‑
cooper变换和在线梯度上升算法进行在线鲁棒的分布式多小区大规模mimo波束域功率分配，能够显著降低优化问题求解和物理层实现的复杂度，加快了运算速度，并且对信道估计误差有一定鲁棒性。
附图说明
[0040]
图1为本发明的整体流程示意图；
[0041]
图2为本发明中基于确定性等同方法的算法流程图；
[0042]
图3为本发明中基于charnes
‑
cooper变换的算法流程图；
[0043]
图4为本发明中基于在线梯度上升算法的算法流程图。
具体实施方式
[0044]
下面结合附图对本发明做更进一步的解释。
[0045]
如图1所示，一种在线鲁棒的分布式多小区大规模mimo预编码方法，包括如下步骤：
[0046]
1)构建小区簇并配置天线阵列，通过波束赋形生成波束集合：
[0047]
将若干个小区构建为一个小区簇，且小区簇内各小区的基站配置大规模天线阵
列，其中天线数为102或103数量级大小，并通过波束赋形生成能够覆盖全部小区的大规模波束集合；
[0048]
其中，基站通过相同的酉变换生成能够覆盖整个小区的大规模波束，每个波束对空间资源进行划分，且基站在生成的波束域与该小区的用户进行通信。
[0049]
根据多小区大规模mimo下行链路通信场景，各基站侧配置m根发送天线的大规模天线阵列，m为102或103数量级大小。若共有u个小区，考虑每个小区有一个用户，每个用户配置n根接收天线。
[0050]
2)各小区基站动态获取本小区内用户的统计信道状态信息：
[0051]
统计信道状态信息的获取途径包括用户反馈、上行探测信号，或基站直接估计。当用户在通信过程中发生移动时，基站与用户之间的统计信道状态信息发生变化，各基站动态获取统计信道状态信息；
[0052]
3)各小区基站根据动态获取的统计信道状态信息进行在线鲁棒的分布式多小区大规模mimo预编码设计：
[0053]
在每一个统计信道状态信息变化的阶段，各基站利用确定性等同方法、charnes
‑
cooper变换和在线梯度上升算法，根据这一阶段及之前获取的所有信道状态信息预测下一阶段的功率分配矩阵。
[0054]
其中，统计信道状态信息为波束域特征模式能量耦合矩阵，其获取包括以下步骤：
[0055]
各个小区的用户在上行信道探测阶段发送上行探测信号至基站；
[0056]
基站根据接收到的探测信号估计波束域统计信道状态信息，第ν个小区的基站到第u个小区的用户的统计信道状态信息为：
[0057][0058]
其中，g
u,v
为第ν个小区的基站到第u个小区的用户的波束域信道矩阵，表示期望运算，
⊙
表示矩阵的阿达玛积，为矩阵g
u,v
的共轭矩阵。
[0059]
各基站将发送给各用户的空间域信号通过统一的酉变换变换到波束域，向本小区的每个用户发送信号。小区间的协作是分布式的。特别的，各基站进行功率分配时，暂时将来自其他小区的干扰看作常数。功率分配完成，基站之间交互信息时，更新干扰信号。
[0060]
将第u个小区的用户定义为用户u，假设基站发送给用户u的信号为x
u
，功率分配矩阵为λ
u
。特别的，干扰加噪声协方差矩阵为将其与波束域信道矩阵合并，生成等效信道矩阵，表示为且等效统计信道状态信息表示为：其中σ2为噪声方差，i
n
为n
×
n的单位矩阵，表示对矩阵g
u,j
进行共轭转置。用户u的遍历速率可以表示为：
[0061][0062]
考虑到在多小区大规模mimo通信中的能效问题，为了优化能效，因此需要对发送信号的功率分配矩阵进行优化，即在每个基站侧对发射波束进行功率分配。具体的，对于第u个小区，其功耗模型为：
[0063]
p
u
(λ
u
)＝ξ
u tr{λ
u
}+mp
c,u
+p
s,u
[0064]
其中，tr{
·
}表示取矩阵迹的运算，log{.}表示对数运算，det{.}表示取矩阵的行列式运算，ξ
u
(＞1)为基站放大系数，p
c,u
为每根天线上耗散的动态功率，p
s,u
为在硬件中耗散的静态电路功率。
[0065]
当用户在通信过程中发生移动时，信道状态信息发生变化，基站根据变化的信道状态信息实行在线鲁棒的分布式多小区大规模mimo预编码方法，其优化目标为各小区的能效值，优化变量为各个基站发送信号的协方差矩阵，约束条件为各基站发送信号的协方差矩阵满足功率约束。设第s个阶段的等效信道为则第u个小区的能效最大化在线鲁棒大规模mimo预编码问题表示为：
[0066][0067]
s.t.tr{λ
u
}≤p
max,u
,
[0068]
λ
u
≥0,λ
u diagonal,
[0069]
且，
[0070][0071]
其中，m为小区的基站天线数，n为每个用户的天线数，为小区u的能效，p
max,u
为基站的发送功率约束。
[0072]
由于此能效优化问题的目标函数非凸，很难得到全局最优解，实现复杂度很高。因此步骤3)中通过利用波束域统计信道状态信息的在线多小区分布式大规模mimo波束域功率分配问题，包括了基于确定性等同方法、charnes
‑
cooper变换和在线梯度上升算法的优化算法。
[0073]
具体的，在计算能效表达式中的分子项，即小区所有用户速率和时，需要对信道进行遍历，计算期望值。由于该期望没有闭式表达式，因而需要monte
‑
carlo仿真计算。为了避免高复杂度的求期望运算，本发明利用大维矩阵随机理论计算目标函数的确定性等同表达，降低计算复杂度。确定性等同方法仅利用统计信道状态信息，通过迭代计算确定性等同辅助变量，即可获得目标函数值的逼近结果。同时，由于确定性等同的结果可以很好逼近速率项的准确表达。
[0074]
如图2所示为基于确定性等同方法的实现过程示意图，其实现的具体过程包括：
[0075]
步骤i：首先对于每一个用户引入确定性等同辅助变量，其中第u个用户的辅助变量为：
[0076][0077][0078]
其中，辅助变量和由下式确定：
[0079][0080]
ψ
u,n
＝[λ
u
(i
m
+ξ
u
λ
u
)
‑1]
n,n
,n＝1,...,n,
[0081]
其中u
u,m
是u
u,u
的第m列，u
u,u
是确定性酉矩阵，i
m
为m
×
m的单位矩阵；
[0082]
步骤ⅱ：固定辅助变量ψ
u
的值，根据步骤i给出的计算式求出ξ
u
，并由求得的ξ
u
计算出新的ψ
u
的值，迭代计算直到辅助变量的值收敛；
[0083]
步骤ⅲ：根据收敛后的辅助变量的值，计算出系统频谱效率的确定性等同值，表示为：
[0084][0085]
如图3所示为charnes
‑
cooper变换的实现过程示意图，该算法包括以下步骤：
[0086]
步骤i：首先对能效优化问题进行变量代换。引入一个辅助变量：
[0087][0088]
其中，p
b,u
＝mp
c,u
+p
s,u
。根据λ
u
的可行域，得到变量x
u
的可行域为：
[0089][0090]
步骤ⅱ：对优化问题的目标函数和约束条件进行变量代换，变换后的等效问题以x
u
为变量：
[0091][0092][0093]
其中目标函数为：
[0094][0095]
步骤ⅲ：利用在线梯度上升算法对变换后的等效问题进行求解。解出x
u
后，利用以下的变换式即可求得原问题的解λ
u
：
[0096][0097]
如图4所示为基于在线梯度上升算法的实现过程示意图，其实现的具体过程包括以下步骤：
[0098]
步骤i：首先计算问题目标函数对求解变量的导数：
[0099][0100]
其中，
[0101]
步骤ⅱ：设统计信道状态信息的变化的指标为s。在第s个阶段，基站由下式计算第s+1阶段的功率分配矩阵：
[0102][0103]
其中为步长，通常设定为π
u
(y)为基于对角矩阵的映射函数，表示
为：
[0104][0105][0106]
其中β＞0并且β的值可以用二分法计算。
[0107]
步骤ⅲ：当统计信道状态信息变化时，即由阶段s变到阶段s+1时，重复上述过程。
[0108]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。