实验三Matlab的数字调制系统仿真实验(参考)概论
成都理工大学实验报告
课程名称：数字通信原理
姓名：__________________学号：______________ 成绩：____ ___ 实验三Matlab的数字调制系统仿真实验（参考）
1 数字调制系统的相关原理
数字调制可以分为二进制调制和多进制调制，多进制调制是二进制调制的推广，主要讨论二进制的调制与解调，简单讨论一下多进制调制中的差分相位键控调制（M-DPSK）。
最常见的二进制数字调制方式有二进制振幅键控（2-ASK）、移频键控（2-FSK）和移相键控（2-PSK 和2-DPSK）。下面是这几种调制方式的相关原理。
1.1 二进制幅度键控（2-ASK）
幅度键控可以通过乘法器和开关电路来实现。载波在数字信号1 或0 的控制下通或断，在信号为1 的状态载波接通，此时传输信道上有载波出现；在信号为0 的状态下，载波被关断，此时传输信道上无载波传送。那么在接收端我们就可以根据载波的有无还原出数字信号的1 和0。
幅移键控法(ASK)的载波幅度是随着调制信号而变化的，其最简单的形式是，载波在二进制调制信号控制下通断，此时又可称作开关键控法(OOK)。多电平MASK调制方式是一种比较高效的传输方式，但由于它的抗噪声能力较差，尤其是抗衰落的能力不强，因而一般只适宜在恒参信道下采用。
2-ASK 信号功率谱密度的特点如下：
（1）由连续谱和离散谱两部分构成；连续谱由传号的波形g(t)经线性调制后决定，离散谱由载波分量决定；
（2）已调信号的带宽是基带脉冲波形带宽的二倍。
1.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.1（1）；第二种方法是用以个受基带脉冲控制的开关电路去选择两个独立频率源的振荡作为输出，如图1.1（2）。两者的区别是前者的相位是连续的，后者由于两个独立的频率源产生的两个不同频率的信号，故相邻码元的相位不一定是连续的。
2-FSK 信号的接受也分为相干和非相干接受两类。最常用的解调方法是采用的相干检测法，相干检测的具体解调电路是同步检波器，原理方框图如图8-8所示。图中两个带通滤波器的作用同于包络检波法，起分路作用。它们的输出分别与相应的同步相干载波相乘，再分别经低通滤波器滤掉二倍频信号，取出含基带数字信息的低频信号，抽样判决器在抽样脉冲到来时对两个低频信号的抽样值V 0（t ）和V 1(t)进行比较判决（判决规则同于包络检波法），即可还原出基带数字信号。如图1.2 所示。
图1.1（2）
图1.1（1）
1.3 二进制相移键控（2-PSK ）
在相移键控中，载波相位受数字基带信号的控制，如在二进制基带信号中为0 时，载波相位为0 或π，为1 时载波相位为π 或0。载波相位和基带信号有一一对应的关系，从而达到调制的目的。2-PSK 信号的功率密度有如下特点：
(1) 由连续谱与离散谱两部分组成； (2) 带宽是绝对脉冲序列的二倍；
(3) 与2ASK 功率谱的区别是当P ＝1/2 时，2PSK 无离散谱，而2ASK 存在离散谱。
2PSK 信号的产生方法主要也是两种。第一种是相乘法，用二进制不归零矩形脉冲信号与载波相乘，得到相位反相的两种码元，如图1.3（1）；第二种方法叫选择法，是用此基带信号控制一个开关电路，以选择输入信号，开关电路的输入信号是相位相差π的同频载波，如图1.3（2）。
由于2PSK 信号实际上是以一个固定初相的末调载波为参考的，因此，解调时必须有与此同频同相的同步载波。如果同步载波的相位发生变化，如0相位变为相
位或
相位变为0相位，则恢复的数字信息就会发生“0”变“1”或“1”变“0”，从而造
带通滤波f 0
带通滤波f 1
相乘
相乘
低通滤波
低通滤波 抽样判决 输入 输出V o (t)
V 1(t)
cos 0ω(t) cos 1ω(t)
y 0(t)
y 1(t)
定时脉冲
图1.2 2FSK 信号的相干接收原理方框
A(t)
s(t)
载波
移相π
图1.3（2）
相乘
s(t)
A(t)
图1.3（1）
本地载波
成错误的恢复。这种因为本地参考载波倒相，而在接收端发生错误恢复的现象称为“倒
”现象或“反向工作”现象。绝对移相的主要缺点是容易产生相位模糊，造成反
向工作。
2PSK 信号的解调方法是相干接受法。由于PSK 信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号，如图1.4所示。
1.4 二进制差分相移键控（2DPSK ）
二进制差分相移键控(2DPSK)二进制差分相移键控常简称为二相相对调相，记作2DPSK 。它不是利用载波相位的绝对数值传送数字信息，而是用前后码元的相对载波相位值传送数字信息。所谓相对载波相位是指本码元初相与前一码元初相之差。
与2PSK 的波形不同，2DPSK 波形的同一相位并不对应相同的数字信息符号，而前后码元的相对相位才唯一确定信息符号。这说明解调2DPSK 信号时，并不依赖于某一固定的载波相位参考值，只要前后码元的相对相位关系不破坏，则鉴别这个相位关系就可正确恢复数字信息。这就避免了2PSK 方式中的“倒π”现象发生。
单从波形上看，2DPSK 与2PSK 是无法分辩的，一方面，只有已知移相键控方式是绝对的还是相对的，才能正确判定原信息；另一方面，相对移相信号可以看作是把数字信息序列（绝对码）变换成相对码，然后再根据相对码进行绝对移相而形成。这就为2DPSK 信号的调制与解调指出了一种借助绝对移相途径实现的方法。
2DPSK 信号的解调有两种解调方式，一种是差分相干解调，另一种是相干解调-码变换法。后者又称为极性比较－码变换法。
带通滤波
相乘
低通滤波
抽样判决
本地载波 提取
V(t)
图1.4 2PSK 信号相干接收原理方框
2 数字调制系统各个环节分析
2.1 仿真框图
典型的数字通信系统由信源、编码解码、调制解调、信道及信宿等环节构成，其框图如图2.1 所示：
数字调制是数字通信系统的重要组成部分，数字调制系统的输入端是经编码器编码后适合在信道中传输的基带信号。对数字调制系统进行仿真时，我们并不关心基带信号的码型，因此，我们在仿真的时候可以给数字调制系统直接输入数字基带信号，不用在经过编码器。
图2.1 数字通信系统模型
MATLAB 提供的图形界面仿真工具Simulink 由一系列模型库组成,包括Sources(信源模块)，Sinks(显示模块)，Discrete(离散系统模块)，Linear(线性环节)，Nonlinear(非线性环节)，Connections(连接),Blocksets&Toolboxes(其他环节)。特别是在Blocksets&Toolboxes 中还提供了用于通信系统分析设计和仿真的专业化模型库CommTbxLibrary。在这里，整个通信系统的流程被概括为：信号的产生与输出、编码与解码、调制与解调、滤波器以及传输介质的模型。在每个设计模块中还包含有大量的子模块，它们基本上覆盖了目前通信系统中所应用到的各种模块模型。通信系统一般都可以建立数学模型。根据所需仿真的通信系统的数学模型(或数学表达式)，用户只要从上述各个模型库中找出所需的模块,用鼠标器拖到模型窗口中组合在一起,并设定好各个模块参数, 就可方便地进行动态仿真.从输出模块可实时看到仿真结果,如时域波形图、频谱图等。每次仿真结束后还可以更改各参数,以便观察仿真结果的变化情况。另外，对Simulink 中没有的模块，可运用S 函数生成所需的子模块，并且可以封装和自定义模块库，以便随时调用。
根据Simulink 提供的仿真模块，数字调制系统的仿真可以简化成如图2.2 所示的模型：
图2.2 数字调制系统仿真框图
2.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 保持一致。
2.3 调制与解调模块
Simulink 通信工具箱中提供了数字信号各种调制方式的模块，如AM、CPM、FM 及PM 等。虽然不同的调制模块，参数设置有所不同，但很多参数在各种调制中是一致的，下面我们以DPSK 调制模块为例介绍一下调制模块的参数及其设置，其余模块将在下面仿真模型的建立过程中详细介绍。
M-DPSK Modulator Baseband和M-DPSK Demodulator Baseband 分别是数字信号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)]
2.4 信道
在分析通信系统时通常选择高斯噪声作为系统的噪声来考查，因为这种噪声在现实中比较常见而且容易分析。Simulink 中提供了带有加性高斯白噪声的信道：AWGN Channe。仿真时可以用该模块模拟现实中的信道，该模块的主要参数有：Initial seed：控制随机数产生的参数，要求不小于30，且与前面信号源中的Initial seed 设置不同的值；
Es/No (dB)：信号每个符号的能量与噪声的功率谱密度的比值；
SNR (dB)：信号功率与噪声功率的比值；
注：Es/No (dB) 和SNR (dB)是表征信号与噪声关系的两种方法，在一次仿真中只能选择其中一个。
2.5 误码计算仪
信号经过信道后，由于噪声的干扰，在接收端可能出现误码，Simulink中提供了Error Rate Calculation 模块来计算误码率，其主要参数的设置为：
Receive delay：接收延迟，表明在计算误码率时接收到的信号比源信号延迟的码元数，便于准确计算
Output data：数据输出，将误码率、误码数及码元总数输出，有两个选项可选择：Work space 和Port。将数据输出到Work space就是将误码率等数据存在内存中，以便下一步使用，而输出到Port中，则是在误码计算仪后面再接一个模块（比如结果显示模块），将数据传到该模块中（显示出来）；
Variable name：变量名称，该参数只有在前面选择了Work space后才有用，它决定数据输出到Wok space后的名称，默认值为ErrorVec。
2.6 示波器
在仿真过程中，可以观察各个环节的时域和频域波形，因此，在各个环节加上示波器以观察波形。另外，还可将示波器的数据输出到Work space 中存储，以便对仿真结果做进一步处理，比如将各个环节的波形对比显示和做频域变换等。
3 数字调制系统仿真及结果分析
3.1 仿真模型的设计及结果分析
了解了仿真所需的主要模块后，下一步就是设计和仿真各种数字调制模型，并对仿真结果在时域和频域进行分析。运用Matlab/Simulink软件进行各种数字通信系统的仿真，并对结果进行分析比较。Matlab/Simulink中有各种调制解调方式的现成模块可以直接应用，但是必须得设定正确的参数才能取得理想的仿真结果。通过仿真比较各种系统的占用带宽和误码率等性能指标。为了正确比较误码率，必须保证误码率比较器中比较的为发送码元和对应的接收码元，因此根据系统不同需要进行不同的延时；为了便于比较，建立了三个条件基本一致的仿真系统，即相同的信号源（伯努利随机二进制发生器），相同的传输环境（加性高斯白噪声环境信噪比为-5dB），都是基带调制、解调模块，最后显示出误码率、传输错误个数、传输总数。
3.2 二进制频移键控（2FSK）
利用伯努利二进制信号发生器（Bernoulli Binary Generator）成生数字信号，利用Simulink 通信工具箱中提供了专门的FSK 调制和解调模块应用FSK 调制模块能可方便地产生2FSK 信号，根据图3.1所示仿真框图。
模型中运用了Simulink 工具箱中的现成调制解调模块和信道模块，然后用示波器观察各环节波形，最后由误码计算仪计算误码。
图3.1 2FSK仿真模型
重要模块参数设置如下：
信源：probability of a zero : 0.5
Initial seed : 12345
Sample time : 0.001
调制解调器：m-ary number ：2
Frequency separation(hz) : 1000;2000;4000 Sample per symbol : 16 信道：initial seed : 12345
Mode : SNR SNR : -5 Input signal power : 1 误码器：receive delay : 0 Computation delay : 0 Mode : entire frame Output data : port 延时器：initial condition : 0 Sample time : -1
频谱：buffer input buffer size : 1024 buffer overlap : 512 window type : boxcar
Fft length : 1024 number of spectral average : 64 Frequency units : herts frequency range : [-fs/2…fs/2]
Amplitude acaling : db minimum y-limit : -30 maximum y-limit : 30
y-axis title : madnitude-squared,db
根据上述参数设置各模块，运行后，得出如下结果，分别设置FSK 调制解调器的频率设为1kHz 、2kHz 、4kHz 。系统的功率谱密度曲线如图3.2至图3.4所示。
图3.2 kHz f 1=?时2FSK 功率谱密度曲线
信源的码元速率为1000Baud ，对于2FSK 已调信号功率谱密度曲线来说，当
1Hz f k ?=时为单峰；当2Hz f k ?=时刚刚分离为双峰；当4Hz f k ?=为完全分离的双
峰。即，当f ?小于码元速率两倍时，2FSK 已调信号功率谱密度曲线为单峰，出现混叠现象；当f ?等于两倍码元速率时，恰好出现双峰；当f ?远远大于两倍码元速率时，功率谱密度曲线为双峰，频谱基本不再混叠。
图3.3 kHz f 2=?时2FSK
功率谱密度
图3.4 kHz f 4=?时2FSK 功率谱密度
3.3 最小频移键控（MSK ）
最小频移键控MSK 二进制频移键控（2FSK ）的改进，其Simulink 仿真模型如图3.5所示
将信号源的Sample time 设为1/2，其余参数可参照2-FSK ，两者参数类似。
从上图可以看出，与2FSK 调制频谱相比，MSK 信号的频谱比较紧凑，在主瓣之外，频谱旁瓣的下降非常迅速。这说明MSK 信号的功率主要包含在主瓣之内。因此，MSK 信号比较适合在窄带信道中传输，对邻道的干扰也较小。
3.4 二进制相移键控（2PSK ）
利用随机整数信号发生器（Random Integer Generator ）成生数字信号，利用Simulink 通信工具箱中提供了专门的BPSK 调制和解调模块，应用BPSK 调制模块能可方便地产生2PSK 信号，根据图3.7所示仿真框图。
图3.5 MSK 的仿真模型
图3.6 MSK 的仿真模型
重要模块参数设置如下： 信源：m-ary number ：2 Initial seed :12345 Sample time : 0.001 调制解调器：m-ary number ：2
Inputtype ：integer Phase rotation :pi/4
Output data type :double
信道、误码器、延时器、频谱器参数设置参考2FSK 。
根据上述参数设置各模块，运行后，得出如下结果，系统的功率谱密度曲线如图3.8所示。
在数字调制中，2PSK （后面将会看到2DPSK 也同样）的频谱特性与2ASK 十分相似。相位调制和频率调制一样，本质上是一种非线性调制，但在数字调相中，由于表征信息的相位变化只有有限的离散取值，因此，可以把相位变化归结为幅度
变化。这样一来，数字调相同线性调制的数字调幅就联系起来了，为此可以把数字
图3.8 2PSK 功率谱密度曲线
图3.7 2PSK 仿真模型
基于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上机实验答案
Matlab上机实验答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/180)/(1+exp(2)) z1 = >> x=[2 1+2i; 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = - + + -
>> a=::; >> z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2) (>> z33=(exp*a)-exp*a))/2.*sin(a++log(+a)/2)可以验证z3==z33，是否都为1） z3 = Columns 1 through 5 + + + + + Columns 6 through 10 + + + + + Columns 11 through 15 + + + + + Columns 16 through 20 + + + + +
Columns 21 through 25 + + + + + Columns 26 through 30 + + + + + Columns 31 through 35 + + + + + Columns 36 through 40 + + + + + Columns 41 through 45 + + + + + Columns 46 through 50
+ + + + + Columns 51 through 55 + + + + + Columns 56 through 60 + + + + + Column 61 + (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:: >> t=0::; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^ 2-2.*t+1) z4 =
基于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全部实验及答案
MATLAB全部实验及答案 实验一、MATLAB基本操作 实验内容及步骤 4、有关向量、矩阵或数组的一些运算 （1）设A=15；B=20；求C=A+B与c=a+b？ （2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B与 A.*B？ A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？ （4）设a=[1 -2 3;4 5 -4;5 -6 7] 请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全 下标的形式），并将其单下标转换成全下标。 clear,clc a=[1 -2 3;4 5 -4;5 -6 7]; [x,y]=find(a<0); c=[]; for i=1:length(x) c(i,1)=a(x(i),y(i)); c(i,2)=x(i); c(i,3)=y(i); c(i,4)=(y(i)-1)*size(a,2)+x(i); end c
（5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行A=[1,2;3,4]+i[5,6;7,8]，结果又如何？前面那 个是虚数矩阵，后面那个出错 （6）请写出完成下列计算的指令： a=[1 2 3;3 4 2;5 2 3]，求a^2=？,a.^2=？ a^2= 22 16 16 25 26 23 26 24 28 a.^2= 1 4 9 9 16 4 25 4 9 （7）有一段指令如下，请思考并说明运行结果及其原因 clear X=[1 2;8 9;3 6]; X( : ) 转化为列向量 （8）使用三元组方法，创建下列稀疏矩阵 2 0 8 0 0 0 0 1 0 4 0 0 6 0 0 0 方法一： clear,clc
数通实验报告二.数字调制
中南大学 通信原理实验报告书 题目：实验二 专业： 姓名： 学号： 时间：2014-12-13
通信原理实验报告（实验二） 实验名称：数字调制 一.实验目的 1、掌握绝对码、相对码概念及它们之间的变换关系。 2、掌握用键控法产生2ASK、2FSK、2DPSK信号的方法。 3、掌握相对码波形与2PSK信号波形之间的关系、绝对码波形与2DPSK信号波形之间的关系。 1、了解2ASK、2FSK、2DPSK信号的频谱与数字基带信号频谱之间的关系。 二.实验内容 1、用示波器观察绝对码波形、相对码波形。 2、用示波器观察2ASK、2FSK、2PSK、2DPSK信号波形。 3、用频谱仪观察数字基带信号频谱及2ASK、2FSK、2DPSK信号的频谱。 三.实验步骤 本实验使用数字信源单元及数字调制单元。 1、熟悉数字调制单元的工作原理。接通电源，打开实验箱电源开关。将数字调制单元单刀双掷开关K7置于左方N（NRZ）端。 2、用数字信源单元的FS信号作为示波器的外同步信号，示波器CH1接信源单元的(NRZ-OUT)AK（即调制器的输入），CH2接数字调制单元的BK，信源单元的K1、K2、K3置于任意状态（非全0），观察AK、BK波形，总结绝对码至相对码变换规律以及从相对码至绝对码的变换规律。 3、示波器CH1接2DPSK，CH2分别接AK及BK，观察并总结2DPSK信号相位变化与绝对码的关系以及2DPSK信号相位变化与相对码的关系（此关系即是2PSK信号相位变化与信源代码的关系）。注意：2DPSK信号的幅度比较小，要调节示波器的幅度旋钮，而且信号本身幅度可能不一致，但这并不影响信息的正确传输。 4、示波器CH1接AK、CH2依次接2FSK和2ASK；观察这两个信号与AK的关系（注意“1”码与“0”码对应的2FSK信号幅度可能不相等，这对传输信息是没有影响的）。 5、用频谱议观察AK、2ASK、2FSK、2DPSK信号频谱（条件不具备时不进行此项观察）。 四.实验过程及结果： 按实验步骤连接，得到全零码的AK、BK波形如下：
基于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)课后实验答案
实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012 2sin 851z e =+ (2) 21ln(2 z x =，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022 a a e e a z a a --+= ++=--L (4) 22 42011 122123t t z t t t t t ?≤
4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。解：(1) 结果： (2). 建立一个字符串向量例如： ch='ABC123d4e56Fg9';则要求结果是：
实验二 MATLAB 矩阵分析与处理 1. 设有分块矩阵33322322E R A O S ?????? =? ??? ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22 E R RS A O S +?? =???? 。 解: M 文件如下； 5. 下面是一个线性方程组： 1 231 1 12340.951110.673450.5211145 6x x x ?? ??????????????=??? ?????????????????? ? ch = 123d4e56g9
(1) 求方程的解。 (2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。 (3) 计算系数矩阵A 的条件数并分析结论。 解： M 文件如下： 实验三 选择结构程序设计 1. 求分段函数的值。 2226035605231x x x x y x x x x x x x ?+-<≠-? =-+≤<≠≠??--? 且且及其他 用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。 解：M 文件如下：
基于MATLAB模拟调制系统的仿真设计
1 线性模拟调制 1.1模拟调制原理 模拟调制是指用来自信源的基带模拟信号去调制某个载波，而载波是一个确知的周期性波形。模拟调制可分为线性调制和非线性调制，本文主要研究线性调制。 线性调制的原理模型如图1.1所示。设c(t)=Acos2t f o π,调制信号为m(t),已调信号为s(t)。 图1.1 线性调制的远离模型 调制信号m(t)和载波在乘法器中相乘的结果为：t A t m t s w o cos )()('=，然后通过一个传输函数为H(f)的带通滤波器，得出已调信号为。 从图1.1中可得已调信号的时域和频域表达式为： （1-1） 式（1-1）中，M(f)为调制信号m(t)的频谱。 由于调制信号m(t)和乘法器输出信号之间是线性关系，所以成为线性调制。带通滤波器H(f)可以有不同的设计，从而得到不同的调制种类。 1.2双边带调制DSB 的基本原理 在幅度调制的一般模型中，若假设滤波器为全通网络，调制信号m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，或称抑制载波双边带（DSB ）调制信号，简称双边带（DSB ）信号。 设正弦型载波c(t)=Acos( t) ，式中：A 为载波幅度， 为载波角频率。 根据调制定义，幅度调制信号（已调信号）一般可表示为： (t)=Am(t)cos(t) （1-2） ?? ???-++==) ()]()([21)()(*]cos )([)(f H f f M f f M f s t h t t m t s o o o w m(t) H(t) A os t w o c s(t) )(' t s
其中，m(t)为基带调制信号。 设调制信号m(t)的频谱为M()，则由公式2-2不难得到已调信号 (t)的频谱： )]()([2 )(c c m M M A s ωωωωω-++= （1-3） 由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域的简单搬移。 标准振幅就是常规双边带调制，简称调幅（AM ）。假设调制信号m(t)的平均值为0，将其叠加一个直流偏量 后与载波相乘，即可形成调幅信号。其时域表达式为: )cos())(()(0t t m t c AM A s ω+= （1-4） 式中： 为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。 若为确知信号，则AM 信号的频谱为: （1-5） AM 信号的频谱由载频分量、上边带、下边带三部分组成。AM 信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关，也就是说，载波分量并不携带信息。因此，AM 信号的功率利用率比较低。 AM 调制器模型如下图所示。 图1.2 AM 调制器模型 AM 信号的时域和频域表达式分别为 （1-6） （1-7） 式中，A o 为外加的直流分量；m(t)可以是确知信号也可以是随机信号，但通常认为其平均值为0，即0)(=t m — 。 由频谱可以看出，AM 信号的频谱由载波分量、上边带、下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像。因此，AM 信号是带有载波 分量的双边带信号，他的带宽是基带信号带宽 的2倍，即 ) (cos )()(cos ) (cos )]([)(t w c t m t w c A t w c t m A o t s o AM +=+=)]()([2 1)]()([)(w c w M w c w M w c w w c w A o t s AM -+++-++=δδπ)] ()([2 1)]()([)(0 ω ω ω ω ωωωδωδπωc c c c m M M A s -+++-++=f H
MATLAB实验题答案
1、求以下变量的值，并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A
matlab所有实验及答案
实验二 习题 1、 矩阵Y= ???? ? ???? ???3472123100451150425 ，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。取出子矩阵?? ? ? ??21301，并求该矩阵的维数。 解：命令为： Y=[5,2,4;0,15,1;45,100,23;21,47,3] Y(2,3) Y(10) sub2ind([4 3],2,3) [i,j]=ind2sub([4 3],10) find(Y==100) sub2ind([4 3],3,2) B=Y(2:2:4,3:-2:1) 或 B=Y([2 4],[3 1]) [m n]=size(Y) 2、 建立一个数值范围为0—100内4*5的整数随机矩阵。 求出大于50的元素的位置。 解：命令为： G=int8(100*rand(4,5)) find(G>50) 3、 已知矩阵A=[1 0 -1 ；2 4 1； -2 0 5]，B=[0 -1 0；2 1 3；1 1 2] 求2A+B 、A 2-3B 、A*B 、B*A 、A .*B ，A/B 、A\B 解：命令为： A=[1 0 -1 ；2 4 1； -2 0 5] B=[0 -1 0；2 1 3；1 1 2] E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B 4、 利用函数产生3*4阶单位矩阵和全部元素都为8的 4*4阶矩阵，并计算两者的乘积。 解：命令为： A=eye(3,4) B=8*ones(4) C=A*B 5、 创建矩阵a=????? ???????------7023021.5003.120498601 ,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b 构成矩阵d ，计算a*b 、c 实验二数字调制 一、实验目的 1、掌握绝对码、相对码概念及它们之间的变换关系。 2、掌握用键控法产生2ASK、2FSK、2DPSK信号的方法。 3、掌握相对码波形与2PSK信号波形之间的关系、绝对码波形与2DPSK信号波形之间的关系。 1、了解2ASK、2FSK、2DPSK信号的频谱与数字基带信号频谱之间的关系。 二、实验内容 1、用示波器观察绝对码波形、相对码波形。 2、用示波器观察2ASK、2FSK、2PSK、2DPSK信号波形。 3、用频谱仪观察数字基带信号频谱及2ASK、2FSK、2DPSK信号的频谱。 三、基本原理 本实验用到数字信源模块和数字调制模块。信源模块向调制模块提供数字基带信号（NRZ码）和位同步信号BS（已在实验电路板上连通，不必手工接线）。调制模块将输入的绝对码AK（NRZ码）变为相对码BK、用键控法产生2ASK、2FSK、2DPSK信号。调制模块内部只用+5V电压。 数字调制单元的原理方框图如图2-1所示，电原理图如图2-2所示（见附录）。 图2-1 数字调制方框图 本单元有以下测试点及输入输出点： ? CAR 2DPSK信号载波测试点 ? BK 相对码测试点 ? 2DPSK 2DPSK信号测试点/输出点，V P-P>0.5V ? 2FSK 2FSK信号测试点/输出点，V P-P>0.5V ? 2ASK 2ASK信号测试点，V P-P>0.5V 用2-1中晶体振荡器与信源共用，位于信源单元，其它各部分与电路板上主要元器件对应关系如下： ?÷2（A）U8：双D触发器74LS74 ?÷2（B）U9：双D触发器74LS74 ?滤波器A V6：三极管9013，调谐回路 ?滤波器B V1：三极管9013，调谐回路 ?码变换U18：双D触发器74LS74；U19：异或门74LS86 ? 2ASK调制U22：三路二选一模拟开关4053 ? 2FSK调制U22：三路二选一模拟开关4053 ? 2PSK调制U21：八选一模拟开关4051 ?放大器V5：三极管9013 ?射随器V3：三极管9013 将晶振信号进行2分频、滤波后，得到2ASK的载频2.2165MHZ。放大器的发射极和集电极输出两个频率相等、相位相反的信号，这两个信号就是2PSK、2DPSK的两个载波，2FSK 信号的两个载波频率分别为晶振频率的1/2和1/4，也是通过分频和滤波得到的。 下面重点介绍2PSK、2DPSK。2PSK、2DPSK波形与信息代码的关系如图2-3所示。 图2-3 2PSK、2DPSK波形 图中假设码元宽度等于载波周期的1.5倍。2PSK信号的相位与信息代码的关系是：前后码元相异时，2PSK信号相位变化180?，相同时2PSK信号相位不变，可简称为“异变同不变”。2DPSK信号的相位与信息代码的关系是：码元为“1”时，2DPSK信号的相位变化180?。码元为“0”时，2DPSK信号的相位不变，可简称为“1变0不变”。 应该说明的是，此处所说的相位变或不变，是指将本码元内信号的初相与上一码元内信号的末相进行比较，而不是将相邻码元信号的初相进行比较。实际工程中，2PSK或2DPSK 信号载波频率与码速率之间可能是整数倍关系也可能是非整数倍关系。但不管是那种关系，上述结论总是成立的。 本单元用码变换——2PSK调制方法产生2DPSK信号，原理框图及波形图如图2-4所示。相对于绝对码AK、2PSK调制器的输出就是2DPSK信号，相对于相对码、2PSK调制器的输出是2PSK信号。图中设码元宽度等于载波周期，已调信号的相位变化与AK、BK的关系当然也是符合上述规律的，即对于AK来说是“1变0不变”关系，对于BK来说是“异变同不变”关系，由AK到BK的变换也符合“1变0不变”规律。 图2-4中调制后的信号波形也可能具有相反的相位，BK也可能具有相反的序列即00100，这取决于载波的参考相位以及异或门电路的初始状态。 2DPSK通信系统可以克服上述2PSK系统的相位模糊现象，故实际通信中采用2DPSK而不用2PSK（多进制下亦如此，采用多进制差分相位调制MDPSK），此问题将在数字解调实验中再详细介绍。 实验四GUI 设计专业 学号姓名成绩电气工程及其自 动化201409140305 杨诚1.创建GUI 绘制方程c bx ax y ++=2图形，需要显示绘图结果的坐标系窗口，还能够输入 a ， b ， c 的值和x 取值范围（最大值和最小值）。 回调函数的编写： a=str2num(get(handles.edit1,'String'));b=str2num(get(handles.edit2,'String'));c=str2num(get(handles.edit3,'String'));xmin=str2num(get(handles.edit4,'String'));xmax=str2num(get(handles.edit5,'String'));x=xmin:0.1:xmax;y=a*x.^2+b*x+c;plot(x,y); 设计的运行界面截图： 2.设计一个GUI，完成画出y=sin(x)、y=cos(x)和y=x 3.的波形图。回调函数的编写： x=-5:0.1:5 plot(x,sin(x)) x=-5:0.1:5 plot(x,cos(x)) x=-5:0.1:5 y=x.^3 plot(x,y) 设计的运行界面截图： 3.创建一个GUI，含有下拉菜单，下拉菜单中有背景颜色选择。回调函数的编写： yanse=get(handles.popupmenu1,'value'); switch yanse case1 set(gcf,'color','r'); case2 set(gcf,'color','y'); case3 set(gcf,'color','g'); end 设计的运行界面截图： 实验一 模拟通信的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） matlab BPSK 调制与解调 1、调制 clear all; g=[1 0 1 0 1 0 0 1];%基带信号 f=100; %载波频率 t=0:2*pi/99:2*pi; cp=[];sp=[]; mod=[];mod1=[];bit=[]; for n=1:length(g); if g(n)==0; die=-ones(1,100); %Modulante se=zeros(1,100); % else g(n)==1; die=ones(1,100); %Modulante se=ones(1,100); % end c=sin(f*t); cp=[cp die]; mod=[mod c]; bit=[bit se]; end bpsk=cp.*mod; subplot(2,1,1);plot(bit,'LineWidth',1.5);grid on; title('Binary Signal'); axis([0 100*length(g) -2.5 2.5]); subplot(2,1,2);plot(bpsk,'LineWidth',1.5);grid on; title('ASK modulation'); axis([0 100*length(g) -2.5 2.5]); 2、调制解调加噪声 clc; close all; clear; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % 假定： % 2倍载波频率采样的bpsk信号 % 调制速率为在波频率的 N/2m % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% m=128; N=512; n=1:1:N; N0=0.5*randn(1,N) %噪声 h0=zeros(1,N); % 30阶低通滤波器 h0 f = [0 0.3 0.3 1]; w0 = [1 1 0 0]; b = fir2(30,f,w0); [h,w] = freqz(b,1,N/2); h0(1,1:N/2)=abs(h'); for i=1:N/2 h0(1,N-i+1)=h0(1,i); end; %%%%%%%%% 随机序列 a=rand(1,m); for i=1:m if(a(1,i)>0.5) a(1,i)=1; else a(1,i)=-1; end; end; %%% 生成BPSK信号 bpsk_m=zeros(1,N); j=1;k=1; for i=1:N if(j==(N/m+1)) j=1; k=k+1; end; % 0.05*pi 为初始相位，可以任意改变 bpsk_m(1,i)=a(1,k)*sin(2*pi*0.5*i+0.05*pi)+a(1,k)*cos(2*pi*0.5*i+ 0.05*pi); j=j+1; end; bpsk_m=bpsk_m+N0;% 信号加噪声，模拟过信道 % 接收处理用正交本振与信号相乘，变频 bpsk_m1=bpsk_m.*sin(2*pi*0.5*n); bpsk_m2=bpsk_m.*cos(2*pi*0.5*n); %滤波 tempx=fft(bpsk_m1); 实验二一维二维数组的创建和寻访 一、实验目的 1、掌握一维数组、二维数组创建和寻访的几种方法。 2、区别数组运算和矩阵运算的差别。 3、熟悉执行数组运算的常用数组操作函数。 4、掌握数组运算中的关系和逻辑操作及常用的关系、逻辑函数。 5、掌握“非数”、“空”数组在MA TLAB中的应用。 二、实验主要仪器与设备 装配有MA TLAB7.6软件的计算机 三、预习要求 做实验前必须认真复习第三章MATLAB的数值数组及向量化运算功能。 四、实验内容及实验步骤 1、一维数组的创建方法有哪几种？举例说明。 答：一维数组的创建方法有： ①递增/递减型一维数组的创建：冒号生成法：x=a：inc：b 线性（或对数）定点法：x=linspace(a,b,n)，x=logspace(a,b,n) ②逐个元素输入法：如x=[0.1，sin(pi/5)，-exp(-3)，-2*pi] ③运用MA TLAB函数生成法：例ones，rand等。 2、输入以下指令，并写出运行结果。本例演示：数组元素及子数组的各种标识和寻访格式；冒号的使用；end的作用。 A=zeros(2,6) %创建(2×6)的全零数组 A(:)=1:12 %赋值号左边：单下标寻访(2×6) 数组A的全部12个元素 %赋值号右边：拥有12个元素的一维数组 A(2,4) %双下标：A数组的第2行第4列元素 A(8) %单下标：数组A的第8个元素 A(: , [1,3]) %双下标：显示A的“第1列和第3列上全部行的元素” A([1, 2, 5, 6]') %单下标：把A数组第1,2,5,6个元素排成列向量 A(: , 4:end) %双下标：显示A的“从第4起到最后一列上全部行的元素” %在此end用于“列标识”，它表示“最后一列” A(2,1:2:5)=[-1, -3, -5] %把右边的3个数分别赋向A数组第2行的第1,3,5个元素位置 B=A([1, 2, 2, 2], [1, 3, 5]) %取A数组的1,3,5列的第1行元素作为B的第1行 %取A数组的1,3,5列的第2行分别作为B的第2,3,4行 L=A<3 %产生与A维数相同的“0，1”逻辑数组 A(L)=NaN %把逻辑1标识的位置上的元素赋为“非数” 运行结果： A = 0 0 0 0 0 0 0 0 0 0 0 0 Matlab FM调制仿真 目录 引言.................................................................................. 一．课程设计的目的与要求 .............................................. 1.1课程设计的目的.................................................... 1.2课程设计的要求.................................................... 二．FM调制解调系统设计............................................... 2.1FM调制模型的建立............................................. 2.2调制过程分析........................................................ 2.3FM解调模型的建立............................................. 2.4解调过程分析........................................................ 2.5高斯白噪声信道特性 ............................................ 2.6调频系统的抗噪声性能分析 ................................ 三．仿真实现...................................................................... 3.1MATLAB源代码.................................................. 3.2仿真结果................................................................ 四．心得体会...................................................................... 五．参考文献...................................................................... 引言 本课程设计用于实现DSB信号的调制解调过程。信号的调制与解调在通信系统中具有重要的作用。调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位 实验一：Matlab操作环境熟悉 一、实验目的 1．初步了解Matlab操作环境。 2．学习使用图形函数计算器命令funtool及其环境。 二、实验内容 熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算： 1．单函数运算操作。 求下列函数的符号导数 (1) y=sin(x); (2) y=(1+x)^3*(2-x); 求下列函数的符号积分 (1) y=cos(x); (2) y=1/(1+x^2); (3) y=1/sqrt(1-x^2); (4) y=(x-1)/(x+1)/(x+2); 求反函数 (1) y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2)); 代数式的化简 (1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4); (2) sin(x)^2+cos(x)^2; (3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x); 2．函数与参数的运算操作。 从y=x^2通过参数的选择去观察下列函数的图形变化 (1) y1=(x+1)^2 (2) y2=(x+2)^2 (3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6) y6=x^2/2 3．两个函数之间的操作 求和 (1) sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5 乘积 (1) exp(-x)*sin(x)实验二 数字调制
matlab实验四及其答案
实验一 模拟通信的MATLAB仿真
MATLAB仿真 BPSK调制
MATLAB实验题目及答案
matlabFM调制仿真
Matlab实验指导书(含答案)汇总