BPSK的matlab仿真简洁代码-基带频带
MATLAB仿真应用-第5章(2)
键控）的仿真试验框图，图5-41是基带FSK（移频键控）调制
信号的频谱，图5-42是MSK（最小移频键控）的仿真试验框
图，图5-43是基带MSK（最小移频键控）调制信号的频谱，
图5-44是GMSK（高斯滚降最小移频键控）的仿真试验框图，
图5-45是基带GMSK（高斯滚降最小移频键控）调制信号的频
谱。表5-31~表5-37分别给出了三个仿真系统中模块的主要参
15
图5-44 GMSK基带调制仿真系统
2021/2/4
1
16
图5-45 GMSK基带调制信号频谱
2021/2/4
1
17
表5-31BernoulliRandomBinaryGenerator （伯努利二进制随机数产生器）的主要参数
2021/2/4
1
18
表5-32 SpectrumScope（频谱仪）的主要参数
2021/2/4
1
19
表5-33M FSKModulatorBaseband （基带M FSK
2021/2/4
1
20
表5-34MSKModulatorBaseband （基带MSK
2021/2/4
1
21
表5-35 GMSKModulatorBaseband （基带GMSK调制器）的主要参数
2021/2/4
2021/2/4
1
41
表5-45ErrorRateCalculation （误码率计算）的主要参数
2021/2/4
1
42
仿真结果说明在相同传输条件下，QPSK、OQPSK 调制以比BPSK调制高1倍的速率传输信息，QPSK误码 率高于BPSK约30倍，OQPSK误码率高于BPSK约1.4倍。 由于是相位调制，频谱反映了具有不同相位特性的同 一载波，因此频谱特性一致。正因为是相位调制，每 符号采样取16（大于1）可以得到更好的系统性能，此 时在误码表中应该正确地选取与之相应的接收延迟参 数。星座图显示的是传输环境SNR（信噪比）为6dB时 的图形。
MATLABPSK调制与解调及仿真
学士学位毕业设计（论文）基于MATLAB的PSK调制和解调及仿真摘要Psk调制是通信系统中最为重要的环节之一，Psk调制技术的改进也是通信系统性能提高的重要途径。
本文首先分析了数字调制系统的基本调制解调方法，然后，运用Matlab及附带的图形仿真工具——Simulink设计了这几种数字调制方法的仿真模型。
通过仿真，观察了调制解调过程中各环节时域和频域的波形，并结合这几种调制方法的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。
最后，在仿真的基础上分析比较了各种调制方法的性能，并通过比较仿真模型与理论计算的性能，证明了仿真模型的可行性。
另外，本文还利用Matlab的图形用户界面（GUI）功能为仿真系统设计了一个便于操作的人机交互界面，使仿真系统更加完整，操作更加方便。
关键词：数字调制；分析与仿真；Matlab；Simulink；GUI图形界面ABSTRACTIn this paper, methods of psk modulation are introduced firstly. Then their simulation models are bu ilt 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。
--BPSK通信系统的计算机性能分析与MATLAB仿真.
淮海工学院课程设计报告书课程名称：通信系统的计算机仿真设计题目:BPSK通信系统性能分析与MATLAB仿真系（院）：学期：专业班级：姓名：学号：BPSK通信系统性能分析与MATLAB仿真1绪论随着通信技术的发展，信号处理方面硬件设计与专业软件设计结合日趋紧密，已经有一些公司开付出专业数字信号处理软件。
比较优秀的而且得到广大技术人员认可的有MATLAB。
MATLAB等优秀软件使仿真技术得到很好的应用。
通过对通信过程的仿真，我们就可以在低成本的条件下检测某一个方案是否能够实现，是否有更好的方案可以代替原来的方案，这样对通信的研究就站在了一个更高的起点，使通信技术的发展日新月异，近几年手机的普及率的迅速提高就从侧面反映移动通信技术的发展。
现代移动通信系统的发展是以多种先进的通信技术为基础发展起来的。
移动通信的主要基本技术包括调制技术、移动信道中颠簸的传播特性、多址方式、抗干扰技术以及组网技术。
在移动通信中，数字调制解调技术是关键技术，其中数字调相信号具有数字通信的诸多优点，在数字移动通信中广泛使用它来传送各种控制信息。
1.1 研究背景与研究意义随着通信系统复杂性不断增加,传统设计已不能适应发展的需要,通信系统的模拟仿真技术越来越受到重视,因此在设计新系统时,要对原有的系统做出修改或者进行相关研究,通常要进行建模和仿真,通过仿真结果来衡量方案的可行性,从中选择合理的系统配置和参数设置,然后进行实际应用。
MATLAB 作为一种功能强大的数据分析和工程计算高级语言,已被广泛应用于现代科学技术研究和工程设计的各个领域。
调制解调技术在通信系统中不可或缺，因此，基于MATLAB的调制解调模块仿真设计对通信系统的教学和科研都具有积极的意义。
1.2 课程设计的目的和任务本次课程设计是根据“通信工程专业培养计划”要求而制定的。
通信系统的计算机仿真设计课程设计是通信工程专业的学生在学完通信工程专业基础课、通信工程专业主干课及科学计算机仿真专业课后进行的综合性课程设计。
完整版PSK原理及MATLAB仿真
S(t)
相对
码元1 0 1 10 0
调相
t
图4 相对移相的波形示意图
一般情况下，相对移相可通过对信码进行变换 和绝对移相来实现。将信码经过差分编码变换成新 的码组 ——相对码，再利用相对码对载波进行绝对 移相，使输出的已调载波相位满足相对移相的相位 关系。
3. 2PSK的调制原理
? 数字调相：如果两个频率相同的载波同时开始 振荡，这两个频率同时达到正最大值，同时达 到零值，同时达到负最大值，它们应处于“同 相”状态；如果其中一个开始得迟了一点，就 可能不相同了。如果一个达到正最大值时，另 一个达到负最大值，则称为“反相”。
PSK 原理及MATLAB 仿真
前言
? 移动通信的迅速发展，离不开很多关键技术的支
持与应用 ，数字调制在通信领域中就发挥着重大
的作用
? 为了使数字信号在带通信道中传输，必须使用数 字基带信号对载波进行调制，以使信号与信道的 特性匹配
? 由于PSK在生活中有着广泛的应用，本课题主要 介绍了 PSK 波形的产生和仿真过程
形');
所得波形图如下：
for i=1:7 if a(i)==0 for j=1:50 out(j+50*(i-1))=f(j); end else for j=1:50 out(j+50*(i-1))=-f(j); end end
end for i=1:7
for j=1:50 m=0;n=0; if out(j+50*(i-1))-f(j)==0 m=m+1; end if out(j+50*(i-1))+f(j)==0 n=n+1; end
码元定时与载波的关系可以是同步的，以 便清晰观察码元变化时对应调制载波的相应 变化；也可以是异步的，因为实际的系统都 是异步的，码元速率约为1Mbt/s。
基于matlab的QPSK与BPSK信号性能比较仿真
┊┊┊┊┊┊┊┊┊┊┊装┊┊┊┊┊订┊┊┊┊┊线┊┊┊┊┊┊┊┊┊┊┊┊┊目录第一章概述 (1)第二章QPSK通信系统原理与仿真 (1)2.1 QPSK系统框图介绍 (1)2.2QPSK信号的调制原理 (2)2.2.1QPSK信号产生方法 (2)2.2.2QPSK星座图 (2)2.3QPSK解调原理及误码率分析 (3)2.3.1QPSK解调方法 (3)2.3.2QPSK系统误码率 (3)2.4QPSK信号在AWGN信道下仿真 (4)第三章BPSK通信系统原理与仿真 (4)3.1BPSK信号的调制原理 (4)3.2BPSK解调原理及误码率分析 (4)第四章QPSK与BPSK性能比较 (5)4.1QPSK与BPSK在多信道下比较仿真 (5)4.1.1纵向比较分析 (5)4.1.2横向比较分析 (7)4.2仿真结果分析 (7)4.2.1误码率分析 (7)4.2.2频带利用率比较 (7)附录 (8)代码1 (8)代码2 (8)代码3 (10)代码4 (12)┊┊┊┊┊┊┊┊┊┊┊装┊┊┊┊┊订┊┊┊┊┊线┊┊┊┊┊┊┊┊┊┊┊┊┊第一章概述QPSK是英文Quadrature Phase Shift Keying的缩略语简称，意为正交相移键控，是一种数字调制方式。
它以其抗干扰性能强、误码性能好、频谱利用率高等优点，广泛应用于数字微波通信系统、数字卫星通信系统、宽带接人、移动通信及有线电视系统之中。
BPSK是英文Binary Phase Shift Keying的缩略语简称，意为二相相移键控，是利用偏离相位的复数波浪组合来表现信息键控移相方式的一种。
它使用了基准的正弦波和相位反转的波浪，使一方为0，另一方为1，从而可以同时传送接受2值(1比特)的信息。
本文所研究的QPSK系统与二进制的BPSK系统相比，具有以下特点:1.在传码率相同的情况下，四进制数字调制系统的信息速率是二进制系统的2倍。
2.在相同信息速率条件下，四进制数字调制系统的传码率是二进制系统的1/4倍，这一特点使得四进制码元宽度是二进制码元宽度的2倍，码元宽度的加大，可增加每个码元的能量，也可减小码间串扰的影响。
通信原理实验数字基带传输仿真实验
通信原理实验数字基带传输仿真实验本文记录的是一次通信原理实验，具体实验内容是数字基带传输仿真实验。
这个实验旨在让学生了解并掌握数字基带传输的基本原理、信号调制和调制解调的方法，并通过仿真实验加深对数字基带传输的理解。
实验步骤：第一步：实现数字基带信号的产生。
我们采用MATLAB编写代码来产生数字基带信号。
具体而言，我们可以选择产生脉冲振幅调制（PAM）、脉冲宽度调制（PWM）、脉冲频率调制（PFM）等各种调制方式。
第二步：实现数字基带信号的传输。
我们可以通过MATLAB编写代码，将数字基带信号在传输媒介中进行仿真。
具体而言，我们可以选择传输介质为AWGN信道、多径信道等，通过加入信噪比、码元传输速率、波特率等参数来模拟不同的传输环境。
第三步：实现数字基带信号的调制。
我们采用调制器进行数字信号的调制。
常见的数字调制方式有AM调制、FM调制、PM调制等。
此处我们选择了二进制相移键控（BPSK）调制来进行数字基带信号的调制。
第四步：实现数字基带信号的解调。
我们采用解调器来实现数字基带信号的解调。
常见的数字解调方式有包络检测法、抑制互调法等。
此处我们选择了直接判决法来进行数字基带信号的解调。
第五步：实现数字基带信号的重构。
我们通过将数字基带信号解调后还原成原始信号进行数字信号的重构。
此处我们需要通过MATLAB代码将解调后的数字信号还原成原始信号，并绘制出波形图进行对比分析。
实验结果：通过对仿真实验的分析，我们得出了一些结论。
首先，不同的数字基带信号相对应不同的调制方式，比如我们可以选择PAM调制来实现计算机通讯中的以太网传输。
其次，数字基带信号的传输受到了多种因素的影响，包括信道的噪声、信噪比、码元传输速率、波特率等。
第三，数字基带信号的解调方式有很多种，我们需要根据传输环境的不同来选择最适宜的解调方式。
最后，数字基带信号的重构是一个非常重要的环节，它能够让我们了解数字基带信号在传输过程中所带来的信息损失和失真情况。
基于Matlab的调制解调技术仿真
实验报告书------基于Matlab的调制解调技术仿真基于Matlab的调制解调技术仿真班级：姓名：学号：一、设计原理数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。
为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。
这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。
数字调制技术的两种方法：①利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理；②利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。
这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（PSK）基本的调制方式。
二、实验仪器1、电脑-MATLAB 一台三、实验目的1、掌握数字带通BPSK调制解调相关知识2、运用MATLAB进行编程实现BPSK的调制解调过程3、仿真输出调制前的基信号、调制后的BPSK信号，解调器在接收到信号后解调的各点的信号波形4、对仿真结果进行分析四、实验报告（填写相应原理，用MATLAB 实现仿真，列出仿真源程序，分析仿真输出结果，总结相关内容）1、BPSK 的调制原理如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于"同相"状态；如果其中一个开始得迟了一点，就可能不相同了。
如果一个达到正最大值时，另一个达到负最大值，则称为"反相"。
一般把信号振荡一次（一周）作为360度。
如果一个波比另一个波相差半个周期，我们说两个波的相位差180度，也就是反相。
当传输数字信号时，"1"码控制发0度相位，"0"码控制发180度相位。
载波的初始相位就有了移动，也就带上了信息。
相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。
BPSK调制的MATLAB仿真课程设计_实验报告
北京邮电大学移动通信课程设计实验报告目录一、背景 (4)二、基本要求 (4)三、设计概述 (4)四、Matlab设计流程图 (5)五、Matlab程序及仿真结果图 (6)1、生成m序列及m序列性质 (6)2、生成50位随机待发送二进制比特序列，并进行扩频编码 (7)3、对扩频前后信号进行BPSK调制，观察其时域波形 (9)4、计算并观察扩频前后BPSK调制信号地频谱 (10)5、仿真经awgn信道传输后，扩频前后信号时域及频域地变化 (11)6、对比经信道前后两种信号地频谱变化 (12)7、接收机与本地恢复载波相乘，观察仿真时域波形 (14)8、与恢复载波相乘后，观察其频谱变化 (15)9、仿真观察信号经凯萨尔窗低通滤波后地频谱 (16)10、观察经过低通滤波器后无扩频与扩频系统地时域波形 (17)11、对扩频系统进行解扩，观察其时域频域 (18)12、比较扩频系统解扩前后信号带宽 (19)13、比较解扩前后信号功率谱密度 (20)14、对解扩信号进行采样、判决 (21)15、在信道中加入2040~2050Hz窄带强干扰并乘以恢复载波 (24)16、对加窄带干扰地信号进行低通滤波并解扩 (25)17、比较解扩后信号与窄带强干扰地功率谱 (27)六、误码率simulink仿真 (28)1、直接扩频系统信道模型 (28)2、加窄带干扰地直扩系统建模 (29)3、用示波器观察发送码字及解扩后码字 (30)4、直接扩频系统与无扩频系统地误码率比较 (31)5、不同扩频序列长度下地误码率比较 (32)6、扩频序列长度N=7时，不同强度窄带干扰下地误码率比较 (33)七、利用Walsh码实现码分多址技术 (34)1、产生改善地walsh码 (35)2、产生两路不同地信息序列 (36)3、用两个沃尔什码分别调制两路信号 (38)4、两路信号相加，并进行BPSK调制 (39)5、观察调制信号频谱，并经awgn信道加高斯白噪和窄带强干扰 (40)6、接收机信号乘以恢复载波，观察时域和频域 (42)7、信号经凯萨尔窗低通滤波器 (43)8、对滤波后信号分别用m1和m2进行解扩 (44)9、对两路信号分别采样，判决 (45)八、产生随机序列Gold码和正交Gold码 (47)1、产生Gold码并仿真其自相关函数 (48)2、产生正交Gold码并仿真其互相关函数 (50)九、实验心得体会 (51)直接序列扩频系统仿真一、背景直接序列扩频通信系统(DSSS)是目前应用最为广泛地系统.在发送端，直扩系统将发送序列用伪随机序列扩展到一个很宽地频带上去，在接受端又用相同地扩频序列进行解扩，回复出原有信息.由于干扰信息与伪随机序列不相关，扩频后能够使窄带干扰得到有效地抑制，提高输出信噪比.系统框图如下图所示：二、基本要求：通过matlab建模，对直扩系统进行仿真，数据调制方式可以自由选择，可以使用基带信号，但最好能使用频带信号，信道为高斯白噪信道.要仿真出扩频前地信号地频偏，扩频后地信号频谱，过信道之后地频谱以及解扩之后地频谱.2.研究并仿真产生m序列，写出生成m序列地算法.验证直扩系统对窄带干扰地抑制能力，在信道中加入一个窄带强干扰，仿真出加了干扰后地频谱图和解扩后地频谱图，给出误码率等仿真图.4.在以上基础上仿真实现码分多址技术，使用Walsh码进行复用，实现多个信号同时传输.(选做)可选项：1.在信道中加入多径，使用rake接收来抗多径效应.2.产生除m序列之外地其他随机序列，如Gold码，正交Gold码等等.3.对比无扩频地系统地误码率.三、设计概述本次课设完成基本要求，并选作了可选项码分多址，Gold码及误码率对比.通过matlab 建模仿真了直扩系统BPSK调制地各点频偏及时域信号，并仿真了窄带强干扰对直扩系统地影响以及利用改善地W ALSH码实现码分多址技术.另外，通过matlab地simulink工具盒bertool工具仿真对比了直扩系统和无扩频系统地误码率.四、matlab设计流程图基本扩频系统仿真流程图五、matlab程序及仿真结果图1、生成m序列及m序列性质实验产生7位m序列，频率100Hz，模拟线性反馈移位寄存器序列，原理图如下：clear all。
基于MATLAB的BPSK通信系统仿真
基于MATLAB的BPSK通信系统仿真摘要：BPSK就是数字带通传输系统中一种非常重要的方式。
对于它的研究有助于我们掌握更加复杂的数字通信系统的传输方式。
本文介绍了BPSK的调制原理，以及使用Simulink对其调制过程进行仿真。
Simulink是MATLAB中最重要的组件之一，它向用户提供一个动态系统建模、仿真和综合分析的集成环境，在这种环境中，用户无须写大量的程序，而只要通过简单直观的鼠标操作，选取适当的库模块，就可以构造出复杂的仿真模型。
文章对于Simulink模块的库作了比较详细的介绍，如何进行模块参数的设置也作了说明。
使用Simulink仿真，直观简便，极大地减轻了工作量，简化了复杂的通信过程分析，也极大的方便了我们对于通信系统的研究与掌握。
关键词：BPSK；Simulink；仿真BPSK Communication System Simulation Based on MATAB Abstract:BPSK is the digital transmission system in a very important way. For its research helps us to grasp more complex transmission of a digital communication system. This article introduces the BPSK modulation principle, as well as its modulation process simulation using Simulink. Simulink is one of the most important components in MATLAB, it offers the user a comprehensive analysis of dynamic systems modeling, simulation and integration environments, in this environment, users do not have to write a lot of programs, but merely through simple, intuitive mouse operation, select the appropriate library module, you can construct sophisticated simulation models. Article for a Simulink block library for a more detailed introduction, described how to set module parameters. Using Simulink simulation, intuitive and easy, greatly reducing the workload, simplifying analysis of complex communication process, was also very handy for communication system of research and at our disposal.Key words:BPSK; Simulink; simulation目录1 序言 (1)1.1 数字通信系统研究的背景及意义 (1)1.2 数字通信的应用现状和未来的发展 (1)2 通信的基本知识 (4)2.1通信的基本概念 (4)2.2通信系统的通信过程与原理 (4)2.2.1通信系统的一般模型 (4)2.2.2模拟和数字通信系统 (6)3 BPSK的调制原理 (8)3.1数字调制的基础 (8)3.2 二进制相移键控（BPSK） (9)4 MATLAB的Simulink工具仿真 (11)4.1 MATLAB的简介 (11)4.2 Simulink的主要优点[7] (12)4.3 Simulink的模块库介绍 (12)4.4 Simulink简单模型的建立 (15)4.5 2PSK调制的仿真 (15)4.5.1 键控调制法 (15)4.5.2 模拟调制法 (19)4.6 仿真结果分析 (20)5 结论 (21)参考文献 (22)致谢 (23)1 序言1.1数字通信系统研究的背景及意义古代通信利用自然界的基本规律和人的基础感官（视觉，听觉等）可达性建立通信系统，是人类基于需求的最原始通信方式。
(完整版)PSK原理及MATLAB仿真
2PSK调制
2PSK系统的调制部分框图如下图所示
M序列 发生器
差分 编码
调相
PSK 频
外载波
放大器
外内
2分频
10MHZ
晶振
图2 2PSK调制部分框图
1. M序列发生器
实际的数字基带信号是随机的，为了实验 和测试方便，一般都是用M序列发生器产生 一个伪随机序列来充当数字基带信号源。按 照本原多项式f(x)=X5+X3+1组成的五级线性 移位寄存器，就可得到31位码长的M序列。
e w Acos( t )
2 psk
c
n
其中， n表示第n个符号的绝对相位： n
0发送‘1’时
发送‘0’时
• 2PSK信号波形如图5所示
图5
解调
2PSK系统的解调部分框图如下图所示
鉴相器
Ud1
低通
再生码
PSK入
Um2
压控
Ud
振荡器
环路
Ud
滤波器
模拟 相乘器
差分 译码
Um1
鉴相器
Ud2
低通
相对移相的调制规律是：每一个码元的载波相
位不是以固定的未调载波相位作基准的，而是以相 邻的前一个码元的载波相位来确定其相位的取值。 例如，当某一码元取“1”时，它的载波相位与前 一码元的载波同相；码元取“0”时，它的载波相 位与前一码元的载波反相。相对移相的波形如下图 所示。
U(t )
基带 信号
t
相应的信号波形如图1所示：
图1 正弦载波的三种键控波形
2. 三种键控的比较：
在大多数情况下，数字调制是利用数 字信号的离散值去键控载波。对载波的幅 度、频率或相位进行键控，便可获ASK、 FSK、PSK等。这三种数字调制方式在抗 干扰噪声能力和信号频谱利用率等方面， 以相干PSK的性能最好，目前已在中、高 速传输数据时得到广泛应用。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
clc
clear
num = 1e5;
fs = 1e4;
Rb = 1e3;
tnum = fs/Rb;
fc = Rb*8;
int = randint(1,num)*2-1;
int_c = upsample(int,tnum);
g = ones(1,tnum);
int_c = filter(g,1,int_c);
%»ù´ø
for snr_db = -5:10;
snr = 10^(snr_db/10);
sigma = sqrt(1/snr);
noise = randn(1,num*tnum)*sigma;
reciever_in = int_c+ noise;
reciever_panjue = reciever_in(tnum/2:tnum:end);
panjue = im2bw(reciever_panjue,0)*2-1;
pe(snr_db+6) = length(find(int~=panjue))/num;
end
semilogy(-5:10,pe,'*r');grid on;hold on;
snr_db = -5:10;
snr = 10.^(snr_db/10);
pe_lilun = 0.5*erfc(sqrt(snr/2));
semilogy(-5:10,pe_lilun,'-b');
%%%%Ƶ´ø´«Êä
t = 0:1/fs:num/Rb-1/fs;
carrier = cos(2*pi*fc*t);
int_c = int_c.*carrier;
[a1,b1] = butter(2,0.1);
for snr_db = -5:10;
snr = 10^(snr_db/10);
sigma = sqrt(0.5/snr);
noise = randn(1,num*tnum)*sigma;
reciever_in = int_c+ noise;
reciever_c = reciever_in.*carrier;
reciever_c = filter(a1,b1,reciever_c);
reciever_panjue = reciever_in(1:tnum:end);
panjue = im2bw(reciever_panjue,0)*2-1;
pe(snr_db+6) = length(find(int~=panjue))/num;
end
semilogy(-5:10,pe,'*g');
snr_db = -5:10;
snr = 10.^(snr_db/10);
pe_lilun = 0.5*erfc(sqrt(snr));
semilogy(-5:10,pe_lilun,'-k');