基于simulink的gmsk仿真

1、专业课程设计报告 题 目： 基于Matlab的基带GMSK的仿真研究 姓 名：陈国强专 业：通信工程班级学号：09042112同 组 人：指导老师：夏思满 南昌航空大学信息工程学院 20 12 年 6 月 20 日 专业 课程设计任务书20 11 20 12 学年 第 2 学期第 17 周 20 周 题目基于Matlab的基带GMSK的仿真研究内容及要求设计要求：（1）掌握GMSK的原理和Simulink仿真基本方法；（2）通过SIMULINK对BT=0.3的GMSK调制解调系统进行仿真；（3）观察调制信号和已调信号波形；（4）改变BT参数，分析调制性能和BT参数的关系。进度安排 17周：查找

2、资料，进行系统软件方案设计； 18周：软件的分模块调试； 19周：系统联调； 20周：设计结果验收，报告初稿的撰写。学生姓名：指导时间 2011.62011.7指导地点： E楼 610 室任务下达2011 年 6 月 13 日任务完成2011 年 7 月 8 日考核方式1.评阅 2.答辩 3.实际操作 4.其它指导教师夏思满系（部）主任付崇芳基于Matlab的基带GMSK的仿真研究摘要：随着现代通信技术的发展，移动通信技术得到快速发展，许多优秀的调制技术应运而生，其中高斯最小频移键控（GMSK）技术是无线通信中比较突出的一种二进制调制方法，它具有良好的功率谱特性和较好的抗干扰性能，特别适用于无

3、线通信和卫星通信，目前，很多通信标准都采用了GMSK技术，例如，GSM，DECT等。本文首先介绍了MSK的一般原理，接着对GMSK的调制原理和几种调制方法进行了阐述，然后，重点研究了GMSK的几种差分解调方法并进行了比较，最后用Matlab软件中的simulink进行仿真，结果表明GMSK具有包络恒定、相位连续、频道干扰小、误码率较低等优点。关键词：高斯最小频移键控；调制；差分解调；Matlab；simulink目 录第一章 设计要求11.1 设计内容11.2 设计要求1第二章 系统的组成及设计原理2第三章 系统功能模块设计53.1 信号发生模块53.2 调制、解调模块53.3 误码率计算器6

4、3.4 波形观察模块7调制、解调信号观察模块73.4.2 调制信号频谱观察模块7眼图观察模块8第四章 系统调试与结果分析94.1 实验调试94.2结果分析104.2.1 GMSK调制与解调波形104.2.2 GMSK调制信号眼图13结论16参考文献17附录一 ：程序18附录二：GMSK调制解调建模图18第一章 设计要求1.1 设计内容：通过SIMULINK对BT=0.3的GMSK调制系统进行仿真。1.2 设计要求：（1）掌握GMSK的原理和Simulink仿真基本方法；（2）通过SIMULINK对BT=0.3的GMSK调制解调系统进行仿真；（3）观察调制信号和已调信号波形；（4）改变BT参数，

5、分析调制性能和BT参数的关系。第二章 系统的组成及设计原理GMSK系统主要由信号产生模块、信号调制模块、信道、信号解调模块、误码率计算模块组成。在图形观察方面还包含频谱仪、示波器和眼图绘制模块。本系统由信号产生模块产生一个二进制序列，再经过调制器进行调制，之后便将调制信号送入信道，经过解调器解调得到解调信号。为计算系统误码率，则在调制器后加一误码率计算模块，计算误码率。信道解调模块误码率计算模块频谱仪示波器调制模块信号产生模块图2.1 系统原理框图在设计中，选用贝努力二进制序列产生器来产生器（Bernoulli Binary Generator）产生一个二进制序列，将序列送入GMSK基带调制器

6、模块（GMSK Modulator Baseband）中得到已调信号，再将已调信号送入一个加性高斯白噪声信道，将信噪比设为一个变量，用于绘制信噪比误码率曲线。解调阶段则将通过加性高斯白噪声信道的信号输入GMSK基带解调器模块（GMSK Demodulator Baseband）中，其后接一误码率统计模块（Error Rate Calculation），且误码率统计模块另一输入端接至源信号处。而用示波器观察解调波形并与源信号波形进行比较。因为已调信号是一复合信号，所以要用complex to Magnitude-Angle 模块，再用示波器分别观察其幅度与相角。另外还用频谱仪观察了已调信号的频谱

7、。GMSK调制调制原理图如图2.2，图中滤波器是高斯低通滤波器，它的输出直接对VCO进行调制，以保持已调包络恒定和相位连续。非归零数字序列高斯低通滤波器频率调制器（VCO）GMSK已调信号图2.2 GMSK调制原理图为了使输出频谱密集，前段滤波器必须具有以下待性：1.窄带和尖锐的截止特性，以抑制FM调制器输入信号中的高频分量；2.脉冲响应过冲量小，以防止FM调制器瞬时频偏过大；3.保持滤波器输出脉冲响应曲线下的面积对应丁pi2的相移。以使调制指数为12。前置滤波器以高斯型最能满足上述条件，这也是高斯滤波器最小移频键控(GMSK)的由来。GMSK解调GMSK本是MSK的一种，而MSK又是是FSK

8、的一种，因此，GMSK检波也可以采用FSK检波器，即包络检波及同步检波。而GMSK还可以采用时延检波，但每种检波器的误码率不同。GMSK非相干解调原理图如图2.3，图中是采用FM鉴频器（斜率鉴频器或相位鉴频器）再加判别电路，实现GMSK数据的解调输出。带通滤波器限幅器鉴频器判决器数据GMSK信号图2.3 GMSK解调原理图如图2.4为GMSK调制解调系统的SimuLink仿真模型，整个系统主要包括五大模块：随机信号发生模块、GMSK调制模块、信道、GMSK解调模块、误码率统计模块。所选库模块如图2.4中所示。图2.4 系统SimuLink仿真模型图第三章 系统功能模块设计3.1 信号发生模块

9、因为GMSK信号只需满足非归零数字信号即可，本设计中选用（Bernoulli Binary Generator）来产生一个二进制序列作为输入信号。GMSK信号输出图3.1 GMSK信号产生器该模块的参数设计这只主要包括以下几个。其中probability of a zero 设置为0.5表示产生的二进制序列中0出现的概率为0.5；Initial seed 为61表示随机数种子为61；sample time为1/1000表示抽样时间即每个符号的持续时间为0.001s。当仿真时间固定时，可以通过改变sample time参数来改变码元个数。例如仿真时间为10s，若sample time为1/100

10、0，则码元个数为10000。GMSK解调信号3.2 调制、解调模块GMSK信号图3.2 GMSK调制解调模块GMSK Modulator Baseband为GMSK基带调制模块，其input type参数设为Bit表示表示模块的输入信号时二进制信号（0或1）。BT product为0.3表示带宽和码元宽度的乘积。其中B是高斯低通滤波器的归一化3dB带宽，T是码元长度。当B·T=时，GMSK调制信号就变成MSK调制信号。BT=0.3是GSM采用的调制方式。Plush length则是脉冲长度即GMSK调制器中高斯低通滤波器的周期，设为4。Symbol prehistory表示GMSK调

11、制器在仿真开始前的输入符号，设为1。Phase offset 设为0，表示GMSK基带调制信号的初始相位为0。Sample per symbol为1表示每一个输入符号对应的GMSK调制器产生的输出信号的抽样点数为1。AWGN Channel为加性高斯白噪声模块，高斯白噪声信道的Mode参数（操作模式）设置为Signal to noise(SNR),表示信道模块是根据信噪比SNR确定高斯白噪声的功率，这时需要确定两个参数：信噪比和周期。而将SNR参数设为一个变量xSNR是为了在m文件中编程，计算不同信噪比下的误码率，改变SNR即改变信道信噪比。GMSK Demodulator Baseband是

12、GMSK基带解调器。其前六项参数与GMSK调制器相同，并设置的值也相同。最后一项为回溯长度Traceback Length,设为变量Tracebacklength，在m文件通过改变其值，可以观察回溯长度对调制性能的影响。3.3 误码率计算模块基带信号GMSK解调信号图3.3误码率计算模块Receive dely(接收端时延)设置为回溯长度加一，表示接收端输入的数据滞后发送端数据TracebackLength+1个输入数据；Computation delay(计算时延)设为0，表示错误率统计模块不忽略最初的任何输入数据。Computation mode(计算模式)设置为Entire frame(

13、帧计算模块)，表示错误率统计模块对发送端和接收端的所有数据进行统计。Output data(输出数据)设为workspace，表示竟统计数据输出到工作区。Variable name (变量名)则是设置m文件中要返回的参数的名称，设为xErrorRate。3.4 波形观察模块调制、解调信号观察模块因为GMSK调制信号是一个复合信号，所以只用示波器（Scope）无法观察到调制波形，所以在调制信号和示波器间加一转换模块Complex to magnitude-angle将调制信号分别在幅度和相角两方面来观察。GMSK调制信号图3.4 调制信号观察模块将Complex to magnitude-ang

14、leoutput的output参数设为magnitude and angle,表示同时输出调制信号的幅度和相角。示波器scope1的number of axes 为2表明有纵坐标个数为2；time range表示时间轴的显示范围，设为auto,表示时间轴的显示范围为整个仿真时间段。Tick Tabels 设为bottom axis only时，只显示各个纵坐标以及最下面的横坐标的标签。图3.5 解调信号观察模块 调制信号频谱观察模块图3.6 GMSK调制信号频谱观察模块设置了坐标Y的范围为0到7，X的范围为-FS,FS，Amplitude scaling表示幅度计算，选择一般模式即以V为单位进

15、行计算。但Y坐标标记Y-axis title设为magnitude，dB转换为dB形式。眼图观察模块图3.7 GMSK调制信号眼图观察模块Offset(sample)参数表示MATLAB在开始绘制眼图之前应该忽略的抽样点的个数。Symbols per trace表示每径符号数，每条曲线即成为一个“径”。Traces displayed 则是要显示的径数。New traces per display 是每次重新显示的径的数目。在系统中要求通过m文件编程绘制误码率曲线。其程序流程图如图3.8：开始X=0:10; xTracebackLength=4;Y=xi=1:length(x)执行SimuLi

16、nk 仿真模型Y(i)=xErrorRate(1)返回返回Y(i)=xErrorRate(1)图3.8 程序流程图第四章 系统调试与结果分析4.1 实验调试调试过程中主要通过MATLAB自带的Help功能来进行调试，在Help中查找所需函数的定义及形式和使用方法。通过报错信息找出相应的错误，翻阅相关资料，与同组人经过讨论后进行修改。在最终解决不了的情况下，请教老师，最终改正所有错误。设计模块、参数设置及程序代码编写完成后。先将高斯白噪声信道信噪比xSNR和GMSK解调模块的回溯长度参数设为常数，运行实验模型，观察示波器，发现没有出现基带与解调信号波形。先检查示波器参数，发现并无问题，编译Sim

17、uLink的.mdl文件时信号发生器报错，错误信息为：For integer inputs, the input values must be in the range +/- (M-2i-1), i=0,1, ., (M/2)-1，检查GMSK调制模块参数input type与GMSK解调模块out put 参数均设置为integer，但实际上贝努力二进制序列产生器产生的是一个由0和1组成二进制序列，与integer产生冲突，将上述两参数就改为bit，再编译mdl文件，无错误显示。进而运行m文件，mdl文件界面弹出，说明无法执行mdl模型。检查程序，发现xSimulation Time在m文

18、件中有设置，而此参数在SimuLink中的simulation/simulation parameters中已根据start time 和stop time 设定，删除m文件中的xSimulation Time=10，再运行，观察示波器，示波器显示波形。误码率曲线也能画出。署名系统基本功能已经实现。在执行瑞丽信道模块时，多径道瑞丽信道模块报错，报错信息为：egal rate transition found involving block 'ruili_error/Multipath Rayleigh Fading Channel/Multipath Fading Channel/Mu

19、ltiply with back propagation/S-Function' at input port 1. A Rate Transition block must be inserted between the two blocks；在资料上查找多径道瑞丽信道模块的参数，发现其Sample time参数必须设置为1/BitRate/SampleperSymbol，前面二进制序列发生器的sample time 为1/1000，而多径道瑞丽信道模块SampleperSymbol参数为1，故多径道瑞丽信道模块的Sample time参数应为1/1000。改正后，运行文件，无错。各调

20、制信号观察时，频谱仪显示的图形都与理论频谱形状相差很大，尤其GMSK的频谱，都没有出现主瓣与旁瓣的明显区分，重新修改频谱仪的参数，将Amplitude scaling 参数由dB改为magnitude，情况就好很多了。但是无论怎样改变，都不能得到理想的状态，估计是其他模块的一些参数对频谱仪的图形观察有影响。最后执行总的文件，各模块都能顺利执行，说明软件调试基本完成。4.结果分析4.2.1 GMSK调制与解调波形:图4.1 GMSK调制信号幅度和相角波形由于调制信号时一个复合信号，不能直接由示波器观察，通过一complex to magnitude-angle模块将调制信号分为幅度和相角两个变量

21、来观察。通过幅度的波形（上）验证了GMSK的幅度不变，由相角波形（下）来看，相角连续，与理论符合。所以图形基本正确。图4.2 GMSK基带信号与解调信号由图4.2中基带信号（上）与解调信号波形（下）比较可得，其由起始码元到最后一个码元，发现调制信号波形从第四个码元开始与基带信号完全符合，说明系统的调制性能较好，基本实现了解调的目的将调制信号还原为基带信号。图4.3 BT=0.3的GMSK调制信号频谱图4.4 GMSK等理论调制频谱对比图4.3和图4.4，实验所得频谱图的主瓣与理论频谱近似，只是顶端稍显尖锐，不够圆滑，可能的频谱仪的参数或去其他模块参数设置不恰当。图4.5 BT=0.9的GMSK

22、调制信号频谱比较图4.3和图4.5中频谱，发现BT=0.3与BT=0.9得GMSK调制频谱，并无明显差异，与GMSK调制信号的频谱随着B.T的减小而变得紧凑起来的理论结果不符合，从而验证可能是系统的某些参数设置不太合理，导致得不到正确的结果。图5.6MSK调制信号频谱比较图4.3和图4.5，发现GMSK的旁瓣衰减比MSK明显，也充分说明了GMSK频谱特性较MSK更好。4.2.2 GMSK调制信号眼图图4.7 BT=0.1分析：由图中混乱的线条可知，BT=0.1时，眼图“眼睛”睁开很小，失真严重，系统码间串扰较大。图4.8 BT=0.3分析：由图中混乱的线条可知，BT=0.3时，眼图“眼睛”睁开

23、比图5.8中大，但存在过零点失真，仍然存在码间串扰，但比BT=0.1时好得多。图4.9 BT=0.9分析：与图4.7，4.8相比较，图4.9中眼图最为清晰，眼睛睁开程度也较大，且眼图端正，说明码间串扰较小。综合上述分析，可知BT值越小，码间串扰越大，这也是GMSK体制的缺点。图4.11 不同BT值时的GMSK误码率曲线在BT=0.2、0.3、0.7时，对系统误码率进行仿真。比较三条曲线，可以看到其差别并不大。结果表明:不同BT值的信号调制性能差别不大.随着信噪比的增大,BT=0.2与BT=0.3的系统性能基本一致。当BT=0.3时,既可以使频域带宽很窄,时域持续时间适当,又使时域信号容易实现。

24、结论在本次专业课程设计中第一次接触到SimuLink ，刚开始是一头雾水毫无头绪。后来经过资料查阅逐渐了解了SimuLink是MATLAB提供的用于对动态系统进行建模、仿真和分析的工具包。设计中要求用SimuLink搭建GMSK调制与解调模块、计算误码率，并且绘制信噪比误码率曲线。第一星期内，主要进行了资料查询及建模任务，先将各模块参数均设为常数，对调制解调波形进行观察，分析器实际波形与理论波形间的差距以及产生误差的原因。依次采用分步进行的方式，逐步实现系统所需各个功能。调试过程采用的也是相同的方法，每搭建一个功能模块就先进行仿真，调试待得到满意结果后再进行下一功能模块的搭建和调试。并在不断出现错误的过程中学会了应用MATLAB的系统帮助。整个专业课程设计中遇到的最大问题就是FFT频谱仪的参数设置及仿真参数的设置，总是solver options得选择时出问题，把握不好固定步长和可变步长的选择，以及固定步长时连续求解器的选择。经实践证明，