根据Matlab的MIMO通信系统仿真

《根据Matlab的MIMO通信系统仿真》由会员分享，可在线阅读，更多相关《根据Matlab的MIMO通信系统仿真（14页珍藏版）》请在装配图网上搜索。

1、北京邮电大学基于M atlab的MIMO通信系统仿真专业：信息工程 班级：2011211126 姓名：学号：目录一、概述11、课题的研究背景12、课程设计的研究目的13、MIMO 系统1【1】MIMO 的三种主要技术1【2】MIMO 系统的概述2【3】MIMO 系统的信道模型2二、基本原理31、基本流程32、MIMO 原理33、空时块码4三、仿真设计51、流程图52、主要模块及参数53、信源产生54、信道编码65、调制66、AWGN 信道67、输出统计7四、程序块设计71、代码7五、仿真结果分析111、仿真图112、结果分析12六、重点研究的问题12七、心得与体会12八、参考文献12一、概述1

2、、背景MIMO表示多输入多输出。在第四代移动通信技术标准中被广泛采用，例如 IEEE 802.16e (Wimax),长期演进(LIE)。在新一代无线局域网(WLAN )标准中,通常用于 IEEE 802.11n,但也可以用于其他802.11技术。MIMO有时被称作空间分集，因为它使用 多空间通道传送和接收数据。只有站点(移动设备)或接入点AP)支持MIMO时才能部署 MIMO。MIMO技术可以显著克服信道的衰落，降低误码率。该技术的应用，使空间成为一种可 以用于提高性能的资源，并能够增加无线系统的覆盖范围。通常，多径要引起衰落，因而被视为有害因素。然而研究结果表明，对于MIMO系统来 说，多

3、径可以作为一个有利因素加以利用。MIMO系统在发射端和接收端均采用多天线(或阵 列天线)和多通道，MIMO的多入多出是针对多径无线信道来说的。传输信息流s(k)经过空时 编码形成N个信息子流ci(k)，I=1，N。这N个子流由N个天线发射出去，经空间信 道后由M个接收天线接收。多天线接收机利用先进的空时编码处理能够分开并解码这些数据 子流，从而实现最佳的处理。2、课程设计的研究目的设计一个基于Mat lab的MIMO通信系统仿真。3、MIMO系统【1】MIMO的三种主要技术当前，MIMO技术主要通过3种方式来提升无线传输速率及品质：(1) 空间复用(Spatial Multiplexing):

4、系统将数据分割成多份，分别在发射端的多根天线上发射出去，接收端接收到多个数据 的混合信号后，利用不同空间信道间独立的衰落特性，区分出这些并行的数据流。从而达到 在相同的频率资源内获取更高数据速率的目的。(2) 传输分集技术，以空时编码(Space Time Coding)为代表： 在发射端对数据流进行联合编码以减小由于信道衰落和噪声所导致的符号错误率。空时编码通过在发射端增加信号的冗余度，使信号在接收端获得分集增益。(3) 波束成型(Beam Forming)： 系统通过多根天线产生一个具有指向性的波束，将信号能量集中在欲传输的方向，从而提升信号质量，并减少对其他用户的干扰。【2】MIMO系统

5、的概述多输入多输出(Multi-inputMulti-output； MIMO)是一种用来描述多天线无线通信系 统的抽象数学模型，能利用发射端的多个天线各自独立发送信号，同时在接收端用多个天线 接收并恢复原信息。该技术最早是由马可尼于1908年提出的，他利用多天线来抑制信道衰 落(fading)。根据收发两端天线数量，相对于普通的单输入单输出系统(SingleTnpu t Single-Output, SISO)，MIMO此类多天线技术尚包含早期所谓的“智能型天线”，亦即单 输入多输出系统(Single-Input Multi-Output, SIMO)和多输入单输出系统(Multiple-I

6、nput Single-Output, MISO)。由于MIMO可以在不需要增加带宽或总发送功率耗损(t ransmi t power expend it ure) 的情况下大幅地增加系统的数据吞吐量(throughput)及发送距离，使得此技术于近几年受 到许多瞩目MIMO的核心概念为利用多根发射天线与多根接收天线所提供之空间自由度来 有效提升无线通信系统之频谱效率，以提升传输速率并改善通信质量。【3】MIMO系统的信道模型采用多发送多接收的天线系统，n副发送天线，n副接收天线的MIMO的系统信道模 tr型如图1.3所示。用h表示第i根发送天线到第j根接收天线之间的信道衰落系数，x表i ,

7、j i示i副发送天线发送的信号，y为第j副接收天线接收到的信号，考虑平坦的衰落信道则j它们之间关系写成矩阵的形式为：y(t)二Hx(t) +n (t),其中定义x(t) = (x (t), x (t), A A , x (t )T表示n发射天线在发射的信号向量;12nk接收天线在时间t接收的信号向量;y(t) = (y (t), y (t),A A , y (t)T 表示 n12nkn (t)=们(t),n (t),aa，耳(t)t表t时刻的加性高斯白噪声向量;-h1,1hh1,2hAAH=2,12,2MMAhhAL 厂1nj,212n分布。knj,nkh，H为信道衰落系数矩阵。其中h的幅度服

8、从Rayleighj,k二、基本原理1、基本流程2、MIMO 原理 核心思想：时间上空时信号处理同空间上分集结合。 时间上空时通过在发送端采用空时码实现： 空时分组、空时格码，分层空时码。 空间上分集通过增加空间上天线分布实现。此举可以把原来对用户来说是有害的无线电波多 径传播转变为对用户有利。MIMO 系统模型：ElVhh.22hV空时解码空时编码MIMO 信道模型弓发送天线，叫接收天线 信道矩阵：接收信号：y = Hx + n分集与复用 根据各根天线上发送信息的差别，MIMO可以分为发射分集技术和空间复用技术。 发射分集:在不同的天线上发射包含同样信息的信号（信号的具体形式不一定完全相同）

9、, 达到空间分集的效果，起到抗衰落的作用典型代表：空时块码（STBC） 空间复用:在不同的天线上发射不同的信息，获得空间复用增益，从而大大提高系统的 容量和频谱利用率典型代表：分层空时码3、空时块码 Alamouti 提出了采用两个发射天线和一个接收天线的系统可以得到采用一个发射天线 两个接收天线系统同样的分集增益。 将每k个输入字符映射为一个矩阵，矩阵的每行对应在p个不同的时间间隔里不同天线 上所发送的符号。 2*2 举例：单天线接收，用G与表示第2个发射符号间隔接收天 线的接收信号： z .人=Vl+Vl+Lh二沟兀;+息坊+电多天线接收？采用最优极大似然译码可得：x=arg(irim|r

10、-HA|3为了要满足各根天线上发送数据的丘交，它的编码矩阵 需要满足如下条件，% W=W十|工十L十圍)1弓三、仿真设计1、流程图2、主要模块及参数3、信源产生 产生独立等概二进制信源 Mat lab 函数：randsrc() 编码块长度：40-6114(3GPP TS36.212 Table 5.1.3-3) 推荐值：320,1024,2560,5120利用randsrc ()随机产生0与1的数，利用公式将其变成1与-1，即产生信源。4、信道编码Option 1:卷积码(3GPP TS25.212 4.2.3.1)Output QGo = 557 (ottal)y i” Output 2Gj

11、 -711 (Octal)根据要求，卷积码编码器为(3, 1,8)编码器，G0= (001101111 )转换为八进制为557； G1= (010110011)转换为八进制为663； G1= (011001001)转换为八进制为711.卷积码记忆长 度为9.5、调制QPSK调制的函数可以利用matlab函数库自带的QPSK函数进行编写。但是由于QPSK是M=4 进行映射的，所以将随机产生的信源利用reshape函数变成两列，即00映射成，01映射成， 10映射成，11映射成。解调时利用matlab自带的demodulate函数。6、AWGN信道Matlab 实现：randn()信噪比与AWGN

12、噪声的关系： SNRa = y/2RM lO-由于信源经过QPSK调制，相当于在两条坐标轴的平面（一条实数轴，一条虚数轴）每 一个方向都有随机的高斯白噪声分布，然后各自乘以相应的幅度，将得出的噪声加入QPSK 调制的信号里，得到的便是经过AWGN信道的信号。7、输出统计 误码率统计 输出：横坐标-SNR Eb/NO （Eb/NO vs Es/NO?）纵坐标-BER （误比特率），BLER （误块率） 画图：semilogy（x,y）四、程序块设计1、代码1 、MIMO:SNR_dB=O:1:3O; %信噪比范围 LOOP=5OO;%信源产生%for i=1:length(SNR_dB)snr=

13、10(SNR_dB(i)/10);%信噪比转换sgma=sqrt(8/snr/3);%噪声方差errorO=O;error1=O;error2=O;error3=O;for a=1:LOOPN=256O;%仿真序列长度S=randsrc(l,N,O,l); %二进制序列信源产生n个随机数%有信道编码%SOO=convolutionalcode(S);%QPSK调制s=QPSK(SOO);s0=zeros(1,3*N/4);%将调制后的信息奇偶拆分成两个信息以进行MIMO发射s1=zeros(1,3*N/4);for b1=1:2:(3*N/2) sO(ceil(b1/2)=s(b1);s1(c

14、eil(b1/2)=s(b1+1);endN1=size(s0,2);%有信道编码(2发2收)h0=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1); %设置信道冲激响应 h1=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1);h2=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1); h3=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1);%接收到的信号(冲激响应与发送序列卷积)-R0=h0.*s0+h1.*s1;%接收信号矢量(不含高斯噪声)R1二-conj(sl).*hO

15、+hl.*conj(sO);%移动通信 P166R2=h2.*s0+h3.*s1;R3=-conj(s1).*h2+h3.*conj(s0);%通过AWGN信道加入噪声后的冲击序列 rO=(RO+sgma.*(randn(l,Nl)+li*randn(l,Nl);% 接收信号矢量(有噪声) r1=(R1+sgma.*(randn(1,N1)+1i*randn(1,N1);r2=(R2+sgma.*(randn(1,N1)+1i*randn(1,N1); r3=(R3+sgma.*(randn(1,N1)+1i*randn(1,N1);s_0二conj(hO).*rO+hl.*conj(rl)+

16、conj(h2).*r2+h3.*conj(r3);%STBC 空时译码 s_1=conj(h1).*r0-h0.*conj(r1)+conj(h3).*r2-h2.*conj(r3);for l0=1:3*N/4 %将接收到的信号合并Y0(2*l0-1)=s_0(l0);Y0(2*l0)=s_1(l0);endyO=DEQPSK(YO); %QPSK 解调m0=viterbi(y0) ; %维特比译码nO,eO=bi terr(mO,S);% 计算误比特率error0=error0+n0;%有信道编码2发1收hr0=sqrt (0. 5) . *(randn(1, N1) +1i. *ran

17、dn(1, N1) ) ; %瑞利信道函数 hr1=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1);Rr0=hr0.*s0+hr1.*s1;%接收信号矢量(不含高斯噪声)Rr1=-conj(s1).*hr0+hr1.*conj(s0);rNO=(RrO+sgma.*(randn(l,Nl)+li*randn(l,Nl);% 接收信号矢量(含复高斯噪声) rN1=(Rr1+sgma.*(randn(1,N1)+1i*randn(1,N1);sr_O二conj(hrO).*rNO+hrl.*conj(rNl);%STBC 空时译码 sr_1=conj(hr1).*r

18、N0-hr0.*conj(rN1);for l1=1:3*N/4 %将接收到的信号合并 Y1(2*l1-1)=sr_0(l1);Y1(2*l1)=sr_1(l1);endyl=DEQPSK(Yl);%QPSK 解调m1=viterbi(y1); %维特比译码nl,el=bi terr(ml,S);% 计算误比特率 error1=error1+n1;%无信道编码%sn=QPSK(S); %QPSK调制for b2=l:2:(N/2)%将调制后的信息奇偶拆分成两个信息以进行MIMO发射sn0(ceil(b2/2)=sn(b2);sn1(ceil(b2/2)=sn(b2+1);endN2=size(

19、sn0,2);%无信道编码(2发2收)hn0=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2); %瑞利信道函数 hn1=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2);hn2=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2); hn3=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2);Rn0=hn0.*sn0+hn1.*sn1;%接收信号矢量(无噪声)Rn1=-conj(sn1).*hn0+hn1.*conj(sn0);Rn2=hn2.*sn0+hn3.*sn1; Rn3

20、=-conj(sn1).*hn2+hn3.*conj(sn0);rn0= (Rn0+sgma.* (randn ( 1 , N2) +1 i*randn ( 1 , N2) ) ) ; %接收信号矢量(有噪声) rn1=(Rn1+sgma.*(randn(1,N2)+1i*randn(1,N2);rn2=(Rn2+sgma.*(randn(1,N2)+1i*randn(1,N2); rn3=(Rn3+sgma.*(randn(1,N2)+1i*randn(1,N2);sn_0=conj(hn0).*rn0+hn1.*conj(rn1)+conj(hn2).*rn2+hn3.*conj(rn3)

21、; sn_1=conj(hn1).*rn0-hn0.*conj(rn1)+conj(hn3).*rn2-hn2.*conj(rn3);for l2=1:N/4 %将接收到的信号合并Y2(2*l2-1)=sn_0(l2); Y2(2*l2)=sn_1(l2);endy2=DEQPSK(Y2); %QPSK解调y2=y2;n2,e2=bi terr(y2,S);% 计算误比特率 error2=error2+n2;%无信道编码(2发1收) hN0=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2); %瑞利信道函数 hN1=sqrt(0.5).*(randn(1,N2)+

22、1i.*randn(1,N2);RN0=hN0.*sn0+hN1.*sn1;%发射序列与接收序列的卷积RN1=-conj(sn1).*hN0+hN1.*conj(sn0);rN0=(RN0+sgma.*(randn(1,N2)+1i*randn(1,N2); %接收信号矢量(有噪声) rN1=(RN1+sgma.*(randn(1,N2)+1i*randn(1,N2);sN_0=conj(hN0).*rN0+hN1.*conj(rN1); sN_1=conj(hN1).*rN0-hN0.*conj(rN1);for l3=1:N/4%将接收到的信号合并Y3(2*l3-1)=sN_0(l3);

23、Y3(2*l3)=sN_1(l3);endy3=DEQPSK(Y3);%QPSK 解调y3=y3;n3,e3=bi terr(y3,S);% 计算误比特率 error3=error3+n3;endBER0(i)=error0/(LOOP*N);BER1(i)=error1/(LOOP*N);BER2(i)=error2/(LOOP*N);BER3(i)=error3/(LOOP*N);endsemilogy(SNR_dB,BER0,-b+);hold on;semilogy(SNR_dB,BER1,r);gridon;semilogy(SNR_dB,BER2,-g);gridon;semilo

24、gy(SNR_dB,BER3,-.k);gridon;legend(2Tx2R:卷积信道编码,2TxlR有信道编码,2Tx2R无信道编码,2TxlR无信道编 码);xlabel(SNR/dB);ylabel(BER);2 、viterbifunction decodingsequence=viterbi(m) trel = poly2trellis(9,557,663,711);tblen=7;%卷积码记忆深度decodingsequence=vitdec(m,trel,tblen,trunc,hard); decodingsequence=decodingsequence;3 、QPSKfu

25、nction modulatingsignal=QPSK(s)s=s;H=modem.pskmod(M,4,InputType, Bit); modulatingsignal=modulate(H,s); modulatingsignal=modulatingsignal;4 、DEQPSKfunction demodulatingsignal=DEQPSK(s)h= modem.pskdemod(M,4,OutputType, Bit);s=sdemodulatingsignal=demodulate(h,s);5 、convolutionalcode %卷积码编码器%msg信号 %cons

26、trainlength 约束长度(寄存器数目加一) %codegenerator冲击响应g (转换成八进制) function code=convolutionalcode(msg) constrainlength=9;trel = poly2trellis(constrainlength,557,663,711); code = convenc(msg,trel);五、仿真结果分析1、仿真图10I2Tx2R:卷积信道编码ZbdR有信道编玛2Fx2只无信道编码2T)dR无信道编阳DC山mIO71015SNR/dB202、结果分析 由误码率仿真图看出在四种信道编码的方法中：横坐标：信噪比 纵坐标

27、：误码率 信噪比在零到五之间时四种方式的误码率都差不多。 在信噪比大于五之后2x2卷积信道编码的误码率最低、仿真性能最好。 信噪比大于8左右时2x1无信道编码的误码率最咼、仿真性能最差。 2x1有信道编码和2x2无信道编码的仿真效果差不多。 有卷积码的误码率比无卷积码的低，因为卷积码有纠错功能。六、重点研究的问题七、心得与体会这次是我上大学以来觉得最为困难的一次小学期，这一次的任务综合了信道系统编码和 通信原理的知识与Matlab软件的应用。以前几乎没用过matlab，所以刚看到这个题目时我 简直觉得我从下手。Ppt看了好多遍也不知道到底做什么。通过翻阅通原课本、上网查阅资 料，了解ppt上给

28、出的每个函数的功能与用途，经过两个星期终于完成了本次任务。在这次 任务重，我碰到了许许多多的问题，其中最复杂的部分就是调制部分的代码与瑞利信道函数 部分的代码。借助网上找的资料和在同学的帮助下才完成了这次的任务。这次小学期让我获 益匪浅，首先我温习了通信原理的知识，而且学习了 matlab的一些的用法。最开始的迷茫 到最后的完成，感觉万事开头难，只要有足够的决心与毅力，没有什么是完成不了的。希望 下一次我可以做的更好。八、参考文献1 吴伟陵，牛凯，移动通信原理，电子工业出版社。2 张梧桐.MIMO通信技术导论J.科学周刊，2013(5).3 周炯槃、庞沁华、续大我、吴伟陵、杨鸿文通信原理，北京邮电大学出版社4 啜钢、王文博、常永宇、全庆一移动通信原理与系统，北京邮电大学出版社