通信系统综合课程设计

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

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

3、实验原理框图 22.2 PCM原理介绍32.2.1 扌由样（Sampling） 32.2.2 量化（quantizing ） 32.2.3 A压缩律：42.2.4 编码（Coding ） 92.2基带传输PST码92.3循环码10循环码介绍10循环码编码原理12循环码的纠错原理 132.4二进制频移键控（2PSK调制与解调原理 152.4.1 数字调制技术15二进制相移键控（2PSK）基本原理 16实现2PSK调制与解调 182.5衰落信道183 MATLAB软件仿真及结果 193.1 MATLAB 简介193.2 发送端仿真结果与分析 19仿真结果图193.2.2 实验结果分析 244.小结

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

5、 MATLAB件编程，完成了整个通信系统编码仿真分析。关键词：MATLAB件编程、数字通信系统、PCM PST码、循环码、PSKAbstractCommunication system is the role of the information from the source sent to one or more of the desti nati on. A digital com muni cati on system is to use digital sig nal to tran sfer in formatio n com muni cati on system, basica

6、lly have source codi ng and decoding, channel coding and decoding, digital modulation and demodulation. The specific process is the source of the an alog sig nal after a/d convert into nu mber son signal and then, again carries on the baseband coding, channel coding, carrier modulation and then sent

7、 to the channel, then made on the receiving end corresponding demodulation, decoding, d/a conversion to get the final receiving sig nal.This design with PCM sampling and quantization and coding principle to an alog sig nal into the digital sig nal, ham ming code kn itt ing decod ing theory to realiz

8、e channel knitting decoding, PST knitting decoding principle used to implement the source knitting decoding, PSK modulation demodulation principle used to implement the digital modulati on and demodulati on, Fadi ng cha nnel and the prin ciple for no ise approximation, using MATLAB software programm

9、ing, complete the whole com muni cati on system cod ing simulatio n an alysis.Keywords: MATLAB software program ming, digital com muni catio n system, PCM, PST, Cyclic code, PSK1设计任务1.1设计目的通过课程设计，使学生加强对电子电路的理解，学会对电路分析计算以及设计。进一 步提高分析解决实际问题的能力，通过完成综合设计型和创新性实验及训练，创造一个动 脑动手、独立开展电路实验的机会，锻炼分析、解决电子电路问题的实际本

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

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

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

13、以下原理框图图2.2实验原理框图2.2 PCM原理介绍模拟信号数字化必须经过三个过程，即抽样、量化和编码，以实现话音数字化的脉冲编码调制(PCM，Pulse Coding Modulation)技术。2.2.1 抽样(Sampli ng)离散时间信号通常是有连续时间信号经周期采样得到的。完成采样功能的器件称为采样器，下图所示为采样器的示意图。图中Xa(t)表示模拟信号，Xa(nt)表示采样信号，T为采样周期，n=0，1，2，。一般可以把采样器视为一个每隔 T秒闭合一次的电子开关So 在理想情况下，开关闭合时间 t满足t v<T。实际采样过程可视为脉冲调幅过程，Xa(t)为调制信号，被调脉

14、冲载波p(t)是周期为T、脉宽为t的周期脉冲串。当T 0时的理想 采样情况是实际采样的一种科学的、本质的抽象，同时可使数学推导得到简化。下面主要 讨论理想采样。2.2.2 量化(quantizing )抽样信号虽然是时间轴上离散的信号，但仍然是模拟信号，其样值在一定的取值范围 内，可有无限多个值。显然，对无限个样值给出数字码组来对应是不可能的。为了实现以数字码表示样值，必须采用“四舍五入”的方法把样值分级“取整”，分的级数越多,即量化级差或间隔越小，量化噪声也越小。223 A压缩律：所谓A压缩律也就是压缩器具有如下特性的压缩律:1 十 In / y - <1十血川穴.1 + ln A其中

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

16、中k为比例常数。当量化级数很大时，可以将它看成连续曲线，因而式(2.2)成为线性微分方程解此微分方程其中C为常数。为了满足归一化要求，当x=1时,k H-c 0y=1,代入式(2.3)可得故所得结果为In工二炒一氏如果压缩特性满足上式，就可获得理想的压缩效果,(2.2)(2.3)(2.4)其量化信噪比和信号幅度无关。满图2.5理想压缩特性曲线足上式的曲线如下图所示,定的修改。A律压缩特性就是对式(2.4)修改后的函数。在上图中，通过原点作理想压缩特性曲线 的切线0C,将oc、cd作为实际的压缩特性。修改以后，必须用两个不同的方程来描述这 段曲线，以切点c为分界点，线段0C的方程:设切点c的坐标

17、为(x1，y1)斜率为ck则由式(2.4)可得(2.5)所以线段oc的方程为1所以当x=x1时，y1=1/k时，有1 , L=1 + 一n x, k k因此有芯=严宀所以，切点坐标为(exp-(k-1),1/k)，令可二丄二严1丄则上三1+ In将它代入式(2.5),就可得到以切点c为边界的段的方程为y =, 0 c x < (2.6)(2.7)l+lii JA因cd段的方程，满足式(3.4),所以由该式可得1.1 + In 貝 xV = 1 +In 芹匚1+1hj4 Uln A由以上分析可见，经过修改以后的理想压缩特性与图5中所示的曲线近似，而式(2.6)式(2.7)和式(2.4)完全

18、一样13折线：实际中，A压缩律通常采用13折线来近似，13折线法如图2.4所示，图中 先把丫轴的0,1区间分为8个不均匀段。7/86/84/fl2/8/81/41/2*二 u图2.6折线示意图1/64 i/32其具体分法如下：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再一分

19、为二，其中点为1/32,取区间1/32,1/16作为第四段;f. 将剩下的区间0,1/32再一分为二，其中点为1/64,取区间1/64,1/32作为第三段；g. 将剩下的区间0,1/64再一分为二，其中点为1/128取区间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

20、.4中的八段线段的斜率分别为：表1各段落的斜率段落12345678斜率161684211/21/4从上表中可以看出，除一、二段外，其他各段折线的斜率都不相同。图2.4中只画出了第一象限的压缩特性，第三象限的压缩特性的形状与第一象限的压缩特性的形状相同，且它们以原点为奇对称，所以负方向也有八段直线，总共有16个线段。但由于正向一、二两段和负向一、二两段的斜率相同，所以这四段实际上为一条直线，因此，正、负双向 的折线总共由13条直线段构成，这就是13折线的由来。从A律压缩特性中可以看出，取 A=87.6主要基于下述两个原因：1使压缩特性曲线在原点附近的斜率为16;2当用13折线逼近时，的八段量化分

21、界点近似为1/2八n(n=0,1,2,7)。从表1可以看出，当要求满足x=1/2An时，相应有y=1-n/8代入式中，有1 _ 1 _ 1 _ 1f =冏亠=严=用f因此有2 =将上式代入式(2.7),就可以得到对应A=94.4时的压缩特性(2.8)1 + ln (j1)In(256a)l + lii Z山宓) In 256此压缩特性如果用13折线逼近，除了第一段落起始点外，其余各段落的分界点的 x、 y都应满足式(2.8)。在13折线中，第一段落起始点要求的 x、y都应该为零，而若按照式 (2.8)计算时，当x=0时，Ls;而当y=0, 乂=1/2八8。因此，需要对式(2.8)的压缩特性曲线

22、 作适当的修正，我们可以在原点和点(1/2A7,1/8)之间用一段直线代替原来的曲线，这段直线 的斜率是1/8- 1/2八7=16。为了找到一个能够表示修正后的整个压缩特性曲线的方程，将式(2.8)变成_ ln(1 + 255z)-(2.9)从上式中可以看出，它满足x=0时，y=0;x=1时，y=1。虽然式(2.9)在其他点上会有误 差，但x在区间(1/128,1 内，1+255X都能和原来的256x比较接近。所以，在绝大部分范 围内的压缩特性仍和A律压缩特性非常接近，只有在x-0的小信号部分和A律压缩特性 有些差别若在式(2.9)中，令 卩=255，则式(2.9)可写成（2.10）ha 1

23、+ 占）血1+山式(2.10)的压缩特性与 卩律压缩特性完全一致。(2)按照量化的维数分，量化分为标量量化和矢量量化。标量量化是一维的量化， 一个幅度对应一个量化结果。而矢量量化是二维甚至多维的量化，两个或两个以上的幅度 决定一个量化结果。以二维情况为例，两个幅度决定了平面上的一点。而这个平面事先按照概率已经划分 为N个小区域，每个区域对应着一个输出结果(码数，codebook)。由输入确定的那一点落在了哪个区域内，矢量量化器就会输出那个区域对应的码字(codeword)。矢量量化的好处是引入了多个决定输出的因素，并且使用了概率的方法，一般会比标量量化效率更高。编码（Coding）量化后的抽样

24、信号在一定的取值范围内仅有有限个可取的样值，且信号正、负幅度分布的对称性使正、负样值的个数相等，正、负向的量化级对称分布。若将有限个量化样值的绝对值从小到大依次排列，并对应地依次赋予一个十进制数字代码(例如，赋予样值 0的十进制数字代码为0),在码前以“ +”、“一”号为前缀，来 区分样值的正、负，则 量化后的抽样信号就转化为按抽样时序排列的一串十进制数字码流，即十进制数字信号。简单高效的数据系统是二进制码系统，因此，应将十进制数字代码变换成二进制编码。根据十进制数字代码的总个数，可以确定所需二进制编码的位数，即字长。这种把量化的 抽样信号变换成给定字长的二进制码流的 过程称为编码。2.2基带

25、传输 PST码PST:Pair Selected Ternary成 对选择三进码。在实际的基带传输系统中，并不是所有代码的电波形都能在信道中传输。含有直流分 量和较丰富低频分量的单极性基带波形就不适宜在低频传输特性差的信道中传输。PST码是成对选择三进码。其编码过程是：先将二进制代码两两分组，然后再把每一码组编码成 两个三进制数字(+、-、0)。因为两位三进制数字共有9种状态，故可灵活地选择其中的 4种状态。为防止PST码的直流漂移，当在一个码组中仅发送单个脉冲时，两个模式应交 替变换。例如：表1 PST码二进制代码+模式-模式0 0-+-+0 10 +0 -1 0+ 0-01 1+ -+ -

26、PST码能提供足够的定时分量，且无直流成分，编码过程也较简单。但这种码在识别时需要提供 分组”信息，即需要建立帧同步。 每位二进制信码都被变换成1位三电平取值 （+1、0、-1）的码，因而有时把这类码称为IB/ 1T码。2.3循环码231循环码介绍循环码是线性分组码的一种，所以它具有线性分组码的一般特性，此外还具有循环性。 循环码的编码和解码设备都不太复杂，且检 （纠）错能力强。它不但可以检测随机的错误， 还可以检错突发的错误。（n,k）循环码可以检测长为n-k或更短的任何突发错误，包括首 尾相接突发错误。循环码是一种无权码，循环码编排的特点是相邻两个数码之间符合卡诺图中的邻接条 件，即相邻两

27、个数码之间只有一位码元不同，码元就是组成数码的单元。符合这个特点的 有多种方案，但循环码只能是表中的那种。循环码的优点是没有瞬时错误，因为在数码变 换过程中，在速度上会有快有慢，中间经过其它一些数码形式，称它们为瞬时错误。这在 某些数字系统中是不允许的，为此希望相邻两个数码之间仅有一位码元不同，即满足邻接 条件，这样就不会产生瞬时错误。循环码就是这样一种编码，它可以在卡诺图中依次循环 得到。循环码最大的特点就是码字的循环特性，所谓循环特性是指：循环码中任一许用码组 经过循环移位后，所得到的码组仍然是许用码组。若（an 1 an 2aj a。）为一循环码组，则（an 2 an 3a。a. 1）、

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

29、码多项式,则xkg(x)是一个幕次为n的码多项式。按模(xn 1)运算，此时:Q ( x) xkg(x)R (x)nR(x)(2.12)(2.13)且因xk g(x)也是n阶幕，故Q(x)=1。由于它是循环码，故xkg(x)按模(xn 1)运算后的 余 式”也是循环码的一个码字，它必能被g(x)整除，即：3 F(x)(214)由以上两式可以得到：g(x)xkg(x)Q(x)(xn1) R(x) (xn 1) f(x)g(x)( 2.15)和xn1 xk f (x) g(x) h(x)g(x)(2.16)从上式中可以看出，生成多项式g(x)应该是xn 1的一个因式，即循环码多项式应该是xn 1的

30、一个n-k次因式。(3)循环码的生成矩阵和一致校验矩阵对所有的i=0,1,2, -1,用生成多项式g(x)除xn ki，有：xn k ：a： (x)g(x) bi(x)(2.17)式中bi(x)是余式，表示为：bi (x)b,n k 1Xn k 1b,1Xbi,0(2.18)因此，xn k i b(x)是g(x)的倍式，即xn k 1bi(x)是码多项式，由此得到系统形式的生成矩阵为：1oobk 1,n k 1bk 1,1bk 1,oGo1obk 2,n k 1bk 2,1bkx,ooooobo, n k 1bo,1bo,o(2.19)它是一个 k n 阶的矩阵。同样，由 G HT=0 可以得

31、到系统形式的一致校验矩阵为：bk 1, n k 1bk 2,n k 1bo,n k 11ooH( 2.2o)bk 1,1bk 2,1bo,1o1obk 1,obk 2,obo,ooo1g(x) x3 x 1 ，42h(x) xxx 1 。写得其生成矩阵和校验矩阵分别为1o11oooo1o11ooGoo1o11oooo1o112.21)如已知( 7， 4)循环码的生成多项式和校验多项式分别为：111o1ooHo111o1ooo111o12.22)2.3.2 循环码编码原理有信息码构成信息多项式 m(x) mk 1xk 1mo，其中最高幕次为k-1 ；用xn k乘以信息多项式m(x)，得到的xn

32、km(x)，最高幕次为n-1 ，该过程相当于把信息码mk 1 ，mk 2,m!, mo)移位到了码字德前k个信息位，其后是r个全为零的监督位；用g(x)除xn km(x)得到余式r(x),其次数必小于g(x)的次数，即小于(n-k)，将此r(x)加于信息位后做监督位，即将r(x)于xn km(x)相加，得到的多项式必为一码多项式。(7, 4)循环码编码的程序框图如图 2.5 所示：图2.7编码程序框图循环码的纠错原理纠错码的译码是该编码能否得到实际应用的关键所在。译码器往往比编码较难 实现，对于纠错能力强的纠错码更复杂。根据不同的纠错或检错目的，循环码译码器 可分为用于纠错目的和用于检错目的的

33、循环码译码器。通常，将接收到的循环码组进行除法运算，如果除尽，则说明正确传输；如果未 除尽，则在寄存器中的内容就是错误图样，根据错误图样可以确定一种逻辑，来确定 差错的位置，从而达到纠错的目的。用于纠错目的的循环码的译码算法比较复杂，感 兴趣的话可以参考一些参考书。而用于检错目的循环码，一般使用ARQ通信方式。检测过程也是将接受到的码组进行除法运算，如果除尽，则说明传输无误；如果未除尽，则表明传输出现差错， 要求发送端重发。用于这种目的的循环码经常被成为循环冗余 校验码，即CRC校验码。CRC校验码由于编码电路、检错电路简单且易于实现，因此 得到广泛的应用。在通过MODEMS输文件的协议如 Z

34、MODEMXMODE协议中均用到了CRC校验技术。在磁盘、光盘介质存储技术中也使用该方法。当码字c通过噪声信道传送时，会受到干扰而产生错误。如果信道产生的错误图样是e,译码器收到的n重接受矢量是y,则表示为:yc e(2.23)上式也可以写成多项式形式：y(x)c(x)e(x)(2.24)译码器的任务就是从y(x)中得到Ae(x)，然后求的估值码字Ac(x)y(x)Ae(x)(2.25)并从中得到信息组m(x)循环码译码可按以下三个步骤进行：(1) 有接收到的y(x)计算伴随式s(x);(2) 根据伴随式s(x)找出对应的估值错误图样e(x)；AAA(3) 计算c(x) y(x) e(x)，得

35、到估计码字c(x)。若c(x) c(x)，则译码正确，否则，A若c(x)c(x)，则译码错误。由于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)循环码错误图样表错误图样错误图样码字伴随式S(x)伴随式E6(X)=X 61000000x2100E5(X)=X 501000002 X +x110E4(X)=X

36、 400100002x2+x+1111E3(X)=X 30001000x+1011E2(X)=X 20000100x2+1101E1 (x)=x 10000010X010Eo(x)=x 000000011001E(x)=000000000000上式指出了系统循环码的译码方法：将收到的码字R(x)用g(x)去除，如果除尽则无错；否则有错。如果有错，可由余式 S(x) 一一找出对应图样，然后将错误图样E(x)与R(x)模2和，即为所求码字C(x)，从而实现纠错目的。根据前面的讨论，可得(7, 4)循环码译码的程序框图如图2.6所示：图2.8译码程序框图2.4二进制频移键控(2PSK调制与解调原理数

37、字调制技术为了使数字信号在带通信道中传输，必须用数字基带信号对载波进行调制以使信号与 信道的特性相匹配。这种用数字信号控制载波，把数字基带信号变换为数字带通信号(已 调信号)的过程成为数字调制。在接收端通过解调器把带通信号还原成数字基带信号的过 程成为数字解调。通常把包括调制和解调过程的数字传输系统叫做数字带通传输系统。一般来说，数字调制与模拟调制的基本原理相同，但是数字调制有离散取值的特点。 因此数字调制技术有两种方法： 利用模拟调制的方法去实现数字式调制； 通过开关键控载波，通常称为键控法。基本键控方式有振幅键控、频移键控、相移 键控振幅键控频移键控相移键控图2.9正弦载波的三种键控波形二

38、进制相移键控(2PSK基本原理2PSK数字调制可分为二进制调制和多进制调制相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在中，通常用初始相位0和n分别表示二进制“T和“ 0”。因此，2PSK信号的时域表达式为:e2PSK（t）Acos（ctn）式中，n表示第n个符号的绝对相位：0，发送,发送“因此，上式可以改写为0”1”A cos ct, A cos ct, 由于表示信号的两种码元的波形相同，极性相反,e2PSK（t）概率为概率为1故2PSK信号可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘:e2psK(t) st cos ct式中s(t)ang(t nTs)n

39、这里，g(t)是脉宽为Ts的单个矩形脉冲，而an的统计特性为1, 概率为P1,概率为1 P即发送二进制符号“ 0”时，e2psk(t)取0相位；发送二进制符号“ 1”时，Qpsk(t)取n相位。 这种以载波的不同相位直接去表示相应二进制数字信号的调制方式，称为二进制绝对相移 方式。双极性(a)模拟调制的方法开关电路(b)键控法图2.10 2PSK信号调制器原理方框图2PSK信号的解调通常采用相干解调法，解调器原理框图如下图2.9所示。假设数字信息为“ 1001T,相干解调各点波形如下图2.9所示。图2.11 2PSK信号的解调原理框图abcde图2.12 2PSK信号相干解调时各点时间波形24

40、3 MATLAB实现2PSK调制与解调1、2PSK调制流程图图2.13 2psk调制2、2PSK解调流程图图2.14 2psk解调2.5衰落信道瑞利衰落信道(Rayleigh fading channe)是一种无线电信号传播环境的统计模型。这种模型假设信号通过无线信道之后，其信号幅度是随机的，表现为衰落”特性，并且多径衰落的信号包络服从瑞利分布。由此，这种多径衰落也称为瑞利衰落。这一信道模型能够描述由电离层和对流层反射的短波信道，以及建筑物密集的城市环境。瑞利衰落只适用 于从发射机到接收机不存在直射信号的情况，否则应使用莱斯衰落信道作为信道模型。假设经反射(或散射)到达接收天线的信号为N个幅值

41、和相位均随机的且统计独立的信号之和。信号振幅为r,相位为，则其包络概率密度函数为(2.26)(2.27)rtP(r)= 2 e 2 (r 0)相位概率密度函数为：P( )=1/2( 02 )3 MATLAB软件仿真及结果3.1 MATLAB 简介MATLAB是一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输 出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编 写好一个较大的复杂的应用程序（M文件）后再一起运行，可移植性好、可拓展性极强。MATLAB是一个包含大量计算算法的集合。 其拥有600多个工程中要用到的数学运算 函数，可以方便的实现用户所需的

42、各种计算功能。MATLAB的这些函数所能解决的问题大 致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅 立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角 函数和其他初等数学运算、多维数组操作以及建模动态仿真等。3.2发送端仿真结果与分析3.2.1 仿真结果图图3. 1原始信号与PCk编码输出图3.2 PST编码图3.3二进制转三进制后循环码信号输出P图3.4 PSK调制信号带迫遊液后国埶2F2D01002003(10400500600 7DC別 0900n&叵1” 图3.6衰落信道信号图3.9 PST解码后的信号图解吗信号图3

43、.10 PCM解调后的信号322实验结果分析实验最后的结果与预想的还是差不多的， 调制，编码和解调，解码的过程都是正确的 有几点可能对本设计存在较大影响：1.是取样的点数我取了 50个不是很多会带来误差。 2.再就是信道噪声会造成误差。3.PSK解调的过程中通过低通滤波器和抽样判决的过程也 都是有误差的。4.最后用A律曲线逼近过程逼近时以一小段一小段的曲线去近似逼近。4.小结在本次课程设计中我学习了整个通信系统的工作原理，了解了PCM编解码、循环码编解码、PST编解码、PSK调制解调的原理和 MATLAB®程，也学习很多 MATLA的函数。虽然之前有接触MATLA这个软件，但是也只是

44、会调用一些简单的函数来做计算， 或用 一些简单的画图命令来作图，一直觉得这个仿真软件的功能很强大，但没有怎么深入地去 运用。这次的课设让我知道，我之前对它的了解简直只是冰山一角。而要用它来实现自己 的专业目的，更是增加了难度。将在课堂上学到的理论知识用它来仿真实现，这就要求我 们对理论知识有全面深入透彻的理解， 并且很熟悉仿真软件，但我在这两个方面都有欠缺， 所以实际操作起来也没有刚开始想象的那样容易。这次的课程设计让我收获颇丰，同时也让我发现了自身的不足。在此间发现的不足， 我将努力改善，通过学习、实践等方式不断提高，克服那些不应成为学习、获得知识的障 碍。这里要感谢帮助我解决编程过程中遇到

45、问题的同学和老师。参考文献1 樊昌信 曹丽娜，通信原理第六版，国防工业出版社，20072 周开利，邓春晖主编 MATLAB基础及其应用教程，北京大学出版社，20073 董振海，精通MATLAB 7编程与数据库应用，电子工业出版社,20074 陈怀琛，MATLAB在电子信息课程中的应用（第二版），电子工业出版社，2006 张德丰主编，MATLAB通信工程仿真，机械工业出版社，2010附录：% 主函数 %clc close all clear all% 生成时间函数 %t=1:0.1:5; s=sin(2*3.14*1000*t);%pcm 调制 %pcmcode,sMax=pcm(s);figu

46、re(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);% 三进制转为二进制 %code=ZH_C(pstcode);figure(3);stairs(code);zoom xon;axis(0 750 -1 2);% 循环码编码 %xh=XH_C(code)

47、;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'); % 衰落信道 %chan el=rayleighcha n(1/64000,100);%生产瑞丽衰落信道 noisig0=filter(chanel,psk);% 把处理过的信号送

48、入信道 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 2);% 循环码解码 %xm=XH_D(pskdec);figure(11);stairs(xm);zoom xon;axis(0 750 -1 2);% 二进制转为三进

49、制 %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(' 还原的正弦信号 '); %PCM 调制 %判断 signal 的正负%求 signal 的最大值%归一化量化%存储矩阵(全零)function code_lin,sM

50、ax=pcm(signal) sSig=sign(signal); sMax=max(abs(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;endif (Q(m)>32 && Q(m)<128) | (Q(m)>512 && Q(m)<=2048)cod

51、e(m,3)=1;%在第三四七八段内，段位码第二位为“ 1endif (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%符号位的判断if sSig(m)>0code(m,1)=1;elseif sSig(m)<0code(m,1)=0;endend%段内码判断程序N=zeros(len

52、gth(signal),1);for m=1:length(signal)N(m)=code(m,2)*4+code(m,3)*2+code(m,4)+1;% 找到 Y 位于第几段enda=0,16,32,64,128,256,512,1024;%量化间隔b=1,1,2,4,8,16,32,64;%除以 16，得到每段的最小量化间隔for m=1:length(signal)q=ceil(Q(m)-a(N(m)/b(N(m);if q=0code(m,(5:8)=0,0,0,0;else k=num2str(dec2bin(q-1,4);code(m,5)=str2num(k(1);code(m,6)=str2num(k(2);code(m,7)=str2num(k(3);code(m,8)=str2num(k(4);end%求出在段内的位置%如果输入为零则输出“ 0”%编码段内码为二进制end%将 N 行 8 列矩阵转换为 1行 8*N 列的矩阵 code_lin=zeros(size(code,1)*8,1);c=1;for i=1:size(code,1)for j=1:8code_lin(c)=code(i,j); c=c+1;endendcode_lin=code_lin' %code_lin End%PST 编码 % functiony=pst_c(