毕业设计（论文）平原GMSK仿真设计

1、目 录第一章 matlab简介11.1matlab概念11.2基本功能11.3主要应用2第二章 数字通信系统32.1系统组成及模型32.2道编码的基本思想和基本方式32.2.1信道编码的基本思想32.2.2信道编码的基本方式32.3常用的信道编码32.3.1线性分组码32.3.2卷积码4第三章 调制解调技术53.1调制解调的概念53.2数字调制的基本类型53.2.1线性调制技术53.2.2恒定包络调制技术53.3几种常见的数字调制方式63.3.1二进制频移键控（2fsk）63.3.2最小移频键控（msk）63.3.3二进制移相键控（bpsk）73.3.4多进制频移键控mfsk83.3.5高斯滤

2、波的最小移频键控gmsk9第四章 信道的仿真设计与研究104.1课题分析和设计104.2数字通信系统程序设计104.2.1仿真程序设计104.2.2仿真程序设计概述114.3 仿真框图设计144.3.1设计原理144.3.2平原中gmsk调制的仿真框图设计的过程15第五章 仿真结果及结论225.1数字通信系统仿真框图比较结果与分析225.1.1数字通信系统框图仿真结果225.1.2数字通信系统框图仿真结果分析225.2数字通信系统仿真程序比较结果与分析225.2.1数字通信系统程序仿真结果225.2.2数字通信系统程序仿真结果分析22致 谢22参考文献24附录 matlab通讯系统工具箱库函数

3、及仿真模块28第一章 matlab简介1.1matlab概念matlab是矩阵实验室（matrix laboratory）的简称，是美国mathworks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括matlab和simulink两大部分。1.2基本功能matlab是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提

4、供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如c、fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 matlab和mathematica、maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。matlab可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 图1-1 matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 matlab的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用matlab来解算问题要比用

5、c，fortran等语言完成相同的事情简捷得多，并且matlab也吸收了像maple等软件的优点,使matlab成为一个强大的数学软件。在新的版本中也加入了对c，fortran，c+ ，java的支持。可以直接调用,用户也可以将自己编写的实用程序导入到matlab函数库中方便自己以后调用，此外许多的matlab爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。1.3主要应用matlab 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理技术 数字信号处理 技术 通讯系统设计与仿真 图1-2 matlab在通讯系统设计与仿真的

6、应用 财务与金融工程 matlab 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 matlab 函数集）扩展了 matlab 环境，以解决这些应用领域内特定类型的问题。第二章 数字通信系统2.1系统组成及模型最简单的数字通信系统模型由信源、信道和信宿三个基本部分组成。实际的数字通信系统模型要比简单的数字通信系统模型复杂得多。一般的数通信系统模型如图1-1所示。图2-1 数字通信系统模型2.2道编码的基本思想和基本方式2.2.1信道编码的基本思想在发端被传输的信息序列上附加一些多余的检验码元，这些监

7、督码元与信息码元间以某种特定的规则相互关联。接收端按照既定的规则检验信息码元与监督码元间的关系，一但传输出错，则两者间的关系将会受到破坏，从而可以发现错误，乃至纠正错误。2.2.2信道编码的基本方式前向纠错（fec）、检错重传（arq）、混合差错控制（hec）和信息反馈（irq）。2.3常用的信道编码2.3.1线性分组码在分组码中，编码后的码元序列每n位为一组，其中k是信息码元，r（r=n-k）个是附加的监督码元。如果信息元与监督码元之间呈线性关系，则为线性分组码。奇偶校验码和汉明码是典型的线性分组码。循环码是线性分组码，既可纠错又可检错。特点是任一码组的每一次循环移位得到的是码中的另一码组。

8、bch码是一类纠正多个随机错误的循环码。2.3.2卷积码卷积码是另一类信道编码，它也是把k位信息编成n（n大于等于k）位，但k和n都很小，适于串行传输。特点是编码后的n个码元不仅与当前段的k位信息有关，而且与前n-1段的信息有关，编码过程中相互关联的码元有nn个。纠错能力随n的增加而增加，而差错率随n的增加而指数下降。编码过程可以看成是输入信息序列与由移位寄存器和模2加连接所决定的的另一序列的卷积，因此称为卷积没码。n称为约束长度，m=n-1称为编码存储。第三章 调制解调技术3.1调制解调的概念调制解调是数字通信系统的重要组成部分。调制解调的目的是使已调信号具有高的频谱利用率和较强的抗干扰和抗

9、衰落的能力。下面就对数字调制解调技术进行介绍。所谓调制，就是用调制信号（基带）去控制或改变载波的一个或几个参数，使调制后的信号（已调信号）含有原来调制信号的全部信息，但载波的某些参数按调制信号的规律变化。调制的目的是使传输的信号与信道相匹配，从而有效传输信号。解调是调制的逆过程，它是从已调信号中恢复出原来调制信号的过程。从广义上讲，调制与解调属于信道编/译码范畴，但调制与解调的目的是实现载波传输，而信道编/译码的主要目的是实现差错控制3.2数字调制的基本类型数字信号调制的基本类型分振幅键控（ask）、频移键控（fsk）和相移键控（psk）。此外，还有许多由基本调制类型改进或综合而获得的新型调制

10、技术。表1-1给出调制方式以及主要用途。表3-1 调制方式及用途调制方式主要用途振幅键控ask数据传输频率键控fsk数据传输相位键控psk、dpsk、qpsk等数据传输、数字微波、空间通信高效数字调制qam、msk等（提高频带利用率）数字微波、空间通信3.2.1线性调制技术这里包括psk、qpsk、dqpsk、ok-qpsk和多电平的psk等。所谓的“线性”，是指这类调制技术要求设备从频率变换到放大和发射的过程中保持充分的线性。这种要求在制造移动设备中增加了难度和成本，但这种方式可获得较高的频谱利用率。3.2.2恒定包络调制技术主要包括msk、gmsk和tfm等。这类调制技术的优点是已调信号具

11、有相对窄的功率谱和对放大设备没有线性要求，不足之处是其频谱利用率通常低于线性调制技术。另一种获得迅速发展的数字调制技术是振幅和相位联合调制（qam）.3.3几种常见的数字调制方式3.3.1二进制频移键控（2fsk）调制原理：在二进制移频键控中，调制信号1或0，分别对应载波的两个频率f1 或 f2 。因此，其调制实现方法就是用输入的二进制信号去控制两个独立的载波发生器，如图1-2所示。图3-1 2fsk调制解调也有相干和非相干两种。更简单的方法是过零检测器法：根据移频键控的过零率的大小来检测已调信号的频率变化，其组成及各点波形如图1-3所示。图3-2 2fsk信号过零检测法解调波形示意图其中a是

12、已调信号，b是限幅后已调信号，c是经过微分的信号，d是整流后的已调信号，e是经过宽脉冲发生后的已调信号，f是经过低通的输出信号。3.3.2最小移频键控（msk）调制原理：msk是fsk的一种特殊情况，其特殊之处在于msk在相邻符号交界处相位保持连续，是一种连续相位fsk。若fsk看作非正交2fsk，则msk调制方法如图1-4。图3-3 msk非正交调制式解调原理：msk可用2fsk方法进行相干解调，并每隔tb时刻作出判决；也可用鉴频器的方法进行非相干解调。3.3.3二进制移相键控（bpsk）调制原理：载波的相位（通常为00或1800）随调制信号“1”或“0”而改变，这种调制称为二进制移频键控。

13、bpsk信号是双极性非归零码的双边带调制。bpsk调制有直接相乘方法和相位选择法两种，见图1-5。图3-4 bpsk调制方法解调原理：对于bpsk信号的解调必须采用相干解调的方法。由于bpsk解调器中的本地载波的相位有00、1800的模糊，通常采用在调制输入的数字基带信号中进行差分编码的方法来解决。这种方法称为二进制差分移相键控（2dpsk）。解调时利用延迟电路将其前一码元延迟一个码元时间ts作为参考相位，并与后一码元相乘，再进行低通滤波，最后经取样判决后恢复出原二进制码。差分相干解调法见图1-6。图3-5 bpsk差分相干解调法3.3.4多进制频移键控mfsk调制原理：在mfsk方式中，mf

14、sk信号常用频率选择法产生。调制器框图如图1-7所示，先将数据比特流经过一个串/并转换器变成多路并行数据，将这多路并行数据分别用不同的频率通过线性门进行选择，再将各线性门输出的数据进行叠加，输出的即为mfsk信号。图3-6 mfsk调制器框图解调原理：mfsk的解调方法常用非相干解调，如图1-8所示，将接受到的mfsk信号通过滤波器进行滤波，然后将滤波后的信号分别送到个带通滤波器进行滤波，再将输出信号进行包络检波，最后将这个信号进行抽样和判决，最后经过一个逻辑电路，将多进制信号转化为二进制信号即可。图3-7 mfsk信号非相干解调3.3.5高斯滤波的最小移频键控gmsk调制原理：gmsk信号是

15、通过在fm调制器前加入高斯低通滤波器（称为预调制滤波器）而产生的。原理图如图1-9所示。图3-8 gmsk信号产生原理图解调原理：gmsk信号的解调可以采用msk一样的正交相干解调电路。在相干解调中最为重要的是相干载波的提取，在移动通信的环境中是比较困难的，因而通常采用差分解调和鉴频器解调等非相干解调。第四章 信道的仿真设计与研究4.1课题分析和设计本章内容是仿真设计的重要组成部分。这一章全面、系统的介绍了仿真软件的整个设计过程。全章共分为：数字通信系统程序和框图仿真设计、数字通信系统界面设计和需要说明的问题三部分。仿真设计部分又分仿真程序设计和仿真框图设计两个方面介绍。在仿真设计和界面设计的

16、介绍中均给出了一般的操作方法和详尽的设计流程并付有图片。相信，通过本章的介绍，可以使您对设计过程有一个清楚明了的认识。本章中的所有使用的仿真程序、m函数和模块可参见附录。4.2数字通信系统程序设计4.2.1仿真程序设计图4-1 程序设计流程图4.2.2仿真程序设计概述（i）m文件简介以matlab程序代码所编写的文件通常以“m”为扩展名，所以这些文件又称为“m文件”（m-files）。m文件是一个文本文件，所以可以用各种文本编辑器去编辑它，但记得在储存文件时必须以文本模式储存。此外，matlab在win95/98/nt及mac的平台上，也提供了一个内建的m文件编辑器（m-file editor

17、）,可以单击matlab命令窗口的file/open下拉式菜单，来打开m文件编辑器；另一种更快的方法，则是在matlab命令窗口直接输入“edit filename.m”。m文件编辑器会用不同的颜色来区分注释、关键词、字符串及一般程序代码。（ii）程序仿真中相关matlab库函数（m函数）的介绍*randint功能：归一化分布随机整数、矩阵产生器。说明：out= randint(n,m)中 n代表矩阵的行数，m代表矩阵的列数。此命令可以产生一个 n*m均匀分布的二进制矩阵。*encode功能：差错控制编码（信道编码）说明：code=encode(msg,n,k,method,opt)中 msg

18、是信息， n 是码字长度，k是信息位长度，method注明编码方式，opt是有些编码方式需要的参数。此命令可以完成汉明码、线性分组码、cyclic码、bch码、r-s码、卷积码六种主要的差错控制编码。*decode功能：差错控制译码（信道译码）说明：msg=decode(code,n,k,method)中 code是指接收到的码字，n是码字长度，k是信息位长度，method注明译码的方式。此命令可以对接收到的码长为n,信息位为k的码字进行译码，恢复出原始的信息，译码方式必须和编码时采用的严格相同。同样具有六种主要的差错控制译码。*dmod功能：数字调制（通带信号）说明：modu=dmod(x,

19、fc,fd,fs,method,m)中x是输入信号，fc是输入信号x的载波频率，fd是输入信号x的抽样频率， fs是仿真是的抽样频率，method注明了调制方式，m是进制数，此命令可以对数字信号x进行调制。method对应的调制方法有ask,fsk,msk,psk,qask。注意fs/fc=2, fs/fd必须是整数。*ddemod功能：数字解调（通带信号）说明：x=dmod(y,fc,fd,fs,method,m) 中y是经过dmod 调制的信号，fc是输入信号x的载波频率，fd是输入信号x的抽样频率， fs是仿真是的抽样频率，method注明了解调方式，m是进制数，此命令可以把用dmod调

20、制的信号y解调出。method也有五种调制方法。*biterr功能：计算误比特数、误比特率说明：rate= biterr(x,y,k)中x 和y都是矩阵，k是指定x,y矩阵中每个元素都用k位二进制数表示。此命令可以计算两个输入矩阵x,y的不同元素的个数和之比。x,y中的元素必须是非负的二进制数式整数。此命令可以用来计算误比特率，并且，x,y 中的一个为接收数据，一个为正确的参考数据。*fopen功能:打开文件说明：fid=fopen(filename,permission)其中filename表示要读写的文件名称，permission则表示要对文件进行的处理方法，处理方法可以是下列任一字符串：

21、 r:只读文件 w: 只写文件 a: 只加入文件 a+: 可读取及加入文件*fread 功能: 可从文件中读取二进制数据，将每一个字节看成一个整数，并将结果以一矩阵返回。说明：fid=fopen(filename,r);% 读取filename文件 data=fread(fid,n m) %fid 是要读取文件的标志码（由fopen 产生），一般是一个非负的整数。若返回的文件标志码为-1，则代表该文件无法打开，n m是指把读出的文件转化成n行m列的矩阵。此命令可读以取filename文件并且将其转化为n行m列的矩阵，然后返回到data。*de2bi功能:十二进制转化说明：msg1=de2bi(

22、data)把data中的数据进行十二进制转换。*reshape功能：可将一矩阵保持其元素个数不变，变成所须矩阵形式。说明：msg=reshape(msg1,n,m)n是现矩阵的行数，m是现矩阵的列数，msgl是需变换的矩阵，注意n*m的值必须与msg1矩阵的元素个数相同，否则就会发生错误。此命令可将msg1变换成n行m列的矩阵msg。4.2.3仿真程序设计过程（i）设计gmsk调制程序如下：clear all; （1）clc; （2）m=2;n=15;k=11; （3）fc=1000;fd=100;fs=2000; （4）msg=randint(k*100,1); （5）code=encode

23、(msg, n, k, cyclic ); （6）modu=dmod(code, fc, fd, fs, gmsk,m); （7）std_value=0.1 ; （8）modu_noise =modu+randn(length(modu), 1 )*std_value; （9）demo =ddemod(modu_noise, fc, fd, fs, gmsk , m); （10）msg_r =decode(demo, n, k, cyclic ); （11）rate=biterr(msg, msg_r, m) （12）下面对以上程序设计进行详细说明：程序第三行先定义进制m ,码长n和信息长度k

24、； 程序第四行定义载波频率fc ,数字传输频率fd和模拟频率(仿真采样频率)fs； 程序第五行用randint命令来产生一个1100*1包含待传输信息的随机整数矩阵，作为信号源；程序第六行使用cyclic差错控制编码方式对信道进行编码，码长为 n,信息长度为 k；程序第七行对已编码的信号源进行二进制fsk调制，fc =1000， fd =100 ， fs =2000；程序第八行和第九行在信道加入加性高斯白噪声，假定加性高斯白噪声的 std（噪声系数）为 0 .1，从而获得加入加性高斯白噪声的已调信道modu_noise；程序第十行是fsk解调过程；fc =1000， fd =100 ， fs

25、=2000；程序第十一行进行cyclic解码，码长为 n,信息长度为 k；程序第十二行对接收信号进行误码计算。（ii）fsk调制与msk 调制/文本传输仿真程序设计比较：a. fsk调制与msk 调制仿真程序设计的过程和思路大致相同，但在信道调制方法上不同，前者使用的是fsk调制，后者使用的是msk调制。并且二者在相同的信源和信道的情况下误码率不同。b.fsk调制仿真程序设计与文本传输仿真程序设计的主要不同在于：第一，前者的信源是一个随机的信源。后者是通过对一个file.txt文本进行打开（第1行）、读（第2行）、十-二进制转换（第3行）、重新构造矩阵形式（第4行）的操作，来生成一个18行1列

26、的信号源。程序如下：fid=fopen(file.txt,r); （1）data=fread(fid,k*100,1); （2）msg1=de2bi(data); （3）msg=reshape(msg1,18,1); （4）第二，在信号调制上文本传输仿真程序使用的是qask。其它程序设计过程与gmsk调制仿真程序设计过程相同。（iii）仿真界面中下拉式菜单程序设计下面详细介绍仿真界面中地形选择中bpsk下拉式菜单程序设计过程：在kuangtu.m文件第98行后，插入以下程序：h2 = findobj(0, tag, popupmenu1); （1）value=get(h2,value); （2

27、）switch get(h2, value) （3）case 1 （4）bpskpy; （5）case 2 （6）bpskcs; （7）case 3 （8）bpskql; （9）otherwise （10）disp(unknown option); （11）end （12）程序第一行以系统默认的方式寻找该弹出框的tag属性值，作为一个对象（即句柄操作）赋予变量h2，建立h2对该控件的控制关联；程序第二行的value值代表下拉条的对应值，并通过get函数将该值动态赋值给h2，至此h2这一变量获得对该控件的控制权；程序第四行到第十一行是进行对象选择。当value值为1时，系统调用bpskpy.md

28、l; 当value值为2时，系统调用bpskcs.mdl; 当value值为3时，系统调用bpskql.mdl。mfsk下拉式菜单和gmsk下拉式菜单的程序设计与bpsk下拉式菜单的程序设计相似。4.3 仿真框图设计4.3.1设计原理由于在相同的调制方式下不同的地形首先反映在信号源模块参数的不同，即信号源中“0”（无用信息）的概率的大小。例如：平原地形信号源中“0”（无用信息）的概率为0.4,城市地形信号源中“0”（无用信息）的概率为0.4,丘陵地形信号源中“0”（无用信息）的概率为0.2；其次，根据具体情况，不同的地形也就意味着不同的信道，因为信道中所受噪声干扰的种类不同，信噪比不同，这样我

29、们根据对信道参数的设置就可以仿真现实中不同的地形。例如：平原地形的信道噪声为瑞利噪声，信噪比为4，城市地形的信道噪声为高斯白噪声，信噪比为14，丘陵地形的信道噪声为高斯白噪声，信噪比为4。4.3.2平原中gmsk调制的仿真框图设计的过程仿真框图如下（1） （2） （3） （4） （8） （7） （6） （5）图4-2 平原中gmsk调制的仿真框图现在对仿真框图中各个模块的参数设置进行说明：（1） 信号源参数设置说明双击模块出现如图4-3所示参数设置框，更改第一个设置框“0信号出现概率”为0.4，更改最后一个设置框“每帧采样码元个数”为8。其它设置采用系统默认值。图4-3 信号源参数设置（2）

30、crc编码参数设置说明crc编码参数设置如图4-4所示，设置采用系统默认值。图4-4 crc编码参数设置（3）gmsk调制参数设置说明将输入类型改为bit，其它设置采用系统默认值。图4-5 gmsk调制参数设置（3） 信道调制参数设置说明将es/no改为，且initial seed项改成1237，表示开始的种子数为1237。其它设置采用系统默认值。图4-6信道参数设置（4） gmsk解调参数设置说明该参数与gmsk调制参数设置相对应。图4-7 gmsk解调参数设置（5） crc译码参数设置说明crc译码参数设置与crc编码参数设置相对应。图4-8 crc译码参数设置（6） 误码率统计参数设置说

31、明输出端口“output data”改为port,这样就增加一个输出端口。其它设置采用系统默认值。图4-9 误码率统计参数设置（7） 误码仪参数设置说明所有设置采用系统默认值。图4-10 误码仪参数设置其它调制方式的框图设计与其类似，主要参数差别见表4-1表4-1 框图设计主要参数差别地 形信号源中“0”的概率噪声类型及信噪比bpsk平原mfsk平原gmsk平原0.4瑞利噪声信噪比es/no=4bpsk城市mfsk城市gmsk城市0.4高斯白噪声信噪比es/no=4bpsk丘陵mfsk丘陵gmsk丘陵0.2高斯白噪声信噪比es/no=14第五章 仿真结果及结论5.1数字通信系统仿真框图比较结果

32、与分析5.1.1数字通信系统框图仿真结果（结果格式为“调制方式地形：误码率-出错码元数-码元总数”）平原地形： bpskpy:0.4681-2344-5008 mfskpy:0.4155-4.073e+007-9.804e+007 gmskpy:0.3942-1974-5008城市地形： bpskcs:0.4145-2076-5008 mfskcs:0.3942-1974-5008 gmskcs:0.3942-1974-5008丘陵地形： bpskql:0.8055-4034-5008 mfskql:0.5643-2826-5008gmskql:0.1945-2826-50085.1.2数字通

33、信系统框图仿真结果分析通过对以上仿真结果进行比较，分析得出：gmsk调制方式在三种地形情况下的误码率均为最底、失真最小。在相同条件下，较另外两种调制方式的性能优异。从理论实践方面，说明了现实生活中的数字通信系统选择gmsk调制方式的正确性。5.2数字通信系统仿真程序比较结果与分析5.2.1数字通信系统程序仿真结果为0程序运行结果为0 ，即说明误码率为0。5.2.2数字通信系统程序仿真结果分析结果说明了编制的仿真调制程序是正确可执行的。致 谢这次毕业论文能够得以顺利完成，是所有曾经指导过我的老师、帮助过我的同学和一直支持着我的家人对我的教诲、帮助和鼓励的结果。我要在这里对他们表示深深的谢意！首先

34、，要特别感谢我的指导老师。老师在我毕业论文的撰写过程中，给我提供了极大的帮助和指导。从开始选题到中期修正，再到最终定稿，老师给我提供了许多宝贵建议。老师渊博的专业知识，严谨的治学态度，精益求精的工作作风，诲人不倦的高尚师德，朴实无华、平易近人的人格魅力对我影响深远。使我树立了远大的学术目标、掌握了基本的研究方法和许多待人接物和为人处世的道理。 下来，要感谢江西理工大学南昌校区所有曾经给我们08应用电子技术班任课的老师，老师们教会我的不仅仅是专业知识，更多的是对待学习、对待生活的态度。 其次，感谢我的父母亲，你们是我力量的源泉，只要有你们，不管面对什么样的困难，我都不会害怕，谢谢你们对我的支持与

35、鼓励！ 第三，感谢我的室友及其他好友，因为有你们的帮助，我的论文得以顺利完成。感谢你们，大学三年给我了那么多的帮助与鼓励，在我不开心的时候，总能让我开心起来。不会忘记，大学四年里我们一起度过的欢乐时光，那些开心的日子，总是那么令人难以忘怀。 最后对老师，同学和家人再次致以我最衷心的感谢！教导过我的老师，你们的人格魅力永记我心间；身边的同学和朋友，有你们，我的大学才算完整；寝室的好友，你们的天赋犹如上天恩赐，有了你们我的生活更加精彩！参考文献1 william h. tranter/ k. sam shanmugan/ theodore s. rappaport 美, 通信系统仿真原理与无线应用

36、：机械工业出版社，2005.2 刘敏，魏玲编著.matlab通信仿真与技术应用：国防工业出版社，2001.3 赵静，张瑾，高新科编著.基于matlab的通信系统仿真：北京航空航天大学出版社，2007.4 张智星等编著.matlab程序设计与应用.北京：清华大学出版社，2002.5 张森 张正亮等编著.matlab仿真技术与实力应用教程.北京：机械工业出版社，2004.6 樊昌信等编著.通信原理.北京：国防工业出版社，2003.7 张葛祥，李娜等编著.matlab仿真技术与应用.北京：清华大学出版社，2003.8 韦岗，季飞，傅娟编著.通信系统建模与仿真：北京航空航天大学出版社，2003.9 樊昌信等，通信原理北京：国防工业出版社，200310 张志涌等，精通matlab6.5版北京：北京航空航天大学出版社，20031