16QAM调制与解调的MATLAB实现及调制性能分析.docx

《16QAM调制与解调的MATLAB实现及调制性能分析.docx》由会员分享，可在线阅读，更多相关《16QAM调制与解调的MATLAB实现及调制性能分析.docx（12页珍藏版）》请在第一文库网上搜索。

1、通信原理课程设计报告书16QAM调制与解调课题名称的MAT1AB实现及调制性能分析姓名学号学院通信与电子工程学院专业通信工程指导教师一、设计任务及要求：设计任务：利用MAT1AB设计一个16QAM调制与解调系统，并对其进行性能分析。要求：1 .设封一一个16QAM调制与解调系统。2 .设计程序时必须使得程序尽可能的简单。3 .利用MAT1AB进行程序编写并对系统进行仿真分析。指导教师签名：20年月日二、指导教师评语:指导教师签名：2010年月日三、成绩验收盖章2010年月日16QAM调制与解调的MAT1AB实现及调制性能分析1设计目的(1)掌握16QAM调制与解调的原理。(2)掌握星座图的原理

2、并能熟悉星座图的应用。(3)熟悉并掌握MAT1AB的使用方法。(4)通过对16QAM调制性能的分析了解16QAM调制相对于其它调制方式的优缺点。2设计原理正交振幅调制(QUadratUreAmPIitUdeModUIation,QAM)是一种振幅和相位联合键控。虽然MPSK和MDPSK等相移键控的带宽和功率方面都具有优势，即带宽占用小和比特噪声比要求低。但是由图1可见，在MPSK体制中，随着M的增大，相邻相位的距离逐渐减小，使噪声容限随之减小，误码率难于保证。为了改善在M大时的噪声容限，发展出了QAM体制。在QAM体制中，信号的振幅和相位作为两个独立的参量同时受到调制。这种信号的一个码元可以表

3、示为sk(Z)=Akcos(0tk)kTt(k+Y)T(21)式中：k=整数；4和4分别可以取多个离散值。式(2-1)可以展开为(22)sk(Z)=Akcoskcost-AksinksinOyXk=AkCOSek，Yk=-AkSinok则式(2-1)变为(23)sk(t)=Xkcosj+YksinCOJXK和匕也是可以取多个离散的变量。从式(23)看出，与Q)可以看作是两个正交的振幅键控信号之和。在式(21)中，若纵值仅可以取4和-4,A4值仅可以取+A和则此QAM信号就成为QPSK信号，如图2所示：所以，QPSK信号就是一种最简单的QAM信号。有代表性的QAM信号是16进制的，记为16QAM

4、,它的矢量图示于下图中：图316QAM信号矢量图图中用黑点表示每个码元的位置，并且示出它是由两个正交矢量合成的。类似地,有64QAM和256QAM等QAM信号，如图4、图5所示。它们总称为MQAM图464QAM信号矢量图图5256QAM信号矢量图调制。由于从其矢量图看像是星座，故又称星座调制。16QAM信号的产生方法主要有两种。第一种是正交调幅法，即用两路独立的正交4ASK信号叠加，形成16QAM信号，如图6所示。第二种方法是复合相移法，它用两路独立的QPSK信号叠加，形成16QAM信号，如图7所示。图中图7复合相移法虚线大圆上的4个大黑点表示一个QPSK信号矢量的位置。在这4个位置上可以叠加

5、上第二个QPSK矢量，后者的位置用虚线小圆上的4个小黑点表示。3设计过程3.1设计思路由设计原理中可知MQAM调制又称为星座调制，故我们在设计16QAM调制系统时就可以星座图来进行编程。下面我们就借用如图8所示的星座图设计一个16QAM调制系统。1011IOObIO1Ot1QOO(X)O100。(W1o(Mmw00110010-01010111图816QAM星座在图中共有16个点，每个点用4个比特表示，代表调制以后的一个矢量位置(这个点拥有唯一的振幅与相位)。因此我们可以把横轴看作是实轴，纵轴看作虚轴。由于每个点与跟它相邻的四个点是等距，且设为2,则每个点都可用一个虚数进行表示。例如点0000

6、可用T-j表示，这个虚数的模就是相当于16QAM信号的振幅，相角就相当于16QAM信号的相位。所以16QAM的调制过程就可以用如下语句进行描述：ifA(1,bx)=0000B(k)=-1-1i;e1seifA(1,bx)=0001B(k)=-3-1i;e1seifA(1,bx)=0010B(k)=-1-3i;e1seifA(1,b:c)=O011B(k)=-3-3i;e1seifA(1,bx)=0100B(k)=1-1i;e1seifA(1,b:c)=O1O1B(k)=1-3i;e1seifA(1,bx)=0110B(k)=3-1i;e1seifA(1,b:c)=O111B(k)=3-3i;e

7、1seifA(1,b：c)=1000B(k)=-1+1i;e1seifA(1,bx)=1001B(k)=-1+3i;e1seifA(1,b:c)=1O10B(k)=-3+1i;e1seifA(1,b:c)=1O11B(k)=-3+3i;e1seifA(1,bx)=11O0B(k)=1+1i;e1seifA(1,bx)=1101B(k)=3+1i;e1seifA(1,b:c)=1110B(k)=1+3i;e1seifA(1,bx)=1111B(k)=3+3i;end当调制以后的信号经过信道加噪以后，我们必须对其进行解调。由于加噪了的缘故，调制以后的信号不再是原来的信号，而应该有不同。因此在解调时

8、不能简单的将上述过程逆转，即不能由-1-j就判断为0000。而应该对虚数的实部和虚部设定一个范围后再进行判断。这个范围边界的选取原理我们可以借用量化的概念,取相邻虚数的实部的平均数和虚部的平均数。以下为16QAM解调过程的程序语句：if(rea1(D(n)-2)&(imag(D(n)-2)C(1,dx)=0011;e1seif(rea1(D(n)-2)&(imag(D(n)0)C(1,d:e)=OOO1;e1seif(rea1(D(n)-2)&(imag(D(n)2)C(1,d:e)=1010;e1seif(rea1(D(n)=2)C(1,dx)=1011;e1seif(rea1(D(n)0)

9、&(imag(D(n)-2)C(1,d:e)=OO10;e1seif(rea1(D(n)0)&(imag(D(n)0)C(1,d:e)=OOOO;e1seif(rea1(D(n)0)&(imag(D(n)2)C(1,d:e)=1OOO;e1seif(rea1(D(n)=2)C(1,d:e)=1OO1;e1seif(rea1(D(n)2)&(imag(D(n)-2)C(1,d:e)=O1O1;e1seif(rea1(D(n)2)&(imag(D(n)0)C(1,d:e)=O1OO;e1seif(rea1(D(n)2)&(imag(D(n)2)C(1,d:e)=11OO;e1seif(rea1(D

10、(n)=2)C(1,d:e)=1110;e1seif(rea1(D(n)=2)&(imag(D(n)=2)&(imag(D(n)=2)&(imag(D(n)=2)&(imag(D(n)=2)C(1,d:e)=1111;end3.2设计总程序c1eara11;c1osea11;N=40000;K=4*N;%信息长度1=7*N;W=7*N4;E=randsrc(1,K,0,1);%信源B=ZerOS(I,W);%16QAM调制后的信号C=zeros(1,1);%16QAM解调后的信号num=zeros(20,1);%误比特数ber=zeros(20,1);%误比特率forSNR=I:1:20%汉明

11、编码A=encode(E,7,4,hammingbinary,);%16QAM调制过程fork=1:Wb=4*k-3;c=4*k;ifA(1,bx)=0000B(k)=-1-1i;e1seifA(1,bx)=OOOIB(k)=-3-1i;e1seifA(1,b:c)=OO10B(k)=-1-3i;e1seifA(1,bx)=0O11B(k)=-3-3i;e1seifA(1,bx)=01O0B(k)=1-1i;e1seifA(1,bx)=01O1B(k)=1-3i;e1seifA(1,bx)=0110B(k)=3-1i;e1seifA(1,b:c)=O111B(k)=3-3i;e1seifA(1

12、tbx)=1000B(k)=-1+1i;e1seifA(1,bx)=1001B(k)=-1+3i;e1seifA(1,bx)=1O10B(k)=-3+1i;e1seifA(1,b:c)=(1O11B(k)=-3+3i;e1seifA(1,bx)=11O0B(k)=1+1i;e1seifA(1,bx)=11O1B(k)=3+1i;e1seifA(1,bx)=1110B(k)=1+3i;e1seifA(1,b:c)=1111B(k)=3+3i;endend%信道加噪D=awgn(B,SNR);%16QAM解调过程forn=1:Wd=4*n-3;e=4*n;if(rea1(D(n)-2)&(imag

13、(D(n)-2)C(1,dx)=0011;e1seif(rea1(D(n)-2)&(imag(D(n)0)C(1,dx)=0001;e1seif(rea1(D(n)-2)&(imag(D(n)2)C(1,d:e)=1O1O;e1seif(rea1(D(n)=2)C(1,d:e)=1O11;e1seif(rea1(D(n)0)&(imag(D(n)-2)C(1,d:e)=OO1O;e1seif(rea1(D(n)0)&(imag(D(n)0)C(1,dx)=0000;e1seif(rea1(D(n)0)&(imag(D(n)2)C(1,d:e)=1OOO;e1seif(rea1(D(n)=2)C(1,d:e)=1001;e1seif(rea1(D(n)2)&(imag(D(n)-2)C(1,dx)=01O1;e1seif(rea1(D(n)2)&(imag(D(n)0)C(1,d:e)=O1OO;e1seif(rea1(D(n)2)&(imag(D(n)2)C(1,d:e)=11OO;e1seif(rea1(D(