通信原理课程设计BPSK调制与解调1

1、摘要 数字通信系统是当代通信领域的主流，在社会生活各个方面占据重要地位。bpsk作为数字通信系统中的一种简单基础的调制解调方法，抗干扰能力强，容易仿真实现。本文通过bpsk的仿真，希望学习到数字通信的基础知识，为以后的学习打下基础。本文介绍了数字化调制解调技术的现状发展及其应用，通信系统仿真软件matlab中的一种可视化仿真工具simulink；然后介绍了bpsk数字调制解调的理论基础，包括数字带通传输分类以及重点分析了bpsk数字调制和解调的原理。本文在深刻理解通信系统理论的基础上，利用matlab强大的仿真功能,在simulink仿真环境下设计了bpsk调制解调系统仿真模型，给出各路观察波

2、形，证实了解调算法的可行性。 关键词：bpsk；调制解调器；matlab ；蒙特卡洛分析；目录一、课程设计目的及内容31.1、课程设计的目的31.2课程设计的内容3二、bpsk仿真设计思路42.1 相移键控系统概述42.2数字带通传输分类42.3 bpsk信号调制/解调原理42.3.1 bpsk信号调制原理42.3.2 bpsk 信号解调原理6三、matlab软件简介8四、bpsk调制解调的matlab仿真94.1 bpsk调制的数学模型94.2 bpsk解调的原理94.3 实验程序94.4 仿真波形图：15五、 总体系能分析 19六、设计总结20七、参考文献21致谢22一、课程设计目的及内容

3、1.1、课程设计的目的通过本课程的学习我们不仅能加深理解和巩固理论课上所学的有关 pcm编码和解码的基本概念、基本理论和基本方法，而且能锻炼我们分析问题和解决问题的能力；同时对我们进行良好的独立工作习惯和科学素质的培养，为今后参加科学工作打下良好的基础。本课程设计主要研究8psk信号的调制解调原理性能分析。通过完成本课题的设计，拟达到以下目的：1学习如何利用计算机仿真方法和技术对通信系统的理论知识进行验证，并学会搭建简单的系统模型；2掌握matlab7.0的基础知识，熟悉matlab进行通信系统仿真中各个常用模块的使用方法；3通过系统仿真加深对通信课程理论知识的理解。通过该课题的设计与仿真，可

4、以提高学生综合应用所学基础知识的能力和计算机编程的能力，为今后的学习和工作积累经验。1.2课程设计的内容根据题目要求，查阅相关资料，掌握数字带通的 bpsk 调制解调的相关知识。学习 matlab软件，掌握 matlab各种函数的使用。在此基础上，完成以下实验要求：1) 设计系统整体框图及数学模型。2) 运用 matlab进行编程，实现 bpsk 的调制解调过程的仿真。其中包括信源、bpsk信号的产生，信道噪声的加入，bpsk信号的载波提取和相干解调。3) 系统性能的分析包括信号带宽，波形对比以及误码率的计算。二、bpsk仿真设计思路2.1 相移键控系统概述相移键控是目前扩频系统中大量使用的调

5、制方式，也是和扩频技术结合最成熟的调制技术，原则上看是一种线性调制。从基带变换到中频以及射频，中间的频谱搬移和信号放大需要一个要求较高的线性信道，因而，设计要求较高。相移键控系统中，有待传输的基带数字脉冲控制着载波相位的变化，从而形成振幅与频率不变，而相位取离散值变化的已调波。2.2数字带通传输分类数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。数字带通传输中一般利用

6、数字信号的离散取值特点通过开关键控载波，从而实现数字调制，比如对载波的振幅、频率和相位进行键控可获得振幅键控（ask）、频移键控（fsk）和相移键控（psk）。这三种数字调制方式在抗干扰噪声能力和信号频谱利用率等方面，以相干psk的性能最好，目前已在中、高速传输数据时得到广泛的应用。2.3 bpsk信号调制/解调原理2.3.1 bpsk信号调制原理二进制相移键控 bpsk（binary phase shift keying）方式一般是键控的载波相位按基带脉冲序列的规律而改变的数字调制方式，也就是说，二进制的数字基带信号 0 与 1 分别用相干调制的载波的 0 与相位的波形来表示。其表达式由公式

7、（1-1）给出： （1-1）其中为双极性的二进制数字序列，的取值为 1，为二进制的符号间隔，基带的发送成形滤波器的冲激响应，通常具有升余弦特性；是调制载波的频率，是调制载波的初始相位。用 bpsk 调制方式时，因为发送端以某一个相位作为基准，所以在接收端也一定有这样一个固定的基准相位作为参考。假如参考相位发生变化了，那么接收端恢复的信息也会出错，也就是存在“倒”现象。因此需要在接收端使用载波同步，才能够正确恢复出基带的信号。bpsk信号的调制原理框图如图2-1所示，典型波形如图2-2所示。图2-1 bpsk调制原理图图图2-2 发送码元为1 0 0 1 1的bpsk波形bpsk信号的频谱如图2

8、-3所示，可以计算频谱效率，所谓频谱效率是指信号传输速率与所占带宽之比。在bpsk中，信号码元为，故信号传输速率为，以频谱的主瓣宽度为传输带宽，忽略旁瓣的影响，则射频带宽为2/，频谱效率为： （每赫） 即每赫兹带宽传输0.5b/s。注意，这里是以射频带宽计算的，若以基带带宽来计算，那就是每赫兹1 b/s。图2-3 bpsk的频谱bpsk的调制器非常简单，只要把数字信号与载波相乘即可。不过这里数字信号的“0”要用“-1”来表示（在数字通信中，符号“1”用“+1”来表示，“0”则用“-1”来表示）。由图2-3可见，bpsk波形与信息代码之间的关系是“异变同不变”，即：若本码元与前一码元相异，则本码

9、元内bpsk信号的初相相对于前一码元内bpsk信号末相变化180；否则不变。2.3.2 bpsk 信号解调原理因为bpsk信号的幅度与基带信号无关，故不能用包络检波法而只能用相干解调法解调bpsk信号，在相干解调过程中需要用到与接收的bpsk信号同频同相的相干载波,相干接收机模型如图2-4所示： 图2-4 bpsk相干接收机模型具体的bpsk信号解调原理框图如图2-5所示。带通滤波器相乘器 低通滤波器抽样判决图2-5 bpsk解调原理框图如图2-5给出了一种bpsk信号相干解调原理框图，图中经过带通滤波的信号在相乘器与本地载波相乘，在相干解调中，如何得到与接收的bpsk信号同频同相的相干载波是

10、关键，然后用低通滤波器去除高频分量，再进行积分采样判决，判决器是按极性进行判决，得到最终的二进制信息。假设相干载波的基准相位于bpsk信号的调制载波的基准相位一致。但是，由于在bpsk信号的载波恢复过程中存在的相位迷糊（phase ambiguity），即恢复的本地载波与所需的相干载波可能同相，也可能反相，这种相位关系的不确定性将会造成解调出数字基带信号与发送的数字基带信号正好相反，即1变为0，0变为1，判决器输出数字信号全部出错。这种现象称为bpsk方式的倒现象。载波同步器从bpsk信号中提取的相干载波可能与接收信号的载波同相，也可能反相，称此为相干载波的相位模糊现象。如果收到的信号与载波信

11、号同相，则相乘为正值，积分采样后必为一大于0的值，即可判决为“1”。如果收到的信号与参考信号相反，则相乘之后必为负值，积分采样后判决为“0”，因此解调完成。具体波形如图2-6所示。图2-6 bpsk解调信号示意图三、matlab简介 matlab软件是美国math works公司的产品，matlab是英文matrixlaborat -ory(矩阵实验室)的缩写。matlab软件系列产品是一套高效强大的工程技术数值运算和系统仿真软件，广泛应用于当今的航空航天、汽车制造、半导体制造、电子通信、医学研究、财经研究和高等教育等领域，被誉为“巨人肩膀上的工具”。研发人员借助matlab软件能迅速测试设想

12、构想，综合评测系统性能，快速设计更好方案来确保更高技术要求。同时matlab也是国家教委重点提倡的一种计算工具。 matlab主要由c语言编写而成，采用lapack 为底层支持软件包。 matlab的编程非常简单，它有着比其他任何计算机高级语言更高的编程效率、更好的代码可读性和移植性，以致被誉为“第四代”计算机语言，matlab是所有mathworks公司产品的数值分析和图形基础环境。此外matlab 还拥有强大的2d和3d甚至动态图形的绘制功能，这样用户可以更直观、更迅速的进行多种算法的比较，从中找出最好的方案。 从通信系统分析与设计、滤波器设计、信号处理、小波分析、神经网络到控制系统、模糊

13、控制等方面来看，matlab提供了大量的面向专业领域的工具箱。通过工具箱，以往需要复杂编程的算法开发任务往往只需一个函数就能实现，而且工具箱是开放的可扩展集，用户可以查看或修改其中的算法，甚至开发自己的算法。目前， matlab已经广泛地应用于工程设计的各个领域，如电子、通信等领域；它已成为国际上最流行的计算机仿真软件设计工具。现在的matlab不再仅仅是一个矩阵实验室，而是一种实用的、功能强大的、不断更新的高级计算机编程语言。现在从电子通信、自动控制图形分析处理到航天工业、汽车工业，甚至是财务工程。matlab都凭借其强大的功能获得了极大的用武之地。广大学生可以使用matlab来帮助进行信号

14、处理、通信原理、线性系统、自动控制等课程的学习；科研工作者可以使用matlab进行理论研究和算法开发；工程师可以使用matlab进行系统级的设计与仿真。四、bpsk调制解调的matlab仿真4.1 bpsk调制的数学模型由于bpsk的两种码元的波形相同，极性相反，故bpsk信号可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘： 其中，即s(t)为双极性全占空(非归零)矩形脉冲序列。4.2 bpsk解调的原理由于psk信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息采用相干解调法来解调信号。4.3 实验程序clear allclose allclcnum=10; %码

15、元个数tnum=200;%码元长度n=num*tnum;%10个码元整体长度a=randint(1,num,2); %产生1行num列的矩阵，矩阵内0和1随机出现fc=0.5; %载波频率为0.5t=0:0.05:9.99;%t从0到9.99，间隔为0.05s=;c=;for i=1:num %i从1到10循环 if(a(i)=0) a=zeros(1,tnum); %i=0时，产生一个码元长度为tnum（200）的0码元 else a=ones(1,tnum); %i=1时，产生一个码元长度为tnum（200）的1码元 end s=s a; %s为随机基带信号 cs=sin(2*pi*fc*

16、t); c=c cs; %c为载波信号end%采用模拟调制方法得到调制信号s_nrz=;for i=1:num %i从1到num（10）循环 if(a(i)=0) a=ones(1,tnum); %i=0时，产生一个码元长度为tnum（200）的1码元 else a=-1*ones(1,tnum); %i非0时，产生一个码元长度为tnum（200）的-1码元 end s_nrz=s_nrz,a; %s_nrz为双极性非归零码ende=s_nrz.*c; %e为bpsk调制信号figure(1); %图1subplot(3,2,1); %图1分为32部分的第一部分plot(s); %作s（基带信

17、号）的波形图grid on;axis(0 n -2 2); %横轴长度为0到n，纵轴范围为-2到+2xlabel(基带信号s(t); %x轴的注释ylabel(基带信号幅值); %y轴的注释subplot(323);plot(c);grid on;axis(0 n -2 2);xlabel(bpsk载波信号);ylabel(bpsk载波信号幅值); %作c（bpsk载波信号）的波形图subplot(325);plot(e);grid on;axis(0 n -2 2);xlabel(bpsk调制信号);ylabel(bpsk调制信号幅值); %作e（bpsk调制信号）的波形图%信号的频谱fs=

18、200; %采样频率n=length(s); %基带信号长度f=0:fs/n:fs-fs/n-fs/2; %修正频率f的范围s=fft(s); %基带信号s的快速傅里叶变换e=fft(e); %基带信号e的快速傅里叶变换c=fft(c); %基带信号c的快速傅里叶变换subplot(322);plot(f,abs(fftshift(s); %基带信号的频谱title(基带信号频谱);xlabel(f/hz);ylabel(s(w); grid on;subplot(324);plot(f,abs(fftshift(c); %载波信号的频谱title(载波信号频谱);xlabel(f/hz);y

19、label(c(w); grid on;subplot(326);plot(f,abs(fftshift(e); %调制信号的频谱title(调制信号频谱);xlabel(f/hz);ylabel(e(w); grid on;%加高斯噪声am=0.7; %输入信号经信道后振幅由1衰减为0.7snr=5; %输入信噪比snr=10(snr/10);n0=(am*am)/2/snr; %计算噪声功率n0_db=10*log10(n0); %将噪声功率转换为dbwni=wgn(1,n,n0_db); % 产生1行n列的高斯噪声yi=e+ni; %bsk已调信号中加入白噪声，输入信噪比为snrfigu

20、re(2);subplot(2,1,1);plot(yi);grid on;xlabel(加入高斯白噪声的已调信号yi(t);%带通滤波器b1,a1 = butter(3,2*pi*0.0001,2*pi*0.01); %计算带通滤波器的h(z)系数y=filter(b1,a1,yi); %对信号yi进行滤波，得到信号yfigure(2);subplot(2,1,2);plot(y);grid on;xlabel(经带通滤波器后信号);%与恢复载波相乘x1=2*c.*y;figure(3);subplot(2,1,1);plot(x1);grid onxlabel(与恢复载波相乘后的信号x1(

21、t);%低通滤波器b2,a2=butter(2,0.005); %计算h(z)系数 ，频率为(1/200)x=filter(b2,a2,x1); %对信号x1滤波，得到信号xfigure(3);subplot(2,1,2);plot(x);grid onaxis(0 n -2 2);xlabel(经低通滤波器后信号波形)%抽样判决x=fun_panjue(x);%调用函数，进行抽样判决figure(4);subplot(2,1,1);plot(x);grid on;xlabel(加噪后解调信号x(t);axis(0 n -2 2);%消除延迟x=fun_yanc(x); %调用函数，进行消除延

22、迟figure(4);subplot(2,1,2);plot(x);grid on;xlabel(加噪后去掉延迟的解调信号x(t);axis(0 n -2 2);%误码率计算err1=length(find(x=s) %计算解调信号中错误码元个数pe_test1=err1/n %计算实际误码率pe1=(1/2)*erfc(sqrt(snr) %计算系统理论误码率% 理论误码率曲线pe=;for snr=1:10 am=0.7; %输入信号经信道后振幅由1衰减为0.7 e=am*am/2; snr=10(snr/10); n0=(am*am)/(2*snr); no=n0/(2*200); %计

23、算噪声功率 n0_db=10*log10(n0);%将噪声功率转换为dbw ni=wgn(1,n,n0_db);% 产生1行n列的高斯噪声 yi=e+ni; %bsk已调信号中加入白噪声，输入信噪比为snr y=filter(b1,a1,yi);%对yi进行滤波(带通滤波器)，得到信号y x1=2*c.*y; %与恢复载波相乘 xx=filter(b2,a2,x1); %经低通滤波器滤波 xx=fun_panjue(xx);%抽样判决 xx=fun_yanc(xx); %消除延迟 snr=10(snr/10); pe=pe,(1/2)*erfc(sqrt(snr); %计算理论误码率endpe

24、;figure;snr=1:10;semilogy(snr,pe,b-);hold on %以log10(pe)为纵坐标画图grid on子程序一：抽样判决（由于存在“倒pi”现象，故以0为基准进行判决，大于0，判为0；小于0，判为1）function w=fun_panjue(w)n=length(w);if w(100)0 w(1:100)=0;else w(1:100)=1;endfor i=101:n if w(i)0; w(i)=0; else w(i)=1; endend子函数二：消除延迟（由于抽样判决后会产生延迟，故认为对解调信号进行时移，以减小甚至消除延迟，便于与基带信号进行比

25、较，计算实际误码率）function m=fun_yanc(m)n=length(m);leng=0;if m(1)=0 for i=1:n if m(i)=1 leng=i; break; end endelse for i=1:n if m(i)=0 leng=i; break; end endendleng1=leng-(floor(leng/200)*200;for i=1:(n-leng1) m(i)=m(i+leng1);endfor i=(n-leng1):n m(i)=m(n-200+10);end4.4仿真波形图：图7 基带信号、载波信号和调制信号时域及频域波形图图8 加入白

26、噪声及带通滤波后的波形图9 图10图11 误码率曲线err1 = 1910pe_test1 = 0.9550pe1 = 0.00605、 总体系能分析 信道噪声对系统性能的主要影响是在接收信号中引入了比特差错。在二进制系统中,比特差错率表现为将符号1误认为0,或将符号0误认为符号1。很明显比特差错的频率越高,接收机的输出信号与原始信息之间的差异就越大。在存在信道噪声的情况下,可以用平均符号差错概率来衡量二进制信息传输的逼真度。平均符号差错概率的定义为,接收机输出的重构符号与所传输的二进制不相同的平均概率。在原始二进制波形中的所有比特均具有相同重要性的条件下,平均符号差错概率又称为误比特率(be

27、r)。但是,在重构原始消息信号的模拟波形时,不同的符号差错可能需要区别对待。例如码字(表示消息信号的量化抽样值) 中重要的比特发生的错误要比不重要的比特发生的错误有害得多。六、设计总结 bpsk，把模拟信号转换成数据值的转换方式之一，利用偏离相位的复数波浪组合来表现信息键控移相方式。bpsk使用了基准的正弦波和相位反转的波浪，使一方为0，另一方为1，从而可以同时传送接受2值得信息。由于最单纯的键控移相方式虽抗噪较强但传送效率差，所以常常使用利用4个相位的qpsk和利用8个相位的8psk。 总体来说这次实验的原理比较困难，但是转换为编程还是费了不少时间和精力的，特别是里边的一些算法，在c语言中要

28、通过若干行代码才能实现，在matlab中一个函数就实现了，我为此还请教了同学，matlab的真的是又实用又方便，具有巨大的魅力，以后还要努力学习，争取连考试过了，然后再好好学习matlab，争取在matlab上取得一些成绩。 通信原理是通信专业的一门主干技术基础课，通过该门课程的学习，使我们掌握确定信号调制解调原理的特性，信号的传输的特性，确定信号的发送与接收的基本方法以及某些典型系统引出的一些重要的基本概念。 本设计要求采用matlab实现对8psk通信系统的蒙特卡罗仿真并且绘制相关的图形，此题比较难做涉及到的方面很多，对8psk及蒙特卡罗仿算法很不了解查阅了很多资料然而涉及的有较少，但查资料当中学到了很多不知道的东西，加深了对psk的了解。 在本次课程设计