移动通信中的数字调制技术与毕业论文

1、毕业设计说明书移动通信中的数字调制技术与仿真实现信息与通信工程学院学生姓名： 学号： 通信工程学 院： 专 业： 指导教师： 2012年 6 月移动通信中的数字调制技术与仿真实现摘要调制技术是任何频带通信中最重要的一项技术。现代移动通信系统都使用数字调制技术，数字调制技术的改进也是通信系统性能提高的重要途径。本文首先研究了几种基本的数字调制方式和广泛应用在现代移动通信中的新型数字调制技术。然后，运用仿真软件MATLAB对其进行编程仿真，通过仿真图形观察了调制解调过程中各环节的时域波形，并结合数字调制技术的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。最后，在仿真的基础上分析比

2、较了各种调制技术的性能。关键词：移动通信；数字调制；分析与仿真；MATLAB移动通信中的数字调制技术与仿真实现AbstractModulation of any frequency band communications is the most important technology. Modern mobile communication systems use digital modulation, using digital modulation techniques to improve communication system performance is also an impo

3、rtant way to improve. This paper studies several basic digital modulation methods and which widely used in modern mobile communication in the new digital modulation techniques. Then, using simulation software MATLAB be programmed simulation to observe the modulation and demodulation process of drawi

4、ng each part of the time domain waveforms, and modulation combined with digital modulation principle, track and analyze the performance of each part of the impact on the modulation and simulation reliability of the model. Finally, analysising and comparing the performance of various modulation techn

5、iques based on the simulation.Key words: mobile communications; digital modulation; analysis and simulation; MATLAB目录 TOC o 1-3 h z u HYPERLINK l _Toc327097824 1 绪论 PAGEREF _Toc327097824 h 1 HYPERLINK l _Toc327097825 1.1 移动通信技术概述 PAGEREF _Toc327097825 h 1 HYPERLINK l _Toc327097826 1.2 调制技术 PAGEREF _

6、Toc327097826 h 1 HYPERLINK l _Toc327097827 1.2.1 调制技术的概念 PAGEREF _Toc327097827 h 1 HYPERLINK l _Toc327097828 1.2.2 调制技术的分类 PAGEREF _Toc327097828 h 1 HYPERLINK l _Toc327097829 1.3 数字调制的意义 PAGEREF _Toc327097829 h 2 HYPERLINK l _Toc327097830 1.4 MATLAB在移动通信系统仿真中的应用 PAGEREF _Toc327097830 h 2 HYPERLINK l

7、 _Toc327097831 2 基本数字调制系统的原理 PAGEREF _Toc327097831 h 4 HYPERLINK l _Toc327097832 2.1 二进制数字调制的原理 PAGEREF _Toc327097832 h 4 HYPERLINK l _Toc327097833 2.1.1 二进制幅度键控（2ASK） PAGEREF _Toc327097833 h 4 HYPERLINK l _Toc327097834 2.1.2 二进制频移键控（2FSK） PAGEREF _Toc327097834 h 5 HYPERLINK l _Toc327097835 2.1.3 二进

8、制相移键控（2PSK） PAGEREF _Toc327097835 h 5 HYPERLINK l _Toc327097836 2.1.4 二进制差分相移键控（2DPSK） PAGEREF _Toc327097836 h 5 HYPERLINK l _Toc327097837 2.2 多进制数字调制 PAGEREF _Toc327097837 h 6 HYPERLINK l _Toc327097838 2.3 二进制数字调制方式的性能比较 PAGEREF _Toc327097838 h 6 HYPERLINK l _Toc327097839 3 新型调制系统的原理 PAGEREF _Toc32

9、7097839 h 8 HYPERLINK l _Toc327097840 3.1 最小频移键控（MSK） PAGEREF _Toc327097840 h 8 HYPERLINK l _Toc327097841 3.1.1 MSK信号的基本原理 PAGEREF _Toc327097841 h 8 HYPERLINK l _Toc327097842 3.1.2 MSK信号的产生 PAGEREF _Toc327097842 h 10 HYPERLINK l _Toc327097843 3.1.3 MSK信号的解调 PAGEREF _Toc327097843 h 11 HYPERLINK l _To

10、c327097844 3.2 高斯滤波最小频移键控（GMSK） PAGEREF _Toc327097844 h 11 HYPERLINK l _Toc327097845 3.2.1 GMSK调制原理 PAGEREF _Toc327097845 h 12 HYPERLINK l _Toc327097846 3.3 四相相移键控（QPSK） PAGEREF _Toc327097846 h 15 HYPERLINK l _Toc327097847 3.3.1 QPSK的基本原理 PAGEREF _Toc327097847 h 15 HYPERLINK l _Toc327097848 3.3.2 QP

11、SK的调制原理 PAGEREF _Toc327097848 h 17 HYPERLINK l _Toc327097849 QPSK解调原理 PAGEREF _Toc327097849 h 18 HYPERLINK l _Toc327097850 3.4 交错正交相移键控（OQPSK） PAGEREF _Toc327097850 h 19 HYPERLINK l _Toc327097851 3.4.1 OQPSK基本原理 PAGEREF _Toc327097851 h 19 HYPERLINK l _Toc327097852 3.4.2 OQPSK的调制原理 PAGEREF _Toc327097

12、852 h 20 HYPERLINK l _Toc327097853 3.5 正交频分复用（OFDM） PAGEREF _Toc327097853 h 20 HYPERLINK l _Toc327097854 3.5.1 OFDM概述 PAGEREF _Toc327097854 h 20 HYPERLINK l _Toc327097855 OFDM的基本原理 PAGEREF _Toc327097855 h 21 HYPERLINK l _Toc327097856 3.5.3 OFDM的实现 PAGEREF _Toc327097856 h 23 HYPERLINK l _Toc327097857

13、 3.6 数字调制技术的应用 PAGEREF _Toc327097857 h 25 HYPERLINK l _Toc327097858 4 数字调制系统的仿真和结果分析 PAGEREF _Toc327097858 h 27 HYPERLINK l _Toc327097859 4.1 数字调制系统的仿真分析 PAGEREF _Toc327097859 h 27 HYPERLINK l _Toc327097860 4.2 MSK信号仿真 PAGEREF _Toc327097860 h 28 HYPERLINK l _Toc327097861 MSK信号仿真思路 PAGEREF _Toc327097

14、861 h 28 HYPERLINK l _Toc327097862 4.2.2 MSK信号的仿真结果分析 PAGEREF _Toc327097862 h 30 HYPERLINK l _Toc327097863 QPSK信号仿真 PAGEREF _Toc327097863 h 31 HYPERLINK l _Toc327097864 4.3.1 QPSK信号的仿真思路 PAGEREF _Toc327097864 h 31 HYPERLINK l _Toc327097865 QPSK信号仿真结果分析 PAGEREF _Toc327097865 h 33 HYPERLINK l _Toc3270

15、97866 5 结 论 PAGEREF _Toc327097866 h 34 HYPERLINK l _Toc327097867 附录1 PAGEREF _Toc327097867 h 35 HYPERLINK l _Toc327097868 附录2 PAGEREF _Toc327097868 h 44 HYPERLINK l _Toc327097869 参考文献 PAGEREF _Toc327097869 h 52 HYPERLINK l _Toc327097870 致 谢 PAGEREF _Toc327097870 h 541 绪论 移动通信技术概述移动通信技术是在20世纪80年代开始发展

16、起来的，移动通信技术的发展速度远远超过固定网络技术，普及范围相当广泛。ITU预计2010年全球移动蜂窝用户数量将达到50亿，人们对移动通信的需求推动了移动通信技术的快速发展，至今，移动通信已经历了20世纪80年代的第一代模拟技术（1G）和90年代的第二代窄带数字技术（2G）这两个发展阶段。近些年来，随着无线通信宽带化技术的突破，移动通信正在向以CDMA为基础，以宽带化通信为特征的第三代3G技术发展，伴随着第三代移动通信技术（3G）逐步实施，移动通信未来的发展及演进问题成了研究热点，因此第四代移动通信技术（4G）被提出。移动通信经历了1G和2G，完成了从模拟技术向数字技术的过渡，现正在向3G过渡

17、和走向更远的4G，把移动通信从窄带推向宽。1.2 调制技术1.2.1 调制技术的概念在移动通信中，信源产生的原始信号绝大部分需要经过调制，变换为适合于在信道内传输的信号，才能在线路中传输。把输入信号变换为适合于通过信道传输的波形，这一变换过程成为调制。通常把原始信号称为调制信号，也称基带信号；被调制的高频周期性脉冲起运载原始信号的作用，因此称载波1。调制技术其实也就是实现了信源与信道的频带匹配。 调制技术主要有一下3个方面的功能：1.频率变换：为了采用无线传送方式，如将（0.3MHz3.4kHz）有效带宽内的信号调制到高频段上去。2.实现信道复用：例如将多路型号互不干扰地安排在同一物理信道中传

18、输。3.提高抗干扰性：抗干扰性（即可靠性）与有效性互相制约，通常可通过牺牲有效性来提高抗干扰性，如FM替代AM。1.2.2 调制技术的分类数字信号调制的基本类型分为振幅键控（ASK）、频移键控（FSK）、和相移键控（PSK）此外许多基本 调制类型改进或综合而获得的新型调制技术。在实际应用中，有两类用得最多的数字调制技术方式11： (1) 线性调制技术，主要包括：移相键控（PSK）,四相移键控（QPSK）,差分四相移相键控（DQPSK）,偏置键控-四相移相键控（OK-QPSK）,4-DQPSK和多电平（移相键控）PSK等。 (2) 恒定包络(连续相位)调制技术,主要包括：最小移频键控（MSK）、

19、高斯滤波最小移频键控（GMSK）、高斯滤波的移频键控（GFSK）和平滑调频（TFM）等。1.3 数字调制的意义现代移动通信系统都使用数字调制技术。现有的通信系统都在由模拟方式向数字方式过渡，数字通信具有很多模拟通信不可比拟的优势，数字通信技术采用数字技术进行加密和差错控制，便于集成。因此这里我们重点讨论数字调制技术。 数字调制是指用数字基带信号对载波的某些参量进行控制，使载波的这些参量随基带信号的变化而变化。根据控制的载波参量的不同，数字调制有调幅、调相和调频三种基本形式，并可以派生出多种其他形式。由于传输失真、传输损耗以及保证带内特性的原因，基带信号不适合在各种信道上进行长距离传输。为了进行

20、长途传输，必须对数字信号进行载波调制，将信号频谱搬移到高频处才能在信道中传输。因此，大部分现代通信系统都使用数字调制技术1、3。另外，由于数字通信具有建网灵活，容易采用数字差错控制技术和数字加密，便于集成化，并能够进入综合业务数字网（ISDN网），所以通信系统都有由模拟方式向数字方式过渡的趋势。因此，对数字通信系统的分析与研究越来越重要，数字调制作为数字通信系统的重要部分之一，对它的研究也是有必要的。通过对调制系统的仿真，我们可以更加直观的了解数字调制系统的性能及影响性能的因素，从而便于改进系统,获得更佳的传输性能。 MATLAB在移动通信系统仿真中的应用随着通信系统复杂性的增加，传统手工分析

21、与电路板试验等分析设计方法已经不能适应发展的需要，通信系统计算机模拟技术日益显示出其巨大的优越性。计算机仿真是根据被研究的真实系统的模型。利用计算机进行实验研究的一种方法。它具有利用模型进行仿真的一系列的模型优点如费用低，易于进行真实系统难于实现各种试验，以及易于实现完全相同条件下重复试验等。Matlab仿真软件就是分析通信系统常用的工具之一。MATLAB是一种交互式的，以矩阵为基础的软件开发环境,它用于科学和工程的计算与可视化。其编程功能简单,并且很容易扩展和有创造新的命令与函数。应用可方便地解决复杂数值计算问题。具有强大的Simulink动态仿真环境,可以实现可视化建模和多工作环境间文件互

22、用和数据交换。Simulink支持连续、离散及两者混合的线性和非线性系统,也支持多种采样速率的多速率系统;Simulink为用户提供了用方框图进行建模的图形接口,它与传统的仿真软件包用差分方程和微分方程建模相比,更直观、方便和灵活。用户可以在和Simulink两种环境下对自己的模型进行仿真、分析和修改。用于实现通信仿真的通信工具包(Communication toolbox,也叫Commlib,通信工具箱)是语言中的一个科学性工具包,提供通信领域中计算、研究模拟发展、系统设计和分析的功能,可以在环境下独立使用,也可以配合Simulink使用。另外，MATLAB的图形界面功能GUI（Graphi

23、cal User Interface）能为仿真系统生成一个人机交互界面，便于仿真系统的操作。因此，在通信系统仿真中得到了广泛应用，本文也选用该工具对数字调制系统进行仿真7。2 基本数字调制系统的原理 调制的目的的是把要传输的模拟信号或数字信号变换成适合信道传输的高频信号。该信号成为已调信号。调制过程用于通信系统的发端。在接收端将已调信号，该过程称为解调。移动通信是数字移动通信，其中的关键技术之一是数字调制技术。对数字调制的主要要求是已调信号的频谱窄和带外衰减快（即所占频带窄，或者说频谱利用率高）；易于采用相干或非相干解调；抗噪声和抗干扰的能力强；以及适宜在衰落信道中的传输。 数字调制的目的是使

24、所传送的信息能够很好地适应于信道特性，以达到最有效、最可靠的传播。在移动通信中，由于颠簸传输的条件极其恶劣，是接收信号幅度发生很大的变化，衰减幅度达到最小。因此，在移动通信中必须采用抗干扰能力强的调制方式。调频制在抗干扰和抗衰落性能上优于调幅制，但调频制也存在着固有的缺点，需要占用较大的带宽，同时还存在着门限效应。当然要同时实现这些最佳的特性是不可能的，因为每种特性都有其局限性，且互相之间会有矛盾。例如，要获得较高的带宽效率，可选用多电平调制，但它要求线性放大，因此会使功率效率降低，而且已调波的包络变化大。如果采用恒包络调制，因要求非线性放大，所以它具有高的功率效率，但又会引起较大的带外辐射。

25、因此，只能折中考虑上述要求。总之，采用调制技术的最终目的，就是使调制以后的调制信号对于干扰有较强的抵抗作用，同时对相邻的信道信号干扰较小，解调方便且易于集成。数字调制可以分为二进制调制和多进制调制，多进制调制是二进制调制的推广。近年来随着移动通信的快速发展，调制技术也随之快速发展，基础的调制技术已不能满足现代移动通信的要求，因此在原有调制技术的基础之上发展而来的调制技术有MSK、GMSK、QPSK、OQPSK、OFDM等。下面介绍二进制调制方式以及新型的调制方式。2.1 二进制数字调制的原理2.1.1 二进制幅度键控（2ASK）振幅键控是利用载波的幅度变化来传递数字信息，而其频率和初始相位保持

26、不变。载波在数字信号1或0的控制下通或断，在信号为1的状态载波接通，此时传输信道上有载波出现；在信号为0的状态下，载波被关断，此时传输信道上无载波传送。那么在接收端我们就可以根据载波的有无还原出数字信号的1和0。2ASK信号功率谱密度的特点如下5：（1）由连续谱和离散谱两部分构成；连续谱由传号的波形g(t)经线性调制后决定，离散谱由载波分量决定；（2）已调信号的带宽是基带脉冲波形带宽的二倍。2.1.2 二进制频移键控（2FSK） 频移键控是利用两个不同频率f1和f2的振荡源来代表信号1和0，用数字信号的1和0去控制两个独立的振荡源交替输出。在2FSK中，载波的频率随二进制基带信号在f1和f2两

27、个频率点间变化。对二进制的频移键控调制方式，其有效带宽为B=2xF+2Fb,xF是二进制基带信号的带宽也是FSK信号的最大频偏，由于数字信号的带宽即Fb值大，所以二进制频移键控的信号带宽B较大，频带利用率小。2FSK功率谱密度的特点如下5： (1) 2FSK信号的功率谱由连续谱和离散谱两部分构成,离散谱出现在和位置；(2) 功率谱密度中的连续谱部分一般出现双峰。若两个载频之差，则出现单峰。2.1.3 二进制相移键控（2PSK）在相移键控中，载波相位受数字基带信号的控制，如在二进制基带信号中为0时，载波相位为0或，为1时载波相位为或0，从而达到调制的目的。2PSK信号的功率密度有如下特点1、5：

28、 (1) 由连续谱与离散谱两部分组成；(2) 带宽是绝对脉冲序列的二倍； (3) 与2ASK功率谱的区别是当P1/2时，2PSK无离散谱，而2ASK存在离散谱。 2.1.4 二进制差分相移键控（2DPSK）前面讨论的2PSK信号中，相位是以未调载波的相位作为参考基准的。由于它利用载波相位的绝对数值表示数字信息，所以又称为绝对相移。2PSK在进行相干解调时，由于载波恢复中相位有0、模糊性，导致解调过程中出现“反向工作”现象，恢复出的数字信号“1”和“0”倒置，从而使2PSK难以实际应用。为了克服此缺点，提出了二进制差分数字相移键控（2DPSK）方式1。2.2 多进制数字调制上面所讨论的都是在二进

29、制数字基带信号的情况，在实际应用中，我们常常用一种称为多进制（如4进制，8进制，16进制等）的基带信号。多进制数字调制载波参数有M种不同的取值，多进制数字调制比二进制数字调制有两个突出的优点：一是由于多进制数字信号含有更多的信息使频带利用率更高；二是在相同的信息速率下持续时间长，可以提高码元的能量，从而减小由于信道特性引起的码间干扰。现实中用得最多的一种调制方式是多进制相移键控（MPSK）。多进制相移键控又称为多相制，因为基带信号有M种不同的状态，所以它的载波相位有M种不同的取值，这些取值一般为等间隔。多进制相移键控有绝对移相和相对移相两种，实际中大多采用四相绝对移相键控（4PSK也称QPSK

30、），四相制的相位有0、/2、3/2四种，分别对应四种状态11、01、00、10。2.3 二进制数字调制方式的性能比较2ASK和2PAK所需要的带宽是码元速率的2倍；2FSK所需的带宽比2ASK和2PAK都要高。各种二进制数字调制系统的误码率取决于解调器输入信噪比r。在抗加性高斯白噪声方面，相干2PSK性能最好，2FSK次之,2ASK最差。ASK是一种应用最早的基本调制方式。其优点是设备简单，频带利用率较高；缺点是抗噪声性能差，并且对信道特性变化敏感，不易是抽样判决器工作在最佳判决门限状态。FSK是数字通信中不可或缺的一种调制方式。其优点是抗干扰能力较强,不受信道参数变化的影响，因此FSK特别适

31、合应用于衰落信道；缺点是占用频带较宽，尤其是MFSK，频带利用率较低。目前，调频体制主要应用于中，低速数据传输中。PSK和DPSK是一种高传输效率的调制方式，其抗噪声能力比ASK和FSK都强，且不易受信道特性变化的影响，因此在高、中速数据传输中得到了广泛的应用。绝对相移（PSK）在相干解调时存在载波相位模糊度的问题，在实际中很少采用于直接传输，MDPSK应用更为广泛5。和ASK、FSK、PSK、和DPSK对应，分别有MASK、MFSK、MPSK和MDPSK。这些多进制数字键控的一个码元中包括更多的信息量。但是，为了得到相同的误比特率，它们需要使用更大的功率或占用更宽的频率。3 新型调制系统的原

32、理 最小频移键控（MSK） MSK信号的基本原理MSK波形有多种表示形式。下面是其中一种 St=Acos2ft+at (3.1)f为载波频率，A为振幅，信号的功率E与A2成正比，相位a(t)携带了所有的信息，其中 at=a0+20.5dqt-nT，KTtk+1T (3.2)a0为初始相位，我们认为它是已知的。0.5为调制指数，它决定了一个符号带来的相位变化, ,qt为相位平滑函数，它很大程度上决定了信号的形状继而影响到性能5。给定输入序列at，MSK的相位轨迹如图3.1所示。各种可能的输入序列所对应的所有可能路径如图3.2所示11。图3.1 MSK的相位序列 图3.2 MSK可能的相位轨迹qt

33、为一个分段函数:当t0时，qt=0，当tLT时，qt=0.5。其中L可以被看作调制方法的记忆长度，它决定了每一个符号究竟影响到该符号以后的多少个符号间隔。实际上MSK属于连续相位调制(CPM)的一种，在CPM中L=1时被称作全响应调制，当L时被称作部分响应调制。MSK属于全响应调制3，即L=1。从MSK的表达式可以得知，MSK的相位是由两部分组成的，一部分是载波随时间连续增加的相位2ft，另外一部分是携带信息的附加相位，它与原始数据息息相关，可以被称为基带相位。一般移频键控(2FSK)两个信号波形具有以下的相关系数8: R=sin2Fk-FtT2Fk-FtT+sin4FcTsin4FcT (3

34、.3)其中Fc=Fk+Ft2 (3.4)因为MSK是一种正交调制，其信号波形的相关系数等于零，所以上式等号右侧的第一项和第二项均应为零。第一项等于零的条件是sin2Fk-FtT=k,（K=1，2，3）令k等于其最小值1，则得到 Fk-Ft=12T=12Fb (3.5)即传号频率和空号频率在一个符号周期内的相位累计严格的相差180。式(3.3)中等号右侧第二项等于零的条件是4FcT=nFc=n4Fb。综上所述得到的频率约束关系: Fc=n4Fb，n4；Fk=n+14Fb，n4； Ft=n4Fb，n4 (3.6)在一个符号周期内必须包含四分之一载波周期的整数倍。载波频率和传号频率相差四分之一符号速

35、率，与空号频率也相差四分之一的符号速率:Fk-Ft=14Fb；Fc-Ft=14Fb；Fk-Ft=14Fb (3.7)从(3.2 )式可以看出，在一个码元周期内，其基带相位总是线性累积2，因此码元终止位与起始相位之差也是2。如果一个码元是“1”那么在该码元周期内，基带相位均匀增加2，在码元末尾处基带相位比码元开始处基带相位要大2。相反如果一个码元是“0”，那么在该码元周期内，基带相位均匀减小2，即在码元末尾处基带相位比码元开始处基带相位要小2，这是MSK的一个重要特征，也是差分解调的依据。所谓“连续”是指当前所要讨论的码元。ak范围kTb-kTb+1内，其起始相位等于与ak相邻的前一个码元的终一

36、止相位(对应于t=kTb时的相位)。对于任何一个码元来说，它在一个码元间隔内，相对于载波相位差虽然只变化2，但在这个码元内，相对于载波相位的实际数值却是千变万化的，这与它前面己经发送过的码流有关。相对于载波相位来说，由式(3.2)可知基带相位值与时间t之间存在着一定的关系。at又称为附加相位函数，它是MSK信号的总相位减去随时间线性增长的载波相位而得到的剩余相位。At的尾的基带相位。此外，随着k值的不同，dk是取值为1的随机数，所以dkTb也是分段线表达式（3.2 )是一直线方程式，直线的斜率是dkTb，截距是上个码元末的相位函数以码元宽度Tb为段。在任意码元期间，此函数的变化量总是2。当dk

37、=1，增大2；当dk=-1时，减少2。 MSK信号的产生MSK是一种在无线移动通信中很有吸引力的数字调制方式，是由2FSK信号的改进而来，因为它有以下两种主要的特点:（1）信号能量的99.5%被限制在数据传输速率的倍的带宽内。谱密度随频率(远离信号带宽中心)倒数的四次幂而下降，而通常的离散相位FSK信号的谱密度却随频率倒数的平方下降。因此，MSK信号在带外产生的干扰非常小。这正是限带工作情况下所希望有的宝贵特点。（2）信号包络是恒定的，系统可以使用廉价高效的非线性器件。从相位路径的角度来看，MSK属于线性连续相位路径数字调制，是连续相位频移键控(CPFSK)的一种特殊情况，有时也叫做最小频移键

38、控(MSK)，MSK的“最小(Minimum )”二字指的是这种调制方式能以最小的调制指数)获得正交的调制信号。MSK信号表达式可正交展开为下式，其调制和解调框图如3.3和3.4所示11。St=cosct+2Tbakt+xk=cosxkcost2Tbcosct-akcosxksint2Tbsinct (3.8)图3.3 MSK信号的产生方法之一 MSK信号的解调 差分译码并串转换90积分清洗判决相关载波判决积分清洗MSK输入数据输出实际解调器往往需要解决载波恢复时的相位模糊问题，因此在编码器中，采用差分编码的预编码是必要的，同时在接收端必须在正交相干解调器输出段也要附加一个差分译码器。MSK解

39、调器的原理框图如图3.4所示11。图中，Xt=bItcost2Tbcos0t+bQtsint2Tbsin0t。定时时钟速率为12Tb，需要一个专门的同步电路来提取，如用平方环、判决反馈环、逆调制环等。图3.4 MSK相干解调框图3.2 高斯滤波最小频移键控（GMSK）GSKM作为一种高效的调制技术，是从OQPSK，MSK调制的基础上发展起来的一种数字调制方式，GMSK的很多方面都优于OQPSK和MSK，比如频带更窄，实现起来更简单，抗干扰能力更强。其特点是在数据流送交频率调制器前先通过一个Gauss滤波器(预调制滤波器)进行预调制滤波，以减小两个不同频率的载波切换时的跳变能量，使得在相同的数据

40、传输速率时频道间距可以变得更紧密，因此GMSK信号比MSK信号具有更窄的带宽。由于数字信号在调制前进行了Gauss预调制滤波，调制信号在交越零点不但相位连续，而且平滑过滤。GMSK调制的信号频谱紧凑、误码特性好，在数字移动通信中得到了广泛使用12。GMSK信号是在MSK调制信号的基础上发展起来的，MSKFSK信号。尽管MSK它具有包络恒定、相位连续、相对较窄的带宽和能相干解调的优点，但它不能满足某些通信系统对带外辐射的严格要求。为了压缩MSK信号的功率谱，在MSK调制前增加一级预调制滤波器，从而有效的抑制了信号的带外辐射。预调制滤波器应具有的特性：带宽窄而带外截止尖锐，以抑制不需要的高频分量；

41、脉冲响应的过冲量较小，防止调制器产生不必要的瞬时频偏；输出脉冲响应曲线的面积应对应于1/2的相移量，使调制指数为1/2。因此，GMSK采用满足以上条件的高斯滤波器作为脉冲形成的滤波器。数据通过高斯滤波器，然后进行MSK调制，滤波器的带宽由时间带宽常数BT决定。在没有载波漂移以及邻道的带外辐射功率相对与总功率小于60dB的情况下，选择BT0.28比较适合于常规的(IEEE定义频段为3001000MHz)移动无线通信系统。预制滤波器的引入使得信号的频谱更为紧凑，但是它同时在时域上展宽了信号脉冲，引入了码间干扰(ISI)，具体的说，预调制滤波器使得脉冲展宽，使得波形在时域上大于码元时间T。因此，有时

42、候将GMSK信号归入部分响应信号。 GMSK调制原理高斯低通滤波器的脉冲响应h(t)可以表示为ht=12Texp-t222T2 （3.9）其中=ln22BT，B是滤波器的3dB带宽。GMSK调制信号为11：st=2EbTbcoswot+-dn-g-nTbd+ dn=-+1, 0tLTb （3.10）gt=12TbrecttTb*hg=12TbQ2Bt-Tb2ln2-Q2Bt+Tbln2-Q2Bt+Tb2ln2Qt=t12e-22d （3.11）其中dn为发送信号序列，Eb为码元能量，Tb为符号周期，L为高斯滤波冲击响应长度，hgt为预高斯成形函数，B为高斯滤波器的3dB带宽，wo为载波频率，为

43、载波相位。对于BTbsbt=N=-expjhn=-Ndn+Cot-NTb +n=-NK=17expjhn=-NAk,N+ckt-NTb （3.12）其中：Cot=n=03sint+nTbsinht=-tgd，tLTbh-t-LTbgd，tLTb （3.13）其中Ak,N为系数。Cot的能量占GMSK信号能量的99.38%，对于更大的BTb，Cot项所占的比重更大，故可以忽略Sbt的后半部分，GMSK基带信号近似表示为：sbtN=-expjhn=-Ndn+c0t-NTb=k=0wkc0t-kTb+jk=0zkc0t-kTsbtUfC0f （3.14）预编码高斯低通滤波器MSK调制器带输入数据序列

44、GMSK调制的IF信号其中wk=cos2n=0kdn+,zk=-sin2n=0kdn+,C0f为C0t的频域响应协议GSM05.04V8.0.0中定义了GMSK调制方式，如图3.5所示。 图3.5 GMSK信号产生原理高斯滤波器的输出脉冲经MSK调制得到GMSK信号，其相位路径由脉冲的形状决定。由于高斯滤波后的脉冲无陡峭沿，也无拐点，因此，相位路径得到进一步平滑，如图3.6所示11。图3.6 GMSK的相位路径实现GMSK调制的关键是滤波器的设计。为了方便GMSK的解调，在调制之前需要对输入数据进行差分预编码。设输入数据为di0,1，di_diff=didi-1 (di0,1) ai=1-2d

45、i-diff (ai-1,1) （3.15）其中代表模2加，将差分编码之后的不归零数据ai，通过高斯低通滤波器和VCO，即可输出GMSK调制信号。高斯预调制滤波器的传递函数为：H2f=exp-2f2 （3.16）式中*Bb=ln22=0.5887。Bb是滤波器的3dB带宽，BbTb为系统中可变参数，BbTb取的小，能够使调制后的带宽变窄，但会引起码间干扰。BbTb=时即为MSK调制。高斯预调制滤波器的冲击响应函数为： ht=exp-t222Tb22*Tb （3.17） 式中=ln22BbTb。 高斯预调制滤波器的矩形脉冲响应为:gt=ht*rectt其中函数rectt=1,fortTb20,o

46、therwise,其中*代表卷积。将上式代入上上式中，得到： gt=Q2Bbln2t-Tb2-Q2Bbln2t+Tb2 （3.18）式中，Qt=t12e-22d。在欧洲GSM标准中，信道传输速率为:1Tb=16256ksymb/s=270.833ksymb/s=16256kbit/s BbTb=0.3 （3.19）当BbTb=0.3时，横坐标每格表示1个码元宽度。gt的积分满足：-+gtdt=Tb据式(3.18)可得GMSK信号的表达式为xt=cosct+2Tb-taig-iTb-Tb2d （3-20）式中：c为载波角频率；Tb为比特宽度；ai为输入的不归零数据3。 四相相移键控（QPSK）

47、四相相移键控信号简称“QPSK”， 意为正交相移键控，是一种数字调制方式。它分为绝对相移和相对相移两种。由于绝对相移方式存在相位模糊问题，所以在实际中主要采用相对移相方式QDPSK。它具有一系列独特的优点，目前已经广泛应用于无线通信中，成为现代通信中一种十分重要的调制解调试4。 QPSK的基本原理在数字相位调制中，M个信号波形可表示为 Smt=Regtej2m-1Mej2ft =gtcos2fct+2Mm-1 =gtcos2Mm-1cos2fct-gtsin2Mm-1sin2fct （m=1,2,，M，0t=0.5) r(i)=1; else r(i)=-1; end;end; %将这9999

48、个数整数化为 1，-1for i=1:n-count % 取9900个 R(i)=r(i-1)-mod(i-1),count)/count+1);end; %a(1)=1;for i=2:N if r(i-1)=1 a(i)=-a(i-1); else a(i)=a(i-1); end;end;for i=1:n A(i)=a(i-1)-mod(i-1),count)/count+1);end;for i=1:N if(mod(i,2) l(i)=a(i); l(i+1)=a(i); else q(i)=a(i); q(i-1)=a(i); end;end;for i=1:n I(i)=l(i

49、-1)-mod(i-1),count)/count+1);end;for i=1:n Q(i)=q(i-1)-mod(i-1),count)/count+1);end; ts=0.0001; Ts=ts*count; fs=1/ts; fc=5/(2*Ts); t0=Ts*N-ts; f=1/(4*Ts); df=0.3; tI=-Ts:ts:t0-Ts; tQ=0:ts:t0; tQ_R=0:ts:t0-Ts; c=2*pi*f; u=I.*cos(c*tI); v=Q.*sin(c*tQ); U=u.*cos(2*pi*fc*tI); V=v.*sin(2*pi*fc*tQ); for i

50、=count+1:n W(i)=U(i)+V(i-count); end; figure subplot(4,1,1); plot(tQ_R,R,k);axis(-Ts,0.1,-2,2); title(数字信号); grid on subplot(4,1,2); plot(tI,A,k);axis(-Ts,0.1,-2,2); title(差分编码后的信号); grid on subplot(4,1,3); plot(tI,I,k);axis(-Ts,0.1,-2,2); title(同相调制信号I); grid on subplot(4,1,4); plot(tQ,Q,k);axis(-T

51、s,0.1,-2,2); title(正交调制信号Q); grid onfigure subplot(3,1,1); plot(tI,U,k);axis(-Ts,0.1,-2,2); title(Icos(wct)cos(wt);grid on subplot(3,1,2); plot(tQ,V,k);axis(-Ts,0.1,-2,2); title(Qsin(wct)sin(wt);grid on subplot(3,1,3); plot(tI,W,k);axis(-Ts,0.1,-2,2); title(MSK信号);grid on noise=0.15*randn(1,n); S=W+

52、noise; S1,s,df1=fftseq(S,ts,df); S1=S1/fs; f=0:df1:df1*(length(s)-1)-fs/2; n_cutoff=100; H=zeros(size(f); for i=floor(fc-n_cutoff)/df1):floor(fc+n_cutoff)/df1) H(i)=1; end; for i=(65536-floor(fc-n_cutoff)/df1):-1:(65536-floor(fc+n_cutoff)/df1), H(i)=1; end; DEM=S1.*H; dem=real(ifft(DEM)*fs; dm=dem(1

53、:length(tQ); ds1=dm.*cos(2*pi*fc*tQ+pi*tQ/(2*Ts); ds2=dm.*cos(2*pi*fc*tQ-pi*tQ/(2*Ts); DW1,ww,df1=fftseq(ds1,ts,df); DW2,ww,df1=fftseq(ds2,ts,df); DW1=DW1/fs; DW2=DW2/fs; n_cutoff1=floor(200/df1); f1=0:df1:df1*(length(ww)-1)-fs/2; H1=zeros(size(f1); H1(1:n_cutoff1)=2*ones(1,n_cutoff1); H1(length(f1)

54、-n_cutoff1+1:length(f1)=2*ones(1,n_cutoff1); DEM1=H1.*DW1; DEM2=H1.*DW2; s1=real(ifft(DEM1)*fs; s11=s1(1:length(tQ); s2=real(ifft(DEM2)*fs; s22=s2(1:length(tQ); for i=1:n if (s11(i)=0.95)|(s11(i)=0.95)|(s22(i)=0.5) r(i)=1; else r(i)=-1; end;end;r(1:10)=1 -1 -1 1 1 1 -1 1 -1 -1;for i=1:n-count, R(i)=

55、r(i-1)-mod(i-1),count)/count+1);end;a(1)=1;for i=2:N if r(i-1)=1 a(i)=-a(i-1); else a(i)=a(i-1); end;end;for i=1:n A(i)=a(i-1)-mod(i-1),count)/count+1);end;for i=1:N if(mod(i,2) l(i)=a(i); l(i+1)=a(i); else q(i)=a(i); q(i-1)=a(i); end;end;for i=1:n I(i)=l(i-1)-mod(i-1),count)/count+1);end;for i=1:n

56、Q(i)=q(i-1)-mod(i-1),count)/count+1);end; ts=0.0001; Ts=ts*count; fs=1/ts; fc=3/(2*Ts); t0=Ts*N-ts; f=1/(4*Ts); df=0.3; tI=-Ts:ts:t0-Ts; tQ=0:ts:t0; tQ_R=0:ts:t0-Ts; c=2*pi*f; u=I.*cos(c*tI); v=Q.*sin(c*tQ); U=u.*cos(2*pi*fc*tI); V=v.*sin(2*pi*fc*tQ); for i=count+1:n W(i)=U(i)+V(i-count); end; figur

57、e subplot(3,1,1); plot(tQ_R,R,k);axis(0,0.1,-2,2); %title(数字信号); grid on subplot(3,1,2); plot(tI,W,k);axis(0,0.1,-2,2); %title(MSK信号);grid on附录2 QPSK仿真程序%信源为双极性不归零码（NRZ），取样点数自己定义QPSKglobal dt df t f N T %定义全局变量close all %关闭以前的应用窗口clear Eb_N0 Pe% % k=input(取样点数2k, k=13:); %采样点数自定义，默认为213 % if k=, k=1

58、3; end% alpha=input(滚降系数);% db=input(请选择输出眼图信噪比0-14:);% if db= | db14 | db0% db=10;% end;% fc=input(调制频率 in MHz:);%调制信号频率% A=input(输入载波幅度);k=13;alpha=0.5;db=10;fc=10;A=1;N=2k; %总的取样点数L=64;%L为每个码元的取样点数M=N/L;%M码元总数Rb=2;%码元速率Rb=2Mb/sTs=1/Rb;%码元宽度Tsdt=Ts/L;%时域取样间隔df=1/(N*dt);%频域取样间隔T=N*dt; %时间截断长度Bs=N*d

59、f/2; %带宽Na=4; %示波器扫描宽度为4个码元%alpha=0; %升余弦滚降系数t=linspace(-T/2,T/2,N); %频域横坐标f=linspace(-Bs,Bs,N)+eps;%时域横坐标%升余弦hr1=sin(pi*t/Ts)./(pi*t/Ts);hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).2);hr=hr1.*hr2;HR=abs(t2f(hr); %变为H(w),取模 GT=sqrt(HR);GR=GT;%。两个载波m=A*cos(2*pi*fc*t);n=-A*sin(2*pi*fc*t); EP=zeros(size(

60、f); %发送功率谱EPa=zeros(size(f); EPt=zeros(size(f);EPtq=zeros(size(f);EPti=zeros(size(f);EPtb=zeros(size(f);for loop1=1:15 %作误码率曲线时用到 Eb_N0(loop1)=2*(loop1-1); %Eb/N0 in dB eb_n0(loop1)=10(Eb_N0(loop1)/10); Eb=0.25;%单位比特能量（相当于能量谱） n0=Eb/(eb_n0(loop1); %信道的噪声谱密度 sita=n0*Bs; %信道中噪声功率 n_err=0; %误码计数 for lo