基于Matlab的QPSK调制解调仿真设计与研究.docx

1、基于Matlab的QPSK调制解调仿真设计与研究天津理工大学计算机与通信工程学院通信工程专业设计说明书基于Matlab/Simulink 的QPSK调制解调仿真设计与研究 姓 名 韩双年 学 号 20092226 班 级 09通信3班 指导老师 白媛 日 期 2012-12-16 摘 要正交相移键控（QPSK），是一种数字调制方式。QPSK技术具有抗干扰能力好、误码率低、频谱利用效率高等一系列优点。论文主要介绍了正交相移键控(QPSK)的概况，以及正交相移键控QPSK的调制解调概念和原理，利用Matlab中M文件和Simulink模块对QPSK的调制解调系统进行了仿真，对QPSK在高斯白噪声信

2、道中的性能进行了，分析了解Simulink中涉及到QPSK的各种模块的功能。【关键词】Matlab QPSK Simulnk 仿真第一章 前 言1.1 专业设计任务及要求1了解并掌握QPSK调制与解调的基本原理；2在通信原理课程的基础上设计与分析简单的通信系统；3学会利用MATLAB7.0编写程序进行仿真，根据实验结果能分析所设计系统的性能。4学习MATLAB的基本知识，熟悉MATLAB集成环境下的Simulink的仿真平台。5利用通信原理相关知识在仿真平台中设计QPSK调制与解调仿真系统并用示波器观察解调后的波形6在指导老师的指导下，独立完成课程设计的全部内容，能正确的阐述和分析设计和实验结

3、果。1.2 Matlab简介 MATLAB是MATrix LABoratory的缩写，是一款由美国Math Works公司出品的商业数学软件。MATLAB 是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。除了矩阵运算、绘制函数/数据图像等常用功能外，MATLAB还可以用来创建用户界面及与调用其它语言(包括C，C+和FORTRAN)编写的程序。尽管MATLAB主要用于数值计算,但是因为大量的额外的工具箱它也适合于不同领域的应用,如控制系统设计与分析、图像处理和信号处理和通信、金融建模和分析等。除了一个完整的Simulink包,提供了一个可视化的开发环境,通常用

4、于系统仿真、动态/嵌入式系统开发等。1.3 Matlab下的simulink简介Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中只要通过简单的鼠标操作，就可以构造出复杂的系统。Simulink提供了一个建立模型方块图的图形用户接口，这个创建过程只需单击和拖动鼠标操作就能完成，它提供了一种更快捷、直接明了的方式，而且用户可以立即看到系统的仿真结果。Simulink具有适应面广、结构和流程清晰及仿真精细、效率高、贴近实际、等优点，基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件

5、应用于Simulink。1.4 通信系统模型通信系统就是传递信息所需要的一切技术设备和传输媒质的总和，包括信息源、发送设备、信道、接收设备和信宿(受信者) ，它的一般模型如图1.4.1所示。 图1.4.1 通信系统一般模型模拟通信系统是利用模拟信号来传递消息的通信系统，其模型如图1.4.2所示。 图1.4.2 模拟通信系统模型第二章 QPSK调制2.1 QPSK介绍 Quadrature Phase Shift Keying通过使用载波的四个各不相同的相位差来表示输入的信息，是具有四进制的相移键控。QPSK是在M=4时的数字的调相技术，它通过约定的四种载波相位，分别为45，135，225，27

6、5，输入数据为二进制的数字序列，因为载波相位是四进制的，所有我们需要把二进制的数据变为四进制的，即把二进制序列中每两个比特分成一组，四种排列组合，即00，01，10，11，双比特码元即为一组。每两位二进制信息比特构成每一组，它们分别表示着着四个符号中的某一个符号。2.2 QPSK调制原理QPSK的调制有两种产生方法相乘电路法和选择法。2.2.1 相乘法输入信号是二进制不归零的双极性码元，它通过“串并变换”电路变成了两路码元。变成并行码元后，每个码元的持续时间是输入码元的两倍。用两路正交载波去调制并行码元。图2.2.1选择法QPSK的调制中，QPSK信号可以看成是两个载波正交的2PSK信号调制器

7、构成。原理分析如下：基本原理和系统结构QPSK与二进制PSK一样，传输信号包含的信息都存在于相位中。个别的载波相位取四个等间隔值之一，如/4、3/4、5/4、7/4。相应的，可将发射信号定义为：其中，i1，2，3，4；E是发射信号的每个符号的能量，T为符号的持续时间，载波频率f等于nc/T，nc为固定整数。每一个可能的相位值对应于一个特定的二位组。下面介绍QPSK信号的产生和检测。如图为典型的QPSK发射机框图。输入的二进制数据序列首先被不归零(NRZ)电平编码转换器转换为极性形式，即负号1和0分别用和-表示。该二进制波形被分接器分成两个分别由输入序列的奇数位偶数位组成的彼此独立的二进制波形，

8、这两个二进制波形分别用a1(t)和a2(t)表示。此时，在任何一信号时间间隔内a1(t)，和a2(t)的幅度恰好分别等于Si1和 Si2，即由发送的二位组决定。这两个二进制波形a1(t)和a2(t)被用来调制一对正交载波：,。这样就得到一对二进制PSK信号。和的正交性使这两个信号可以被独立地检测。最后，将这两个二进制PSK信号相加，从而得期望的QPSK。2.2.2 选择法输入基带信号经过串并变换后用于控制一个相位选择电路，按照当时的输入双比特ab，决定选择哪个相位的载波输出。图2.2.2选择法2.3 QPSK调制原理框图 图2.3 调制原理框图2.4 QPSK调制方式的Matlab仿真I 路信

9、号是用余弦载波，由2进制数据流的奇数序列组成；Q路信号用正弦载波，由2进制数据流的偶数序列组成。下面的a是Idata，b就是Qdata，它们分布与各自的载波相乘分别输出 I 路信号和 Q 路信号。I 路信号加上Q路信号就是QPSK输出信号。当 I 路载波信号是0相位时为1，是180相位时为0；当Q路载波信号是0相位时为1，是180相位时为0。2.4 matlab调制仿真图2.5 QPSK调制方式Matlab-simulink仿真2.5.1 simulink调制建模 图2.5.1调制框图（1)产生需要的信号源在搭建QPSK调制解调系统中使用伯努力信号发生器产生随机的01比特序列，每两比特代表就一

10、个符号。Bernoulli Binary Generator模块利用伯努利分布的原理，相应得到参数为p的伯努利分布。伯努利分布的均值1 - p和方差p(1 p)的。一个零概率参数指定p。本次实验中的p设置为0.5，即0和1等概。采样时间可根据需要进行设置，例如测误码率时采样时间设为0.01s。图2.5.2信号源参数设置（2）串并变换我们先通过使用buffer 这个模块来实现将信号源信号转变为两路信号。Buffer 模块可以重新分配缓冲区块的输入样本，用到了Demux，可以将一个复合输入转化为多个单一输出，即可以输出多个采样率较低的帧信号。但会产生与缓冲区容量相同的时延。所以，我们可以设置buf

11、fer的参数容量为2。图2.5.3 Buffer的参数设置（3）单极性信号转化为双极性信号因为QPSK的调制信号要求的是双极性信号，所以用伯努利随机生成二进制Generator模块产生的信号必须经过转化才能够被使用。利用加法模块和常数产生模块将1和0的序列各自减去1/2，再利用比例运算模块乘以2，就得到了1和-1 的双极性序列。（4）调制模块分别将两路信号乘以相位相差 的载波，然后相加。载波由正弦信号发生器产生。正弦波模块的参数设置为可基于时间的模式，时间设为使用仿真时间，我们设载波信号的幅度为1，载波频率可根据需要来进行设置，两路载波同频正交，相位相差 / 2。我们设上支路的相位为0，下支路

12、的相位为 。图2.5.4上支路载波参数图2.5.5下支路载波参数2.5.2 simulink调制仿真结果图2.5.6信号源和转变后的双极性信号图2.5.7上支路载波图2.5.8 下支路载波图2.5.9调制信号第三章 QPSK解调3.1QPSK解调原理QPSK接收机由一对共输入地相关器组成。这两个相关器分别提供本地产生地相干参考信号和。相关器接收信号x(t),相关器输出地x1和x2被用来与门限值0进行比较。如果x10,则判决同相信道地输出为符号1；如果x1=-1 & t(i)=5& t(i)=0 & t1(i)=4& t1(i)0 %20点为同一个值，20点数据叠加后与阈值0比较 data_re

13、cover_a(i:i+19)=1; %data_recover_a是并/串转换后的20000点 bit_recover=bit_recover 1; %bit_recover是1000点数据 else data_recover_a(i:i+19)=-1; bit_recover=bit_recover -1; endenderror=0;dd = -2*bit_in+1; % 将bit_in中的1变成-1，0变成1ddd=dd; %ddd是1表示0，-1表示1的原始序列，1000个点ddd1=repmat(ddd,20,1); %ddd1是20*1000的矩阵for i=1:2e4 ddd2

14、(i)=ddd1(i); %将ddd1拉直为1*20000的行向量ddd2endfor i=1:1e3 if bit_recover(i)=ddd(i) error=error+1; endendp=error/1000;figure(1)subplot(4,1,1);plot(t2,ddd2);axis(0 100 -2 2);title(原序列);subplot(4,1,2);plot(t1,I_demo);axis(0 100 -2 2);title(I 支路解调);subplot(4,1,3);plot(t1,Q_demo);axis(0 100 -2 2);title(Q 支路解调);subplot(4,1,4);plot(t2,data_recover_a);axis(0 100 -2 2);title(解调后序列);