通信系统综合课程设计

1、 武汉理工大学通信系统课群综合训练与设计报告 课程设计任务书学生姓名： 专业班级： 指导教师： 工作单位： 题 目: 通信系统课群综合训练与设计1课程设计的目的通过课程设计，使学生加强对电子电路的理解，学会对电路分析计算以及设计。进一步提高分析解决实际问题的能力，通过完成综合设计型和创新性实验及训练，创造一个动脑动手独立开展电路实验的机会，锻炼分析解决电子电路问题的实际本领，实现由课本知识向实际能力的转化；加深对通信原理的理解，提高学生对现代通信系统的全面认识，增强学生的实践能力。2 课程设计要求要求：掌握以上各种电路与通信技术的基本原理，掌握实验的设计、电路调试与测量的方法。1.培养学生根据

2、需要选学参考书，查阅手册，图表和文献资料的自学能力，通过独立思考深入钻研有关问题，学会自己分析解决问题的方法。2.通过对实验电路的分析计算，了解简单实用电路的分析方法和工程设计方法。3.掌握示波器，频谱仪，失真度仪的正确使用方法，学会简单电路的实验调试和整机指标测试方法，提高动手能力。3 课程设计进度安排序号设 计 内 容所用时间1根据设计任务，分析电路原理，确定实验方案2天2根据实验条件进行电路的测试，并对结果进行分析7天3撰写课程设计报告1天合 计2周指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日 目 录摘 要IAbstractI1设计任务11.1设计目的11.2课程设计

3、要求11.3课程设计任务12实验原理分析22.1实验原理框图22.2 PCM原理介绍32.2.1 抽样(Sampling)32.2.2量化（quantizing）32.2.3 A压缩律：42.2.4 编码（Coding）92.2基带传输PST码92.3循环码102.3.1循环码介绍102.3.2 循环码编码原理122.3.3 循环码的纠错原理132.4二进制频移键控（2PSK）调制与解调原理152.4.1 数字调制技术152.4.2 二进制相移键控（2PSK）基本原理162.4.3 MATLAB实现2PSK调制与解调182.5衰落信道183 MATLAB软件仿真及结果193.1 MATLAB简

4、介193.2 发送端仿真结果与分析193.2.1仿真结果图193.2.2 实验结果分析244.小结25参考文献26附录：27 摘 要 通信系统的作用就是将信息从信源发送到一个或多个目的地。一个数字通信系统是利用数字信号来传递信息的通信系统，主要有信源编码与译码、信道编码与译码、数字调制与解调。具体过程是把信息源的模拟信号经过模数转换变成数子信号然后，再进行基带编码、信道编码、载波调制然后发送到信道中，然后造接收端进行对应的解调、译码、数模转换得到最终的接收信号。本设计结合PCM的抽样、量化、编码原理把模拟信号转换数字信号、汉明码编译码原理实现信道编译码、PST编译码原理用于实现信源编译码、PS

5、K调制解调原理用于实现数字调制与解调、衰落信道原理用于信道和噪声源的近似，利用MATLAB软件编程，完成了整个通信系统编码仿真分析。关键词：MATLAB软件编程、数字通信系统、PCM、PST码、循环码、PSK Abstract Communication system is the role of the information from the source sent to one or more of the destination. A digital communication system is to use digital signal to transfer informatio

6、n communication system, basically have source coding and decoding, channel coding and decoding, digital modulation and demodulation. The specific process is the source of the analog signal after a/d convert into number son signal and then, again carries on the baseband coding, channel coding, carrie

7、r modulation and then sent to the channel, then made on the receiving end corresponding demodulation, decoding, d/a conversion to get the final receiving signal. This design with PCM sampling and quantization and coding principle to analog signal into the digital signal, hamming code knitting decodi

8、ng theory to realize channel knitting decoding, PST knitting decoding principle used to implement the source knitting decoding, PSK modulation demodulation principle used to implement the digital modulation and demodulation, Fading channel and the principle for noise approximation, using MATLAB soft

9、ware programming, complete the whole communication system coding simulation analysis.Keywords: MATLAB software programming, digital communication system, PCM, PST, Cyclic code, PSKII1设计任务1.1设计目的通过课程设计，使学生加强对电子电路的理解，学会对电路分析计算以及设计。进一步提高分析解决实际问题的能力，通过完成综合设计型和创新性实验及训练，创造一个动脑动手独立开展电路实验的机会，锻炼分析解决电子电路问题的实际

10、本领，实现由课本知识向实际能力的转化；加深对通信原理的理解，提高学生对现代通信系统的全面认识，增强学生的实践能力。1.2课程设计要求要求：掌握以上各种电路与通信技术的基本原理，掌握实验的设计、电路调试与测量的方法。1.培养学生根据需要选学参考书，查阅手册，图表和文献资料的自学能力，通过独立思考深入钻研有关问题，学会自己分析解决问题的方法。2.通过对实验电路的分析计算，了解简单实用电路的分析方法和工程设计方法。3.掌握示波器，频谱仪，失真度仪的正确使用方法，学会简单电路的实验调试和整机指标测试方法，提高动手能力。1.3课程设计任务利用仿真软件（如Matlab或SystemView），或硬件实验系

11、统平台上设计完成一个典型的通信系统。学生要完成整个系统各环节以及整个系统的仿真，最终在接收端或者精确或者近似地再现输入（信源），计算失真度，并且分析原因。（信源为模拟信源，数字化方式为增量调制，基带码为AMI码，信道码为循环码，调制方式为ASK调制，信道类型为衰落信道）。完成整个系统各环节以及整个系统的仿真，最终在接收端或者精确或者近似地再现输入（信源），计算失真度，并且分析原因。信源：自己构造一时间函数，数字化方式:PCM，基带码:PST，信道码:循环码，调制方式:PSK，信道类型:衰落信道 ；解调，信道解码，基带解码，数模转换的过程与输入端对应。2实验原理分析2.1实验原理框图数字调制信道

12、编码基带编码模数转换模拟信源信道数字解调信道解码模拟信源数模转换基带解码图2.1通信系统的原理框图模数转换中有三个基本过程：抽样、量化、编码。根据编码的不同有几种调制方法，如脉冲振幅调制（PAM），脉冲编码调制（PCM），差分脉冲编码调制（DPCM），增量调制。基带传输的常用码型有AMI码，HDB3码，PST码，曼彻斯特码，密勒码，CMI码。信道编码的常用编码有汉明码，卷积码，循环码，BCH码。常见的数字调制方式有振幅键控（ASK），频移键控（FSK），相移键控（PSK）。常见的信道有加性高斯白噪声信道，多径衰落信道。根据题目要求，数字化方式为PCM调制，基带码为PST码，信道码为循环码，数字

13、调制方式为PSK调制，信道为衰落信道，则有以下原理框图循环编码PSK调制PST编码PCM调制模拟信源衰落信道PSK解调循环解码模拟信源PCM解调PST解码图2.2 实验原理框图2.2 PCM原理介绍 模拟信号数字化必须经过三个过程，即抽样、量化和编码，以实现话音数字化的脉冲编码调制（PCM，Pulse Coding Modulation）技术。2.2.1 抽样(Sampling) 图2.3 抽样原理图离散时间信号通常是有连续时间信号经周期采样得到的。完成采样功能的器件称为采样器，下图所示为采样器的示意图。图中Xa(t)表示模拟信号，Xa(nt)表示采样信号，T为采样周期，n=0，1，2，。一般

14、可以把采样器视为一个每隔T秒闭合一次的电子开关S。在理想情况下，开关闭合时间满足<<T。实际采样过程可视为脉冲调幅过程，Xa(t)为调制信号，被调脉冲载波p(t)是周期为T、脉宽为的周期脉冲串。当0时的理想采样情况是实际采样的一种科学的、本质的抽象，同时可使数学推导得到简化。下面主要讨论理想采样。2.2.2量化（quantizing）抽样信号虽然是时间轴上离散的信号，但仍然是模拟信号，其样值在一定的取值范围内，可有无限多个值。显然，对无限个样值一一给出数字码组来对应是不可能的。为了实现以数字码表示样值，必须采用“四舍五入”的方法把样值分级“取整”，分的级数越多，即量化级差或间隔越小

15、，量化噪声也越小。2.2.3 A压缩律：所谓A压缩律也就是压缩器具有如下特性的压缩律：其中，A为压缩系数；y为归一化的压缩器输出电压；x为归一化的压缩器输入电压。图画出了A为某一取值的归一化压缩特性。A律压缩特性是以原点奇对称的，为了简便，图中只给出了正半轴部分。图2.4 A压缩律特性上图中，x和y都在-1和+1之间，取量化级数为N(在y方向上从-1到+1被均匀划分为N个量化级)，则量化间隔为当N很大时，在每一量化级中压缩特性曲线可看作是直线，因此有式中，xi为第i个量化级间隔的中间值。因此(2.1)为了使量化信噪比不随信号x变化，也就是说在小信号时的量化信噪比不因x的减小而变小，即应使各量化

16、级间隔与x成线性关系，即则式（2.1）可写成(2.2)即其中k为比例常数。当量化级数很大时，可以将它看成连续曲线，因而式(2.2)成为线性微分方程解此微分方程(2.3)其中c为常数。为了满足归一化要求，当x=1时，y=1，代入式(2.3)可得故所得结果为即(2.4)如果压缩特性满足上式,就可获得理想的压缩效果，其量化信噪比和信号幅度无关。满足上式的曲线如下图所示，由于其没有通过坐标原点，所以还需要对它作一定的修改。图2.5 理想压缩特性曲线A律压缩特性就是对式(2.4)修改后的函数。在上图中，通过原点作理想压缩特性曲线的切线oc，将oc、cd作为实际的压缩特性。修改以后，必须用两个不同的方程来

17、描述这段曲线，以切点c为分界点，线段oc的方程:设切点c的坐标为(x1，y1)斜率为则由式(2.4)可得(2.5)所以线段oc的方程为所以当x=x1时，y1=1/k时，有因此有所以，切点坐标为 (exp-(k-1),1/k) ，令则将它代入式(2.5),就可得到以切点c为边界的段的方程为(2.6)因cd段的方程，满足式(3.4),所以由该式可得(2.7)由以上分析可见，经过修改以后的理想压缩特性与图5中所示的曲线近似，而式(2.6)式(2.7)和式(2.4)完全一样。13折线：实际中，A压缩律通常采用13折线来近似，13折线法如图2.4所示，图中先把轴的0,1区间分为8个不均匀段。图2.6折线

18、示意图其具体分法如下:a.将区间0，1一分为二，其中点为1/2,取区间1/2,1作为第八段;b.将剩下的区间0,1/2再一分为二，其中点为1/4,取区间1/4,1/2作为第七段;c.将剩下的区间0,1/4再一分为二，其中点为1/8,取区间1/8,1/4作为第六段;d.将剩下的区间0,1/8再一分为二，其中点为1/16,取区间1/16,1/8作为第五段;e.将剩下的区间0,1/16再一分为二，其中点为1/32,取区间1/32,1/16作为第四段; f.将剩下的区间0,1/32再一分为二，其中点为1/64,取区间1/64,1/32作为第三段;g.将剩下的区间0,1/64再一分为二，其中点为1/12

19、8,取区间1/128,1/64作为第二段;h.最后剩下的区间0,1/128作为第一段。然后将y轴的0,1区间均匀地分成八段，从第一段到第八段分别为0,1/8,(1/8,2/8,(2/8,3/8,(3/8,4/8,(4/8,5/8,(5/8,6/8,(6/8,7/8,(7/8,1。分别与x轴的八段一一对应。采用上述的方法就可以作出由八段直线构成的一条折线，该折线和A压缩律近似，图2.4中的八段线段的斜率分别为:表1 各段落的斜率段落12345678斜率161684211/21/4从上表中可以看出，除一、二段外，其他各段折线的斜率都不相同。图2.4中只画出了第一象限的压缩特性，第三象限的压缩特性的

20、形状与第一象限的压缩特性的形状相同，且它们以原点为奇对称，所以负方向也有八段直线，总共有16个线段。但由于正向一、二两段和负向一、二两段的斜率相同，所以这四段实际上为一条直线，因此，正、负双向的折线总共由13条直线段构成，这就是13折线的由来。从A律压缩特性中可以看出，取A=87.6主要基于下述两个原因: 1 使压缩特性曲线在原点附近的斜率为16;2 当用13折线逼近时，的八段量化分界点近似为1/2n(n=0,1,2,7)。从表1可以看出，当要求满足x=1/2n时，相应有y=1-n/8代入式中，有因此有将上式代入式(2.7),就可以得到对应A=94.4时的压缩特性(2.8)此压缩特性如果用13

21、折线逼近，除了第一段落起始点外，其余各段落的分界点的x、y都应满足式(2.8)。在13折线中，第一段落起始点要求的x、y都应该为零，而若按照式(2.8)计算时，当x=0时，y-;而当y=0，x=1/28。因此，需要对式(2.8)的压缩特性曲线作适当的修正，我们可以在原点和点(1/27,1/8)之间用一段直线代替原来的曲线，这段直线的斜率是1/8÷1/27=16。为了找到一个能够表示修正后的整个压缩特性曲线的方程，将式(2.8)变成(2.9)从上式中可以看出，它满足x=0时，y=0;x=1时，y=1。虽然式(2.9)在其他点上会有误差，但x在区间(1/128,1内，1+255x都能和原

22、来的256x比较接近。所以，在绝大部分范围内的压缩特性仍和A律压缩特性非常接近，只有在x0的小信号部分和A律压缩特性有些差别。若在式(2.9)中，令=255，则式(2.9)可写成(2.10) 式(2.10)的压缩特性与律压缩特性完全一致。（2）按照量化的维数分，量化分为标量量化和矢量量化。标量量化是一维的量化，一个幅度对应一个量化结果。而矢量量化是二维甚至多维的量化，两个或两个以上的幅度决定一个量化结果。以二维情况为例，两个幅度决定了平面上的一点。而这个平面事先按照概率已经划分为N个小区域，每个区域对应着一个输出结果（码数，codebook）。由输入确定的那一点落在了哪个区域内，矢量量化器就会

23、输出那个区域对应的码字（codeword）。矢量量化的好处是引入了多个决定输出的因素，并且使用了概率的方法，一般会比标量量化效率更高。2.2.4 编码（Coding）量化后的抽样信号在一定的取值范围内仅有有限个可取的样值，且信号正、负幅度分布的对称性使正、负样值的个数相等，正、负向的量化级对称分布。若将有限个 量化样值的绝对值从小到大依次排列，并对应地依次赋予一个十进制数字代码（例如，赋予样值0的十进制数字代码为0），在码前以“+”、“”号为前缀，来 区分样值的正、负，则量化后的抽样信号就转化为按抽样时序排列的一串十进制数字码流，即十进制数字信号。简单高效的数据系统是二进制码系统，因此，应将十

24、 进制数字代码变换成二进制编码。根据十进制数字代码的总个数，可以确定所需二进制编码的位数，即字长。这种把量化的抽样信号变换成给定字长的二进制码流的 过程称为编码。2.2基带传输PST码PST: Pair Selected Ternary成对选择三进码。在实际的基带传输系统中，并不是所有代码的电波形都能在信道中传输。含有直流分量和较丰富低频分量的单极性基带波形就不适宜在低频传输特性差的信道中传输。PST码是成对选择三进码。其编码过程是：先将二进制代码两两分组，然后再把每一码组编码成两个三进制数字(+ 、- 、0)。因为两位三进制数字共有9种状态，故可灵活地选择其中的4种状态。为防止PST码的直流

25、漂移，当在一个码组中仅发送单个脉冲时，两个模式应交替变换。例如：表1 PST码二进制代码+模式-模式0 0- +- +0 10 +0 -1 0+ 0- 01 1+ -+ -PST码能提供足够的定时分量，且无直流成分，编码过程也较简单。但这种码在识别时需要提供“分组”信息，即需要建立帧同步。 每位二进制信码都被变换成1位三电平取值（+1、0、-1）的码，因而有时把这类码称为lB1T码。2.3循环码2.3.1循环码介绍循环码是线性分组码的一种，所以它具有线性分组码的一般特性，此外还具有循环性。循环码的编码和解码设备都不太复杂，且检(纠)错能力强。它不但可以检测随机的错误，还可以检错突发的错误。（n

26、,k）循环码可以检测长为n-k或更短的任何突发错误，包括首尾相接突发错误。循环码是一种无权码，循环码编排的特点是相邻两个数码之间符合卡诺图中的邻接条件，即相邻两个数码之间只有一位码元不同，码元就是组成数码的单元。符合这个特点的有多种方案，但循环码只能是表中的那种。循环码的优点是没有瞬时错误，因为在数码变换过程中，在速度上会有快有慢，中间经过其它一些数码形式，称它们为瞬时错误。这在某些数字系统中是不允许的，为此希望相邻两个数码之间仅有一位码元不同，即满足邻接条件，这样就不会产生瞬时错误。循环码就是这样一种编码，它可以在卡诺图中依次循环得到。循环码最大的特点就是码字的循环特性，所谓循环特性是指：循

27、环码中任一许用码组经过循环移位后，所得到的码组仍然是许用码组。若（  ）为一循环码组，则（）、（）还是许用码组。也就是说，不论是左移还是右移，也不论移多少位，仍然是许用的循环码组。（1） 循环码的多项式表示设码长为n的循环码表示为（），其中为二进制数，通常把码组中各码元当做二进制的系数，即把上式中长为n的各个分量看做多项式： （2.11)的各项系数，则码字与码多项式一一对应，这种多项式中，x仅表示码元位置的标记,因此我们并不关心x的取值，这种多项式称为码多项式。（2）(n,k)循环码的生成多项式 (n,k)循环码的生成多项式写为g(x),它是(n,k)循环码码集中唯一的，幂次为n-k

28、的码多项式，则是一个幂次为n的码多项式。按模()运算，此时： （2.12）即 （2.13）且因 g(x)也是n阶幂，故Q(x)=1。由于它是循环码，故按模()运算后的“余式”也是循环码的一个码字，它必能被g(x)整除，即： （2.14）由以上两式可以得到： （2.15）和 （2.16）从上式中可以看出，生成多项式g(x)应该是的一个因式，即循环码多项式应该是的一个n-k次因式。（3） 循环码的生成矩阵和一致校验矩阵对所有的i=0,1,2,k-1，用生成多项式g(x)除，有： （2.17）式中是余式，表示为： （2.18）因此，是g(x)的倍式，即是码多项式，由此得到系统形式的生成矩阵为： (2

29、.19)它是一个kn阶的矩阵。同样，由G=0可以得到系统形式的一致校验矩阵为： （2.20） 如已知（7，4）循环码的生成多项式和校验多项式分别为：，。写得其生成矩阵和校验矩阵分别为： （2.21） （2.22） 2.3.2 循环码编码原理有信息码构成信息多项式，其中最高幂次为k-1；用乘以信息多项式m(x)，得到的，最高幂次为n-1，该过程相当于把信息码（，）移位到了码字德前k个信息位，其后是r个全为零的监督位；用g(x)除得到余式r(x),其次数必小于g(x)的次数，即小于（n-k），将此r(x)加于信息位后做监督位，即将r(x)于相加，得到的多项式必为一码多项式。（7，4）循环码编码的程

30、序框图如图2.5所示：初始化确定余式r(x): 确定c(x): 存储c(x)图2.7 编码程序框图2.3.3 循环码的纠错原理 纠错码的译码是该编码能否得到实际应用的关键所在。译码器往往比编码较难实现，对于纠错能力强的纠错码更复杂。根据不同的纠错或检错目的，循环码译码器可分为用于纠错目的和用于检错目的的循环码译码器。 通常，将接收到的循环码组进行除法运算，如果除尽，则说明正确传输；如果未除尽，则在寄存器中的内容就是错误图样，根据错误图样可以确定一种逻辑，来确定差错的位置，从而达到纠错的目的。用于纠错目的的循环码的译码算法比较复杂，感兴趣的话可以参考一些参考书。而用于检错目的循环码，一般使用AR

31、Q通信方式。检测过程也是将接受到的码组进行除法运算，如果除尽，则说明传输无误；如果未除尽，则表明传输出现差错，要求发送端重发。用于这种目的的循环码经常被成为循环冗余校验码，即CRC校验码。CRC校验码由于编码电路、检错电路简单且易于实现，因此得到广泛的应用。在通过MODEM传输文件的协议如ZMODEM、XMODEM协议中均用到了CRC校验技术。在磁盘、光盘介质存储技术中也使用该方法。 当码字c通过噪声信道传送时，会受到干扰而产生错误。如果信道产生的错误图样是e，译码器收到的n重接受矢量是y,则表示为： （2.23）上式也可以写成多项式形式： （2.24）译码器的任务就是从y(x)中得到，然后求

32、的估值码字 （2.25）并从中得到信息组。循环码译码可按以下三个步骤进行： （1）有接收到的y(x)计算伴随式s(x); （2）根据伴随式s(x)找出对应的估值错误图样； （3）计算，得到估计码字。若,则译码正确，否则，若，则译码错误。由于g(x) 的次数为n - k 次，g(x) 除E(x) 后得余式（即伴随式）的最高次数为n-k-1次，故S(x) 共有2n-k 个可能的表达式，每一个表达式对应一个错误格式。可以知道(7,4）循环码的S(x) 共有2(7-4) = 8个可能的表达式，可根据错误图样表来纠正(7,4）循环码中的一位错误，其伴随式如表2所示。表2 BCH（7，4）循环码错误图样表

33、错误图样错误图样码字伴随式S(x)伴随式E6(x)=x61000000x2100E5(x)=x50100000x2+x110E4(x)=x40010000x2+x+1111E3(x)=x30001000x+1011E2(x)=x20000100x2+1101E1(x)=x10000010x010E0(x)=x000000011001E(x)=000000000000上式指出了系统循环码的译码方法：将收到的码字R(x) 用g(x) 去除，如果除尽则无错；否则有错。如果有错，可由余式S(x) 一一找出对应图样，然后将错误图样E(x) 与R(x) 模2 和，即为所求码字C(x) ，从而实现纠错目的。

34、根据前面的讨论，可得（7，4）循环码译码的程序框图如图2.6所示：初始化图2.8译码程序框图否存储c(x)由S(x)确定错误图样E(x)S(x)=0，无误码误码由R(x)确定S(x)：纠错2.4二进制频移键控（2PSK）调制与解调原理2.4.1 数字调制技术为了使数字信号在带通信道中传输，必须用数字基带信号对载波进行调制以使信号与信道的特性相匹配。这种用数字信号控制载波，把数字基带信号变换为数字带通信号（已调信号）的过程成为数字调制。在接收端通过解调器把带通信号还原成数字基带信号的过程成为数字解调。通常把包括调制和解调过程的数字传输系统叫做数字带通传输系统。一般来说，数字调制与模拟调制的基本原

35、理相同，但是数字调制有离散取值的特点。因此数字调制技术有两种方法：利用模拟调制的方法去实现数字式调制；通过开关键控载波，通常称为键控法。基本键控方式有振幅键控、频移键控、相移键控。振幅键控 频移键控 相移键控图2.9正弦载波的三种键控波形2.4.2 二进制相移键控（2PSK）基本原理数字调制可分为二进制调制和多进制调制。相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中，通常用初始相位0和分别表示二进制“1”和“0”。因此，2PSK信号的时域表达式为：式中，j jn表示第n个符号的绝对相位：因此，上式可以改写为由于表示信号的两种码元的波形相同，极性相反，故2PSK信

36、号可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘：式中这里，g(t)是脉宽为Ts的单个矩形脉冲，而an的统计特性为即发送二进制符号“0”时，e2psk(t)取0相位；发送二进制符号“1”时，e2psk(t)取相位。这种以载波的不同相位直接去表示相应二进制数字信号的调制方式，称为二进制绝对相移方式。(a) 模拟调制的方法 （b）键控法图2.10 2PSK信号调制器原理方框图2PSK信号的解调通常采用相干解调法，解调器原理框图如下图2.9所示。假设数字信息为“10011”，相干解调各点波形如下图2.9所示。图2.11 2PSK信号的解调原理框图图2.12 2PSK信号相干解调时各点时间

37、波形2.4.3 MATLAB实现2PSK调制与解调1、2PSK调制流程图 2、2PSK解调流程图带通滤波器设计与载波相乘低通滤波器设计抽样判决2psk信号相干解调输出载波频率初始化数字基带信号输入符号1=>相位jn:符号0=>相位jn: 0控制载波产生2psk信号e2psk(t)=cos(2*fc*t+ jn)2psk信号输出图2.13 2psk调制 图2.14 2psk解调2.5衰落信道瑞利衰落信道（Rayleigh fading channel）是一种无线电信号传播环境的统计模型。这种模型假设信号通过无线信道之后，其信号幅度是随机的，表现为“衰落”特性，并且多径衰落的信号包络服

38、从瑞利分布。由此，这种多径衰落也称为瑞利衰落。 这一信道模型能够描述由电离层和对流层反射的短波信道，以及建筑物密集的城市环境。瑞利衰落只适用于从发射机到接收机不存在直射信号的情况，否则应使用莱斯衰落信道作为信道模型。假设经反射（或散射）到达接收天线的信号为N个幅值和相位均随机的且统计独立的信号之和。信号振幅为r,相位为,则其包络概率密度函数为 P(r)= (r0) （2.26）相位概率密度函数为： P()=1/2 （） （2.27）3 MATLAB软件仿真及结果3.1 MATLAB简介MATLAB是一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以

39、在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序（M文件）后再一起运行，可移植性好、可拓展性极强。MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。MATLAB的这些函数所能解决的问题大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等。3.2 发送端仿真结果与分析3.2.1 仿真结果图图3. 1 原始信号与PCM编码输出图3.2 PS

40、T编码图3.3二进制转三进制后循环码信号输出图3.4 PSK调制信号图3.5 滤波载波函数信号图3.6 衰落信道信号图3.7 PSK解调后图像 图3.8 循环码解码后二进制转三进制输出图像 图3.9 PST解码后的信号图 图3.10 PCM解调后的信号3.2.2 实验结果分析实验最后的结果与预想的还是差不多的，调制，编码和解调，解码的过程都是正确的。有几点可能对本设计存在较大影响：1.是取样的点数我取了50个不是很多会带来误差。2.再就是信道噪声会造成误差。3.PSK解调的过程中通过低通滤波器和抽样判决的过程也都是有误差的。4.最后用A律曲线逼近过程逼近时以一小段一小段的曲线去近似逼近。4.小

41、结在本次课程设计中我学习了整个通信系统的工作原理，了解了PCM编解码、循环码编解码、PST编解码、PSK调制解调的原理和MATLAB编程，也学习很多MATLAB的函数。虽然之前有接触MATLAB这个软件，但是也只是会调用一些简单的函数来做计算，或用一些简单的画图命令来作图，一直觉得这个仿真软件的功能很强大，但没有怎么深入地去运用。这次的课设让我知道，我之前对它的了解简直只是冰山一角。而要用它来实现自己的专业目的，更是增加了难度。将在课堂上学到的理论知识用它来仿真实现，这就要求我们对理论知识有全面深入透彻的理解，并且很熟悉仿真软件，但我在这两个方面都有欠缺，所以实际操作起来也没有刚开始想象的那样

42、容易。这次的课程设计让我收获颇丰，同时也让我发现了自身的不足。在此间发现的不足，我将努力改善，通过学习、实践等方式不断提高，克服那些不应成为学习、获得知识的障碍。这里要感谢帮助我解决编程过程中遇到问题的同学和老师。参考文献1 樊昌信 曹丽娜，通信原理第六版，国防工业出版社，20072 周开利，邓春晖主编 MATLAB基础及其应用教程，北京大学出版社，20073 董振海,精通MATLAB 7 编程与数据库应用，电子工业出版社,20074 陈怀琛，MATLAB在电子信息课程中的应用（第二版），电子工业出版社，20065 张德丰主编，MATLAB通信工程仿真，机械工业出版社，2010附录：%-主函数

43、-%clcclose allclear all%-生成时间函数-%t=1:0.1:5; s=sin(2*3.14*1000*t);%-pcm调制-%pcmcode,sMax=pcm(s);figure(1)subplot(2,1,1);plot(s); title('原始正弦信号');subplot(2,1,2);stairs(pcmcode);zoom xon;axis(0 750 -1 2);%-pst编码-%pstcode=pst_c(pcmcode);figure(2);stairs(pstcode);zoom xon;axis(0 750 -1 2);%-三进制转为二

44、进制-%code=ZH_C(pstcode);figure(3);stairs(code);zoom xon;axis(0 750 -1 2);%-循环码编码-%xh=XH_C(code);figure(4);stairs(xh);zoom xon;axis(0 750 -1 2);%-psk调制-%psk,bit,mod=psk_encode(xh);figure(5); plot(psk); zoom xon;axis(0 750 -1 2);% figure(12);plot(psk);axis(0 1000 -1.2 1.2);title('PSK Modulation'

45、;);%-衰落信道-%chanel=rayleighchan(1/64000,100);%生产瑞丽衰落信道noisig0=filter(chanel,psk);%把处理过的信号送入信道noisig=abs(noisig0);figure(6); plot(noisig);zoom xon;title('通过衰落信道后的波形');% figure;plot(noisig);axis(0 1000 -1.2 1.2)%-psk解调-%pskdec=psk_decode(psk,mod);figure(10);stairs(pskdec);zoom xon;axis(0 750 -1

46、 2);%-循环码解码-%xm=XH_D(pskdec);figure(11);stairs(xm);zoom xon;axis(0 750 -1 2);%-二进制转为三进制-%dcode=ZH_D(xm);figure(12);stairs(dcode);zoom xon;axis(0 750 -1 2);%-pst解码-%pstdec=repst(dcode);figure(13);stairs(pstdec);zoom xon;axis(0 750 -1 2);%-pcm解调-%pcmdec=pcm_de(pstdec,sMax);plot(pcmdec); title('还原的

47、正弦信号');%-PCM调制-%function code_lin,sMax=pcm(signal)sSig=sign(signal); %判断signal的正负sMax=max(abs(signal); %求signal的最大值Q=abs(signal/sMax)*2048; %归一化量化code=zeros(length(signal),8); %存储矩阵（全零）%段落码判断%for m=1:length(signal) if Q(m)>128 && Q(m)<=2048 %在第五段与第八段之间，段位码第一位都为“1” code(m,2)=1; end

48、if (Q(m)>32 && Q(m)<128) | (Q(m)>512 && Q(m)<=2048) code(m,3)=1; %在第三四七八段内，段位码第二位为“1” end if (Q(m)>16&&Q(m)<32)|(Q(m)>64&&Q(m)<128)|(Q(m)>256&&Q(m)<512)|(Q(m)>1024&&Q(m)<=2048) code(m,4)=1; %在二四六八段内，段位码第三位为“1” end %符号位的判断