基于MATLAB的线性常系数差分方程求解

1、数字信号处理课程设计题目： 基于MATLAB的线性常系数差分方程求解 学院： 专业： 班级： 学号： 姓名： 指导教师： 目 录摘 要1第一章 背景31.1 背景知识31.2 数字信号课程特点31.3 软件介绍41.4 MATLAB及数字信号处理4第二章 设计目的及要求62.1 设计目的62.2 课程设计的内容要求72.2.1 设计要求7第三章 设计任务8第四章 设计原理94.1 差分与差分方程94.2 线性常系数差分方程144.3 线性常系数差分方程的求解15第五章 设计过程165.1 用MATLAB求解差分方程16第六章 设计代码及结果186.1 MATLAB源程序186.2 程序运行结果

2、206.3 比较结果总结24第七章 收获与体会25致 谢27参考文献28摘 要数字信号处理分析了数字信号处理课程的重要性及特点，为了帮助学生理解与掌握课程中的基本概念、基本原理、基本分析方法，提出了用MATLAB进行数字信号处理课程设计的思路，并阐述了课程设计的具体方法、步骤和内容。MATLAB语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，MATLAB功能强大、简单易学、编成效率高，深受广大科技工作者的喜爱，特别是MATLAB还具有信号分析工具箱，不需具备很强的编程能力，就可以很方便地进行语音信号分析、处理和设计。线性常系数差分方程求解是数字信号处理课程中常出现的课题，也是现代科学

3、中值得深入研究的一个课题 本文介绍了线性常系数差分方程的基本概念，论述了其求解方法，并用MATLAB具体实现了线性常系数差分方程的求解。基于MATLAB的线性常系数差分方程求解主要是用MATLAB作为工具平台，设计中涉及到差分方程的递推求解以及用filter对系数向量的归一化等等。通过数字信号处理课程的理论知识的综合运用，从实践上初步实现对数字信号的处理。关键字：MATLAB，线性常系数差分方程，数字信号处理。Abstract Digital signal processing analysis of digital signal processing course of the import

4、ance and features, in order to help the students to understand and grasp basic concepts, basic principles, basic analysis method, is put forward with the MATLAB for digital signal processing curriculum design, curriculum design and describes the specific methods, steps and content.The MATLAB languag

5、e is widely used in engineering calculation and the numerical analysis in the field of advanced language, MATLAB powerful, easy to learn, a high efficiency, by the vast number of scientific workers favorite, especially the MATLAB also has a signal analysis toolbox, does not need to have very strong

6、ability of programming, can be very convenient for the analysis of speech signal, processing and design.Linear constant coefficient differential equation is a digital signal processing program that often appear in the topic, as well as modern science and worthy of in-depth study of a topicThis paper

7、 introduces the linear constant coefficient differential equation basic concept, discussed the solution method, and MATLAB the specific realization of linear constant coefficient difference equation.MATLAB based on the linear constant coefficient difference equation MATLAB is mainly used as a tool p

8、latform, design relate to differential equation recursive solution and the use of filter on the coefficient vector is normalized and so on. Through the course of digital signal processing theory, the integrated use of knowledge, from the practice of preliminary implementation of digital signal proce

9、ssing.Key words: MATLAB, Linear constant coefficient differential equation,Digital signal processing.第一章 背景1.1 背景知识数字信号处理（Digital Signal Processing,简称DSP）是一门设计许多学科而又广泛应用于许多领域的新兴学科。DSP有两种含义：Digital Signal Processing（数字信号处理）、Digital Signal Processor（数字信号处理器）。我们常说的DSP指的是数字信号处理器。数字信号处理器是一种适合完成数字信号处理运算的处

10、理器。20世纪60年代以来，随着计算机和信息技术的飞速发展，数字信号处理技术应运而生并得到迅速的发展。在过去的二十多年的时间里，数字信号处理已经在通信等领域得到极为广泛的应用。数字信号处理是利用计算机专用处理设备，以数字形式对信号进行采集、变换、滤波、估值增强、压缩、识别等处理，以得到符合人们需要的信号形式。它是以众多学科为理论基础的，它所涉及的范围极其广泛。例如，在数学领域，微积分、概率统计、随即过程、数值分析等都是数字信号处理的基本工具，与网络理论、信号与系统、控制论、通信理论、故障诊断等也密切相关，近来新兴的一些学科，如人工智能、模式识别、神经网络等，都与数字信号处理密不可分。可以说，数

11、字信号处理是把许多经典的理论体系作为自己的理论基础，同时又使自己成为一系列新兴学科的理论基础。1.2 数字信号课程特点数字信号处理课程是一门理论和技术发展十分迅速、应用非常广泛的前沿性学科，他的理论性和实践性都很强，他的特点是：(1)要求的数学知识多，包括高等代数、数值分析、概率统计、随机过程等。(2)要求掌握的基础知识强，网络理论、信号与系统是本课程的理论基础。(3)与其他学科密切相关，即与通信理论、计算机、微电子技术不可分，又 是人工智能、模式识别、神经网络等新兴学科的理论基础之一。选择用MATLAB进行课程设计:MATLAB语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，MAT

12、LAB功能强大、简单易学、编程效率高，深受广大科技工作者的欢迎。特别是MATLAB还具有信号分析工具箱，不需具备很强的编程能力，就可以很方便地进行信号分析、处理和设计。1.3 软件介绍 MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数字数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现

13、算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、经融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解答问题要比用C,FORTRAN等语言完成相同的事情简捷的多，并且mathworks也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C,FORTRAN，C+，JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用

14、户可以直接进行下载就可以用。1.4 MATLAB及数字信号处理MATLAB是矩阵实验室之意。除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。MATLAB的基本数据单位是矩阵，它的指令表达式与数学，工程中常用的形式十分相似，故用MATLAB来解决问题要比用C,FORTRAN等语言完全相同的事情简捷得多，可以将自己编写的实用程序导入到MATLAB函数库中方便在新的版本中也加入了对C，FORTRAN，c+，JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。信号时数字信号处理领域中最基本、最重要的

15、概念。简单地说，信号就是信息的载体，是信息的物理体现。信号既可以分为时间连续、幅度也连续的模拟信号和时间和幅度上都经过量化的数字信号，也可以划分为连续时间信号和离散时间信号。几乎在科学技术的每一领域，为了信号的提取，都要进行信号处理，就是以数值计算的方法对信号进行采集，变换，综合，估计，与识别的加工处理过程，借以达到提取信息和便与应用的目的。随着计算机和信息科学的飞速发展，信号处理已经逐渐发展为一门独立的学科，是信息科学的重要组成部分。在语音处理，雷达，航空，航天，地质勘探，通信，生物医学工程等众多领域得到了广泛应用。MATLAB软件，在数字信号处理方面具有得天独厚的优势。利用目录下的系统函数

16、，用户可以实现波形的产生，信号的变换，滤波，功率谱估计，系统设计与稳定性分析，小波信号分析等众多功能。本文即是以数字信号处理的理论基础，应用MATLAB软件求解线性常系数差分方程的一个具体事例。第二章 设计目的及要求2.1 设计目的 1. 全面复习课程所学理论知识，巩固所学知识重点和难点，将理论与实践很好的结合起来。2. 掌握信号分析与处理的基本方法与实现3. 提高综合运用所学知识独立分析和解决问题的能力4. 熟练使用一种高级语言进行编程实现课程设计是教学的最后一个步骤，课程设计有利于基础知识的理解，我们掌握了基础知识和基本技能，但是要真正接触才能真正理解课程的深入部分；还有礼物逻辑思维的锻炼

17、，在许多常规学科的日常教学中，我们不难发现这样一个现象，不少学生的思维常常处于混乱的状态，写起作文来前言不搭后语，解起数学题来步骤混乱，这些都是缺乏思维训练的结果，所以我们可以通过实践来分析问题、解决问题、预测目标等；同时也有利于其他学科的整合，例如我们这次的课程设计就是要运用MATLAB软件的帮助才能实现；最重要的有利于治学态度的培养，在课程设计中，我们可能经常犯很多小错误，可能要通过好几次的反复修改、调试才能成功，但这种现象会随着学习的深入而慢慢改观。这当中就有一个严谨治学、一丝不苟的科学精神的培养，又有一个不拍失败、百折不挠品格的锻炼。数字信号处理课程设计实在学生完成数字信号处理和MAT

18、LAB的结合后基于实验以后开设的。本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现线性常系数差分方程求解。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序编写与调试、写课程设计报告等步骤，使学生进一步掌握数字信号处理课程的基本理论、基本方法和基本技术；使学生增进对MATLAB的认识，利用MATLAB加深对理论知识的理解；使学生了解和掌握使用MATLAB的应用过程和方法，为以后的设计打下良好基础培养学生能根据设计要求，进行理论知识分析、设计方法总结、典型实例设计等方面的设计综合能力；使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。随

19、着信息科学和计算机技术的迅速发展，现代信号处理的理论与应用得到飞跃式的发展，形成一门极其重要的学科。特别是数字信号处理也已经成为高等学校相关专业的必修课程。为了更好地将数字信号处理的理论付诸实践，此次课程设计是一个很好的契机，通过仿真实验，让大家初步了解信号处理的分析法和实现方法。现代信号处理是一门一算法为核心，理论和实践性较强的学科。是电子信息工程、通信工程专业、电子信息科学与技术专业的一门重要的专业基础课。数字信号处理课程是在学习完数字信号处理的相关理论后，进行的综合性训练课程，其目的是：使学生进一步巩固数字信号处理的基本盖帘、理论、分析方法和实现方法；增强学生应用MATLAB语言编写数字

20、信号处理的应用程序及分析、解决实际问题的能力。2.2 课程设计的内容要求 设计要求 1. 巩固所学的专业技术知识；2. 提高综合运用所学理论知识独立分析和解决问题的能力；3. 进一步提高程序设计及调试能力；4. 更好地将理论与实践相结合；5. 学习和掌握科学研究资料检索的方法，学习对已有资料进行消化总结的方法6. 学习撰写科学报告的基本方法；7. 本设计要求分组合作完成；8. 上机前提前熟悉使用课程设计实验平台，掌握其基本的操作方法；9. 上机前了解课程设计实验平台的源代码，掌握其程序结构及在此平台上添加处理程序的方法；10. 设计过程中详细记录产生的图形、参数、数据等，用于编写课程设计报告。

21、本次课程设计的部分课题用到了本科数字信号处理的几乎所有知识，既能帮助学生对知识点的融会贯通，又使学生感到学有所用，培养进一步深入学习信号处理的兴趣。 设计内容 1. 课程设计题目和题目涉及要求2. 设计思想和系统分析功能分析3. 设计中关键部分的理论分析与计算，关键模块的设计思路4. 程序代码清单5. 测试数据、测试输出结果，及必要的理论分析和比较6. 总结，包括设计过程中遇到的问题和解决办法，设计心得与体会等7. 致谢8. 参考文献第三章 设计任务题目：基于MATLAB的线性常系数差分方程求解1、自行产生一个序列，要求：（1）对序列进行差分运算，并画出差分序列的时域波形图；2、已知一个二阶线

22、性常系数差分方程用下式表示：y(n)+a1y(n-1)+a2y(n-2)= b0x(n)+b1x(n-1)+b2x(n-2)，要求：（1） 参数a1、a2、b0、b1、b2由运行时输入；（2） 已知输入，画出x(n)的时域波形图；（3） 求出x(n)的共轭对称分量xe(n)和共轭反对称分量xo(n)，并分别画出时域波形图；（4） 初始条件由运行时输入，求输出y(n)，并画出其波形；（5） 对于不同的初始条件，分析其输出是否一致，从中得出什么结论第四章 设计原理 4.1 差分与差分方程与连续时间信号的微分及积分运算相对应，离散时间信号有差分及序列求和运算。设有序列f(k),则称,f(k+2),f

23、(k+1),f(k1),f(k2),为f(k)的移位序列。序列的差分可以分为前向差分和后向差分。一阶前向差分定义为 （3.11）一阶后向差分定义为 （3.12）式中和称为差分算子。由式（3.11）和式（3.12）可见，前向差分与后向差分的关系为 （3.13）二者仅移位不同，没有原则上的差别，因而它们的性质也相同。此处主要采用后向差分，并简称其为差分。由查分的定义，若有序列、和常数，则 （3.14）这表明差分运算具有线性性质。 二阶差分可定义为 （3.15）类似的，可定义三阶、四阶、n阶差分。一般地，n阶差分 （3.16）式中 （3.17）为二项式系数序列f(k)的求和运算为 （3.18）差分方

24、程是包含关于变量k的未知序列y(k)及其各阶差分的方程式，它的一般形式可写为 （3.19a）式中差分的最高阶为n阶，称为n阶差分方程。由式（3.16）可知，各阶差分均可写为y(k)及其各移位序列的线性组合，故上式常写为 （3.19b）通常所说的差分方程是指式（3.19b）形式的方程。若式（3.19b）中，y(k)及其各移位序列均为常数，就称其为常系数差分方程；如果某些系数是变量k的函数，就称其为变系数差分方程。描述LTI离散系统的是常系数线性差分方程。差分方程是具有递推关系的代数方程，若一直初始条件和激励，利用迭代法渴求的差分方程的数值解。 差分方程的经典解一般而言，如果但输入单输出的LTI系

25、统的激励f(k)，其全响应为y(k)，那么，描述该系统激励f(k)与响应y(k)之间关系的数学模型式n阶常系数线性差分方程，它可写为 （3.110a）式中、都是常数。上式可缩写为 （3.110b）与微分方程的经典解类似，上述差分方程的解由齐次解和特解两部分组成。齐次解用表示，特解用表示，即 （3.111）a.齐次解当式（3.110)中的f(k)及其各移位项均为零时，齐次方程 （3.112）的解称为齐次解。首先分析最简单的一阶差分方程。若一阶差分方程的齐次方程为 （3.113）它可改写为y(k)与y(k1)之比等于a表明，序列y(k)是一个公比为a的等比级数，因此y(k)应有如下形式 （3.11

26、4）式中C式常数，有初始条件确定。对于n阶齐次差分方程，它的齐次解由形式为的序列组合而成，将代入到式（3.112），得由于C0，消去C；且0，以除上式，得（3.115）上式称为差分方程式（3.110）和式（3.112）的特征方程，它有n个根，称为差分方程的特征根。显然，形式为的序列都满足式（3.112），因而它们是式（3.110）方程的齐次解。依特征根取值的不同，差分方程齐次解的形式见表31,其中、等为待定常数表31 不同特征根所对应的齐次解特征根齐次解单实根重实根一对共轭复根重共轭复跟b.特解特解的函数形式与激励的函数形式有关，表32列出了集中典型的激励f(k)所对应的特解。选定特解后代入原

27、差分方程，求出其待定系数等，就得出方程的特解。表32 不同激励所对应的特解激励特解 所有特征根均不等于1时 当有重等于1时的特征根时 当不等于特征根时 当是特征单根时 当是重特征根时或所有特征根均不等于c.全解 式（3.110）的线性差分方程的全解是齐次解与特解之和。如果方程的特征根均为单根，则差分方程的全解为（3.116）如果特征根为重根，而其余nr个特征根为单根时，差分方程的全解为（3.117）式中各系数由初始条件确定。如果激励信号是在k=0时接入的，差分方程的解适合于k0。对于n阶差分方程，用给定的n个初始条件y(0),y(1),y(n1)就可确定全部待定系数。如果差分方程的特解都是单根

28、，则方程的全解为式（3.116），将给定的初始条件y(0),y(1),y(n1)分别代入到式（3.116），可得（3.118）由以上方程可求得全部待定系数。 零输入响应系统的激励为零，仅由系统的初始状态引起的响应，称为零输入响应，用表示。在零输入条件下，式（3.110）等号右端为零，化为齐次方程，即 （3.125）一般设定激励是在k=0时接入系统的，在k0时，激励尚未接入，故式（3.125）的几个初始状态满足 （3.126）式（3.126）中的y(1),y(2),y(n)为系数的初始状态，由式（3.125）和式（3.126）可求得零输入响应。 零状态响应当系统的初始状态为零，仅由激励f(k)所

29、产生的响应，称为零状态响应，用 表示。在零状态情况下，式（3.110）仍是非齐次方程，其初始状态为零，即零状态响应满足（3.130）的解。若其特征根均为单根，则其零状态响应为 （3.131）式中为待定常数，为特解。需要指出，零状态响应的初始状态为零，但其初始值不一定等于零。4.2 线性常系数差分方程4.21 一个N 阶线性常系数差分方程可用下式表示： （）或者 （）式中，x（n）和y(n)分别是系统的输入序列和输出序列，ai和bi均为常系数，式中y(n-i)和x(n-i)项只有一次幂，也没有相互交叉相乘项，故称为线性常系数差分方程。差分方程的阶数是用方程y(n-i)项中i的最大取值与最小取值之

30、差确定的。在（）式中，y(n-i)项i最大的取值N，i的最小取值为零，因此称为N阶差分方程。4.3 线性常系数差分方程的求解已知系统的输入序列，通过求解差分方程可以求出输出序列。求解差分方程的基本方法有以下三种：（1） 经典解法。这种方法类似于模拟系统中求解微分方程的方法，它包括齐次解与特解，由边界条件求待定系数，上节已作简单介绍，这里不作介绍。（2） 递推解法。这种方法简单，且适合用计算机求解，但只能得到数值解，对于阶次较高的线性常系数差分方程不容易得到封闭式（公式）解答。（3） 变换域方法。这种方法是将差分方程变换到z域进行求解，方法简便有效。当然还可以不直接求解差分方程，而是先由差分方程

31、求出系统的单位脉冲响应，再与已知的输入序列进行卷积运算，得到系统输出。但是系统的单位脉冲响应如果不是预先知道，仍然需要求解差分方程，求其零状态响应解。（4） 卷积法：由差分方程求出系统的h(n)，再与已知的x(n) 进行卷积，得到y(n)。观察（）式，求n时刻的输出，要知道n时刻以及n时刻以前的输入序列值，还要知道n时刻以前的N个输出序列值。因此求解差分方程在给定输入序列的条件下，还需要确定N个初始条件。如果求n0时刻以后的输出，n0时刻以前N个输出值y（n0-1）、y（n0-2）、y（n0-N）就构成了初始条件。（）式表明，已知输入序列和N个初始条件，则可以求出n时刻的输出；如果将该公式中的

32、n用n+1代替，可以求出n+1时刻的输出，因此（1.4.1）式表示的差分方程本身就是一个适合递推法求解的方程。第五章 设计过程5.1 用MATLAB求解差分方程 MATLAB信号处理工具箱提供的filter函数实现线性常系数差分方程的递推求解，调用格式如下： yn=filter(B,A.xn) 计算系统对输入信号向量xn的零状态响应输出信号向量yn，yn与xn长度相等，其中，B和A是（）式所给差分方程的系数向量，即 B=b0,b1,bM, A=a0,a1,aN其中a0=1，如果a01，则filter用a0对系数向量B和A归一化。 yn=filter(B,A.xn,xi) 计算系统对输入信号向量

33、xn的全响应输出信号yn。所谓全响应，就是由初始状态引起的零输入响应和由输入信号xn引起的零状态响应之和。其中，xi是等效初始条件的输入序列，所以xi是由初始条件确定的。MATLAB信号处理工具箱提供的filtic就是由初始条件计算xi的函数，其调用格式如下：xi=filtic（B,A,ys,xs）其中，ys和xs是初始条件向量：ys=y(-1),y(-2),y(-3),y(-N),xs=x(-1),x(-2),x(-3),x(-M)。如果xn是因果序列，则xs=0.调用时可缺省xs。 例的MATLAB求解程序ep141.m如下：%1.4.1.m:调用MATLAB解差分方程y(n)-0.8y(

34、n-1)=x(n)a=0.8；ys=1; %设差分方程系数a=0.8，初始状态：y(-1)=1xn=1,zeros(1,30); %x(n)=单位脉冲序列，长度N=31B=1;A=1,-0.8; %差分方程系数xi=filtic(B,A,ys); %由初始条件计算等效初始条件的输入序列xiyn=filter(B,A,xn,xi); %调用fiter解差分方程，求系统输出信号y(n)n=0:length(yn)-1;stem(n,yn,.)title(时域波形图);xlabel(n);ylabel(y(n)程序中取查分方程系数a=0.8时，得到系统输出y（n）如图（a）所示，与例1.4.1的解析

35、递推结果完全相同。如果令初始条件y（-1）=0（仅修改程序中ys=0），则得到系统输出y（n）=h（n），如图1.4.1（b）所示。(a) (b)图（a）为a=0.8，y（-1）=1时，系统输出时域波形图，图（b）为a=0.8，y(-1)=0时，系统输出时域波形图。 第六章 设计代码及结果6.1 MATLAB源程序 源程序如下：%1.m:调用MATLAB解差分方程y(n)-0.8y(n-1)=x(n)ys=1; %初始状态：y(-1)=1xn=1,zeros(1,30); %x(n)=单位脉冲序列，长度N=31B=1;A=1,-0.8; %差分方程系数xi=filtic(B,A,ys); %由

36、初始条件计算等效初始条件的输入序列xiyn=filter(B,A,xn,xi); %调用filter解差分方程，求系统输出信号y(n)n=0:length(yn)-1; %n的取值范围stem(n,yn,.) %画出时域波形图title(时域波形图);xlabel(n);ylabel(y(n) %x轴、y轴分别代表n，x(n) -2 源程序如下：n=-5:5; %n的取值范围xn=0.5.n; %xn=0.5.nstem(n,xn,fill),grid on %画出时域波形图xlabel(n),ylabel(x(n), title(时域波形图) %x轴、y轴分别代表n，x(n)-3 源程序如下

37、：n=-5:5; %n的取值范围a=0.5; %设a=0.5xen=a*a.n+a.(-n); %xen=a*a.n+a.(-n)xon=a*a.n-a.(-n); %xon=a*a.n-a.(-n)figure(1); stem(n,xen,filled),grid on %画出xe(n)时域波形图 title(时域波形图);xlabel(n);ylabel(xe(n) %x轴、y轴分别代表n，xe(n)figure(2); stem(n,xon,filled),grid on %画出xo(n)时域波形图 title(时域波形图);xlabel(n);ylabel(xo(n) %x轴、y轴分

38、别代表n，xo(n)-4 源程序如下：b0=2;b2=0;b2=-1;a1=-0.7;a2=0.1;ys=0; %设差分方程系数，初始状态：y(-1)=1B=2,0,-1;A=1,-0.7,0.1; %差分方程系数n=-5:5; %n的取值范围xn=0.5.n; %x(n)=0.5.nxi=filtic(B,A,ys); %由初始条件计算等效初始条件的输入序列xiyn=filter(B,A,xn,xi); %调用fiter解差分方程，求系统输出信号y(n)stem(n,yn,.) %画出时域波形图title(a);xlabel(n);ylabel(y(n) %x轴、y轴分别代表n，y(n)6.

39、2 程序运行结果 结果如下图 差分序列时域波形图-1 输入参数a1=-0.7,a2=0.1,b0=2,b1=0,b2=-1 得到二阶线性常系数差分方程为 y(n)-0.7y(n-1)+0.1y(n-2)=2x(n)-x(n-2)-2 结果如下：图-2 x（n）的时域波形图-3 结果如下 图-3 共轭对称分量xe（n）的时域波形图 图-3 共轭对称分量xo（n）的时域波形图 -4 输入初始条件ys=y(-1)=1,得到结果如下：图-4 输入初始条件ys=y(-1)=1时，输出y(n)波形图-5（1） 改变初始条件ys=y(-1)=0,得到结果如下： 图-5（1） 输入初始条件ys=y(-1)=0

40、时，输出y(n)波形图-5（2） 改变初始条件ys=y(-1)=100,得到结果如下: 图-5（2） 输入初始条件ys=y(-1)=100时，输出y(n)波形图 6.3 比较结果总结由上实验可知，通过改变初始条件ys的值，得到的输出波形大小并不一致，即输出信号y(n)是不相同的；从而我们可以得出，对于同一个差分方程和同一个输入信号，因为初始条件不同，得到的输出信号是不相同的。第七章 收获与体会 本次MATLAB课程设计让我熟悉了该软件的一些功能，但是对于灵活应用MATLAB，以及掌握各方面的设计思维以及技巧，还需要投入更多的时间。在熟悉MATLAB程序和操作的同时培养了我的独立思考能力，专研精

41、神，解决问题能力和动手能力。在此之前了解到MATLAB是一个很重要很有用的工具，但我并没有完全理解，本课程设计中，通过查阅资料，阅读网上程序并读写程序，对于MATLAB的应用有了更深的了解，同时也认识到MATLAB功能非常的强大，有着很多方面的应用，如绘制函数，处理音频，图像数据，创建用户界面等功能，实为一个功能强大的软件。本次课程设计我完成了基于MATLAB的线性常系数差分方程求解的题目，通过实际操作回顾所学的内容，强化基础，实践理论知识。相信在以后的学习中，还会更加深入的了解MATLAB，应用它。随着课程设计报告的基本完成，本次课程设计终于接近了尾声。本次课程设计要求我们利用上学期所学的信号与线性系统分析的知识结合MATLAB编程工具，完成差分方程求解设计的题目，通过实际操作，回顾所学内容，务实基础，强化理论知识，并体验理论与实际相结合的过程。设计过程中遇到的第一个问题便是对于MATLAB语言的不熟悉，其实现在想想这个问题不应该成为问题。毕竟本专业曾开设过MATLAB程序设计这门课，而且老师还特