扩频信号的BPSK调制仿真.docx

扩频信号的BPSK调制仿真.docx

《扩频信号的BPSK调制仿真.docx》由会员分享，可在线阅读，更多相关《扩频信号的BPSK调制仿真.docx（15页珍藏版）》请在冰豆网上搜索。

扩频信号的BPSK调制仿真

《控制系统CAD与数字仿真》

2014～2015学年第一学期

学院（部）

电子电气工程学院

学号

0

姓名

纪辰

授课教师

陈剑雪

扩频信号的BPSK调制仿真

摘要：

随着数字信号处理技术的不断发展，数字化软件无线电接收机已经成为趋势，调制/解调技术是数字通信系统中的核心技术。

现代计算机科学技术快速发展，使得通信系统仿真的设计和分析过程变得相对直观和便捷，由此也使得通信系统仿真技术得到了更快的发展。

通信系统仿真具有广泛的适应性和极好的灵活性,有助于更好地研究通信系统性能。

本文介绍了数字化调制解调技术的现状发展及其应用，然后介绍了BPSK数字调制解调的理论基础，重点分析了BPSK数字调制和解调的原理。

本文利用MATLAB强大的仿真功能,在其仿真环境下建立了直扩系统和BPSK调制解调系统仿真模型，给出各路观察波形,证实了解调算法的可行性。

最后，本文对所做研究工作进行了总结，并且提出了今后的工作和研究方向。

关键词：

BPSK；调制解调器；MATLAB；直扩系统

ThesimulationofspreadspectrumsignalwithBPSKdigitalmodulation

Abstract：

Withdigitalsignalprocessingtechnologycontinuestoevolve,digitalsoftwareradioreceiverhasbecomeatrend,modulation/demodulationtechnologyisthecoretechnologyindigitalcommunicationsystem。

Therapiddevelopmentofmoderncomputerscienceandtechnology,makesthecommunicationssystemsimulationdesignandanalysisprocessbecomerelativelyintuitiveandconvenient,whichalsomakesthecommunicationsystemsimulationtechnologyhasbeenfasterdevelopment.Communicationsystemsimulationwithwideadaptabilityandexcellentflexibility,helpstobetterstudythecommunicationsystemperformance.

ThispaperintroducesthedigitalmodemtechnologysituationandintroduceddigitalmodulationanddemodulationofBPSKtheoreticalfoundation,includinganalysisoftheBPSKdigitalmodulationanddemodulationprinciple.

Inthispaper,basedonthepowerfulsimulationusingMATLABfunctioninitsenvironmentdesignedthespreadspectrumsignalsystemandtheBPSKmodulationdemodulationsystemsimulationmodel,andthroughtheconstellationconfirmedthatthedemodulationalgorithm.

Finally,thispapermadeasummaryoftheresearchwork,andproposedfutureworkandresearchdirections.

Keywords:

BPSK,Modem,MATLAB,Spreadspectrumsignalsystem

一、引言

研究背景及其意义

随着数字信号处理技术的不断发展，数字化软件无线电接收机已经成为趋势。

在卫星通信、移动通信、微波通信、光纤通信等现代通信系统中，信道中传输的都是数字已调信号，称它们为数字调制系统。

数字调制技术是数字通信系统中的核心技术，随着数字调制技术的出现，在有限的带宽内传输高速的数据已经成为可能，并且与过去使用的模拟调制相比，如调幅（AM）和调频（FM）、频移键控（FSK）、开关键控（OOK）、脉宽调制（PWM）、脉位调制（PPM）、脉幅调制（PAM）等技术相比有更高的可靠性和抗干扰性。

数字调制是正交幅度调制（QAM）、正交相移键控（QPSK）、二进制相移键控（BPSK）以及由这些技术派生的调制方法。

BPSK是最简单的二进制相移键控调制方法，其它更先进的调制方法大都由BPSK改进和增强，BPSK作为一种多进制的调相技术获得了广泛的应用。

直接序列扩频通信系统（DSSS）是目前应用最为广泛的系统。

在发送端，直扩系统将发送序列用伪随机序列扩展到一个很宽的频带上去，在接受端又用相同的扩频序列进行解扩，回复出原有信息。

由于干扰信息与伪随机序列不相关，扩频后能够使窄带干扰得到有效的抑制，提高输出信噪比。

数字化调制解调技术的现状发展及其应用

调制和解调是现代通信的重要手段。

调制就是用基带信号对载波波形的某些参数进行控制，使这些参量随基带信号的变化而变化；解调是调制的逆过程，它是从己调制的信号中恢复出原来调制信号的过程。

根据被调制的是模拟还是数字信号，调制技术分为模拟调制和数字调制两类。

模拟调制是对载波信号的参量进行连续调制，在接收端对载波信号的参量连续地进行估值；而数字调制都是用载波信号的某些离散状态来表征所传送的信息，在接收端也主要对载波信号的离散调制参量进行检测。

本论文主要是对数字调制解调技术的研究。

数字调制技术主要有三种方式：

振幅键控（ASK）、移频键控（FSK）和相移键控（PSK）。

五十年代末就已经出现了二相相移键控，继而为了提高频谱利用率，又出现了四相相移键控（QPSK），但这两种方式的已调波在码元转换时刻都可能产生150°的相位突跳，这使得频谱高频滚降慢，边带起伏大。

为了消除这种相位突跳，六十年代末,在QPSK的基础上又出现了交错正交移相键控（OQPSK）,这种调制方式虽然克服了180°的相位突跳，但仍然存在着90°的相位突跳。

为了彻底解决相位突跳问题，在七十年代提出了最小频移键控技术（MSK）,但这种技术和后来在此基础上出现的几种调制方式的相位特性仅局限于一个码元内进行,这大大限制了选择不同相位路径的可能性。

所以调制技术发展的过程，其实也正是己调波相位路径不断得到改进和完善的过程。

数字调制解调方式尽管种类繁多，但是随着数字技术的快速发展，将来定会出现更多先进的调制解调技术，谁都无法准确描述未来调制解调技术的发展。

总之，未来的调制解调技术都是朝着一个方向发展：

将不断地更新改进，使通信更高速和更可靠，传输速率更快，误码率更低。

MATLAB的简介

目前，在国际流行的科技应用软件中，数学类（区别于文字处理类和图像处理类）软件共有几十款之多。

从它们的数学处理的原始内核来看，不外乎两种类型：

数值计算型和数学分析型。

前者如MATLAB、Xmath等，它们对大量数据具有较强的管理、计算和可视化能力，运行效率高；后者如Mathematica、Maple等，它们擅长于符号计算，可以得到问题的解析符号解和任意精度解，但处理大量数据时速度较慢。

MATLAB是英文MatrixLaboratory（矩阵实验室）的缩写。

1980年前后，时任美国墨西哥大学计算机科学系主任的CleveMoler教授在给学习讲授线性代数课程时，想教学生使用当时流行的线性代数软件包（Linpack）和基于特征值计算的软件包（Eispack），但发现用其他高级语言编程极为不便，于是CleveMoler教授为学生编写了方便使用Linpack和Eispack的借口程序并命名为MATLAB,这就是MATLAB最早的雏形。

从MATLAB诞生开始，由于其高度的集成性及应用的方便性，在高校中受到了极大的欢迎。

由于它使用方便，能非常快的实现科研人员的设想，极大的节约了科研人员的时间，受到了大多数科研人员的支持，经过一代代人的努力，目前已发展到了版本。

MATLAB是一种解释性执行语言，具有强大的计算、仿真、绘图等功能。

由于它使用简单，扩充方便，尤其是世界上有成千上万的不同领域的科研工作者不停的在自己的科研过程中扩充MATLAB的功能，使其成为了巨大的知识宝库。

可以毫不夸张的说，哪怕是你真正理解了一个工具箱，那么就是理解了一门非常重要的科学知识。

科研工作者通常可以通过MATLAB来学习某个领域的科学知识，这就是MATLAB真正在全世界推广开来的原因。

目前的MATLAB版本已经可以方便的设计漂亮的界面，它可以像VB等语言一样设计漂亮的用户接口，同时因为有最丰富的函数库（工具箱），所以计算的功能实现也很简单，进一步受到了科研工作者的欢迎。

另外，MATLAB和其他高级语言也具有良好的接口，可以方便的实现与其他语言的混合编程，进一步拓宽了MATLAB的应用潜力。

可以说，MATLAB已经也很有必要成为大学生的必修课之一，掌握这门工具对学习各门学科有非常重要的推进作用。

二、BPSK调制解调原理

BPSK信号调制原理

二进制相移键控BPSK（BinaryPhaseShiftKeying）方式一般是键控的载波相位按基带脉冲序列的规律而改变的数字调制方式，也就是说，二进制的数字基带信号0与1分别用相干调制的载波的0与π相位的波形来表示。

其表达式由公式（2-1）给出：

（2-1）

其中{

}为双极性的二进制数字序列，

的取值为±1，

为二进制的符号间隔，

基带的发送成形滤波器的冲激响应，通常具有升余弦特性；

是调制载波的频率，

是调制载波的初始相位。

用BPSK调制方式时，因为发送端以某一个相位作为基准，所以在接收端也一定有这样一个固定的基准相位作为参考。

假如参考相位发生变化了，那么接收端恢复的信息也会出错，也就是存在“倒π”现象。

因此需要在接收端使用载波同步，才能够正确恢复出基带的信号。

BPSK信号的调制原理框图如图2-2所示，典型波形如图2-3所示。

图2-2BPSK调制原理图

图2-3发送码元为10011的BPSK波形

BPSK信号的频谱如图2-4所示，可以计算频谱效率，所谓频谱效率是指信号传输速率与所占带宽之比。

在BPSK中，信号码元为

，故信号传输速率为

，以频谱的主瓣宽度为传输带宽，忽略旁瓣的影响，则射频带宽为2/

，频谱效率为：

（每赫）

即每赫兹带宽传输s。

注意，这里是以射频带宽计算的，若以基带带宽来计算，那就是每赫兹1b/s。

图2-4BPSK的频谱

BPSK的调制器非常简单，只要把数字信号与载波相乘即可。

不过这里数字信号的“0”要用“-1”来表示（在数字通信中，符号“1”用“+1”来表示，“0”则用“-1”来表示）。

由图2-4可见，BPSK波形与信息代码之间的关系是“异变同不变”，即：

若本码元与前一码元相异，则本码元内BPSK信号的初相相对于前一码元内BPSK信号末相变化180°；否则不变。

BPSK信号解调原理

因为BPSK信号的幅度与基带信号无关，故不能用包络检波法而只能用相干解调法解调BPSK信号，在相干解调过程中需要用到与接收的BPSK信号同频同相的相干载波,相干接收机模型如图2-5所示：

图2-5BPSK相干接收机模型

具体的BPSK信号解调原理框图如图2-6所示。

图2-6BPSK解调原理框图

如图2-6给出了一种BPSK信号相干解调原理框图，图中经过带通滤波的信号在相乘器与本地载波相乘，在相干解调中，如何得到与接收的BPSK信号同频同相的相干载波是关键，然后用低通滤波器去除高频分量，再进行积分采样判决，判决器是按极性进行判决，得到最终的二进制信息。

假设相干载波的基准相位于BPSK信号的调制载波的基准相位一致。

但是，由于在BPSK信号的载波恢复过程中存在

o的相位迷糊（phaseambiguity），即恢复的本地载波与所需的相干载波可能同相，也可能反相，这种相位关系的不确定性将会造成解调出数字基带信号与发送的数字基带信号正好相反，即‘1’变为‘0’，‘0’变为‘1’，判决器输出数字信号全部出错。

这种现象称为BPSK方式的‘倒

’现象。

三、扩频信号的BPSK调制系统设计仿真

基本扩频系统仿真流程图

生成m序列及m序列性质

产生8位m序列，频率100Hz，模拟线性反馈移位寄存器序列，原理图如下：

由上图可以看出，8位m序列为1，-1，-1，-1，1，-1，1，1。

另外，自相关函数的图形比较尖锐，最大值为1，最小值为-1/8，符合理论结果。

生成50位随机待发送二进制比特序列及进行扩频编码

生成的信息码频率为100/8Hz，利用m序列编码后，频率变为100Hz。

对扩频前后信号进行BPSK调制并观察其时域波形

BPSK调制采用2kHz信号cos（2*2000*t）作为载波

可以看出，100/8Hz的无扩频信号每时由于序列极性变换产生相位变换，100Hz的扩频后调制信号每由于序列极性变换产生相位变换。

计算并观察扩频前后BPSK调制信号的频谱

对信号采用400000点fft计算，得到频谱

如图，扩频前信号主瓣宽度约为2*100/8=25Hz，扩频后，信号频谱展宽，主瓣1900~2100Hz约为200Hz，为无扩频信号频谱宽度的N=8倍，符合理论推算。

四、总结与展望

总结

从以上结果可以看出BPSK调制解调系统的抗噪声能力较强。

但是，在操作中发现Matlab仿真系统运行时响应时间相对较长，需要计算机运行很长时间才能完成仿真运算过程。

经过分析，我认为，除了算法导致了系统运算量大以外，BPSK系统本身的信号传输效率也不高，这也印证了通信原理教材中提到的BPSK系统的特点，即抗噪声干扰能力强，但数据传输效率低下。

基于MATLAB的BPSK调制解调系统仿真模型的设计与分析,为实际系统的构建提供了很好的依据。

由于BPSK优越的性能，随着EDA技术的进一步发展。

该模型有着一定的技术价值。

该模型也可供其他设计方式的数字通信系统所使用，有着广泛的实践性意义。

展望

论文主要在对BPSK调制解调系统进行了设计和仿真,但是基于时间比较仓促,自己研究的内容在很多地方还需要不断改进,有以下几个方面是可以继续改进或深入分析的:

（1）本文在深刻理解通信系统理论的基础上设计了大量的仿真模型,但是这些模型大部分是针对通信系统某一方面的性能而设计的,虽然能说明一定问题,但整体规模比较小,今后应该研究相对复杂的通信系统,按照协议的规范,构建较大规模的通信系统仿真模型;

（2）由于网络的发展,频带资源越来越紧张,因此需要寻找频谱利用率更高的调制解调方案,因此有必要研究QPSK、64QAM,256QAM,1024QAM,这是以后要努力的工作。

参考文献

[1]孔艳芳.16QAM解调模块的设计实现与性能验证[J].重庆大学，.

[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.

[6]王俊胜.全数字BPSK／QPSK解调器原理和应用[J].无线电通信技术，1992，18（4）；225—263.

[7]薛定宇，陈阳泉.基于MATLAB/Simulink的系统仿真技术与应用[M]．北京：

清华大学出版社，.

[8]何一.全数字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序列性质

clearall;

clc;

X1=0;X2=0;X3=1;

m=400;%重复50遍的8位单极性m序列

fori=1:

m

Y3=X3;Y2=X2;Y1=X1;

X3=Y2;X2=Y1;

X1=xor（Y3,Y1）;

L（i）=Y1;

end

fori=1:

m

M（i）=1-2*L（i）;%将单极性m序列变为双极性m序列

end

k=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=abs（ym）;%求双极性m序列频谱

fm=（1:

2048）*200/2048;

plot（fm,magm（1:

2048）*2/4096）;

title（'双极性8位M序列的频谱'）

axis（[90,140,0,]）;

[a,b]=xcorr（M,'unbiased'）;

subplot（3,1,3）%求双极性m序列自相关函数

plot（b,a）;

axis（[-20,20,,]）;

title（'双极性8位M序列的自相关函数'）;

生成50位随机待发送二进制比特序列及进行扩频编码

N=50;a=0;

x_rand=rand（1,N）;%产生50个0与1之间随机数

fori=1:

N

ifx_rand（i）>=%大于等于的取1，小于的取0

x（i）=1;a=a+1;

elsex（i）=0;

end

end

t=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;

fori=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）;

end

s（l）=0;

fori=1:

400%扩频后，码率变为100/8*8=100Hz

s（i）=xor（L（i）,y（i））;

end

tt=0:

8*N-1;

stem（tt,s）;

axis（[0,400,0,1]）;

title（'扩频后的待发送序列码'）;

对扩频前后信号进行BPSK调制并观察其时域波形

figure（3）

subplot（2,1,2）

fs=2000;

ts=0:

:

;%为了使信号看起来更光滑，作图时采样频率为100kHz

%ps=cos（2*pi*fs*ts）;

s_b=rectpulse（s,1000）;%将冲激信号补成矩形信号

s_bpsk=（1-2.*s_b）.*cos（2*pi*fs*ts）;%扩频后信号BPSK调制时域波形

plot（ts,s_bpsk）;

xlabel（'s'）;

axis（[,,,]）

title（'扩频后bpsk信号时域波形'）;

subplot（2,1,1）

s_bb=rectpulse（x,8000）;

s_bpskb=（1-2.*s_bb）.*cos（2*pi*fs*ts）;%无扩频信号BPSK调制时域波形

plot（ts,s_bpskb）;

xlabel（'s'）;

axis（[,,,]）;

title（'扩频前bpsk信号时域波形'）

计算并观察扩频前后BPSK调制信号的频谱

figure（4）

N=400000;

ybb=fft（s_bpskb,N）;%无扩频信号BPSK调制频谱

magb=abs（ybb）;

fbb=（1:

N/2）*100000/N;

subplot（2,1,1）

plot（fbb,magb（1:

N/2）*2/N）;

axis（[1700,2300,0,]）;

title（'扩频前调制信号频谱'）;

xlabel（'Hz'）;

subplot（2,1,2）

yb=fft（s_bpsk,N）;%扩频信号BPSK调制频谱

mag=abs（yb）;

fb=（1:

N/2）*100000/N;

plot（fb,mag（1:

N/2）*2/N）;

axis（[1700,2300,0,]）;

title（'扩频后调制信号频谱'）;

xlabel（'Hz'）;