1、 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 1 页 共 25 页基 于 MATLAB/Simulink 的 PSK 频 带 传 输 系统 仿 真 与 性 能 分 析摘 要 本课程设计主要是设计一个 PSK 频带传输系统并对其进行仿真与性能分析。本课程设计仿真平台为 MATLAB/Simulink。在设计此频带传输系统时，首先对信号进行 PSK 调制，再通过加入高斯白噪声传输信道，接着在接收端对信号进行 PSK 解调，最后把输出的信号和输入的信号进行比较。通过最后仿真结果可知，该模拟信号频带传输通信系统已初步实现了设计指标并可用于解决一些实际性的问题。关键
2、词 PSK 调制解调；高斯白噪声；MATLAB/Simulink1 引 言通 信 (Communication)就 是 信 息 的 传 递 ， 是 指 由 一 地 向 另 一 地 进 行 信 息 的 传 输 与交 换 ， 其 目 的 是 传 输 消 息 。 然 而 ， 随 着 社 会 生 产 力 的 发 展 ， 人 们 对 传 递 消 息 的 要 求也 越 来 越 高 。 在 各 种 各 样 的 通 信 方 式 中 ， 利 用 “电 ”来 传 递 消 息 的 通 信 方 法 称 为 电信 (Telecommunication)， 这 种 通 信 具 有 迅 速 、 准 确 、 可 靠 等 特
3、点 ， 且 几 乎 不 受 时 间 、地 点 、 空 间 、 距 离 的 限 制 ， 因 而 得 到 了 飞 速 发 展 和 广 泛 应 用 。目前，无论是模拟通信还是数字通信，在不同的通信业务中都得到了广泛的应用。但是，数字通信的发展速度已明显超过了模拟通信，成为当代通信技术的主流。与模拟通信相比，数字通信具有以下一些优点：抗干扰能力强，且噪声不积累；传输差错可控；便于用现代数字信号处理技术对数字信息进行处理、变换、存储；易于集成，使通信设备微型化，重量轻；易于加密处理，且保密性好。数字通信的缺点是，一般需要较大的带宽。另外，由于数字通信对同步要求高，因而系统设备复杂。但是，随着微电子技术、
4、计算机技术的广泛应用以及超大规模集成电路的出现，数字系统的设备复杂程度大大降低。同时高效的数据压缩技术以及光纤等大容量传输媒质的使用正逐步使带宽问题得到解决。因此，数字通信的应用必将越来越广泛。本课程设计主要是设计一个 PSK 频带传输系统并对其进行仿真与性能分析 1。在设计此频带传输系统 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 2 页 共 25 页时，首先对输入信号利用相关的模块进行 PSK 调制，再通过加入高斯白噪声传输信道，接着在接收端对信号进行 PSK 解调，最后把输出的信号和输入的信号进行比较。1.1 设计目的本课程设计的目的主要是仿真通信系
5、统中频带传输技术中的 PSK 调制。首先产生一段随机的二进制非归零码的基带信号，对其进行 PSK 调制后再送入加性高斯白噪声（AWGN）信道传输，在接收端对其进行 PSK 解调以恢复原信号，观察还原是否成功，改变 AWGN 信道的信噪比，计算传输前后的误码率，绘制信噪比-误码率曲线，并与理论曲线比较进行说明。1.2 设计要求（1）本设计开发平台为 MATLAB 中的 Simulink。（2）模型设计应该符合工程实际，模块参数设置必须与原理相符合。（3）处理结果和分析结论应该一致，而且应符合理论。（4）独立完成课程设计并按要求编写课程设计报告书。1.3 设计平台MATLAB/Simulink1.
6、4 设计步骤（1）根据课题任务书要求结合教材上讲授的原理，在纸质文档上画出课题的电路级框图。（2）从设计平台上找到实现框图中每个方框对应的功能实现模块，并按图连接方式进行连接。（3）按照通信原理教材介绍的原理进行每个模块的参数设置。（4）运行模块，根据模型的运行结果或出错提示，不断地修改模型和仿真参数，直至仿真成功。 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 3 页 共 25 页2 基本原理2.1 设计平台简介MATLAB 集成环境下的 Simulink：MATLAB 是一种功能强大的科学计算和工程仿真软件，它的交互式集成界面能够帮助用户快速地完成数值分
7、析、矩阵运算、数字信号处理、仿真建模、系统控制和优化等功能。MATLAB 语言采用与数学表达相同的形式，不需要传统的程序设计语言，由于 MATLAB 的这些特性，它已经成为科研工作和工程仿真中的高效助手。Simulink 是 MATLAB 最 重 要 的 组 件 之 一 ， 它 提 供 一 个 动 态 系 统 建 模 、 仿 真 和综 合 分 析 的 集 成 环 境 。 在 该 环 境 中 ， 无 需 大 量 书 写 程 序 ， 而 只 需 要 通 过 简 单 直 观 的鼠 标 操 作 ， 就 可 构 造 出 复 杂 的 系 统 。 Simulink 具 有 适 应 面 广 、 结 构 和 流
8、 程 清 晰 及仿 真 精 细 、 贴 近 实 际 、 效 率 高 、 灵 活 等 优 点 ， 并 基 于 以 上 优 点 Simulink 已 被 广 泛应 用 于 控 制 理 论 和 数 字 信 号 处 理 的 复 杂 仿 真 和 设 计 。 同 时 有 大 量 的 第 三 方 软 件 和 硬件 可 应 用 于 或 被 要 求 应 用 于 Simulink。Simulink 的 储 存 模 块 ， 仿 真 数 据 可 以 方 便 地 以 各 种 形 式 保 存 到 工 作 区 或 文 件 中 ，供 用 户 在 仿 真 结 束 之 后 对 数 据 进 行 分 析 和 处 理 。 另 外 Si
9、mulink 把 具 有 特 定 功 能 的代 码 组 织 成 模 块 的 方 式 ， 并 且 这 些 模 块 可 以 组 织 成 具 有 等 级 结 构 的 子 系 统 ， 因 此 具有 内 在 的 模 块 化 设 计 要 求 。Simulink 是 MATLAB 中 的 一 种 可 视 化 仿 真 工 具 ， 是 一 种 基 于 MATLAB 的框 图 设 计 环 境 ， 是 实 现 动 态 系 统 建 模 、 仿 真 和 分 析 的 一 个 软 件 包 ， 被 广 泛 应 用 于 线性 系 统 、 非 线 性 系 统 、 数 字 控 制 及 数 字 信 号 处 理 的 建 模 和 仿 真
10、 中 。 Simulink 可 以用 连 续 采 样 时 间 、 离 散 采 样 时 间 或 两 种 混 合 的 采 样 时 间 进 行 建 模 ， 它 也 支 持 多 速 率系 统 ， 也 就 是 系 统 中 的 不 同 部 分 具 有 不 同 的 采 样 速 率 。 为 了 创 建 动 态 系 统 模 型 ，Simulink 提 供 了 一 个 建 立 模 型 方 块 图 的 图 形 用 户 接 口 (GUI) ， 这 个 创 建 过 程 只 需单 击 和 拖 动 鼠 标 操 作 就 能 完 成 ， 它 提 供 了 一 种 更 快 捷 、 直 接 明 了 的 方 式 ， 而 且 用 户可
11、以 立 即 看 到 系 统 的 仿 真 结 果 2。 Simulink 工 作 对 话 框 如图 2-1 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 4 页 共 25 页图 2-1 Simulink 工 作 对 话 框2.2 PSK 调制原理相移键控是利用载波的相位变化来传递数字信息的，而振幅和频率保持不变。在PSK 中，通常用初始相位 0 和 分别表示二进制“1” 和“0”。因此，PSK 信号的时域表达式为式（2-1） 。（2-()cosPSKetAt1）其中， 表示第 n 个符号大的绝对相位。当发送“0”时， 为 0; 当发送“1” 时， 为n
12、 nn。因此式（2-2）可改写为式（2-5） 。（2-cos()PSKcAtet2）其中 的概率为 P， 的概率为-P 。由于表示信号的两种码元的波形相cosAtcost同，极性相反，故 PSK 信号一般可以表述为一个双极性全占空矩形脉冲序列与一个正 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 5 页 共 25 页弦载波的相乘，即式（2-3） 。（2-()cosPSKettt3）其中 的表达式为式（2-4） 。()st（2-1()()nisistagtT4）这里， 是脉宽为 的单个矩形脉冲。发送二进制符号“0”时（ 取 1） ， 取 0()gtsT ia(
13、)PSKet相位；发送二进制符号“1” 时（ 取-1） ， 取 相位。这种以载波的不同相位直接ia()PSKet去表示二进制数字信号的调制方式，称为二进制绝对相移方式。PSK 信号模拟调制方法图如下图 2-2 所示：码型变换S(t)乘法器载波epsk(t)双极性不归零图 2-2 PSK 信号模拟调制方法图PSK 信号键控法调制图如下图 2-3 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 6 页 共 25 页载波180 度移相相开关电路0180度S(t)epsk(t)图 2-3 PSK 信号键控法调制图PSK 信号的解调通常采用相干解调法，如图 2-
14、4 所示：BPF 相乘器 LPF 抽样判决器载波 定时脉冲epsk(t) 输出图 2-4 PSK 信号相干解调法3 系统设计3.1 PSK 调制与解调在本设计中采用的是模拟调制法和相干解调法。整个 PSK 的仿真系统的调制与解调过程为：首先将信号源的输出信号经过码型变换后与载波通过相乘器进行相乘，送入加性高斯白噪声（AWGN）信道中传输。在接收端通过带通滤波器后再次与载波相乘，接着通过低通滤波器、抽样判决器，最后由示波器显示出各阶段波形，并用误码器观察误码率。PSK 调制与解调仿真电路图如图 3-1 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 7
15、页 共 25 页图 3-1 PSK 调制与解调仿真电路图载波模块参数设置如下图 3-2 所示：图 3-2 载波模块参数设置说明：将载波频率设置为 5Hz，由于在载波参数设置里，频率的单位是 rad/sec， 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 8 页 共 25 页所以这里 5Hz 即为 10*pi，对话框中“Frequency(rad/sec):”中应填“10*pi”。信号源参数设置如下图 3-3 所示：图 3-3 信号源的参数设置说明：对话框中“Sample time：”为抽样时间，其倒数就是信号源的频率，设信号源的数码率为 1bit/s，所以“
16、Sample time：” 中应填“1”。码变换器参数设置如下图 3-4 所示：图 3-4 码变换器参数设置 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 9 页 共 25 页说明：根据 PSK 调制解调原理，要将输入的单极性码转换为双极性码，且是二进制非归零码。对话框中“M-ary number：”即为进制数，所以应填 “2”。低通滤波器参数设置如图 3-5 所示：图 3-5 低通滤波器参数设置说明：低通滤波器的截止频率等于信号源频率。前面设置信号源频率为 1Hz，所以对话框中“Passband edge frequency (rads/sec)：”应填“
17、2*pi”。抽样判决器参数设置如下图 3-6 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 10 页 共 25 页图 3-6 抽样判决器参数设置说明：在 PSK 系统中，当发送不同符号的概率相等时，判决器的最佳判决门限为零，与接收机输入信号的幅度无关。所以对话框中“Quantization partition” 为“0”。量化器参数设置如图 3-7 所示：图 3-7 量化器参数设置说明：量化器抽样频率等于信号源频率。前面已经设置信号源频率为 1Hz，即抽样频率为 1Hz，所以对话框中“Sample time (-1 for inherited):”应
18、填“1”。误码器的参数设置如下图 3-8 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 11 页 共 25 页图 3-8 误码器的参数设置说明：根据示波器上的波形图可以看出，最后输出的波形入输入波形之间存在一个码元的延时，所以对话框中“Receive delay ：”应填“1”。PSK 调制与解调波形如图 3-9 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 12 页 共 25 页图 3-9 PSK 调制与解调波形由上图可知：第一路为信号源模块波形图，第二路为双极性码波形图，第三路为PSK 调制后波形图，第
19、四路为调制信号与载波相乘后波形图，第五路为经过低通滤波器后波形图，第六路为 PSK 解调波形图。由各波形可看出该 PSK 调制解调系统符合设计要求。3.2 加入高斯白噪声后的 PSK 调制与解调加入高斯白噪声后的 PSK 调制与解调仿真电路图如图 3-10 所示：图 3-10 PSK 调制与解调中加入高斯白噪声仿真图 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 13 页 共 25 页带通滤波器参数设置如图 3-11 所示：图 3-11 带通滤波器参数设置说明：带通滤波器的上截止频率和下截止频率分别等于载波频率与信号源频率的和与差。前面已设置信号源频率为 1
20、Hz，载波频率为 5Hz,计算得上、下截止频率分别为 6Hz、4Hz，转换成以 rads/sec 为单位即为 12*pi 、8*pi。所以“Lower passband edge frequency (rads/sec) Upper passband edge frequency (rads/sec)”应填“8*pi 、12*pi”。高斯白噪声的参数设置如下图 3-12 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 14 页 共 25 页图 3-12 高斯白噪声的参数设置加入高斯白噪声后误码器的参数设置如下图 3-13 所示：图 3-13 加入高斯白
21、噪声后误码器的参数设置说明：根据示波器上的波形图可以看出，最后输出的波形入输入波形之间存在两个码元的延时，所以对话框中“Receive delay ：”应填“2”。PSK 调制与解调中加入高斯白噪声后示波波形如图 3-14 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 15 页 共 25 页图 3-14 PSK 调制与解调中加入高斯白噪声仿真图由上图可知：第一路为信号源模块波形图，第二路为双极性码波形图，第三路为PSK 调制后的波形图，第四路为加入高斯白噪声后的波形图，第五路为经过带通滤波器后的波形图，第六路为经过带通滤波器后与载波相乘后的波形图，第
22、七路为经过低通滤波器后的波形图，第八路为 PSK 解调后的波形图。在 PSK 调制与解调中加入高斯白噪声后，波形出现了失真，解调也有误码存在，系统基本符合设计要求。3.3 误码率的计算误码率是衡量一个数字通信系统性能的重要指标。在信道高斯白噪声的干扰下，各种二进制数字调制系统的误码率取决于解调器输入信噪比，而误码率表达示的形式则取决于解调方式。理论上误码率计算 M 文件如图 3-15 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 16 页 共 25 页图 3-15 理论上误码率计算 M 文件计算误码率时，高斯白噪声的参数设置如图 3-16 所示：图
23、3-16 高斯白噪声参数设置图说明：将对话框中“Variance(vector or matriz)”的参数改为“var”,对应于绘制信噪比-误码率图程序中的参数设置。程序内容见附录。误码器参数设置如图 3-17 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 17 页 共 25 页图 3-17 误码器参数设置说明：将对话框中的“Output data”的参数改为“Workspace”。理论上误码率与信噪比关系曲线图如图 3-18 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 18 页 共 25 页图 3-1
24、8 理论上误码率与信噪比关系曲线图实际上误码率计算 M 文件如图 3-19 所示：图 3-19 实际上误码率计算 M 文件实际上误码率与信噪比关系曲线图如图 3-20 所示： 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 19 页 共 25 页图 3-20 实际上误码率与信噪比关系曲线图由上图可知，信噪比越大，误码率越小，与理论知识吻合。系统基本符合设计要求。 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 20 页 共 25 页4 出现的问题及解决方法4.1 出现的问题第一：示波器图形显示不全。第二：信号在经过所设计的带通
25、通信仿真系统后波形出现失真。第三：在没有加入高斯白噪声的情况下解调误码率不为 0。第四：加入噪声之后，误码率显著提高。第五：解调波形时无失真，但解码后波形严重失真。4.2 解决办法第一：在示波器图中修改 data history 中的 limit data points to last 的数据再重新运行 simulink 观察示波器即可看到准确图形。第二：这是因为信号进行抽样时没有满足奈奎斯特定理。为了能从取样信号中恢复原信号，抽样必须满足奈奎斯特定理，即抽样频率应大于或等于两倍的原始信号频率 3。第三：在出现误码数据时，我们可以通过修改 error rate calculation 中 re
26、ceive delay 的数据这种手段直到误码数据为零。在本设计中我们将 receive delay 设为 2 即可解决此种情况 4。第四：在本课程设计中通过提高载波的频率，并相应的改变带通滤波器的参数设置来解决此种情况的发生。第五：这是由于信号经过低通滤波器后会产生时延，而本次课程设计中信号是以帧的形式进行传输，因而在解调输出端若直接使用解调信号，将会产生波形失真。因此我们在解调输出端加入相对应的延时模块即可解决此种问题 5。 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 21 页 共 25 页结束语此次通信原理课程设计只有短短的两周时间，在这两周时间里，
27、我学会了用MATLAB/Simulink 仿真系统初步的设计方法。初步了解了如何用 MATLAB/Simulink这个仿真软件进行 PSK 系统的调制和解调的设计。并且明白了如何实现信号的频带传输。虽然在课程设计过程中遇到了不少的麻烦和困难，与此同时也深刻认识到自己存在的许多不足，理论知识的学习有待加强，实践的能力有所欠缺，遇到困难自己独立处理有待努力。在课程设计之前，我到图书馆借阅了相关书籍，但是理论的学习并没有让我体会到课程设计的真谛，光凭理论实在无法实现书中所要表达的内容，因此上机实践很关键。通过理论知识的初步学习，加上上机的实践，我渐渐地对 Simulink 有了一定的认识，但是在将书
28、上的理论方框图转换为 Simulink 模块的过程中出现了很大的困难，因为对软件不熟悉的缘故，许多模块都没有找到，因此也大大延误了设计的完成时间。但是在同学和老师的帮助指点下，许多问题得到了解决，同时也让我对此次设计产生了浓厚的兴趣，遇到问题也有了钻研的精神。这次课程设计不仅巩固我们在书上学习的基本内容，还在一定程度上提高了我们的动手能力。同时也让我懂得依靠团队的力量才能使我们更快更好的将各种工作和任务完成好。在这次课程设计中，我不仅巩固了过去学过的知识，同时还学到了许多书本上没有学到的知识培养了我对系统的分析能力。通过这次课程设计，我深刻的了解到，理论与实际的结合是十分重要的，单纯的理论知识
29、学习无法培养我们实际动手能力。只有把理论和实践结合起来，把理论作为实践的基础，把实践作为理论的延伸，使两者有机的结合在一起，才能提高自己的实际动手能力和独立思考的能力， 两周的时间很快就过去了，但是这段时间里的收获和感触却是我的一大笔财富。在此我衷心感谢同学们在课程设计中给予我的无私帮助，共同学习。也衷心感谢指导老师不辞辛苦指导我们课程设计的支持和辅导，感谢大家！ 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 22 页 共 25 页致 谢经过这次的课程设计，我了解到了做任何事都要有耐心、更是要细心做事。这次的课程设计让意识到自己的原理知识还是不够好，在今后的
30、学习中我们需要更努力的学习课本的专业知识，才能更好的服务于实践中。课程设计不仅是对前面所学知识的一种检验，而且也是对自己能力的一种提高。在设计中遇到了很多问题，最后在老师的辛勤的指导下，终于游逆而解，有点小小的成就感，终于觉得平时所学的知识有了实用的价值，达到了理论与实际相结合的目的，不仅学到了不少知识，而且锻炼了自己的能力，使自己对以后的路有了更加清楚的认识。在这次课程设计中也使我们的同学关系更进一步了，同学之间互相帮助，有什么不懂的大家在一起商量，听听不同的看法对我们更好的理解知识，所以在这里非常感谢帮助我的同学。最后,在这里，我要感谢老师悉心指导和无私地帮助，感谢我们的老师，同时也感谢同
31、学们对我的耐心指导，在今后的学习中我也一定倍加努力。 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 23 页 共 25 页参考文献1 孙屹,戴妍峰 . SystemView 通信仿真开发手册 . 北京:国防工业出版社,20042 桑林,郝建军 ,刘丹谱.数字通信. 北京:北京邮电大学出版社,20023 樊昌信,曹丽娜 . 通信原理(第版). 北京:国防工业出版社,20084 康华光，邹寿彬 .电子技术基础（数字部分) .北京：高等教育出版社，20045 樊昌信,曹丽娜 .通信原理.北京：国防工业出版社,2007 基于 MATLAB/Simulink 的 PS
32、K 频带传输系统仿真与性能分析 第 24 页 共 25 页附录：误码率计算（理论值）% 程序名称 ：zhenque.m% 程序功能：绘制理论信噪比-误码率曲线。% 程序作者 ：贺健% 最后修改日期：2011-1-6x=-6:18; % 定义信噪比范围for i=1:length(x); snr=x(i); % 每次运行信噪比增加 1dBvar=0.5/(10.(0.1*snr); % 计算对应的噪声功率sim(zhenque.mdl); % 运行模型文件ebr(i)=ErrorVec(1); % 保存当次运行误码率iendebr_th=erfc(sqrt(-1/3).10);semilogy(
33、x,ebr,x,ebr_th) % 绘制信噪比 -误码率曲线，纵坐标为对数坐标xlabel(信噪比 r/dB);ylabel(误码率 Pe);title(误码率与信噪比关系曲线); 基于 MATLAB/Simulink 的 PSK 频带传输系统仿真与性能分析 第 25 页 共 25 页附录：误码率计算（实际值）% 程序名称 ：hj.m% 程序功能：绘制实际信噪比-误码率曲线。% 程序作者 ：贺健% 最后修改日期：2011-1-6x=-6:18; % 定义信噪比范围for i=1:length(x); snr=x(i); % 每次运行信噪比增加 1dBvar=0.5/(10.(0.1*snr); % 计算对应的噪声功率sim(zhenque.mdl); % 运行模型文件ebr(i)=ErrorVec(1); % 保存当次运行误码率endisemilogy(x,ebr) % 绘制信噪比-误码率曲线，纵坐标为对数坐标xlabel(信噪比 r/dB);ylabel(误码率 Pe);title(误码率与信噪比关系曲线);grid on;
