大规模MIMO传输关键技术——下行链路预编码（1）

在MIMO系统中，如果发送端能够获知信道的某些信息（预编码的前提条件），就可以利用该信息对发射信号进行预处理以提高系统的传输速率和链路可靠性。利用发送端信道状态信息（Channel State Information at the Transmitter，CSIT）对发送信号进行预处理的技术称为预编码技术。

预编码技术是为了解决解码检测（Decoding）技术中的缺陷。解码检测技术的目的是在接收端利用多种信号处理来消除由信道和接收端噪声引起的信号干扰。但是，这种在接收端进行信号处理是存在着巨大的缺陷的，因为接收端必然需要有足够的硬件基础来支撑从而达到信号处理所要求的条件。首先，接收端要进行对接收信号的解码检测，就必须获得相应的及时信道状态信息，这就总味着接收端必须要通过信道估计（Channel Estimate）手段来获取相应的数据信息，但是信道估计技术现阶段仍然存在不够精确或者性能不稳定等缺陷。其次，信号的检测需要大量的复运算，而实际中接收端（即用户）的运算处理能力十分有限。
　　
对于多用户MIMO系统的下行链路，如果不采用预编码技术，基站在同一时频资源上发送多个用户的信号就会造成用户之间的干扰。各个用户受其接收天线数目的限制，很难独自消除来自其他用户的干扰，恢复出所需要的信号。为了解决复用用户之间的干扰问题，基站需要根据 CSIT 对发送信号进行预编码。除此之外，在多用户MIMO系统下行链路中采用预编码技术，通过将大量的复杂计算放置在拥有计算性能较好的发送端，可以降低接收机的复杂度，解决移动台的功耗问题。

预编码技术可以按照编码设计方案分为线性预编码和非线性预编码两类。在多用户MIMO系统中，在基站发送端进行线性处理的发送信号方法称为线性预编码。主要的线性预编码方案有：迫零( Zero Forcing，ZF)预编码、块对角化( Block Diagonalization，BD)预编码、最小均方误差( Minimum Mean Square Error，MMSE)预编码。非线性预编码是指在发送信号上添加设计好的预编码矢量来达到抑制用户间干扰的预编码方法。非线性处理方法主要是基于“脏纸编码( Dirty Paper Coding，DPC) ”思想的各种预处理方案，例如Tomlinson－Harashima Precoding 预编码(THP)以及矢量预编码(Vector Precoding，VP)。尽管这些非线性预编码技术性能指标足够优秀，但是需要的实现条件较高，并且有着极高的计算复杂度，而这些缺陷的存在始终对于硬件实现和运维管理来说是不切实际的。

在常规场景下，非线性预编码性能优于线性预编码，但是线性预编码有复杂度低的优势，更加容易应用于实际场景中。并且大规模MIMO的物理特性（基站的天线数趋近于无穷大时）使得线性预编码技术性能可以无限逼近性能近乎最优的非线性预编码技术，所以大规模MIMO的线性预编码技术己经成为5G通信技术中的研究重点和热点。

线性预编码技术可以划分为基于码本的预编码技术和非码本的预编码技术。基于码本的预编码技术的预编码矩阵产生于用户，即用户通过收到的导频信息来估计获得CSI，并按照一定的准则从预先设计过的预编码矩阵码本中选取最优或者次优的预编码矩阵，最后通过反馈链路传递所选预编码矩阵的索引信息。一套码本的适用场景相对具体和固定，且不同场景下的码本设计差异较大。相对于基于码本设计的预编码技术而言，非码本的线性预编码算法就简单许多。在CSI可用的基础上，基站端负责全部的预编码矩阵的计算，生成及运用等工作。这样就降低了用户端的要求。

基于CSIT，发送端可以对需要发送的数据向量s进行预处理后发送，即发送向量可以表示为：
x=g(s)(1) \boldsymbol x = g(\boldsymbol s)\tag {1} x=g(s)(1)

其中，g(⋅)g(⋅)g(⋅)是预处理函数，如果g(⋅)g(⋅)g(⋅)是非线性函数，则称该处理方法为非线性预编码；如果g(⋅)g(⋅)g(⋅)是线性函数，则称为线性预编码。

CSIT的获取：时分双工（TDD）的通信系统中，信道链路具有良好的互易性，即上下行信道链路的状态信息在一个相干时间内几乎一致，并且上下行链路使用在同一频率内，节省了频谱资源，所以用户只需要通过上行信道向基站端发送导频信号，基站通过从上行信道获得导频信息从而估计得到此时的上行信道的状态信息，再利用信道互易性便可得到下行信道的状态信息。在基站得到其到每个用户的信道状态矩阵后，将根据该矩阵计算出预编码矩阵，最后将编码后的数据映射到天线上进行传输。相比之下，频分双工（FDD）通信系统获得信道状态信息的方式比较复杂，因此较少讨论。

考虑由一个基站和KKK个用户组成的多用户MIMO系统，其中基站配置MMM根天线，用户kkk配置NkN_kNk​根天线（k=1,2,⋯ ,Kk=1, 2,\cdots, Kk=1,2,⋯,K）。基于平坦衰落信道的假设，将用户kkk到基站的信道矩阵记为Gk∈CM×NkG_k\in \mathbb{C} ^ {M\times N_k}Gk​∈CM×Nk​，基站到用户kkk的信道矩阵记为Hk∈CNk×MH_k\in \mathbb{C} ^ {N_k\times M}Hk​∈CNk​×M。（注：当上行链路和下行链路使用相同的信道时，HkH_kHk​为GkG_kGk​的共轭转置）

基站向用户kkk发送Lk(Lk≤Nk,∑k=1KLk=M)L_k(L_k \leq N_k,\sum \limits_{k=1} ^ {K} L_k =M)Lk​(Lk​≤Nk​,k=1∑K​Lk​=M) 路数据流，即sk∈CLk×1\boldsymbol s_k\in\mathbb {C}^{L_k\times1}sk​∈CLk​×1。如果不采用预编码方案，基站的发送信号向量为x=[s1T,s2T,⋯ ,sKT]T(∈CM×1)\boldsymbol x=[\boldsymbol s_1^T,\boldsymbol s_2^T,\cdots,\boldsymbol s_K^T]^T(\in \mathbb{C} ^ {M\times 1})x=[s1T​,s2T​,⋯,sKT​]T(∈CM×1)，用户kkk收到的信号向量为：
yk=Hkx+nk(2) \boldsymbol{y}_k=\boldsymbol {H}_k\boldsymbol x+\boldsymbol{n}_k\tag{2} yk​=Hk​x+nk​(2)
　　显然，用户kkk收到的信号中包含了来自其他用户的干扰，现有的移动通信系统中，由于各用户在地理位置上是分散的，因此一般采用基站集中控制的模式，用户之间并不协同工作。而各个用户的接收天线数目一般小于基站的发送天线数，难以独自消除来自其他用户的干扰。由此可见，在发送端采用预编码对用户间的干扰进行预先消除是非常必要的。

采用线性预编码，基站将sk\boldsymbol s_ksk​经预编码矩阵Wk\boldsymbol {W}_kWk​后发送，发送向量
x=∑k=1KWksk=Ws(3) \boldsymbol x =\sum\limits_ {k=1}^K\boldsymbol {W}_k\boldsymbol s_k=\boldsymbol {Ws}\tag{3} x=k=1∑K​Wk​sk​=Ws(3)
其中，W≜[W1,W2,⋯ ,WK]\boldsymbol W\triangleq[\boldsymbol W_1,\boldsymbol W_2,\cdots,\boldsymbol W_K]W≜[W1​,W2​,⋯,WK​]，s≜[s1T,s2T,⋯ ,sKT]\boldsymbol s\triangleq [\boldsymbol s_1^T,\boldsymbol s_2^T,\cdots,\boldsymbol s_K^T]s≜[s1T​,s2T​,⋯,sKT​]。且向量x\boldsymbol xx满足基站平均发射功率的约束：
E[xHx]=P(4) \mathbb{E}[\boldsymbol x^H \boldsymbol x]=P\tag{4} E[xHx]=P(4)
根据式(1)(2)，用户kkk的接收向量为：
yk=HkWksk+Hk∑j=1,j≠kKWjsj+nk(5) \boldsymbol{y}_k=\boldsymbol {H}_k\boldsymbol {W}_k\boldsymbol s_k+\boldsymbol {H}_k\sum\limits_{j=1,j\neq k}^K\boldsymbol {W}_j\boldsymbol s_j+\boldsymbol{n}_k\tag{5} yk​=Hk​Wk​sk​+Hk​j=1,j​=k∑K​Wj​sj​+nk​(5)
上式中等号右边第一项为用户kkk需要的信号，第二项为其他用户造成的干扰，最后一项是加性高斯白噪声。用矩阵形式表示总的输入-输出关系为：
y=HWs+n(6) \boldsymbol y =\boldsymbol {HWs}+\boldsymbol{n}\tag{6} y=HWs+n(6)

对于下行大规模MIMO信道，系统模型：有LLL个小区，每个小区1个基站，配备MMM根天线，小区中有KKK个用户，每个终端配备NNN根天线。s=[s1,s2,⋯ ,sK]Ts=[\boldsymbol s_1,\boldsymbol s_2,\cdots,\boldsymbol s_K]^Ts=[s1​,s2​,⋯,sK​]T表示发送端发送给用户的信息矢量。预编码矩阵为W=[W1,W2,⋯ ,WK]\boldsymbol W= [\boldsymbol W_1,\boldsymbol W_2,\cdots,\boldsymbol W_K]W=[W1​,W2​,⋯,WK​]，WK\boldsymbol W_KWK​表示第KKK个用户的预编码矩阵。 H≜[h1T,h2T,⋯ ,hKT]T\boldsymbol H\triangleq [\boldsymbol h_1^T,\boldsymbol h_2^T,\cdots,\boldsymbol h_K^T]^TH≜[h1T​,h2T​,⋯,hKT​]T为（基站到）用户的信道矩阵。PPP表示功率矩阵。

ZF（Zero Forcing）预编码，这是一种信道求逆方法，发送端乘以信道矩阵的逆（伪逆）矩阵，使得等效信道是一单位阵。ZF预编码技术的思路就是利用已知信道状态信息，并且基于最小二乘估计，将各独立信号之间的干扰实施线性迫零，从而得到所需信号。在理想CSIT下，预编码矩阵WWW设计为：
W=cHH(HHH)−1(7) \boldsymbol W=c\boldsymbol H^H(\boldsymbol H\boldsymbol H^H)^{-1}\tag{7} W=cHH(HHH)−1(7)
其中，c=P/∥(HHH)−1∥F2c=\sqrt {P/\left\|(\boldsymbol H\boldsymbol H^H)^{-1}\right\|_F^2}c=P/∥∥∥​(HHH)−1∥∥∥​F2​​是功率归一化因子。将（6）带入（5）可得，这种预编码算法可以完全消除用户之间的干扰和用户内部各数据流之间的干扰，即[HW]k,j=0,k≠j[HW]_{k,j}=0,k\neq j[HW]k,j​=0,k​=j，因此称为迫零预编码。迫零预编码的缺点是会使加性噪声被加权放大。⟹\Longrightarrow⟹规则化逆信道（Regularized Channel Inversion，RCI）方法避免了迫零预编码中的噪声放大问题⟹\Longrightarrow⟹规则化迫零预编码（RZF）。

区别于ZF 预编码方案将其他用户的干扰化为零，MMSE 预编码方案以接收端的SINR（信干噪比）达到最大为目的。对于秩亏的信道，迫零预编码下的性能可以通过正则化伪逆获得提升，即：
W=HH(HHH+βI)−1(8) \boldsymbol W=\boldsymbol H^H(\boldsymbol H\boldsymbol H^H+\beta I)^{-1}\tag{8} W=HH(HHH+βI)−1(8)
　　其中，β\betaβ是正则化因子。这种预编码算法常被称为MMSE预编码，因为它与白噪声情况下MMSE波束成形设计原则相类似。通过引入参数β\betaβ ，使得系统的和容量随着$min( M, N ) $线性增长，并且MMSE 预编码方案在低信噪比的情况下性能更好。

在预编码信号处理中，基于SVD准则的预编码充分地利用信道增益和信道资源信息，将整个信道拆分为多个互不干扰的且具有不同增益的子信道。基于这样的思想，SVD准则的预编码可以有效的消除ISI。

首先对信道矩阵H∈CN×M\boldsymbol H\in \mathbb{C} ^ {N \times M}H∈CN×M进行奇异值分解：H=UΣVH\boldsymbol H=\boldsymbol U\boldsymbol \Sigma \boldsymbol V^HH=UΣVH；式中，Σ\boldsymbol \SigmaΣ是N×MN×MN×M非负对角矩阵，U\boldsymbol UU和V\boldsymbol VV分别是N×NN×NN×N和M×MM×MM×M的酉矩阵。此时，该技术对应的最优预编码矩阵为矩阵V\boldsymbol VV中那些奇异值不为0的列的构成的矩阵：
W=V~P1/2(9) \boldsymbol W= \boldsymbol {\widetilde{V} P^{1/2}\tag{9}} W=VP1/2(9)
其中，V~∈CM×N\boldsymbol {\widetilde{V} }\in \mathbb{C} ^ {M\times N}V∈CM×N是矩阵中V\boldsymbol VV非零奇异值对应的奇异值向量构成的矩阵。P=diag(p1,p2,…,pN)\boldsymbol P=diag(p_1,p_2,\dots,p_N)P=diag(p1​,p2​,…,pN​)是功率分配矩阵，并且矩阵元素满足功率限定条件tr(P)=Ptrtr(\boldsymbol P)=P_{tr}tr(P)=Ptr​。

BD线性预编码方案：在MU-MIMO系统中，引入BD预编码。它采用收发端协调波束成形技术，同时消除用户间和同一个用户各天线间的干扰，其主要思想是将下行信道分解为多个独立的平行或正交的单用户MIMO信道。

通常，BD预编码算法分两个步骤实现。第一步为除第k个用户信道外的总信道矩阵H~k\widetilde H_kHk​的零空间计算一个正交基Wa\boldsymbol W ^aWa，可用于消除MUI，获得等效的并行SU-MIMO信道。 第二步用于并行化每个用户流，并通过Wb\boldsymbol W^bWb获得最大的预编码增益。因此，预编码矩阵W\boldsymbol WW可以表示为：W=WaWb\boldsymbol W=\boldsymbol W ^a\boldsymbol W^bW=WaWb

波束成形：波束成形是天线技术与数字信号处理技术的结合，目的用于定向信号传输或接收。波束成形，源于自适应天线的一个概念。接收端的信号处理，可以通过对多天线阵元接收到的各路信号进行加权合成，形成所需的理想信号；同样原理也适用用于发射端。
  如果要采用波束成形技术， 前提是必须采用多天线系统。例如，多入多出(MIMO)，不仅采用多接收天线，还可用多发射天线。由于采用了多组天线，从发射端到接收端无线信号对应同一条空间流(spatial streams)， 是通过多条路径传输的。在接收端采用一定的算法对多个天线收到信号进行处理，就可以明显改善接收端的信噪比。即使在接收端较远时，也能获得较好的信号质量。

秩亏：Rank defect，如果一个矩阵的秩是具有相同大小的矩阵能达到的最高秩，则该矩阵为满秩；如果矩阵不具有满秩，则该矩阵为秩亏。通常引起系数阵秩亏的原因有两种：由于起算数据不足而引起系数阵秩亏，称之为数亏；另一种就是当必要观测数据不足而引起的系数阵秩亏，称为形亏。但是在测量中第二种情况实际上不存在，我们在测量中观测总个数总会大于必要观测个数，因而很少会出现形亏的现象，因此，主要讨论有关数亏的情况。

正则化：Regularization，代数几何中的一个概念。求解不适定问题的普遍方法是：用一组与原不适定问题相"邻近"的适定问题的解去逼近原问题的解,这种方法称为正则化方法。
  不适定问题：适定问题是指满足下列三个要求的问题:①解是存在的(存在性);②解是唯一的(唯一性);③解连续依赖于初始值条件(稳定性)。这三个要求中，只要有一个不满足，则称之为不适定问题。特别，如果条件③不满足，那么就称为阿达马意义下的不适定问题。一般地说不适定问题，常常是指阿达马意义下的不适定问题。

伪逆：Pseudo inverse，伪逆矩阵是逆矩阵的广义形式。由于奇异矩阵或非方阵的矩阵不存在逆矩阵，假设存在一个与AAA的转置矩阵ATA^TAT同型的矩阵XXX，并且满足：AXA=A,XAX=XAXA=A,XAX=XAXA=A,XAX=X，此时，称矩阵XXX为矩阵AAA的伪逆，也称为广义逆矩阵。

[1]金石, 温朝凯等. 大规模MIMO传输理论与关键技术[M]. 北京: 电子工业出版社, 2018
[2]秦舒雅, 杨龙祥. 面向5G的大规模MIMO预编码算法比较研究[J]. 计算机技术与发展, 2015, 25(07): 150-154.
[3]徐悦. Massive MIMO系统中预编码技术的研究[D]. 北京邮电大学, 2017.
[4]Jian Wu,Shu Fang,Lei Li,Yang Yang. QR decomposition and gram Schmidt orthogonalization based low-complexity multi-user MIMO precoding[P]. 10th International Conference on Wireless Communications, Networking and Mobile Computing (WiCOM 2014),2014.

mimo 编码 关键