基于MATLAB的PSK调制和解调及仿真

基于MATLAB的PSK调制和

解调及仿真

—I—

摘要 摘要

Psk调制是通信系统中最为重要的环节之一，Psk调制技术的改进也是通信系统性能提高的重要途径。本文首先分析了数字调制系统的基本调制解调方法，然后，运用Matlab及附带的图形仿真工具 ——Simulink设计了这几种数字调制方法的仿真模型。通过仿真，观察了调制解调过程中各环节时域和频域的波形，并结合这几种调制方法的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。最后，在仿真的基础上分析比较了各种调制方法的性能，并通过比较仿真模型与理论计算的性能，证明了仿真模型的可行性。另外，本文还利用Matlab的图形用户界面（GUI）功能为仿真系统设计了一个便于操作的人机交互界面，使仿真系统更加完整，操作更加方便。

关键词：数字调制；分析与仿真；Matlab；Simulink；GUI图形界面

—II—

ABSTRACT ABSTRACT

In this paper, methods of psk modulation are introduced

firstly. Then their simulation models are built by using MATLAB’s simulation tool, SIMULINK. Through observing the results of simulation, the factors that affect the capability of the psk modulation system and the reliability of the simulation models are analyzed. And then, the capability of three digital modulation simulation models, 2-PSK, 4-PSK and , have been compared, as well as comparing the results of simulation and theory. At last, the conclusion is gotten: The simulation models are reasonable. In addition, an operation interface is designed, which can simplify the manipulation of the simulation system, by mean of the Graphical User Interface, which short for GUI.

Keywords: PSK modulation; analysis; simulation; MATLAB; SIMULINK; GUI

—III—

目录 目录

摘要 ..................................................... II ABSTRACT ................................................ III 目录 ..................................................... IV 前言 ...................................................... 1 1绪论 .................................................... 2

1.1通信技术的历史和发展 ............................... 2 1.2数字调制技术 ....................................... 3 1.3数字调制的发展现状和趋势 ........................... 4 1.4本章小结 ........................................... 5 2 MATLAB仿真技术 ......................................... 6

2.1通信仿真 ........................................... 6 2.2 MATLAB简介 ........................................ 9 2.3 Simulink简介 ..................................... 12 2.4 本章小结 .......................................... 14 3 PSK 调制系统 ........................................... 15

3.1 2PSK数字调制原理 ................................. 15 3.2 4PSK的调制和解调 ................................. 19 3.3 本章小结 .......................................... 23 4 PSK调制解调系统的仿真 ................................. 24

4.1 2PSK调制解调系统的仿真 ........................... 24 4.2 4PSK调制解调系统的仿真 ........................... 25 4.3利用MATLAB研究4PSK信号 .......................... 27 4.4 本章小结 .......................................... 29 结论 ..................................................... 30 参考文献 ................................................. 31 致谢 ..................................................... 32 附录 ..................................................... 33

—IV—

黑龙江八一农垦大学毕业设计（论文） 前言

现代通信的发展趋势为数字化，随着现代通信技术的不断开发，数字调制技术已日趋成熟，在各个领域都得到了广泛的应用和认同。因此本文对PSK的数字调制和解调进行了仿真。现代社会发展要求通信系统功能越来越强，性能越来越高，构成越来越复杂；这就要借助于功能强大的计算机辅助分析设计技术和工具才能实现。现代计算机科学技术快速发展，已经研发出了新一代的可视化的仿真软件。这些功能强大的仿真软件，使得通信系统仿真的设计和分析过程变得相对直观和便捷，由此也使得通信系统仿真技术得到了更快的发展。本文使用的是功能强大的MATLAB软件。

MATLAB是一种使用简便的、特别适用于科学研究和工程计算的高级语言，与其他计算机语言相比，它的特点是简洁和智能化，具有极高的编程和调试效率。通过使用MATLAB工具箱函数对数字调制进行仿真，更能直观彻底的掌握数字通信，数字调制的原理。有助于我们的学习和研究，加深对知识的理解和运用。 MATLAB的便利性还体现在它的仿真结果还可以存放到MATLAB 的工作空间里做事后处理。方便我们修改参数对不同情况下的输出结果进行对比。

由于MATLAB和SIMULINK是集成在一起的，因此用户可以在这两种环境下对自己的模型进行仿真、分析和修改。

—1—

黑龙江八一农垦大学毕业设计（论文） 1绪论

1.1通信技术的历史和发展

1.1.1通信的概念

通信就是克服距离上的障碍，从一地向另一地传递和交换消息。消息是信息源所产生的，是信息的物理表现，例如，语音、文字、数据、图形和图像等都是消息（Message）。消息由模拟消息（如语音、图像等）以及数字消息（如数据、文字等）之分。所有消息必须在转换成电信号（通常简称为信号）后才能在通信系统中传输。所以，信号（Signal）是传输消息的手段，信号是消息的物资载体。

相应的信号可以分为模拟信号和数字信号，模拟信号的自变量可以是连续的或离散的，但幅度是连续的，如电话机、电视摄像机输出的信号就是模拟信号。数字信号的自变量可以是连续的或离散的，但幅度是离散的，如计算机等各种数字终端设备输出的信号就是数字信号。

通信的目的是传递消息，但对受信者有用的是消息中包含的有效内容，即信息（Information）。消息是具体的、表面的，而信息是抽象的、本质的，且消息中包含的信息的多少可以用信息量来度量。通信技术，特别是数字通信技术近年来发展非常迅速，它的应用越来越广泛。

数字通信系统较模拟通信系统而言，具有抗干扰能力强、便于加密、易于实现集成化、便于与计算机连接等优点。因而，数字通信更能适应对通信技术的高要求。

1.1.2通信的发展史简介

远古时代，远距离的传递消息是以书信的形式来完成的，这种通信方式明显具有传递时间长的缺点。为了在尽量短的时间内传递尽量多的消息，人们不断地尝试所能找到的各种最新技术手段。1837年发明的莫尔斯电磁式电报标志着电通信的开始。之后，利用电进行通信的研究取得了长足的进步。1866年利用海底电缆实现了跨大西洋的越洋电报通信。1876年贝耳发明了电话，利用电信号实现了语音信号的有线传递，使信息的传递变得既迅速又准确，这标志着模拟通信的开始，由于它比电报更便于交流使用，所以直到20世纪前半叶这种采用模拟技术的电话通信技术比电报得到了更为迅速和广泛的发展。1937年瑞威斯发明的脉冲编码调制标志数字通信的开始。20世纪60年代以后集成电路、电子计算机的出现，使得数字通信迅速发展。在70年代末在全球发展起来的模拟移动电话在90年代中期被数字移动电话所代替，现有的模拟

—2—

黑龙江八一农垦大学毕业设计（论文） 电视也正在被数字电视所代替。数字通信的高速率和大容量等各方面的优越性也使人们看到了它的发展前途

1.2数字调制技术

通信按照传统的理解就是信息的传输与交换。在当今信息社会，通信则与遥感，计算技术紧密结合，成为整个社会的高级“神经中枢”。没有通信，人类社会是不可想象的。一般来说，社会生产力水平要求社会通信水平与之相适应。若通信水平跟不上，社会成员之间的合作程度就受到限制。可见，通信是十分重要的。

通信传输的消息是多种多样的，可以是符号的，文字的，数据和图像的等等。各种不同的消息可以分为两类：一类称为离散消息；另一类称为连续消息。离散消息的状态是可数的或离散的，比如符号，文字或数据等。离散消息也称数字消息。而连续消息则是其状态连续变化的消息，例如，连续变化的语音，图像等。连续消息也称模拟消息。因此按照信道中传输的是模拟信号还是数字信号可以将通信系统分为模拟通信系统和数字通信系统。

数字通信有以下突出的特点：第一，数字信号传输时，信道噪声或干扰所造成的差错，原则上是可以控制的。第二，当需要保密的时候，可以有效的对基带信号进行人为的“扰乱”，即加上密码。

数字通信系统可以用下图表示：

信息源 加密器 编码器 调制器 信 道 解调器 译码器 解密器 受信者 噪 声 源 数字通信在近20年来得到了迅速的发展，其原因是： （1） 抗干扰能力强

（2） 便于进行各种数字信号处理

—3—

黑龙江八一农垦大学毕业设计（论文） （3） 易于实现集成化

（4） 经济效益正赶上或超过模拟通信

（5） 传输与交换可结合起来，传输电话与传输数据也可结合起来，成为一个

统一整体，有利于实现综合业务通信网。

本设计主要研究数字通信过程中的调制解调过程。从原理上说受调载波可以是任意的，只要已调信号适合心动的传输就可以了，但是实际上，大多数通信系统中，都选择正弦信号作为载波。这是因为正弦信号简单，便于产生和接收。

1.3数字调制的发展现状和趋势

进入20世纪以来，随着晶体管、集成电路的出现与普及、无线通信迅速发展。特别是在20世纪后半叶，随着人造地球卫星的发射，大规模集成电路、电子计算机和光导纤维等现代技术成果的问世，通信技术在以下几个不同方向都取得了巨大的成功。

（1） 微波中继通信使长距离、大容量的通信成为了现实。

（2） 移动通信和卫星通信的出现，使人们随时随地可通信的愿望可以实现。 （3） 光导纤维的出现更是将通信容量提高到了以前无法想象的地步。 （4） 电子计算机的出现将通信技术推上了更高的层次，借助现代电信网和计

算机的融合，人们将世界变成了地球村。

（5） 微电子技术的发展，使通信终端的体积越来越小，成本越来越低，范围

越来越广。例如2003年我国的移动电话用户首次超过了固定电话用户。根据国家信息产业部的统计数据，到2005年底移动电话用户近4亿。 随着现代电子技术的发展，通信技术正向着数字化、网络化、智能化和宽带化的方向发展。随着科学技术的进步，人们对通信的要求越来越高，各种技术会不断地应用于通信领域，各种新的通信业务将不断地被开发出来。到那时人们的生活将越来越离不开通信。

本文中提到的调制方式大都是可以实用的，已经采用多年，并且至今仍然被采用着。但是，这些调制方法还不是很完善，有许多值得改进之处。因此，在这些基本的数字调制方法基础上，多年来不断研究出新的或改进的调制方法。实际上，在基本的和先进的调制方法之间并没有明确的界限。这些方法都是不间断地发展出来的，后来者自然比原有者更先进。

此外，随着技术的进步，特别是超大规模集成电路和数字信号处理技术的发展，使得复杂的电路设计得以用少量的几块即成电路模块实现，有些硬件电路的功能还可以用软件代替实现。因此使得一些较复杂的调制技术能够容易地实现并投入使用。这方面的条件使得新的更复杂的调制体制迅速地不断涌现。

—4—

黑龙江八一农垦大学毕业设计（论文） 目前，改进的数字调制方式主要有偏置正交相移键控， /4正交差分相移键控，最小频移键控，高斯最小频移键控，正交频分复用，网格编码调制等，这里对最小频移键控作一介绍

1.4本章小结

本章主要讲述了通信、调制、数字通信的概念，发展以及在生活中的应用价值。从大的方面到小的方面，层层递进一直到本文所要讲的主旨内容数字调制。

—5—

黑龙江八一农垦大学毕业设计（论文） 2 MATLAB仿真技术

2.1通信仿真

2.1.1通信仿真简介

随着科学技术、仿真理论及计算机的不断发展，仿真技术不断提高。在如今的科学研究中，仿真技术提高了科学研究水平，缩短了科学研究周期、降低了科学研究成本及风险、促进了各不同领域学科融合、加速了科研成果转化为生产力。可以说仿真技术已成为科学研究中必不可少的实用技术。因此在现代科学研究及应用中，仿真技术被广泛应用于数学、物理、电子、通信、医学、生物等众多领域。

所谓系统仿真，通俗的说就是模型试验，它是指通过系统模型试验去研究一个已经存在的或者正在设计的系统的过程。系统仿真，它不是对原型的简单再现，而是按照研究的侧重点对系统进行提炼，以利于研究者抓住问题的本质，这种建立在模型系统上的实验技术，称为仿真技术。仿真方法可以分为3类：实物仿真、数学仿真和半实物仿真。系统仿真的研究重点在于仿真环节，即在模型建立之后，设计适当的算法，并编制成计算机程序。因此，便产生了很多仿真算法和仿真软件，其中以MATLAB提供的动态仿真工具Simulink最为耀眼，它不仅具有强大的功能，并且具有很好的使用性。MATLAB就是大量的计算机仿真软件中的优秀代表，它在科学研究特别是电子信息科学中有着极为广泛的应用。MATLAB现已被广泛应用于数学、通信、信号处理、自动控制、神经网络、图形处理等许多不同学科的研究中。

仿真是衡量系统性能的工具，它通过仿真模型的仿真结果来推断原系统的性能，从而为新系统的建立或原系统的改造提供可靠的参考。实际的通信系统是一个功能结构相当复杂的系统，对这个系统做出的任何改变都可能影响到整个系统的性能和稳定。因此，在对原有的通信系统做出改进或建立一个新系统之前，通常需要对这个系统进行建模和仿真，通过仿真结果衡量方案的可行性，从中选择最合理的系统配置和参数设置，然后再应用于实际系统中。这个过程就是通信仿真。通信系统仿真可以分成离散事件仿真和连续事件仿真。在离散事件仿真中，仿真系统只对离散事件做出反应，而在连续事件仿真中，仿真系统对输入信号产生连续的输出信号。离散事件仿真是对实际通信系统的一种简化，它的仿真建模比较简单，整个仿真过程需要花费的时间也比连续仿真少。虽然离散事件仿真舍弃了一些仿真细节，在有些场合显得不够具体，但仍然是通信系统仿真的主要形式。

与一般的仿真过程类似，在对通信系统实施仿真之前，首先需要研究通信系统的特性，通过归纳和抽象建立通信系统的仿真模型。

—6—

黑龙江八一农垦大学毕业设计（论文） 2.1.2 通信仿真的一般步骤

通信系统仿真一般分成3个步骤，即仿真建模、仿真实验和仿真分析。应该注意的是，通信仿真是一个螺旋式发展的过程，因此这3个步骤可能需要循环执行多次之后才能够获得令人满意的仿真结果。 （1）仿真建模

仿真建模是根据实际通信系统建立仿真模型的过程，它是整个通信系统仿

真过程中的一个关键步骤，因为仿真模型的好坏直接影响着仿真结果的真实性和可靠性。

仿真模型一般是一个数学模型。数学模型有多种分类方式，包括确定性模型和随机性模型，静态模型和动态模型。确定性模型的输入变量和输出变量都有固定数值，而在随机模型中，至少有一个输入变量是随机的。静态模型不需要考虑时间变化因素，动态模型的输入输出变量则需要考虑时间变化因素。一般情况下，通信系统模型是一个随机动态系统。 在仿真建模过程中，首先需要分析实际系统存在的问题或设立系统改造的目标，并且把这些问题和目标转化成数学变量和公式。有了这些具体的仿真目标之后，下一步是获取实际通信系统的各种参数，如通信系统占用的带宽及其频率分布，系统对于特定的输入信号产生的输出等。同时，对于通信系统中的各个随机变量，可以采集这些变量的数据，然后通过数学工具来确定随机变量的分布特性。

有了上面的准备工作，下一步就可以通过仿真软件来建立仿真模型了。 （2）仿真实验

仿真实验是一个或一系列针对仿真模型的测试。在仿真实验中，通常需要

多次改变仿真模型输入信号的数据，以观察和分析仿真模型对这些输入信号的反应，以及仿真系统在这个过程中表现出来的性能。需要强调的一点是，仿真过程中使用的输入数据必须具有一定的代表性，即能够从各个角度显著地改变仿真输出信号的数值。

实施仿真之前需要确定的另外一个因素是性能尺度。性能尺度指的是能够衡量仿真过程中系统性能的输出信号的数值（或根据输出信号计算得到的数值），因此，在实施仿真之前，首先需要确定仿真过程中应该收集哪些仿真数据，这些数据以什么样的格式存在，以及收集多少数据。在明确了仿真系统对输入信号的要求之前，最好把这些设置整理成一份简单的文档。编写文档是一个好习惯，它能够帮助我们回忆起仿真过程的一些细节。当然，文档的编写不一定要求很规范，并且文档的大小应视仿真设计的规模而定。最后，还应该明确各个输入信号的初始设置以及仿真系统内部各个状态的初始值。

仿真的运行实际上是计算机的计算过程，这个过程一般不需要人工干预，

—7—

黑龙江八一农垦大学毕业设计（论文） 花费的时间由仿真的复杂度确定。如果需要比较仿真系统在不同参数设置下的性能，应该使仿真系统在取不同参数值时具有相同的输入信号，这样才能保证分析和比较的客观性和可靠性。 （3）仿真分析

仿真分析是一个通信系统仿真流程中的最后一个步骤。在仿真分析过程

中，用户已经从仿真过程中获得了足够多的关于系统性能的信息，但是这些信息只是一个原始的数据，一般还需要经过数值分析和处理才能够获得衡量系统性能的尺度，从而获得对仿真系统的一个总体评价。常用的系统性能尺度包括平均值，方差，标准差，最大值和最小值等，它们从不同的角度描绘了仿真系统的性能。

图表是最简洁的说明工具，它具有很强的直观性，便于分析和比较，因此，仿真分析的结果一般都绘制成图表形式。我们使用的仿真工具一般都具有很强的绘图功能，能够便捷地绘制各种类型的图表。

2.1.3 仿真技术在通信系统设计中的作用

仿真技术的主要应用之一就是辅助系统的设计过程，也就是计划、设计和实现通信系统的过程。实际系统的设计不同于传统的分析问题，它有自己的特点。在设计过程中，不能依靠对系统各单元特性的详细而精密的了解，这里面有现实而又重要的原因。

系统设计师的基本目的是制造出一个系统使它在预计的使用寿命内满足一定的性能指标。假如我们已经有了一个能满足一定性能的具体设计方案，所谓具体设计方案，就是说系统中的每个部件都唯一地用一定的函数来描述.

2.1.4仿真的应用

我们主要通过以下几种仿真方法，讲述目前正迅猛发展的仿真技术的许多的新的思想。

（1） 针对光纤系统，研究的重点是光检测器的特性。光纤系统中误码率是很小的，这样我们就可能在该例子采用运行时间减小技术配合蒙特卡罗仿真方法。

（2） 主要研究多径衰落环境下64QAM无线中继系统的错误概率性能（BER）的估计。在该系统中，不同的QA输出可表示为单个基本脉冲时延加权的和，仿真就是产生这些基本脉冲。进而计算BER。

（3） 主要是检测数字卫星通信系统的环境。这个系统中含有非线性，这是一个难以处理的问题。另一个问题是模型描述越详细，仿真就变得复杂费时。在这个研究中，我们应当掌握复杂性与准确性的平衡，在降低复杂性时，采用一些技巧避免准确性的过分降低。

（4） 野战网初始拓扑的优化问题。由于网络节点和用户都在运动之中，节

—8—

黑龙江八一农垦大学毕业设计（论文） 点数和用户数都可能十分大，要保证野战网覆盖制定的地域范围，并在很恶劣的战场环境下具有很好的抗毁能力，网络的设计必须考虑很多参数，这可能使网络仿真的运行时间过长，启发算法利用迭代逼近，使得设计结果与仿真运行时间能取得较好的折衷。

（5） 关于自适应滤波器的性能仿真.

2.2 MATLAB简介

2.2.1 MATLAB 的介绍及发展

MATLAB软件系列产品是一套高效强大的工程技术数值运算和系统仿真软件，广泛应用于当今的航空航天、汽车制造、半导体制造、电子通信、医学研究、财经研究和高等教育等领域，被誉为“巨人肩膀上的工具”。研发人员借助MATLAB软件能迅速测试设想构想，综合评测系统性能，快速设计更好方案来确保更高技术要求。同时MATLAB也是国家教委重点提倡的一种计算工具。

MATLAB 软件是美国 Math works 公司的产品，MATLAB 是英文 MATrix LABoratory(矩阵实验室)的缩写。从1984年推出了它的第一个DOS版本至今，一经推出了6.5版。Matrix Laboratory意为“矩阵实验室”，从它的本意可以知道，最初的MATLAB只是一个数学计算工具。但现在的MATLAB已经远不仅仅是一个“矩阵实验室”，它已经成为一个集概念设计、算法开发、建模仿真、实时实现于一体的集成环境，它拥有许多衍生的子集工具[9]。

新的版本集成了日常数学处理中的各种功能，包括高效的数值计算、矩阵运算、信号处理和图形生成等等的常用功能。在 MATLAB 环境下，用户可以集成地进行程序设计、数值计算、图形绘制、输入输出、文件管理等各项操作。 MATLAB 提供了一个人机交互的数学系统环境，该系统的基本数据结构是矩阵，在生成矩阵对象时，不要求作明确的维数说明，所谓交互式语言，是指人们给出一条命令，立即就可以得出该命令的结果。该语言无需像 C 和 Fortran 语言那样，首先要求使用者去编写源程序，然后对之进行编译、连接，最终形成可执行文件。这无疑会给使用者带来了极大的方便，因此，利用 MATLAB可以节省大量的编程时间。2002年6月Mathworks公司正式推出MATLAB Release 13，即MATLAB 6.5／Simulink 5.0 这是目前应用最广的版本。

MATLAB主要由C语言编写而成，采用LAPACK 为底层支持软件包。 MATLAB的编程非常简单，它有着比其他任何计算机高级语言更高的编程效率、更好的代码可读性和移植性，以致被誉为“第四代”计算机语言，MATLAB是所有MathWorks公司产品的数值分析和图形基础环境。此外MATLAB 还拥有强大的2D和3D甚至动态图形的绘制功能，这样用户可以更直观、更迅速的进行多种算法的比较，从中找出最好的方案。

—9—

黑龙江八一农垦大学毕业设计（论文） 从通信系统分析与设计、滤波器设计、信号处理、小波分析、神经网络到控制系统、模糊控制等方面来看，MATLAB提供了大量的面向专业领域的工具箱。通过工具箱，以往需要复杂编程的算法开发任务往往只需一个函数就能实现，而且工具箱是开放的可扩展集，用户可以查看或修改其中的算法，甚至开发自己的算法。

目前， MATLAB已经广泛地应用于工程设计的各个领域，如电子、通信等领域；它已成为国际上最流行的计算机仿真软件设计工具。现在的MATLAB不再仅仅是一个矩阵实验室，而是一种实用的、功能强大的、不断更新的高级计算机编程语言。

2.2.2 MATLAB的主要特点

（1） 科学计算

MATLAB是一种面向科学与工程计算的高级语言，允许采用数学形式的语言编写程序。用 MATLAB编写程序犹如在演算纸上排列出公式与求解问题，因此MATLAB语言被称为“演算纸式”科学计算语言。MATLAB拥有五百多种数学、统计及工程函数，可使用户立刻实现所需强大的数学计算功能。 （2） 先进的可视化工具

MATLAB提供功能强大的、交互式的二维绘图和三维绘图功能。可使用户创建富有表现力的彩色图形。 （3） 直观灵活的语言

MATLAB不仅仅是是一套打好包的函数库，同时也是一种高级的、面向对象的编程语言。

（4） 开放性、可扩展性强

M文件是可见的MATLAB程序，所以用户可以查看源代码。开放的系统设计使用户能够检查算法的正确性，修改已存在的函数，或者加入自己的新部件。 （5） 特殊应用工具箱

MATLAB的工具箱加强了对工程及科学种特殊应用的支持。工具箱也和MATLAB一样是完全用户化的，可扩展性强。将某个或某几个工具箱与MATLAB联合使用，可以得到一个功能强大的计算组合包，满足用户的特殊要求。 （6） 使用方便

MATLAB语言是一种解释型语言，执行之前不需要进行专门的编译。一般情况下，在采用任何高级语言编写和调试程序时需要经历4个阶段，即编辑、编译、链接以及执行调试，并且这四个步骤之间是顺次执行的。MATLAB语言与其他语言相比，较好地解决了上述的问题。 （7） 强大的图形处理功能[10]

2.2.3 MATLAB的功能

MATLAB 之所以成为世界顶级的科学计算与数学应用软件,是因为它随着版本的升级与不断完善而具有愈来愈强大的功能。

—10—

黑龙江八一农垦大学毕业设计（论文） （1） 数值计算功能。MATLAB 出色的数值计算功能是使之优于其他数学应用软件的决定性因素之一。

（2） 符号计算功能。MATLAB 符号运算的独特之处:无须事先对变量赋值,而所得的结果以标准的符号形式表达,符号计算的整个过程以字符进行。

（3） 数据分析功能。MATLAB 可以给计算数据以二维、三维乃至四维的图形表现。这不仅使数据间的关系清晰明了,而且对于揭示其内在本质有着非常重要的作用。

（4） 动态仿真功能。MATLAB 提供了一个模拟动态系统的交互程序SIMULINK用户通过简单的鼠标操作,就可建立起直观的系统模型,并进行仿真[14]

。

（5） 图形文字统一处理功能。MATLAB Notebook成功地将Microsoft Word 与MATLAB 集成为一个整体,为文字处理、科学计算、工程设计营造了一个完美统一的工作环境。它既拥有Word 强大的文字处理功能,又能从Word 访问MATLAB 的数据计算和可视化结果。

2.2.4 MATLAB系统的组成

（1） 编程语言：它是以矩阵和数组为基本单位的编程语言。具有条件控制、函数调用、数据结构、输入输出、面向对象等程序语言特性。利用它既可以进行小规模编程，完成算法设计和算法实验的基本任务，也可以进行大规模编程，开发复杂的应用程序。 （2） 工作环境： 包括了一系列的应用工具，提供编写和调试程序的环境。 （3） 图形处理：这是MATLAB 图形系统的基础，包括绘制二维和三维图形，也包括用户对图形图像等对象进行特性控制的低层 MATLAB 命令，以及开发 GUI 应用程序的各种工具。

（4） 数学函数库：包括了大量的数学函数，也包括复杂的功能。这是对 MATLAB 使用的各种数学算法的总称。

（5） 应用程序接口(API)：提供接口程序， 使MATLAB可以与用其他语言

编写的程序进行交互。

2.2.5 MATLAB的仿真应用及应用领域

现在从电子通信、自动控制图形分析处理到航天工业、汽车工业，甚至是财务工程。MATLAB都凭借其强大的功能获得了极大的用武之地。广大学生可以使用MATLAB来帮助进行信号处理、通信原理、线性系统、自动控制等课程的学习；科研工作者可以使用MATLAB进行理论研究和算法开发；工程师可以使用MATLAB进行系统级的设计与仿真。如今MATLAB在系统级设计和仿真方面的强大能力已经获得了包括Texas Instruments和Motorola这样的世界著名公司的认可，Texas Instruments公司的工程师们使用MATLAB直接在可复用高层系统模型下面提炼实现细节，同时生成实时软件原型；Motorola公司的IC设计师甚至认为Simulink的使用大大加快了产品设计的仿

—11—

黑龙江八一农垦大学毕业设计（论文） 真速度，使仿真时间由原来的几小时、几天缩短到了几分钟！ MATLAB的主要应用领域包括： （1） 仿真和建模 （2） 实时仿真 （3） 自动控制

（4） 信号处理与通信 （5） 数据分析/科学计算 （6） 算法开发

（7） 图形和可视化法 （8） 独立应用开发

2.3 Simulink简介

2.3.1 Simulink相关内容

近几年，在学术界和工业领域，Simulink已成为在动态系统领域建模和仿真方面分，Simulink具有相对独立的功能和使用方法。确切的说，它是一个用来对动态系统进行建模、应用最广泛的软件包之一 。它的魅力在于强大的功能和简便的操作。作为MATLAB的重要组成部仿真和分析的软件包[5]。它支持连续、离散及两者混合的线性和非线性系统，也支持具有多种采样频率的系统，而且系统可以是多进程的。在Simulink环境中，利用鼠标就可以在模型窗口中直观地\"画\"出系统模型，然后直接进行仿真。它为用户提供了方框图进行建模的图形接口，采用这种结构画模型就像你用手和纸来画一样容易。它与传统的仿真软件包微分方程和差分方程建模相比，具有更直观、方便、灵活的优点。

Simulink包含有SINKS（输出方式）、SOURCE（输入源）、LINEAR（线性环节）、NONLINEAR（非线性环节）、CONNECTIONS（连接与接口）和EXTRA（其他环节）子模型库，而且每个子模型库中包含有相应的功能模，用户也可以定制和创建用户自己的模块[5]。用Simulink创建的模型可以具有递阶结构，因此用户可以采用从上到下或从下到上的结构创建模型。用户可以从最高级开始观看模型，然后用鼠标双击其中的子系统模块，来查看其下一级的内容，以此类推，从而可以看到整个模型的细节，帮助用户理解模型的结构和各模块之间的相互关系。在定义完一个模型后，用户可以通过Simulink的菜单或MATLAB的命令窗口键入命令来对它进行仿真。菜单方式对于交互工作非常方便，而命令行方式对于运行一大类仿真非常有用。采用SCOPE模块和其他的画图模块，在仿真进行的同时，就可观看到仿真结果。除此之外，用户还可以在改变参数后来迅速观看系统中发生的变化情况。仿真的结果还可以存放到MATLAB的工作空间里做事后处理。

—12—

黑龙江八一农垦大学毕业设计（论文） 模型分析工具包括线性化和平衡点分析工具、MATLAB的许多工具及MATLAB的应用工具箱。由于MATLAB和SIMULINK的集成在一起的，因此用户可以在这两种环境下对自己的模型进行仿真、分析和修改。

2.3.2 Smulink仿真原理

Simulink通过利用模块组合的方法可以方便用户快速、准确地创建动态系统的计算机模型。它可以用来模拟线性与非线性系统，连续与非连续系统，或者这些混合的系统，是强大的系统仿真工具[27]。 1. 典型的模型结构

Simulink的每一个模块实际上都是一个系统、一个典型的Simulink模块包括输入、状态和输出三个部分：

（1） 输入模块：即信号源模块，包括常数信号源、函数信号发生器和用户自定义信号；

（2） 状态模块：即被模拟的系统模块，它是Simulink的中心模块，是系统建模的核心和主要部分；

输出模块：即信号显示模块，它能够以图形方式、文件格式进行显示， 也可以在MATLAB的工作空间显示，输出模块主要集中在Sinks库。在Simulink中，模块都是用矢量来表示这三个部分的，如图1 所示

u (input)

x (states)

y (output)

2. Simulink仿真过程

Simulink仿真分为两个阶段：初始化和模型执行。 （1） 初始化阶段

初始化阶段需要完成的主要工作及其步骤如下：

 对模型的参数进行估计，得到它们实际计算的值。  展开模型的各个层次。

 按照更新的次序对模块进行排序。

 确定那些非显示化的信号属性，并检查每个模块是否能够接受连接

到它们输入端的信号。

 确定所有非显示化的信号采样时间模块的采样时间。

 分配和初始化存储空间，以便存储每个模块的状态和当前值的输

出。

（2） 模型执行阶段

—13—

黑龙江八一农垦大学毕业设计（论文） 对于一般的仿真模型是通过采用数值积分来来进行仿真的，计算数积分可以采用以下两步来进行：

 按照秩序计算每个模块的积分。

 根据当前输入和状态来决定状态的微分，得到微分矢量，然后把它返回给解法器，以计算下一个采样点的状态矢量。

2.4 本章小结

本章主要介绍MATLAB和Simulink两种仿真技术的概念，仿真原理。在通信系统中的应用。两种仿真的技术的强大功能以及今后的发展。

—14—

黑龙江八一农垦大学毕业设计（论文） 3 PSK 调制系统

在通信和信息传输系统、工业自动化或电子工程技术中，调制和解调应用最为广泛。而调制和解调的基本原理是利用信号与系统的频域分析和傅里叶变换的基本性质，将信号的频谱进行搬移，使之满足一定需要，从而完成信号的传输或处理。调制与解调又分模拟和数字两种，在现代通信中，调制器的载波信号几乎都是正弦信号，数字基带信号通过调制器改变正弦载波信号的幅度、频率或相位，产生幅度键控（ASK）、相位键控（PSK）、频率键控（FSK）信号，或同时改变正弦载波信号的几个参数，产生复合调制信号。本课程设计主要介绍基于Matlab对2PSK 和4PSK进制的调制仿真实现.

3.1 2PSK数字调制原理

3.1.1 2PSK数字调制

2PSK信号用载波相位的变化来表征被传输信息的状态，通常规定0相位载波和π相位载波分别表示传“1”和传“0”。

设二进制单极性码为an,其对应的双极性二进制码为bn,则2PSK信号的一般时域信号可以表示为：

S2psk（t）= [bn g(t-nTs)]cosωct

n式中 bn=-1（当an=0时，概率为P） bn=1（当an=1时，概率为1-P）’ 则时域信号可以变为:

S2psk（t）= [g(t-nTs)]cos（ωct+π）， 当an=0时

n S2psk（t）= [g(t-nTs)]cos（ωct+0） 当an=1时

n由此可知2PSK信号是一种双边带信号，功率谱为：

P2PSK（ƒ）= ƒs=P（1-P）[|G（ƒ+ ƒs）|+|G（ƒ- ƒs）|] +

1 ƒs2（1-P）2|G（0）|2[δ(ƒ+ ƒs)+ δ(ƒ- ƒs)] 42

2

2PSK信号的带宽为B2PSK=（ƒc+Rs）-（ƒc-Rs）= 2Rs

式中Rs为码元速率。

值得注意的是，2PSK码元序列的波形与载频和码元持续时间之间的关系有关。当一个码元中包含有整数个载波周期时，在相邻码元的边界处波形是不连续的，或者说相位是不连续的。当一个码元中包含的载波周期数比整数个周

—15—

黑龙江八一农垦大学毕业设计（论文） 期多半个周期时，则相位连续。当载波的初始相位差90度时，即余弦波改为正弦波时，结果类似。以上说明，相邻码元的相位是否连续与相邻码元的初始相位是否相同不可混为一谈。只有当一个码元中包含有整数个载波周期时，相邻码元边界处的相位跳变才是由调制引起的相位变化[16]。

2PSK信号的产生方法主要有两种。第一种叫相乘法，是用二进制基带不归零矩形脉冲信号与载波相乘，得到相位反相的两种码元。第二种方法叫选择法，是用此基带信号控制一个开关电路，以选择输入信号，开关电路的输入信 号是相位相差的同频载波。这两种方法的复杂程度差不多，并且都可以用数字信号处理器实现。

S（t） 双极性 不归零 eo（t） 码变换 相 乘 载波 2PSK及2DPSK的调制方框图

3.1.1用M文件编程实现

程序见附录一

—16—

黑龙江八一农垦大学毕业设计（论文）

3.1.3 用Simulink实现PSK调制

首先构造如图20所示的Simulink模型：

这里载波的参数设置见图下图，为了便于观察将载波频率设置为：

—17—

黑龙江八一农垦大学毕业设计（论文）

—18—

黑龙江八一农垦大学毕业设计（论文） 这样，所需模块找到以后，排列好，用简单的直线连接起来，点击运行，然后双击示波器scope，便出现所需的PSK调制波形。这里对示波器进行参数设置，使之同时显示二进制序列，载波波形和调制信号波形，见图

3.2 4PSK的调制和解调

3.2.1 4PSK信号

四进制绝对相移键控(4PSK)直接利用载波的四种不同相位来表示数字信息。4PSK信号相位φn矢量图如下

—19—

黑龙江八一农垦大学毕业设计（论文） 135o 01 10 90o 45o 11 00 0o 11 180o 参考相位 00 225o 01 270o

10 315o 参考相位 由于每一种相位代表两个比特信息，因此每个四进制码元可以用两个二进制码元的组合来表示。两个二进制码元中的前一比特用a来表示，后一比特用b表示，则双比特ab与载波相位的关系入下图： 双比特ab与载波相位的关系

双比特码元 载波相位（φn） A方式 B方式 0o

a 0 1 1 0 b 0 0 1 1 225o 90o 315 o 45 o 135 o 180o 270o 4PSK信号可以表示为

e4PSK(t)=

ng(t-nTs)cos(ωct+φn)

式中，g(t) 为信号包络波形，通常为矩形波，幅度为1：Ts为码元时间宽度；ωc 为角频率；φn 为第n个码元对应的相位，可取π/4,3π/4,5π/4,7π/4这四种值。

4PSK信号可以表示为 正交形式：

e4PSK(t)= [g(t-nTs)cosφn] cosωct - [g(t-nTs)sinφn ]sinωc

nn

= [ang(t-nTs)cosωct-bng(t-nTs) sinωct]

nn =I(t) cosωct-Q(t) sinωct

—20—

黑龙江八一农垦大学毕业设计（论文） 式中I(t)= Q(t)=

nang(t-nTs) bng(t-nTs)

n其中an ,bn分别是cosφn ，sinφn的值。an ,bn可取+1和-1。

可见，四进制信号可等效为两个正交载波进行双边带调制所得信号之和。这样，就把数字调相和线性调制联系起来，为四相波形的产生提供依据。

3.2.1 4PSK信号调制和解调

（1）4PSK调制原理：

4PSK的调制方法有正交调制方式（双路二相调制合成法或直接调相法）、相位选择法、插入脉冲法等。这里我们采用正交调制方式。

4PSK的正交调制原理如图。

a 单/双极性载波震荡 输入 × cosωct 串/并变移相π/2 + + - sinωct 单/双极性b 4PSK正交调制原理方框图 4PSK输出 ×

它可以看成是由两个载波正交的2PSK调制器构成的。图中串/并变换器将输入的二进制序列分为速度减半的两个并行双极性序列a和b（a,b码元在事件上是对齐的），再分别进行极性变换，把极性码变为双极性码（0→-1，1→+1）然后分别调制到cosωct和sinωct两个载波上，两路相乘器输出的信号是相互正交的抑制载波的双边带调制（DSB）信号，其相位与各路码元的极性有关，分别由a和b码元决定。经相加电路后输出两路的合成波形，即是4PSK信号。图中两个乘法器，其中一个用于产生0o与180o两种相位状态，另一个用于产

—21—

黑龙江八一农垦大学毕业设计（论文） 生90o与270o两种相位状态，相加后就可以得到45o，135o，225o，和315o四种相位状

（2）4PSK解调原理

4PSK信号是两个载波正交的2PSK信号的合成。所以，可以仿照2PSK相干检测法，用两个正交的相干载波分别检测两个分量 a和b，然后还原成二进制双比特串行数字信号。此法称作极性比较法（相干解调加码反变换器方式或相干正交解调发）。

4PSK信号解调器原理方图 yA(t) × 4PSK 输入 zA(t) 低通滤cosωct 正交载波源 sinωct × xA(t) 抽样判a 带通滤波yi(t) 位定并/串变器 b 低通滤zB(t) xB(t) 抽样判yB(t)

在不考虑噪声及传输畸变时，接收机输入的4PSK信号码元可表示为 yi(t)=A cos(ωct+φn)

式中φn为45o，135o，225o，315o四个相位值。 带通滤波器输出的两路信号 yA(t)= yB(t)= yi(t) 两路相乘器输出分别为

ZA（t）=A cos(ωct+φn) cosωct=

AAcos(2ωct+φn)+ cosφn 22AA sin(2ωct+φn)+ sinφn 22ZB（t）= A cos(ωct+φn)（-sinωct）=-

—22—

黑龙江八一农垦大学毕业设计（论文） 低通滤波器输出为

XA（t）=

Acosφn 2Asinφn 2XB（t）=

抽样判决器的判决准则如下表：

输入相位 φn 性 45o cosφn 的极性 n sinφ的极 判决器输出 a b

+ - - + + + - - 1 0 0 1 1 1 0 0 135 o o225 o 315 判决器是按

极性来判决的。即正抽样值判为1，负抽样值判为0.两路抽样判决器输出a、b，经并/串变换器就可将并行数据恢复成串行数据。

3.3 本章小结

本章讲述了数字调制的原理及二种调制方式相位键控（PSK、2PSK）的调制和解调方法和结果，通过仿真过程和结果来加深对三种键控的理解。

—23—

黑龙江八一农垦大学毕业设计（论文） 4 PSK调制解调系统的仿真

4.1 2PSK调制解调系统的仿真

系统方框图

各个模块的参数见附录（二） 运行结果如下：

—24—

黑龙江八一农垦大学毕业设计（论文）

4.2 4PSK调制解调系统的仿真

系统方框图

各个模块的参数见附录（三） 运行结果如下：

—25—

黑龙江八一农垦大学毕业设计（论文）

—26—

黑龙江八一农垦大学毕业设计（论文） 4.3利用MATLAB研究4PSK信号

系统方框图：

各个模块的参数见附录（四）

运行结果如下：

—27—

黑龙江八一农垦大学毕业设计（论文）

实际的传输系统不可能完全的满足无码间串扰传输条件，评价传输系统的一种定性而且简单的方法就是观察信号通过加有噪声的信道后的眼图上右图是信号没有通过任何的系统的眼图，下面是信号通过信噪比为26dB的高斯白噪声后的眼图。最下面是信号通过信噪比为16dB后的眼图。通过对比可以知道信噪比越大，信号的眼图越清晰，信噪比越小，眼图越模糊，没有噪声的话，眼图最清晰，近似线状。

眼图不仅反映信噪比对信号的影响，而且反映信号之间串扰的大小。眼图中改变参数设置，可以发现最佳判决时刻的改变。

—28—

黑龙江八一农垦大学毕业设计（论文）

4.4 本章小结

本章利用MATLAB中的SIMULINK研究了2PSK和４PSK通信系统，以及对其进行了仿真，对4PSK信号进行了研究。

—29—

结论 结论

本设计研究了2PSK和4PSK的调制和解调原理，以及利用MATLAB对其调制和解调进行了编程和编译仿真，得到的结论和理论上是一致的。简单而且快捷。同时利用MATLAB中的SIMYULINK对2PSK和4PSK的通信系统进行了仿真研究了其传输的特性，及传输中噪声对系统的影响。

本研究具有可对比性，对比2PSK和4PSK的通信原理和星座图可发现其中的不同点，但是频谱图近似相同。通信中信道的信噪比设置越大信噪传输越理想，与理论上是相符合的。2PSK和4PSK的传输系统也具有对比性，本研究在文中列出了仿真过程中每个元件的仿真参数的设置。比较其中不同点我们发现其中参数基本相似。也说明了他们的传输原理基本相同，都利用了相位的不同表示了不同的码元传输。

我们知道，2PSK信号是用载波的不同相位直接表示相应的数字信号而得出的，在这种绝对的移相过程中由于发端是以某一个相位作为基准的，因而在接收系统中也必须有这样一个固定的基准作为参考。如果这个参考相位发生变化，则回复的数字信息就会与发送的信息完全相反，从而造成错误的恢复。这种现象常称为2PSK的“倒”现象，因此，实际中一般不采用2PSK方式，而采用差分相移（2DPSK）方式。2DPSK是利用前后相邻码元的相对载波相位值去表示数字信息的一种方式。

2DPSK信号波形与2PSK的不同。2DPSK波形的 同一相位并不 对应相同的数字信息符号，而前后码元的的相对相位差才唯一的决定信息符号。这说明，解调2DPSK信号时并不依赖于某一固定的载波相位参考值。只要前后码元的相对相位关系不破坏则鉴别这个关系就可以正确恢复数字信息，避免了2PSK的方式中的“倒”现象的发生。但从波形上我们是无法分辨两者的不同，说明，一方面，只有已知移相键控方式是绝对还是相对的，才能正确的判断原信息；另一方面，相对移相信号可以看成是把数字信息序列（绝对码）变换成相对码，然后再根据相对码进行绝对相移而形成的。 2DPSK和2PSK的解调方式相同。

综上所述，用计算机仿真电子通信系统，具有广泛的适应性和极高的灵活性。在硬件试验中改变硬件的参数设置就意味着重做硬件，而在软件中只需对相应的参数进行重新设置，同时利用Simulink的可视化建模仿真和MATLAB简单编程的特点，可以实现较为复杂的系统，因此MATLAB/Simulink在通信系统仿真方面具有强大的功能和优越性！

—30—

附录 参考文献

[1]樊昌信 徐炳祥等.通信原理(第5版).北京.国防工业出版社.2005 [2]王福昌 熊兆飞 黄本雄.通信原理.清华大学出版社.2006

[3]钟麟 王峰.MATLAB仿真技术与应用教程.国防工业出版社.2003

[4]刘敏 魏玲.《MATLAB通信仿真与技术应用》.国防工业出版社.2001 [5]孙屹 吴磊.《Simulink通信仿真开发手册》.国防工业出版社，2004 [6]王兴亮.《数字通信原理与技术》（第二版）.西安电子科技大学出版社.2000 [7]沈越泓 高媛媛 魏以民.《通信原理》.机械工业出版社.2004

[8]郭文彬 桑森.通信原理（基于MATLAB的计算机仿真）.北京邮电大学出版社.2006

[9]张志涌等.精通MATLAB6.5版.北京航空航天大学出版社.2003

[10]邓华等.MATLAB通信仿真及应用实例祥解.人民邮电出版社.2003 [11]李建新.现代通信系统分析与仿真—MATLAB通信工具箱.西安.西安电

子科技大学出版社.2000

[12]李仲令 曹世文 葛造坤.现代通信系统仿真及应用.电子科技大学出

版社.1998.

[13]孙屹.MATLAB通信仿真开发手册.北京.国防工业出版社.2005 [14]薛定宇 陈阳泉.基于MATLAB/Simulink的系统仿真技术与应用[M].北

京.清华大学出版社.2002

[15]韩利竹 王华.MATLAB电子仿真与应用.北京.国防工业出版社.2003 [16]黎洪松 张卫刚等.数字通信原理.西安电子科技大学出版社.2005 [17]Peter M.Grant.Digital Communications(Second Edition).China Machine

Press

[18]H.S.Black.Modulation Theory.1993,Chap.4 [19]Emmanuel C.Ifeachor Barrie W.Jervis.Digital Sigital

Processing.Publishing House Of Electronics Industry

[20]Richard G.Lyons.Understanding Digital Siginal Processing.Science Press [21]W.R.Bennett.Spectra of Quantized Signals.BSTJ,1984,27(3)

[22]M.R.Winkler.Chirp signals for Communications.WESCON Convention

Record ,1998

[23]Adam Lender. Correlative Digital Communication Techniques. IEEE

Transactions on Communication Technology, 1996

[24] H.M.Boettinger.The Telephone Book.New York:Riverwood publishers,

1997

—31—

附录 致谢

我选择了田老师的课题，作设计的过程中，我有许多不懂得地方，在老师的指导下我一步步的解决问题完成论文，在完成过程中老师指导我去怎么选择资料，如何去利用网络资源，在这个学习的过程中，我了解到MATLAB的实用价值，更深的理解数字调制技术的调制原理。对三种键控方式的调制原理也理解的更加透彻。这都是老师的功劳。老师不厌其烦的一次次为我解答思路及程序上方方面面的问题。对论文的要求及答辩时间的改动老师都及时地联系我，通知我。在作毕业设计中，没有老师的帮助我是不可能完成的。感谢老师对我的关心和帮助。经过我的努力和索老师的耐心指导毕业设计顺利按时完成，它是对我们把本科四年所学的理论知识运用到实践中的一次系统的检验。，通过这段时间的亲身经历，我感觉自己学到了：收集、整理资料、共同协作、分析及处理问题等许多方面的知识。我真诚感谢这期间老师给予我的全力帮助，细心指导以及对我的严格要求，是她在我遇到问题时，不辞辛苦帮我解决,感谢她在设计和任务安排上长时间的指导。在设计的过程中，我还从老师身上学到好多东西：索老师热心工作的精神感动了我们，她还用实际行动告诉我们在工作中要脚踏实地，在学术上要严谨，在思维上要活跃，在学业上要勤奋刻苦。还有，在设计过程中感谢给我提供微机等便利的条件，使得我们能够顺利的完成毕业设计。最后感谢各位评委给予批评指正。

—32—

附录 附录

附录一

clear all; close all;

fs=8e5;%抽样频率 fm=20e3;%基带频率 n=2*(6*fs/fm); final=(1/fs)*(n-1); fc=2e5; % 载波频率 t=0:1/fs:(final); Fn=fs/2;%耐奎斯特频率

%用正弦波产生方波

%========================================== twopi_fc_t=2*pi*fm*t; A=1; phi=0;

x = A * cos(twopi_fc_t + phi); % 方波 am=1; x(x>0)=am; x(x<0)=-1;

figure(1) subplot(321); plot(t,x);

axis([0 2e-4 -2 2]); title('基带信号'); grid on

car=sin(2*pi*fc*t);%载波 ask=x.*car;%载波调制

—33—

附录

subplot(322); plot(t,ask);

axis([0 200e-6 -2 2]); title('PSK信号'); grid on;

%===================================================== vn=0.1;

noise=vn*(randn(size(t)));%产生噪音 subplot(323); plot(t,noise); grid on;

title('噪音信号'); axis([0 .2e-3 -1 1]);

askn=(ask+noise);%调制后加噪 subplot(324); plot(t,askn);

axis([0 200e-6 -2 2]); title('加噪后信号'); grid on;

%带通滤波

%====================================================================== fBW=40e3; f=[0:3e3:4e5]; w=2*pi*f/fs; z=exp(w*j); BW=2*pi*fBW/fs;

a=.8547;%BW=2(1-a)/sqrt(a) p=(j^2*a^2); gain=.135;

Hz=gain*(z+1).*(z-1)./(z.^2-(p));

subplot(325); plot(f,abs(Hz)); title('带通滤波器');

—34—

附录 grid on;

Hz(Hz==0)=10^(8);%avoid log(0) subplot(326);

plot(f,20*log10(abs(Hz))); grid on;

title('Receiver -3dB Filter Response'); axis([1e5 3e5 -3 1]);

%滤波器系数

a=[1 0 0.7305];%[1 0 p]

b=[0.135 0 -0.135];%gain*[1 0 -1] faskn=filter(b,a,askn); figure(2) subplot(321); plot(t,faskn);

axis([0 100e-6 -2 2]); title('通过带通滤波后输出'); grid on;

cm=faskn.*car;%解调 subplot(322); plot(t,cm);

axis([0 100e-6 -2 2]); grid on;

title('通过相乘器后输出'); %低通滤波器

%================================================================== p=0.72;

gain1=0.14;%gain=(1-p)/2 Hz1=gain1*(z+1)./(z-(p));

subplot(323);

Hz1(Hz1==0)=10^(-8);%avoid log(0) plot(f,20*log10(abs(Hz1))); grid on;

title('LPF -3dB response'); axis([0 5e4 -3 1]);

—35—

附录

%滤波器系数

a1=[1 -0.72];%(z-(p)) b1=[0.14 0.14];%gain*[1 1]

so=filter(b1,a1,cm); so=so*10;%add gain

so=so-mean(so);%removes DC component subplot(324); plot(t,so);

axis([0 8e-4 -3.5 3.5]); title('通过低通滤波器后输出'); grid on; %Comparator

%====================================================== High=2.5; Low=-2.5;

vt=0;%设立比较标准 error=0;

len1=length(so); for ii=1:len1 if so(ii) >= 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)'),

附录二：

Random-Integer Genrrator （随机整数发生器）的主要参数 模块名称 Random-Integer Genrrator

位置 Communications Blockset\\Comm Sources

—36—

附录 参数名称 M-ry number（元数） Initial seed（初始化种子） Sample time(采样时间) Spectrum Scope(频谱仪)的主要参数 模块名称 Spectrum Scope 位置 DSP Blockset\\DSP sinks

参 数 名 称 Buffer size（缓存长度） Buffer overlap（缓存交叠） FFT length（FFT长度） （点）数） Scope position（显示器位置） Frequency units（频率单位） Frequency range （频率范围） Amplitude scaling（幅度刻度） 号采样时间一致） MinmumY-limit(Y轴最小刻度) MaxmumY-limit(Y轴最大刻度) 模块名称 Discrete-Time Scatter Diagram 位置 CommunicationsBlockset\\ Comm Sinks 参数名称 Sanmples per symbol Offset（samples） Points display New traces per display 模块名称 BPSK Modulator Baseband

位置 Communications Blockset\\Modulation\\Digital Bseband Modulation 参数名称 Phase offset Samples per symbol 参数值 Pi/4 16 —37—

参数值 2PSK 2 12345 0.001 参 数 值 1024 512 1024 Number ofspectral averages（谱（计算）平均64 Get(0，‘defaultfigureposition’) Hertz [-Fs/2...Fs/2] dB Inherit sample increment from input（与输入信使能 -30 15 Discrete-Time Scatter Diagram（离散使劲星座图仪）的主要参数

参数值 1 0 1000 100 BPSK Modulator Baseband（基带BPSK调制器）的主要参数

附录

AWGN Channel（加性高斯白噪声信道）的主要参数 模块名称 AWGN Channel

位置Communications Blockset\\Channels 参 数 名 称 Initial seed（初始化种子） Mode（模式） SNR（dB）（信噪比） Input signal power

Error Rate Calculation（误码率计算）的主要参数 模块名称 Error Rate Calculation

位置 CommunicationsBlockset\\ Comm Sinks 参数名称 Receive delay（接收延迟） Communication delay（计算延迟） Computational mode（计算模式） Output data（输出数据）

附录三：

Random-Integer Genrrator （随机整数发生器）的主要参数 模块名称 Random-Integer Genrrator

位置 Communications Blockset\\Comm Sources 参数名称 M-ry number（元数） Initial seed（初始化种子） Sample time(采样时间)

Spectrum Scope(频谱仪)的主要参数 模块名称 Spectrum Scope

—38—

参数值 BPSK 18233 Signal to noise ration（SNR）（信噪比） 16 1 参数值 BPSK 1 0 Entire frame Port 参数值 QPSK 4 12345 0.001 附录 位置 DSP Blockset\\DSP sinks

参 数 名 称 Buffer size（缓存长度） Buffer overlap（缓存交叠） FFT length（FFT长度） （点）数） Scope position（显示器位置） Frequency units（频率单位） Frequency range （频率范围） Amplitude scaling（幅度刻度） 号采样时间一致） MinmumY-limit(Y轴最小刻度) MaxmumY-limit(Y轴最大刻度) 模块名称 Discrete-Time Scatter Diagram 位置 CommunicationsBlockset\\ Comm Sinks 参数名称 Sanmples per symbol Offset（samples） Points display New traces per display 模块名称 QPSK Modulator Baseband

位置 Communications Blockset\\Modulation\\Digital Bseband Modulation 参数名称 Input type（输入类型） Phase offset(相位偏移) 参数值 Integer Pi/4 参数值 1 0 1000 100 -30 15 Get(0，‘defaultfigureposition’) Hertz [-Fs/2...Fs/2] dB 参 数 值 1024 512 1024 Number ofspectral averages（谱（计算）平均64 Inherit sample increment from input（与输入信使能 Discrete-Time Scatter Diagram（离散使劲星座图仪）的主要参数

QPSK Modulator Baseband（基带QPSK调制器）的主要参数

Samples per symbol（每符号采样数） 16 AWGN Channel（加性高斯白噪声信道）的主要参数 模块名称 AWGN Channel

位置Communications Blockset\\Channels 参 数 名 称 Initial seed（初始化种子） Mode（模式） 参数值 QPSK 18233 Signal to noise ration（SNR）（信噪比） —39—

附录 SNR（dB）（信噪比） 模块名称 Error Rate Calculation

16 Error Rate Calculation（误码率计算）的主要参数 位置 CommunicationsBlockset\\ Comm Sinks 参数名称 Receive delay（接收延迟） Communication delay（计算延迟） Computational mode（计算模式） Output data（输出数据）

附录四：

Random-Integer Genrrator （随机整数发生器）的主要参数 模块名称 Random-Integer Genrrator

位置 Communications Blockset\\Comm Sources 参数名称 M-ry number（元数） Initial seed（初始化种子） Sample time(采样时间) 模块名称 QPSK Modulator Baseband

位置 Communications Blockset\\Modulation\\Digital Bseband Modulation 参数名称 Input type（输入类型） Phase offset(相位偏移) Samples per time 模块名称 FIR Interpolation 位置 DSP Blockset\\ Filters 参数名称 数） Interpolation（插补系数） Framing（分帧） 缓存初始条件） —40—

参数值 BPSK 1 0 Entire frame Port 参数值 4 12345 0.001 QPSK Modulator Baseband（基带QPSK调制器）的主要参数

参数值 Integer Pi/4 1 FIR Interpolation（FIR插补器）的主要参数

参数值 FIR filter coefficients （FIR滤波器系rcosine(1,8,[],0.3,1)(升余弦滤波器) 8 Maintain Input frame size Output Buffer initial conditions（输出NaN 附录 AWGN Channel（加性高斯白噪声信道）的主要参数 模块名称 AWGN Channel

位置Communications Blockset\\Channels 参 数 名 称 SNR（dB）（信噪比） Input signal power Eye Diagram（眼图仪）的主要参数 模块名称Eye Diagram

位置Communications Blockset\\Comm Sinks 参数名称 参数值 8 0 2 40 Sanmples per symbol Offset（samples）（预置） Symbols per trace（每迹符号数） Traces displayed（显示的迹数） 数） Eye diagram to display （显示的眼图） In-phase and Quadrature （同相及正交） 参数值 26 1 New traces per display（每显示新迹10 —41—