连续时间系统的时域分析——求零输入响应和零状态响应

1、成绩评定表学生姓名丛蕾蕾班级学号16专业通信工程课程设计题目求零输入响应和零状态响应评语组长签字：成绩日期20 年 月日课程设计任务书学院信息科学与工程专业通信工程学生姓名丛蕾蕾班级学号1103060106课程设计题目连续时间系统的时域分析求零输入响应和零状态响应内容与要求：1s学习Matlab软件知识与应用2、学习并研究信号可视化与时域运算3、利用Matlab编程，完成常用连续信号与信号的时移、反褶、尺度 变换4、写出课程设计报告打印程序，给出运行结果 进度安排：周一、周二：1、布置课程设计任务、要求2、学习Matlab软件知识与应用周三、周四：1、利用Matlab编程，完成相应的信号分析与

2、处理课题2、上机编程、调试3、撰写课程设计报告书周五：答辩，上交报告指导教师：20：1年 月 日专业负责人：201 年 月日学院教学副院长：201 年 月日1 .弓I言12 Matlab 入门22.1 Matlab7.0 介绍22.2 利用Matlab7.0编程完成习题设计33 Matlab7.0实现连续时间系统的时域分析3.1 常用连续时间信号的类别与原理43.2 编程设计与实现43.3 运行结果与其分析7结论20参考文献211 / 141.引言人们之间的交流是通过消息的传播来实现的，信号则是消息的表现形式，消 息是信号的具体内容。本文概述了信号仿真系统的需求、总体结构、基本功能。重点介绍了

3、利用 Matlab软件设计实现信号仿真系统的基本原理与功能，以与利用Matlab软件提 供的图形用户界面(Graphical User Interfaces ,GUI)设计具有人机交互、界面友 好的用户界面。本文采用Matlab的图形用户界面设计功能，开发出了各个实验界 面。在该实验软件中，集成了信号处理中的多个实验，应用效果良好。本系统是 一种演示型软件用可视化的仿真工具,以图形和动态仿真的方式演示部分基本信 号的传输波形和变换，使学习人员直观、感性地了解和掌握信号与系统的基本知 识。近年来，计算机多媒体教育手段的运用逐步普与，大量优秀的科学计算和系 统仿真软件不断涌现，为我们实现计算机辅助

4、教学和学生上机实验提供了很好的 平台。通过对这些软件的分析和对比，我们选择MATLAB语言作为辅助教学工具, 借助MATLAB强大的计算能力和图形表现能力，将信号与系统中的概念、方 法和相应的结果，以图形的形式直观地展现给我们，大大的方便我们迅速掌握和 理解老师上课教的有关信号与系统的知识。MATLAB 7.0是一个包含大量计算算法的集合。其拥有600多个工程中要用 到的数学运算函数,可以方便的实现用户所需的各种计算功能。函数中所使用的 算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。 在通常情况下，可以用它来代替底层编程语言，如C和C+ o在计算要求相同 的情况下，使用

5、MATLAB的编程工作量会大大减少。MATLAB 7.0的这些函数集 包括从最简单最基本的函数到诸如矩阵特征向量、快速傅立叶变换的复杂函数。 函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程与偏微 分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问 题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组 操作以与建模动态仿真等。作为信号与系统的基本分析软件之一，利用MATLAB进行信号与系统的分析 与设计是通信以与信息工程学科的学生所要掌握的必要技能之一。通过学习并使 用MATLAB语言进行编程实现课题的要求，对学生能力的培养极为重要。尤其会

6、 提高综合运用所学理论知识进行分析问题、解决问题的能力，也便于将理论知识 与实践相结合，并得以更好地掌握信号分析与处理的基本方法与实现。这也将为 后续相关的课程学习打下一定的基础,从而在以后相关课程设计与分析的时候达 到对MATLAB的熟练应用与融会贯通。2.Matlab 入门2.1 Matlab7.0 介绍MATLAB7.02009版本，高级技术计算语言和交互式环境可以较使用传统的 编程语言（如C、C+和Fortran）更快地解决技术计算问题。MATLAB 7.0主要功能为交互式工具可以按迭代的方式探查、设计与求解问 题此高级语言可用于技术计算此开发环境可对代码、文件和数据进行管理各种工具

7、可用于构建自定义的图形用户界面各种函数可将基于MATLAB的算法与外部应 用程序和语言（如 用C+ +、FortranJava. COM以与Microsoft Excel）集成 数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以与数值积分等二维 和三维图形函数可用于可视化数据；特点：1）高效的数值计算与符号计算功能， 能使用户从繁杂的数学运算分析中解脱出来；2）具有完备的图形处理功能，实现 计算结果和编程的可视化友好的用户界面与接近数学表达式的自然化语言， 使学者易于学习和掌握；4）功能丰富的应用工具箱（如信号处理工具箱、通信工 具箱等），为用户提供了大量方便实用的处理工具。MATLAB

8、7.0由一系列工具组成。这些工具方便用户使用MATLAB的函数和 文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历 史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件 的浏览器。随着MATLAB 7.0的商业化以与软件本身的不断升级，MATLAB 7.0的 用户界面也越来越精致，更加接近Windows的标准界面人机交互性更强，操 作更简单。而且新版本的MATLAB 7.0提供了完整的联机查询、帮助系统，极大 的方便了用户的使用。简单的编程环境提供了比较完备的调试系统，程序不必经 过编译就可以直接运行，而且能够与时地报告出现的错误与进行出错原因分析。

9、综上，在进行信号的分析与仿真时，MATLAB7.0无疑是一个强大而实用的 工具。尤其对于信号的分析起到了直观而形象的作用，非常适合与相关课题的研 究与分析。3. Matlab7.0实现连续时间系统的时域分析3.1 连续时间系统的零输入响应3.1.1 零输入响应的求解方法描述n阶线性时不变(LTI)连续系统的微分方程为：dny dxydydmuducii :h %：h+ ci -:F ci x y = bx! +bF b /1 dt 2 dt 41 dtm m dt 向已知y与各阶导数的初始值为y(0),/(0),/，求系统的零输入响应。建模如下：当lit系统的输入为零时，其零输入响应为微分方程

10、的齐次解(即令微分方程的等号右端为零)，其形式为(设特征根均为单根)y=+ C2/M +.+/即其中,pn是特征方程aAn+a2Xn l+-i+anX+an=0的根,它们可以用root(a)语句求得。各系数 由y与其各阶导数的初始值来确定。对此有下列方程：a + G + , +C = %p + P2c2 + ,+ P“C =力0川一匕+审一匕+。:匕=。1%写成矩阵形式为：PJ匕+ P2n-1C2+-+ PE成kD11，。PiPTPl P2”T即 V，C=Y。其解为：OVY。式中c = G=% Dy. DtPiPl_ n-1PlPnV为范德蒙矩阵,在matlab的特殊矩阵库中有vandero3

11、.1.2 零输入响应的编程设计与实现以下面式子为例：y + 2y(r) + 3) = 4/(r)- 6/(r)初始条件为y(0_)= 2,y(0_)= 7MATLAB 程序：a=input(输入分母系数 a=al1a2,.=);n=length(a)-l;Y0=input(输入初始条件向量 Y0=y0,Dy0,D2y0.=);p=roots(a);V=rot90(vander(p);c=VY0,;dt=input(dt=);te=input(,te=1);t=0;dtte;y=zeros(l,length(t);for k=l:n y=y+c(k)*exp(p(k)*t);end-5 - /

12、14plot(t,y):gnd xlabel(t) ;ylabel(y);title。零输入响应)程序运行结果：用这个通用程序来解一个三阶系统，运行此程序并输入a=1,2,3 Y0=2,7 dt=0.002 te=7结果如图3.1所示。图3.1零输入响应的波形图3.1.3 零输入响应的结果分析根据图可以分析零输入响应，它的起始值与输入函数无关，只与它的初始状态值有关，其起始值等于y(0_)的值。随着时间的推移，最后零输入响应的值无限 的趋近于0。3.2 连续时间系统的零状态响应3.2.1 零状态响应的求解方法我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述，NM(，)+40y)=

13、4/(，)+4/ )+4/ )+o/(f)可 用4=。3，。2吗，4乃=回也，也,输入函数=/),得出它的冲激响应八, 再根据LTI系统的零状态响应y (t)是激励u (t)与冲激响应h (t)的卷积积分。 注意，如果微分方程的左端或右端表达式中有缺项，则其向量a或b中的对应元 素应为零，不能省略不写，否则出错。例如：求函数的零状态响应，已知：)，”(/) + 2/(0 + 3X0 = 4/() - 6/(r)与初始状态 %(一)= a(.)= 。输入函数 y = Sin(2 *f) + cos * t) o 建模：先求出系统的冲激响应，写出其特征方程万+2几+ 3 = 0求出其特征根为pl和

14、P2,与相应的留数口，2 ；则冲激响应为h(t) = tep1 + r2ep-输入y (t)可用输入u (t)与冲激响应h (t)的卷积求得。3.2.2 零状态响应的编程设计与实现MATLAB 程序：a=input，输入分母系数 a=al,a2,.b=input(,输入输入信号系数 b=bl,b2,.=1);dt=input(dt=);te=input(,te=1);t=O:dtte;u二input(输入函数 u=1);te=t(end);dt=te/(length(t)-l);rtp,k=residue(bla);h=r(l)*exp(p(l)*t)+r(2)*exp(p(l)*t);su

15、b plot(2,l, l),plot(t, h);gridEe(冲激函数)y = conv(u,h)7t;subplot(2,12)tplot(ty(l：length(t);gridEe(零状态响应)程序运行结果执行这个程序，取 a=123b=4,6dt=0.001 te=7得出图形如图3.2所示。-7- / 14图3.2零状态响应的波形图3.2.3 零状态响应的结果分析由于初始状态为零，所以零状态的起始值也为零，即h(t)包含了连续系统的固 有特性，与系统的输入无关。只要知道了系统的冲激响应，即可求得系统在不同 输入时产生的输出。因此，求解系统的冲激响应h对进行连续时间系统的分析具 有非常

16、重要的意义-9 - / 14总结本次的信号与系统课程设计，我的课设题目是求连续时间系统的零输入与零 状态响应，需要用到MATLAB进行波形仿真。由于我对MATLAB这个仿真软件并不 是很熟悉，所以在进行课设的过程中，遇到了很多问题，不过最后还是顺利完成, 虽然做的程序并不是很完美,但是我还是从中学会了很多。设计过程中，出现了各种各样的问题，有些是由对软件的使用不熟悉引起的, 有的是由对零状态响应和零输入响应的求解方法不熟悉弓I起的。但是最后我还是 找到了解决的办法。比如说,在我用系统求解法求解零状态响应和零输入响应时, 由于求解过程复杂.用MATLAB软件进行仿真时，编程很难，所以我采用了矩阵

17、 求解的方法，可以在编程时直接调用特殊矩阵库中的Vander进行系统求解，这 样编程就容易了许多。在整个设计过程中我懂得了许多东西，也培养了独立思考和设计的能力，树 立了对知识应用的信心，相信会对今后的学习工作和生活有非常大的帮助，并且 学习到了 MATLAB软件主要功能交互式工具可以按迭代的方式探查、设计与求解 问题此高级语言可用于技术计算此开发环境可对代码、文件和数据进行管理各种 工具可用于构建自定义的图形用户界面各种函数可将基于MATLAB的算法与外 部应用程序和语言（如 C、C+ +、Fortrans Java COM 以与 Microsoft Excel） 集成数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以与数值积分等 二维和三维图形函数可用于可视化数据。课程设计不仅是对前面所学知识的一种检验，也是对自己能力的一种提高， 通过这次课程设计使自己明白了原来的光靠书面知识是非常欠缺的，还是需要不 断的实践巩固。因此在以后的工作和生活中还是应该不断的学习，努力提高自己 的知识和综合素质。本次的课程设计让我对信号与系统这一门学科有了更深层次的理解，在分析 并解决问题的过程中