扩频信号的BPSK调制仿真

1、控制系统CAD与数字仿真20142015 学年 第 一 学期学院（部）电子电气工程学院学号021211229姓名纪辰授课教师陈剑雪扩频信号的BPSK调制仿真摘要：随着数字信号处理技术的不断发展， 数字化软件无线电接收机已经成为 趋势，调制/解调技术是数字通信系统中的核心技术。现代计算机科学技术快 速发展，使得通信系统仿真的设计和分析过程变得相对直观和便捷，由此也使得通信系统仿真技术得到了更快的发展。 通信系统仿真具有广泛的适应性和极 好的灵活性，有助于更好地研究通信系统性能。本文介绍了数字化调制解调技术的现状发展及其应用，然后介绍了BPSK数字调制解调的理论基础，重点分析了 BPSK 数字调制

2、和解调的原理。本文利用 MATLAB 强大的仿真功能，在其仿真环境下建立了直扩系统和BPSK 调制解调系统仿真模型，给出各路观察波形，证实了解调算法的可行性。最后，本文对所做研究工作进行了总结，并且提出了今后的工作和研究方 向。关键词：BPSK 调制解调器；MATLAB;直扩系统The simulation of spread spectrum signalwith BPSK digital modulationAbstract: With digital sig nal process ing tech no logy continues to evolve, digitalsoftware

3、 radio receiver has become a trend, modulati on /demodulati on tech no logyis the core tech no logy in digital com muni cati on system 。 The rapid developme ntof moder n computer scie nee and tech no logy, makes the com muni cati ons systemsimulati on desig n and an alysis process become relatively

4、in tuitive and convenient,which also makes the com muni cati on system simulatio n tech no logy has bee nfaster development. Communication system simulation with wide adaptability andexcelle nt flexibility, helps to better study the com muni catio n system performa nee.This paper in troduces the dig

5、ital modem tech no logy situati on anddevelopment.Then introduced digital modulation and demodulation of BPSKtheoretical foundation, including analysis of the BPSK digital modulation anddemodulati on prin ciple.In this paper, based on the powerful simulatio n using MATLAB fun cti on in its envir onm

6、ent desig ned the spread spectrum sig nal system and the BPSK modulatio ndemodulati on system simulatio n model, and through the con stellatio n con firmedthat the demodulatio n algorithm.Fin ally, this paper made a summary of the research work, and proposed futurework and research directi ons.Key w

7、ords: BPSK, Modem, MATLAB, Spread spectrum sig nal system目录一、弓丨言. -3 -1.1 研究背景及其意义 . -3 -1.2 数字化调制解调技术的现状发展及其应用 . -4 -1.3 MATLAB 的简介. -4 -二、BPSK 调制解调原理 . -5 -2.1 BPSK 信号调制原理 . -5 -2.2 BPSK 信号解调原理 . -7 -三、 扩频信号的 BPSK 调制系统设计仿真 . -7 -3.1 基本扩频系统仿真流程图 . -7 -3.2 生成 m 序列及 m 序列性质. 83.3 生成 50 位随机待发送二进制比特序列及进

8、行扩频编码 83.4 对扩频前后信号进行 BPSK 调制并观察其时域波形 .83.5 计算并观察扩频前后 BPSK 调制信号的频谱.8四、总结与展望. 84.1 总结. 84.2 展望. 错误！未定义书签。参考文献. 错误！未定义书签。附录. 10、弓I言1.1 研究背景及其意义随着数字信号处理技术的不断发展，数字化软件无线电接收机已经成为趋势。在卫星 通信、移动通信、微波通信、光纤通信等现代通信系统中，信道中传输的都是数字已调信 号，称它们为数字调制系统。数字调制技术是数字通信系统中的核心技术，随着数字调制 技术的出现，在有限的带宽内传输高速的数据已经成为可能，并且与过去使用的模拟调制 相比

9、，如调幅(AM)和调频(FM)、频移键控(FSK)、开关键控(00K)、脉宽调制(PWM)、脉 位调制(PPM)、脉幅调制(PAM)等技术相比有更高的可靠性和抗干扰性。数字调制是正交幅度调制(QAM)、正交相移键控(QPSK)、二进制相移键控(BPSK)以 及由这些技术派生的调制方法。BPSK是最简单的二进制相移键控调制方法，其它更先进 的调制方法大都由BPSK改进和增强，BPSK作为一种多进制的调相技术获得了广泛的应 用。直接序列扩频通信系统(DSSS)是目前应用最为广泛的系统。在发送端，直扩系统将发 送序列用伪随机序列扩展到一个很宽的频带上去，在接受端又用相同的扩频序列进行解 扩，回复出原

10、有信息。由于干扰信息与伪随机序列不相关，扩频后能够使窄带干扰得到有 效的抑制，提高输出信噪比。1.2 数字化调制解调技术的现状发展及其应用调制和解调是现代通信的重要手段。调制就是用基带信号对载波波形的某些参数进行 控制，使这些参量随基带信号的变化而变化；解调是调制的逆过程，它是从己调制的信号 中恢复出原来调制信号的过程。根据被调制的是模拟还是数字信号，调制技术分为模拟调 制和数字调制两类。模拟调制是对载波信号的参量进行连续调制，在接收端对载波信号的 参量连续地进行估值；而数字调制都是用载波信号的某些离散状态来表征所传送的信息，在接收端也主要对载波信号的离散调制参量进行检测。本论文主要是对数字调

11、制解调技术 的研究。数字调制技术主要有三种方式：振幅键控(ASK)、移频键控(FSK)和相移键控(PSK)。五十年代末就已经出现了二相相移键控，继而为了提高频谱利用率，又出现了四相相 移键控(QPSK)，但这两种方式的已调波在码元转换时刻都可能产生150的相位突跳，这使得频谱高频滚降慢，边带起伏大。为了消除这种相位突跳，六十年代末，在QPSK的基础上又出现了交错正交移相键控(OQPSK),这种调制方式虽然克服了180的相位突跳，但仍 然存在着90的相位突跳。为了彻底解决相位突跳问题，在七十年代提出了最小频移键控 技术(MSK),但这种技术和后来在此基础上出现的几种调制方式的相位特性仅局限于一个

12、码 元内进行，这大大限制了选择不同相位路径的可能性。所以调制技术发展的过程，其实也 正是己调波相位路径不断得到改进和完善的过程。数字调制解调方式尽管种类繁多，但是随着数字技术的快速发展，将来定会出现更多 先进的调制解调技术，谁都无法准确描述未来调制解调技术的发展。总之，未来的调制解 调技术都是朝着一个方向发展：将不断地更新改进，使通信更高速和更可靠，传输速率更 快，误码率更低。1.3 MATLAB 的简介目前，在国际流行的科技应用软件中，数学类(区别于文字处理类和图像处理类)软 件共有几十款之多。从它们的数学处理的原始内核来看，不外乎两种类型：数值计算型和 数学分析型。前者如MATLAB、Xm

13、ath等，它们对大量数据具有较强的管理、计算和可视 化能力，运行效率高；后者如Mathematica Maple等，它们擅长于符号计算，可以得到问 题的解析符号解和任意精度解，但处理大量数据时速度较慢。MATLAB是英文Matrix Laboratory（矩阵实验室）的缩写。1980年前后，时任美国墨 西哥大学计算机科学系主任的Cleve Moler教授在给学习讲授线性代数课程时， 想教学生 使用当时流行的线性代数软件包 （Linpack）和基于特征值计算的软件包（Eispack），但发 现用其他高级语言编程极为不便，于是Cleve Moler教授为学生编写了方便使用Lin pack和Eisp

14、ack的借口程序并命名为MATLAB,这就是MATLAB最早的雏形。从MATLAB诞生开始，由于其高度的集成性及应用的方便性，在高校中受到了极大 的欢迎。由于它使用方便，能非常快的实现科研人员的设想，极大的节约了科研人员的时 间，受到了大多数科研人员的支持，经过一代代人的努力，目前已发展到了7.X版本。MATLAB是一种解释性执行语言，具有强大的计算、仿真、绘图等功能。由于它使用简 单，扩充方便，尤其是世界上有成千上万的不同领域的科研工作者不停的在自己的科研过 程中扩充MATLAB的功能，使其成为了巨大的知识宝库。可以毫不夸张的说，哪怕是你 真正理解了一个工具箱，那么就是理解了一门非常重要的科

15、学知识。科研工作者通常可以 通过MATLAB来学习某个领域的科学知识，这就是MATLAB真正在全世界推广开来的原 因。目前的MATLAB版本已经可以方便的设计漂亮的界面，它可以像VB等语言一样设计 漂亮的用户接口，同时因为有最丰富的函数库（工具箱），所以计算的功能实现也很简单，进一步受到了科研工作者的欢迎。另外，MATLAB和其他高级语言也具有良好的接口，可以方便的实现与其他语言的混合编程，进一步拓宽了MATLAB的应用潜力。可以说，MATLAB已经也很有必要成为大学生的必修课之一，掌握这门工具对学习各门学科有非 常重要的推进作用。二、BPSK调制解调原理2.1 BPSK 信号调制原理二进制相

16、移键控BPSK（ Bi nary Phase Shift Keyi ng）方式一般是键控的载波相位按基带脉冲序列的规律而改变的数字调制方式，也就是说，二进制的数字基带信号0与1分别用相干调制的载波的0与n相位的波形来表示。其表达式由公式（2-1）给出：s（t）=a.gT（t -nTb）cos（讥R）（2-1）n wj其中 An为双极性的二进制数字序列，An的取值为土1，Tb为二进制的符号间隔，gT（t）基带的发送成形滤波器的冲激响应，通常具有升余弦特性；,是调制载波的频率，E 是调制载波的初始相位。用BPSK调制方式时，因为发送端以某一个相位作为基准，所以在接收端也一定有这 样一个固定的基准相

17、位作为参考。假如参考相位发生变化了，那么接收端恢复的信息也会 出错，也就是存在“倒n”现象。因此需要在接收端使用载波同步，才能够正确恢复出基 带的信号。BPSK言号的调制原理框图如图2-2所示，典型波形如图2-3所示。图2-2 BPSK调制原理图的频谱如图2-4所示，可以计算频谱效率，所谓频谱效率是指信号传输速率与所占带宽之 比。在BPS!中，信号码元为Tb，故信号传输速率为fb=1/仏，以频谱的主瓣宽度为传输带宽，忽略旁瓣的影响，则射频带宽为2/Tb，频谱效率为:信号传输速率即每赫兹带宽传输0.5b/s。注意，这里是以射频带宽计算的，若以基带带宽来计算，那就 是每赫兹1 b/s。图2-4 B

18、PSK的频谱BPSK勺调制器非常简单，只要把数字信号与载波相乘即可。不过这里数字信号的“0”要用“-1”来表示(在数字通信中，符号“1”用“+T来表示，“0”则用“-1”来 表示)。由图2-4可见，BPSK波形与信息代码之间的关系是“异变同不变”，即：若本码 元与前一码元相异，则本码元内BPSK言号的初相相对于前一码元内BPSK言号末相变化180;否则不变。2.2 BPSK 信号解调原理因为BPSK言号的幅度与基带信号无关，故不能用包络检波法而只能用相干解调法解调BPSK言号， 在相干解调过程中需要用到与接收的BPSK言号同频同相的相干载波,相干接 收机模型如图2-5所示：图2-5 BPSK相

19、干接收机模型具体的BPSK信号解调原理框图如图2-6所示。键，然后用低通滤波器去除高频分量，再进行积分采样判决，判决器是按极性进行判决,得到最终的二进制信息。假设相干载波的基准相位于BPSK信号的调制载波的基准相位一致。但是，由于在BPSK言号的载波恢复过程中存在180o的相位迷糊(phase ambiguity)， 即恢复的本地载波与所需的相干载波可能同相，也可能反相，这种相位关系的不确定性将 会造成解调出数字基带信号与发送的数字基带信号正好相反，即1变为0,0变为1,判决器输出数字信号全部出错。这种现象称为BPSK方式的倒二现象。三、扩频信号的BPSK调制系统设计仿真eBP SK（t）如图

20、2-6给出了一种BPSKI号相干解调原理框图与本地载波相乘，在相干解调中，如何得到与接收的图中经过带通滤波的信号在相乘器BPSK信号同频同相的相干载波是关3.1基 本 扩 频 系 统 仿 真 流 程 图e输出卜对信号采用400000点fft计算，得到频谱如图，扩频前信号主瓣宽度约为2*100/8=25Hz，扩频后，信号频谱展宽,主瓣19002100Hz约为200Hz,为无扩频信号频谱宽度的N=8倍，符合理论推算。四、总结与展望4.1 总结从以上结果可以看出BPSKS制解调系统的抗噪声能力较强。但是，在操作 中发现Matlab仿真系统运行时响应时间相对较长，需要计算机运行很长时间才 能完成仿真运

21、算过程。经过分析，我认为，除了算法导致了系统运算量大以 外，BPSK系统本身的信号传输效率也不高，这也印证了通信原理教材中提到的BPSK系统的特点，即抗噪声干扰能力强，但数据传输效率低下。基于MATLA的BPSK调制解调系统仿真模型的设计与分析，为实际系统的构建提供了很好的依 据。由于BPSK优越的100/8HZ 二进制比特信息100 Hz 8 位双极性 m 序列3.2生成m 序列及- 由上图可以看出，8位m序列为1,-1，-2000Hz 载波 COS4000nt关函数的图形比较尖锐，最大值为1,最小值为,1,-1,1, 1。另外，自相3.3 生成 50 位随机待发送二进制比特序列及进行扩频编

22、码一 BPSK 调制信号生成的信息码频率为100/8Hz|，利用m序列编码后，频率变为100Hzo3.4 对扩频前后信号进行 BPRK 调制并观察其时域波形、BPSKM制采用2kHz信号cos（2*2000*t）观察波形波换,可以看出，100/8HZ的无扩频信号每0.08s时由于序列极性变换产生相位变100Hz的扩频后调制信号每0.01s由于序列极性变换产生相位变换。3.5 计算并观察扩频前后 BPSK 调制信号的频谱100Hz1 !，符合理论结果。性能，随着EDA技术的进一步发展。该模型有着一定的技 术价值。该模型也可供其他设计方式的数字通信系统所使用，有着广泛的实践 性意义。4.2 展望论

23、文主要在对BPSK调制解调系统进行了设计和仿真，但是基于时间比较仓 促，自己研究的内容在很多地方还需要不断改进，有以下几个方面是可以继续改 进或深入分析的：(1)本文在深刻理解通信系统理论的基础上设计了大量的仿真模型，但是这 些模型大部分是针对通信系统某一方面的性能而设计的，虽然能说明一定问题，但整体规模比较小，今后应该研究相对复杂的通信系统，按照协议的规范，构建较 大规模的通信系统仿真模型；(2)由于网络的发展，频带资源越来越紧张，因此需要寻找频谱利用率更高的 调制解调方案,因此有必要研究QPSK 64QAM,256QAM,1024QA这是以后要努力 的工作。参考文献1 孔艳芳.16QAM

24、解调模块的设计实现与性能验证J.重庆大学，2011.04.2 常红霞基于 FPGA 的 QDPSK 调制解调技术的研究及实现J.南京理工大学，2007 . 06.3 昆仑，郭黎利.全数字 BPSK 调制解调器J.哈尔滨工程大学学报.2000(4):13-19.4 沈兰荪，李智群，张文忠等.调制解调的数字实现J.电信科学，1993,9(6):27-32.5 陈皓.基于 FPGA 的调制解调器的研究和设计J.华中师范大学，2007 . 05.王俊胜.全数字 BPSK / QPSK 解调器原理和应用J.无线电通信技术，1992，18；225263.7薛定宇，陈阳泉.基于 MATLAB/Simulin

25、k的系统仿真技术与应用 M 北京：清华大学岀版社，2011.02.何一.全数字 BPSK/ QPSK 解调机同步技术研究J.解放军信息工程大学，2003 . 05.9杨丽.基于 MATLAB 的通信系统仿真研究J.南京信息工程大学，2006 . 05.10 李贺冰.Simulink 通信仿真教程M.第 1 版.北京：国防工业岀版社，2006.11 李建新，刘乃安.现代通信系统分析与仿真-MATLAB 通信工具箱M.西安：西安电子科技大学岀版社,2000.12 苏伟.基于软件无线电的 BPSK/ QPSK 解调器的研究J.北京化工大学，2008 . 06.附录生成 m 序列及 m 序列性质cle

26、ar all;clc;X1=0;X2=0;X3=1;m=400;%重复50遍的8位单极性m序列for i=1:mY3=X3; Y2=X2; Y1=X1;X3=Y2; X2=Y1;X1=xor(Y3,Y1);L(i)=Y1;endfor i=1:mM(i)=1-2*L(i);%将单极性m序列变为双极性m序列endk=1:1:m;figure(1)subplot(3,1,1)%做m序列图stem(k-1,M);axis(0,8,-1,1);xlabel(k);ylabel(M序列)；title(移位寄存器产生的双极性8位M序列);subplot(3,1,2)ym=fft(M,4096);magm=

27、abs(ym);%求双极性m序列频谱fm=(1:2048)*200/2048;plot(fm,magm(1:2048)*2/4096);title(双极性8位M序列的频谱)axis(90,140,0,0.1);a,b=xcorr(M, un biased);subplot(3,1,3)%求双极性m序列自相关函数plot(b,a);axis(-20,20,-0.5,1.2);title(双极性8位M序列的自相关函数);生成 50 位随机待发送二进制比特序列及进行扩频编码泸生50个0与1之间随机数%大于等于0.5的取1，小于0.5的取0endt=0:N-1;figure(2)%故信息码图subplot(2,1,1)stem(t,x);title(扩频前待发送二进制信息序列);tt=0:349;subplot(2,1,2)l=1:8*N;y(l)=0; for i=1:N k=8*i-7;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1 ;y(k)=x(i);k=k+1;y(k)=x(i);ends(l)=0;for i=1:400频后，码率变为100/8*8=100Hzs(i)=xor(L(i),y(i);end tt=0