利用矩阵理论详细推导MIMO信道容量

利用矩阵理论详细推导MIMO信道容量摘要多输入多输出(MIMO)技术被认为是现代通信技术中的重大突破之一，以其能极大增加系统容量与改善无线链路质量的优点而受到了越来越多的重视与关注。通信信道容量是信道进行无失真传输速率的上界，因此研究MIMO的信道容量具有巨大的指导意义。本文把矩阵理论知识与MIMO技术信道容量中的应用紧密结合，首先建立了MIMO信道模型，利用信息论理论和矩阵理论详细推导出MIMO信道容量。并得出重要结论。关键词：MIMO；信道容量；奇异值分解一、引言MIMOMultipleInput-MultipleOutput)是指在通信链路的发送端与接收端均使用多个天线元的传输系统,它能够将传统通信系统中存在的多径因素变成对用户通信性能有利的因素,从而成倍地提高业务传输速率。矩阵理论在通信，自动控制等工程领域里应用广泛，而通信的难点在于信道的处理，因此，矩阵理论与无线信道的研究是一个很好的切入点。目前，MIMO技术的信道容量和空时编码，空时复用等技术都离不开矩阵理论的应用。二、利用矩阵理论详细推导MIMO信道容量1)MIMO信道介绍MIMO是多输入多输出系统，当发送信号所占用的带宽足够小的时候，信道可以被认为是平坦的，这样，MIMO系统的信道用一个RTnn的复数矩阵H描述，H的子元素,jih表示从第(1,2,...)Rjjn根发射天线到第(1,2,...)Tiin根接收天线之间的空间信道衰落系数[1]。如下图所示：111212122212TTRTRRnnnnnnHhhhhhhhhh(2.1)每个符号周期内，发送信号可以用一个1Tn的列向量12[]TTinxxxxx表示，其中ix表示在第i个天线上发送的数据。同时，用一个1Rn的列向量12[]RTinyyyyy表示，其中iy表示在第i个天线上发送的数据。对于高斯信道,发射信号的最佳分布也是高斯分布[1]。因此,x的元素是零均值独立同分布的高斯变量。发送信号的协方差可以表示为：{}HxxRExx(2.2)发送信号的功率可以表示为()xxPtrR(2.3)接收信号和噪声可以分别用两个1Rn的列向量y和n表示。其中信道噪声是加性噪声，服从循环对称复高斯分布，并且与发射信号x不相关，假设n均值为0，功率为2。噪声的协方差为：2RHnnnREnnI(2.4)通过这样一个线性模型，接收信号可以表示为yHxn(2.5)??????????????接收信号的协方差可表示为2[]RHyyHHHHHHHxxnREyyEHxnHxnHExxHEnnHExnHEnxHRHI(2.6)因为x与噪声n不相关，所以0HExnEnx。2)MIMO信道容量一般公式推导下面根据信息论知识，我们对MIMO信道容量做一般性推导。在下面的推导过程中我们假设信道矩阵H在接收端已经完全已知，但是它是随机的，因此我们可以得到瞬时信道容量为：[1]()max,XxCHIxy(2.7)其中，,Ixy是在已知信道H的情况下输入x与输出y之间的互信息量，有：,|IxyHyHyx(2.8)其中，Hy是y的信息熵（微分熵），定义：2()()log()Hypypy，其中()py是y的概率（概率密度）。Hy是y的差分嫡，(|)Hyx是给定x条件下y的差分嫡，由于发送信号与噪声之间是独立的，因此有(|)()HyxHn[1]，所以上式可以重新写为：,IxyHyHn(2.9)由于噪声概率密度函数确定，所以Hn为定值，当信道为加性高斯信道时，信源x服从高斯分布时此时接收信号y也服从高斯分布，根据信息论理论，此时(,)Ixy取最大，即为信道容量。此时y和n的信息熵分别为：[1]212()logdetyybitHyeR(2.10)2212()logdetRnbitHneI(2.11)所以我们可以得到信道瞬时交互信息(,)Ixy，也即信息容量为：222222221()logdet/det21logdet()/det21logdet2RRRRRyynHxxnnnHxxnCHeReIHRHeIIeIHRHIbit(2.12)工程中一般定义信道容量为单位时间内的平均互信息，故定义MIMO的信道容量：[2]1CCHT(2.13)其中T为一个符号周期，根据采样定理，(1/)2TB，其中B为信号带宽，取(1/)2TB，代入（2.13）式，得：[2]22logdet/RHxxnHRHCBIbits(2.14)这即是MIMO的信道容量一般公式。3)奇异值分解计算MIMO信道容量对于MIMO无线信道，信道是极其复杂的。因此原始的信道矩阵也就显得复杂，不便于分析，而且一般矩阵不经过处理计算行列式很困难。这就自然想到在信源端对发射信号做某种预处理，使得经过预处理的信号经过的信道变得简单易分析，而且具体实现也变得简单。对于信道矩阵来说，对角矩阵是最简单的，所以自然就想到把信道矩阵分解，利用矩阵理论中的奇异值分解可以达到这种目的。下面就矩阵的奇异值分解来计算MIMO的信道容量。由奇异值分解理论,任何一个RTnn的矩阵H可以写成：[3]HHUDV(2.15)式中，D是RTnn的对角阵，其对角非零元素模值为H的正奇异值，U和V分别为RRnn和TTnn的酉矩阵。把公式（2.15）代入公式（2.5），得：HyUDVxn(2.16)公式（2.16）两边同时左乘酉矩阵HU，利用HRnUUI，同时令HyUy，HxVx即对信源信号x进行预编码，HnUn，得：[3]yDxn(2.17)根据矩阵理论，矩阵HHH的特征值为非负数，非零特征值的个数等于矩阵H的秩，用r表示，对于RTnn的矩阵H，秩的最大值min(,)RTmnn，用i表示矩阵H的奇异值，将i代入公式（2.17）得：(1,2,...)iiiiyxnir(2.18)(1,...)iiRynirn(2.19)式（2.19）显示，接收元素(1,...)iRyirn并不依赖于发射信号，即信道增益为零。另一方面,接收元素(1,2...,)iyir仅仅取决于发射元素ix，因此,可以认为,通过式（2.18）和（2.19）得到的等效MIMO信道是由个去耦平行子信道组成的[4]。其信道增益为矩阵H的奇异值。可以进一步推导出iy，ix，和in的协方差和迹：,,HHHyyxxnnyyxxnnRURURVRVRURU(2.20)()(),()(),()()yyxxnnyyxxnntrRtrRtrRtrRtrRtrR(2.21)对于接收端已知信道参数，发射端未知，发射端平均分配功率的MIMO系统。设发射端发射总功率为TP，则每个天线发射功率为/TTPn，此时{}HTxxTTnPRExxIn(2.22)根据式（2.20）得：HTxxxxTTnPRVRVIn(2.23)然后把根据（2.15）算出HHH后和（2.23）代入MIMO信道一般公式（2.14）得：221log(1)/riTiTPCBbitsn(2.24)式（2.24）可以看出,MIMO链路的信道容量很大程度上取决于H的秩r。矩阵的秩越大，容量也越大。所以，MIMO正是利用无线信道的多径效应使相距超过半个波长的天线尽量不相关，从而使信道矩阵秩越大，进而在不增加带宽和发射功率的情况下增加系统容量。[4]下面我们说明信道容量是如何与信道矩阵H相关的。令min(,)RTmnn为矩阵H的秩，定义：,,HRTHRTHHnnQHHnn(2.25)设(1,2,...)iim是Q的特征值，则：1det()()mmiiIQ(2.26)用2TTnP代替（2.26）中的，得：221det()(1)miTTmiTTpPIQnn(2.27)把式(2.27)代入式（2.24）得平均功率分配的MIMO信道容量：22log[det()]/TmTPCBIQbitsn(2.28)当Tn逐渐增大，使得1mTQIn，此时，可以得到信道容量的近似表达式：22log(1)/TPCBmbits(2.29)由此可以看出当发射天线数目很多时,信道容量随m的增大而线性增大,也就是说可以在不增加带宽和发射功率的情况下成倍的提高信道容量。参考文献：1梁毓锋.MIMO无线通信系统的信道容量分析:[D].大连:大连海事大学，20082ThomasM.Cover，JoyA.Thomas（著），阮吉寿，张华（译）.信息论基础[M].北京:机械工业出版社,20093孙丹，张晓光.MIMO系统信道容量研究[J].现代电子技术,2006(19):4-64刘冰.MIMO系统的信道容量分析[J].微计算机信息,2005,21(12-3):129-131