通信原理课程设计报告数字信号频带传输系统设计.doc

《通信原理课程设计报告数字信号频带传输系统设计.doc》由会员分享，可在线阅读，更多相关《通信原理课程设计报告数字信号频带传输系统设计.doc（38页珍藏版）》请在三一办公上搜索。

1、课 程 设 计 报 告课程名称 通信原理课程设计 课题名称 数字信号频带传输系统设计 专 业 通信工程 班 级 学 号 姓 名 指导教师 2011年12月 12 日湖南工程学院课 程 设 计 任 务 书课程名称 通信原理课程设计 课 题 数字信号频带传输系统设计 专业班级 学生姓名 学 号 指导老师 审 批 任务书下达日期 2011 年 12月 12 日任务完成日期 2011 年 12 月 26日一、设计目的与设计内容结合通信原理课程内容，培养实际动手能力，增强把书本知识转化为实践设计的能力。在强调基本原理的同时，更突出设计过程的锻炼。强化学生的实践创新能力和独立进行科研工作的能力。要求学生经

2、过课程设计这一教学环节学会应用MATLAB软件来实现信号传输中的各个基本环节。1设计内容：具体设计题目包括：模拟调制、数字调制、PCM编解码、循环码的编译码、m序列伪随机码、最佳接收等内容。学生还可以根据自己的兴趣，围绕通信原理课程知识自由选题。题目1 信息论基本计算要求：（1）编程实现信源平均信息量的计算(以高斯分布的信源为例)（2）编程实现离散信道容量的计算(以输入符号等概分布为例)（3）编程实现信源编码过程(以Huffman编码为例) 题目2 PCM编译码器设计及应用要求：（1）用simulink对系统建模设计一个PCM编译码器（2）输入模拟话音信号观察其输出波形（3）输入数字信号观察其

3、输出波形（4）对所设计的系统性能进行仿真分析（5）对其应用举例阐述题目3 数字信号基带传输系统设计要求：（1）编程实现常见基带信号的波形、码型转换：非归零/归零码、AMI码（2）并绘出每种波形、码型的功率谱分布，给出与在理论课上所学相符合的分析与理解（3）编程实现基带传输系统的误码率计算：二电平误码率计算（4）编程实现基带信号传输的扰码与解扰 题目4 数字信号频带传输系统设计（1）2ASK调制/相干解调器和非相干解调系统（包络检波法）设计要求：a编程实现或是用simulink对系统建模 b输入数字信号序列并进行接收判决c通过多次输入输出对两种系统性能进行分析比较d绘出信号的时域波形和频谱图（2

4、）2FSK调制/相干解调器和非相干解调系统（包络检波法）设计要求：a编程实现或是用simulink对系统建模 b输入数字信号序列并进行接收判决 c通过多次输入输出对两种系统性能进行分析比较d绘出信号的时域波形和频谱图（3）2PSK、2DPSK调制/相干解调系统设计要求：a编程实现或是用simulink对系统建模 b输入数字信号序列并进行接收判决 c通过多次输入输出对两种系统性能进行分析比较d绘出信号的时域波形和频谱图（4）16/64QAM调制/解调系统设计要求：a编程实现或是用simulink对系统建模 b输入数字信号序列并进行接收判决 c通过多次输入输出对两种系统性能进行分析比较d绘出信号的

5、时域波形和频谱图题目5 设计一个先验等概的2ASK最佳接收机 要求：（1）编程实现或是用simulink对系统建模 （2）输入数字信号序列并进行接收判决 （3）通过多次输入输出对所设计的系统性能进行分析题目6 设计一个DSB调制/解调（用平方环和科斯塔斯环实现载波提取）电路要求： （1）编程实现或是用simulink对系统建模（2）输入模拟话音信号观察其输出波形（3）对所设计的系统性能进行仿真分析（4）对其应用举例阐述题目7 数字通信频带传输系统综合设计必做题目：题目1，2，3；选作题目：4（1）、4（2）、4(3)、4(4)，5，6中的一个；综合题目7：每班同学中，做2、3、4、5题目的各选

6、2个同学，共计6人合作来完成题目7。2设计要求：1、在一周内学生须上机16小时以上，程序调试完后，须由指导老师在机器上检查运行结果，经教师认可后的源程序可通过打印机输出，并请教师在程序清单上签字。2、课程设计报告内容和格式：设计题目，设计的详细步骤，设计过程中的结果、图形等，设计总结。3、每组每人必须独立完成，成绩的考核按设计结果、答辩成绩及课程设计报告来综合评定。成绩分为优、良、中、及格、不及格五级分评定。4、指导教师：彭祯，谭小兰。二、进度安排第十六周 星期三 8：0012：00 星期五 8：0012：00 星期日 8：0012：00 星期日 14：0018：00第十七周 星期二 14：0

7、018：00 星期四 14：0018：00附：课程设计报告装订顺序：封面、任务书、目录、正文、评分、附件（A4大小的图纸及程序清单）。 正文的格式:一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。正文的内容:一、课题的主要功能；二、课题的功能模块的划分（要求画出模块图）；三、主要功能的实现（至少要有一个主要模块的流程图）；四、程序调试；五、总结；六、附件（所有程序的原代码，要求对程序写出必要的注释）。正文总字数要求在5000字以上（不含程序原代码）。 目录一、 主要功能.7二、 功能模块的划分.71、 信源平均信息量的计算(以高斯分布的信源为例).72、 离散信道容量

8、的计算(以输入符号等概分布为例).73、 信源编码过程(以Huffman编码为例).74、 PCM编译码器设计及应用.85、 数字信号基带传输系统设计.86、 数字信号频带传输系统设计.9三、 主要功能的实现1、信源平均信息量的计算(以高斯分布的信源为例).102、离散信道容量的计算(以输入符号等概分布为例).103、信源编码过程(以Huffman编码为例).114、PCM编译码器设计及应用.115、编程实现常见基带信号的波形、码型转换包括：单极性非归零码、单极性归零码、双极性非归零码、双极性归零码、AMI码.126、波形、码型的功率谱分布的绘图，包括：单极性非归零码、单极性归零码、双极性非归

9、零码、双极性归零码、AMI码.157、编程实现基带传输系统的误码率计算：二电平误码率计算.18四、功能模块的分析与流程图.211、信息论基本运算模块.212、PCM编译码的设计及运用模块.21五、程序及Simulink调试.23六、心得体会.23七、评分表.25八、附件（程序清单）.26一、 主要功能 MATLAB是由MathWorks公司开发的一种主要用于数值计算及可视化图形处理的工程语言，是当今最优秀的科技应有软件之一。它将集数值分析、矩阵运算、图形图像处理、信号处理和仿真等诸多强大的功能集成在较易使用的交互计算机环境之中，为科学研究、工程应用提供了一种功能强、效率高的编程工具。本次课程设

10、计学生通过自己设计及建立通信系统，掌握通信系统的构成原理、信号传输的概念，加深对通信系统及信号的认识，提高自身的实际应用、动手能力。最终使学生通过课程设计学会应用MATLAB软件来实现信号传输中的各个基本环节。二、 功能模块的划分1、信源平均信息量的计算(以高斯分布的信源为例) 根据题目要求，用高斯过程X(t)d的一维概率密度函数服从正态分布的表达式f=1/(u*sqrt(2*pi)*exp(-(x-o)2/2*u2来完成信源平均信息量的计算。平均信息量：平均每个符号所能提供的信息量，也叫平均自信息量。H(X)= ；高斯分布函数：；2、离散信道容量的计算(以输入符号等概分布为例) 我们利用函数

11、dmessage来求信源的熵，利用函数hemssage来求平均互信息量，并最终得到信道的容量。离散信道容量：信道容量是信道所能传送的最大的信息量。C=maxI(X;Y) （比特/码元）I(X;Y)=H(Y)H(Y/X)；3、信源编码过程(以Huffman编码为例) 对于所要求的信源，可以对其进行作为一元信源的哈夫曼编码并得到编码效率，相应的二元信源的哈夫曼编码及其编码效率。huffman编码方法如下：（1） 将信源信息呼号按其出现的概率大小依次排列；（2） 取两个概率最小的字母分别配以0和1两个码元，并将这两个概率相加作为一个新字母的概率，与未分配的二进制符号的字母重新排队；（3） 对重排后的

12、两个概率最小符号重复步骤（2）的过程；（4） 不断重复上述过程，直到最后两个符号配以0和1为止；（5） 从最后一级开始，向前返回得到各个信源符号所对应的码元序列，及相应的码字。Huffman编码的意义：将概率大的信息符号编以短的码字，概率小的符号配以长的码字，使得平均码字长度最短，冗余度减小。4、PCM编译码器设计及应用（1）PCM编码器系统 Simulink实现PCM技术的三个过程：采样、量化与编码。其中以saturation作为限幅器，将输入信号幅度值限制在PCM编码定义的范围内，Relay模块的门限设置为0，其输出即可作为PCM编码输出的最高位极性码。样值取绝对值后，以LookUp Ta

13、ble模块进行13折线压缩，并用增益模块将样值范围放大到0127,然后用间距为1的Quantizer进行四舍五入取整，最后将整数编码为7位二进制序列，作为PCM编码的低7位。（2）PCM译码器系统 PCM解码器中首先分离并行数据中的最高位（极性码）和7位数据，然后将7位数据转换为整数值，再进行归一化、扩张后与双极性码相乘得出解码值。（3）PCM编译码器整体系统5、数字信号基带传输系统设计（1） 编程实现常见基带信号的波形、码型转换包括：单极性非归零码、单极性归零码、双极性非归零码、双极性归零码、AMI码。单极性非归零码：用电平1来表示二元信息中的“1”，用电平0来表示二元信息中的“0”，电平在

14、整个码元的时间里不变，记作NRZ码。单极性归零码：它与单极性非归零码不同处在于输入二元信息为1时，给出得码元前半时间为1，后半时间为0，输入0则完全相同。双极性非归零码：它与单极性非归零码类似，区别在于双极性使用电平-1来表示信息0.双极性归零码：此种码型比较特殊，它使用前半时间1,、后半时间0来表示1；采用前半时间-1、后半时间0来表示信息0.因此它具有3个电平，严格的说是一种三元码（电平1,0，-1）。AMI码：全称是传号交替反转码。其编码规则是将消息码中的“1”交替变成“+1”和“-1”，将消息码中的“0”仍保持为“0”。（2） 波形、码型的功率谱分布的绘图，包括：单极性非归零码、单极性

15、归零码、双极性非归零码、双极性归零码、AMI码。数字基带信号一般是随机信号，因此分析随随机信号的频谱特性要用功率谱密度来分析。一般来说求解功率谱是一件相当困难的事情，但由于上述几种码型比较简单，因此可以求出功率谱。（3）编程实现基带 传输系统的误码率计算：二电平误码率计算。6、数字信号频带传输系统设计（1）2ASK调制/相干解调系统设计（2）2ASK调制/非相干解调系统设计（包络检波法）产生二进制振幅键控信号的方法，或称调制方法，主要有两种。第一种采用相乘电路法，用基带信号和载波相乘就得到了已调信号输出。第二种是采用开关电路法。在振幅键控中载波幅度是随着基带信号的变化而变化的。使载波在二进制基

16、带信号1或0的控制下通或断，即用载波幅度的有或无来代表信号中的“1”或“0”，这样就可以得到2ASK信号，这种二进制振幅键控方式称为通断键控（OOK）。其时域数学表达式为： （1）式中，A为未调载波幅度，为载波角频率，为符合下列关系的二进制序列的第n个码元： （2）综合式1和式2，令A1，则2ASK信号的一般时域表达式为：式中，Ts为码元间隔，为持续时间 Ts/2，Ts/2 内任意波形形状的脉冲（分析时一般设为归一化矩形脉冲），而就是代表二进制信息的随机单极性脉冲序列。采用相乘电路产生二进制ASK信号时，用基带信号A（t）与载波coswt相乘就得到已调信号输出。2ASK信号解调时可采用相干解调

17、法和包络检波法，包络检波法解调器是用整流器和低通滤波器构成的；相干解调中相乘电路需要有相干载波coswt,它必须从接收信号中提取，并且和接收信号的载波同频同相，所以这种方法比包络检波法要复杂得多。2ASK信号调制器方框图如下：2ASK信号的解调原理图如下：三、 主要功能的实现1、信源平均信息量的计算(以高斯分布的信源为例)（1）程序代码见附件（程序清单）（2）运行结果截图如下：2、离散信道容量的计算(以输入符号等概分布为例)（1）程序代码见附件（程序清单）（2）运行结果截图如下：3、信源编码过程(以Huffman编码为例)（1）程序代码见附件（程序清单）（2）运行结果截图如下：4、 PCM编译

18、码器设计及应用13折线近似的PCM编码器测试模型及仿真结果如下：13折线近似的PCM解码器测试模型及仿真结果如下：PCM编码器系统（子系统）设计如下：PCM译码器系统（子系统）设计如下：PCM编译码器整体设计如下：PCM编译码器运行结果如下：（1） 输入模拟语音信号观察其输出波形如下：（2） 输入数字信号观察其输出波形如下：5.1编程实现常见基带信号的波形、码型转换包括：单极性非归零码、单极性归零码、双极性非归零码、双极性归零码、AMI码。（1）单极性非归零码程序代码见附件（程序清单）运行结果截图如下：（2）单极性归零码程序代码见附件（程序清单）运行结果截图如下：（3）双极性非归零码程序代码见

19、附件（程序清单）运行结果截图如下：（4）双极性归零码程序代码见附件（程序清单）运行结果截图如下： （5）AMI码程序代码见附件（程序清单）运行结果截图如下：5.2波形、码型的功率谱分布的绘图，包括：单极性非归零码、单极性归零码、双极性非归零码、双极性归零码、AMI码。（1）单极性非归零码程序代码见附件（程序清单）运行结果截图如下：（2）单极性归零码程序代码见附件（程序清单）运行结果截图如下：（3）双极性非归零码程序代码见附件（程序清单）运行结果截图如下：（4）双极性归零码程序代码见附件（程序清单）运行结果截图如下：（5）AMI码程序代码见附件（程序清单）运行结果截图如下：5.3编程实现基带传输

20、系统的误码率计算：二电平误码率计算。程序代码见附件（程序清单）运行结果截图如下：6、数字信号频带传输系统设计（1）2ASK调制/相干解调系统设计程序代码见附件（程序清单）运行结果截图如下：a、 各信号波形b、 各信号频谱 （2）2ASK调制/非相干解调系统设计（包络检波法）程序代码见附件（程序清单）运行结果截图如下：a、各信号波形b、各信号频谱四、 功能模块的分析与流程图1、信息论基本运算模块（1）通过系统产生一个高斯随机信源，再求出它的平均信息量。（2）编写hmessage(The average mutual information)函数求出平均互信息，编写message(Discrete

21、 Entropy)函数，求出离散信息熵，调用hmessage函数和message函数得出hf和hx,最后信道容量c=hx-hf。（3）先编写huffman函数，对系列排序并求出huffman编码。调用huffman函数得出编码后的码字2、PCM编译码的设计及运用模块PCM编译码器的步骤框图如下:低通压缩抽 样量 化编 码低通扩张解 调解 码信道接 收输入输出PCM即脉冲编码调制，在通信系统中完成将语音信号数字化功能。PCM的实现主要包括三个步骤完成：抽样、量化、编码。分别完成时间上离散、幅度上离散、及量化信号的二进制表示。根据CCITT的建议，为改善小信号量化性能，采用压扩非均匀量化，有两种建

22、议方式，分别为A律和律方式，我国采用了A律方式，由于A律压缩实现复杂，常使用 13 折线法编码,采用非均匀量化PCM编码。A律压扩特性是连续曲线，A值不同压扩特性亦不同，在电路上实现这样的函数规律是相当复杂的。实际中，往往都采用近似于A律函数规律的13折线（A=87.6）的压扩特性。这样，它基本上保持了连续压扩特性曲线的优点，又便于用数字电路实现，本设计中所用到的PCM编码正是采用这种压扩特性来进行编码的。下表列出了13折线时的值与计算值的比较。0101按折线分段时的01段落12345678斜率16168421表中第二行的值是根据时计算得到的，第三行的值是13折线分段时的值。可见，13折线各段

23、落的分界点与曲线十分逼近，同时按2的幂次分割有利于数字化。 编码是把量化后的信号变换成代码，其相反的过程称为译码。当然，这里的编码和译码与差错控制编码和译码是完全不同的，前者是属于信源编码的范畴。 在现有的编码方法中，若按编码的速度来分，大致可分为两大类：低速编码和高速编码。通信中一般都采用第二类。编码器的种类大体上可以归结为三类：逐次比较型、折叠级联型、混合型。在逐次比较型编码方式中，无论采用几位码，一般均按极性码、段落码、段内码的顺序排列。下面结合13折线的量化来加以说明。 段落码 段内码段落序号段落码量化级段内码81111511111411107110131101121100610111

24、10111010105100910018100040117011160110301050101401002001300112001010001000100000 PCM编码输出的二进制序列中，每个样值用8位二进制码来表示，其中最高比特位是样值的正负极性，规定负极性为用0表示，正极性用1表示，接下来三位为段落码，最后四位是样值处于段落内16个均匀间隔上的间隔序号。五、 程序及Simulink调试1、 在离散信道容量的计算(以输入符号等概分布为例)引用程序文件hmessage.m文件中误将function r=hmessage(x,f,nx,my)搞成了function r=message(x,f

25、,nx,my)结果出现如下报错：? Input argument nx is undefined.Error in = hmessage at 3for i=1:nx改正后运行正常，出现平均互信息量hf 、此离散信源的平均信息量r、信道的信道容量c 的正确求值结果。2、 在波形、码型的功率谱分布的绘图时把相应引用程序文件中的grid 值设置错误（以单极性非归零码威力）导致无法正常显示波形、码型的功率谱分布。例如在单极性非归零码的求波形、码型的功率谱分布时snrz.m文件代码中grid 设置为300而在单极性非归零码波形、码型的功率谱分布绘图的主程序中M设置为64结果出现以下报错：? Error

26、 using = plotVectors must be the same lengths.Error in = Untitleddanfeigonglv at 14plot(f,10*(log10(ys);将snrz.m文件代码中grid 设置为64后程序运行正常出现正确的功率谱分布图像。3、在解码子系统的Block Parameter:Mux1中Number of inputs设置参数时将其设置为8，结果出现报错，分析其原因为Mux1只能接收7位二进制序列应将其设置为7。将其设置为7后报错消失，系统运行正常能够显示出输入模拟话音信号的正确输出波形和输入数字信号的正确输出波形。六、心得体会这

27、次通信原理的课程设计本来计划是两周完成，但是由于一些因素，我们压缩了一半的时间，这对我们来说不能不说是一个挑战。时间紧任务重，无时无刻不被这个课程设计搞得神经紧绷。特别是刚开始的时候，彭老师让我们去图书馆查资料，说真的我真的不知道如何下手了。逼着自己去想，去找，结果发现图书馆的资料真是“神似”啊！车到山前必有路，我们从图书馆、网上找到一些资料，然后仔细研究讨论，渐渐地有了一些眉目。俗话说万世开头难。接下来，总该做的顺利一些了吧。结果可以引用方志敏的一句话：“你们想错了。”面对资料里的各种信息，我不断的甄别、调整。按资料上的不行，我就考虑到底是哪些地方出错了，然后一步步调。程序不对就改，数据不对

28、就调。就这样一个个题目慢慢的有了结果，虽然速度慢了点，但还是很高兴。不仅如此，这次课程设计让我再一次感受到掌握matlab这个软件对我们通信专业学生的重要性。因为我们通信原理的实验是在通信实验室用通信实验箱做的，本以为课程设计会和实验一样，结果课程设计却是用matlab设计通信系统。虽然这学期我们学了matlab但用它来设计通信系统却是头一次。在接下来的几天里我充分体会到matlab这个软件的功能之强大，运用之宽广。尤其是PCM编译码器系统用Simulink仿真，有许多模块是以前没有接触过的，所以得慢慢摸索，在摸索中学着如何运用。除此之外，从这次课程设计中我还得到了一些其它的体会。比如在编写波

29、形、码型的功率谱分布的绘图的程序时，我就纠结于是在时域显示还是在频域显示。最终我还是选择了在频域显示。这让我体会到有时到达一个目的地可以选择不同的路径。在后来的2ASK调制/相干解调器和非相干解调系统（包络检波法）设计中我同样面临着不同的选择，选择编程实现或是用simulink对系统建模。为了简化问题我选择了编程。有时就是这样，我们要选择最简的路劲去完成任务。通过本次课程设计我明白了，我们掌握的理论知识还远远不足以满足我们实际操作的需要。我们要加强理论知识的学习与拓展，由点到线，由线到面，不断地学习理论知识。但是仅仅依靠理论知识，那只是纸上谈兵。就像这次课设中的一些问题，咋一看和书上的差不多，

30、但如果直接把书上的拿来就用，会发现根本就不是那回事。所以我们必须及时的把理论知识运用到我们的实践中去。在实践中发展和完善理论。最终把理论和实践结合起来，把自己锻炼成一个合格实用的通信人才。七、评分表课程设计评分表课程名称： 通信原理 项 目评 价设计方案的合理性与创造性设计与调试结果设计说明书的质量答辩陈述与回答问题情况课程设计周表现情况综合成绩 教师签名： 日 期： 八、附件（程序清单）1、信源平均信息量的计算(以高斯分布的信源为例)syms x u ou=3;o=4;f=1/(u*sqrt(2*pi)*exp(-(x-o)2/2*u2);t=-f*log(f)/log(2);r=int(t

31、,-inf,inf);disp(平均信息量为)r=double(r)2、离散信道容量的计算(以输入符号等概分布为例)（1）message.m文件代码如下：function r=message(x,n)r=0;for i=1:n r=r-x(i)*log(x(i)/log(2);end;disp(此离散信源的平均信息量为:)r（2） hmessage.m文件代码如下：function r=hmessage(x,f,nx,my)sum=0;for i=1:nx for j=1:my t=f(i,j)*x(i); if t=0 sum=sum-t*log(f(i,j)/log(2); end; en

32、d;end;r=sum;disp(平均互信息量为:)double(r)（3） 主程序代码如下： clear;x=0.25,0.25,0.25,0.25;f1=1/2 1/2 0 0;0 1/2 1/2 0;0 0 1/2 1/2;1/2 0 0 1/2;hf=hmessage(x,f1,4,4)hx=message(x,4)disp(信道的信道容量为:) c=hx-hf3、 信源编码过程(以Huffman编码为例)(1) huffman.m文件代码如下：function h,l=huffman (p)if (length(find(p10e-10) error(Not a prob.vecto

33、r,component do not add to 1) endn=length(p); q=p; m=zeros(n-1,n);for i=1:n-1 q,l=sort(q); m(i,:)=l(1:n-i+1),zeros(1,i-1); q=q(1)+q(2),q(3:n),1; endfor i=1:n-1 c(i,:)=blanks(n*n);endc(n-1,n)=0; c(n-1,2*n)=1; for i=2:n-1 c(n-i,1:n-1)=c(n-i+1,n*(find(m(n-i+1,:)=1)-(n-2):n*(find(m(n-i+1,:)=1); c(n-i,n)=

34、0; c(n-i,n+1:2*n-1)=c(n-i,1:n-1); c(n-i,2*n)=1; for j=1:i-1 c(n-i,(j+1)*n+1:(j+2)*n)=c(n-i+1,n*(find(m(n-i+1,:)=j+1)-1)+1:n*find(m(n-i+1,:)=j+1); end; end for i=1:n h(i,1:n)=c(1,n*(find(m(1,:)=i)-1)+1:find(m(1,:)=i)*n); ll(i)=length(find(abs(h(i,:)=32); end l=sum(p.*ll);（2）主程序如下：p=1/2,1/4,1/8,1/16,1

35、/16 h,l=huffman(p)4、编程实现常见基带信号的波形、码型转换：（1） 单极性非归零码snrz.m文件代码如下：function y=snrz(x)grid =300;t=0:1/grid:length(x);for i=1:length(x); if(x(i)=1), for j=1:grid, y(i-1)*grid+j)=1; end else for j=1:grid, y(i-1)*grid+j)=0; end endendy=y,x(i);M=max(y);m=min(y)plot(t,y);axis(0,i,m-0.1,M+0.1);title(1 0 0 1 1

36、0 0 0 0 1 0 1);主程序如下：t=1 0 0 1 1 0 0 0 0 1 0 1;snrz(t);(2)单极性归零码srz文件代码如下：function y=srz(x)grid =300;t=0:1/grid:length(x);for i=1:length(x); if(x(i)=1), for j=1:grid/2, y(grid/2*(2*i-2)+j)=1; y(grid/2*(2*i-1)+j)=0; end else for j=1:grid/2, y(i-1)*grid+j)=0; end endendy=y,x(i);M=max(y);axis(0 12 0 1.

37、2)subplot(2,1,1);plot(t,y);title(1 0 0 1 1 0 0 0 0 1 0 1);主程序如下：t=1 0 0 1 1 0 0 0 0 1 0 1;srz(t);（3） 双极性非归零码dnrz函数代码如下：function y=dnrz(x)grid =300;t=0:1/grid:length(x);for i=1:length(x); if(x(i)=1), for j=1:grid, y(i-1)*grid+j)=1; end else for j=1:grid, y(i-1)*grid+j)=-1; end endendy=y,x(i);M=max(y);m=min(y)plot(t,y);axis(0,i,m-0.1,M+0.1);title(1 0 0 1 1 0 0 0 0 1 0 1);主程序如下：t=1 0 0 1 1 0 0 0 0 1 0 1;dnrz(t);(4) 双极性归零码drz.m文件代码如下：function y=drz(x)grid=300;t=0:1/grid:length