大规模mimo系统中低复杂度信道估计方法

文档序号:10675475阅读:469来源:国知局
大规模mimo系统中低复杂度信道估计方法
【专利摘要】本发明提供的大规模MIMO系统中低复杂度信道估计方法,该方法包括:对基站接收的信号进行向量化,得到基于克罗内克模型的接收信号模型;对上述接收信号模型,根据MMSE估计原理,得到包含有对信号协方差矩阵和干扰协方差矩阵进行求逆运算的MMSE估计结果;对MMSE估计结果中的求逆运算转换为求解线性方程组的问题,通过共轭梯度法求解方程组的近似最优解,将近似最优解应用与MMSE的估算结果中,获取修正的MMSE估算结果,降低MMSE信道估计的复杂度。本发明通过将求逆运算转换为求解方程组的问题,然后利用共轭梯度法求得方程组的近似最优解,以达到降低计算复杂度的目的。
【专利说明】
大规模ΜI MO系统中低复杂度信道估计方法
技术领域
[0001] 本发明设及通信技术领域,尤其是大规模MIM0(Massive Multiple Input Multiple Output,Massive-MIMO)系统中低复杂度信道估计算法。
【背景技术】
[0002] 由于大规模ΜΙΜΟ系统在基站端配置有数百根天线,并且在上行链路进行信道估计 时需要对协方差矩阵进行求逆的操作,运样就导致的信道估计的计算复杂度为0(λ/ι'),μ 为协方差矩阵的维度,运个在硬件实现过程中将是一个极其复杂度过程。

【发明内容】

[0003] 本发明要解决的技术问题在于提供一种大规模ΜΙΜΟ系统中低复杂度信道估计方 法,W解决现有技术存在的问题。
[0004] 本发明采用W下技术方案:
[0005] 大规模ΜΙΜΟ系统中低复杂度信道估计方法,该方法包括W下步骤:
[0006] 步骤1:对基站接收的信号进行向量化,得到基于克罗内克模型的接收信号模型;
[0007] 步骤2:对上述接收信号模型,根据MMSE估计原理,得到包含有对信号协方差矩阵 和干扰协方差矩阵进行求逆运算的MMSE估计结果;
[000引步骤3:对MMSE估计结果中的求逆运算转换为求解线性方程组的问题,通过共辆梯 度法求解方程组的近似最优解,将近似最优解应用与MMSE的估算结果中,获取修正的MMSE 估算结果,降低MMSE信道估计的复杂度。
[0009] 所述步骤3中,通过共辆梯度法求解方程组的近似最优解时,通过对方程组的系数 矩阵进行分裂,分裂为对角线上的元素构成的对角阵、严格下Ξ角矩阵W及严格上Ξ角矩 阵,设计由所述对角线上的元素构成的对角阵、严格下Ξ角矩阵W及严格上Ξ角矩阵组成 的预处理矩阵,根据设定的预处理矩阵通过预处理共辆梯度法求解方程组的近似最优解, 得到修正后的MMSE估算结果。
[0010] 对于所述大规模的ΜΙΜΟ系统,接收端配置化根天线,发射端配置有化根天线,通过 发射长度为Β的预定义导频序列Ρ,得到接收信号Υ:
[0011] Υ = ΗΡ+Ν
[0012] 对上述的Η、Ν、Υ、Ρ进行向量化W后,得到步骤1中的基于克罗内克模型的接收信号 模型为:
[0013] y 二 ph + n
[0014] 其中,索为Y的向量化表示,?为Η的向量化表示,fi为N的向量化表示,g为P的向量 化表示,且P二护ΦI,I为一个化X化的单位阵,@素示Kronecker乘积,N为包含噪声和导 频污染的干扰项。
[0015] 对基于克罗内克模型的接收信号模型,根据MMSE估计原理,得到的步骤2中的包含 有对信号协方差矩阵和干扰协方差矩阵进行求逆运算的MMSE估计结果为:
[0016]
[0017]其中,S为干扰协方差矩阵,R为信道的协方差矩阵。
[001引记所述MMSE估计结果中的求逆运算为A,即:
[0019]
[0020] 此时,A-1歹的求逆运算转化为求解线性方程组Αχ = y的解,通过步骤3所述的共 辆梯度法通过L次迭代得到方程组的近似最优解进而得到修正的MMSE估算结果:
[0021]
[0022] 所述步骤3中,通过预处理共辆梯度法求解方程组的近似最优解用W加快共辆梯 度收敛速度的过程为:
[0023] 通过对系数矩阵A进行分裂:
[0024] A = D-L-U
[0025] 其中,D,-L,-U分别是A对角线上的元素构成的对角阵,严格下Ξ角矩阵W及严格 上Ξ角矩阵;
[00%] 设计预处理矩阵W为:
[0038] 然后通过L次迭代得到上述方程组的近似最优解X*,得到修正的MMSE估计结果:
[0039]
[0040] 所述通过共辆梯度法或预处理共辆梯度法获取修正的MMSE估算结果的方法包括 W下步骤:
[0041 ] (1)选取任意初始向量XI和精度要求ε,并置k= 1;
[0042] (2)计算觀=Αχ* -f。.若豁; < 么',则停止计算,此时x* = xk,否则转下一步;
[004引(3)构造捜索方向:令dk 二-gk+0k-ldk-l ,其中,当 k = 1 时,0k-l 二 0 , dk 二-gl ,否则
[0044] (4)-维捜索:沿着dk方向进行捜索,通过迭代步长
:得到Xk+i = Xk +Akdk,k = k+1;然后转(2);
[0045] 对上述过程经过L次迭代得到方程组的近似最优解然后可W得到修正的MMSE 估计结果:
[0046]
[0047] 其中,XI为X的初始向量;ε为共辆梯度法截止精度;K为迭代次数;阱为误差向量;F 为F范数;dk为共辆梯度法捜索方向;β为信道干扰系数;λ为共辆梯度法迭代步长。
[0048] 本发明的有益效果:本发明首先通过对基站接收信号进行向量化,得到基于克罗 内克模型的MMSE信道估计结果。分析信道估计结果表明,由于MMSE设及到协方差矩阵的求 逆运算导致信道估计的计算复杂度为本发明通过将求逆运算转换为求解方程组的 问题,然后利用共辆梯度法求得方程组的近似最优解,W达到降低计算复杂度的目的。同时 为了加快共辆梯度法的收敛速度,设计预处理矩阵,依此来更大幅度的降低信道估计的计 算复杂度,来达到估计性能和计算复杂度之间的平衡。
【附图说明】
[0049] 图1为本发明的流程图。
[0050] 图2为不同估计器的MMSE在不同干扰情况下随着迭代次数的变化情况,其中,(a) 为干扰为零(即β = 〇)的情况,(b)为干扰为0.1(即β = 0.1的情况),(c)为干扰为1(即β=1) 的情况。
[0051] 图3为不同估计器的MMSE在不同干扰情况下随着SNR的变化情况,其中(a)为CG估 计器的情况,(b)为PCG估计器的情况。
【具体实施方式】
[0052] 下面结合附图1~3和【具体实施方式】对本发明作进一步详细说明。
[0053] 如图1所示,本发明包括W下步骤:
[0054] S1:通过对基站端接收信号进行向量化,得到基于克罗内克模型的接收信号模型;
[0055] 对于一个大规模ΜΙΜΟ系统,接收端配置有化根天线,发射端配置有化根天线,通过 发射长度为Β的预定义导频序列Ρ,得到接收信号Υ:
[0056] Υ = ΗΡ+Ν
[0057] 其中Ν € (C^'xs为循环对称复高斯干扰:齡(Ν) e £W(0,S),干扰协方差矩阵 SeCAWxAvs是正定的,它包括常规的非相关接收机噪声和不同类型的来自其他系统的干扰。 引入向量化算子:对H,N,Y分别进行向量化得到:衣=v'er(Y),&二. η二vv【-:(N),哀和 好均为一个化Β X 1的向量,?为一个化化X 1的向量。用同样的方法,我们定义一个导频序列 矩阵:
[0化引
[0化9] I为一个化X化的单位阵,@表示Kronecker乘积。则由W上几个式子可已将信道 模型表示为:
[0060]
[0061] S2:基于上述接收信号模型,根据最小均方误差(MMSE)估计原理得到MMSE估计结 果:
[0062]
[0063] 然后得到MMSE的估计误差:
[0064]
[0065] 由MMSE的估计结果可W看到由于需要对信道协方差矩阵和干扰协方差矩阵进行 求逆运算,而导致MMSE信道估计的计算复杂度为〇(/V/3)。
[0066] S3:将求逆运算转化为求解线性方程组的问题,通过共辆梯度法求得方程组的近 似最优解;
[0067] 记
[006引
[0069] A-1家就可W转化为求解线性方程组Ax二家的问题。然后通过共辆梯度法通过L次 迭代得到方程组的近似最优解
[0070] 算法一为使用共辆梯度法求解方程组的近似最优解的过程:
[0071 ] (1).选取任意初始向量XI和精度要求ε,并置k = 1。
[0072] (2).计算
,:则停止计算,此时x* = xk,否则转下一步。
[007;3] (3).构造捜索方向:令 dk = -gk+0k-idk-i,其中,当 k=l 时,0k-i = 〇,dk = -gi,否则
D.
[0074] (4). -维捜索:沿着dk方向进行捜索,通过迭代步长
,得到Xk+i = Xk +Akdk,k = k+1。然后转第2步。
[0075] 通过算法一经过L次迭代得到方程组的近似最优解然后可W得到修正的MMSE 估计结果:
[0076]
[0077] S4:为了加快共辆梯度法的收敛速度,设计预处理矩阵,在方程组两端同时左乘预 处理矩阵的逆;
[0078] 对系数矩阵A进行分裂:
[0079] A = D-L-U
[0080] 其中,D,-L,-U分别是A对角线上的元素构成的对角阵,严格下Ξ角矩阵W及严格 上Ξ角矩阵。在此设计预处理矩阵W:
[0081 ] W=(D-L)D-i(D-L)h
[0082]同时在方程组两端左乘
[0090] S5:通过预处理共辆梯度法求得方程组的近似最优解,得到修正之后的MMSE估计 结果。
[0091 ]算法二为使用预处理共辆梯度法求解方程组的近似最优解的过程:
[0092] (1).选取任意初始向量XI和精度要求ε,并置k = 1。
[0093] (2).计莫
,则停止计算,此时/ = xk,否则转下一步。
[0094] (3).构造捜索方向:令 dk = -gk+0k-idk-i,其中,当 k=l 时,0k-i = 〇,dk = -gi,否则
[0095] (4). 一维捜索:沿着dk方向进行捜索,通过迭代步长
,得到xk+i = xk+Akdk,k = k+1。然后转第2步。
[0096] 通过算法二经过L次迭代得到方程组的近似最优解然后可W得到修正的MMSE 估计结果:
[0097]
[0098] 综上所述,本发明提出了一种在大规模ΜΙΜΟ系统中低复杂度的信道估计算法,可 W使得传统的MMSE估计器的计算复杂度降低一个量级,并且在存在导频污染的情况下,随 着迭代次数的增加,其估计精度完全可W达到MMSE的估计精度,达到了估计性能和计算复 杂度之间的平衡。
[0099] 本发明通过将信道估计中的求逆运算转化为解线性方程组的问题,然后通过共辆 梯度法求得方程组的近似最优解,并且为了加快共辆梯度法的收敛速度,对系数矩阵A进行 分裂,设计预处理矩阵W,同时在方程组两端左乘Fi,然后在通过L次迭代求得方程组的近似 最优解
[0100] 传统的MMSE估计器计算复杂度为Μ的立方级。对于基于预处理共辆梯度法的信道 估计方法,其计算复杂度大小主要是在计算得到近似最优解X巧寸产生的计算量,并且X*是通 过L次迭代得到;在每一次的迭代过程中,由于λ谢化都是一个标量,那么在计算xk+1和dk+拥 的计算复杂度就只设及到标量与矩阵相乘运算W及矩阵相加运算,其计算量是可W忽略不 计的。也就意味着只有在计算求得gk,Ak和化的过程中才会设及到计算复杂度的问题,但是 由于gk,Ak和故公式中设及到都是向量进行运算,只有一小部分的矩阵与向量的相乘,其计 算复杂度为〇?\<'ν/2)。虽然在计算Ck的过程中包含W-1,但是对于c=W-ig,其相当于Ξ角分 解法中的前代与回代过程,计算复杂度为因此,在利用基于预处理共辆梯度法的 估计器整体计算复杂度为.其中Μ为经过克罗内克模型变换之后信道协方差矩阵 的维度,L为迭代次数,在高效集成电路模型中,运是一个极其容易实现并执行的求解过程。 与此相反,传统的MMSE估计器中设及到矩阵求逆的过程,运个在硬件实现中是一个极其复 杂的过程。由此可见基于预处理共辆梯度法的信道估计算法,可W使传统的MMSE信道估计 的计算复杂度降低一个量级。
[0101] 图2,图3为本发明的一种大规模ΜΙΜΟ系统的实施例,在本实施例中,接收天线数化 = 100,发送天线数化= 10,导频序列长度Β=10。为了不失一般性,设定零均值的信道和干 扰。为了更好地体现出信道的相关特性,遵循克罗内克模型来描述目标信道和干扰信道天 线之间的相关性:
[0102]
[0103] Hw是独立同分布的随机矩阵,矩阵中所有元素均服从期望为0,方差为1分布。在性 能仿真中,大规模ΜΙΜΟ系统中所有的协方差矩阵被模拟为:/?=保,Θ/?,,干扰小区和目标 小区模型一样,第i个干扰小区的协方差矩阵天
3>〇,;'€?,讀为干 扰小区集合,β因子表示该污染小区的污染严重程度。β = 〇时说明没有受到相邻小区的干 扰,表示是噪声受限的情况,β=1说明干扰小区的污染程度比较严重,在此假设一共有两个 干扰小区。
[0104] 定义归一化的导频信噪比
[0112]图2反映了不同估计器的MSE与迭代次数L之间的关系。在此设置信噪比为5地。从 该图中可W看出,无论是噪声受限情况下还是存在导频污染的情况下,随着迭代次数的增 加,CG和PCG估计器的MSE会逐渐降低,最终都能够达到MMSE估计器的性能。而且PCG估计器 的Μ沈渐近于MMSE估计器的速度非常快,从图中就可W看出,无论是噪声受限还是存在导频 污染的情况下,PCG估计器最多经过4次迭代,其估计精度就可W达到MMSE的估计精度,而CG 估计器基本需要经过16次迭代,也可W完全达到MMSE的估计精度。
[0113] 图3中主要描述了不同估计器在迭代次数固定的情况下,不同干扰情况下的MSE随 着信噪比的变化情况。在图3(a)中,设置迭代次数为固定值L=15,图3(b)中,设置迭代次数 为固定值L = 3,从图中可W发现,在噪声受限的情况下,CG和PCG估计器的MSE不能够趋近于 MMSE估计器。但是在存在导频污染的情况下,CG估计器和PCG估计器的MSE随着信噪比的增 加,会逐渐趋于稳定,但是相比于MMSE估计器,估计精度会稍微有所下降,不过基本可W认 为与匪SE估计性能相当。虽然在迭代次数固定的情况下,CG和PCG估计性能会稍微有所下 降,但是却可W使信道估计的计算复杂度降低一个量级,达到了估计性能和计算复杂度之 间的平衡。
[0114] 本发明中各字符含义如下:
[0115] 化-接收端天线数目;
[0116] Nt-发射端天线数目;
[0117] P-预定义导频序列;
[011引B-预定义导频序列P的长度;
[0119] Y-接收端获取的接收信号;
[0120] H-信道矩阵;
[0121] N-包括噪声和导频污染的干扰项;
[0122] S-干扰协方差矩阵;
[0123] y-Y的向量表示;
[0124] h-H的向量表示;
[0125] fi-N的向量表不;
[0126] P-P的向量表示;
[0127] I-NrX化的单位阵;
[012引 ρΤ-Ρ的转置矩阵;
[01巧]良"_ρ的共辆转置矩阵;
[0130] R-信道的协方差矩阵;
[0131] Α-系数矩阵;
[0132] X*-Αχ二f的近似最优解;
[0133] X1-X的初始向量;
[0134] ε-共辆梯度法截止精度;
[0135] Κ-迭代次数;
[01%] gk-误差向量;
[0137] F-F 范数;
[0138] 也-共辆梯度法捜索方向;
[0139] β-信道干扰系数;
[0140] λ-共辆梯度法迭代步长;
[0141] D-A的对角线上的元素构成的对角阵;
[0142] -L-A的严格下Ξ角矩阵;
[0143] -U-A的严格上Ξ角矩阵;
[0144] w-预处理矩阵;
[0145] R广接收相关矩阵;
[0146] Hw-元素服从(0,1)独立同分布的随机矩阵;
[0147] Rt-发射相关矩阵。
【主权项】
1. 大规模ΜΙΜΟ系统中低复杂度信道估计方法,其特征在于,该方法包括以下步骤: 步骤1:对基站接收的信号进行向量化,得到基于克罗内克模型的接收信号模型; 步骤2:对上述接收信号模型,根据MMSE估计原理,得到包含有对信号协方差矩阵和干 扰协方差矩阵进行求逆运算的MMSE估计结果; 步骤3:对MMSE估计结果中的求逆运算转换为求解线性方程组的问题,通过共辄梯度法 求解方程组的近似最优解,将近似最优解应用与MMSE的估算结果中,获取修正的MMSE估算 结果,降低MMSE信道估计的复杂度。2. 根据权利要求1所述的大规模ΜΙΜΟ系统中低复杂度信道估计方法,其特征在于:所述 步骤3中,通过共辄梯度法求解方程组的近似最优解时,通过对方程组的系数矩阵进行分 裂,分裂为对角线上的元素构成的对角阵、严格下三角矩阵以及严格上三角矩阵,设计由所 述对角线上的元素构成的对角阵、严格下三角矩阵以及严格上三角矩阵组成的预处理矩 阵,根据设定的预处理矩阵通过预处理共辄梯度法求解方程组的近似最优解,得到修正后 的MMSE估算结果。3. 根据权利要求1所述的大规模ΜΙΜΟ系统中低复杂度信道估计方法,其特征在于:对于 所述大规模的ΜΜ0系统,接收端配置Nr根天线,发射端配置有Nt根天线,通过发射长度为Β 的预定义导频序列P,得到接收信号Y: Y = HP+N 对上述的H、N、Y、P进行向量化以后,得到步骤1中的基于克罗内克模型的接收信号模型 为:其中,f为Y的向量化表示,ft为Η的向量化表示,fi为N的向量化表示,豆为P的向量化表 示,且p_ = P7?I .,I为一个Nr XNr的单位阵,?表示Kronecker乘积,N为包含噪声和导频污 染的干扰项。4. 根据权利要求3所述的大规模ΜΙΜΟ系统中低复杂度信道估计方法,其特征在于:对基 于克罗内克模型的接收信号模型,根据MMSE估计原理,得到的步骤2中的包含有对信号协方 差矩阵和干扰协方差矩阵进行求逆运算的MMSE估计结果为:其中,S为干扰协方差矩阵,R为信道的协方差矩阵。5. 根据权利要求4所述的大规模MM0系统中低复杂度信道估计方法,其特征在于:记所 述MMSE估计结果中的求逆运算为A,即: A-pRp 4 S 此时,的求逆运算转化为求解线性方程组Ax = f的解,通过步骤3所述的共辄梯 度法通过L次迭代得到方程组的近似最优解X'进而得到修正的MMSE估算结果: h = Rp7 Xs .6. 根据权利要求2或5所述的大规模ΜΙΜΟ系统中低复杂度信道估计方法,其特征在于: 所述步骤3中,通过预处理共辄梯度法求解方程组的近似最优解用以加快共辄梯度收敛速 度的过程为: 通过对系数矩阵A进行分裂: A = D-L-U 其中,D,-L,-U分别是A对角线上的元素构成的对角阵,严格下三角矩阵以及严格上三 角矩阵; 设计预处理矩阵W为: w=(d-l)d_1(d-l)h 在方程组两端左乘F1后 ff_1Ax = ff_1Aff_HffHx=ff_1b 记 ff_1Ax = ff_1Aff_HffHx=ff_1b 记然后通过L次迭代得到上述方程组的近似最优解,,得到修正的MMSE估计结果:6. = 。7.根据权利要5或6所述的大规模ΜΙΜΟ系统中低复杂度信道估计方法,其特征在于: 所述通过共辄梯度法或预处理共辄梯度法获取修正的MMSE估算结果的方法包括以下 步骤: (1) 选取任意初始向量XI和精度要求ε,并置k= 1; (2) 计算gA = AXa- -f。若H ,则停止计算,此时x* = xk,否则转下一步; (3 )构造搜索方向:令dk = -gk+0k-ldk-1,其中,当k = 1时,0k-1 = 0,dk = -gi,否则(4)一维搜索:沿着dk方向进行搜索,通过迭代步长得到xk+i = xk+X kdk,k = k+l;然后转(2);对上述过程经过L次迭代得到方程组的近似最优解X'然后可以得到修正的MMSE估计结 果: h-Rp'x 其中,XI为X的初始向量;ε为共辄梯度法截止精度;K为迭代次数;gk为误差向量;F为F范 数;dk为共辄梯度法搜索方向;β为信道干扰系数;λ为共辄梯度法迭代步长。
【文档编号】H04B7/08GK106059645SQ201610332896
【公开日】2016年10月26日
【申请日】2016年5月18日
【发明人】高向川, 李臣阳, 王树坤, 张卫党, 王法松, 朱政宇, 李青
【申请人】郑州大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1