基于Matlab的数字调制系统仿真与分析本科生毕业论文

本科毕业学员毕业实践（论文、设计）报告论文题目：基于Mat lab的数字调制系统仿真与分析

毕业论文（设计）原创性声明

本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经发表或撰写过的研究成果。对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。

作者签名：日期：

毕业论文（设计）授权使用说明

本论文（设计）作者完全了解**学院有关保留、使用毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。学校可以公布论文（设计）的全部或部分内容。保密的论文（设计）在解密后适用本规定。

作者签名：指导教师签名：

日期：日期：

注意事项

1.设计（论文）的内容包括：

1）封面（按教务处制定的标准封面格式制作）

2）原创性声明

3）中文摘要（300字左右）、关键词

4）外文摘要、关键词

5）目次页（附件不统一编入）

6）论文主体部分：引言（或绪论）、正文、结论

7）参考文献

8）致谢

9）附录（对论文支持必要时）

2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。

3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。

4.文字、图表要求：

1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写

2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印

4）图表应绘制于无格子的页面上

5）软件工程类课题应有程序清单，并提供电子文档

5.装订顺序

1）设计（论文）

2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订

3）其它

摘要

数字调制是通信系统中最为重要的环节之一，数字调制技术的改进也是通信系统性能提高的重要途径。本文首先分析了数字调制系统的五种基本调制解调方法，然后，运用Matlab及附带的图形仿真工具——Simulink设计了这几种数字调制方法的仿真模型。通过仿真，观察了调制解调过程中各环节时域和频域的波形，并结合这几种调制方法的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。最后，在仿真的基础上分析比较了各种调制方法的性能，并通过比较仿真模型与理论计算的性能，证明了仿真模型的可行性。

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

ABSTRACT

In this paper, five usual methods of digital 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 digital modulation system and the reliability of the simulation models are analyzed. And then, the capability of three digital modulation simulation models, 2-FSK, 2-DPSK and MSK, have been compared, as well as comparing the results of simulation and theory.

Keywords: Digital modulation; analysis; simulation; MATLAB; SIMULINK.

目录

1引言 (1)

1.1数字调制系统概述 (1)

1.1.1数字通信系统的组成 (1)

1.1.2数字通信系统的特点 (2)

1.2数字调制的意义 (5)

1.3 Matlab在通信系统仿真中的应用 (6)

2数字调制系统的相关原理 (7)

2.1二进制幅度键控(2-ASK) (7)

2.2二进制频移键控(2-FSK) (7)

2.3二进制相移键控(2-PSK) (8)

2.4多进制数字调制 (8)

3数字调制系统的仿真设计 (9)

3.1数字调制系统各个环节分析 (9)

3.1.1仿真框图 (10)

3.1.2信号源仿真及参数设置 (11)

3.1.3调制与解调模块 (11)

3.1.4信道 (12)

3.2仿真模型的设计及结果分析 (13)

3.2.1 2-ASK (13)

3.2.2 2-FSK (18)

3.2.3 2-DPSK (22)

3.2.4 2-MSK (25)

3.2.5 M-DPSK (27)

3.3数字调制的性能比较 (30)

3.3.1各种仿真模型的性能比较 (30)

3.3.2仿真模型性能与理论性能的比较 (32)

4结论 (33)

致谢 (34)

参考文献 (35)

基于Matlab的数字调制系统仿真与分析

1引言

1.1数字调制系统概述

数字载波调制(简称数字调制)与模拟调制没有本质上的区别，它是用数字基带信号作为原始信号，去控制高频正弦载波信号的振幅、频率和相位，相应的有三种基本的调制方式：数字振幅调制（ASK）、数字频率调制（FSK）、数字相位调制（PSK）。

1.1.1数字通信系统的组成

数字通信系统就是利用数字信号传递消息的通信系统。而数字信号指的是不仅在时间上是离散的且在幅度上也是离散的信号。数字通信系统的形式各种各样，但从数字通信的特点以及所完成的功能上来看，可把它概括成图1-1所示的系统模型。

图1.1数字通信系统的组成框图

信源是信息的发出源，接其性质可分为离散信息源，如电报、数据等;另一种是连续的模拟信息源，如电路等。信宿是信息的归宿。

信源编码广义地说包含两个方面:(1)将输入信号变换成适合于数字通信系统处理和传输的数字信号。如果信源是模拟信号，应首先进行模拟/数字变换，使经过编码后的输出信号成为时间上离散、幅度取值有限，且

按一定规律组合的数字脉冲串:(2)通过信源编码提高数字信号的有效性，尽可能地减少原信号中的多余度，进行压缩信号的带宽的编码，使单位时间、单位系统频带上所传输的信息量最大.这两个方面是在信源编码的过程中同时完成的。信源译码则是信源编码的逆过程。

加密与解密是为了实现保密通信，通过加密人为地把待传输的数字序列扰乱。这种编码可采用周期非常长的伪随机码序列等，在接收端根据己知的解密方法，对接收序列进行解密。

信道编解码主要是为解决可靠性问题而设置的。由于信源编码后的数字信号是要通过信道来传输，而信号又不可避免地要受到各种噪声的干扰，因此可能会导致接收端数字信号的判决错误。信道编码就是采用一种对传输的原始信息按一定规则人为地加入一些数据，在接收端通过信道译码以达到自身发现和纠正误码的目的，这种技术称为“差错控制技术”。

一般来讲经信道编码的二元数字信号不适合在信道上直接传输，调制器的任务是把数字信号变为适合于信道传输的信号:而解调器的过程正好相反。通常对数字信号的频带调制有ASK, FSK和PSK等.调制与解调方式对通信质量的影响比较大，因此应合理选择。

定时同步系统使数字信号序列按节拍一步一步地工作，收、发两端的节拍一定要一致，否则将出现混乱。另外发送的数字信号序列常常是编组的，收端必须知道这些编组的头尾，否则就无法恢复原始信息。要保证收、发两端的节拍一致，必须有同步系统的控制。

1.1.2数字通信的特点

从数字通信的过程来看，很容易发现数字通信有着许多模拟通信无法比拟的优点。

(1) 抗干扰能力强

信号在传输过程中不可避免地要受到各种噪声的干扰。对于模拟信号来说，叠加在信号上的噪声难以与信号分开。同样，叠加在数字信号波形上的噪声也是难以去掉，但由于数字通信系统传送的数字信号，其信息并不包含在信号脉冲的波形之中，而是包含在脉冲的有无之中，因而只有当噪声在判决时超过某个范围，才有可能改变信号的值，产生错误判决，造成误码。因此，数字信号比模拟信号的抗干扰能力强，且数字信号还可进行纠错编码，进一步提高其抗干扰能力。由于数字信号的抗千扰能力强，在类似的信道条件下，数字通信的传输精度比模拟通信高的多。

(2) 采用再生中继可实现高质量远距离传输

远距离模拟通信系统中的噪声是积累的，因而随着通信距离的增加，传输质量也随之下降。而在数字通信系统中传送的数字信号大多是二元或三元信号。例如:二元数字信号只有两个状态“0”和“1"，在传输信道中也会受到噪声干扰，当干扰达到一定程度，进入再生中继器，再生中继器中的幅度识别电路对收到的二元波形信号进行判决，这些判决值通过波形形成电路，以没有噪声千扰的“纯净”脉冲波形向下一站发送，这样就消除了噪声的千扰。在理想情况下，噪声可全部清除，不会产生积累。正因为数字信号可以再生，所以可通过多个再生中继器实现高质量的远距离传输。

(3) 灵活性强适应各种业务要求

在数字通信中，各种消息(电报、电话、图像和数据等)均可变为统一的数字信号进行传输。在系统中对数字信号传输情况的监视信号、控制信号及业务信号均可采用数字信号。数字传输与数字交换技术结合起来组成的综合业务数字通信网(ISDN)，对于来自不同信息源的信号自动地进行变换、综合、传输、处理、存储和分离，实现各种综合业务，这给实际应用带来了极大的方便。

(4) 便于自动化、智能化

由于数字通信所传输的信号与数字电子计算机所采用的数字信号完全一致，所以可以方便地与计算机进行接口连接，构成复杂的、远距离的、大规模的、灵活多样的系统，使数字通信系统趋向于自动化、智能化。

(5) 易于加密

数字信号可用各种极其复杂规律的密码进行加密，只要用简单的逻辑电路就能实现，从而使通信具有高度的保密性。

(6) 设备易于集成化、微型化

由于数字通信系统中大都采用数字电路，数字电路比模拟电路易于集成化，因此通信设备可以采用中、大规模、超大规模集成电路制成体积小、功耗低、成本低、可靠性高的设备。

数字通信较模拟通信有如此多的优点是有代价的。传输数字信号所需带宽远比模拟信号的宽。但是随着宽频带传输媒质(光纤等)的广泛使用和频带压缩编码技术(如ADPCM、矢量编码、预测编码等)的日趋成熟和

实用化，数字通信占用频带宽的问题己得到解决，不再是数字通信发展的一个障碍.另一方面，数字通信的实现过程要复杂于模拟通信的实现过程，若是没有集成工艺作基础，很难在设备的体积、功耗、可靠性和经济性方面与发展已十分完美的模拟系统竞争，数字通信这种通信手段早在20世纪30年代就己提出，但它的真正发展和使用是从晶体管和集成电路的发明之后才开始的。特别是现在，大规模、超大规模集成电路的广泛使用，使数字通信得到了突飞猛进的发展和广泛的使用。

1.2数字调制的意义

数字调制是指用数字基带信号对载波的某些参量进行控制，使载波的这些参量随基带信号的变化而变化。根据控制的载波参量的不同，数字调制有调幅、调相和调频三种基本形式，并可以派生出多种其他形式。由于传输失真、传输损耗以及保证带内特性的原因，基带信号不适合在各种信道上进行长距离传输。为了进行长途传输，必须对数字信号进行载波调制，将信号频谱搬移到高频处才能在信道中传输。因此，大部分现代通信系统都使用数字调制技术。另外，由于数字通信具有建网灵活，容易采用数字差错控制技术和数字加密，便于集成化，并能够进入综合业务数字网（ISDN网），所以通信系统都有由模拟方式向数字方式过渡的趋势。因此，对数字通信系统的分析与研究越来越重要，数字调制作为数字通信系统的重要部分之一，对它的研究也是有必要的。通过对调制系统的仿真，我们可以更加直观的了解数字调制系统的性能及影响性能的因素，从而便于改进系统，获得更佳的传输性能。

1.3 Matlab在通信系统仿真中的应用

随着通信系统复杂性的增加,传统的手工分析与电路板试验等分析设计方法已经不能适应发展的需要,通信系统计算机模拟仿真技术日益显示出其巨大的优越性.。计算机仿真是根据被研究的真实系统的模型,利用计算机进行实验研究的一种方法.它具有利用模型进行仿真的一系列优点,如费用低,易于进行真实系统难于实现的各种试验,以及易于实现完全相同条件下的重复试验等。Matlab仿真软件就是分析通信系统常用的工具之一。

Matlab是一种交互式的、以矩阵为基础的软件开发环境,它用于科学和工程的计算与可视化。Matlab的编程功能简单,并且很容易扩展和创造新的命令与函数。应用Matlab可方便地解决复杂数值计算问题。Matlab 具有强大的Simulink动态仿真环境,可以实现可视化建模和多工作环境间文件互用和数据交换。Simulink支持连续、离散及两者混合的线性和非线性系统,也支持多种采样速率的多速率系统;Simulink为用户提供了用方框图进行建模的图形接口,它与传统的仿真软件包用差分方程和微分方程建模相比,更直观、方便和灵活。用户可以在Matlab和Simulink 两种环境下对自己的模型进行仿真、分析和修改。用于实现通信仿真的通信工具包(Communication toolbox,通信工具箱)是Matlab语言中的一个科学性工具包,提供通信领域中计算、研究模拟发展、系统设计和分析的功能,可以在Matlab环境下独立使用,也可以配合Simulink使用。另外，Matlab的图形界面功能GUI（Graphical User Interface）能为仿真系统生成一个人机交互界面，便于仿真系统的操作。因此，Matlab在通信系统仿真中得到了广泛应用，本文也选用该工具对数字调制系统进

行仿真。

2 数字调制系统的相关原理

数字调制可以分为二进制调制和多进制调制，多进制调制是二进制调制的推广，所以本文主要讨论二进制的调制与解调，最后简单讨论一下多进制调制中的差分相位键控调制（M-DPSK）。

最常见的二进制数字调制方式有二进制振幅键控（2-ASK）、移频键控（2-FSK）和移相键控（2-PSK和2-DPSK）。下面是这几种调制方式的相关原理。

2.1 二进制幅度键控（2-ASK）

幅度键控可以通过乘法器和开关电路来实现。载波在数字信号1或0的控制下通或断，在信号为1的状态载波接通，此时传输信道上有载波出现；在信号为0的状态下，载波被关断，此时传输信道上无载波传送。那么在接收端我们就可以根据载波的有无还原出数字信号的1和0。

2-ASK信号功率谱密度的特点如下：

（1）由连续谱和离散谱两部分构成；连续谱由传号的波形g(t)经线性调制后决定，离散谱由载波分量决定；

（2）已调信号的带宽是基带脉冲波形带宽的二倍。

2.2 二进制频移键控（2-FSK）

频移键控是利用两个不同频率f1和f2的振荡源来代表信号1和0，用数字信号的1和0去控制两个独立的振荡源交替输出。对二进制的频移键控调制方式，其有效带宽为B=2xF+2Fb,xF是二进制基带信号的带宽也是FSK信号的最大频偏，由于数字信号的带宽即Fb值大，所以二进制

频移键控的信号带宽B较大，频带利用率小。2-FSK功率谱密度的特点如下：

(1) 2FSK信号的功率谱由连续谱和离散谱两部分构成,?离散谱出现在f1和f2位置；

(2) 功率谱密度中的连续谱部分一般出现双峰。若两个载频之差|f1 -f2|≤fs，则出现单峰。

2.3二进制相移键控（2-PSK）

在相移键控中，载波相位受数字基带信号的控制，如在二进制基带信号中为0时，载波相位为0或π，为1时载波相位为π或0。载波相位和基带信号有一一对应的关系，从而达到调制的目的。2-PSK信号的功率密度有如下特点：

(1) 由连续谱与离散谱两部分组成；

(2) 带宽是绝对脉冲序列的二倍；

(3) 与2ASK功率谱的区别是当P＝1/2时，2PSK无离散谱，而2ASK存在离散谱。

2.4 多进制数字调制

上面所讨论的都是在二进制数字基带信号的情况，在实际应用中，我们常常用一种称为多进制（如4进制，8进制，16进制等）的基带信号。多进制数字调制载波参数有M种不同的取值，多进制数字调制比二进制数字调制有两个突出的优点：一是有于多进制数字信号含有更多的信息使频带利用率更高；二是在相同的信息速率下持续时间长，可以提高码元的能量，从而减小由于信道特性引起的码间干扰。现实中用得最多的

一种调制方式是多进制相移键控（MPSK）。

多进制相移键控又称为多相制，因为基带信号有M种不同的状态，所以它的载波相位有M种不同的取值，这些取值一般为等间隔。在多相制移键控有绝对移相和相对移相两种，实际中大多采用四相绝对移相键控（4PSK，有称QPSK），四相制的相位有0、π/2、π、3π/2四种，分别对应四种状态11、01、00、10。

3数字调制系统的仿真设计

3.1 数字调制系统各个环节分析

典型的数字通信系统由信源、编码解码、调制解调、信道及信宿等环节构成，其框图如图3.1所示：

数字调制是数字通信系统的重要组成部分，数字调制系统的输入端是经编码器编码后适合在信道中传输的基带信号。对数字调制系统进行仿真时，我们并不关心基带信号的码型，因此，我们在仿真的时候可以给数字调制系统直接输入数字基带信号，不用在经过编码器。

图3.1数字通信系统框图

3.1.1 仿真框图

MATLAB提供的图形界面仿真工具Simulink由一系列模型库组成,包括Sources(信源模块)，Sinks(显示模块)，Discrete(离散系统模块)，

Linear(线性环节)，Nonlinear(非线性环节)，Connections(连接),Blocksets&Toolboxes(其他环节)。特别是在Blocksets&Toolboxes 中还提供了用于通信系统分析设计和仿真的专业化模型库CommTbxLibrary。在这里，整个通信系统的流程被概括为：信号的产生与输出、编码与解码、调制与解调、滤波器以及传输介质的模型。在每个设计模块中还包含有大量的子模块，它们基本上覆盖了目前通信系统中所应用到的各种模块模型。通信系统一般都可以建立数学模型。根据所需仿真的通信系统的数学模型(或数学表达式)，用户只要从上述各个模型库中找出所需的模块,用鼠标器拖到模型窗口中组合在一起,并设定好各个模块参数, 就可方便地进行动态仿真.从输出模块可实时看到仿真结果,如时域波形图、频谱图等。每次仿真结束后还可以更改各参数,以便观察仿真结果的变化情况。另外，对Simulink中没有的模块，可运用S 函数生成所需的子模块，并且可以封装和自定义模块库，以便随时调用。

根据Simulink提供的仿真模块，数字调制系统的仿真可以简化成如图3.2所示的模型：

图3.2 数字调制系统仿真框图

3.1.2 信号源仿真及参数设置

Simulink通信工具箱中的Comm Sources/Data Sources提供了数字信号源Bernoulli Binary Generator，这是一个按Bernoulli分布提供随机二进制数字信号的通用信号发生器。在现实中，对受信者而言，发送端的信号是不可预测的随机信号。因此，我们在仿真中可以用Bernoulli Binary Generator来模拟基带信号发生器。

其中主要参数的含义为：

Probability of a zero ：产生的信号中0符号的概率，在仿真的时候一般设成0.5，这样便于频谱的计算；

Initial seed ：控制随机数产生的参数，要求不小于30，而且与后面信道中的Initial seed设置不同的值；

Sample time：抽样时间，这里指一个二进制符号所占的时间，用来控制号发生的速率，这个参数必须与后面调制和解调模块的Symbol period保持一致。

3.1.3 调制与解调模块

Simulink通信工具箱中提供了数字信号各种调制方式的模块，如AM、CPM、FM及PM等。虽然不同的调制模块，参数设置有所不同，但很多参数在各种调制中是一致的，下面我们以DPSK调制模块为例介绍一下调制模块的参数及其设置，其余模块将在下面仿真模型的建立过程中详细介绍。

M-DPSK Modulator Passband和M-DPSK Demodulator Passband 分别是数字信号DPSK调制和解调的专用模块，其中主要参数有：

M-ary number：输入信号的阶次数，比如2-DPSK就是2阶的；

Symbol period：符号周期，即，一个符号所占的时间，这必须与信号源的Sample time保持一致；

Carrier frequency：载波频率；

Carrier initial phase：载波的初始相位；

Input sample time：输入信号的抽样时间；

Output sample time：输出信号的抽样时间。

其中，各参数要满足以下关系：

Symbol period > 1/(Carrier frequency)

Input sample time < 1/[2*Carrier frequency + 2/(Symbol period) Output sample time <1/[2*Carrier frequency + 2/(Symbol period)]

3.1.4 信道

在分析通信系统时通常选择高斯噪声作为系统的噪声来考查，因为这种噪声在现实中比较常见而且容易分析。Simulink 中提供了带有加性高斯白噪声的信道：AWGN Channe。仿真时可以用该模块模拟现实中的信道，该模块的主要参数有：

Initial seed：控制随机数产生的参数，要求不小于30，且与前面信号源中的Initial seed设置不同的值；

Es/No (dB)：信号每个符号的能量与噪声的功率谱密度的比值；

SNR (dB)：信号功率与噪声功率的比值；

注：Es/No (dB) 和SNR (dB)是表征信号与噪声关系的两种方法，在一次仿真中只能选择其中一个。

3.2 仿真模型的设计及结果分析

了解了仿真所需的主要模块后，下一步就是设计和仿真各种数字调制模型，并对仿真结果在时域和频域进行分析。

3.2.1 2-ASK

通常，二进制振幅键控信号（2-ASK）的产生方法（调制方法）有两种，如图3.3所示：

图3.3 2-ASK信号产生的两种方法

2-ASK解调的方法也有两种相应的接收系统组成方框如图3.4所示：

图3.4 2-ASK信号接收系统组成框图

根据图3.3（a）所示方框图产生2-ASK信号，并用图3.4（b）所示的相干解调法来解调，设计2-ASK仿真模型如图3.5所示：

图3.5 2-ASK模型

在该模型中，调制和解调使用了同一个载波，目的是为了保证相干解调的同频同相，虽然这在实际运用中是不可能实现的，但是作为仿真，这样能获得更理想的结果。

主要模块参数设置如下：

1.Bernoulli Binary Generator的参数设置为：

Probability of a zero ：0.5

Initial seed ：67

Sample time：1

2.载波频率设为：50（可调）

3. Sample and Decide 模块是一个子系统，其内部结构由抽样和判决两部分组成，其中，抽样由同步冲激信号（Sychronizing signal）完成，其参数period（sec）设置和信号源的参数Sample time保持一致。判决模块是一个由M文件编写的S函数，S函数是Simulnk中用以功能扩展

的一个功能，用S函数可以自己编制Simulink库中没有的Simulink模块，从而使Simulink的功能大大加强，本模型中使用的判决模块就是这样一个应用。Sample and Decide 模块内部结构如图3.6所示：

图3.6 Sample and Decide 子系统内部结构

仿真结果时域分析

设信息源发出的是由二进制符号0、1组成的序列，且假定0符号出现的概率为P，1符号出现的概率为1-P，他们彼此独立。则，2ASK信号的时间表示式为

:

（3.2.）

将图3.5中各示波器的值输出到Work space中做统一处理，各环节波形如图3.7所示；

基于Matlab的FM仿真实现

摘要 本次设计主要是以Matlab为基础平台，对FM信号进行仿真。介绍了FM信号，及其调制和解调的基本原理，并设计M文件，分析在混入噪声环境下的波形失真，以及分析FM的抗噪声性能。本设计的主要目的是对Matlab的熟悉和对模拟通信理论的更深化理解。 关键词：Matlab；FM；噪声

前言 (2) 1 设计基础 (3) 1.1 Matlab及M文件的简介 (3) 1.2模拟调制概述 (4) 1.2.1模拟调制系统各个环节分析 (5) 1.2.2 模拟调制的意义 (6) 2 FM基本原理与实现 (7) 2.1 FM的基本原理 (7) 2.1.1调制 (7) 2.1.2解调 (8) 2.2 FM的实现 (8) 2.2.1 FM调制的实现 (8) 2.2.2 FM解调的实现 (9) 2.3 调频系统的抗噪声性能 (10) 2.3.1 高斯白噪声信道特性 (10) 3 FM的仿真实现与分析 (14) 3.1 未加噪声的FM解调实现 (14) 3.2 叠加噪声时的 FM解调 (16) 总结 (20) 致谢 (21) 参考文献 (22) 附录 (23)

通信按照传统的理解就是信息的传输。在当今高度信息化的社会，信息和通信已成为现代社会的命脉。信息作为一种资源，只有通过广泛传播与交流，才能产生利用价值，促进社会成员之间的合作，推动社会生产力的发展，创造出巨大的经济效益。而通信作为传输信息的手段或方式，与传感技术、计算机技术相融合，已成为21世纪国际社会和世界经济发展的强大动力。可以预见，未来的通信对人们的生活方式和社会的发展将会产生更加重大和意义深远的影响。 在通信系统中，从消息变换过来的原始信号所占的有效频带往往具有频率较低的频谱分量(例如语音信号)，如果将这种信号直接在信道中进行传输，则会严重影响信息传送的有效性和可靠性，因此这种信号在许多信道中均是不适宜直接进行传输的。在通信系统的发射端通常需要有调制过程，将调制信号的频谱搬移到所希望的位置上，使之转换成适于信道传输或便于信道多路复用的已调信号；而在接收端则需要有解调过程，以恢复原来有用的信号。调制解调方式常常决定了一个通信系统的性能。随着数字化波形测量技术和计算机技术的发展，可以使用数字化方法实现调制与解调过程。 调制在通信系统中具有重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号，而且它对系统的传输有效性和传输可靠性有着很大的影响。调制方式往往决定了一个通信系统的性能。调制技术是指把基带信号变换成传输信号的技术。基带信号是原始的电信号，一般是指基本的信号波形，在数字通信中则指相应的电脉冲。在无线遥测遥控系统和无线电技术中调制就是用基带信号控制高频载波的参数（振幅、频率和相位），使这些参数随基带信号变化。用来控制高频载波参数的基带信号称为调制信号。未调制的高频电振荡称为载波（可以是正弦波，也可以是非正弦波，如方波、脉冲序列等）。被调制信号调制过的高频电振荡称为已调波或已调信号。已调信号通过信道传送到接收端，在接收端经解调后恢复成原始基带信号。

基于MATLAB的OQPSK调制解调实现

基于MATLAB的OQPSK调制解调实现学生姓名：周翌指导老师：吴志敏 摘要本课程设计的目标在于深切理解OQPSK调制与解调的基本原理，学会使用MATALB软件中的M文件来实现OQPSK的调制与解调以及分析加入不同噪声时对信号的影响程度。首先产生一个数字基带信号，接下来调用MATLAB中的相应函数对这个基带信号进行调制，然后分析调制后的波形:，记录结果后对调制后的信号进行解调，观察解调结果并做好记录，最后在信号中加入噪声并观察其时频图的变化，分析信噪比的噪声对调制结果的影响。本课程设计的实验开发/运行平台为windowsXP/windows7，程序设计使用MATLAB语言。通过调试运行，基本完成设计目标,达到调制与解调的目的。 关键词：MATLAB；M文件；OQPSK；调制与解调；噪声 1 引言 数字调制与解调技术在数字通信中占有非常重要的地位，数字通信技术与MATLAB 的结合是现代通信系统发展的一个必然趋势。在数字信号通信过程中，噪声的影响往往比较大，同时我们都希望有较高的频带利用率和功率利用率，而OQPSK也是一种恒包络调制技术，其频谱特性好，既保留着2PSK的高抗噪声性能、高频带利用率和高功率利用率，又有效地减弱了2PSK的“反相工作”缺陷，在通信研究中有着非常重要的意义，特别是在卫星通信和移动通信的领域有着广泛的应用。MATLAB作为当前国际控制界最流行的面向工程与科学计算的高级语言，在控制系统的分析、仿真与设计方面得到了非常广泛的应用，随着其信号处理专业函数和专业工具箱的成熟，越来

越受到通信领域人士的欢迎，其在通信领域的应用也将更加广泛。 1.1课程设计目的 熟悉OQPSK的基本原理，掌握MATLAB中M文件的使用及相关函数的调用方法，在此基础上通过编程实现OQPSK的调制与解调，并通过加入的噪声来判断所设计的系统性能。这次课程设计不仅让我对OQPSK有了更加深入的了解，而且学会了如何利用MATLAB中的M文件来实现通信系统方面的应用，最重要的是，自己能够独立完成一个小项目了，有了这方面的经验，我在以后的学习中就会有更充足的信心和动力。 1.2课程设计要求 熟悉MATLAB中M文件的使用方法，并在深切理解OQPSK调制解调原理的基础上，编写出OQPSK调制解调程序。绘制出OQPSK信号解调前后在时域和频域中的波形，并观察解调前后频谱有何变化以加深对OQPSK信号解调原理的理解。分别对信号叠加不同噪声，并进行解调，绘制出解调前后信号的时频波形，分析不同噪声对信号传输造成的影响大小。 1.3课程设计步骤 先产生随机信号，然后对信号进行调制和解调，在调制和解调过程中加入高斯白噪声，观察现象。 1、产生四进制数字作为数字基带信号，对其进行调制； 2、将函数调制信号改为相应的时域波形调制信号； 3、在函数调制信号中加入高斯白噪声，生成加入噪声后的时域波形调制信号； 4、分别生成没加或加了噪声的调制信号波形图和频谱图； 5、分别对没加或加了噪声的调制信号进行解调； 6、计算误码率。

基于MATLAB的模拟线性调制

基于MATLAB 的模拟调制实验报告 一、实验目的 1.进一步学习调制的知识，掌握调频与调角两种模拟调制技术。 2.进一步学习MATLAB 的编程，熟练使用MATLAB 进行作图。 二、实验原理 1.调制的概念 调制（modulation ）就是对信号源的信息进行处理加到载波上，使其变为适 合 于信道传输的形式的过程，是使载波随信号而改变的技术。 一般，用来传送消息的信号()t u c 叫作载波或受调信号，代表所欲传送消息的信 号叫作调制信号，调制后的信号()t u 叫作已调信号。用调制信号()t u Ω控制载波的某些参数，使之随()t u Ω而变化，就可实现调制。 2.调制的目的 频谱变换 当所要传送的信号的频率或者太低，或者频带很宽，对直接采用电磁波的形 式进行发送很不利，需要的天线尺寸很大，而且发射和接受短的天线与谐振回路的参数变化范围很大。为了信息有效与可靠传输，往往需要将低频信号的基带频谱搬移到适当的或指定的频段。这样可以提高传输性能，以较小的发送功率与较短的天线来辐射电磁波。 实现信道复用 为了使多个用户的信号共同利用同一个有较大带宽的信道，可以采用各种复用技术。如模拟电话长途传输是通过利用不同频率的载波进行调制。将各用户话音每隔4 kHz 搬移到高频段进行传输。 提高抗干扰能力 不同的调制方式，在提高传输的有效性和可靠性方面各有优势。如调频广播系统，它采用的频率调制技术，付出多倍带宽的代价，由于抗干扰性能强，其音质比只占10 kHz 带宽的调幅广播要好得多。扩频通信就是以大大扩展信号传输带宽，以达到有效抗拒外部干扰和短波信道多径衰落的特殊调制方式。 3.调制的种类 根据()t u Ω和()t u c 的不同类型和完成调制功能的调制器传递函数不同，调制分为以下多种方式： （1）．按调制信号()t u Ω的类型分为： ● 模拟调制：调制信号()t u Ω是连续变化的模拟量，如话音与图像信号。 ● 数字调制：调制信号是数字化编码符号或脉冲编码波形。 （2）．按载波信号()t u c 的类型分： ● 连续波调制：载波信号为连续波形，通常以正弦波作为载波。

基于Matlab的数字水印设计——基于DCT域的水印实现

摘要 数字水印（Digital Watermark）技术是指用信号处理的方法在数字化的多媒体数据中嵌入隐蔽的标记，这种标记通常是不可见的，只有通过专用的检测器或阅读器才能提取。数字水印是信息隐藏技术的一个重要研究方向。随着数字水印技术的发展，数字水印的应用领域也得到了扩展，数字水印的基本应用领域是版权保护、隐藏标识、认证和安全不可见通信。 当数字水印应用于版权保护时，潜在的应用市场在于电子商务、在线或离线地分发多媒体内容以及大规模的广播服务。数字水印用于隐藏标识时，可在医学、制图、数字成像、数字图像监控、多媒体索引和基于内容的检索等领域得到应用。数字水印的认证方面主要ID卡、信用卡、ATM卡等上面数字水印的安全不可见通信将在国防和情报部门得到广泛的应用。 本文主要是根据所学的数字图象处理知识，在MATLAB环境下，通过系统编程的方式，建立并实现基于DCT域的数字水印加密系统。该系统主要包含数字水印的嵌入与提取，仿真结果表明，数字水印算法具有有效性、可靠性、抗攻击性、鲁棒性和不可见性，能够为数字媒体信息在防伪、防篡改、认证、保障数据安全和完整性等方面提供有效的技术保障。 关键词：数字水印；MATLAB；DCT

目录 1 课程设计目的 (1) 2 课程设计要求 (2) 3 数字水印技术基本原理 (3) 3.1 数字水印基本框架 (3) 3.2 算法分类 (3) 3.2.1 DCT法 (4) 3.2.2 其他方法 (4) 3.3 实际需要考虑的问题 (4) 3.3.1 不可见性 (4) 3.3.2 鲁棒性 (5) 3.3.3 水印容量 (5) 3.3.4 安全性 (5) 4 基于DCT变换仿真 (6) 4.1 算法原理 (6) 4.1.1 准备工作 (6) 4.1.2 选取8*8变换块 (7) 4.1.3 边界自适应 (7) 4.1.4 DCT变换与嵌入 (7) 4.1.5 恢复空域 (8) 4.2 嵌入算法扩展 (8) 4.2.1 RGB彩色图像三个矩阵的划分 (8) 4.2.2 八色彩色水印 (8) 4.3 水印的提取 (9) 4.4 仿真程序 (9) 5 结果分析 (14) 结束语 (16) 参考文献 (17)

基于MATLAB的模拟调制系统仿真与测试(AM调制)

闽江学院 《通信原理设计报告》 题目：基于MATLAB的模拟调制系统仿真与测试学院：计算机科学系 专业：12通信工程 组长：曾锴（3121102220） 组员：薛兰兰（3121102236） 项施旭（3121102222） 施敏（3121102121） 杨帆（3121102106） 冯铭坚（3121102230） 叶少群（3121102203） 张浩（3121102226） 指导教师：余根坚 日期：2014年12月29日——2015年1月4日

摘要在通信技术的发展中，通信系统的仿真是一个重点技术，通过调制能够将信号转化成适用于无线信道传输的信号。 在模拟调制系统中最常用最重要的调制方式是用正弦波作为载波的幅度调制和角度调制。在幅度调制中，文中以调幅、双边带和单边带调制为研究对象，从原理等方面阐述并进行仿真分析；在角度调制中，以常用的调频和调相为研究对象，说明其调制原理，并进行仿真分析。利用MATLAB下的Simulink工具箱对模拟调制系统进行仿真，并对仿真结果进行时域及频域分析，比较各个调制方式的优缺点，从而更深入地掌握模拟调制系统的相关知识，通过研究发现调制方式的选取通常决定了一个通信系统的性能。 关键词模拟调制；仿真；Simulink 目录 第一章绪论 (1) 1.1 引言 (1) 1.2 关键技术 (1) 1.3 研究目的及意义 (2) 1.4 本文工作及内容安排 (2) 第二章模拟调制原理 (3) 2.1 幅度调制原理 (3) 2.1.1 AM调制 (4) 第三章基于Simulink的模拟调制系统仿真与分析 (6) 3.1 Simulink工具箱简介 (6) 3.2 幅度调制解调仿真与分析 (8) 3.2.1 AM调制解调仿真及分析 (8) 第四章总结 (12) 4.1 代码 (13) 4.2 总结 (14)

基于MATLAB的QPSK通信系统仿真设计毕业设计论文

基于MATLAB的QPSK通信系统仿真设计 摘要 随着移动通信技术的发展，以前在数字通信系统中采用FSK、ASK、PSK 等调制方式，逐渐被许多优秀的调制技术所替代。本文主要介绍了QPSK调制与解调的实现原理框图，用MATLAB软件中的SIMULINK仿真功能对QPSK调制与解调这一过程如何建立仿真模型，通过对仿真模型的运行，得到信号在QPSK 调制与解调过程中的信号时域变化图。通过该软件实现方式，可以大大提高设计的灵活性，节约设计时间，提高设计效率，从而缩小硬件电路设计的工作量，缩短开发周期。 关键词 QPSK，数字通信，调制，解调，SIMULINK -I-

Abstract As mobile communications technology, and previously in the adoption of digital cellular system, ASK, FSK PSK modulation, etc. Gradually been many excellent mod ulation technology substitution, where four phase-shift keying QPSK technology is a wireless communications technology in a binary modulation method. This article prim arily describes QPSK modulation and demodulation of the implementation of the prin ciple of block diagrams, focuses on the MATLAB SIMULINK software emulation in on QPSK modulation and demodulation the process how to build a simulation model, through the operation of simulation model, I get signal in QPSK modulation and dem odulation adjustment process domain change figure. The software implementation, ca n dramatically improve the design flexibility, saving design time, increase efficiency, design to reduce the workload of hardware circuit design, and shorten the developmen t cycle. Keywords QPSK, Digital Communication,modulation,demodulation,SIMULINK -II-

实验三matlab的数字调制系统仿真实验(参考)

成都理工大学实验报告 课程名称：数字通信原理 姓名：__________________学号：______________ 成绩：____ ___实验三 Matlab的数字调制系统仿真实验（参考） 1 数字调制系统的相关原理 数字调制可以分为二进制调制和多进制调制，多进制调制是二进制调制的推广，主要讨论二进制的调制与解调，简单讨论一下多进制调制中的差分相位键控调制（M-DPSK）。 最常见的二进制数字调制方式有二进制振幅键控（2-ASK）、移频键控（2-FSK）和移相键控（2-PSK 和2-DPSK）。下面是这几种调制方式的相关原理。 二进制幅度键控（2-ASK） 幅度键控可以通过乘法器和开关电路来实现。载波在数字信号1 或0 的控制下通或断，在信号为1 的状态载波接通，此时传输信道上有载波出现；在信号为0 的状态下，载波被关断，此时传输信道上无载波传送。那么在接收端我们就可以根据载波的有无还原出数字信号的1 和0。 幅移键控法(ASK)的载波幅度是随着调制信号而变化的，其最简单的形式是，载波在二进制调制信号控制下通断，此时又可称作开关键控法(OOK)。多电平MASK 调制方式是一种比较高效的传输方式，但由于它的抗噪声能力较差，尤其是抗衰落的能力不强，因而一般只适宜在恒参信道下采用。 2-ASK 信号功率谱密度的特点如下： （1）由连续谱和离散谱两部分构成；连续谱由传号的波形g(t)经线性调制后决

定，离散谱由载波分量决定； （2）已调信号的带宽是基带脉冲波形带宽的二倍。 二进制频移键控（2-FSK） 数字频率调制又称频移键控（FSK），二进制频移键控记作2FSK。数字频移键控是用载波的频率来传送数字消息，即用所传送的数字消息控制载波的频率。2FSK信号便是符号“1”对应于载频f1，而符号“0”对应于载频f2（与f1不同的另一载频）的已调波形，而且f1与f2之间的改变是瞬间完成的。 从原理上讲，数字调频可用模拟调频法来实现，也可用键控法来实现。模拟调频法是利用一个矩形脉冲序列对一个载波进行调频，是频移键控通信方式早期采用的实现方法。2FSK键控法则是利用受矩形脉冲序列控制的开关电路对两个不同的独立频率源进行选通。键控法的特点是转换速度快、波形好、稳定度高且易于实现，故应用广泛。 频移键控是利用两个不同频率f1 和f2 的振荡源来代表信号1 和0，用数字信号的1 和0 去控制两个独立的振荡源交替输出。对二进制的频移键控调制方式，其有效带宽为B=2xF+2Fb,xF 是二进制基带信号的带宽也是FSK 信号的最大频偏，由于数字信号的带宽即Fb 值大，所以二进制频移键控的信号带宽B 较大，频带利用率小。2-FSK 功率谱密度的特点如下： (1) 2FSK 信号的功率谱由连续谱和离散谱两部分构成,?离散谱出现在f1 和f2 位置； (2) 功率谱密度中的连续谱部分一般出现双峰。若两个载频之差|f1 -f2|≤fs，则出现单峰。2FSK信号的产生方法主要是两种。第一种是用二进制基带矩形脉冲信号区调制一个调频器，使其能够输出两个不同的频率的码元，如图（1）；第二种方法是用以个受基带脉冲控制的开关电路去选择两个独立频率源的振荡作为输出，如图（2）。两者的区别是前者的相位是连续的，后者由于两个独立的频率源产生的两个不同频率的信号，故相邻码元的相位不一定是连续的。

数字水印技术及基于MATLAB的快速实现

2011年3月刊计算机工程应用技术信息与电脑 China Computer&Communication 1. 引言 多媒体及网络的迅速发展使得多媒体信息的交流和传输变得更 加简单和快捷，然而，这也使盗版者能以低廉的成本复制及传播未经 授权的数字产品，这种对数字产品保护和信息安全的迫切需求，导致 了数字水印技术成为多媒体信息安全领域的一个热点问题。数字水印 技术是在不影响宿主媒体主观质量的情况下，在宿主媒体（文本、图 像、视频、音频）中嵌入不易被人察觉的标识信息，用以证明原创作 者对其作品的所有权，并作为鉴定、起诉非法侵权的证据。 2. 数字水印的特征 一般认为数字水印应具有以下特征： (1) 安全性。数字水印应该是安全、难以被篡改的。当数字作品 发生变化时，数字水印应当也相应发生变化；同时，未经授权的个人 不得修改水印，理论上是未经授权的用户不能检测到产品中是否含有 水印。 (2) 鲁棒性。当被保护的数据在经过攻击后，嵌入的水印信息仍 能保持好的完整性并能以一定的正确概率被检测到。这些可能的攻击 包括噪声、滤波、剪切、旋转和编码等。 (3) 不可感知性。数字水印的嵌入不应使得原始作品发生可以感 知的变化，也不能使得被保护数据在质量上发生可以感觉到的失真。 (4) 可证明性。在多媒体作品的实际应用过程中可能需要多次加 入水印，这时水印系统必须能够允许水印被多次嵌入到被保护的数 据，而且每个水印均能独立地被证明。 (5) 无歧义性。恢复出的水印或对水印判决结果能够表明版权的 惟一，不会发生多重版权纠纷问题。 3. 数字水印的基本原理 通用的水印技术包含两个方面：水印的嵌入和水印的提取或检 测，如图1和图2所示。 图1 水印信号嵌入 图2 水印信号提取或检测 4. 数字水印的研究现状 4.1 文本水印 文本水印就是将代表著作人身份的信息(水印)嵌入到电子出版物 中，在产生版权纠纷时来验证版权的归属。其主要分为三大类：基于 文档结构的水印方法、基于自然语言处理技术的水印方法、基于传统 图像的水印方法。 基于文档结构的各种水印方法都只是提留在文本的表层，无法抵 抗对于文本结构和格式的攻击，简单的重新录入攻击就能使之失效， 因此这些水印方法普遍存在鲁棒性差的缺点。自然语言文本水印方法 相对提高了抗攻击的能力，但普遍存在容量不足的问题。基于传统图 像的文本水印普遍存在鲁棒性不高、操作复杂的缺点。 4.2 图像水印 根据水印的实现过程，图像水印算法可分为空域算法和变换域算 法。空域算法是通过直接改变原始图像的像素值来嵌入水印，通常具 有较快的速度，但鲁棒性差，且水印容量也会受到限制；变换域算法 是通过改变某些变换系数来嵌入水印，通常具有很好的鲁棒性和不可 见性。其实现一般是基于图像变换，如DCT、DFT、DWT等。重点介 绍一下变换域算法。 4.2.1 离散傅里叶变换 (DFT) 该方法是利用图像的DFT来嵌入信息。通信理论中调相信号的抗 干扰能力比调幅信号的抗干扰能力强，同样在图像中利用相位信息嵌 入的水印也比用幅值信息嵌入的水印更稳健。实验表明该方法的抗压 缩能力比较弱。 4.2.2 离散余弦变换 (DCT) DCT能把空间域的图像转换到变换域上进行研究，从而能很容易 了解到图像的各空间频域成分，进行相应处理。基于DCT的水印方法 与基于DFT的水印方法相比有较好的鲁棒性，但是无法做到对图像信 号内容的自适应，因此往往会造成对图像特征的明显损害，不可感知 性不是最佳。 4.2.3 离散小波变换 (DWT) DWT是一种时间---频率信号的多分辨率分析方法，在时频两域 都具有表征信号局部特征的能力。实验表明，与DCT、DFT变换相比 较，基于DWT的水印算法的鲁棒性最优，且与JPEG2000、MPEG4压 缩标准兼容，利用DWT产生的水印具有良好的视觉效果和抵抗多种 攻击的能力，且不可感知性最好。 4.3 音频水印 音频水印利用音频文件的冗余信息和人耳听觉系统的特点来嵌入 水印，其可以保护声音数字产品不被随意复制和篡改，如CD唱片， 广播电台的节目内容等。音频水印的三种基本方法：扩频嵌入方 法、回声隐藏方法和相位编码方法。 4.4 视频水印 视频水印是通过对视频载体的时间和空间冗余来嵌入水印，其既 不影响视频质量，又能达到保护节目制作者的合法权益和控制数字产 品的复制。视频水印从算法要求上同图像水印有许多相似之处，但视 频水印也有一些独特之处，如能够在压缩和未压缩的格式下实时完成 水印的检测，对MPEG压缩、A/D和D/A转换等都有较好的稳健性。 数字水印技术涉及到通信理论、编码理论、噪声理论、视听觉 感知理论、扩频技术、信号处理技术、数字图像处理技术、多媒体技 术、模式识别技术、算法设计等理论，用到经典的DFT离散傅立叶变数字水印技术及基于MATLAB的快速实现 张 巍1 时宏伟2 （1.78179部队，四川成都 610011；2. 川大智胜，四川成都 610045） 摘要：数字水印是近几年来出现的数字产品版权保护技术，是当前国际学术界的研究热点．该文论述了数字水印的提出及研究现状、水印的基本原理和算法、水印的分类等情况，并介绍了一种可以快速上手的高效的实用语言——MATLAB，同时给出了一个用MATLAB工具在静止图像上嵌入水印的实例。 关键词：数字水印；MATLAB；DCT 中图分类号：TP39 文献标识码：A 文章编号：1003-9767（2011）03-0130-02

基于MATLAB的模拟信号频率调制(FM)与解调分析

课程设计任务书 学生姓名：杨刚专业班级：电信1302 指导教师：工作单位：武汉理工大学 题目：信号分析处理课程设计 －基于MATLAB的模拟信号频率调制（FM）与解调分析 初始条件： 1.Matlab6.5以上版本软件； 2.先修课程：通信原理等； 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1、利用MATLAB中的simulink工具箱中的模块进行模拟频率（FM）调制与解调，观 察波形变化 2、画出程序设计框图，编写程序代码，上机运行调试程序，记录实验结果（含计算结 果和图表等），并对实验结果进行分析和总结； 3、课程设计说明书按学校统一规范来撰写，具体包括： ⑴目录；⑵理论分析； ⑶程序设计；⑷程序运行结果及图表分析和总结； ⑸课程设计的心得体会（至少800字，必须手写。）； ⑹参考文献（不少于5篇）。 时间安排： 周一、周二查阅资料，了解设计内容； 周三、周四程序设计，上机调试程序； 周五、整理实验结果，撰写课程设计说明书。 指导教师签名： 2013 年 7月 2 日 系主任（或责任教师）签名： 2013年 7月 2日

目录 1 Simulink简介 (1) 1.1 Matlab简介······················································错误！未定义书签。 1.2 Simulink介绍 ···················································错误！未定义书签。 2 原理分析 ·····························································错误！未定义书签。 2.1通信系统 ·························································错误！未定义书签。 2.1.1通信系统的一般模型 ···································错误！未定义书签。 2.1.2 模拟通信系统 (3) 2.2 FM调制与解调原理···········································错误！未定义书签。 3 基于Matlab方案设计 (6) 3.1 Matlab代码 (6) 3.2 Matlab仿真 (8) 4 基于Simulink方案设计 (12) 4.1 使用Simulink建模和仿真的过程 (12) 4.1.1 Simulink模块库简介 (12) 4.1.2 调制解调模块库简介 (13) 4.2 FM调制与解调电路及仿真 (14) 4.3 仿真结果分析 (17) 5 心得体会 ·····························································错误！未定义书签。 6 参考文献 (20) 本科生课程设计评定表

数值分析的matlab实现

第2章牛顿插值法实现 参考文献：[1]岑宝俊. 牛顿插值法在凸轮曲线修正设计中的应用[J]. 机械工程师,2009,10:54-55. 求牛顿插值多项式和差商的MA TLAB 主程序： function[A,C,L,wcgs,Cw]=newpoly(X,Y) n=length(X);A=zeros(n,n);A(:,1) =Y'; s=0.0;p=1.0;q=1.0;c1=1.0; for j=2:n for i=j:n A(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1)); end b=poly(X(j-1));q1=conv(q,b);c1=c1*j;q=q1; end C=A(n,n);b=poly(X(n));q1=conv(q1,b); for k=(n-1):-1:1 C=conv(C,poly(X(k)));d=length(C);C(d)=C(d)+A(k,k); end L(k,:)=poly2sym(C);Q=poly2sym(q1); syms M wcgs=M*Q/c1;Cw=q1/c1; （1）保存名为newpoly.m 的M 文件 （2）输入MA TLAB 程序 >> X=[242,243,249,250]; >> Y=[13.681,13.526,13.098,13.095]; >> [A,C,L,wcgs,Cw]=newpoly(X,Y) 输出3阶牛顿插值多项式L 及其系数向量C 差商的矩阵A ，插值余项wcgs 及其 ) ()()1(ξ+n n f x R 的系数向量Cw 。 A = 13.6810 0 0 0 13.5260 -0.1550 0 0 13.0980 -0.0713 0.0120 0 13.0950 -0.0030 0.0098 -0.0003 C = 1.0e+003 *

基于Matlab 的QPSK调制解调仿真设计与研究

天津理工大学计算机与通信工程学院通信工程专业设计说明书 基于Matlab/Simulink 的QPSK调制解调仿真设计与研究 姓名韩双年 学号 20092226 班级 09通信3班 指导老师白媛 日期 2012-12-16

目录 摘要 (2) 第一章前言 (2) 1.1 专业设计任务及要求 (2) 1.2 Matlab简介 (2) 1.3 Matlab下的simulink简介 (3) 1.4 通信系统模型 (3) 第二章QPSK调制 (4) 2.1 QPSK介绍 (4) 2.2 QPSK调制原理 (4) 2.2.1 相乘法 (4) 2.2.2 选择法 (5) 2.3 QPSK调制原理框图 (6) 2.4 QPSK调制方式的Matlab仿真 (6) 2.5 QPSK调制方式Matlab-simulink仿真 (7) 2.5.1 simulink调制建模 (7) 2.5.2 simulink调制仿真结果 (8) 第三章QPSK解调 (13) 3.1 QPSK解调原理 (13) 3.2 QPSK解调原理框图 (13) 3.3 QPSK解调方式Matlab仿真 (13) 3.4 QPSK解调方式的Matlab-simulink仿真 (14) 3.4.1 QPSK解调建模 (14) 3.4.2 传输信道 (16) 3.4.3仿真结果 (16) 3.5 仿真结果分析 (18) 第四章QPSK通信系统性能分析 (19) 第五章结论 (19) 参考文献 (20) 附录 (20)

摘要 正交相移键控（QPSK），是一种数字调制方式。QPSK技术具有抗干扰能力好、误码率低、频谱利用效率高等一系列优点。论文主要介绍了正交相移键控(QPSK)的概况，以及正交相移键控QPSK的调制解调概念和原理，利用Matlab中M文件和Simulink模块对QPSK的调制解调系统进行了仿真，对QPSK在高斯白噪声信道中的性能进行了，分析了解Simulink中涉及到QPSK的各种模块的功能。 【关键词】Matlab QPSK Simulnk 仿真 第一章前言 1.1专业设计任务及要求 1了解并掌握QPSK调制与解调的基本原理； 2在通信原理课程的基础上设计与分析简单的通信系统； 3学会利用MATLAB7.0编写程序进行仿真，根据实验结果能分析所设计系统的性能。 4学习MATLAB的基本知识，熟悉MATLAB集成环境下的Simulink的仿真平台。 5利用通信原理相关知识在仿真平台中设计QPSK调制与解调仿真系统并用示波器观察解调后的波形 6在指导老师的指导下，独立完成课程设计的全部内容，能正确的阐述和分析设计和实验结果。 1.2 Matlab简介 MATLAB是MATrix LABoratory的缩写，是一款由美国Math Works公司出品的商业数学软件。MATLAB 是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。除了矩阵运算、绘制函数/数据图像等常用功能外，MATLAB还可以用来创建用户界面及与调用其它语言(包括C，C++和FORTRAN)编写的程序。尽管MATLAB主要用于数值计算,但是因为大量的额外

基于MATLAB的2ASK、2FSK和2PSK的调制仿真

实验报告（一） 一、实验名称：基于MATLAB 的2ASK 、2FSK 和2PSK 的调制仿真 二、实验目的： （1）熟悉2ASK 、2FSK 和2PSK 的调制原理。 （2）学会运用Matlab 编写2ASK 、2FSK 和2PSK 调制程序。 （3）会画出原信号和调制信号的波形图。 （4）掌握数字通信的2ASK 、2FSK 和2PSK 的调制方式。 三、实验原理分析 3.1二进制振幅键控（2ASK ） 振幅键控是利用载波的幅度变化来传递数字信息，而其频率和初始相位保持不变。在2ASK 中，载波的幅度只有两种变化状态，分别对应二进制信息“0”或“1”。二进制振幅键控的表达式为： s(t) = A(t)cos(w 0+θ) 0＜t ≤T 式中，w 0=2πf 0为载波的角频率；A(t)是随基带调制信号变化的时变振幅，即 A(t) = ?? ? 0A 典型波形如图所示： 2ASK 信号的产生方法通常有两种：相乘法和开关法，相应的调制器如图2。图2（a ） 就是一般的模拟幅度调制的方法，用乘法器实现；图2（b ）是一种数字键控法，其中的开关电路受s(t)控制。 在接收端，2ASK 有两种基本的解调方法：非相干解调（包络检波法）和相干解调（同步检测法），相应的接收系统方框图如图：

3.2、二进制频移键控（2FSK ） 二进制频移键控信号码元的“1”和“0”分别用两个不同频率的正弦波形来传送，而其振幅和初始相位不变。故其表达式为： =)(s t ???? ?++时 "0发送“),cos(”时1发送“),cos 21(?ω?ωn n t A t A 图4 2FSK 信号时间波形 由图可见，2FSK 信号的波形（a ）可以分解为波形（b ）和波形（c ）,也就是说，一个2FSK 信号可以看成是两个不同载频的2ASK 信号的叠加。 2FSK 信号的调制方法主要有两种。第一种是用二进制基带矩形脉冲信号去调制一个调频器，使其能够输出两个不同频率的码元。 第二种方法是用一个受基带脉冲控制的开关电路

数值分析的MATLAB程序

列主元法 function lianzhuyuan(A,b) n=input('请输入n:') %选择阶数A=zeros(n,n); %系数矩阵A b=zeros(n,1); %矩阵b X=zeros(n,1); %解X for i=1:n for j=1:n A(i,j)=(1/(i+j-1)); %生成hilbert矩阵A end b(i,1)=sum(A(i,:)); %生成矩阵b end for i=1:n-1 j=i; top=max(abs(A(i:n,j))); %列主元 k=j; while abs(A(k,j))~=top %列主元所在行 k=k+1; end for z=1:n %交换主元所在行a1=A(i,z); A(i,z)=A(k,z); A(k,z)=a1; end a2=b(i,1); b(i,1)=b(k,1); b(k,1)=a2; for s=i+1:n %消去算法开始m=A(s,j)/A(i,j); %化简为上三角矩阵 A(s,j)=0; for p=i+1:n A(s,p)=A(s,p)-m*A(i,p); end b(s,1)=b(s,1)-m*b(i,1); end end X(n,1)=b(n,1)/A(n,n); %回代开始 for i=n-1:-1:1 s=0; %初始化s for j=i+1:n s=s+A(i,j)*X(j,1);

end X(i,1)=(b(i,1)-s)/A(i,i); end X 欧拉法 clc clear % 欧拉法 p=10; %贝塔的取值 T=10; %t取值的上限 y1=1; %y1的初值 r1=1; %y2的初值 %输入步长h的值 h=input('欧拉法please input number(h=1 0.5 0.25 0.125 0.0625):h=') ; if h>1 or h<0 break end S1=0:T/h; S2=0:T/h; S3=0:T/h; S4=0:T/h; i=1; % 迭代过程 for t=0:h:T Y=(exp(-t)); R=(1/(p-1))*exp(-t)+((p-2)/(p-1))*exp(-p*t); y=y1+h*(-y1); y1=y; r=r1+h*(y1-p*r1); r1=r; S1(i)=Y; S2(i)=R; S3(i)=y; S4(i)=r; i=i+1; end t=[0:h:T]; % 红线为解析解，'x'为数值解 plot(t,S1,'r',t,S3,'x')

利用MATLAB实现QPSK调制及解调

郑州轻工业学院 课程设计说明书题目：利用MATLAB实现QPSK调制及解调 姓名： 院系：电气信息工程学院 专业班级：电子信息工程09-1 学号： 540901030154 指导教师：赵红梅 成绩： 时间： 2012 年 6 月 18 日至 2012 年 6 月 22 日

郑州轻工业学院 课程设计任务书 题目利用MATLAB实现QPSK调制及解调 专业班级电子信息工程09级 1班学号 54 姓名 主要内容、基本要求、主要参考资料等： 主要内容： 已知数字信号1011000101101011，码元速率为2400波特，载波频率为1200Hz，利用MATLAB画出QPSK调制波形，并画出调制信号经过高斯信道传输后解调波形及接收误码率，将其与理论值进行比较。 基本要求： 1、通过本课程设计，巩固通信原理QPSK调制的有关知识； 2、熟悉QPSK产生原理； 3、熟悉高斯信道的建模及QPSK解调原理； 4、熟悉误码率的蒙特卡罗仿真； 5、学会用MATLAB来进行通信系统仿真。 主要参考资料： 主要参考资料： 1、王秉钧等. 通信原理[M].北京：清华大学出版社，2006.11 2、陈怀琛.数字信号处理教程----MATLAB释义与实现[M].北京：电子工业出版社,2004. 完成期限：2012.6.18—2012.6.23 指导教师签名： 课程负责人签名： 2012年6月16日

目录 一前言 (4) 1.1QPSK系统的应用背景简介 (4) 1.2 QPSK实验仿真的意义 (4) 1.3 实验平台和实验内容 (5) 1.3.1实验平台 (5) 1.3.2实验内容 (5) 二、系统实现框图和分析 (5) 2.1、QPSK调制部分， (5) 2.2、QPSK解调部分 (7) 三、实验结果及分析 (7) 3.1、理想信道下的仿真 (7) 3.2、高斯信道下的仿真 (8) 3.3、先通过瑞利衰落信道再通过高斯信道的仿真 (9) 参考文献： (11) 附录 (12)

同济大学数值分析matlab编程题汇编

MATLAB 编程题库 1.下面的数据表近似地满足函数2 1cx b ax y ++=，请适当变换成为线性最小二乘问题，编程求最好的系数c b a ,,，并在同一个图上画出所有数据和函数图像. 625 .0718.0801.0823.0802.0687.0606.0356.0995 .0628.0544.0008.0213.0362.0586.0931.0i i y x ---- 解： x=[-0.931 -0.586 -0.362 -0.213 0.008 0.544 0.628 0.995]'; y=[0.356 0.606 0.687 0.802 0.823 0.801 0.718 0.625]'; A=[x ones(8,1) -x.^2.*y]; z=A\y; a=z(1); b=z(2); c=z(3); xh=-1:0.1:1; yh=(a.*xh+b)./(1+c.*xh.^2); plot(x,y,'r+',xh,yh,'b*')

2.若在Matlab工作目录下已经有如下两个函数文件，写一个割线法程序，求出这两个函数 10 的近似根，并写出调用方式： 精度为10 解： >> edit gexianfa.m function [x iter]=gexianfa(f,x0,x1,tol) iter=0; while(norm(x1-x0)>tol) iter=iter+1; x=x1-feval(f,x1).*(x1-x0)./(feval(f,x1)-feval(f,x0)); x0=x1;x1=x; end >> edit f.m function v=f(x) v=x.*log(x)-1; >> edit g.m function z=g(y) z=y.^5+y-1; >> [x1 iter1]=gexianfa('f',1,3,1e-10) x1 = 1.7632 iter1 = 6 >> [x2 iter2]=gexianfa('g',0,1,1e-10) x2 = 0.7549 iter2 = 8

实验一 模拟通信的MATLAB仿真

实验一 模拟通信的MATLAB 仿真 姓名：左立刚 学号：031040522 简要说明： 实验报告注意包括AM ，DSB ，SSB ，VSB ，FM 五种调制与解调方式的实验原理，程序流程图，程序运行波形图，simulink 仿真模型及波形，心得体会，最后在附录中给出了m 语言的源程序代码。 一．实验原理 1．幅度调制（AM ） 幅度调制（AM ）是指用调制信号去控制高频载波的幅度，使其随调制信号呈线性变化的过程。AM 信号的数学模型如图3-1所示。 图2-1 AM 信号的数学模型 为了分析问题的方便，令 δ =0， 1.1 AM 信号的时域和频域表达式 ()t S AM =[A 0 +m ()t ]cos t c ω （2-1） ()t S AM =A 0 π[()()ωωωωδC C ++-]+()()[]ωωωωc c M M ++-2 1 （2-2）

AM 信号的带宽 2 =B AM f H （2-3） 式中， f H 为调制信号的最高频率。 2.1.3 AM 信号的功率P AM 与调制效率 η AM P AM =()222 2 t m A +=P P m c + （2-4） 式中，P C =2 A 为不携带信息的载波功率；()2 2 t m P m =为携带信息的边带 功率。 ()() t t m A m P P AM C AM 2 2 2+= = η （2-5） AM 调制的优点是可用包络检波法解调，不需要本地同步载波信号，设备简单。AM 调制的最大缺点是调制效率低。 2.2、双边带调制（DSB ） 如果将在AM 信号中载波抑制，只需在图3-1中将直流 A 0 去掉，即可输出 抑制载波双边带信号。 2.2.1 DSB 信号的时域和频域表达式 ()()t t m t c DSB S ωcos = （2-6） ()()()[]ωωωωωC C DSB M M S ++-=2 1 （2-7） DSB 信号的带宽 f B B H AM DSB 2 == （2-8）