通信系统毕业论文-基于MATLAB的数字基带传输系统及数字调制的仿真

1、(此文档为word格式，下载后您可任意编辑修改！) 课 程 设 计 说 明 书 设 计 及 说 明主要结果目 录前言 2MATLAB软件简介 2数字基带传输系统仿真 51.1数字基带传输系统简介 51.2 数字基带信号 71.2.1 数字基带信号的要求 71.2.2 数字基带信号81.2.3 常用的基带传输码型912系统模型12传输系统模型12131.4.1 余弦滚降基带传输系统131.4.2 眼图 1417数字调制系统仿真前言17移键控2FSK 172FSK信号的产生172FSK解调原理192FSK调制与解调基于MATLAB仿真 20幅键控(2ASK232ASK调制与解调基于MATLAB仿真

2、25移键控(2PSK272PSK原理272PSK基于MATLAB的仿真29移键控(4PSK334PSK信号的产生与解调334PSK基于MATLAB的仿真34正交幅度调制(16QAM42正交调制及相干解调原理框图4216QAM调制与解调基于MATLAB的仿真43 48参考文献：49前言 随着通信系统的规模和复杂度不断增加， 统的设计方法已经不能适应开展传的需要，通信系统的模拟仿真技术越来越受到重视。传统的通信仿真技术主要分 可以 得到与真实环境十分接近的结果，为手工分析与电路试验2种， 但耗时长方法比较繁杂，而通信系统的计算机模拟仿真技术是介于上述2种方法的一种系统设计方法，它可以让用户 在很

3、短的时间内建立整个通信系统模型，并对其进行模拟仿真。通信原理计算机仿真实验，是对数字基带传输系统的仿真。仿真工具是MATLAB程序设计语 言。 MATLAB是一种先进的高技术程序设计语言，主要用于数值计算及可视化图形处理。特点 是将数值分析、矩阵计算、图形、图像处理和仿 真等诸多强大功能集成在一个极易使用的交互式环境中伪科学研究、工程设计以及必须进行有效数值计算的众多学科提供了一种高效率的编程 工具。运用MATLAB，可以对数字基带传输系统 进行较为全面地研究。为了使本科类学生学好通信课程，我们进行了试点，通过课程设计的方式针对通信原理的很多内容进行了仿真。MATLAB软件简介 MATLAB是

4、由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言如C、Fortran的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创立用户界面、连接其

5、他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的根本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也参加了对C，FORTRAN，C+ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。

6、MATLAB 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 技术 数字信号处理 技术 通讯系统设计与仿真 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱单独提供的专用 MATLAB 函数集扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。 20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担，用FORTRAN编写了最早的MATLAB。1984年由Little、Moler、Steve

7、 Bangert合作成立了的MathWorks公司正式把MATLAB推向市场。到20世纪90年代，MATLAB已成为国际控制界的标准计算软件。综上所述，Matlab 语言有如下特点：编程效率高 它是一种面向科学与工程计算的高级语言，允许用数学形式的语言编写程序，且 比 Basic、 Fortran 和 C 等语言更加接近我们书写计算公式的思维方式， Matlab 用 编写程序犹如在演算纸上排列出公式与求解问题。因此，Matlab 语言也可通俗 地称为演算纸式科学算法语言由于它编写简单，所以编程效率高，易学易懂。 2用户使用方便 Matlab 语言是一种解释执行的语言在没被专门的工具编译之前，它

8、灵活、 方便，其调试程序手段丰富，调试速度快，需要学习时间少。人们用任何一种语 言编写程序和调试程序一般都要经过四个步骤：编辑、编译、连接以及执行和调 试。各个步骤之间是顺序关系，编程的过程就是在它们之间作瀑布型的循环。 Matlab 语言与其它语言相比，较好地解决了上述问题，把编辑、编译、连接和 执行融为一体。它能在同一画面上进行灵活操作快速排除输入程序中的书写错 误、语法错误以至语意错误，从而加快了用户编写、修改和调试程序的速度，可 以说在编程和调试过程中它是一种比 VB 还要简单的语言。 具体地说，Matlab 运行时，如直接在命令行输入 Mailab 语句命令，包括调 用 M 文件的语

9、句，每输入一条语句，就立即对其进行处理，完成绩译、连接和运 行的全过程。又如，将 Matlab 源程序编辑为 M 文件，由于 Mat1ab 磁盘文件也是 M 文件，所以编辑后的源文件就可直接运行，而不需进行编译和连接。在运行 M 文件时， 如果有错， 计算机屏幕上会给出详细的出锗信息， 用户经修改后再执行， 直到正确为止。所以可以说，Mat1ab 语言不仅是一种语言，广义上讲是一种该 语言开发系统，即语言调试系统。 3扩充能力强 Matlab 语言有丰富的库函数，在进行复杂的数学运算时可以直接调用，而且 Matlab 的库函数同用户文件在形成上一样， 所以用户文件也可作为 Matlab 的库

10、函数来调用。因而，用户可以根据自己的需要方便地建立和扩充新的库函数，以 便提高 Matlab 使用效率和扩充它的功能。另外，为了充分利用 Fortran、C 等语 言的资源， 包括用户已编好的 Fortran， 语言程序， C 通过建立 Me 调文件的形式， 混合编程，方便地调用有关的 Fortran，C 语言的子程序。 4语句简单，内涵丰富 Mat1ab 语言中最基本最重要的成分是函数，其一般形式为a，6，c = fun d，e，f，即一个函数由函数名，输入变量 d，e，f,和输出变量 a，b，c组成，同一函数名 F，不同数目的输入变量包括无输入变量及 不同数目的输出变量，代表着不同的含义有

11、点像面向对象中的多态性。这不仅 使 Matlab 的库函数功能更丰富，而大大减少了需要的磁盘空间，使得 Matlab 编写的 M 文件简单、短小而高效。 5高效方便的矩阵和数组运算 Matlab 语言象 Basic、Fortran 和 C 语言一样规定了矩阵的算术运算符、关系运 算符、逻辑运算符、条件运算符及赋值运算符，而且这些运算符大局部可以毫无 改变地照搬到数组间的运算，有些如算术运算符只要增加“就可用于数组间 的运算，另外，它不需定义数组的维数，并给出矩阵函数、特殊矩阵专门的库函 数，使之在求解诸如信号处理、建模、系统识别、控制、优化等领域的问题时， 显得大为简捷、高效、方便，这是其它高

12、级语言所不能比较的。在此根底上，高 版本的 Matlab 已逐步扩展到科学及工程计算的其它领域。因此，不久的将来， 它一定能名符其实地成为“万能演算纸式的科学算法语言。 6方便的绘图功能 Matlab 的绘图是十分方便的，它有一系列绘图函数命令，例如线性坐标、 对数坐标，半对数坐标及极坐标，均只需调用不同的绘图函数命令，在图上 标出图题、XY 轴标注，格栅绘制也只需调用相应的命令，简单易行。另外， 在调用绘图函数时调整自变量可绘出不变颜色的点、线、复线或多重线。这种为 科学研究着想的设计是通用的编程语言所不及的。总之，Matlab 语言的设计思 想可以说代表了当前计算机高级语言的开展方向。 数

13、字基带传输系统仿真1.1数字基带传输系统 数字基带传输系统的介绍 在数字传输系统中，其传输的对象通常是二进制数字信号，它可能是来自计算机、电传打字机或其它数字设备的各种数字脉冲，也可能是来自数字 终端的脉冲编码调制PCM信号。这些二进制数字信号的频带范围通常从直流和低频开始，直到某一频率m f ，我们称这种信号为数字基带信号。在某些有线信道中，特别是在传输距离不太远的情况下，数字基带信号可以不经过调制和解调过程在信道中直接传送，这种不使用调制和解调设备而直接传输基带信号的通信系统，我们称它为基带传输系统。而在另外一些信道，特别是无线信道和光信道中，数字基带信号那么必须经过调制过程，将信号频谱搬

14、移到高频处才能在信道中传输，相应地，在接收端必须经过解调过程，才能恢复数字基带信号。我们把这种包括了调制和解调过程的传输系统称为数字载波传输系统。数字基带传输系统的模型如图 1-1 所示，它主要包括码型变换器、发送滤波器、信道、接收滤波器、均衡器和取样判决器等局部。图1-1 数字基带传输系统模型数字 基带传输系统 的输入信号是由终端设备或编码设备产生的二进制脉 冲序列，通常是单极 性的矩形脉冲信号NRZ 码。为了使这种信号适合 于信道的传输，一般 要经过码形变换器，把单极性的二进制脉冲变成双极 性脉冲如 AMI 码或 3 HDB 码。发送滤波器对码脉冲进行波形变换，以减 小信号在基带 传输系

15、统中传输时产生的码间串扰。信号在传输过程中，由于 信道特性不理想 及加性噪声的影响，会使接收到的信号波形产生失真 ，为了减小失 真对信 号的影响，接收信号首先进入接收滤波器滤波，然后 再经均衡器对 失真信 号进行校正，最后由取样判决器恢复数字基带脉冲序 列。 目前，虽 然在实际使用的数字通信系统中，基带传输方式不如数字载波传 输方式那样应 用广泛 ，但由于数字基带传输系统是数字通信系统中最根本 的传输方式，而且从 理论上来说，任何一种线性载波传输系统都可以等效 为基带传输系 统，因 此理解数字信号的基带传输过程十分重要。 数字基带信号有二元码和三元码，有归零码和非归零码等，有的具有直流分量，在

16、波形上具有不同的特点，他们有不同的特点，有的低频成份多，有的高频成份多，有的具有直流分量，有的占有带宽等，所有这些在波形处理时会对一些学生产生模糊的概念，针对本科类的学生要求，他们如何理解、区分、掌握这些信号波形的特点，同时可以让学生在仿真过程中对通信原理的各种概念加深理解。另外，此仿真实验只需在计算机的虚拟实验室即可，不受实验场地、环境的限制。软件的功能主要有： 实现各种常用码型的数字基带信号仿真；能产生随机的数字信号序列，具有普遍性； 3)能绘制直观、清晰、准确、可靠的数字基带 信号仿真图形；4)要对相应的码型的特点进行相应的描述。图 1-2 仿真结构图 在通信中，数字基带信号有多种码型表

17、示，它们在传输过程中有随机性，为 用让这种波形描述具有普遍性， m 序列伪随机码来作为码型的仿真数字序列。利用 MATLAB 软件仿真出每一种码型，让学生通过仿真软件的使用，加深对码和波形的理解。1.2 数字基带信号1.2.1 数字基带信号的要求不同形式的数字基带信号又称为码型具有不同的频谱结构，为适应信道的传输特性及接收端再生、恢复数字基带信号的需要，必须合理地设计数字基带信号，即选择适宜的信号码型。适合于在有线信道中传输的数字基带信号形式称为线路传输码型。一般来说，选择数字基带信号码型时，应遵循以下根本原那么：数字基带信号应不含有直流分量，且低频及高频分量也应尽量的少。在基带传输系统中，往

18、往存在着隔直电容及耦合变压器，不利于直流及低频分量的传输。此外，高频分量的衰减随传输距离的增加会快速地增大，另一方面，过多的高频分量还会引起话路之间的串扰，因此希望数字基带信号中的高频分量也要尽量的少。数字基带信号中应含有足够大的定时信息分量。基带传输系统在接收端进行取样、判决、再生原始数字基带信号时，必须有取样定时脉冲。一般来说，这种定时脉冲信号是从数字基带信号中直接提取的。这就要求数字基带信号中含有或经过简单处理后含有定时脉冲信号的线谱分量，以便同步电路提取。实际经验告诉我们，所传输的信号中不仅要有定时分量，而且定时分量还必须具有足够大的能量，才能保证同步提取电路稳定可靠的工作。基带传输的

19、信号码型应对任何信源具有透明性，即与信源的统计特性无关。这一点也是为了便于定时信息的提取而提出的。信源的编码序列中，有时候会出现长时间连“0的情况，这使接收端在较长的时间段内无信号，因而同步提取电路无法工作。为防止出现这种现象，基带传输码型必须保证在任何情况下都能使序列中“1和“0出现的概率根本相同，且不出现长连“1或“0的情况。当然，这要通过码型变换过程来实现。码型变换实际上是把数字信息用电脉冲信号重新表示的过程。此外，选择的基带传输信号码型还应有利于提高系统的传输效率；具有较强的抗噪声和码间串扰的能力及自检能力。实际系统中常常根据通信距离和传输方式等不同的要求，选择适宜的基带码型。1.2.

20、2 数字基带信号 对不同的数字基带传输系统，应根据不同的信道特性及系统指标要求，选择不同的数字脉冲波形。原那么上可选择任意形状的脉冲作为基带信号波形，如矩形脉冲、三角波、高斯脉冲及升余弦脉冲等。但实际系统常用的数字波形是矩形脉冲，这是由于矩形脉冲易于产生和处理。下面我们就以矩形脉冲为例，介绍常用的几种数字基带信号波形。1.单极性波形NRZ 这是一种最简单的二进制数字基带信号波形。这种波形用正或负电平和零电平分别表示二进制码元的“1码和“0码，也就是用脉冲的有无来表示码元的“1和“0，这种波形的特点是脉冲的极性单一，有直流分量，且脉冲之间无空隙，即脉冲的宽度等于码元宽度。故这种脉冲又称为不归零码

21、NRZ -NonReturn to ZeroNRZ波形一般用于近距离的电传机之间的信号传输。双极性波形 在双极性波形中，用正电平和负电平分别表示二进制码元的“1码和“0码，这种波形的脉冲之间也无空隙。此外，从信源的统计规律来看，“1码和“0码出现的概率相等，所以这种波形无直流分量。同时这种波形具有较强的抗干扰能力。故双极性波形在基带传输系统中应用广泛。单极性归零波形RZ 这种波形的特点是脉冲的宽度 小于码元的宽度T ，每个电脉冲在小于码元宽度的时间内总要回到零电平，故这种波形又称为归零波RZ-Return to Zero。归零波形由于码元间隔明显，因此有利于定时信息的提取。但单极性RZ波形中仍

22、含有直流分量，且由于脉冲变窄，码元能量减小，因而在匹配接收时，输出信噪比较不归零波形的低。双极性归零波形 这种波形是用正电平和负电平分别表示二进制码元的“1码和“0码，但每个电脉冲在小于码元宽度的时间内都要回到零电平，这种波形兼有双极性波形和归零波形的特点。 差分波形相对码波形 信息码元与脉冲电平之间的对应关系是固定不变的绝对的，故称这些波形为绝对码波形，信息码也称为绝对码。所谓差分波形是一种把信息码元“1和“0反映在相邻信号码元的相对电平变化上的波形，差分波形中，码元“1和“0分别用电平的跳变和不变来表示，即用相邻信号码元的相对电平来表示码元“1和“0，故差分波形也称为相对码波形。差分波形也

23、可以看成是差分码序列bn 对应的绝对码波形，差分码bn 与绝对码an之间的关系可用以下的编码方程表示 bnbn1 an 1.1式中，为模2和运算符号。 由上式看出，当绝对码an 每出现一个“1码时，差分码bn电平变化一次；当 an 出现“0码时，差分码bn 电平与前一码元bn-1 相同。可见， bn 前后码元取值的变化代表了原信码n a 中的“1和“0。由式1.1可以导出译码方程为 an bn-1bn 1.2 由上式可看出，译码时只要检查前后码元电平是否有变化就可以判决发送的是“1码还是“0码。多电平脉冲波形多进制波形 上述各种波形都是二进制波形，实际上还存在多电平脉冲波形，也称为多进制波形。

24、这种波形的取值不是两值而是多值的。例如，代表四种状态的四电平脉冲波形，每种电平可用两位二进制码元来表示，如00 代表-3E， 代表-E， 代 01 10表E，11 代表3E，这种波形一般在高速数据传输系统中用来压缩码元速率，提高系统的频带利用率。但在相同信号功率的条件下，多进制传输系统的抗干扰性能不如二进制系统。1.2.3 常用的基带传输码型 前面提到，为满足基带传输系统的特性要求，必须选择合适的传输码型。基带传输系统中常用的线路传输型码主要有：传号交替反转码-AMI 码、三阶高密度双极性码- 3 HDB码、分相码-Manchester 码、传号反转码-CMI 码以及4B3T 码等。下面我们详

25、细地介绍这些码型。1、传号交替反转码-AMI 码 AMI Alternate Mark Inversion码又称为平衡对称码。这种码的编码规那么是：把码元序列中的“1码变为极性交替变化的传输码1、-1、1、-1、，而码元序列中的“0码保持不变。例如： 码元序列： 1 00 1 1 0 1 0 1 1 1 1 00 由AMI 码的编码规那么可以看出，由于1和-1各占一半，因此，这种码中无直流分量，且其低频和高频分量也较少，信号的能量主要集中在2 T f 处，其中Tf 为码元速率。此外，AMI 码编码过程中，将一个二进制符号变成了一个三进制符号，即这种码脉冲有三种电平，因此我们把这种码称为伪三电平

26、码，也称为1B1T 码型。AMI码除了上述特点外，还有编译码电路简单及便于观察误码情况等优点。但是AMI码有一个重要的缺陷，就是当码元序列中出现长连“0时，会造成提取定时信号的困难，因而实际系统中常采用AMI 码的改良型 HDB3码。HDB3 码HDB3 High Density Bipolar 3是三阶高密度双极性码，它是为了克服传输波形中出现长连“0码情况而设计的AMI 码的改良型。HDB3 码的编码规那么是：1把码元序列进行AMI 编码，然后去检查AMI 码中连0 的个数，如果没有四个以上包括四个连0 串时，那么这时的AMI 码就是3 HDB 码。2如果出现四个以上连0 串时，那么将每4

27、 个连0 小段的第4 个0变成与其前一个非0 码1 或-1相同的码。显然，这个码破坏了“极性交替反转的规那么，因而称其为破坏码，用符号V 表示即1 记为V， 记为-V -1 。3为了使附加V 码后的序列中仍不含直流分量，必须保证相邻的V 码极性交替。这一点，当相邻的V 码之间有奇数个非0 码时，是能得到保证的；但当相邻的V 码之间有偶数个非0 码时，那么得不到保证。这时再将该连0 小段中的第1 个0 变成B 或-B，B 的极性与其前一个非0 码相反，并让后面的非零码从V 码后开始再极性交替变化。例如： 码元序列： 1 0000 1 0 1 0 0 0 0 1 000 0 1 1 AMI 码：

28、1 0000 -1 0 1 0 0 0 0 1 000 0 11 HDB3 码： 1 000V -1 0 1 -B00-V 1 000V -1 1 上例中，第1个V码和第2个V码之间，有2个非0 码偶数，故将第2个4 连0小段中的第1个0变成-B；第2个V码和第3个V码之间，有1个非0码奇数，不需变化。最后可看出， HDB3 码中，V码与其前一个非0码1 或-1极性相同，起破坏作用；相邻的V码极性交替；除V码外，包括B码在内的所有非0码极性交替。虽然HDB3 码的编码规那么比较复杂，但译码却比较简单。从编码过程中可以看出，每一个V码总是与其前一个非0码包括B码在内同极性，因此从收到的码序列中可

29、以很容易地找到破坏点V码，于是可断定V码及其前3个码都为0码，再将所有的-1变为1后，便可恢复原始信息代码。HDB3码的特点是明显的，它既保存AMI码无直流分量，便于直接传输的优点，又克服了长连0串连0的个数最多3个的出现，HDB3 码的频谱中既消除了直流和甚低频分量，又消除了方波中的高频分量，非常适合基带传输系统的特性要求。因此，HDB3码是目前实际系统中应用最广泛的码型。虽然HDB3码比AMI 码的性能更好，但它仍属于1B1T 码型。曼彻斯特Manchester码曼彻斯特码又称数字双相码或分相码，曼彻斯特码用一个周期的方波来代表码元“1，而用它的反相波形来代表码元“0。这种码在每个码元的中

30、心部位都发生电平跳变，因此有利于定时同步信号的提取，而且定时分量的大小不受信源统计特性的影响。曼彻斯特码中，由于正负脉冲各占一半，因此无直流分量，但这种码占用的频带增加了一倍。曼彻斯特码适合在较短距离的同轴电缆信道上传输。CMI 码 CMI 码称为传号反转码。在CMI 码中，“1码传号交替地用正、负电平脉冲来表示，而“0码那么用固定相位的一个周期方波表示，CMI 码和曼彻斯特码相似，不含有直流分量，且易于提取同步信号。CMI 码的另一个特点是具有一定的误码检测能力。这是因为，CMI 码中的“1码相当于用交替的“00和“11两位码组表示，而“0码那么固定地用“01 码组表示。正常情况下，序列中不

31、会出现“10码组，且“00和“11码组连续出现的情况也不会发生，这种相关性可以用来检测因干扰而产生的局部错码。根据原CCITT 的建议，CMI码可用作脉冲编码调制四次群的接口码型以及速率低于8448 kb s的光纤数字传输系统中的线路传输码型。此外，CMI 码和曼彻斯特码一样都是将一位二进制码用一组两位二进制码表示，因此称其为1B2B 码。 54 B3T 码4B3T 码是1B1T 码的改良型它把4 个二进制码元变换为3 个三进制码元。显然，在相同信息速率的条件下，4B3T 码的码元传输速率要比1B1T 码的低，因而提高了系统的传输效率。 4B3T 码的变换过程中需要同步信号，变换电路比较复杂，

32、故一般较少采用。1.系统模型传输系统模型1.4.1 余弦滚降基带传输系统升余弦滚降传输特性H()可表示为 H()是对截止频率b的理想低通特性H0()按H()的滚降特性进行“圆滑得到的，H1()对于b具有奇对称的幅度特性，其上、下截止角频率分别为b+1、b-1。它的选取可根据需要选择，升余弦滚降传输特性H1()采用余弦函数， 此时H()为 称为滚降系数。余弦滚降系统基于matlab仿真源程序如下:% 数字基带信号传输 码间串扰 升余弦滚降系统的频谱及其时域波形Ts=1;N=17;dt=TsN;df=1.0(20.0*Ts);t=-10*Ts:dt:10*Ts;f=-2Ts:df:2Ts;a=0,

33、0.5,1;for n=1:length(a) for k=1:length(f) if abs(f(k)0.5*(1+a(n)Ts Xf(n,k)=0; elseif abs(f(k)=st2(i) st(i)=0; else st(i)=st2(i); endend 仿真结果1仿真结果2 仿真结果3 仿真结果4 仿真结果5 幅键控(2ASK数字幅度调制又称幅度键控ASK，二进制幅度键控记作2ASK。2ASK是利用代表数字信息“0”或“1”的基带矩形脉冲去键控一个连续的载波，使载波时断时续地输出。有载波输出时表示发送“1”，无载波输出时表示发送“0”。借助于第3章幅度调制的原理，2ASK信号

34、可表示为 8-1式中，为载波角频率，为单极性NRZ矩形脉冲序列 8-2其中，是持续时间为、高度为1的矩形脉冲，常称为门函数；为二进制数字 8-32ASK信号的产生方法调制方法有两种，如图8-1所示。图a是一般的模拟幅度调制方法，不过这里的由式8-2规定；图b是一种键控方法，这里的开关电路受控制。图c给出了及的波形例如。二进制幅度键控信号，由于一个信号状态始终为0，相当于处于断开状态，故又常称为通断键控信号OOK信号。图8-1 2ASK产生原理框图和示意波形2ASK信号解调的常用方法主要有两种：包络检波法和相干检测法。包络检波法的原理方框图如图8-2所示。带通滤波器BPF恰好使2ASK信号完整地

35、通过，经包络检测后，输出其包络。低通滤波器LPF的作用是滤除高频杂波，使基带信号包络通过。抽样判决器包括抽样、判决及码元形成器。定时抽样脉冲位同步信号是很窄的脉冲，通常位于每个码元的中央位置，其重复周期等于码元的宽度。不计噪声影响时，带通滤波器输出为2ASK信号，即，包络检波器输出为。经抽样、判决后将码元再生，即可恢复出数字序列。图8-2 2ASK的包络非相干解调原理框图相干检测法原理方框图如图8-3所示。相干检测就是同步解调，要求接收机产生一个与发送载波同频同相的本地载波信号，称其为同步载波或相干载波。利用此载波与收到的已调信号相乘，输出为图8-3 2ASK的相干解调原理框图经低通滤波滤除第

36、二项高频分量后，即可输出信号。低通滤波器的截止频率与基带数字信号的最高频率相等。由于噪声影响及传输特性的不理想，低通滤波器输出波形有失真，经抽样判决、整形后再生数字基带脉冲。虽然2ASK信号中确实存在着载波分量，原那么上可以通过窄带滤波器或锁相环来提取同步载波，但这会给接收设备增加复杂性。因此，实际中很少采用相干解调法来解调2ASK信号。2ASK调制与解调基于MATLAB仿真：2ASK基于MATLAB的程序代码：function ask1()sim(dm);t1=ScopeData5(:,1);y1=ScopeData5(:,2);t2=ScopeData2(:,1);y2=ScopeData

37、2(:,2);t3=ScopeData(:,1);y3=ScopeData(:,2);t5=ScopeData1(:,1);y5=ScopeData1(:,2);t6=ScopeData7(:,1);y6=ScopeData7(:,2);t7=ScopeData3(:,1);y7=ScopeData3(:,2); %引入示波器的数据subplot(3,2,5);plot(t1,y1),grid;xlabel(滤波后的波形)subplot(3,2,4);plot(t2,y2),grid;xlabel(相干解调后的波形)subplot(3,2,1);stairs(t3,y3),axis(0 20

38、-0.5 2.5),grid;title( ASK调制与解调各环节波形)xlabel(源信号)subplot(3,2,2);plot(t5,y5),grid;xlabel(调制后的波形)subplot(3,2,3);plot(t6,y6),grid;xlabel(经过高斯信道后的波形)subplot(3,2,6);stairs(t7,y7),axis(0 20 -0.5 2.5),grid;xlabel(抽样判决后的波形)figure(2)n1=1024;f1=(-n12:1:n12-1);s=fft(y3,n1);subplot(3,2,1);plot(f1,abs(fftshift(s),

39、%axis(-30 30 0 7)xlabel(源信号频谱)n2=1024;f2=(-n22:1:n22-1);s=fft(y5,n2);subplot(3,2,2);plot(f2,abs(fftshift(s),%axis(-400 400 0 300)xlabel(调制后的频谱)s=fft(y2,n2);subplot(3,2,3);plot(f2,abs(fftshift(s),%axis(-400 400 0 300)xlabel(相干解调后的频谱)s=fft(y1,n2);subplot(3,2,4);plot(f2,abs(fftshift(s),%axis(-400 400 0

40、 300)xlabel(滤波后的频谱)s=fft(y7,n2);subplot(3,2,5);plot(f2,abs(fftshift(s),%axis(-400 400 0 300)xlabel(抽样判决后的频谱)仿真结果：移键控(2PSK2PSK原理2PSK的一般原理及实现方法绝对相移是利用载波的相位指初相直接表示数字信号的相移方式。二进制相移键控中，通常用相位0和来分别表示“0或“1。2PSK已调信号的时域表达式为 8-20这里，与2ASK及2FSK时不同，为双极性数字基带信号，即 8-21式中，是高度为1，宽度为的门函数； 8-22因此，在某一个码元持续时间内观察时，有，或8-23当码

41、元宽度为载波周期的整数倍时，2PSK信号的典型波形如图8-11所示图8-112PSK信号的典型波形2PSK信号的调制方框图如图8-12示。图a是产生2PSK信号的模拟调制法框图；图b是产生2PSK信号的键控法框图。图8-122PSK调制器框图就模拟调制法而言，与产生2ASK信号的方法比较，只是对要求不同，因此2PSK信号可以看作是双极性基带信号作用下的DSB调幅信号。而就键控法来说，用数字基带信号控制开关电路，选择不同相位的载波输出，这时为单极性NRZ或双极性NRZ脉冲序列信号均可。2PSK信号属于DSB信号，它的解调，不再能采用包络检测的方法，只能进行相干解调，其方框图如图8-13。工作原理

42、简要分析如下。图8-13 2PSK信号接收系统方框图不考虑噪声时，带通滤波器输出可表示为 8-24式中为2PSK信号某一码元的初相。时，代表数字“0；时，代表数字“1。与同步载波相乘后，输出为8-25经低通滤波器滤除高频分量，得解调器输出为 8-26根据发端产生2PSK信号时0或代表数字信息“1或“0的规定，以及收端与的关系的特性，抽样判决器的判决准那么为 8-28其中为在抽样时刻的值。2PSK接收系统各点波形如图8-14所示。图8-14 2PSK解调各点波形可见，2PSK信号相干解调的过程实际上是输入已调信号与本地载波信号进行极性比较的过程，故常称为极性比较法解调。由于2PSK信号实际上是以

43、一个固定初相的末调载波为参考的，因此，解调时必须有与此同频同相的同步载波。如果同步载波的相位发生变化，如0相位变为相位或相位变为0相位，那么恢复的数字信息就会发生“0变“1或“1变“0，从而造成错误的恢复。这种因为本地参考载波倒相，而在接收端发生错误恢复的现象称为“倒现象或“反向工作现象。绝对移相的主要缺点是容易产生相位模糊，造成反向工作。2PSK基于MATLAB的仿真：2PSK基于MATLAB的程序代码：clear all; close all; fs=8e5;%抽样频率 fm=20e3;%基带频率 n=2*(6*fsfm); final=(1fs)*(n-1); fc=2e5; % 载波频

44、率 t=0:1fs:(final); Fn=fs2;%耐奎斯特频率 %用正弦波产生方波 twopi_fc_t=2*pi*fm*t; A=1; phi=0; x = A * cos(twopi_fc_t + phi); % 方波 am=1; x(x0)=am; x(x= vt Vs(ii)=High; else Vs(ii)=Low; end end Vo=Vs; subplot(325); plot (t,Vo), title(解调后输出信号), axis(0 2e-4 -5 5) grid on; xlabel(时间 (s), ylabel(幅度(V)仿真结果： 移键控(4PSK4PSK信号

45、的产生与解调在进制数字相位调制中，四进制绝对移相键控4PSK，又称QPSK和四进制差分相位键控4DPSK，又称QDPSK用的最为广泛。下面着重介绍多进制数字相位调制的这两种形式。4PSK利用载波的四种不同相位来表征数字信息。由于每一种载波相位代表两个比特信息，故每个四进制码元又被称为双比特码元，习惯上把双比特的前一位用代表，后一位用代表。14PSK信号的产生多相制信号常用的产生方法有：直接调相法及相位选择法。1相位选择法由式8-54可以看出，在一个码元持续时间内，4PSK信号为载波四个相位中的某一个。因此，可以用相位选择法产生4PSK信号，其原理如图8-27所示。图中，四相载波发生器产生4PS

46、K信号所需的四种不同相位的载波。输入的二进制数码经串并变换器输出双比特码元。按照输入的双比特码元的不同，逻辑选相电路输出相应相位的载波。例如，B方式情况下，双比特码元为11时，输出相位为45?的载波；双比特码元为01时，输出相位为135?的载波等。图8-27相位选择法产生4PSK信号B方式方框图图中产生的是B方式的4PSK信号。要想形成A方式的4PSK信号，只需调整四相载波发生器输出的载波相位即可。2直接调相法由式8-56可以看出，4PSK信号也可以采用正交调制的方式产生。B方式4PSK时的原理方框图如图8-28a所示。它可以看成是由两个载波正交的2PSK调制器构成，分别形成图8-28b中的虚

47、线矢量，再经加法器合成后，得图b中实线矢量图。显然其为B方式4PSK相位配置情况。图8-28直接调相法产生4PSK信号方框图假设要产生4PSK的A方式波形，只需适当改变振荡载波相位就可实现。24PSK信号的解调由于4PSK信号可以看作是两个载波正交的2PSK信号的合成，因此，对4PSK信号的解调可以采用与2PSK信号类似的解调方法进行。图8-29是B方式4PSK信号相干解调器的组成方框图。图中两个相互正交的相干载波分别检测出两个分量和，然后，经并串变换器复原成二进制双比特串行数字信号，从而实现二进制信息恢复。此法也称为极性比较法。图8-294PSK信号的相干解调假设解调4PSK信号A方式，只需

48、适当改变相移网络。在2PSK信号相干解调过程中会产生“倒即“180相位模糊现象。同样，对于4PSK信号相干解调也会产生相位模糊问题，并且是0、90、180和270四个相位模糊。4PSK基于MATLAB的仿真% 调相法clear allclose allt=-1:0.01:7-0.01;tt=length(t);x1=ones(1,800);for i=1:tt if (t(i)=-1 & t(i)=5& t(i)=0 & t1(i)=4& t1(i)0 data_recover_a(i:i+19)=1; bit_recover=bit_recover 1; else data_recover_

49、a(i:i+19)=-1; bit_recover=bit_recover -1; endenderror=0;dd = -2*bit_in+1;ddd=dd;ddd1=repmat(ddd,20,1);for i=1:2e4 ddd2(i)=ddd1(i);endfor i=1:1e3 if bit_recover(i)=ddd(i) error=error+1; endendp=error1000;figure(1)subplot(2,1,1);plot(t2,ddd2);axis(0 100 -2 2);title(原序列);subplot(2,1,2);plot(t2,data_rec

50、over_a);axis(0 100 -2 2);title(解调后序列); 仿真结果2：% 设定 T=1, 不加噪声clear allclose all% 调制bit_in = randint(1e3, 1, 0 1);bit_I = bit_in(1:2:1e3);bit_Q = bit_in(2:2:1e3);data_I = -2*bit_I+1;data_Q = -2*bit_Q+1;data_I1=repmat(data_I,20,1);data_Q1=repmat(data_Q,20,1);for i=1:1e4 data_I2(i)=data_I1(i); data_Q2(i)

51、=data_Q1(i);end;t=0:0.1:1e3-0.1;f=0:0.1:1;xrc=0.5+0.5*cos(pi*f);data_I2_rc=conv(data_I2,xrc)5.5;data_Q2_rc=conv(data_Q2,xrc)5.5;f1=1;t1=0:0.1:1e3+0.9;I_rc=data_I2_rc.*cos(2*pi*f1*t1);Q_rc=data_Q2_rc.*sin(2*pi*f1*t1);QPSK_rc=(sqrt(12).*I_rc+sqrt(12).*Q_rc);% 解调I_demo=QPSK_rc.*cos(2*pi*f1*t1);Q_demo=Q

52、PSK_rc.*sin(2*pi*f1*t1);I_recover=conv(I_demo,xrc);Q_recover=conv(Q_demo,xrc);I=I_recover(11:10010);Q=Q_recover(11:10010);t2=0:0.05:1e3-0.05;t3=0:0.1:1e3-0.1;data_recover=;for i=1:20:10000 data_recover=data_recover I(i:1:i+19) Q(i:1:i+19);end;ddd = -2*bit_in+1;ddd1=repmat(ddd,10,1);for i=1:1e4 ddd2(

53、i)=ddd1(i);endfigure(1)subplot(4,1,1);plot(t3,I);axis(0 20 -6 6);subplot(4,1,2);plot(t3,Q);axis(0 20 -6 6);subplot(4,1,3);plot(t2,data_recover);axis(0 20 -6 6);subplot(4,1,4);plot(t,ddd2);axis(0 20 -6 6); % QPSK误码率分析SNRindB1=0:2:10;SNRindB2=0:0.1:10;for i=1:length(SNRindB1) pb,ps=cm_sm32(SNRindB1(i)

54、; smld_bit_err_prb(i)=pb; smld_symbol_err_prb(i)=ps;end;for i=1:length(SNRindB2) SNR=exp(SNRindB2(i)*log(10)10); theo_err_prb(i)=Qfunct(sqrt(2*SNR);end;title(QPSK误码率分析);semilogy(SNRindB1,smld_bit_err_prb,*);axis(0 10 10e-8 1);% semilogy(SNRindB1,smld_symbol_err_prb,o);仿真结果3：% semilogy(SNRindB1,smld_

55、symbol_err_prb,o);semilogy(SNRindB2,theo_err_prb);legend(仿真比特误码率,理论比特误码率);y=Qfunct(x)y=(12)*erfc(xsqrt(2);functionpb,ps=cm_sm32(SNRindB)N=10000;E=1;SNR=10(SNRindB10);sgma=sqrt(ESNR)2;s00=1 0;s01=0 1;s11=-1 0;s10=0 -1;for i=1:N temp=rand; if (temp0.25) dsource1(i)=0; dsource2(i)=0; elseif (temp0.5) d

56、source1(i)=0; dsource2(i)=1; elseif (temp0.75) dsource1(i)=1; dsource2(i)=0; else dsource1(i)=1; dsource2(i)=1; end;end;numofsymbolerror=0;numofbiterror=0;for i=1:N n=sgma*randn(size(s00); if(dsource1(i)=0)&(dsource2(i)=0) r=s00+n; elseif(dsource1(i)=0)&(dsource2(i)=1) r=s01+n; elseif(dsource1(i)=1)

57、&(dsource2(i)=0) r=s10+n; else r=s11+n; end; c00=dot(r,s00); c01=dot(r,s01); c10=dot(r,s10); c11=dot(r,s11); c_max=max(c00 c01 c10 c11); if (c00=c_max) decis1=0;decis2=0; elseif(c01=c_max) decis1=0;decis2=1; elseif(c10=c_max) decis1=1;decis2=0; else decis1=1;decis2=1; end; symbolerror=0; if(decis1=d

58、source1(i) numofbiterror=numofbiterror+1; symbolerror=1; end; if(decis2=dsource2(i) numofbiterror=numofbiterror+1; symbolerror=1; end; if(symbolerror=1) numofsymbolerror=numofsymbolerror+1; end;end;ps=numofsymbolerrorN;pb=numofbiterror(2*N);仿真结果4：正交幅度调制(16QAM正交调制及相干解调原理框图正交调制原理框图相干解调原理框图16QAM调制与解调基于

59、MATLAB的仿真：仿真程序：clc;close all;N=60; %设定码元数量 fb=1; %基带信号频率 fs=32; %抽样频率 fc=4; %载波频率,为便于观察已调信号，我们把载波频率设的较低 Kbase=2; % Kbase=2,基带经成形滤波器滤波后，再进行调制 info=random_binary(N); %产生二进制信号序列 y,I,Q=qam(info,Kbase,fs,fb,fc); %对基带信号进行16QAM调制 y1=y; y2=y; %备份信号，供后续仿真用 T=length(info)fb; m=fsfb;nn=length(info); dt=1fs; t=

60、0:dt:T-dt; subplot(221);plot(info);title(原序列)subplot(222); %显示的已调信号 %由于测试信号码元数量为10000个，在这里我们只显示其总数的110 plot(t(1:1000),y(1:1000),t(1:1000),I(1:1000),t(1:1000),Q(1:1000),0 20,0 0); title(已调信号输出波形); %傅里叶变换，求出已调信号的频谱 n=length(y); y=fft(y)n; y=abs(y(1:fix(n2)*2; q=find(y1e-04); y(q)=1e-04; y=20*log10(y);