基于matlab的MIMO通信系统

1、北京邮电大学实验报告基于Matlab的MIMOB信系统仿真目录概述错误!课程设计目标：错误!数字通信系统概述错误!二、基本原理错误!三、仿真设计错误!流程图错误!主要模块错误!主要参数错误!错误!结构性和关键语句错误!状态检验错误!性能测试错误!五、误码性能错误!六、仿真结果与分析错误!七、重点研究的问题错误!八、结论错误!九、参考文献错误!一、概述课程设计目标：（1）了解移动通信关键技术（2）了解数字通信系统仿真流程（3）学会用Matlab实现基本的信道编译码、调制解调等通信模块（4）学习并实现MIMO空时处理技术（5）学习和掌握性能分析的思路和方法。数字通信系统概述信道中传输的是数字信号的

2、通信方式称为数字通信，它包括将基带数字信号直接送往信道传输的数字基带传输和经载波调制后在送往信道传输的数字载波传输。对应的通信系统称为数字通信系统。具体的通信流程如上图所示。对于本次实验中需要完成的MIMO通信系统，流程相应为：产生信源，而后利用卷积码进行信道编码，调制方式使用QPSKS制，再将其通过双天线传输后通过白高斯噪声信道，传播后利用ML进行判决，译码方式采用维特比译码，最终得到传输后的数据。二、基本原理概述及原理MIMO表示多输入多输出。在第四代移动通信技术标准中被广泛采用，例如IEEE（Wimax）,长期演进（LTE）在新一代无线局域网（WLAN肺准中，通常用于IEEE,但也可以用

3、于其他技术。MIMO有时被称作空间分集，因为它使用多空间通道传送和接收数据。只有站点（移动设备）或接入点（AP）支持MIMO时才能部署MIMO。MIMO技术的应用，使空间成为一种可以用于提高性能的资源，并能够增加无线系统的覆盖范围。MIMO允许多个天线同时发送和接收多个空间流，并能够区分发往或来自不同空间方位的信号。多天线系统的应用，使得多达min(Nt,Nr)的并行数据流可以同时传送。同时，在发送端或接收端采用多天线，可以显著克服信道的衰落，降低误码率。利用MIMO技术可以提高信道的容量，同时也可以提高信道的可靠性，降低误码率。前者是利用MIMO信道提供的空间复用增益，后者是利用MIMO信道

4、提供的空间分集增益。实现空间复用增益的算法主要有贝尔实验室的BLAST算法、ZF算法、MMSE算法、ML算法。目前MIMO技术领域另一个研究热点就是空时编码。常见的空时码有空时块码、空时格码。空时码的主要思想是利用空间和时间上的编码实现一定的空间分集和时间分集，从而降低信道误码率。信道模型nT发送天线，nR接收天线信道矩阵:hiihi2L几nTh2ih22Lh2nTMMOM接收信号：yHxn即接收信号为信道衰落系数发射信号+接收端噪声MIMO模型中有一个空时编码器，同时含有有多根天线，其系统模型和上述MIMO系统理论一致。一般来说，移动终端所支持的天线数目总是比基站端要少所以发送天线的数目要大

5、于接收天线。分集与复用根据各根天线上发送信息的差别，MIMO可以分为发射分集技术和空间复用技术。发射分集：在不同的天线上发射包含同样信息的信号（信号的具体形式不一定完全相同），达到空间分集的效果，起到抗衰落的作用典型代表：空时块码（STBC）空间复用：在不同的天线上发射不同的信息，获得空间复用增益，从而大大提高系统的容量和频谱利用率典型代表：分层空时码空时块码（STB。Alamouti提出了采用两个发射天线和一个接收天线的系统可以得到采用一个发射天线两个接收天线系统同样的分集增益。Alamouti发送分集空时编码方案将每k个输入字符映射为一个nTp矩阵，矩阵的每行对应在p个不同的时间问隔里不同

6、天线上所发送的符号。2*2举例:XiX2X2Xi单天线接收，用r1, r2表示第1,2个发射符号间隔接收天线的接收信号:r1'X1h2x2 n1r2h1x2h2x1n2采用最优极大似然译码可得:arg(minr-H?2)为了要满足各根天线上发送数据的正交，它的编码矩阵需要满足如下条件HXnTXnT22x1x2L|Xk 2 InTx1 argminx2argminhj,1hj,2hj,1hj,2j 121 x2d22dx2, x2多天线接收，发送端的编码与传输方案和单接收天线系统一样。只是在接收端的处理变得复杂，需要对不同接收天线上接收到的信号进行合并处理。多接收天线下的判决度量可以通过

7、把各副接收天线上的接受信号得到的判决度量线性合并得到。判决式如下：分层空时码（BLAST空时编码利用多天线组成的天线阵同时发送和接收。在发送端，将数据流分离成多个支流，对每个支流进行空时处理和信号设计（空时编码），然后通过不同天线同时发送；在接收端，利用天线阵接收，并经过空时处理和空时码解码，还原成发送数据流。天线1串并变换信道编码器1信道编码器2空时编码调制器1天线n调制器2图1分层空时码的发射端系统模型在接收端，用多个天线分集接收，信道参数通过信道估计获得，由线性判决反馈均衡器实现分层判决反馈干扰抵消，然后进行分层空时译码，单个信道译码器完成信道译码，分层空时码接收端系统框图如下图所示：天

8、线1天线2|线性判决反馈均衡空时译信道译三、仿真设计流程图信道编码采用卷积码进行QPSK调制STBC双天线发送通过高斯白噪声信道单/多天线接收Ml准则进行判决维特比译码进行解调结束主要模块信源产生要求：产生独立等概的二进制信源实现：matlab自带函数randint()一个信源产生N/2个0,1随机等概数据S1=randint(1,N/2);S2=randint(1,N/2);信道编码方式：卷积码(3GPP )adHHdmdVHd_Dh-r*LPT卷积纠错编码函数convenc()Output QGn=557gHl)Output 1Gt =63 (octal)h_Output2G2=711(o

9、ctal)格式：code=convenc(msg,trellis)功能：利用poly2trellis函数定义的格形trellis结构，对二进制矢量信息msg进行卷积编码。将卷积编码多项式转换成格形(trellis)结构函数poly2trellis()格式：trellis=poly2trellis(constrainlength,codegenerator)功能：将前向反馈卷积编码器的多项式转换成一格形(trellis)结构。Viterbi算法译卷积码函数vitdec()格式：decoded=vitdec(code,trellis,tblen,opmode,dectype)功能：利用Viterb

10、i算法译卷积码。Code为poly2trellis函数或istrellis函数定义的格形trellis结构的卷积码。参数tblen取正整数，表示记忆(traceback麻度。参数opmode代表解码操作模型。调制与解调基本要求QPSK戳破根址为。多工芋的L2K星座，冷栽油相位为于，芋寿和学的qp$k星座在数字信号的调制方式，QPSK四相移键控是目前最常用的一种卫星数字信号调制方式，它具有较高的频谱利用率、较强的抗干扰性、在电路上实现也较为简单。QPSK是在M=4时的调相技术，它规定了四种载波相位，分别为45°,135°,225°,315°,调制器输入的数

11、据是二进制数字序列，为了能和四进制的载波相位配合起来，则需要把二进制数据变换为四进制数据，这就是说需要把二进制数字序列中每两个比特分成一组，共有四种组合，即00,01,10,11,其中每一组称为双比特码元。每一个双比特码元是由两位二进制信息比特组成，它们分别代表四进制四个符号中的一个符号。QPSK中每次调制可传输2个信息比特，这些信息比特是通过载波的四种相位来传递的。解调器根据星座图及接收到的载波信号的相位来判断发送端发送的信息比特。实现方法:调制：Matlab自带函数reshape以及自定义star()利用星座点调制解调：Matlab自带函数reshape以及ML思想解调，过程中使用min(

12、)函数提高要求16QAM:t农16QAM是用两路独立的正交4ASK信号叠加而成，4ASK是用多电平信号去键控载波而得到的信号。它是2ASK体制的推广，和2ASK相比，这种体制的优点在于信息传输速率高。用/并变换器将速率为Rb的二进制码元序列分为两路，速率为Rb/电平变换为Rb/2的二进制码元序列变成4个电平信号,4电平信号与正交载波相乘，完成正交调制，两路信号叠加后产生16QAM信号.在两路速率为Rb/2的二进制码元序列中，经2-4电平变换器输出为4电平信号，即M=16.经4电平正交幅度调制和叠加后，输出16个信号状态，即16QAM.实现方法：调制：Matlab自带函数reshape以及sta

13、r()用星座点进行调制解调：Matlab自带函数reshape以及ML思想解调，过程中使用min()函数空时分组码发送与接收两发一收:r；二ha工十h二匕十n：r；=-h；Xi+h：x3十址(2)两发两收：实现：Matlab自带conj()与.*信道加性高斯白噪声AWGN(AdditiveWhiteGaussianNoise)是最基本的噪声与干扰模型。加性噪声：叠加在信号上的一种噪声，通常记为n(t),而且无论有无信号，噪声n都是始终存在的。因此通常称它为加性噪声或者加性干扰。白噪声：噪声的功率谱密度在所有的频率上均为一常数，则称这样的噪声为白噪声。如果白噪声取值的概率分布服从高斯分布，则称这

14、样的噪声为高斯白噪声。实现:Matlab自带函数randn()主要参数% 数据包个数%包长PACKET=10;N=2560*4;SNR_dB=0:2:30;%信噪比范围(030),图中横坐标间隔为 2tblen=30;%卷积码记忆深度，寄存器的35倍最佳，选择30四、程序块设计结构性和关键语句信源产生ppt中建议使用randsrc()函数，1由于randsrc()产生的信源序列是+1和-1的序列，而进行卷积码编码时输入编码器的序列必须是单极性不归零序列，所以要将序列中的-1修改为0.于是我采用了randi()函数，直接可以产生随机且等概的0,1序列。S1=randi(1,N/2);%产生第1组

15、信号，长度为N/2S2=randi(1,N/2);%产生第2组信号，长度为N/2信道编码卷积码%8进制表示寄存器状态%转换为格行结构%产生随机信号1%产生随机信号2%对 1进行卷积编码%对 2进行卷积编码codegen=557663711;trellis=poly2trellis(constlen,codegen);forbag=1:PACKETS1=randi(1,N/2);S2=randi(1,N/2);S=S1S2;SJ_1=convenc(S1,trellis);SJ_2=convenc(S2,trellis);根据图可得，有9个移位寄存器，1个输入和3个输出。将编码器的3个冲激响应写

16、成二进制序列，最后转化成八进制，最后得到的结果为557,663,711,为八进制的表示形式。维特比译码直接利用公式vitdec()即可decode.l=vitdeg(result_ljtKellisjtblen,'tmndJhard');*2度2收QPSK信号1decode_2=vitdec(Te5uLt_2j,trellT-SjlblenjJtrime'','hatd''】!&2lft!QP£k信号2decode_3=vitdec(result-,IrelliSjtblen,'truncJbar)：*?直llf

17、eQP5K信号1decode_4=vitdec(Tesult_trellis,ttilen.,trunc''hard")涓口发1收QPSR信号2decode_5=vitdec(result.Sjtrellis,tblen,'trune4FJhartf】:第2窟2收1.60/僖W1decQd_6=vitdec(xesult_6jtrellis,tblen,?tr-unc'',；hard')”二发二版15QAM信号二deeode_7=vitdec(result_7rtrellis,tblen,'trune/Jbard*)：K2!SL

18、收L6QAM信号1decode_3-vitdec(iesult_8jtrellis,tblen,Jtrnri,c''；hard"):雷2发】收1gAM信等2冗余校验能被已知g(x)®除即判断为无错，否则产生错误选取g(x)=1+xA5+xA12+xA16加冗余过程：上-sire-input,2):g«n=reros(ljx+i):output二(bk+a);if(output(1)二二1)outputCI:x+11=inodt(output(1ixH-D-bgen)j2):end检查过程forii1:n-xif(input(1)=1)input(l

19、;x+1)=mgd(finput(1;h+L)+5en)?21;endirpvt=input(2:end)input(l);endifsujn''input)=Qindicate-0：elseindicate=1;调制与解调QPSK星座图与调制star_qpsk=ezip(j*-3/4+pi3/4*pi-1/4*pil/4*pi);JfiqPSK星由图(格雷科顺序)%star=eip(j*-3/4*pi3/4*pi-l/4*pi)：A格雷码顺序bits_reshape(bits_in,Z,口)：黑串并变换bits_0_3=2l*bits_2:常Q】序列变k1,&1序列j

20、ipduL_seq=star(bi-ts_0_3-l);需调制点输出QPSK解调fori=l:4iusti,;)=(abs(seQ-star(i),*2;，肛思想，遍历所有星座点选取距离最小end,1)；£0存械最小值的行数(1T)一网存的最小的数字cle3rdis;demothil_seq=a(10,:):16QAM星座图i=l;Ifork=-3:2:3%k为实部,H为虚部forjh=-3:2:3start(i)=(k+j*)/sqrt(10);珀功军归一匕i=i+l;end.：end;star_16qamFstart(013245761213151439111OJ+1)克160A

21、M星座典(格雷码I顺序)16QAM调制bitreshape(bLts_Ln54,)：句串并变换bits_0_l&=842l*bit_4;WT序列费Q-15序列medul_seq=starbits_0_15+1)；与调刮点输出16QAM解调fori=l:16曜醴二）=（whsIwmq-star（i）.,2;里ML思熊，遍历断背星座点选期题离最小endyljilj=TRiTiCdis?jj1）；学工。存放最小值的行数（1-4）下。存放最小的数字cleardis2:deniodul_sq=b1ib：J;瑞利信道信道中的信号hQ-sqrttQ.5L*(rantin.11,lenn-j.tian

22、dnChlen>):%瑞利信道参数hl=sqrt<0,5)+(randnlj4-j.*randn(ljlen)：h2=sqr-t(0.5).*(Tandn.'19len.)+ai.*TajidnCljlen):h3=sqit(0-5)-*(randnCljlen)*i,*iandn(ljlen)：根据瑞利信道衰减参数，加之公式可得两根天线在两个符号周期内在信道中的总信号rij（包括高斯噪声）%蔚一个符号周期KU=hO,*T_seql+hl.*T_seq2;1天线1R21=h2.*I_seql+h3.*I_seq2;看天线E%第二个符号周期K12=-hO.*conj(T_sq

23、2)+hk*conjG_s«ql)：¥天然1K22=-h2.*conj(T_geq2)+h3»*conj(T_scql);篇天战2为描敢信号加入高斯喔声黑第一个符号周期搔收总信旨r11=R11+sgnta.*2*1ranjdn-(13len)+i*t3ndn(.llen);r21=(R21+sgnia.*2*(Esndn(Llez+itrsndndjlen),)：，第二个符母周期撞收总信号rl'(RI2+senta,*"1214SB7B410UU 1T 0100 j 11G1 1T12J4567B910111 卜”El - &|0L707

24、卜UMWl + 0U 心丽十TJEH 中0JW1 十:1»硒 + k|a7OTl + X3OT1 04二冽中Q7OTLOugM -经过16QAM后生成每路信号3840个数据因为16QAM将每路信号分成四路后调制即3840=153604 个数据randn.Ujlen)+*rwdn(ljIcnl):r22-(R22+sgjna.*2*(randn.(ljlen)+i*randnfljlen):误比特率的计算利用matlab自带biterr()函数Biterr()用于比较两组标量序列中不同位数的个数用个数与总位数相除记得到本次循环中的误比特率nO,eO=biterrld«eod&

25、amp;_qpskjS);见本次循环的课:比特率temp.EEO1皿)=nO/W;因为根据不同的系统可能需要多次循环，故而要对每次循环中的误比特率进行求和并与总循环次数相除，等到真正的误比特率。5E0(bagjjj)=Iswn(t2)/loop;用整体循环的限比恃率状态检验每路信源产生N/2=5120个数据田SI<lriM2OdtMbla*经过卷积编码后生成15360个数据因为有三路输出即15360=5120*3个数据SJ.ldoubles234567gID11111011D1一QPSK1座图（即坐标数据）士5tsr_lcomplexdouble>12341回.7071仇-0.70

26、71¥.0.7071国0.0.7071+0.16QAM星座图数据（格雷码顺序排列）Ii|%Lar_2<IxlbKempledeubrL2J4567&91G111-0_-O.=D.94a?-J-C.94B?+_-O_=_31&2-0.-031&2t.-0312-一D.S4B?-口a.MH?0l.口.皿7，£l_EM3I''I-I1-I经过QPSKt生成每路信号7680个数据因为QPSK等每路信号分成两路后调制即7680=15360/2个数据SdpflkjmH7M0CTfflipIndoubh?>12145678gID111回

27、工皿To.03162031醛一鼠031£2+k9W-CL.-DlMBJ*“也数一.444S7-限-0-M«7*a4一-,&«?-J。.帆2经过两发两收的QPSK言号调制后产生的信号每路7680个数据与发送前相同h|£一1七Lt：七即CQE理型4LlCMjbk、工234567a91a111?.2742",1JJ21,L-。跳"417M巾+-2.9640+.3川/*1M51+2-3熨仙*3r395ST-k.2J776*25611-2L60：卜u两发两收未编码情况每路产生2560个数据djL*wmple<dcul>e&g

28、t;123i4S67BiDll1|口湖-0£37第-23-339.丐+280舫+Z.-2.1J074+3535gl+d，L459S+LZ941K)-3-2.48St-2.I.&2SJ+L.0713-3上怔经过两发一收的QPSK言号调制后产生的信号每路7680个数据H|5_1-Ix7bSOcomplexcinubJeL234567a910n113血时.i口wot.uj<BM+14=urar*-ijk7*.xaaw+二117»*Jjaw*j.nw*Ji4aoe经过两发两收的16QAM信号调制后产生的信号每路3840个数据与发送前相同士755,1*1x38+5com

29、plexdcvjblr>12345S78gIQ111loim一:“10194-J-aJ6M-0.1.1565-+L.-19727-Cl.-13时?+-5Jai5-2-1,M54+一-7.4178-6Mm-*经过两发一收的16QAM信号调制后产生的信号每路3840个数据士拶屿138*0ecmpl?*dotib>12土4567sg1011|10?倒-k.(UBM-flJQ2双4叫2447(142.9m*-W1CMr鼻丽d4irr.Mlk-CL.2,M702L州E2解调2发2收QPS所产生序列15360个经过维特比译码后2发2收产生序列5120个已恢复为0,1序列性能测试当采用packe

30、t=10时，每包长2460*4共98400个数据耗时：34s较为快捷，可见代码性能较优此时加入冗余校验，速度大大下降耗时：8'45'约是不添加冗余校验的17倍且由于在此系统中冗余校验无法实现其功能，即实际纠错故而认为冗余校验用处不大增加包长，不加入冗余校验packet=50,每包长为依然为2460*4共492000个数据耗时2'44''且两幅图得到的误比特率图线无明显差别，相较之下选择包长为10进行仿真.五、误码性能包长为packet=10时，数据共有98400个此时得到的曲线与packet=50数据增加5倍即492000个数据时的在QPSKF误码率图线

31、无明显差异，细微来看，packet=50稍平滑一些=10,tblen=30二50当包长都取10时，改变卷积码的记忆深度。查取资料可知，卷积码的记忆深度为寄存器个数的35倍时性能最佳。本次仿真中因为寄存器个数为8个，所以记忆深度在2440中选取比较合适。于是我选取30作为常态。比较大的记忆深度50和比较小的记忆深度3来画图比较，从中可以看出，记忆深度较小（=3）时，未编码曲线较为光滑记忆深度增大（=30,=50）,编码更为曲线光滑packet=10tblen=501010一2发2收卷积编码图线-2发2收未编码图线2发1收卷积编码图线1001015SNRdB202530packet=10tblen

32、=310°10"10"2发2收卷积编码图线2发2收未编码图线2发1收卷积编码图线1(/10J1015SNRdB202530六、仿真结果与分析情况下误码率图线我选取了三个状态的误比特率来进行分析分别是2发2收的卷积编码BER图线2发1收的卷积编码BER图线以及未进行卷积编码直接调制的BER图线。首先我们已知2发2收的性能优于2发1收，观察图线，2发2收图线一直在2发1收图线之下，即相同信噪比时，2发2收的误比特率小于2发1收，正确。对于未进行卷积编码的图线与进行编码的图线相比，在信噪比较小的时候（约小于10db）未编码信号的误比特率比2发1收的进行QPSK的信号误比

33、特率低，因为本次仿真中的误比特率是经过逐次取平均，扩大范围继续平均，数据准确，于是可以得到信噪比小的时候，收发方式对误比特率影响较大。而当超过10db时，未进行卷积码的信号误比特率远远高于经过编码的信号，可见在信噪比大时，编码对通信系统性能的影响更大。可见在信噪比约为15dB以下时，两条曲线几乎交叠在一起，误比特率相当,然而当信噪比大于15dB时，2发2收误比特率大于2发1收的情况。不符合常理，有待解决。七、重点研究的问题、QPSKt16QAM的调制最开始也考虑过使用matlab直接提供的函数pskmod(),qammod()和相应的解调函数来进行QPSKM16QAM的调制与解调。但因为学习过

34、通信原理，理解了调制的本质实际就是映射，将一维标量映射到QPSK的单位圆及16QAM的星座点上。所以我采用了自己设计调制与解调的函数其实实施起来也很简单。首先定义QPSKS16QAM的星座点(6载波相位为手，斗,孥和千的。P5K星座我采用了这种在通信原理中因为错到相邻符号概率小而性能较好的方式，只用直接定义4个点的位置就完成了。16QAM也很简单，同样采用通信原理中较为推荐的格雷码方式排序星座点，用循环的方式定义星座点，以期获得较低的误符号率。之后进行串并变化后，把数据分别变为2行和4行，一列看为一个符号，转换为10进制数字后+1得到在星座点中的序列位置，就完成了调制。、QPSKt16QAM的

35、解调完成了调制，解调的思想也比较简单。最基础的先要应用ML准则，即找到离当前位置最近的星座点，记录它的位置，再回到回复序列中找到相应序号对应的2或4位数据，再次使用reshape即可恢复为原信号。、瑞利信道中的发送与接收弹性要求中有瑞利信道这一项，所以我查阅了相关资料后了解到，只要定义好锐利信道的参数h0h3后就能利用书中的公式，将两根天线上相邻时刻到达的信号与瑞利信道系数想乘，就可以得到天线的发射信号，再加上随机噪声sgma.*2*(randn(1,len)+j*randn(1,len)，便得到了信道中的总函数。译码部分比起发送要复杂一些。发送的时候有相应的公式，直接就可以得到，而对于译码，

36、则需要理解两发两收和一收进行的含义，才能写出相应的算法。两发两收交两发一收复杂，但在理解了具体实现的基础上，需要我们用已有的简单的数学知识(如负数的共轭计算，算式相消技巧)，最终将收端得到的信号恢复为原信号。、冗余校验在本学期的信息处理与编码与编码课中，作为课余知识简单的了解了CRC冗余校验的实现。即与给定的g(x)相除，整除可认为无差错。根据一篇论文尝试编写了冗余校验的程序，虽然程序实现并不太难，但却给程序造成了极大的冗余。在调制前要给两路信号添加冗余，解调译码后要对相应的8个结果解校验。于是程序的运行时间平均提高了近20倍。在之后的调试代码过程中一直将冗余校验注释掉，以免浪费过多的时间。、误比特率的计算比起通信原理中误比特率的计算需要很长的公式，开根乘方不同，在matlabbiterr（