MSK调制解调器及同步性能的仿真分析.doc

2西南交通大学本科毕业设计（论文） 第 页重庆交通大学信息科学与工程学院课程设计报告题 目MSK调制解调器及同步性能的仿真分析专 业班 级学 号姓 名指导教师信息科学与技术学院二一四年九月 15信息科学与工程学院通信工程专业方向课程设计 通信工程专业课程设计成绩鉴定表进行时间2013-2014学年下学期第1、2周成绩鉴定学习内容（20分）与教学任务计划结合程度（10分）与专业培养结合程度（5分）其它（5分）过程评价（20分）实践能力（10分）学习态度（5分）学习纪律（5分）报告鉴定（60分） 报告内容与实践过程紧密结合（15分）报告内容与教学计划内容紧密结合（15分）报告质量(主题、结构、观点、逻辑、资料、字数 30分)指导教师姓名吴仕勋职称 讲师成绩摘 要 最小频移键控又称快速频移键控，他可以解决OQPsK调制方式不能解决包络起伏的问题，从而能够产生恒定包络、相位连续的调制信号。最小频移键控(MSK)是2FSK的改进调制方式，它具有波形连续，相位稳定，带宽最小并且严格正交的特点，其改进型GMSK在无线通信GSM系统中得到了广泛地应用。之所以称为最小频移键控，是由于其频移指数=0.5，即他的两个频率差是最小的正交频率差。本文研究了最小频移键控系统MSK调制与解调的工作原理，并给出了基于Matlab软件环境的仿真实现仿真运行结果。关键词：数字通信；调制；解调；系统仿真目 录摘 要II目 录III第1章 绪论11.1 数字通信的发展11.2 课题背景及意义1第2章 MSK调制解调原理32.1 MSK信号的产生32.2 MSK信号调制解调方法52.3 msk通信系统调制解调原理框图7第三章 仿真分析9结论11参考文献12附录13 信息科学与工程学院通信工程专业方向课程设计 第1章 绪论1.1 数字通信的发展通信按照传统的理解就是信息的传输与交换,为了传递消息，各种消息需要转换成电信号，消息与电信号之间必须建立单一的对应关系，否则在接收端就无法复制出原来的消息。通常，消息被载荷到电信号的某一参量上，如果电信号的该参量携带着离散消息，则该参量必将是离散取值的。这样的信号就称为数字信号。如果电信号的参量连续取值，则称这样的信号为模拟信号。按照信道中传输的是模拟信号还是数字信号，可以相应地把通信系统分为两类：数字通信系统和模拟通信系统。自1844年5月24日莫尔斯在华盛顿和巴尔的摩之间发送世界上斯一份电报以来 ，电报通信已经经历了150多年。但是长期以来，由于电报通信不如电话通信方便，作为数字通信主要形式的电报却比1876年贝尔发明的电话发展缓慢。直到20世纪60年代已后，数字通信才日益兴旺起来，数字通信迅速发展的基本原因是它与模拟通信相比，更能适应对通信技术越来越高的要求。第一数字传输抗干扰能力强，尤其是在中继时，数字信号可以再生而消除噪声的积累；第二，传输差错可以控制，从而改善了传输的质量；第三，便于使用现代数字信号处理技术来对数字信息进行处理；第四，数字信息易于做高保密性的加密处理；第五，数字通信可以综合传递各种消息，使通信系统功能增强。然而，数字通信的许多优点都是用比模拟通信占据更宽的带宽的系统频带而换来的。以电话为例，一路模拟电话只占据4khz的带宽，而一路传输质量相同的数字电话这可能要占用数十千赫兹的带宽。 在系统频带紧张的场合，数字通信这一缺点显得很突出，但是在系统频带富裕的场合，比如毫米波通信，光通信等场合，数字通信几乎成了唯一的选择。随着计算机技木和大规模集成技术的发展，数字通信在其发展过程中表现出了强大的生命力，它冲破了传统模拟通信方式的统治，逐步地发展、完善。可以预言：随着通信事业的发展，特别是各种宽带传输技术(例如光纤传输、数字微波等)、综合业务数字网(ISDN)的实用化，全数字化的通信方式必将逐步取代模拟通信方式而得到蓬勃发展1。1.2 课题背景及意义 当今社会已经步入信息时代，在各种信息技术中信息的传输及通信起着支撑作用。而在频带资源日益紧张的今天，为了提高系统的容量（满足更多的用户）信道间隔已经是一减再减已经由最初的100khz减到了今天的12.5khz甚至更小。数字通信系统因其组网灵活，差错控制和保密性都比较容易，而且能够进入ISDN网所以通信系统已逐步由模拟制式向数字制式过渡，信号的调制方式也逐步由模拟方式持续、广泛地向数字方式转化，数字通信系统成为了信息的传输的一种重要手段。 然而，一般的数字调制技术，如ASK、PSK和FSK因传输效率低和抗干扰能力差而无法满足移动通信的要求，为此，需要专门研究一些抗干扰性强、误码性能好、频谱利用率高的数字调制技术，尽可能地提高单位频谱内传输数据的比特率，以适用于移动通信窄带数据传输的要求。MSK因具有：（1）已调信号振幅是恒定的。（2）信号的频率偏移严格等于1/4TS ，相应的调制指数H=（f2-f1）ts =0.5。(3)以载波相位在一个码元期间内准确地线性变化pi/2；（4）在一个码元期间内，信号应包含四分子一载波周期的整数倍。(5)在码元转换时刻信号的相位是连续的，或者说，信号的波形没有突变。的特点使得msk通信系统抗干扰能力强适用于移动通信等窄带数据传输的要求。在设计新系统或者对原有的通信系统做出修改或者进行相关的研究时，通常要进行建模和仿真，通过仿真结果衡量方案的可行性，从中选择最合理的系统配置和参数设置，然后再应用于实际系统中。通过仿真，可以提高研究开发工作的效率，发现系统中潜在的问题，优化系统整体性能。与一般的仿真过程类似，在对通信系统实施仿真之前，首先需要研究通信系统的特性，通过归纳和抽象建立通信系统的仿真模型。通过对系统的仿真,可以不需要实际的硬件环境就可以分析系统的特点。人们能够通过仿真实验就可以了解msk数字通信系统性能。这样大大的减少实验的开销，对科学技术的发展是很重要的。matrix公司的matlab软件是一套功能非常强大的工程技术数值运算和系统仿真软件。第2章 MSK调制解调原理2.1 MSK信号的产生 频移键控是数字通信中用得较广的一种形式，在衰落信道中传输数据时，它被广泛采用。Fsk信号是0符号对应载频1，而1符号对应于载频2（与1不同的另一载频）的已调波形，而且1与2之间的改变是瞬间完成的。基本调制方法有模拟调频法和键控法.一般来说，键控法得到的得到的调制信号的相位是不连续的（两载波频率相差为pi/2的整数倍时相位连续） 图2-1 二进制移频键控信号的时间波形MSK叫最小移频键控，它是移频键控（FSK）的一种改进型。这里“最小”指的是能以最小的调制指数（即0.5）获得正交信号，它能比PSK传送更高的比特速率。二进制MSK信号的表达式可写为： (2-1) 式中，k称为附加相位函数；c为载波角频率；Tk为第k个输入码元，s为码元宽度；a取值为1；k为第k个码元的相位常数，在时间kTst(k+1)Ts中保持不变，其作用是保证在t=kTs时刻信号相位连续。由 (2-2)可知当1时，信号的频率为：当1时，信号的频率为：由此可得频率之差为：H=Ts=x Ts=0.5那么MSK信号波形如图示： 图2-2 MSK信号波形从图中可以看出，+信号和信号在一个码元期间恰好相差二分之一周，即相差为了保持相位的连续，在t=时间内应有下式成立：=（）（） (2-3) 即：当时，=；当时，=（）； (2-4) 若令0，则0或，此式说明本比特内的相位常数不仅与本比特区间的输入有关，还与前一个比特区间内的输入及相位常数有关。由附加相位函数k(t)的表示式可以看出，k(t)是一直线方程，其斜率为 截距为k。由于ak的取值为1，故是分段线性的相位函数。因此，MSK的整个相位路径是由间隔为Ts的一系列直线段所连成的折线。在任一个码元期间Ts，若ak=+1，则k(t)线性增加；若ak=-1， 则k(t)线性减小 。对于给定的输入信号序列ak，相应的附加相位函数k(t)的波形如图所示。对于各种可能的输入信号序列，k(t)的所有可能路径是一个从-2到+2的网格图。 图3-3 附加相位函数k(t)的波形图 图4-4 附加相位路径网格从以上分析总结得出，MSK信号具有以下特点： （1）MSK信号是恒定包络信号； （2）在码元转换时刻，信号的相位是连续的，以载波相位为基准的信号相位在一个码元期间内线性地变化/2 ； （3） 在一个码元期间内， 信号应包括四分之一载波周期的整数倍，信号的频率偏移等于，相应的调制指数h=0.5。 (4)信号频率偏移严格等于 。 2.2 MSK信号调制解调方法二进制移频键控信号fsk很容易用一个矩形脉冲序列对一个载波进行调频而获得。这正是频率键控通信方式早期采用的实现方法，也是利用模拟调频法实现数字调频的方法。对2fsk信号产生的另一方法是采用键控法，即利用受矩形脉冲序列控制的开关电路对两个不同的独立频率源进行选通4。二进制移频键控信号的解调方法很多，有模拟鉴频法和数字检测法，有非相干解调方法也有相干解调方法。其解调原理是将二进制移频键控信号分解为上下两路二进制振幅键控信号，分别进行解调，通过对上下两路的抽样值进行比较最终判决出输出信号。 图5-5二进制移频键控信号解调器原理图 (a) 非相干解调; (b) 相干解调MSK信号属于数字频率调制信号，可以采用相干解调，也可以采用非相干解调方式。由 = (2-5) 令， 则： (2-6)所以msk信号可以看成是由两个彼此正交的载波和分别进行振幅调制而合成的。根据上面表达式的描述可构成一种MSK调制器由MSK信号的一般表示式可得 (2-7)因为 代入式可得 (2-8) 上式即为MSK信号的正交表示形式。其同相分量为 也称为I支路。 其正交分量为 也称为Q支路。Cos 和sin 称为加权函数。Q支路信号先延迟Ts，经sinct加权调制和正交载波sinct相乘输出正交分量xQ(t)。xI(t)和xQ(t)相减就可得到已调MSK信号。 MSK信号属于数字频率调制信号，可以采用相干解调，也可以采用非相干解调方式。由于MSK信号调制指数较小，采用一般鉴频器方式进行解调误码率性能不太好，因此在对误码率有较高要求时大多采用相干解调方式。本设计模块中采用一种相干解调的方式。已知： (2-9)把该信号进行正交解调可得到：Ik路 (2-10)=+Qk路 (2-11)=+我们需要的是、两路信号，所以必须将其它频率成份、通过低通滤波器滤除掉，然后对、采样即可还原成、两路信号。根据上面描述可构成MSK解调器.2.3 msk通信系统调制解调原理框图 到目前为止，已经把整个设计做了详尽的阐述了，接下来便只是需要把以上的几个设计模汇总了，形成一个完整的原理图。 介绍MSK调制的基本原理，研究MSK信号的特点及一种正交实现MSK信号的方法，并对这方法进行理论的分析和仿真。 图2-6 Msk通信系统调制原理图 图2-7 Msk通信系统解调原理图第三章 仿真分析 打开matlab simulink启动msk通信系统仿真框图程序.首先观察两种系统调制信号波形： 图3-1 Msk的传输数据和差分后的波形 图3-1为MSK传输数据，和经过差分变换后的数据波形。 图3-2 I路Q路波形 下面我们再来分析一下，msk信号可以看成是由两个彼此正交的载波和分别进行振幅调制而合成的。经sinct加权调制和正交载波sinct相乘输出正交分量xQ(t)。xI(t)和xQ(t)相减就可得到已调MSK信号。 图3-3 Msk的已调和加噪后的波形 通过上面的仿真结果，我们可以看到通过信道后，信号上叠加了很大的噪声，使信号产生了很大的畸变。从msk通信系统传输方针结果来看信号在传输过程中高低电平宽度基本一致，在误码上没有很明显的误判现象，在接收端接收信号只是偶尔会出现毛刺。通过比较我们看出在衰落很严重和干扰的窄带通信系统中使用msk调制解调的通信系统更具优势. 图3-4 Msk解调后的波形和误码曲线从3-4中可看出，解调后的波形与传输的数据波形几乎一致，MSK传输数据经过调制后在进行加噪，模拟传输信道的噪声干扰，然而在解调之后，MSK解调数据与传输数据间的误差很小，可看出MSK调制的抗干扰性能确实很好。同时图3-5的误码曲线也证明如此。 结论 通信系统性能评价主要从：1.抗噪声性能（从误码率体现）。2.调制信号频带宽度。3.信道特性变化的敏感性来衡量。MSK调制方式具有相位连续、包络恒定、频谱利用率高等优点。在通信系统中低误码率的MSK调制解调器硬件方案一直为研究的热点。可见msk调制解调数字通信系统更适合应用于信道衰落严重的通信系统。通过实验仿真我们可以设计出合理的系统参数。仿真技术作为一门综合性科学，将随着其相关领域技术的深入发展，继续向纵深快速发展，同时将扩大其综合应用的领域，发挥更大的作用。但是，作为一门综合性技术学科，仿真技术还有许多理论及技术问题需要继续进行深入的研究探讨。我国应大力开展仿真技术的理论研究和技术应用研究，尽快缩短与先进发达国家在技术上的差距。仿真技术的发展，必将推动我国科学技术水平的进一步提高。参考文献 1张志勇等.精通matlab6.5J.北京:北京航空航天大学出版社,2003.3.3樊昌信,张甫詡、徐炳祥、吴成柯 .通信原理M .北京:国防工业出版社,2001.5.4王立宁,乐光新,詹菲等.MATLAB与通信仿真M .北京:人民邮电出版社, 1999 5 陈怀琛，吴大正，高西金 matlab 及在信息课程中的应用J北京：电子工业出版社,20066 楼天顺、刘晓东、李博涵基于matlab7.x的系统分析与设计M.西安：西安电子科技大学出版社, 2005.5 7张肃文主编 .高频电子线路M .北京：高等教育出版社, 2004.11 8程佩青编著 .数字信号处理教程M .北京：清华大学出版社, 2006.3 9张辉、曹丽娜、王勇 编著 .通信原理辅导M 西安.西安电子科技大学出版社, 2004.11 10徐明远 .MATLAB仿真在通信工程中的应用J. 西安. 西安电子科技大学出版社, 2005.6 11田丽华等.编码理论J.西安 西安电子科技大学出版社, 2003.8附录部分代码：%MSK调制，解调 clear all data_len = 30000; %码元个数 sample_number = 8;%采样个数 Rb = 24000; %码元速率fc = 96000; %载波频率 data = rand_binary(data_len); %MSK基带调制 signal_out,I_out,Q_out = mod_msk(data,data_len,sample_number,Rb); %中频搬移 multi = fc/Rb; I_temp=interp(I_out,multi); Q_temp=interp(Q_out,multi);Fs=fc*sample_number; t=1/Fs:1/Fs:length(I_temp)*1/Fs; signal_i=I_temp.*cos(2*pi*fc*t); signal_q=Q_temp.*sin(2*pi*fc*t); signal_mod=I_temp.*cos(2*pi*fc*t)-Q_temp.*sin(2*pi*fc*t); % signal_mod1=real(I_temp+j*Q_temp).*exp(j*2*pi*fc*t); %加噪声 for SNR = 0:8 signal_mod1 = 0.01 * awgn(signal_mod,SNR); N=300; % 滤波器的阶数为(N+1)? ? F=0,fc-1000,fc+1000,Fs/2*2/Fs; A=1,1,0,0; lpf=firls(N,F,A); amp_lpf,w=freqz(lpf); I_dem=signal_mod1.*cos(2*pi*fc*t)*2 .* cos(pi*t*Rb/2); I_d