QPSK调制解调的仿真

通信工程专业综合课程设计第1页共23页引言近年来，软件无线电作为解决通信体制兼容性问题的重要方法受到各方面的注意。它的中心思想是在通用的硬件平台上，用软件来实现各种功能，包括调制解调类型、数据格式、通信协议等。通过软件的增加、修改或升级就可以实现新的功能，充分体现了体制的灵活性、可扩展性等。其中软件的增加、高频谱效率的调制解调模块是移动通信系统的关键技术，它的软件化也是实现软件无线电的重要环节。QPSK是英文QuadraturePhaseShiftKeying的缩略语简称，意为正交相移键控，是一种数字调制方式。在19世纪80年代初期,人们选用恒定包络数字调制。这类数字调制技术的优点是已调信号具有相对窄的功率谱和对放大设备没有线性要求,不足之处是其频谱利用率低于线性调制技术。19世纪80年代中期以后,四相绝对移相键控(QPSK)技术以其抗干扰性能强、误码性能好、频谱利用率高等优点,广泛应用于数字微波通信系统、数字卫星通信系统、宽带接入、移动通信及有线电视系统之中。通过完成设计内容，复习QPSK调制解调的基本原理，同时也要复习通信系统的主要组成部分，了解调制解调方式中最基础的方法。了解QPSK的实现方法及数学原理。并对“通信”这个概念有个整体的理解，学习数字调制中误码率测试的标准及计算方法。同时还要复习随机信号中时域用自相关函数，频域用功率谱密度来描述平稳随机过程的特性等基础知识，来理解高斯信道中噪声的表示方法，以便在编程中使用。理解QPSK调制解调的基本原理，并使用MATLAB编程实现QPSK信号在高斯信道和瑞利衰落信道下传输，以及该方式的误码率测试。复习MATLAB编程的基础知识和编程的常用算法以及使用MATLAB仿真系统的注意事项，并锻炼自己的编程能力，通过编程完成QPSK调制解调系统的仿真，以及误码率测试，并得出响应波形。在完成要求任务的条件下，尝试优化程序。本课设是基于Matlab的软件仿真，只需PC机上安装MATLAB6.0或者以上版本即可。课设的要求是1.构建一个理想信道基本QPSK仿真系统,要求仿真结果有：基带输入波形及其功率谱;QPSK信号及其功率谱；QPSK信号星座图。2.构建一个在AWGN（高斯白噪声）信道条件下的QPSK仿真系统，得出高斯白噪声信道条件下的误码性能以及高斯白噪声的理论曲线，要求所有误码性能曲线在同一坐标比例下绘制。通过本次实验，除了和队友培养了默契学到了知识之外，还可以将次实验作为一种推广，让更多的学生来深入一层的了解QPSK以至其他调制方式的原理和实现方法。可以方便学生进行测试和对比。足不出户便可以做实验。通信工程专业综合课程设计第2页共23页1方案论证本次课设是基于MATLAB环境下对QPSK调制解调系统进行仿真。在Matlab环境下有两种仿真方案，一种是基于simulink对QPSK进行模块化的仿真，要求是不能直接调用软件里面的集成模块，以此来实现QPSK系统的仿真。另一种方案是在Matlab环境下，用软件编程的方法来实现调制解调系统的仿真，并且得出不同信道的误码率，并作以比较。1.1方案一在Matlab环境下，用软件编程的方法来实现QPSK调制解调系统的仿真。Matlab是一款由美国MathWorks公司出品的商业数学软件。MATLAB是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。除了矩阵运算、绘制函数/数据图像等常用功能外，MATLAB还可以用来创建用户界面及与调用其它语言(包括C，C++和FORTRAN)编写的程序。我们此次可以应用它的这种功能，在熟练掌握语言编程以及深刻理解QPSK的基础上，严格按照程序仿真的流程图对各部分进行程序编写并仿真。应用理论知识来验证最终的仿真结果是否正确。1.2方案二在Matlab环境下的Simulink模块库中应用模块来搭建QPSK调制解调系统的仿真电路模块。在深入了解QPSK调制解调系统理论知识的情况下，根据理论知识构建调制解调的系统框图。在此基础上选择合适的电路器件以及设置合理的参数，构成调制解调的总体电路框图。Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中只要通过简单的鼠标操作，就可以构造出复杂的系统。Simulink提供了一个建立模型方块图的图形用户接口，这个创建过程只需单击和拖动鼠标操作就能完成，它提供了一种更快捷、直接明了的方式，而且用户可以立即看到系统的仿真结果。在此基础上完成电路图的搭建，得出仿真结果，并与实际的理论结果相比较。1.3方案确立在经过详细的考虑和比较之后，我确立了方案一为本次课程设计的最终方案。QPSK的调制方法有两种，分别是相乘电路法和选择法，结合本次设计的实际情况，最终选择用程序方式实现相乘电路的调制方式。因为在Matlab中应用程序对QPSK调制解调系统进行仿真，一方面降低了系统设计的复杂性，并且有效的克服了电子瓶颈的问题。因为在电子电路中，总体电路对参数设置是非常敏感的，一旦参数设置出现小的偏差，将会影响到整个电路的结果。另一方面程序的仿真过程中可以将主程序分为很多个子程序，将它们逐个进行仿真，这样不会影响到整体结果，可以逐步进行调试。结合以上的比较论证以及自身的能力，最终我决定用程序法来完成此次对QPSK调制解调系统的仿真。通信工程专业综合课程设计第3页共23页2仿真原理四相相位调制解调是利用载波的四种不同相位差来表征输入的数字信息,是四进制移相键控。QPSK是在M=4时的调相技术,它规定了四种载波相位,分别为45°,135°,225°,275°,调制器输入的数据是二进制数字序列,为了能和四进制的载波相位配合起来,则需要把二进制数据变换为四进制数据,这就是说需要把二进制数字序列中每两个比特分成一组,共有四种组合,即00,01,10,11,其中每一组称为双比特码元。每一个双比特码元是由两位二进制信息比特组成,它们分别代表四进制四个符号中的一个符号。QPSK中每次调制可传输2个信息比特,这些信息比特是通过载波的四种相位来传递的。解调器根据星座图及接收到的载波信号的相位来判断发送端发送的信息比特。数字调制用“星座图”来描述,星座图中定义了一种调制技术的两个基本参数:(1)信号分布;(2)与调制数字比特之间的映射关系。星座图中规定了星座点与传输比特间的对应关系,这种关系称为“映射”,一种调制技术的特性可由信号分布和映射完全定义,即可由星座图来完全定义[3]。在QPSK调制中,QPSK信号可以看作两个载波正交的2PSK调制器构成。串/并变换器将输入的二进制序列分为速率减半的两个并行的双极性序列,然后分别对sinωct和cosωct调制,相加后得到QPSK调制信号。QPSK同相支路和正交支路可分别采用相干解调方式解调,得到I(t)和Q(t)。经抽样判决和并/串变换器,将上、下支路得到的并行数据恢复成串行数据。2．1QPSK调制原理在QPSK调制中,QPSK信号可以看作两个载波正交的2PSK调制器构成。串/并变换器将输入的二进制序列分为速率减半的两个并行的双极性序列,然后分别对sinωct和cosωct调制,相加后得到QPSK调制信号。QPSK同相支路和正交支路可分别采用相干解调方式解调,得到I(t)和Q(t)。经抽样判决和并/串变换器,将上、下支路得到的并行数据恢复成串行数据。调制原理框图如图2.1所示图2.1QPSK调制原理框图通信工程专业综合课程设计第4页共23页原理分析：基本原理及系统结构QPSK与二进制PSK一样，传输信号包含的信息都存在于相位中。的别的载波相位取四个等间隔值之一，如л/4,3л/4,5л/4,和7л/4。相应的，可将发射信号定义为2/cos[2(21)/4]Etfti0≤t≤TSi（t）＝0（2.1.1）其中，i＝1，2，2，4；E为发射信号的每个符号的能量，T为符号持续时间，载波频率f等于nc/T，nc为固定整数。每一个可能的相位值对应于一个特定的二位组。例如，可用前述的一组相位值来表示格雷码的一组二位组：10，00，01，11。下面介绍QPSK信号的产生和检测。如果a为典型的QPSK发射机框图。输入的二进制数据序列首先被不归零（NRZ）电平编码转换器转换为极性形式，即负号1和0分别用bE和－bE表示。接着，该二进制波形被分接器分成两个分别由输入序列的奇数位偶数位组成的彼此独立的二进制波形，这两个二进制波形分别用a1（t），和a2（t）表示。容易注意到，在任何一信号时间间隔内a1（t），和a2（t）的幅度恰好分别等于Si1和Si2，即由发送的二位组决定。这两个二进制波形a1（t），和a2（t）被用来调制一对正交载波或者说正交基本函数：1（t）＝2cos(2)cftT，2（t）＝2sin(2)cftT。这样就得到一对二进制PSK信号。1（t）和2（t）的正交性使这两个信号可以被独立地检测。最后，将这两个二进制PSK信号相加，从而得期望的QPSK。2.2QPSK解调原理在QPSK解调中,正交支路和同相支路分别设置两个相关器(或匹配滤波器),得到I(t)和Q(t),经电平判决和并/串变换后即可恢复原始信息从发射机发射的已调信号经过传输媒质传播到接收端,接收机接收到的己调信号为:SQPSK(t)=I(t)cosωct+Q(t)sinωct(2.2.1)I(t)、Q(t)分别为同相和正交支路,ωc为载波频率,那么相干解调后,同相支路相乘可得:Ii(t)=SQPSK(t)cosωct=[I(t)cosωct+Q(t)sinωct]cosωct=I(t)cos2wct+Q(t)sinwct/2=I(t)/2-I(t)cos2ωct+Q(t)sin2ωct（2.2.2）正交支路相乘可得:Qq(t)=SQPSK(t)sinωct=[I(t)cosωct+Q(t)sinωct]sinωct=I(t)sinωct×cosωct+Q(t)sin2ωct=I(t)sin2ωct/2+Q(t)-Q(t)cos2ωct（2.2.3）经过低通滤波器可得:Ii(t)=0.5I(t),Qq(t)=0.5Q(t)（2.2.4）通信工程专业综合课程设计第5页共23页原理框图如图2.2所示：1（t）同相信道门限＝02（t）正交信道门限＝0图2.2QPSK解调原理框图原理分析：QPSK接收机由一对共输入地相关器组成。这两个相关器分别提供本地产生地相干参考信号1（t）和2（t）。相关器接收信号x（t），相关器输出地x1和x2被用来与门限值0进行比较。如果x10,则判决同相信道地输出为符号1；如果x10,则判决同相信道的输出为符号0。；类似地。如果正交通道也是如此判决输出。最后同相信道和正交信道输出这两个二进制数据序列被复加器合并，重新得到原始的二进制序列。在AWGN信道中，判决结果具有最小的负号差错概率。2.3调制解调系统总框图串/并变换脉冲成型升余弦滤波相干载波coswt相干载波sinwt二进制双极性不归零序列××+AWGN/RAY信道××相干解调coswt相干解调sinwt过采样过采样脉冲成型升余弦滤波低通滤波低通滤波并/串变换抽样判决抽样判决输出图2.3QPSK调制与解调系统总框图判决门限低通filrer判决门限复接器接收信号x（t）低通filrer发送二进制序列的估计通信工程专业综合课程设计第6页共23页3仿真与分析3.1程序仿真流程图程序仿真流程图是程序仿真的步骤及大体框架，基于流程图，我们可以有效的避免程序次序的先后问题，将主程序分割成很多个子程序来逐一实现，严格按照流程图的次序进行编程仿真。最终将各个子程序综合起来实现最终的程序仿真。流程图如图3.1所示开始产生码元temp=rand;if(temp0.25),00参数初始化numofsymbolerror=0;numofbiterror=0;num=0；N=50000;numofbiterror100num=num+1;Yi=0iNYi++Nps=numofsymbolerror/num;pb=numofbiterror