基于MATLAB的通信系统的设计与仿真

1、基于MATLAB的通信系统的设计 与仿真通信系统是用以完成信息传输过程的技术系统的总称，作用是将信息 从信源发送到一个或多个目的地。调制与解调在信息的传输过程中占据着重要 的地位，是不可或缺的，因此研究系统的调制和解调过程就极为重要。MATLAB 是集数值计算、图形绘制、图像处理及系统仿真等强大功能于一体的科学 计算语言，它强大的矩阵运算和图形可视化的功能以及丰富的工具箱，为 通信系统的调制和解调过程的分析提供了极大的方便。本论文首先介绍了通信系统的概念，进而引出调制和解调，然后介绍 了我们常用的几种调制和解调的方法。由于MATLAB具有的强大功能所以详 细介绍了 MATLAB通信系统工具箱，

2、并给出了基于MATLAB的通信系统的调制 与解调的实现，运用MATLAB仿真软件进行仿真。【关键词】：通信系统调制与解调 MATLAB1ABSTRACTCoinmunication system is used to complete the information transmission process of the floorboard of the technology system, is used to send information from the source to one or more destinations. Modulation and demodulation

3、in the process of information transmission occupies the important position, is indispensable, so the system of modulation and demodulation process is extremely important. MATLAB numerical computation9 graphics rendering, Image processing and system simulation, and powerful functions in the integrati

4、on of scientific computing language9 its powerful function of matrix and graph visualization and rich toolkit, for the analysis of the communication system of modulation and demodulation process provides a great convenience.This paper firstly introduces the concept of communication system, which lea

5、d to modulation and demodulation, and then introduces our commonly used several kinds of modulation and demodulation method.【Key words ： Communication system Modulation and demodulation MATLAB1目录前言4第一章通信系统的调制与解调5第一节通信系统的概念5一、模拟通信系统.5二、数字通信系统6三、QAM调制解调原理.14第二节多进制正交幅度调制及相干解调原理框图17一、正交调制原理框图.17二、相干解调原理

6、框图.18第二章MATLAB通信系统工具箱19第三章MATLAB的模拟调制和解调实例21第一节用MATLAB分析双边带幅度调制（DSMAM） 21第二节用MATLAB分析相干解调过程24第三节基于MATLAB的多进制正交幅度调制及相干解调设计与仿真.26一、对系统进行分析与设计26二、随机信号的生成26三星座图映射27四、波形成形（平方根升余弦滤波器）29五、调制.30六、加入高斯白噪声之后解调.31七、误码率曲线33A. 160AM载波调制信号在4GN信道下的性能.33第四节仿真结果及分析35结论39致谢40参考文献通信按照传统的理解就是信息的传输，在当今高度信息化得社会，信息和 通信已经成

7、为现代社会的“命脉”。信息作为一种资源，只有通过广泛的传播与 交流，才能产生利用价值，促进社会成员之间的合作，推动社会生产力的发展, 创造出巨大的经济效益。而通信作为传输信息的手段或方式，与传感技术、计 算机技术相互融合，已经成为21世纪国际社会和世界经济发展的强大推动力， 所以未来的通信对人们的生活方式和社会的发展将会产生更加重大和意义深远 的影响。在信息传输过程中，要求天线的尺寸要和信号的波长相比拟，信号才能有 效的被辐射。对于语音信号来说，相应的天线尺寸要在几十公里以上，实际上 不可能实现，所以需要经过调制将信号频谱搬移到较高的频率范围，如果不进 行调制就把信号直接辐射出去，那么各电台所

8、发出信号的频率就会相同。调制 作用的实质就是使相同频率范围的信号分别依托于不同频率的载波上，接收机 就可以分离出所需的频率信号，不致互相干扰。有时信号过于复杂，人工计算其调制和解调过程较难实现，对其结果的分 析又缺乏可视化的直观表现，影响了所得结果在实际生活中的应用，美国 MathWorks公司开发的MATLAB解决了这一问题。它应用于自动控制、数学计算、 信号分析、信号处理等诸多领域,也是国内高校和研究部门进行许多科学研究的 重要工具。MATLAB的出现给通信系统的分析提供了极大的方便。第一章通信系统的调制与解调第一节通信系统的概念通信是为了传输信息，通信系统就是将信息从信源发送到一个或多个

9、目的 地，对于电通信来说，首先要把消息转变成电信号，然后经过发送设备，将信 号送入信道，在接收端利用接受设备对接收信号作相应的处理后，送给信宿再 转换为原来的消息，这一过程可利用图1.1所示的通信系统一般模型来概括。图1.1通信系统的一般模型通信(Communication)传输的消息是多种多样的，可以分成两大类：一类称 为连续消息；另一类称为离散消息。消息的传递是通过它的物理载体电信号来 实现的，按信号参量的取值不同，可以把信号分为两类：模拟信号和数字信号。 通常按照信道中传输的是模拟信号还是数字信号，相应的可以把通信系统分为 模拟通信系统和数字通信系统。一、模拟通信系统模拟通信系统是利用模

10、拟信号来传递信息的通信系统，其模型如图1.2所 示，其中包含两种重要变换。第一种变换是，在发送端把连续消息变换成原始 电信号，在接收端进行相反的变换，这种变换由信源和信宿来完成，通常称为 原始电信号为基带信号，基带的意思是指信号的频谱从零频附近开始。有些信 道可以直接传输基带信号，而以自由空间作为信道的无线电传输却无法直接传 输这些信号。因此，模拟通信系统中常常需要进行第二种变换：把基带信号变换成适合在信道中传输的信号，并在接收端进行反变换。完成这种变换和反变 换的通常是调制器和解调器。图1.2模拟通信系统模型二、数字通信系统数字通信系统(Digital Communication Syste

11、m,DCS)是利用数字信号来传 递信息的通信系统，如图1.3所示。数字通信所涉及的技术问题很多，其中主要 有信源编码与译码、信道编码与译码、数字调制与解调、同步以及加密与解密 等。图1.3数字通信系统模型与模拟通信相比，数字通信具有许多优良的特性:(1)、抗干扰能力强，且噪声不积累。数字通信系统中传输的是离散取值 的数字波形，接收端的目标不是精确的还原被传输的波形，而是从收到噪声干 扰的信号中判决出发送端所发送的是哪一个波形。(2)、传输差错可控。在数字同喜系统中，可以通过信道编码技术进行检 错与纠错，降低误码率，提高传输质量。(3)、便于用现代数字信号处理技术对数字信息进行处理、变换、存储。

12、 这种数字处理的灵活性表现为可以将来自不同信源的信号综合到一起传输。(4)、易于集成，使通信设备微型化，重量轻。(5)、易于加密处理，且保密性好。数字通信的缺点是，一般需要较大的传输带宽，设备复杂。近年来，随着 大规模集成电路的出现，数字系统的设备复杂程度和技术难度大大降低了，数 字传输方式日益受到欢迎。数字处理的灵活性使得数字传输系统中传输的数字 信息既可以来自计算机、电传机等数据终端的各种数字代码，也可以来自模拟 信号经过数字化处理后的脉冲编码(PCM)信号等。原理上，数字信息可以宜 接用数字代码序列表示和传输，但在世纪传输中，视系统的要求和信道情况， 一般需要进行不同形式的编码，并且选用

13、一组取值有限的离散波形来表示，这 些取值离散的波形可以是未经调制的电信号，也可以是调制后的信号。未经调 制的数字信号所占据的频谱是从零频或很低频率开始，称为数字基带信号。1、调制和解调的概念调制：把信号转换成适合在信道中传输的形式的一种过程，广义的调制分 为基带调制和带通调制(也称为载波调制)。在无线通信中和其他大多数场合， 调制均指载波调制。载波调制，就是用调制信号去控制载波的参数的过程，使载波的某一个或 某几个参数按照调制信号的规律而变化。调制信号是指来自信源的消息信号(基 带信号)，这些信号可以是模拟的，也可以是数字的。未受调制的周期性振荡信 号称为载波，它可以是正弦波，也可以是非正弦波

14、(如周期性脉冲序列)。载波 调制后称为己调信号，它含有调制信号的全部特征。解调：将己调信号中的调制信号恢复出来，是调制的逆过程。调制方式不 同，解调方法也不一样。解调可以分为正弦波解调(有时也称为连续波解调)。 正弦波解调还可再分为幅度解调、频率解调和相位解调，此外还有一些变种如 单边带信号解调、残留边带信号解调等。同样，脉冲波解调也可分为脉冲幅度 解调、脉冲相位解调、脉冲宽度解调和脉冲编码解调等。调制方式有很多，根据调制信号是模拟信号还是数字信号，载波是连续波 （通常是正弦波）还是脉冲序列，相应的调制方式有模拟连续波调制（简称模 拟调制）、数字连续波调制（简称数字调制）、模拟脉冲调制和数字脉

15、冲调制等。 最重要和最常用的模拟调制方式是用正弦波作为载波的幅度调制和角度调制。 常见的调幅、双边带、单边带和残留边带等调制就是幅度调制的几个典型实例。解调的方法可以分为两类：相干解调和非相干解调（包络检波）。解调过程 与采用何种解调方式有关，对于常规幅度调制，一般用包络检波进行解调，由 于在这种解调方式中，接收机对载波频率和相应精度的了解是无关紧要的，所 以解调过程相对简单。对于DSB-AM调制和SSB-AM调制，用相干解调的方法, 它要求在接收机中有一个与载波同频同相的信号，接收机中产生所需要的正弦 波振荡器，为本地振荡器。数字通信系统中，发送信号某一参数的离散取值（如正弦波的振幅、频率、

16、 相位等）于所承载消息的数字信号（即离散符号）之间通常是一种一一对应关 系。如同模拟信号的频带传输时一样，传输数字信号时也有三种基本的调制方 式：振幅键控（ASK）、频移键控（FSK）和相移键控（PSK）。它们分别对应于 利用载波（正弦波）的振幅、频率和相位来承载数字基带信号，可以是模拟线 性调制和角度调制的特殊情况。理论上数字调制与模拟调制在本质上没有什么不同，它们都属于正弦波调 制。但是，数字调制是源信号为离散型的正弦波调制，而模拟调制则是源信号 为连续型的正弦波调制，因而，数字调制具有由数字信号带来的一些特点。这 些特点主要包括两个方面：第一，数字调制信号的产生，除把数字的调制信号 当做

17、模拟信号的特例而直接采用模拟调制方式产生数字调制信号外，还可以采 用键控载波的方法。第二，对于数字调制信号的解调，为提高系统的抗噪声性 能，通常采用与模拟调制系统中不同的解调方式。数字信息由二进制和多进制之分，因此，数字调制可分为二进制调制和多 进制调制。在二进制调制中，信号参量只有两种可能取值；而在多进制调制中, 信号参量可能有M （M2）中取值。2、幅度调制和相干解调的原理1)幅度调制(线性调制)的原理幅度调制是由调制信号去控制高频载波的幅度，使之随调制信号作线性变 化的过程，设正弦型载波为c(t)=Acos( coct(p)(2-1)式中：A为载波幅度；以为载波角频率；8)为载波初始相位

18、(以后可以假 定为夕。为0,而不失讨论的一般性)。根据调制定义，幅度调制信号(已调信号)一般可表示成Sm(t)=Am(t)cos (coct)(2-2)式中：m(t)为基带调制信号。设调制信号m的频谱为则由上式(2.3-2)不难得到已调信号Sm(t) 的频谱S3)S(g) = M(co + 公)+ Af 3 - coc)(2-3)由以上表示式可见，在波形上，幅度已调信号的幅度随基带信号的规律而 呈正比地变化；在频谱结构上，它的频谱完全是基信号频谱在频域内的简单搬 移(精确到常数因子)。由于这种搬移是线性的，因此，幅度调制通常又称为线 性调制，这里的“线性”并不意味着已调信号与调制信号之间符合线

19、性变换关系, 实际上，任何调制过程都是一种非线性的变换过程。标准调幅就是常规双边带调制，简称调幅(AM)。假设调制信号m的平 均值为0,将其叠加一个直流偏量A0后与载波相乘(如图1.4所示)，即可形成 调幅信号。Sm(t)图L4 AM调制模型其时域表示式为(2-4)SAM(t)=A0+m(t) COS6Z7 =A0 cosxt7+m(t) coscoct式中：AO为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。若m(t)是确知信号，则AM信号的频谱为SAM(co )=fAO 5(。+ tzr) + dco- coc) +l/2 M(co+ 公)+ M(0- coc) (2-5)2)

20、相干解调相干解调也叫同步检波，解调与调制的实质一样，均是频谱搬移。调制是 把基带信号的频谱搬到了载波位置，这一过程可以通过一个相乘器与载波相乘 来实现。解调则是调制的反过程，即把在载频位置的已调信号的频谱搬回到原 始基带位置，因此同样可以用相乘器与载波相乘来实现，相干解调的一般模型 如图1.5所示。Sm(t)Sp(t) Sd(t)vCHt lp ccosart图L5相干解调器的一般模型相干解调时，为了无失真的恢复原基带信号，接收端必须提供一个与接收 的已调载波严格同步(同频同相)的本地载波(称为相干载波)，它与接收的已 调信号相乘后，经低通滤波器取出低频分量，即可得到原始的基带调制信号。相干解

21、调器适用于所有的线性调制信号的解调。送入解调器的己调信号的 一般表达式为 Sm(t)=SI(t) cos以7 +SQ(t) sin coct(2-6)与同频同相的相干载波c(t)相乘后，得Sp(t)= Sm(t) coscoct =1/2 SI(t)+l/2 St) cos2yf +1/2SQ(t) sin 2公7(2-7)经低通滤波器(LPF)后，得Sd(t)=1/2 SI(t)(2-8)其中，SI是m通过一全通滤波器Hl (g)后的结果。因此，Sd(t)就是解调输出，即Sd(t)=l/2 SI(t) m(t)由此可见，相干解调器适用于所有线性调制信号的解调，即对于AM、DSB、 SSB和V

22、SB都是适用的。只是AM信号的解调结果中含有直流成分A0,这时 在解调后加上一个简单的隔直流电容即可。3、正交幅度调制解调原理1)正交幅度调制技术正交振幅调制(Quadrature Amplitude Modulation,QAM)是一种振幅和相 位联合键控。虽然MPSK和MDPSK等相移键控的带宽和功率方面都具有优势, 即带宽占用小和比特噪声比要求低。但是由图L6可见，在MPSK体制中，随 着(21)(22)(23)图L6 8PSK信号相位M的增大，相邻相位的距离逐渐减小，使噪声容限随之减小，误码率难于 保证。为了改善在M大时的噪声容限，发展出了 QAM体制。在QAM体制中, 信号的振幅和相

23、位作为两个独立的参量同时受到调制。这种信号的一个码元可 以表示为 (r) = Ak cos(6?(/ +。4) kT t (k + 1)T式中：k=整数；儿和a分别可以取多个离散值。式(2-1)可以展开为sk (r) = Ak cos 功 cos 卬-Ak sin 3k sin coj令Xk = AkcosSk, Yk = -AksinOk则式(2-1)变为sk (r) = Xk cos + Yk sin coQt和丫也是可以取多个离散的变量。从式(2-3)看出，可以看作是两个正交的振幅键控信号之和。在式(21)中，若8k值仅可以取n/4和迅/4, Ak值仅可以取+A和A,则此QAM信号就成为

24、QPSK信号，如图1.7所示:图1.7 4QAM信号矢量图所以，QPSK信号就是一种最简单的QAM信号。有代表性的QAM信号是 16进制的，记为16QAM,它的矢量图示于下图L8中：图1.8 16QAM信号矢量图图中用黑点表示每个码元的位置，并且示出它是由两个正交矢量合成的。 类似地，有64QAM和256QAM等QAM信号，如图1.9、图1.10所示。它们 总称为MQAM调制。由于从其矢量图看像是星座，故又称星座调制。图 L 9 R40AM图 1J 0 2560AM16QAM信号的产生方法主要有两种。第一种是正交调幅法，即用两路独立的正交4ASK信号叠加，形成16QAM信号，如图L11所示。图

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

26、特序列所调制。发送信号波形如图L13所示旭=AmcgT(t)cos27ifct + AmsgT (r)sinm = 1,2,加图1.13 M=16QAM信号星座图式中4,.和是电平集合,这些电平是通过将k比特序列映射为信号振 幅而获得的。例如一个16位正交幅度调制信号的星座图如下图所示，该星座是 通过用M=4PAM信号对每个正交载波进行振幅调制得到的。利用PAM分别调 制两个正交载波可得到矩形信号星座。QAM可以看成是振幅调制和相位调制的结合。因此发送的QAM信号波形 可表示为%加(f)= 4 (f)cos (2相I + 0n),m = 1,2，,根，n = 1,2，,%,如果2=21那么qa

27、m方法就可以达到以符号速率 %/（占+3）同时发送 勺+七=bg2%M个二进制数据。图1.14给出了 QAM 调制器的框图。图1.14 QAM调制器框图2、QAM的解调和判决假设在信号传输中存在载波相位偏移和加性高斯噪声。因此r可以表示为厂。）=4 g? ）cos（2阴.+。） +）sin（2次.+。） + 其中。是载波相位偏移，且n（t） = nc（t）cos2ct-ns（t）2ct将接收信号与下述两个相移函数进行相关=St （r）cos（2G +。），2（f） = g7）sin（2J +。）如图1.15所示，相关器的输出抽样后输入判决器。使用图115中所示的锁 相环估算接收信号的载波相位偏

28、移。，相移内和忆对该相位偏移进行补 偿。图1.15 QAM信号的解调和判决假设图中1.15所示的时钟与接收信号同步，以使相关器的输出在适当的时 刻及时被抽样。在这些条件下两个相关器的输出分别为八=AIJ1C + nc cos。一凡 sin。r2 = Amc + nc sin + nx cos。其中1 rL o T噪声分量是均值为0,方差为M/2的互不相关的高斯随机变量。最佳判决器计算距离量度（二邑“）=卜一。m =3、QAM的误码率性能矩形QAM信号星座最突出的优点就是容易产生PAM信号可直接加到两个 正交载波相位上，此外它们还便于解调。对于例=2下的矩形信号星座图（K为偶数）,QAM信号星座

29、图与正交载 波上的两个PAM信号是等价的，这两个信号中的每一个上都有屈=个信 号点。因为相位正交分量上的信号能被相干判决极好的分离，所以易于通过PAM 的误码率确定QAM的误码率。M进制QAM系统正确判决的概率是:p=1误码率的的严格上限。P” 抽样判决图1.17相干解调原理框图第二章MATLAB通信系统工具箱MATLAB通信系统工具箱主要应用于参数化模型，频谱分析和估计等。通 信系统工具箱实际上是用MATLAB的基本语句变成的各种子程序集，用于解决 某一方面的专门问题或实现某一类的新算法。MATLAB的通信系统工具箱可以 任意增减，不同的工具箱可以给不同领域的用户提供丰富强大的功能。MATL

30、AB中的通信系统工具箱是目前比较成熟的通信系统仿真工具， MATLAB通信系统工具箱提供的函数主要用于处理通信及传输问题，为我们分 析通信系统中每个过程提供了极大的帮助，使我们能够轻松直观的解决问题。常用的MATLAB函数主要有以下几项：1、figure功能：创建新的图形窗口(用于输出图形的窗口)。格式：figure说明：Hgure函数创建一个新的图形窗口，并成为当前图形窗口，所创建的 图形窗口的序号是按同一 MATLAB程序中创建的顺序号。2、plot功能：线型绘图函数。格式：plot(x)plot(x,y)说明：plot(x)是一种最简单的调用方式，x是长度为n的数值向量。plot(x)

31、的作用是在坐标系中顺序地用直接连接顶点i,x(i),i=l,2,n生成一条折线。当 向量元素充分多时，即可生成一条光滑的曲线。3、subplot功能：多坐标设置与定位当前坐标系。格式：subplot(m,n,k)说明：subplot(m,n,k)将图形窗口分成m行n列m*n块子区域，按从上到下，从左到右的顺序，在第k块子区域定义一个坐标系，使其成为当前坐标系，随 后的绘图函数将在该坐标系输出图形。4、 ezplot功能：该命令用来绘制符号表达式的自变量和对应各函数值的二维曲线。格式：ezplot(F, xmin,xmax,flg)说明：其中F是要画的符号函数；xmin,xmax是绘图的自变量范

32、围，tip 是窗口。5、axis功能：该命令用来控制坐标轴的特性。格式：axis(xmin,xmax, ymin,ymax)说明：此为坐标范围，其中xminxmax , ymlnymax。7、angle功能：求复数的相角。格式：P=angle(Z)说明：当Z为复数矩阵时，用来求矩阵Z中每个元素的相角，相角位于 1严。第三章MATLAB的模拟调制和解调实例第一节用MATLAB分析双边带幅度调制(DSM-AM)我们可以利用MATLAB强大的符号运算功能来进行运算,再根据MATLAB的 可视化结果进行分析。在DSB-AM中，己调信号的时域表示为：u (t) =m (t) c (t) =Acm (t)

33、 cosQR fa + /)式中，In(t)是消息信号，c(t)=Ac cos(2n/d +戊)为载波，fc是载波的频 率(单位：HZ),火是初始相位。为了讨论方便取初相化=0 (以下类似)。随u(t)作傅里叶变换，即可得到信号的频域表示：U (f) =Ac/2 M(/-/c)+Ac/2 M(f + fc)传输带宽Bt是消息信号带宽W的两倍，即：Bt=2W0tt0/3Y某消息信号m(t)= -2 t0/3t zeros (n-n2);df=fs/n;运行后得到的信号和调制信号如图3. 1所示：图3.1双边带幅度调制信号第二节用MATLAB分析相干解调过程0ttO/3例某消息信号，(t)= -2

34、 t0/3t 1等概分布的随机信号。源代码如下所示：random_binary.m%产生二进制信源随机序列function info=random_binary(N)if nargin = 0,%如果没有输入参数，则指定信息序列为10000个码元N=10000；end;for 1=1 :N,temp=rand;if (temp0.5),info(i)=0;% 1/2 的概率输出为0elseinfo(i)=l;% 1/2 的概率输出为1endend;三、星座图映射对产生的二进制随机序列进行串并转换，分离出I分量、Q分量，然后再分 别进行电平映射。由于是调用matlab系统函数调制解调，在此将转换

35、后边的序 列进行四进制转换，方便后面的调制，再将转换好的序列通过调用qam()函数进 行16qam调制，具体代码如下：代码如下：%串/并变换分离出I分量、Q分量，然后再分别进行电平映射I=x(l:2:nn-1);I,InJ=two2four(I,4*m)；Q=x(2:2:nn);Q,Qn=two2four(Q,4*m);if Kbase=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);% 调制二进制转换成四进制代码：two2foiir.m%二进制转换成四进制func

36、tion y,yn=two2four(x,m);T=0 1;3 2;n=length(x);ii=l;for i=l:2:n-l;xl=x(i:i+l)+l;yn(ii)=T(xi(l),xi(2);li=ii+l;end;yn=yn-L5;y=yn；for i=l:m-l;y=y;yn；end;y=y(：); %映射电平分别为L5； 0.5； 0.5； 1.5画出星座图代码如下：consteLm%画出星座图function c=constel(x,fs,fb,fc);N=length(x);m=2*fs/fb; n=fs/fc;ll=m-n;i=l；ph0=(il-l)*2*pi/n;whi

37、le i = N/in;xl=x(il:il+n.l);y=2*fft(xi)/n;c(i)=y(2);i=i+l;il=il+m;end;%如果无输出，则作图if nargoutl;cmax=max(abs(c);ph=(0:5:360)*pi/180;plot(L414*cos(ph),L414*sln(ph),C);hold on;for i=l:length(c);ph=phO-angle(c(i);a=abs(c(i)/cmax* 1.414;plot(a*cos(ph),a*sln(ph),r*);end;plot(-1.5 1.5,0 OJ/kfJO 0,-1.5 1.5,k:)

38、;hold off;axis equal;axis(-L5 1.5 -1.5 1.5);end;四、波形成形(平方根升余弦滤波器)为了避免相邻传输信号之间的串扰，多元符号需要有合适的信号波形。方 波是在本地数字信号处理时常见的波形，但在实际传输时这种方波并不合适。 根据奈奎斯特第一准则，在实际通信系统中一般均使接收波形为升余弦滚降信 号。这一过程由发送端的基带成形滤波器和接收端的匹配滤波器两个环节共同 实现，因此每个环节均为平方根升余弦滚降灌波，两个环节合成就实现了一个 升余弦滚降滤波。实现平方根升余弦滚降信号的过程称为“波形成形”，通过 采用合适的滤波器对多元码流进行滤波实现，由于生成的是基

39、带信号，因此这一过程又称“基带成形滤波”。代码如下：bshape.m%基带升余弦成形滤波器function y=bshape(x,fs,fb,N,alfa,delay);%设置默认参数if nargin6; delay=8; end;if nargin5; alfa=0.5; end;if nargin4; N=16; end;b=firrcos(N,fb,2*alfa*fbfs);y=fllter(b,l,x);五、调制调制代码：qam.mfunction yJ,Q=qam(x,Kbase,fs,fb,fc);T=length(x)/fb;m=fs/fb;nn=length(x);dt=l/

40、fs;t=O:dt:T-dt;%串/并变换分离出I分量、Q分量，然后再分别进行电平映射I=x(l:2:nn-1);IJn=two2four(M*m);Q=x(2:2:nn);Q,Qn=two2four(Q,4*m);if Kbase=2;%基带成形漉波I=bsliape(I,fs,fb/4);Q=bshape(Q,fs,fb/4);end;y=I.*cos(2*pi*fc*t)-Q.*sin(2*pi*fc*t);%调制六、加入高斯白噪声之后解调为了简化程序和得到可靠的误码率，我们在解调时并未从已调信号中恢复 载波，而是直接产生与调制时一模一样的载波来进行信号解调。加入不同强度的高斯白噪声代码

41、：SNR_in_dB=8:2:24;%AWGN 信道信噪比for j=l:length(SNR_ln_dB)y_add_noise=awgn(y2,SNR_in_dB(j);% 力口入不同强度的高斯白噪声y_output=qamdet(y_add_iioise,fs5fb,fc);% 对已调信号进行解调解调时先设计一个巴特沃斯滤波器，然后将I分量、Q分量并/串转换，最终恢复成码元序列XI】，然后进行解调。 解调的代码如下：%QAM信号解调function (xn,x=qamclet(y,fsfb,fc);dt=l/fs;t=O:dt:Qength(y)-l)*dt;I=y.*cos(2*pi*

42、fc*t);Q=-y.*sin(2*pl*fc*t);b,a=butter(2,2*fb/fs);I=flltfilt(b,aJ);Q=flltfllt(b,a,Q);m=4*fs/fb;N=length(y)/m;n=(.6:l:N)*m;%设计巴特沃斯滤波器n=flx(n);1In=I(n);Qn=Q(n);xn=four2two(In Qn);%I分量Q分量并/串转换，最终恢复成码元序列XIInn=length(xn);xn=xn(l:nn/2);xn(nn/2+l:nn);xn=xn(:);xn=xnf;四进制转换成二进制代码如下：%四进制转换成二进制function xn=four2

43、two(yn);y=yn; ymin=inin(y);ymax=max(y);yinax=max(ymax abs(ymin);ymin=-abs(ymax);yn=(y-yinin)*3/(ymax-ymln);%设置门限电平，判决I0=flnd(yn=0.5 & yn=L5 & yn=2.5);yn(I3)=ones(slze(I3)*3;%一位四进制码元转换为两位二进制码元T=0 0;0 1;1 1;1 0;n=length(yn);for i=l:n;xn(i,:)=T(yn(i)+h：);end;xn=xn,;ixn=xn(:); xn=xnf;七、误码率曲线1误码率代码如下；num

44、oferr=0;for i=l:Nif (y_output =lnfo )，numoferr=numoferr+1;end;end;Pe(J )=numoferr/N;% 统计误 码率end;figure;semilogy (SNR_in_d B ,Pe； red *)；grid on;xlabel(*SNR in (IB1);ylabeUTe*);调制在不同信道噪声强度下的误码率，)；八、16QAM载波调制信号在AWGN信道下的性能16-QAM载波调制信号在AWGN信道下的性能代码如下： clear allnsymbol= 100000;%每种信噪比下的发送符号数M=16;% 16-QAMg

45、raycode=0 13 2 4 5 7 6 12 13 15 14 8 9 11 10;% Gray编码规则EsN0=5:20;%信噪比，Es/N0snrl=10eA(EsN0/10);%信噪比转换为线性值msg=randint(lnsymbol,M);% 消息数据msgl=graycode(msg+1);% Gray 映射insginod=qainmod(msgl9M);% 基带16-QAM调制spow=norm(msgmod).A2/nsymbol;% 求每个符号的平均功率for indx=l:length(EsNO)sigma=sqrt(spow/(2*snrl(indx);%根据符号功率求噪声功率rx=insgmod+sigina*(raiidn(lUength(msgmod)+J*rancln(lJength(msgmod);y=qamdemocI(rx,M);decinsg=graycode(y+l);err,ber(lndx)=biterr(insg,decmsgUog2(M);% 误比特率err,ser(indx)=symerr(insg,decmsg);% 误符号率endP4=2