通信原理课程设计

1、通信原理课程设计报告题目：基于MATLAB的M-QAM调制及相干解调的设计与仿真班级:姓名：杨仕浩(2014111347)解博文(2014111321)介子豪（201411132?指导老师:成绩：日期：2016年12月21日基于MATLAB勺M-QAh调制及相干解调的设计与仿真摘要：正交幅度调制技术（QAM）是一种功率和带宽相对高效的信道调制技术，因此在自适应信道调制技术中得到了较多应用。本次课程设计主要运用MATLAB软件对/17=16进制正交幅度调制系统进行了仿真，从理论上验证16进制正交幅度调制系统工作原理，为实际应用和科学合理地设计正交幅度调制系统，提供了便捷、高效、直观的重要方法。实

2、验及仿真的结果证明，多进制正交幅度调制解调易于实现，且性能良好，是未来通信技术的主要研究方向之一，并有广阔的应用前景。关键词：正交幅度调制系统；MATLAB;仿真1引言11.1 课程设计的目的11.2 课程设计的基本任务和要求11.3 仿真平台Matlab12QAM系统的介绍22.1 正交幅度调制技术22.2 QAM调制解调原理52.3 QAM的误码率性能73多进制正交幅度（M-QAM）调制及相干解调原理框图94基于MATLAB的多进制正交幅度（M-QAM调制及柜干解调设计与仿真104.1 系统设计104.2 随机信号的生成104.3 星座图映射114.4 波形成形（平方根升余弦滤波器）134

3、.5 调制144.6 加入高斯白噪声之后解调155仿真结果及分析206总Z吉与体会236.1 总结236.2 心得体会24【参考文献】25附录261引言本次课程设计主要运用MATLAB软件进行程序编写。实现模拟基带信号经QAM调制与相干解调的传输过程，通过分析比较调制解调输出波形以及功率谱特征，理解QAM调制解调原理。1.1 课程设计的目的通信原理课程设计的目的是使我们加深对所学的通信原理知识的理解，扎实掌握通信原理的基础知识和基本理论，增强分析问题和解决问题的能力，培养学生专业素质，提高其利用通信原理知识处理通信系统问题的能力，为今后专业课程的学习、毕业设计和工作打下良好的基础。1.2 课程

4、设计的基本任务和要求(1)设计出规定的数字通信系统的结构，包括信源，调制，发送滤波器模块，信道，接受滤波器模块以及信宿；(2)根据通信原理，设计出各个模块的参数(例如码速率，滤波器的截止频率等)；(3)熟悉MATLAB环境下的Simulink仿真平台,用Matlab/Simulink实现该数字通信系统；(4)观察仿真并进行波形分析(波形图和频谱图等)；1.3 仿真平台MatlabMATLAB目前已发展成为由MATLAB语言、MATLAB工作环境、MATLAB图形处理系统、MATLAB数学函数库和MATLAB应用程序接口五大部分组成的集数值计算、图形处理、程序开发为一体的功能强大的系统。本次课程

5、设计则在深入研究连续时间信号傅里叶级数分析理论知识的基础上，利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能,通过MATLAB编程进行图形功能仿真，从而实现连续时间周期信号频域分析的仿真波形，包括以下内容：用MATLAB实现周期信号的傅里叶级数分解与综合的波形；用MATLAB实现周期信号的单边频谱及双边频谱的波形与分析；用MATLAB实现典型周期信号的频谱的波形。102QAM系统的介绍2.1 正交幅度调制技术正交振幅调制(QuadratureAmplitudeModulation,QAM)是一种振幅和相位联合键控。虽然MPSK和MDPSK等相移键控的带宽和功率方面都具有优势,即

6、带宽占用小和比特噪声比要求低。图18PSK信号相位但是由图1可见，在MPSK体制中，随着M的增大，相邻相位的距离逐渐减小，使噪声容限随之减小，误码率难于保证。为了改善在M大时的噪声容限,发展出了QAM体制。在QAM体制中，信号的振幅和相位作为两个独立的参量同时受到调制。这种信号的一个码元可以表示为：SkAkcos(otk)kTt(k1)T(ii)式中：k二整数;Ak和k分别可以取多个离散值。式Q1)可以展开为：Sk(t)Akcoskcos4Aksinksin4(12)令Xk=AkCOSkYk=-Aksink则式(11)变为：Sk(t)XkcosotY;sint(13)和XkYk也是可以取多个离

7、散的变量。从式(1 3)看出，Sk(t)可以看作是两个正交的振幅键控信号之和在式（11）中，若k值仅可以取/4和-/4,Ak值仅可以取+A和A,则此QAM信号就成为QPSK信号，如图2所示：图24QAM信号矢量图所以，QPSK信号就是一种最简单的QAM信号。有代表性的QAM信号是16进制的，记为16QAM,它的矢量图示于下图中：图316QAM信号矢量图图中用黑点表示每个码元的位置，并且示出它是由两个正交矢量合成的。类似地，有64QAM和256QAM等QAM信号，如图4、图5所示。它们总称为MQAM调制。由于从其矢量图看像是星座，故又称星座调制。图464QAM信号矢量图图5256QAM信号矢量图

8、16QAM信号的产生方法主要有两种。第一种是正交调幅法，即用两路独第二种方法是复合相移法，它移用两路独立的QPSK信号叠加，形成16QAM信号，如图7所示。图中虚线大圆上的4个大黑点表示一个QPSK信号矢量的位置。在这4个位置上可以叠加上第二个QPSK矢量，后者的位置用虚线小圆上的4个小黑点表示。2.2 QAM调制解调原理2.2.1 QAM调制正交幅度调制QAM是数字通信中一种经常利用的数字调制技术，尤其是多进制QAM具有很高的频带利用率，在通信业务日益增多使得频带利用率成为主要矛盾的情况下，正交幅度调制方式是一种比较好的选择。正交幅度调制(QAM)信号采用了两个正交载波cos2ftsin2f

9、ct,每个载波都被一个独立的信息比特序列所调制。发送信号波形如图8所示Um(t)AmcgT(t)cos2fctAmsgT(t)sin2fct,m1,2，Ml(1-4)*图8M=16QAM信号星座图式中加和枷是电平集合，这些电平是通过将k比特序列映射为信号振幅而获得的。例如一个16位正交幅度调制信号的星座图如下图所示,该星座是通过用M=4PAM信号对每个正交载波进行振幅调制得到的。利用PAM分别调制两个正交载波可得到矩形信号星座。QAM可以看成是振幅调制和相位调制的结合。因此发送的QAM信号波形可表7F为Umn(t)AmgT(t)COS(2fn),m1,2,.J111,n1,2,.但！2,(i-

10、5)如果町23M2盯那么QAM方法就可以达到以符号速率2同时发送IO32M1个二进制数据。图9给出了QAM调制器的框图。二遗副数篡送滤祓%平衡调制黑Slli27Tftt衷送滤I器式t)图9QAM调制器框图222QAM的解调和判决假设在信号传输中存在载波相位偏移和加性高斯噪声。因此r可以表示为r(t) AmcgT(t)C0S(2 fc)Amsgi(t)sin(2 f)n(t)(1-6)其中是载波相位偏移，且n(t) n/t)cos2 fct n$2 fct将接收信号与下述两个相移函数进行相关运算i(t) gT(t)cos(2 fct )2(t) gi(t)sin(2 fct )d-9)如图10所

11、示，相关器的输出抽样后输入判决器。使用图10中所示的锁相环估算接收信号的载波相位偏移，相移1，和2m对该相位偏移进行补偿。图10 QAM信号的解调和判决输出判注L假设图中所示的时钟与接收信号同步，以使相关器的输出在适当的时刻及时被抽样。在这些条件下两个相关器的输出分别为iAnccosn5sin(1-1。)moA.2八ncsinn5cos(1-11)me其中1T2ncnc20(t)gi(t)dt(1-12)Unns20s(t)gT(t)dtQ-13)噪声分量是均值为0,方差为N。2的互不相关的高斯随机变量。最佳判决器计算距离量度2D(r,Sm)rSmm1,2,M(1-14)2.3QAM的误码率性

12、能矩形QAM言号星座最突出的优点就是容易产生PAM信号可直接加到两个正交载波相位上，此外它们还便于解调。K对于M2下的矩形信号星座图（K为偶数），QAM信号星座图与正交载波上的两个PAM信号是等价的,这两个信号中的每一个上都有.M2s个信号点。因为相位正交分量上的信号能被相干判决极好的分离，所以易于通过PAM的误码率确定QAM的误码率M进制QAM系统正确判决的概率是：2Pc1Pm(1T5)式中P晃.M进制PAM系统的误码率，该PAM系统具有等价QAM系统的每一个正交信号中的一半平均功率。通过适当调整M进制PAM系统的误码率，可得：P_21l/.MQ3/M1Eav/No(1-16)式中Eav/N

13、o是每个符号的平均信噪比。因此，因此M进制QAM的误码率为：Pm11Pm(1T7)可以注意到，当K为偶数时，这个结果对M情形时精确的，而当K为奇数时，就找不到等价的，M进制PAM系统。如果使用最佳距离量度进行判决的最佳判决器,可以求出任意K=1误码率的的严格上限。P.112Qr3Eav/M1NoX)1LPF抽样判决图12相干解调原理框图如图12所示的M-QAM目干解调原理框图首先对收到的已调信号进行正交相干解调，一路与cos3ct相乘，一路与sin3c相乘。然后经过低通滤波器,低通滤波器LPF滤除乘法器法产生的高频分量，获得有用信号，低通滤波器LPF输出经抽样判决可恢复出电平信号。4基于MAT

14、LA口的多进制正交幅度（M-QAM调制及相干解调设计与仿真4.1 系统设计首先进行系统的分析的设计，整个设计分为如下几个部分：随机序列的产生、序列的串/并和并/串转换、16QAM调制、星座图的绘制、16QAM解调、加入噪声、误码率的测量及仿真图形显示。4.2 随机信号的生成利用Matlab中的random_binary函数来产生0、1等概分布的随机信号。源代码如下所示：random_binary.m%产生二进制信源随机序列functioninfo=random_binary(N)ifnargin=0,%如果没有输入参数,则指定信息序列为10000个码元N=10000;end;fori=l:Nz

15、temp=rand;if(temp0.5),info(i)=0;%1/2的概率输出为0elseinfo(i)=l;%1/2的概率输出为1endend;4.3 星座图映射对产生的二进制随机序列进行串并转换，分离出I分量、Q分量，然后再分别进行电平映射。由于是调用matlab系统函数调制解调,在此将转换后边的序列进行四进制转换，方便后面的调制，再将转换好的序列通过调用qam()函数进行16qam调制，具体代码如下：代码如下：%串/并变换分离出I分量、Q分量，然后再分别进行电平映射I=x(l:2:nn-1);I/In=two2four(I,4*m);Q=x(2:2:nn);Q/Qn=two2four

16、(Q,4*m);ifKbase=2;%基带成形滤波I=bshape(I,fs,fb/4);Q=bshape(Q/fs/fb/4);end;y=L*cos(2*pi*fc*t)-Q.*sin(2*pi*fc*t);%调制二进制转换成进制代码：two2four.m%二进制转换成四进制function y,yn =two2four(xzm);T=0 1;3 2;n = length(x);ii=l;for i = l:2:n-l; xi=x(i:i + l) + l; yn(ii)=T(xi(l),xi(2);11 = 11 + 1;end; yn=yn-1.5; y=yn;for i = l:m-

17、l; y=y;yn; end;y=y(：),； %映射电平分别为-1.5 ； 0.5 ； 0.5 ； 1.5星座图代码如下：constel.m%画出星座图function c=constel(x,fs,fb/fc);N=length(x);m=2*fs/fb;n=fs/fc;il=m-n;i=l;ph0=(il-l)*2*pi/n;whilei=N/m;xi=x(il:il+n-l);y=2*fft(xi)/n;c(i)=y；i=i+l;il=il+m;end;%如果无输出，则作图ifnargoutl;cmax=max(abs(c);ph=(0:5:360)*pi/180;plot(1.414

18、*cos(ph)/1.414*sin(ph)/c1)；holdon;fori=l:length(c);ph=phO-angle(c(i);a=abs(c(i)/cmax*1.414;plot(a*cos(ph)za*sin(ph)/r*);end;plot(-1.51.5,00J-1.51.5；k:);holdoff;axisequal;axis(-1.51.5-1.51.5);end;4.4 波形成形（平方根升余弦滤波器）为了避免相邻传输信号之间的串扰，多元符号需要有合适的信号波形。方波是在本地数字信号处理时常见的波形，但在实际传输时这种方波并不合适。根据奈奎斯特第一准则，在实际通信系统中一

19、般均使接收波形为升余弦滚降信号。这一过程由发送端的基带成形滤波器和接收端的匹配滤波器两个环节共同实现，因此每个环节均为平方根升余弦滚降滤波，两个环节合成就实现了一个升余弦滚降滤波。实现平方根升余弦滚降信号的过程称为波形成形，通过采用合适的滤波器对多元码流进行滤波实现，由于生成的是基带信号，因此这一过程又称基带成形滤波。代码如下：bshape.m%基带升余弦成形滤波器functiony=bshape(x,fs/fb/N,alfa/delay);%设置默认参数ifnargin6;delay=8;end;ifnargin5;alfa=0.5;end;ifnargin=0.5&yn=1.5&yn=2.

20、5);yn(I3)=ones(size(I3)*3;%一位四进制码元转换为两位二进制码元T=00;01;11;10;n=length(yn);fori=l:n;xn(i/:)=T(yn(i)+l/:);end;xn=xn;xn二xn(:);xn=xn;误码率代码如下：numoferr=0;fori=l:Nif(y_output(i)=info(i),numoferr=numoferr+l;end;end;Pe(j)=numoferr/N;%统计误码率end;figure;semilogy(SNR_in_dB,Peered*-);gridon;xlabel(SNRindB);ylabel(Pe)

21、;title(16QAM调制在信道噪声强度下的误码率,)；16-QAM载波调制信号在AWGN信道下的性能代码如下：clearallnsymbol=100000;%每种信噪比下的发送符号数M=16;%16-QAMgraycode=0132457612131514891110;%Gray编码规则EsN0=5:20;%信噪比，Es/NOsnrl=10.A(EsN0/10);%信噪比转换为线性值msg=randint(l/nsymbol/M);%消息数据msgl=graycode(msg+l);%Gray映身寸msgmod=qammod(msgl/M);%基带16-QAM调制spow=norm(msg

22、mod)42/nsymbol;%求每个符号的平均功率forindx=l:length(EsNO)sigma=sqrt(spow/(2*snrl(indx);%根据符号功率求噪声功率rx=msgmod4-sigma*(randn(lzlength(msgmod)4-j*randn(ljength(msgmod)y=qamdemod(rx/M);decmsg=graycode(y+1);err,ber(indx)=biterr(msgzdecmsg/log2(M);%误比特率errlser(indx)=symerr(msg/decmsg);%误符号率endP4=2*(l-l/sqrt(M)*qfu

23、nc(sqrt(3*snrl/(M-l);Aserl=l-(l-P4).2;%理论误符号率berl=l/log2(M)*serl;%理论误比特率semilogy(EsNO,ber,-ko,EsNO,ser-k*,EsNO,serLEsNO,berl,-k.);title(16-QAM载波调制信号在AWGN信道下的性能)xlabel(Es/NO)；ylabelC误比特率和误符号率)legend(误比特率，误符号率，理论误符号率，理论误比特率)5仿真结果及分析0.6-0.8已调信号(ln:red,Qn:green)0.40.20-0.2-0.6515102025图13已调信号波形如图13所示的已调

24、信号，程序码元数量设定为10000个,此处显示1/10,即1000个,基带信号频率1HZ,抽样频率32H乙载波频率4HZ。-30-800已调信号频谱-50-55-6526101248 f/fb1416图14已调信号频谱20星座图1匚图1516QAM星座图通过观察图15的16QAM星座图可以看出16QAM的频带利用率比较高,带宽占用比较小。16QAM调制在不同信道噪声强度下的误码率图1616QAM调制在不同信道噪声强度下的误码率3116-QAM载波调制信号在AWGN信道下的性能-610爸S.D夺n天口爸寺匕n天10Es/N01520图1716-QAM载波调制信号在AWGN信道下的性能由图16图1

25、7可看到当信噪比小的情况下，仿真曲线和理论曲线差距略大,而随着信噪比的增大，仿真曲线越来越逼进理论曲线。在同样信噪比时，误符号率比误比特率（误码率）要大。简单分析不难看出，由于理论误码率曲线是建立在误符号率除以4的基础上的，而这一条件的前提是出现误符号的时候，一个符号中只有一个bit位发生了错误，这表明误码率比较低，也就是说明信噪比比较大。所以，当信噪比比较小的时候，理论计算的误码率的值要小于仿真得到的值。6总结与体会6.1 总结在现代通信中，提高频谱利用率一直是人们关注的焦点之一。近年来，随着通信业务需求的增长，寻找频谱利用率高的数字调制方式已成为数字通信系统设计，研究的主要目标之一o正交振

26、幅调制QAM(QuadratureAmplitudeModulation)就是一种频谱利用率很高的调制方式，其在中、大容量数字微波通信系统、有线电视网络高速数据传输、卫星通信系统等领域得到了广泛应用。本设计是基于MATLAB/Simulink的16QAM调制与解调系统进行设计与仿真，得到的结论是：(1)对16QAM调制解调系统基本原理进行了较为深入的了解与分析，并且根据原理图构建MATLAB/Simulink的仿真模型。(2)较为熟悉的掌握MATLAB/Simulink软件在通信系统中的设计与仿真的步骤与方法。(3)通过观察星座图可以看出16QAM的频带利用率比较高，带宽占用比较小。作为一种现

27、代调制技术，QAM相比较传统调制技术，有着很多优于传统调制技术的特性，使得它在中、大容量数字微波通信系统、有线电视网络高速数据传输、卫星通信系统和蜂窝系统等系统中得到大规模应用，大大提高系统通信质量和通信效率。MATLAB在通信仿真中有着重要的应用,MATLAB/Simulink是通信系统计算机仿真的强大工具，本文提供了一个实际仿真的例子。实际的信道是很复杂的，在实际的应用中应根据不同的要求选用不同的调制方式。本文利用MATLAB/Simulink对M=16进制正交幅度调制系统进行了仿真，从理论上验证16进制正交幅度调制系统工作原理，为实际应用和科学合理地设计正交幅度调制系统，提供了便捷、高效

28、、直观的重要方法。6.2 心得体会这次课程设计通过运用MATLAB/Simulink软件进行设计仿真,掌握了运用技巧；理解了M-QAM调制的原理以及其优越性。课程设计是培养学生综合运用所学知识，发现、提出、分析和解决问题，锻炼实践能力的重要环节，是对学生实际工作能力的具体训练和考察过程。通过课程设计系统的了解了一些理论知识，能把理论知识与具体实践相结合，培养了我们动手、动脑能力。课程设计主要是对理论知识的延伸，它的目的主要是要在设计中发现问题，并且自己要能找到解决问题的方案，形成一种独立的意识，还能从设计中检验平时所学的理论知识到底有多少，巩固已经学会的内容，不断学习新知识，把这门课学的扎实。

29、参考文献】邵玉斌.Matlab/Simulink通信系统建模与仿真实例分析.北京:清华大学出版社，20082张化光，刘鑫蕊，孙秋野.MATLAB/SIMUUNK实用教程.北京:人民邮电出版社，20093樊昌信，曹丽娜.通信原理.北京:国防工业出版社,20084刘卫国.MATLAB程序设计教程.北京:中国水利水电出版社，2005附录main_plot.mclear;clc;echooff;closeall;N=10000;%设定码元数量fb=l;%基带信号频率fs=32;%抽样频率%载波频率，为便于观察已调信号，我们把载波频率设fc=4;的较低Kbase=2;%Kbase=l,不经基带成形滤波，

30、直接调制；%Kbase=2基带经成形滤波器滤波后，再进行调制info=random_binary(N);%产生二进制信号序列y/I,Q=qam(info/Kbase/fs/fb/fc);%对基带信号进行16QAM调制yl=y；y2=y；%备份信号，供后续仿真用T=length(info)/fb;m=fs/fb;nn=length(info);dt=l/fs;t=0:dt:T-dt;n=length(y);y=fft(y)/n;y=abs(y(l:fix(n/2)*2;q=find(yle-04);y(q)=le-04;y=20叫oglO(y);fl=m/n;f=O:fl:(length(y)-

31、l)*fl;%subplot(212);plot(f,y,b);gridon;title(,已调信号频谱,)；xlabel(f/fb)；%画出16QAM调制方式对应的星座图%constel(ylrfszfbzfc);title(星座图)；SNRJn_dB=8:2:24;%AWGN信道信噪比forj=l:length(SNR_in_dB)y_add_noise=awgn(y2,SNR_in_dB(j);%力口入不同强度的高斯白噪声y_output=qamdet(y_add_noise/fszfb/fc);%对已调信号进行解调numoferr=0;fori=l:Nif(y_output(i)=in

32、fo(i)znumoferr=numoferr+l;end;end;Pe(j)=numoferr/N;%统计误码率end;figure;semilogy(SNR_in_dB/Pe/,blue*-);gridon;xlabel(SNRindB);ylabelCPe);title(116QAM调制误码率1);bshape.m%基带升余弦成形滤波器functiony=bshape(x,fs/fb/N/alfa/delay);%设置默认参数ifnargin6;delay=8;end;ifnargin5;alfa=0.5;end;ifnargin4;N=16;end;b=firrcos(N,fb/2*a

33、lfa*fb/fs);y=filter(b/l/x);four2two.mfunctionxn=four2two(yn);y=yn;ymin=min(y);ymax=max(y);ymax=max(ymaxabs(ymin);ymin=-abs(ymax);yn=(y-ymin)*3/(ymax-ymin);%设置门限电平，判决I0=find(yn=0.5&yn=1.5&yn=2.5);yn(I3)=ones(size(I3)*3;%一位四进制码元转换为两位二进制码元T=00;01;11;10;n=length(yn);fori=l:n;xn(i/:)=T(yn(i)+l/:);end;xn=

34、xn;xn=xn(:);xn=xn;two2four.m%二进制转换成四进制functiony/yn=two2four(x,m);T=01;32;n=length(x);ii=l;fori=l:2:n-l;xi=x(i:i+l)+l;yn(ii)=T(xi(l)rxi(2);.11=11+1;end;yn=yn-1.5;y=yn;y=y;yn；end;y=y(：),；%映射电平分别为-1.5；0.5；0.5；1.5random_binary.mfunctioninfo=random_binary(N)if nargin = 0,%如果没有输入参数，则指定信息序列为10000个码元N=10000;end;fori=l:N,temp=rand;% 1/2的概率输出为0% 1/2的概率输出为1if(temp0.5)zinfo(i)=0;elseinfo(i)=l;endend;qamde