MSK调制解调课件

1、通信工程专业综合课程设计目录1 绪论11.1 数字通信的发展11.2研究MSK数字通信系统的意义11.3通信系统仿真的意义12 MSK系统调制解调原理32.1 MSK调制原理32.2 MSK解调原理52.3 MSK的抗噪声性62.4 MSK功率谱密度73 MATLAB仿真设计与系统分析93.1 MATLAB简介9.3.2 MSK信号的调制93.3 MSK解调实现103.4叠加噪声的MSK系统104 MSK系统分析114.1 MSK信号的时域调制解调分析114.2 MSK频域分析124.3 加噪声系统分析124.4误码率分析135 总结14致谢15参考文献：16附录:MSK信号调制解调的源程序代

2、码17 1 绪论 1.1 数字通信的发展通信按照传统的理解就是信息的传输与交换,为了传递消息，各种消息需要转换成电信号，消息与电信号之间必须建立单一的对应关系，否则在接收端就无法复制出原来的消息。通常，消息被载荷到电信号的某一参量上，如果电信号的该参量携带着离散消息，则该参量必将是离散取值的。这样的信号就称为数字信号。如果电信号的参量连续取值，则称这样的信号为模拟信号。按照信道中传输的是模拟信号还是数字信号，可以相应地把通信系统分为两类：数字通信系统和模拟通信系统。自1844年5月24日莫尔斯在华盛顿和巴尔的摩之间发送世界上斯一份电报以来 ，电报通信已经经历了150多年。但是长期以来，由于电报

3、通信不如电话通信方便，作为数字通信主要形式的电报却比1876年贝尔发明的电话发展缓慢。直到20世纪60年代已后，数字通信才日益兴旺起来，数字通信迅速发展的基本原因是它与模拟通信相比，更能适应对通信技术越来越高的要求。第一数字传输抗干扰能力强，尤其是在中中继时，数字信号可以再生而消除噪声的积累；第二，传输差错可以控制，从而改善了传输的质量；第三，便于使用现代数字信号处理技术来对数字信息进行处理；第四，数字信息易于做高保密性的加密处理；第五，数字通信可以综合传递各种消息，使通信系统功能增强。然而，数字通信的许多优点都是用比模拟通信占据更宽的带宽的系统频带而换来的。以电话为例，一路模拟电话只占据4k

4、hz的带宽，而一路传输质量相同的数字电话这可能要占用数十千赫兹的带宽。 在系统频带紧张的场合，数字通信这一缺点显得很突出，但是在系统频带富裕的场合，比如毫米波通信，光通信等场合，数字通信几乎成了唯一的选择。随着计算机技木和大规模集成技术的发展，数字通信在其发展过程中表现出了强大的生命力，它冲破了传统模拟通信方式的统治，逐步地发展、完善。可以预言：随着通信事业的发展，特别是各种宽带传输技术(例如光纤传输、数字微波等)、综合业务数字网(ISDN)的实用化，全数字化的通信方式必将逐步取代模拟通信方式而得到蓬勃发展。 1.2研究MSK数字通信系统的意义 当今社会已经步入信息时代，在各种信息技术中信息的

5、传输及通信起着支撑作用。而在频带资源日益紧张的今天，为了提高系统的容量（满足更多的用户）信道间隔已经是一减再减已经由最初的100khz减到了今天的12.5khz甚至更小。数字通信系统因其组网灵活，差错控制和保密性都比较容易，而且能够进入ISDN网所以通信系统已逐步由模拟制式向数字制式过渡，信号的调制方式也逐步由模拟方式持续、广泛地向数字方式转化，数字通信系统成为了信息的传输的一种重要手段。然而，一般的数字调制技术，如ASK、PSK和FSK因传输效率低和抗干扰能力差而无法满足移动通信的要求，为此，需要专门研究一些抗干扰性强、误码性能好、频谱利用率高的数字调制技术，尽可能地提高单位频谱内传输数据的

6、比特率，以适用于移动通信窄带数据传输的要求。MSK因具有：（1） 已调信号振幅是恒定的；（2） 信号的频率偏移严格等于 ，相应的调制指数=0.5；（3） 以载波相位在一个码元期间内准确地线性变化；（4） 在一个码元期间内，信号应包含四分子一载波周期的整数倍；（5） 在码元转换时刻信号的相位是连续的，或者说，信号的波形没有突变。的特点使得MSK通信系统抗干扰能力强适用于移动通信等窄带数据传输的要求。1.3通信系统仿真的意义在设计新系统或者对原有的通信系统做出修改或者进行相关的研究时，通常要进行建模和仿真，通过仿真结果衡量方案的可行性，从中选择最合理的系统配置和参数设置，然后再应用于实际系统中。通

7、过仿真，可以提高研究开发工作的效率，发现系统中潜在的问题，优化系统整体性能。与一般的仿真过程类似，在对通信系统实施仿真之前，首先需要研究通信系统的特性，通过归纳和抽象建立通信系统的仿真模型。通过对系统的仿真,可以不需要实际的硬件环境就可以分析系统的特点。人们能够通过仿真实验就可以了解MSK数字通信系统性能。这样大大的减少实验的开销，对科学技术的发展是很重要的。matrix公司的MATLAB软件是一套功能非常强大的工程技术数值运算和系统仿真软件。MSK通信系统的仿真设计主要就是使用MATLAB的M函数工具箱进行仿真。 本次毕业设计的主要任务是对MSK数字通信系统进行MATLAB环境下的仿真。首先

8、通过收集资料和学习来理解MSK通信系统的工作原理，然后使用MATLAB软件对MSK数字通信系统的调制解调以及传输进行仿真，通过仿真结果分析得出系统的性能和优势。目的就是让人们可以通过本次仿真就可以了解到MSK数字通信系统的特点。2 MSK系统调制解调原理 2.1 MSK调制原理MSK称为最小移频键控，是移频键控（FSK）的一种改进型。这里“最小”指的是能以最小的调制指数（即0.5）获得正交信号，它能比PSK传送更高的比特速率。二进制MSK信号的表达式可写为： （公式 2.1)载波角频率；码元宽度；第k个码元中的信息，其取值为1；第k个码元的相位常数，它在时间中保持不变；当1时，信号的频率为：当

9、1时，信号的频率为：由此可得频率之差为：那么MSK信号波形如图2.1所示：图2.1 MSK信号波形为了保持相位的连续，在t=时间内应有下式成立：=（）（） （公式 2.2)即：当时，=；当时，=（）；若令0，则0或，此式说明本比特内的相位常数不仅与本比特区间的输入有关，还与前一个比特区间内的输入及相位常数有关。 （公式 2.3)= 令， 则： （公式 2.4)为了便于理解如图2.2所示：图2.2 码元变换及成形信号波形图根据上面描述可构成一种MSK调制器，其方框图如图2.3所示：图2.3 MSK调制原理框图 输入数据NRZ，然后通过CPLD电路实现差分编码及串/并转换，得到Ik、Qk两路数据。

10、波形选择地址生成器是根据接受到的数据（Ik或Qk）输出波形选择的地址。EEPROM（各种波形数据存储在其中）根据CPLD输出的地址来输出相应的数据，然后通过DA转换器得到我们需要的基带波形，最后通过乘法器调制，运放求和就得到了我们需要的MSK调制信号。MSK基带波形只有两种波形组成，见图2.4所示：图2.4 MSK成形信号在MSK调制中，成形信号取出原理为：由于成形信号只有两种波形选择，因此当前数据取出的成形信号只与它的前一位数据有关。如果当前数据与前一位数据相同，输出的成形信号就相反（如果前一数据对应波形1，那么当前数据对应波形2）；如果当前数据与前一位数据相反，输出的成形信号就相同（如果前

11、一数据对应波形1，那么当前数据仍对应波形1）。 2.2 MSK解调原理MSK信号的解调与FSK信号相似，可以采用相干解调，也可以采用非相干解调方式。本实验模块中采用一种相干解调的方式。已知： （公式 2.5)把该信号进行正交解调可得到：Ik路 （公式 2.6)=+Qk路 （公式 2.7) =+我们需要的是、两路信号，所以必须将其它频率成份、通过低通滤波器滤除掉，然后对、采样即可还原成、两路信号。根据上面描述可构成一种MSK解调器，其方框图如图2.5所示： 图2.5 MSK解调原理框图 将得到的MSK调制信号正交解调，通过低通滤波器得到基带成形信号，并对由此得到的基带信号的波形进行电平比较得到数

12、据，再将此数据经过CPLD的数字处理，就可解调得到NRZ码。在实际系统中，相干载波是通过载波同步获取的，相干载波的频率和相位只有和调制端载波相同时，才能完成相干解调。由于载波同步不是本实验的研究内容，因此在本模块中的相干载波是直接从调制端引入，因此解调器中的载波与调制器中的载波同频同相。载波同步的实验可在本实验箱的CDMA系统中实现。2.3 MSK的抗噪声性 数字通信系统中，信号的传输过程会受到各种干扰，从而影响对信号的恢复.通信系统的抗噪声性能是指系统克服加性噪声影响的能力。衡量数字通信系统抗噪声性能的重要指标是误码率. 分析二进制数字调制系统的抗噪声性能，得出误码率与信噪比之间的数学关系。

13、在二进制数字调制系统抗噪声性能分析中，假设信道特性是恒参信道，在信号的频带范围内其具有理想矩形的传输特性(可取传输系数为K)。 噪声为等效加性高斯白噪声，其均值为零，方差为25。 对2FSK信号解调同样可以采用同步检测法和包络检波法，下面分对同步检测法的系统性能进行分析。2FSK信号采用同步检测法时在码元时间宽度Ts区间，发送端产生的2FSK信号可表示为 （公式 2.8) 发送“1”符号是为, 发送“0”符号是为. 其中: （公式 2.9) （公式 2.10) 式中，和分别为发送“1”符号和“0”符号的载波角频率，Ts为码元时间间隔。在(0, Ts)时间间隔，信道输出合成波形为 （公式 2.1

14、1) 式中，为加性高斯白噪声，其均值为零，方差为。解调器采用两个带通滤波器来区分中心频率分别为和的信号。中心频率为的带通滤波器只允许中心频率为的信号频谱成分通过，而滤除中心频率为的信号频谱成分；中心频率为的带通滤波器只允许中心频率为2的信号频谱成分通过，而滤除中心频率为的信号频谱成分。 Pe= （公式 2.12) 在大信噪比条件下， （公式 2.13) 对于MSK通信系统设信道特性为恒参信道，噪声为加斯白噪声，MSK解调器输入、信号与噪声的合波为： （公式 2.14)式中 是均值为0，方差为 的窄带高斯噪声。经过相乘、低通滤波器滤波和抽样，在 时刻I支路的样值为在 时刻Q支路的样值为： 式中

15、和 分别为 和 在取样时刻的 样本值。在I支路和Q支路数据等概率的情况下，支路的误码率为： （公式 2.15) 式中， 为信噪比。经过交替门输出和差分译码后，总比特率为： 2.4 MSK功率谱密度 对相位不连续的二进制移频键控信号，可以看成由两个不同载波的二进制振幅键控信号的叠加，其中一个频率为，另一个频率为。因此，相位不连续的二进制移频键控信号的功率谱密度可以近似表示成两个不同载波的二进制振幅键控信号功率谱密度的叠加。根据二进制振幅键控信号的功率谱密度，我们可以得到二进制移频键控信号的功率谱密度为 （公式 2.16) 令概率， 将二进制数字基带信号的功率谱密度公式代入式可得 （公式 2.17

16、) 相位不连续的二进制移频键控信号的功率谱由离散谱和连续谱所组成 离散谱位于两个载频和处；连续谱由两个中心位于和处的双边谱叠加形成；若两个载波频差小于，则连续谱在处出现单峰；若载频差大于，则连续谱出现双峰。若以二进制移频键控信号功率谱第一个零点之间的频率间隔计算二进制移频键控信号的带宽，则该二进制移频键控信号的带宽B2FSK为 ： （公式 2.18）其中 图2.6相位不连续2FSK信号的功率谱示意图下面我们简要讨论一下MSK信号的功率谱。由可以得到MSK 信号的功率谱密度 （公式 2.19) 与2PSK相比，MSK信号的功率谱更加紧凑，其第一个零点出现在处，而2PSK的第一个零点出现在处。这表

17、明，MSK信号功率谱的主瓣所占的频带宽度比2PSK信号的窄；当(-)时，MSK的功率谱衰减速率更为迅速。它要比2PSK的衰减速率快得多，因此对邻道的干扰也较小。因此MSK信号抗干扰性能要优于2psk。 图2.7 MSK与2PSK信号的归一化功率谱3 MATLAB仿真设计与系统分析3.1 MATLAB简介在MATLAB通信工具箱中有SLMULINK仿真模块和MATLAB函数，形成一个运算函数和仿真模块的集合体，用来进行通信领域的研究、开发、系统设计和仿真。通信工具箱中的模块可供直接使用，并允许修改，使用起来十分方便，因而完全可以满足使用者设计和运算的需要。 MATLAB通信工具箱中的系统仿真，分

18、为用Simulink模块框图进行仿真和用MATLAB函数进行的仿真两种。在用SIMULINK模块框图的仿真中，每个模块，在每个时间步长上执行一次，就是说，所有的模块在每个时间步长上同时执行。这种仿真被称为时间流的仿真。而在用MATLAB函数的仿真中，函数按照数据流的顺序依次执行，意味着所处理的数据，首先要经过一个运算阶段，然后再激活下一个阶段，这种仿真被称为数据流仿真。某些特定的应用会要求采用两种仿真方式中的一种，但无论是哪种，仿真的结果是相同而且很方便。近几年来，在学术界和工业领域，M函数已经成为动态系统建模和仿真领域中应用最为广泛的软件之一。M函数可以很方便地创建和维护一个完整地模块，评估

19、不同地算法和结构，并验证系统的性能。由于M函数是采用模块组合方式来建模，从而可以使得用户能够快速、准确地创建动态系统的计算机仿真模型，特别是对复杂的不确定非线性系统，更为方便。M函数模型可以用来模拟线性和非线性、连续和离散或者两者的混合系统，也就是说它可以用来模拟几乎所有可能遇到动态系统。另外M函数还提供一套图形动画的处理方法，使用户可以方便的观察到仿真的整个过程。M函数没有单独的语言，但是它提供了S函数规则。所谓的S函数可以是一个M函数文件、FORTRAN程序、C或C+语言程序等,通过特殊的语法规则使之能够被M函数模型或模块调用。S函数使M函数更加充实、完备，具有更强的处理能力。同Matla

20、b一样，M函数也不是封闭的,他允许用户可以很方便的定制自己的模块和模块库。同时M函数也同样有比较完整的帮助系统，使用户可以随时找到对应模块的说明，便于应用。综上所述，M函数就是一种开放性的，用来模拟线性或非线性的以及连续或离散的或者两者混合的动态系统的强有力的系统级仿真工具。目前，随着软件的升级换代，在软硬件的接口方面有了长足的进步，使用Simulink可以很方便地进行实时的信号控制和处理、信息通信以及DSP的处理。世界上许多知名的大公司已经使用M函数作为他们产品设计和开发的强有力工。.3.2 MSK信号的调制因为MSK信号可以用两个正交的分量表示： （k-1）TstkTs （公式 3.1)

21、式中：右端第1项称作同相分量，其载波为；第2项称作正交分量，其载波为。 在仿真时先设置输入信号的参数：载波频率调制后,数字信号速率,模拟信号采样率。之后产生要调制的随机数字信号x;然后用dmod函数进行调制，产生调制信号。dmod函数是MATLAB中一个常用来进行信号调制的函数，它后面的参数包括被调制信号，载波信号的频率，输出信号的速率以及采样速率和所进行调制的函数。部分程序代码和仿真图如下所示：M=2; %定义输入随机信号的参数Fc=20; %载波信号的频率Fc为20；Fd=10; %输出信号的速率Fd为10；Fs=200; %采样速率Fs为200； x=randint(1,10000);

22、%产生随机的输入序列y=dmod(x,Fc,Fd,Fs,msk); %对用dnod函数对输入的随机序列进行调制3.3 MSK解调实现由于MSK信号是最小二进制FSK信号，所以它可以采用解调FSK信号的相干法和非相干法解调。在进行程序仿真时，用ddemod函数进行调制信号的解调。同时画出解调前后的时域与频域的波形图。在画频域的波形图时先对已调信号与解调信号进行DTFT变换，之后画出相应的波形。ddemod函数是与dmod函数相对应，用来对已调信号进行解调的。其后面的参数与dmod函数后的一模一样。部分程序与仿真波形图如下所示：z=ddemod(y,Fc,Fd,Fs,msk); %对调制后的MSK

23、信号进行解调 在解调后得到了时域的调制解调波形，然后再将时域的信号通过DTFT变换转变到频域，从而得到频域的调制解调图形n1=1:1:length(y);w=0:1:1000*pi/500; n2=1:1:10000;w=0:1:1000*pi/500; %定义频域参数 n1为y的长度，n2为输入原始信 号的长度Y=y*exp(-j*n1*w); %对已调信号进行DTFT变换Z=z*exp(-j*n2*w); %对解调信号进行DTFT变换magY=abs(Y);magZ=abs(Z);3.4叠加噪声的MSK系统由于信号在传输的过程中难免要受到外来噪声的影响，所以在进行通信仿真时务必要在理想的模

24、拟通信系统中加入高斯白噪声对系统进行影响，以此来判断一个通信系统抗噪声性能的好坏。在这次课程设计中，我利用awagn函数对已调信号加入信噪比为0.001的高斯白噪声。然后对加入了噪声的已调信号进行解调，同时画出其时域与频域的波形，将其与没有加噪声的调制信号进行比较，不难发现：加入噪声后对已调信号的影响很大，但对解调信号而言，噪声对解调后的影响还是有一定限度的。这说明MSK解调系统的抗噪声性能很强。仿真的部分程序与相应的仿真波形图如下所示：y1=awgn(y,0.001); %对已调信号加入高斯白噪声z1=ddemod(y,Fc,Fd,Fs,msk); %对加入高斯白噪声的调制信号进行解调num

25、1,ratio1=symerr(x,z2) %对加入噪声后解调的信号进行误码率的测定 再对信号转换到频域进行调制与解调对比分析，其部分程序如下：Y1=y1*exp(-j*n1*w); %对加噪声后的已调信号进行DTFT变换Z1=z1*exp(-j*n2*w); %对加噪声后的解调信号进行DTFT变换magY1=abs(Y1);magZ1=abs(Z1); 下一步通过画图函数显示出时域与频域的图像，进行比较，通过symerr函数对误码率进行判断，方便分析。4 MSK系统分析 4.1 MSK信号的时域调制解调分析 运行程序得到时域的调制解调图，如下图：图4.1 调制前的信号图4.2 调制后的信号

26、由调制波形图可以看出MSK的调制信号特性与2FSK调制信号的特性非常的相似，即：当输入信号为“1”时，调制后的波形比输入信号为“0”时的波形要密。同时MSK信号的包络是恒定的，相位则是连续的。带宽相对于一般的2FSK信号要小，而且正交。图4.3 MSK信号解调前的频域波形图4.4 MSK信号解调后的频域波形 由解调后的频域波形可以看出MSK信号的稳定性很好，说明MSK信号的抗噪声性能很强。 4.2 MSK频域分析 由程序运行得到频域调制前后的波形图，如下： 图4.5 调制前解调后的频域波形 图4.6 调制前解调后的频域波形通过调制前和解调后的频域波形比较，发现解调非常的成功，波形被完整无误的输

27、出来。 4.3 加噪声系统分析 对信号叠加不同信噪比噪声的噪声，其输出结果与对比如下图：图4.7信噪比为0.001的 MSK信号加噪声解调前的时域波形图4.8 信噪比为0.001的 MSK信号加噪声解调后的时域波形图4.9信噪比为1的MSK信号加噪声解调前的时域波形图4.10 信噪比为1的MSK信号加噪声解调后的时域波形图4.11 信噪比为2的 MSK信号加噪声解调前的频域波形图4.12 信噪比为2的 MSK信号加噪声解调后的频域波形在进行编译与运行后可得到频域的信号图，可以看出，加入噪声后的MSK系统依然十分稳定，虽然有些频率干扰，但对主信号的干扰效果并不显著，依然有十分良好的传输环境，即可

28、以看出，MSK通信系统的抗噪声性能还是十分的良好。 4.4误码率分析通过symerr函数对解调信号的误码率进行测定，得到的误码率为0.0024.因为输入的随机序列长度为10000，错了24个。说明MSK信号的抗噪声性能很强。当将信噪比变为1时，输出的误码率很快发生了变化,由原来的0.0024变为了6.0000e-004.图形也同时发生了变化。通过图形的比较，可以发现噪声对调制信号的影响很大，但通过解调后影响相对减小了很多。在多次改变信噪比后，将信噪比与误码率联合做了如图3-8所式的波形图。图4.13 信噪比与产生误码率曲线图通过图形可知，随着系统信噪比的增加，解调后输出的误码率随之减小，当信噪

29、比增加到一定程度时，误码率变为0。说明噪声对信号的影响在一定程度上是可以减小，但很难消除，因为每一个系统都不可能是理想系统，总会受到外界的干扰。 5 总结本次课程设计是频移键控通信系统调制解调的理论基础上，通过MATLAB软件simulink工具箱完成对整个通信系统的仿真设计。通过仿真我们可以不用通过硬件实验就可以得出实验结果。通过仿真结果衡量方案的可行性，从中选择最合理的系统配置和参数设置，然后再应用于实际系统中。通过仿真，可以提高研究开发工作的效率，发现系统中潜在的问题，优化系统整体性能。而且还可以大大的节约研究经费，这对科学技术的发展是很重要的。通过本次仿真结果可以使我们能很直面的了解到

30、msk通信系统性能。通信系统性能评价主要从：1.抗噪声性能（从误码率体现）。2.调制信号频带宽度。3.信道特性变化的敏感性。来衡量。从上面的仿真结果和比较可以得到结论，msk数字通信系统与fsk数字通信系统带宽相差不大.采用fsk调制的系统抗噪声性能相对要差一些两种系统工作过程复杂度相差不大，所以目前在噪声严重的环境下msk系统应用得更为广泛。当信道中的噪声比较小时（比如最大幅度是1），两种系统都能很好的传输数据；当信道中加入随机噪最大幅度为3时，msk通信系统基本没有误码，可见MSK数字通信系统有着很强的抗干扰能力，而采用fsk调制的系统就明显有很大的误码率，可见msk数字通信系统对信道变化

31、敏感程度不如采用msk调制的系统。可见msk调制解调数字通信系统更适合应用于信道衰落严重的通信系统。通过实验仿真我们可以设计出合理的系统参数。仿真技术作为一门综合性科学，将随着其相关领域技术的深入发展，继续向纵深快速发展，同时将扩大其综合应用的领域，发挥更大的作用。但是，作为一门综合性技术学科，仿真技术还有许多理论及技术问题需要继续进行深入的研究探讨。我国应大力开展仿真技术的理论研究和技术应用研究，尽快缩短与先进发达国家在技术上的差距。仿真技术的发展，必将推动我国科学技术水平的进一步提高。 致谢 一个月的课程设计很快就过去了，在这期间我学习和掌握了使用MATLAB对数字通信系统的仿真。回想起这

32、期间的经历，感慨万千。课程设计是对的知识的一次总结、一次检验。同时也是锻炼和进一步提高自己自学能力的一次宝贵机会。在做本次设计的初期阶段，由于自己在很的多知识点上有所遗忘，因此设计的进度相当缓慢，特别是在原理分析和M文件仿真上更是裹足不前。我重点复习并新学习了相关章节的知识,并很快理解了要点，进度由此快速展开。设计后期阶段，在设计思路上卡了壳。老师热心的帮我讲解帮我理顺了设计思路，使我得以进行软件调试从而验证了设计的正确可行性并按时完成设计，在此特别感谢老师的指导，此外还要感谢学校教研室的老师们的关心和培养，也感谢帮助和支持我的同学。在此表示衷心的感谢！由于理论知识有限、文中难免出现不足之处，

33、敬请各位老师批评指正。参考文献：1张志勇等.精通matlab6.5J.北京:北京航空航天大学出版社,2003.3.2John G.Proakis Masoud Salehi Gerhard Bauch.著 刘树棠译comtemporary communication systems using matlab and simulink.北京:电子工业出版社,2005.4.3樊昌信,张甫詡、徐炳祥、吴成柯 .通信原理M .北京:国防工业出版社,2001.5.4王立宁,乐光新,詹菲等.MATLAB与通信仿真M .北京:人民邮电出版社, 1999 5 陈怀琛，吴大正，高西金 matlab 及在信息课程

34、中的应用J北京：电子工业出版社,20066 楼天顺、刘晓东、李博涵基于matlab7.x的系统分析与设计M.西安：西安电子科技大学出版社, 2005.5 7张肃文主编 .高频电子线路M .北京：高等教育出版社, 2004.11 8程佩青编著 .数字信号处理教程M .北京：清华大学出版社, 2006.3 9张辉、曹丽娜、王勇 编著 .通信原理辅导M 西安.西安电子科技大学出版社, 2004.1110徐明远 .MATLAB仿真在通信工程中的应用J. 西安. 西安电子科技大学出版社, 2005.6 11田丽华等.编码理论J.西安 西安电子科技大学出版社, 2003.8附录:MSK信号调制解调的源程序代码M=2; %定义输入随机信号的参数Fc=20;Fd=10;Fs=200;t=0:500;x=randint(1,100); %产生随机的输入序列y=dmod(x,Fc,Fd,Fs,msk); %对输入的随机序列进行调制figure(1); %第一张图画调制前后的波形subplot(2,1,1),stem(x),title(调制前的输入波形);subplot(2,1,2), (t,y(1:501),title(调制后的输出波形); axis(0 500 -1 1)z=ddemod(y,Fc,