CN114726685B

CN114726685B - 一种低复杂度的大规模mimo通信系统下行链路信道估计方法

Info

Publication number: CN114726685B
Application number: CN202210211539.6A
Authority: CN
Inventors: 王顺; 季晨; 戴继生
Original assignee: Jiangsu University
Current assignee: Jiangsu University
Priority date: 2022-03-04
Filing date: 2022-03-04
Publication date: 2024-04-09
Anticipated expiration: 2042-03-04
Also published as: CN114726685A

Abstract

本发明公开了一种低复杂度的大规模MIMO通信系统下行链路信道估计方法，步骤1：基站采用了一个具有M_t根天线的均匀线性阵列，下行链路中的移动用户装配M_r根天线的均匀线性阵列，在T个时刻内，基站发送导频信号矩阵S，则移动用户接收到的信号为Y＝SA_T(β)WB^T+N。步骤2：将接收信号进行向量化。步骤3：初始化参数。步骤4：设置迭代次数计数变量j＝1。步骤5：固定q(α)，q(γ),β，B,更新q(w)。步骤6：固定q(w),q(γ),β，B,更新q(α)。步骤7：固定q(w)，q(α)，β，B,更新q(γ)。步骤8：固定q(w)，q(α)，q(γ)，β,更新B。步骤9：固定q(w),q(α),q(γ),B,更新β。步骤10：判断迭代计数变量j是否达到上限J或γ是否收敛，如果都不满足，则迭代计数变量j＝j+1，并返回步骤5。步骤11：估计最终的信道H。本发明在信道估计性能上有明显的提升，同时复杂度较低。

Description

一种低复杂度的大规模MIMO通信系统下行链路信道估计方法

技术领域

本发明属于无线通信领域，涉及一种多输入多输出(Multi-input Multi-output，MIMO)通信系统的信道估计方法，具体来说是一种基于变分贝叶斯推理(variationalBayesian inference,VBI)的低复杂度大规模MIMO通信系统的下行链路信道估计方法。

背景技术

大规模多输入多输出(Multiple-Input Multiple-Output，MIMO)技术是下一代无线通信的研究热点和关键技术之一。将大规模MIMO技术应用于无线通信中，可获得较高的信道增益，同时其极低的系统功耗可以使得信号的覆盖范围得到有效扩大。相比于传统无线通信系统，大规模MIMO系统具有很高的频谱利用率、能量效率和系统的鲁棒性能。

信道信息的获取在通信系统中占有举足轻重的地位，然而在无线通信系统中信道信息都是不可知的，因此信道估计是无线通信中极为重要的技术，信道估计的精度将直接影响整个通信系统的性能。同时，随着天线规模的扩大，信道矩阵的维度也随之线性增长，高计算复杂度的问题愈发成为信道估计技术的制约因素。为了降低计算复杂度，现有方法在进行信道估计时大多将用户端的阵列信息并入路径增益矩阵中，进而利用复杂度较低的行稀疏特性进行信道估计。目前，人们已经提出了很多利用行稀疏特性的方法解决大规模MIMO信道估计问题，例如在文献J.Dai,A.Liu and V.K.N.Lau,FDD Massive MIMO ChannelEstimation with Arbitrary 2D-Array Geometry,IEEE Transactions on SignalProcessing,vol.66,no.10,pp.2584-2599,15May,2018中提出了一种基于离网稀疏贝叶斯学习的大规模MIMO通信系统的信道估计方法，由于此方法忽略了用户端的阵列反馈信息，因此会损失一定的信道估计性能。

发明内容

针对现有方法的不足，本发明将提出一种低复杂度的基于VBI的大规模MIMO通信系统的下行链路信道估计方法。

用于实现本发明的技术解决方案包括如下步骤：

步骤1：基站采用了一个具有M_t根天线的均匀线性阵列，下行链路中的移动用户装配M_r根天线的均匀线性阵列，在T个时刻内，基站发送导频信号矩阵S，则移动用户接收到的信号为Y＝SA_T(β)WB^T+N。

步骤2：将接收信号进行向量化。

步骤3：初始化参数。

步骤4：设置迭代次数计数变量j＝1。

步骤5：固定q(α),q(γ),β,B,更新q(w)。

步骤6：固定q(w),q(γ),β,B,更新q(α)。

步骤7：固定q(w),q(α),β,B,更新q(γ)。

步骤8：固定q(w),q(α),q(γ),β,更新B。

步骤9：固定q(w),q(α),q(γ),B,更新β。

步骤10：判断迭代计数变量j是否达到上限J或γ是否收敛，如果都不满足，则迭代计数变量j＝j+1，并返回步骤5。

步骤11：估计最终的信道。

本发明的有益效果为：

利用VBI方法，本发明获得了一种低复杂度基于变分贝叶斯推理的大规模MIMO通信系统的下行链路信道估计方法。与现有方法相比，本发明在信道估计性能上有明显的提升，同时本发明的计算复杂度为明显低于现有方法的/>

附图说明

图1是本发明实施流程图。

图2在信噪比为0dB，发送端天线M_t为150，用户天线M_r为8根，网格数设定为200时，导频时刻T由40到90变化时对信道进行500次估计，仿真结果如图2所示。

图3在导频时刻为50，发送端天线M_t为100，用户天线M_r为8根，网格数设定为200时，信噪比从-10dB到10dB变化时对信道进行500次估计，仿真结果如图3所示。

具体实施方式

下面结合附图对本发明作进一步说明。

本发明的实施如图1所示，包括如下：

(1)基站采用了一个具有M_t根天线的均匀线性阵列，下行链路中的移动用户数装配M_r根天线的均匀线性阵列，在T个时刻内，基站发送维度为T×M_t的导频信号矩阵S，移动用户接收到的信号为Y＝SA_T(β)WB^T+N，其中:

A_T(β)＝[a_T(θ₁+β₁),a_T(θ₂+β₂),...,a_T(θ_N+β_N)]表示发送端阵列流型矩阵，

表示导向矢量，n＝1,2,...,N，

(·)^T表示转置，

λ表示载波波长，d为相邻天线阵元之间的间距，

表示对角域/>均匀划分成N个网格点，即/>

中的元素β_i表示θ_i上的角度偏差，

W代表N×M_r维的稀疏矩阵，

B代表M_r×M_r维的未知矩阵，

N代表T×M_r维均值为0，精度为α的高斯白噪声矩阵。

(2)将接收信号重新排列，可得y＝Φ_Bw+n，其中：

y＝vec(Y)，Φ(β)＝SA_T(β)，w＝vec(W)，n＝vec(N)，

vec(·)代表向量化，

代表克罗内科积，

(3)初始化参数：μ_m＝0_N，γ＝0_N，a＝b＝10^-10，μ_m表示w_m的均值，B^T为Y的右奇异向量，其中：

0_N表示维度为N×1的0向量。

(4)设置迭代次数计数变量j＝1。

(5)固定q(α),q(γ),β,B,更新q(w)：

其中：

q(α),q(γ)，q(w)分别表示α，γ，w的近似后验概率密度函数，

表示均值为μ，方差为Σ的复高斯分布，

w_m代表W的第m列向量，

b_m表示B的第m列向量，

(·)^H表示共轭转置，

diag(·)表示对角运算矩阵。

(6)固定q(w),q(γ),β,B,更新q(α)：

q(α)＝Γ(α；a_α,b_α)，

其中：

Γ(·|a,b)表示形状参数为a，速率参数为b的伽马分布，

a_α＝TM_r+a，

||·||₂表示矩阵的2范数，tr(·)表示矩阵的迹。

(7)固定q(w),q(α),β,B,更新q(γ)：

其中：

γ表示w的精度向量，

γ_n表示γ的第n个元素，

[·]_n,n表示矩阵的第n个对角线元素。

(8)固定q(w),q(α),q(γ),β,更新B：

其中：

g_m＝tr(Φ(β)Σ_mΦ^H(β))，m＝1,2,...,M_r。

(9)固定q(w),q(α),q(γ),B,更新β：

其中：

sign(·)表示取正负号运算，

ζ＝[ζ(β₁),ζ(β₂),…ζ(β_N)]^T，

ζ(β_n)＝2Re(a′(θ_n+β_n)^HS^HSa(θ_n+β_n)c₁+a′(θ_n+β_n)^HS^Hc₂)，n＝1,2,...,N，

Re(·)表示取实部运算，

表示U的第n行，/>表示Σ_m的第(i,i)个元素，

(·)^*代表共轭运算，

a′(θ_n+β_n)表示a(θ_n+β_n)在θ_n+β_n处的导数。

(10)判断迭代计数变量j是否达到上限J＝100或γ是否收敛，如都不满足，则迭代计数变量j＝j+1，并返回(5)。

(11)估计最终的信道H：H＝Φ(β)UB^T。

下面通过列举仿真实验结果对本发明的效果做进一步阐述。

为了评估本方法的性能，假设基站装配了一个具有M_t＝100根天线的均匀线性阵列，下行链路的工作频率为2170MHz，无线信道由3GPP spatial channel model(SCM)模型随机产生，基站发送导频信号矩阵S，S的每个元素均为零均值单位方差的独立高斯分布，加性噪声假设为高斯白噪声。

实验条件

采用本发明在信噪比为0dB，发送端天线M_t为150，用户天线M_r为8根，网格数设定为200时，导频时刻T由40到90变化时对信道进行500次估计，仿真结果如图2所示。

采用本发明在在导频时刻为50，发送端天线M_t为100，用户天线M_r为8根，网格数设定为200时，信噪比从-10dB到10dB变化时对信道进行500次估计，仿真结果如图3所示。

实验分析

从图2和图3可以看出，本发明可以精确地估计出大规模MIMO通信系统的下行链路信道信息，其NMSE性能明显优于现有方法。

上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技术所创的等效方式或变更均应包含在本发明的保护范围之内。

Claims

1.一种低复杂度的大规模MIMO通信系统下行链路信道估计方法，其特征在于，包括如下步骤：

步骤1：基站采用了一个具有M_t根天线的均匀线性阵列，下行链路中的移动用户装配M_r根天线的均匀线性阵列，在T个时刻内，基站发送维度为T×M_t的导频信号矩阵S，移动用户接收到的信号为Y＝SA_T(β)WB^T+N；

步骤2：将接收信号进行向量化；

步骤3：固定q(α),q(γ),β,B,更新q(w)；

步骤4：固定q(w),q(γ),β,B,更新q(α)；

步骤5：固定q(w),q(α),β,B,更新q(γ)；

步骤6：固定q(w),q(α),q(γ),β,更新B；

步骤7：固定q(w),q(α),q(γ),B,更新β；

步骤8：判断迭代计数变量j是否达到上限J或γ是否收敛，如果都不满足，则迭代计数变量j＝j+1，并返回步骤3；

步骤9：估计最终的信道；

所述步骤1中：

表示导向矢量，

(·)^T表示转置，

λ表示载波波长，d为相邻天线阵元之间的间距，

表示对角域/>均匀划分成N个网格点，即/>

中的元素β_i表示θ_i上的角度偏差，

W代表N×M_r维的稀疏矩阵，

B代表M_r×M_r维的未知矩阵，

N代表T×M_r维均值为0，精度为α的高斯白噪声矩阵；

所述步骤2中对接收信号进行向量化的方法：

y＝Φ_Bw+n，

其中：y＝vec(Y)，Φ(β)＝SA_T(β)，w＝vec(W)，n＝vec(N)，vec(·)代表向量化，

代表克罗内科积；

所述步骤3中更新q(w)的方法：

其中：

表示均值为μ，方差为Σ的复高斯分布，

w_m代表W的第m列向量，

b_m表示B的第m列向量，

(·)^H表示共轭转置，

diag(·)表示对角运算矩阵；

所述步骤4中更新q(α)的方法：

q(α)＝Γ(α；a_α,b_α)，

其中：

Γ(·|a,b)表示形状参数为a，速率参数为b的伽马分布，

a_α＝TM_r+a，

||·||₂表示矩阵的2范数，tr(·)表示矩阵的迹；

所述步骤5中更新q(γ)的方法：

其中：

γ表示w的精度向量，

γ_n表示γ的第n个元素，

[·]_n,n表示矩阵的第n个对角线元素；所述步骤6中更新B的方法：

其中：

U＝[μ₁,μ₂,...,μ_Mr]，

g_m＝tr(Φ(β)Σ_mΦ^H(β))；

所述步骤7中更新β的方法：

其中：

sign(·)表示取正负号运算，

δ＝[ζ(β₁),ζ(β₂),…ζ(β_N)]^T，

Re(·)表示取实部运算，

表示U的第n行，/>表示Σ_m的第(i,i)个元素，

(·)^*代表共轭运算，

a'(θ_n+β_n)表示a(θ_n+β_n)在θ_n+β_n处的导数；

所述步骤9中估计最终信道的方法：

H＝Φ(β)UB^T；

还包括：初始化参数μ_m＝0_N，γ＝0_N，a＝b＝10^-10，以及设置迭代次数计数变量初始值为j＝1。