法律状态公告日
法律状态信息
法律状态
2020-07-14
授权
授权
2017-06-23
实质审查的生效 IPC(主分类):H04B7/0456 申请日:20161230
实质审查的生效
2017-05-31
公开
公开
技术领域
本发明属于移动通信领域,主要涉及基于改进的牛顿迭代的大规模MIMO预编码方法。
技术背景
大规模MIMO系统是第五代移动通信系统的关键技术之一,通过在基站和用户端配备大量天线,显著地提高系统的信道容量、数据传输速率、频谱效率和通信质量。由于随着用户数目的增加,用户接收信号中存在大量的用户间干扰,因此在发射端引入预编码技术。在下行链路中,基站首先通过信道估计得到信道状态信息(CSI,channel stateinformation),然后基于CSI将调制后的发射信号进行预编码,以便消除用户间干扰,匹配信道衰落。
传统的预编码算法可以根据运算特性分为线性预编码和非线性预编码两类。线性预编码是使用线性运算处理信息,算法实现简单,计算复杂度低。主要有迫零预编码(ZF)、匹配滤波预编码(MF)和最小均方误差预编码(MMSE)等。非线性预编码是使用非线性运算处理信息,可以获得很好的系统性能,但是计算复杂度过高,不适用于大规模MIMO系统。
在线性预编码中,基于规则化迫零(Regularized ZF,RZF)预编码具有很好的性能和较低的复杂度,它是在ZF预编码和MMSE预编码算法的基础上提出来的。但是随着天线数目的增加,RZF预编码存在高维矩阵求逆的过程,具有很高的计算复杂度。
发明内容
本发明为了解决现有技术中的问题,提供了一种基于改进的牛顿迭代的大规模MIMO预编码方法,该方法采用牛顿迭代法对高维矩阵求逆过程进行估计,将矩阵求逆过程转化为矩阵乘法和矩阵加法的迭代过程。并且对牛顿迭代法进行改进,加快收敛速度,减少迭代次数。利用改进的牛顿迭代法对RZF预编码进行估计,算法收敛速度快,在迭代次数较小的情况下得到性能优的预编码矩阵。
本发明所述的基于改进的牛顿迭代的大规模MIMO预编码方法包括:
(1)估计大规模MIMO的信道矩阵,通过得到的信道矩阵计算RZF预编码矩阵;
(2)采用改进的牛顿迭代法对RZF预编码矩阵中的逆矩阵进行估计,将矩阵求逆运算转化成矩阵加法和矩阵乘法运算,其中,改进的牛顿迭代法是构建高阶迭代式,将位于0附近的特征值经过变换,使其更加靠近1,1附近的特征值保持不变,从而加快牛顿迭代的收敛速度;
(3)利用得到的预编码矩阵对发送信号进行预编码。
其中,步骤(1)具体包括:
(1-1)设置大规模MIMO发射天线数量为M,接收天线数量为K,信道为慢衰落信道,则信道向量为:
hk~CN(0M×1,Φ),k=1,...,K
式中,hk表示第k个信道向量,hk~CN(0M×1,Φ)表示hk服从均值为0M×1方差为Φ的分布,0M×1表示M行1列的0矩阵,Φ是信道的相干矩阵,具有有界的谱范数;
(1-2)对信道向量进行估计,得到信道向量的估计值为:
式中,zk是与hk独立同分布的信道估计噪声,zk~CN(0M×1,Φ),τ代表信道估计的质量,如果τ值为0则表明系统拥有高质量的信道估计,τ值为1则表明系统得到的信道估计没有任何意义;
(1-3)根据信道向量的估计值得到信道矩阵的估计值
(1-4)根据信道矩阵的估计值
其中,步骤(2)具体包括:
(2-1)令矩阵
(2-2)按照以下步骤进行迭代:
a)令i=1;
b)计算Ti=(2I-Ti-1X)Ti-1;
c)令Ai=TiX,判断
d)计算Ti+1=(γ3(TiX)2+γ2(TiX)+γ1I)Ti,其中,
e)令Ai+1=Ti+1X,i=i+1;
f)i=i+1,并返回至b),直至迭代次数达到预设次数停止,求得最后的X;
(2-3)将求得的X带入公式
其中,步骤(3)具体包括:
利用得到的预编码矩阵GRZF对发送信号进行预编码,得到编码后的信号为x=GRZFS,式中,S表示发送信号矩阵。
本发明的有益效果:实验结果表明,当迭代次数超过4次时,传统的牛顿迭代法的性能优于基于泰勒级数展开的逆矩阵估计算法。在迭代次数为2时,改进的牛顿迭代优化算法就已经能获得大约95%的RZF预编码平均用户到达率。
附图说明
图1是本发明在信道估计误差为0.1,发送天线数为256,接收天线数为32,截断阶数和迭代次数均为4的条件下,对截短泰勒级数法和传统牛顿迭代法得到的预编码矩阵获得的用户到达率与RZF预编码用户到达率进行比较。
图2是本发明在信道估计误差为0.1,发送天线数为256,接收天线数为32,截断阶数和迭代次数均为2的条件下,对传统牛顿迭代法和改进的牛顿迭代法得到的预编码矩阵获得的用户到达率与RZF预编码用户到达率进行比较。
具体实施方式
下面对具体实施进行详细的描述。
步骤(1)
设信道为慢衰落信道,信道矩阵建模为:
hk~CN(0M×1,Φ)k=1,...,K(1)
式中,
式中,hk为实际信道,zk是与hk独立同分布的信道估计噪声,zk~CN(0M×1,Φ),τ代表信道估计的质量,如果τ值为0则表明系统拥有高质量的信道估计,τ值为1则表明系统得到的信道估计没有任何意义。最终估计的信道矩阵为
第k个用户接收到的信号可以表示为:
式中,nk是循环对称复高斯随机噪声,其均值为0,方差为σ2,记为nk~CN(0,σ2)。
其中β是使等式tr(GGH)=P成立的系数,P为基站的发送功率,ξ是优化系数的最优值。
步骤(2)
令矩阵
如表1所示,按照以下步骤进行迭代:
a)令i=1;
b)计算Ti=(2I-Ti-1X)Ti-1;
c)令Ai=TiX,判断
d)计算Ti+1=(γ3(TiX)2+γ2(TiX)+γ1I)Ti,其中,
e)令Ai+1=Ti+1X,i=i+1;
f)i=i+1,并返回至b),直至迭代次数达到预设次数停止,求得最后的X;
将求得的X带入公式
表1
其中,步骤(2)的推算过程为:
设令矩阵
Ti=(2I-Ti-1X)Ti-1(6)
其中,当T0=α0XT时,上式收敛。
对矩阵X进行奇异分解可得:X=U∑UT,其中
其中
由(8)式可得如下递推关系
化简得
其中
经过推导,α0应该满足
(11)式需要计算矩阵X的特征值,因此下面给出一个次优解。设λ是矩阵
因此可以做如下近似近,
由(10)式可知,当pj0在0和2附近时,pji向1的收敛速度很慢。经过多次的牛顿迭代最终pji值分布范围主要集中在0和1附近,因此使用构建高阶迭代式进一步提高牛顿迭代的收敛速度。将位于0附近的特征值经过变换,使其更加靠近1,而1附近的特征值继续保持不变,通过这种方式进而加速牛顿迭代收敛速度。
令Ai=TiX,
Ti+1=(γ3(TiX)2+γ2(TiX)+γ1I)Ti(13)
为使得上述方程满足如下条件:f(1)=1,f'(1)=0,f(0)=0,f'(0)>>2,方程解得
上述特征值方程能够保证f'(0)>>2,f'(1)=0,所以pj在0附近的值能够被放大,pj在1附近的值能够继续保持收敛。
步骤(3)
利用得到的预编码矩阵GRZF对发送信号进行预编码,得到编码后的信号为x=GRZFS,式中,S表示发送信号矩阵。
机译: 一种无线通信系统中的预编码方法,该方法可获取基于代码簿的预编码矩阵索引和变化的多级预编码解决方案,便携式终端及其基站
机译: 下行大规模MIMO系统的高效预编码装置及其方法
机译: 大规模MIMO稳健预编码传输方法