基于MATLAB的BPSK调制仿真及性能分析概要

1、* 实践教学 * 兰州理工高校运算机与通信学院2022 年春季学期通信系统仿真训练 课程设计题 目： 基于 MATLAB 的 BPSK 调制仿真及性能分析专业班级： _ 姓 名： _ 学 号： _ 指导老师： _ 成 绩摘要 BPSK 是在 21 世纪的现阶段,数字通信系统是现代通信系统的主流,在社会生活各个 方面占据重要位置,调制与解调是一个不行或缺的环节,它大大提高了信号的传输性能和安 全性； BPSK 作为数字通信系统中的一种简洁、基础的调制解调方法,抗干扰才能强,简洁 仿真实现；在明白并把握了 BPSK 调制与解调原理后,学习MATLAB 仿真软件的使用方法；本次课程设计为基于 MAT

2、LAB 的 BPSK 调制仿真；本次课设着重介绍了算法的实现,并 采纳MATLAB 程序仿真测试了 BPSK 过程中单极性不归零编码、脉冲成形、PSK 调制、信号通过 AWGN 信道、载波复原、解调、解码等过程；关键词： BPSK 、调制解调、 MATLAB 仿真、数字通信系统1目录前言 . . 3 一、BPSK 原理 . 4 1.1 MPSK 的介绍 . 4 1.1.1 MPSK 的介绍 . 4 1.2 BPSK 简介 . 5 1.2.1 BPSK 的概念及分类 . 5 1.3 BPSK 信号的产生 . . 5 1.4 BPSK 调制 . . 5 1.4.1 调制的概念 . . 5 1.4.

3、2 调制的种类 . . 6 1.4.3 调制的作用 . . 6 1.4.4 调制方式 . . 6 1.4.5 BPSK 的调制原理 . 7 1.5 BPSK 的解调 . . 8 1.6 高斯噪声 . 8 1.6.1 高斯白噪声 . . 8 二、设计思路 . . 9 三、设计仿真及分析 . . 10 1 设计仿真 . 10 2 结果总体分析 . 13 总结 . . 14 参考文献 . . 15 致谢 . . 162、八 , 、刖言数字通信系统使用的最基本的就是PSK , 实际上由于交调的存在,大都使用PSK , 由于 电波主要在自由空间传播,信道参数比较稳固,信道的主要干扰是高斯白噪声,因而可视

4、为 恒参信道；理论分析指出,在恒参信道中,采纳 PSK 的误码率最低,更有效地利用卫星频带,在接收机噪声的作用下,PSK的误码率最低；因而在实际中得到广泛的应用；特殊是高数数据传输中 FSK 是一种很重演的调制方式；因此本课程设计尝试对 PSK 中常用的 BPSK 的调制技术进行通信系统仿真训练的设计实 践；本次课程设计为基于 MATLAB 的 BPSK 调制仿真及性能分析；本次课设着重介绍了算法 的实现,并采纳MATLAB 程序仿真测试了 BPSK 过程中单极性不归零编码、脉冲成形、PSK调制、信号通过 AWGN 信道、载波复原、解调、解码等过程；本设计要求采纳 Matlab 实现对 BPS

5、K 通信系统的蒙特卡罗仿真并且绘制相关的图形,此 题比较难做涉及到的方面很多,对 BPSK 及蒙特卡罗仿算法很不明白查阅了很多资料然而涉 及的有较少,但查资料当中学到了很多不知道的东西,加深了对 PSK 的明白；以二进制调相为例,取码元为 1时,调制后载波与未调载波同相；取码元为 0时,调 制后载波与未调载波反相； 1和0时调制后载波相位差180 ；就模拟调制法而言, 与产 生 2ASK 信号的方法比较,只是对 st要求不同,因此 BPSK 言号可以看作是双极性基带信号 作用下的 DSB 调幅信号；而就键控法来说,用数字基带信号 st掌握开关电路,挑选不同相 位的载波输出,这时 st为单极性

6、NRZ 或双极性NRZ 脉冲序列信号均可；本设计分析 BPSK 调制和解调系统,并用 BPSK 软件仿真调制系统,最终建立蒙特卡洛 仿真模型,统计系统误码率；该题目概括了MATLAB 技术、通信系统原理等课程的 主要学问点,通过该设计能够培育和提高同学综合设计才能,为今后的学习和工作积存体会；3、BPSK 原理1 MPSK 的介绍1.1.1 MPSK 的介绍MPSK 即多进制相移键控,又称为多相制；这种键控方式是多进制键控的主要方式；在 M 进制的相移键控信号,用 M 个相位不同的载波分别代表 M 个不同的符号；假如载波 有 2n 个相位,它可以代表 n 位二进制码元的的不同组合的码组；多进制

7、相移键控也分为多 进制肯定相移键控和多进制相对相移键控；示为在 MPSK 信号中,载波相位 有 M 种可能取值, n=2n n/Mn=1,2, -M；因此 MPSK 信号可表St = COS 00t+ n = cos ot+2 n/M 如载波频率是基带信号速率的整数倍,就上式可改写为 St= 刀 gtT scos ot+ 0=cos d 工 cos gt-nT s -sin d 工 sin gt-nT 式中 gt是高度为 1、宽度为 Ts 的门函数； Ts 为 M 进制码元的连续时间；亦即 kk=log 2M 比特二进制 码元的连续时间； n 为第 N 个码元对应的相位,共有 M 个不同的取值

8、；上式说明,MPSK 信号可等效为两个正交载波的MASK 信号之和；所以,MPSK 信号的 带宽和 MASK信号的带宽相同；因此,MPSK 系统是一种高效率的信息传输方式；但是,当 M 的取值增加时,载波间的相位差也随之削减,这就使它的抗噪声性能变差；数字通信系统模型如图 2.1 所示；图 1.1.1 数字通信系统模型2 BPSK 简介1.2.1 BPSK 的概念及分类BPSK Binary Phase Shift Keying , 把模拟信号转换成数据值的转换方式之一,禾 U 用偏 离相位的复数波浪组合来表现信息键控移相方式；BPSK 使用了基准的正弦波和相位反转的 波浪, 使一方为0, 另

9、一方为 1, 从而可以同时传送接受 2 值1 比特 的信息；移相键控分为肯定移相和相对移相两种；以未调载波的相位作为基准的相位调制叫作绝 对移相；以二进制调相为例,取码元为“ T时,调制后载波与未调载波同相；取码元为“ 0” 时,调制后载波与未调载波反相；“ T 和“ 0” 时调制后载波相位差 180 ；3 BPSK 信号的产生就模拟调制法而言,与产生 2ASK 信号的方法比较,只是对 st要求不同,因此 BPSK 信号可以看作是双极性基带信号作用下的 DSB 调幅信号；而就键控法来说,用数字基带信 号 st掌握开关电路,挑选不同相位的载波输出,这时 冲序列信号均可；st 为单极性 NRZ 或

10、双极性 NRZ 脉BPSK 信号属于 DSB 信号,它的解调,不再能采纳包络检测的方法,只能进行相干解 调；BPSK 言号相干解调的过程实际上是输入已调信号与本地载波信号进行极性比较的过程,故常称为极性比较法解调；由于 BPSK 言号实际上是以一个固定初相的未调载波为参考的,因此,解调时必需有与此同频同相的同步载波；假如同步载波的相位发生变化,如0 相位变为 n 相位或 n 相位变为 0 相位,就复原的数字信息就会发生“ 0” 变“1” 或“ 1” 变“ 0” ,从而造成错误的恢 复；这种由于本地参考载波倒相,而在接收端发生错误复原的现象称为“ 倒 n 现象或“ 反向工作” 现象；肯定移相的主

11、要缺点是简洁产生相位模糊,造成反向工作；这也是它实际 应用较少的主要缘由；4 BPSK 调制1.4.1 调制的概念调制 modulation 就是对信号源的信息进行处理加到载波上,使其变为适合于信道传 输的形式的过程,就是使载波随信号而转变的技术；一般来说,信号源的信息 的频率重量,称为基带信号；基带信号往往不能作为传输信 也称为信 源 含有直流重量和频率较低 号,因此必需把基带信号转变为一个相对基带频率而言频率特别高的信号以适合于信道传 输；这个信号叫做已调信号,而基带信号叫做调制信号；调制是通过转变高频载波即消息的载体信号的幅度、相位或者频率,使其随着基带信号幅度的变化而变化来实现的； 而

12、解 调就是将基带信号从载波中提取出来以便预定的接收者 也称为信宿 处理和懂得的过5程；调制的种类调制的种类很多,分类方法也不一样；按调制信号的形式可分为模拟调制和数字调制；用模拟信号调制称为模拟调制；用数据或数字信号调制称为数字调制；按被调信号的种类 可分为脉冲调制、正弦波调制和强度调制 如对非相干光调制 等；调制的载波分别是脉冲,正弦波和光波等；正弦波调制有幅度调制、频率调制和相位调制三种基本方式,后两者合称为角度调制；此外仍有一些变异的调制,如单边带调幅、残留边带调幅等；脉冲调制也 可以按类似的方法分类；此外仍有复合调制和多重调制等；不同的调制方式有不同的特点 和性能；调制的作用调制在通信

13、系统中有特别重要的作用；通过调制,不仅可以进行频谱搬移,把调制信 号的频谱搬移到所期望的位置上,从而将调制信号转换成适合于传播的已调信号,而 且它对系统的传输有效性和传输的牢靠性有着很大的影响,调制方式往往打算了一个通信 系统的性能；调制方式在通信中,我们经常采纳的调制方式有以下几种：1 、模拟调制：用连续变化的信号去调制一个高频正弦波,主要有：幅度调制 调幅 AM、双边带 DSB、单边带 SSB 、残留边带 VSB 以及独立边带 ISB; 角度调制 调频 FM,调相 PM 两种；由于相位的变化率就是频率,所以调相波 和调频波是亲密相关的；2、 数字调制：用数字基带信号掌握载波,把数字基带信号

14、变换为数字带通信号 已调信号 的过称成为数字调制,主要有：振幅键控 ASK ；频率键控 FSK ；相位键控 PSK ；3、 脉冲调制：用脉冲序列作为载波,主要有：脉冲幅度调制 PA M ： Pulse Amplitude Modulation ；脉宽调制 PDM ： Pulse Duration Modulation ；脉位调制 PPM ：Pulse Position Modulation ；6145 BPSK 的调制原理二进制移相键控是用二进制数字信号 0 和 1 去掌握载波的两个相位 0 和 n 的方法；在2PSK 中,通常用初始相位 0 和 n 分别表示二进制 1 和 0；因此, 2PS

15、K 信号的时域表达式为：e2PSK t = Acos ct n 式中,;:n 表示第 n 个符号的肯定相位 : ；0,发送“* 发送 “ T时0” 时e2PSK t= A co 的ct, 概率为 PA cos ct, 概率为 1 - P 因此,上式可以改写为由于两种码元的波形相同,极性相反,故BPSK 信号可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘：e 2PSKt二St cos ct式中st 八 angt - nTs n这里 st为双极性全占空 非归零 矩形脉冲序列, gt是脉宽为 Ts 的单个矩形脉冲,而 an 的统计特性：1, 概率为 P a =n -1,概率为 1-P 7

16、图 1.4.5 BPSK 信号的波形示例5 BPSK 的解调PSK 信号本身就是利用相位传递信息的,所以在接收BPSK 信号的解调方法是相干解调法；由于端必需利用信号的相位信息采纳相干解调法来解调信号；给出了一种 2PSK 信号相干接收设备的原理框图；图中经过带通滤波的信号在相乘器中 与本地载波 相乘,然后用低通滤波器滤除高频重量,再进行抽样判决,判决器是按极性进 行判决,得到最终的二进 制信息；详细波形如下图：R 1AK II 1E 1 - 1- Q- 1图 1.5 BPSK 各点时间波形 1.6 高斯噪声1.6.1 高斯白噪声高斯白噪声：假如一个噪声,它的幅度分布听从高斯分布,而它的功率谱

17、密度又是均 匀分布的,就 称它为高斯白噪声；热噪声和散粒噪声是高斯白噪声；所谓高斯白噪声中的 高斯是指概率分布是正态函数,而白噪声是指它的二阶矩不相关,一阶矩为常数,是指先 信号的两个不同方面的问题；后信号在时间上的相关性；这是考查一个、设计思路2.1 系统总体设计2.1.1 总体的系统设计方案如图 2.1.1 所示 : 8图 2.1.1 系统方案图BPSK（B in ary Phase Shift Keyi ng ）是二进制相移键控,它是一种相位调制算法；相位调制（调相）是频率调制（调频）的一种演化,载波的相位被调整用于把数字信息的比特编码到 每一词相位转变（相移）；数字带通传输中一般利用数

18、字信号的离散取值特点通过开关键控载波,从而实现数字调 制,比如对载波的振幅、频率和相位进行键控分别可以获得振幅键控 和相移键控（ PSK ）；（ASK ）、频移键控（ FSK）就总体而言,信道应看作一个线性系统,满意线性叠加原理；信号在信道中传输,存在 衰耗和时延,信道中总是存在噪声,信号在实际信道中传输,将会产生失真,任何信道都有 肯定的频率带宽,信道不可能传送功率无限大的信号；9三、设计仿真及分析3.1 设计仿真3.1.1 仿真框图 在发送端,通过随机函数立即产生十比特二进制比特序列；然后把这十比特序列在频 率 fc=4000HZ 的载波上进行传输,并且采样频率 fs=8000HZ ；经过

19、调制后,调制信号就可以 在信道上 传输；但是在实际的信道中传输时,会叠加很多噪声,因此,程序模拟在实际信 道上传输,产生噪声,叠加到已调信号上；在接收端,通过相干解调的方法,把接收到的叠加有噪声的信号进行解调,但是解调 后的信号仍不 是最先发送的二进制比特流,需要对解调得到的信号进行抽样判决,才能得 到发送的二进制比特流,即发送信号；软件的仿真流程图如图 3.1.1 所示；图 3.1.1 基于 MATLA 的 BPSK 调制解调仿真框图3.1.2 仿真结果及分析 基带信号、 BPSK 调制信号以及载波信号幅值 高斯白噪声10何20D3星芾信寻顎谊l.:E养 皿%0 5050 100 Itm生

20、吨 养 怛 嵋 建1000 栽液信号頻谱T 500 uldm-%调制信号频谙100lllns芥 咀QFE500 100 展 盖 缶X LU-50 0 1000 圻 iz图 3.1.2 基带、载波、调制信号幅值图11加入加性高斯白噪声后的BPSK 言号如下图所示 : BPSK 言号图 3.1.3 加入加性高斯白噪声后的由上图可知信号在加入高斯白噪声时,信号质量明显变差,会造成信号判决时误码率 急剧增加,而信道的噪声功率固定不变,因此在传输信号时,应增加信号功率,来提高信 息的传输效率；解调以后出来的波形如图 3.1.4 所示：-21 - 0 2DG 400 600 800 10DD 1200 1

21、400 1SOO 1800 2022加噪后 解调信号2 - ： - ,：- - - - . - - .- - ” - ; . . I - ;. . :. . : - ; - . . . . -1 . :- . : . . . : . : . . : - 一 : : . . I I ii I 丨 I I i I0 200 400 E00 800 1000 1200 1400 1600 1000 2022加噪后去掉推迟的解调信号炸）图 3.1.4 加躁后解调信号12经过低通滤波器信号波形 : 图 3.1.5 与复原载波相乘后的信号蒙特卡洛分析实际误码率及理论误码率仿真图如下图所示：12346670

22、9 10图 3.1.6 蒙特卡洛分析误码率图码元为 10000 个的时候,蒙特卡洛仿真统计图如图5.6,从上图可以看出实际误码率要比理论误13码率大得多,是由于理论误码率考虑到的影响因素仅仅是主要因素,而各种外界干扰及温度等都没有考虑；当信噪比为 3.2 结果总体分析6db 的时候仿真结果和实际结果有较大误差；BPSK 调制解调系统的抗噪声才能较强；但是,在操作中发觉 Matlab 仿真系统运行时响 应时间很长,需要运算机运行很长时间才能完成仿真运算过程；经过小组成员的分析,我 们认为,除了算法导致了系统运算量大以外,BPSK 系统本身的信号传输效率也不高,这也印证了通信原理教材中提到的 BP

23、SK S统的特点,即抗噪声干扰才能强,但数据传输效率低 下；信道噪 声对系统性能的主要影响是在接收信号中引入了比特差错；在二进制系统中 , 比 特差错率表现为将符号 1 误认为 0,或将符号 0 误认为符号 1；很明显比特差错的频率越高 , 接收机的输出信号与原始信息之 间的差异就越大；在存在信道噪声的情形下 , 可以用平均符 号差错概率来衡量二进制信息传输的逼真 度； 平均符号差错概率的定义为 , 接收机输出的重 构符号与所传输的二进制不相同的平均概率；在原 始二进制波形中的全部比特均具有相同 重要性的条件下 ,平均符号差错概率又称为误比特率（BER ）；但是 ,在重构原始消息信号的模 拟波

24、形时 , 不同的符号差错可能需要区分对待；例如码字（表示消息 信号的量化抽样值） 中 重要的比特发生的错误要比不重要的比特发生的错误有害得多；14总结通过本次设计,让我能够将课本的学问与实践相结合起来,更全面的明白书本的学问,此外,课程设计仅靠课本学问是远远不够的,通过仿真实践补充了大量的课外学问；同时 对 BPSK 调制解调的工作原理有了更加深刻的懂得,在设计之前,通过图书馆、上网浏览并 收集了很多的相关材料,但当进行设计时,却也遇到了很多的问题,让我体会到了设计的要求在于系统性、可行性、精确性,很多问题的显现给我们的设计带来了难度,也同时是 一次庞大的挑战,最终,在老师和同学的帮忙下,克服

25、了种种困难,顺当的完成了本次通信系统仿真训练的课程设计,在此,第一要感谢学校支配此次课程设计,让我有机会对本 次课程设计能够深化明白和设计,再次感谢老师细心地指导,仍有同组同 学的相互沟通和 帮忙,使我完成了本次通信系统仿真课程设计,让我学到了很多课本里学不到的学问；在 课程设计过程中查阅了大量的有关 BPSK 调制解调设计的书籍,巩固了以前所学过的学问,而且 过该门课程的学 学到了很多在书本里没学过的学问；通信原理是通信专业的一门主干技术基础课,通 习,使我们把握确定信号调制解调原理的特性,信号的传输的特性,确定 信号的发送与接收的基本方 法以及某些典型系统引出的一些重要的基本概念；本设计要

26、求采纳 Matlab 实现对 BPSK 1 信系统的蒙特卡罗仿真并且绘制相关的图形,此 题比较难 BPSK 及蒙特卡罗仿算法很不明白查阅了很多资料然而涉 及的有较少,但查 做涉及到的方面很多,对 资料当中学到了很多不知道的东西,加深了对 PSK 的明白；完成这次课程设计让我收成颇多,不但进一步把握了通信原理这一核心课程的学问及 一门专业仿 真软件的基本操作,而且仍提高了自己的设计才能及分析问题才能；通过这次 课程设计使我懂得了理论与实际相结合的必要性,只有理论学问是远远不够的,只有把所 从理论中得出结论,才能真正做到学以致用；15学的理论学问与实践相结合起来,参考文献曹志刚等编,现代通信原理M

27、 ,清华高校出版社,2022程佩青 . 数字信号处理教程 M. 清华高校出版社刘卫国 .MATLAB 程序设计教程 M. 中国水利水电出版社桑林,郝建军,刘丹,数字通信 M. 北京邮电高校出版社,2022苗云长等主编,现代通信原理及应用 M；电子工业出版社,2022吴伟铃,庞沁华,通信原理 M,北京邮电高校出版社,2022张圣勤, MATLAB7.C 有用教程 M,机器工业出版社, 2022 邵玉斌, Matlab/Simulink通信原理建模与仿真实例分析M,清华高校出版社 , 2022 樊昌信,通信原理（第 6 版） M. 北京 :国防工业出版社.2022 16致谢通过本次运算机通信课程设

28、计,让我在除了课本的学问之外的学问有了更好的懂得,对 8PSK 调制解调的工作原理有了更好的懂得,在设计之前,收集了很多的材料,但当真正 深化设计时,却也遇到了诸多的问题,让我体会到了设计的要求在于系统性,可行性,准确性,诸多问题的显现给我们的设计带来了难度,也同时是更大的一次挑战,最终,在老 师以及同学的帮忙下,克服了种种困难,顺利的完成了本次通信仿真训练的课程设计,在此,第一要感谢学校支配此次课程设计,让我有机会对本次课程设计能够深化懂得和设计,再次感谢老师的细心指导和改正,仍有同组同学的相互团结和帮助,使我完成了本次运算 机通信课程设计,使我在求学的道路上有了更多方面学问的获得；设计过程

29、中查阅了大量 的有关 PSK 调制解调设计的书籍,巩固了以前所学过的学问,而且学到了很多在书本上所 没有学到过的学问；通过这次课程设计使我懂得了理论与实际相结合的必要性,只有理论 学问是远远不够的,只有把所学的理论学问与实践相结合起来,从理论中得出结论,才能 真正为社会服务,从而提高自己的实际动手才能和独立摸索的才能；再次感谢学校,感谢17老师,感谢同学；附录clear all close all clc num=10; % 码元个数tnum=200;% 码元长度N=num*tnum ;%10 个码元整体长度a=randint1,num,2; % 产生 1 行 num 列的矩阵,矩阵内 t=0

30、:0.05:9.99;%t 从 0 到 9.99 ,间隔为 0.05 s=;c=; for i=1:num %i 从 1 到 10 循环ifai=0 0 和 1 随机显现 fc=0.25; % 载波频率为 0.5 A=zeros1,tnum; %i=0 时,产生一个码元长度为tnum 200 的 0 码元else end A=ones1,tnum; %i=1 时,产生一个码元长度为 tnum 200 的 1 码元s=s A; %s 为随机基带信号 cs=sin2*pi*fc*t; c=c cs; %c 为载波信号 end %采纳模拟调制方法得到调制信号 s_NRZ=; for i=1:num

31、%i 从 1 到 num 10 循环 ifai=0 else end A=ones1,tnum; %i=0 时,产生一个码元长度为 tnum 200 的 1 码元A=-1*ones1,tnum; %i 非 0 时,产生一个码元长度为 tnum 200 的 -1 码元18s_NRZ=s_NRZ,A; %s_NRZ 为双极性非归零码end e=s_NRZ.*c; %e 为 BPSK 调制信号figure1; %图 1 subplot3,2,1; % 图 1 分为 3X 2 部分的第一部分plots; %作 s 基带信号 的波形图grid on; axis0 N -2 2; % 横轴长度为 0 到

32、N,纵轴范畴为 -2 到+2 xlabel 基带信号 st; %x 轴的注释ylabel 基带信号幅值 ; %y 轴的注释 subplot323;plotc;grid on; axis0 N -2 2; xlabelBPSK 载波信号 ; ylabelBPSK 载波信号幅值 ; subplot325;plote;grid on; axis0 N -2 2; xlabel BPSK 调制信号 ; ylabel BPSK 调制信号幅值 ; %信号的频谱Fs=200; % 采样频率%作 c BPSK 载波信号 的波形图%作 e BPSK 调制信号 的波形图n=lengths; % 基带信号长度f=0

33、:Fs/n:Fs-Fs/n-Fs/2; %修正频率f 的范畴S=ffts; % 基带信号 s 的快速傅里叶变换E=ffte; %基带信号e 的快速傅里叶变换C=fftc; %基带信号c 的快速傅里叶变换subplot322; plotf,absfftshiftS; %基带信号的频谱title 基带信号频谱 ; xlabelf/hz;ylabelSw; grid on; subplot324; plotf,absfftshiftC; % 载波信号的频谱19title 载波信号频谱 ; xlabelf/hz;ylabelCw; grid on; subplot326; plotf,absfftsh

34、iftE; %调制信号的频谱title 调制信号频谱 ; xlabelf/hz;ylabelEw; grid on; %加高斯噪声am=0.7; % 输入信号经信道后振幅由 1 衰减为 0.7SNR=5; % 输入信噪比sn r=10ASNR/10; N0=am*am/2/snr; % 运算噪声功率N0_db=10*log10N0; % 将噪声功率转换为 dBW ni=wgn1,N,N0_db; % 产生 1 行 N 列的高斯噪声yi=e+ni; %BSK 已调信号中加入白噪声,输入信噪比为 SNRfigure2;subplot2,1,1;plotyi;grid on; xlabel 加入高斯

35、白噪声的已调信号 yit; %带通滤波器b1,a1 = BUTTER3,2*pi*0.0001,2*pi*0.01; %运算带通滤波器的 Hz 系数y=filterb1,a1,yi; %对信号 yi 进行滤波,得到信号 y figure2;subplot2,1,2;ploty;grid on; xlabel 经带通滤波器后信号 ; %与复原载波相乘x1=2*c.*y; figure3;subplot2,1,1;plotx1;grid on xlabel 与复原载波相乘后的信号 x1t; %低通滤波器b2,a2=butter2,0.005; % 运算 Hz系数,频率为 1/200 x=filterb2,a2,x1; % 对信号 x1 滤波,得到信号x figure3;subplot2,1,2;plotx;grid on axis0 N -2 2;xlabel 经低通滤波器后信号波形 % 抽样判决20 x=fun_panjuex;% 调用函数,进行抽样判决 xlabel 加噪后解调信号 xt ; axis0 N -2 2; %排除推迟x=fun_yancx; % 调用函数