基于DFT码本的预编码

-.-.可修编-题目：基于DFT码本的预编码学院：通信工程学院专业：电子与通信工程姓名：小瑜学号：1501120442-.-.可修编-基于DFT码本的预编码摘要本文主要研究了在预编码技术中，针对反馈开销十分昂贵，为了节约成本采用基于DFT码本选择的预编码，并对预编码过程进行了性能仿真。仿真结果表明，在QUOTE天线模式下基于DFT码本的预编码简单实用，并且随着不同反馈比特数的增加，基于DFT码本的预编码技术的系统性能随之增加。关键词预编码DFT码本选择引言随着现代移动通信中对吞吐量和频谱效率要求的不断提高，对于基于多天线的多输入多输出的MIMO技术的研究不断深入。而在多用户MIMO系统中，存在的主要问题就是各个用户之间存在干扰。预编码技术的提出则有效的抑制了MIMO信道中的多用户干扰，不仅如此预编码技术在MIMO系统中能显著提升信道容量，并大大降低接收机的复杂度，成为当前的研究的关键技术之一。预编码技术就是在已知信道状态信息（CSI）的情况下,在发射端调整发射策略，接收端进行均衡，从而提高MIMO系统的性能。在多用户MIMO下行链路中，各个用户之间无法相互协作，不能利用上行链路的联合检测来恢复发射信号，因此预编码是多用户MIMO下行链路获得复用增益和分集增益的关键。预编码技术可以分为基于码本的预编码方式和基于非码本的预编码方式。基于码本的预编码方法要求发送端和接收端共享一套码本集合，然后根据具体的信状况从一个确定的矩阵集合中选取一个使系统性能最优的矩阵，再将矩阵在码本集合中的序号反馈给发送端。这样的编码方案是针对反馈信道十分昂贵的开销提出的，并且反馈信道所需传输的数据量很小，只需要几个比特的大小，大大节约了成本。而基于非码本的预编码是在发送端已知信道信息状态的情况下，对信道信息矩阵H进行适当的分解，进而得到相应的发射端预编码矩阵、接收端均衡矩阵和将MIMO信道变换成若干个独立子信道的等效信道矩阵.2.基本原理2.1预编码基本原理预编码技术是在发射端已知或可以获得信道状态信息的情况下，利用信道状态信息对待发射信号进行预处理的技术，从而进一步提高系统的吞吐量。预编码可以表述为：其中，X为发射信号向量，W为预编码矩阵，R为接收向量，H为空间信道矩阵，n为噪声向量，QUOTE为经过预编码之后的等效空间信道矩阵。预编码技术要求将信道划分为N个互不影响的区域，每个区域由M个比特组成的索引值来表示，并且为每个信道区域确定一个合适的发射模式。2.2基于码本选择的预编码预编码技术一般是要通过反馈在接收端获得完整或部分的信道信息状态，每次通过反馈有限个比特的数据信息来传递更多有效的信息数据。为了使反馈资源充分利用，预编码技术采用基于码本与非码本的设方法。目前常见的基于非码本的预编码方式有：基于线性接收机的SVD(奇异值分解)算法、基于ZF-SIC（破零串行干扰消除）接收机的GMD(几何均值分解)算法、改进的UCD(统一信道分解)算法。但是由于反馈信道的开销十分昂贵，本文主要研究的是基于码本的预编码。常用的基于码本设计方式有基于Household变换、基于DFT矩阵以及直接量化等方式。基于码本的预编码就是在接收端和发送端都共享一个已知码本集合，码本集合包含多个预编码矩阵，接收端根据信道估计的信道矩阵以某一性能目标在码本集合中选择使系统性能最优的预编码矩阵，再将码本序号反馈给发送端，发送端根据序号选择预编码矩阵进行预编码。由于反馈信息只需要码本序号，大大减少了反馈量，节约了带宽。预编码预编码数据流信道码本选择信道估计码本序号均衡图2-1基于码本选择的预编码的系统流程框图3.基于DFT阵的码本设计基于DFT矩阵的码本实现起来简单易行，故本文采用的基于码本的预编码选择的是DFT阵。若发送天线数目为M，U为预编码矩阵，um-1为矩阵U的第m个预编码向量。预编码矩阵构造如下：,因为每个预编码向量类似于DFT结构，所以称为DFT码本算法。如果M=2，则预编码矩阵为：如果M=4，则预编码矩阵为：根据上面构造的码本和秩的自适应可以构造出基DFT的码本，其原则是根据秩的大小选取码本中的子集。一般来讲,基于DFT矩阵构造的码本都是以上述矩阵为基础乘以权值进行的。由于上述的DF码本在满秩的情况下，其预编码矩阵只有一个，进一步改进上述的DFT预编码矩阵。若{U(0)…U(G-1)}为预编码矩阵集，其中为第g个预编码矩阵，u(g)m-1为矩阵U(g)的第m个预编码向量。其定义为：其中；其中G的大小决定满秩情况下的码本大小。比如，当发送天线M=2，G=2时，其预编码码本大小为2，码本包括如下两个预编码矩阵，即：，当M=4，G=2时，码本将包含如下两个预编码，即：，4．性能仿真及分析本文主要是针对基于DFT选择码本的预编码过程进行了仿真分析，在这样预编码中，反馈信道所需传输的数据量很小，一般只有几个比特的大小，大大节约了资源。图3-1分别仿真了反馈信道为2比特、4比特和7比特的系统性能，以4为发射天线，2为接收天线，数据经过层映射后分为2层数据流，经过预编码处理后，由4根发射天线发射。接收端天线进行信道估计后进行DFT预编码码本的选择，最后将反馈结果发给发送端。仿真过程是以QPSK调制，接收端采用MMSE均衡算法，以平坦瑞利衰落信道为模型进行的。从图2-1中可以看出随着反馈比特数的增加，基于DFT码本的预编码技术的系统性能也随之增加。图4-1不同反馈比特对基于DFT码本预编码的系统性能的影响5.结束语本文主要针对预编码技术中的基于DFT码本选择的预编码进行了分析研究，并使用MATLAB仿真工具对其预编码过程进行了性能仿真。仿真结果表明，采用基于DFT码本的预编码技术的系统性能随着其反馈比特数的增加而增加，基于DFT码本的预编码过程简单易行，反馈信道所需传输的数据量很小，只需要几个比特的大小。参考文献[1]SAMPATHH,STOICAP,PAULRAJA,Ceneralizedlinearprecoderanddecoderdesignformimochannelsusingtheweightedmmsecriterion.IEEETransCommun,2001,49:2198~2206.[2]SParkvall,AFuruskar,EDahlman.EvolutionofLTEtowardIMT-Advanced[J].IEEECommunicationsMagazine,2011,49(2):84-91.[3]峰.MIMO无线通信系统的预编码技术[D].:交通大学，2009.[4]戴继生.线性MIMO系统预编码技术的研究[D].:中国科学技术大学，2010.[5]TexasInstruments.CodebookDesignforE-UTRAMIMOPre-coding[J].3GPPTSGRANWG146bis,Seoul,Korea,2006.[6]JrRWHeath,SSandhu,APaulraj.Antennaselectionforsoatialmultiplexingsystemswithlinearreceivers[J].CommunicationsLetters,IEEE,2001,5(4):142-144.[7]YuanXJ,GuoQH,WangXD,etal.Evolutionanalysisoflow-costiterativeequalizationincodedlinearsystemswithcyclicpre-fixes.IEEEJournalonSelectedAreasinCommunications,2008,26(2):301~310[8]MJoham,WUtschick,JANossek.LineartransmitprocessinginMIMOcommunicationssystems[J].SignalProcessing,IEEETransactionson,2005,53(8):2700-2712.附录：仿真代码如下：clc;clearall;GG=128;%产生128个DFT码本form=1:4forn=1:2forg=1:GGW3(m,n,g)=exp(i*2*pi*(m-1)*((n-1)+(g-1)/GG)/4)/2;endendendW3load;L=1000;X=rand(1,L);fori=1:LifX(i)<0.5X(i)=0;elseX(i)=1;endend%随机生成1000个信源序列x=QPSK(X);%QPSK调制s=layer(x,2);%信号分层forG=1:100%循环100次GH=(randn(2,4)+j*randn(2,4))/sqrt(2);%随机生成复高斯分布的信道特性矩阵%DFT码本反馈2bitm=1;forSN=-6:2:14error=0;N=1/10^(SN/10);n=sqrt(N)*(randn(2,L/4)+j*randn(2,L/4))/sqrt(2);fori=1:4M=H*W3(:,:,i);[U1,S1,V1]=svd(M);D1=diag(S1);d1(i)=D1(2);end[t1,A]=max(d1);F=W3(:,:,A);y1=H*F*s+n;K=H*F;y2=(K'*K+N*eye(2))\K'*y1;y=reshape(y2,1,L/2);Y=QPSKJT(y);fori=1:LifY(i)~=X(i)error=error+1;endendPe_svd4(G,m)=error/L;m=m+1;SNend%DFT码本反馈4bitm=1;forSN=-6:2:14error=0;N=1/10^(SN/10);n=sqrt(N)*(randn(2,L/4)+j*randn(2,L/4))/sqrt(2);fori=1:16M=H*W3(:,:,i);[U2,S2,V2]=svd(M);D2=diag(S2);d2(i)=D2(2);end[t2,B]=max(d2);F=W3(:,:,B);y1=H*F*s+n;K=H*F;y2=(K'*K+N*eye(2))\K'*y1;y=reshape(y2,1,L/2);Y=QPSKJT(y);fori=1:LifY(i)~=X(i)error=error+1;endendPe_svd5(G,m)=error/L;m=m+1;SNend%DFT码本反馈7比特m=1;forSN=-6:2:14error=0;N=1/10^(SN/10);n=sqrt(N)*(randn(2,L/4)+j*randn(2,L/4))/sqrt(2);fori=1:128M=H*W3(:,:,i);[U3,S3,V3]=svd(M);D3=diag(S3);d3(i)=D3(2);end[t3,C]=max(d3);