数字信号处理课件及答案(西安电子科技版)第1章.ppt

时域离散信号和时域离散系统第 1 章第1章时域离散信号和时域离散系统1.1学习要点与重要公式学习要点与重要公式1.2解线性卷积的方法解线性卷积的方法1.3例题例题1.4习题与上机题解答习题与上机题解答时域离散信号和时域离散系统第 1 章 1.1学习要点与重要公式学习要点与重要公式本章内容是全书的基础。学生从学习模拟信号分析与处理到学习数字信号处理，要建立许多新的概念。数字信号和数字系统与原来的模拟信号和模拟系统不同，尤其是处理方法上有本质的区别。模拟系统用许多模拟器件实现，数字系统则通过运算方法实现。如果读者对本章关于时域离散信号与系统的若干基本概念不清楚，则学到数字滤波器时，会感到“数字信号处理”这门课不好掌握，总觉得学习的不踏实。因此学好本章是极其重要的。时域离散信号和时域离散系统第 1 章1.1.1学习要点学习要点 （1）信号：模拟信号、时域离散信号、数字信号三者之间的区别；常用的时域离散信号；如何判断信号是周期性的，其周期如何计算等。（2）系统：什么是系统的线性、时不变性以及因果性、稳定性；线性、时不变系统输入和输出之间的关系；求解线性卷积的图解法（列表法）、解析法，以及用MATLAB工具箱函数求解；线性常系数差分方程的递推解法。（3）模拟信号的采样与恢复：采样定理；采样前的模拟信号和采样后得到的采样信号之间的频谱关系；如何由采样信号恢复成原来的模拟信号；实际中如何将时域离散信号恢复成模拟信号。时域离散信号和时域离散系统第 1 章1.1.2 重要公式重要公式 （1）这是一个线性卷积公式，注意公式中是在之间对m求和。如果公式中x(n)和h(n)分别是系统的输入和单位脉冲响应，y(n)是系统输出，则该式说明系统的输入、输出和单位脉冲响应之间服从线性卷积关系。时域离散信号和时域离散系统第 1 章（2）x(n)=x(n)*(n)该式说明任何序列与(n)的线性卷积等于原序列。x(nn0)=x(n)*(nn0)（3）这是关于采样定理的重要公式，根据该公式要求对信号的采样频率要大于等于该信号的最高频率的两倍以上，才能得到不失真的采样信号。这是由时域离散信号理想恢复模拟信号的插值公式。时域离散信号和时域离散系统第 1 章1.2解线性卷积的方法解线性卷积的方法解线性卷积是数字信号处理中的重要运算。解线性卷积有三种方法，即图解法（列表法）、解析法和在计算机上用MATLAB语言求解。它们各有特点。图解法（列表法）适合于简单情况，短序列的线性卷积，因此考试中常用，不容易得到封闭解。解析法适合于用公式表示序列的线性卷积，得到的是封闭解，考试中会出现简单情况的解析法求解。解析法求解过程中，关键问题是确定求和限，求和限可以借助于画图确定。第三种方法适合于用计算机求解一些复杂的较难的线性卷积，实验中常用。时域离散信号和时域离散系统第 1 章解线性卷积也可用Z变换法，以及离散傅里叶变换求解，这是后面几章的内容。下面通过例题说明。设x(n)=R4(n),h(n)=R4(n)，求y(n)=x(n)*h(n)。该题是两个短序列的线性卷积，可以用图解法（列表法）或者解析法求解。表1.2.1给出了图解法（列表法），用公式可表示为y(n)=,0,0,1,2,3,4,3,2,1,0,0,时域离散信号和时域离散系统第 1 章时域离散信号和时域离散系统第 1 章下面用解析法求解，写出卷积公式为在该例题中，R4(m)的非零区间为0m3，R4(nm)的非零区间为0nm3，或写成n3mn，这样y(n)的非零区间要求m同时满足下面两个不等式：0m3m3mn上面公式表明m的取值和n的取值有关，需要将n作分段的假设。按照上式，当n变化时，m应该按下式取值：时域离散信号和时域离散系统第 1 章 max0,n3mmin3,n当0n3时，下限应该是0，上限应该是n；当4n6时，下限应该是n3，上限应该是3；当n6时，上面的不等式不成立，因此y(n)=0；这样将n分成三种情况计算：（1）n6时，y(n)=0 （2）0n3时，时域离散信号和时域离散系统第 1 章（3）4n6时，将y(n)写成一个表达式，如下式：y(n)=0n34n6其它时域离散信号和时域离散系统第 1 章在封闭式求解过程中，有时候决定求和的上下限有些麻烦，可借助于非零值区间的示意图确定求和限。在该例题中，非零值区间的示意图如图1.2.1所示。在图1.2.1(b)中，当n0 时,最后得到例例1.3.3设时域离散线性时不变系统的单位脉冲响应h(n)和输入激励信号x(n)分别为时域离散信号和时域离散系统第 1 章x(n)=cos(n)u(n)求系统的稳态响应y(n)。解解 x(n)=cos(n)u(n)=(1)nu(n)时域离散信号和时域离散系统第 1 章当n时，稳态解为 例例1.3.4假设5项滑动平均滤波器的差分方程为y(n)=x(n)+x(n1)+x(n2)+x(n3)+x(n4)输入信号用图1.3.1表示，画出该滤波器输出的前16个序列值的波形，并说明该滤波器对输入信号起什么作用。时域离散信号和时域离散系统第 1 章图1.3.1时域离散信号和时域离散系统第 1 章解解：已知系统的差分方程和输入信号求系统输出，可以用递推法求解，这里采用MATLAB函数filter 计算。调用MATLAB函数filter计算该系统的系统响应的程exp134.m如下：%程序exp134.m%调用conv实现5项滑动平均滤波xn=0.5*ones(1，15)；xn(4)=1；xn(6)=1；xn(10)=1；hn=ones(1，5)；yn=conv(hn，xn)；时域离散信号和时域离散系统第 1 章%以下为绘图部分n=0：length(yn)1；subplot(2，1，1)；stem(n，yn，.)xlabel(n)；ylabel(y(n)程序运行结果如图1.3.2所示。由图形可以看出，5项滑动平均滤波器对输入波形起平滑滤波作用，将信号的第4、8、12、16的序列值平滑去掉。时域离散信号和时域离散系统第 1 章图1.3.2时域离散信号和时域离散系统第 1 章例例1.3.5已知x1(n)=(n)+3(n1)+2(n2)，x2(n)=u(n)u(n3)，试求信号x(n)，它满足x(n)=x1(n)*x2(n),并画出x(n)的波形。解解：这是一个简单的计算线性卷积的题目。x(n)=x1(n)*x2(n)=(n)+3(n1)+2(n2)*u(n)u(n3)=(n)+3(n1)+2(n2)*R3(n)=R3(n)+3R3(n1)+2R3(n2)=(n)+4(n1)+6(n2)+5(n3)+2(n4)画出x(n)的波形如图1.3.3所示。时域离散信号和时域离散系统第 1 章图1.3.3时域离散信号和时域离散系统第 1 章例例1.3.6已知离散信号x(n)如图1.3.4(a)所示，试求y(n)=x(2n)*x(n)，并绘出y(n)的波形。（选自西安交通大学2001年攻读硕士学位研究生入学考试试题）解解：这也是一个计算线性卷积的题目，只不过要先求出x(2n)。解该题适合用列表法(图解法)。x(2n)=1,1,1,0.5y(n)=x(2n)*x(n)=1,2,3,3,3,3,2.75,2,1,0.25绘出y(n)的波形如图1.3.4(b)所示。时域离散信号和时域离散系统第 1 章图1.3.4时域离散信号和时域离散系统第 1 章1.4习题与上机题解答习题与上机题解答1.用单位脉冲序列(n)及其加权和表示题1图所示的序列。题1图时域离散信号和时域离散系统第 1 章解：x(n)=(n+4)+2(n+2)(n+1)+2(n)+(n1)+2(n2)+4(n3)+0.5(n4)+2(n6)2 给定信号：2n+54n160n40 其它(1)画出x(n)序列的波形，标上各序列值；(2)试用延迟的单位脉冲序列及其加权和表示x(n)序列；(x(n)=时域离散信号和时域离散系统第 1 章（3）令x1(n)=2x(n2)，试画出x1(n)波形；（4）令x2(n)=2x(n+2)，试画出x2(n)波形；（5）令x3(n)=x(2n)，试画出x3(n)波形。解解：(1)x(n)序列的波形如题2解图（一）所示。(2)x(n)=3(n+4)(n+3)+(n+2)+3(n+1)+6(n)+6(n1)+6(n2)+6(n3)+6(n4)时域离散信号和时域离散系统第 1 章（3）x1(n)的波形是x(n)的波形右移2位，再乘以2，画出图形如题2解图（二）所示。(4)x2(n)的波形是x(n)的波形左移2位，再乘以2，画出图形如题2解图（三）所示。(5)画x3(n)时，先画x(n)的波形(即将x(n)的波形以纵轴为中心翻转180)，然后再右移2位，x3(n)波形如题2解图（四）所示。时域离散信号和时域离散系统第 1 章题2解图（一）时域离散信号和时域离散系统第 1 章题2解图（二）时域离散信号和时域离散系统第 1 章题2解图（三）时域离散信号和时域离散系统第 1 章题2解图（四）时域离散信号和时域离散系统第 1 章3 判断下面的序列是否是周期的;若是周期的，确定其周期。(1)(2)解解：（1）因为=,所以,这是有理数，因此是周期序列，周期T=14。（2）因为=,所以=16,这是无理数，因此是非周期序列。时域离散信号和时域离散系统第 1 章4 对题1图给出的x(n)要求：(1)画出x(n)的波形；(2)计算xe(n)=x(n)+x(n)，并画出xe(n)波形；(3)计算xo(n)=x(n)x(n)，并画出xo(n)波形;(4)令x1(n)=xe(n)+xo(n),将x1(n)与x(n)进行比较，你能得到什么结论？时域离散信号和时域离散系统第 1 章解解：（1）x(n)的波形如题4解图（一）所示。(2)将x(n)与x(n)的波形对应相加，再除以2，得到xe(n)。毫无疑问，这是一个偶对称序列。xe(n)的波形如题4解图（二）所示。(3)画出xo(n)的波形如题4解图（三）所示。时域离散信号和时域离散系统第 1 章题4解图（一）时域离散信号和时域离散系统第 1 章题4解图（二）时域离散信号和时域离散系统第 1 章题4解图（三）时域离散信号和时域离散系统第 1 章(4)很容易证明：x(n)=x1(n)=xe(n)+xo(n)上面等式说明实序列可以分解成偶对称序列和奇对称序列。偶对称序列可以用题中(2)的公式计算，奇对称序列可以用题中(3)的公式计算。5 设系统分别用下面的差分方程描述，x(n)与y(n)分别表示系统输入和输出，判断系统是否是线性非时变的。（1）y(n)=x(n)+2x(n1)+3x(n2)（2）y(n)=2x(n)+3 （3）y(n)=x(nn0)n0为整常数 （4）y(n)=x(n)时域离散信号和时域离散系统第 1 章（5）y(n)=x2(n)（6）y(n)=x(n2)（7）y(n)=（8）y(n)=x(n)sin(n)解解：（1）令输入为x(nn0)输出为 y(n)=x(nn0)+2x(nn01)+3x(nn02)y(nn0)=x(nn0)+2x(nn01)+3(nn02)=y(n)时域离散信号和时域离散系统第 1 章故该系统是非时变系统。因为 y(n)=Tax1(n)+bx2(n)=ax1(n)+bx2(n)+2ax1(n1)+bx2(n1)+3ax1(n2)+bx2(n2)Tax1(n)=ax1(n)+2ax1(n1)+3ax1(n2)Tbx2(n)=bx2(n)+2bx2(n1)+3bx2(n2)所以 Tax1(n)+bx2(n)=aTx1(n)+bTx2(n)故该系统是线性系统。时域离散信号和时域离散系统第 1 章（2）令输入为x(nn0)输出为y(n)=2x(nn0)+3y(nn0)=2x(nn0)+3=y(n)故该系统是非时变的。由于Tax1(n)+bx2(n)=2ax1(n)+2bx2(n)+3Tax1(n)=2ax1(n)+3Tbx2(n)=2bx2(n)+3Tax1(n)+bx2(n)aTx1(n)+bTx2(n)故该系统是非线性系统。时域离散信号和时域离散系统第 1 章(3)这是一个延时器，延时器是线性非时变系统，下面证明。令输入为x(nn1)输出为y(n)=x(nn1n0)y(nn1)=x(nn1n0)=y(n)故延时器是非时变系统。由于Tax1(n)+bx2(n)=ax1(nn0)+bx2(nn0)=aTx1(n)+bTx2(n)故延时器是线性系统。时域离散信号和时域离散系统第 1 章(4)y(n)=x(n)令输入为x(nn0)输出为y(n)=x(n+n0)y(nn0)=x(n+n0)=y(n)因此系统是线性系统。由于Tax1(n)+bx2(n)=ax1(n)+bx2(n)=aTx1(n)+bTx2(n)因此系统是非时变系统。时域离散信号和时域离散系统第 1 章（5）y(n)=x2(n)令输入为 x(nn0)输出为y(n)=x2(nn0)y(nn0)=x2(nn0)=y(n)故系统是非时变系统。由于 Tax1(n)+bx2(n)=ax1(n)+bx2(n)2 aTx1(n)+bTx2(n)=ax21(n)+bx22(n)因此系统是非线性系统。时域离散信号和时域离散系统第 1 章（6）y(n)=x(n2)令输入为x(nn0)输出为y(n)=x(nn0)2)y(nn0)=x(nn0)2)=y(n)故系统是非时变系统。由于Tax1(n)+bx2(n)=ax1(n2)+bx2(n2)=aTx1(n)+bTx2(n)故系统是线性系统。时域离散信号和时域离散系统第 1 章（7）y(n)=x(m)令输入为x(nn0)输出为 y(n)=0DD)x(m-n0)y(nn0)=x(m)y(n)故系统是时变系统。由于Tax1(n)+bx2(n)=ax1(m)+bx2(m)=aTx1(n)+bTx2(n)故系统是线性系统。时域离散信号和时域离散系统第 1 章（8）y(n)=x(n)sin(n)令输入为x(nn0)输出为y(n)=x(nn0)sin(n)y(nn0)=x(nn0)sin(nn0)y(n)故系统不是非时变系统。由于Tax1(n)+bx2(n)=ax1(n)sin(n)+bx2(n)sin(n)=aTx1(n)+bTx2(n)故系统是线性系统。时域离散信号和时域离散系统第 1 章6 给定下述系统的差分方程，试判定系统是否是因果稳定系统，并说明理由。（1）y(n)=x(nk)（2）y(n)=x(n)+x(n+1)（3）y(n)=x(k)（4）y(n)=x(nn0)（5）y(n)=ex(n)时域离散信号和时域离散系统第 1 章解解：（1）只要N1，该系统就是因果系统，因为输出只与n时刻的和n时刻以前的输入有关。如果|x(n)|M，则|y(n)|M，因此系统是稳定系统。（2）该系统是非因果系统，因为n时间的输出还和n时间以后（n+1）时间）的输入有关。如果|x(n)|M,则|y(n)|x(n)|+|x(n+1)|2M,因此系统是稳定系统。（3）如果|x(n)|M，则|y(n)|x(k)|2n0+1|M，因此系统是稳定的；假设n00，系统是非因果的，因为输出还和x(n)的将来值有关。时域离散信号和时域离散系统第 1 章（4）假设n00，系统是因果系统，因为n时刻输出只和n时刻以后的输入有关。如果|x(n)|M,则|y(n)|M，因此系统是稳定的。（5）系统是因果系统，因为系统的输出不取决于x(n)的未来值。如果|x(n)|M,则|y(n)|=|ex(n)|e|x(n)|eM，因此系统是稳定的。7 设线性时不变系统的单位脉冲响应h(n)和输入序列x(n)如题7图所示，要求画出y(n)输出的波形。解解：解法（一）采用列表法。y(n)=x(n)*h(n)=x(m)h(nm)时域离散信号和时域离散系统第 1 章题7图时域离散信号和时域离散系统第 1 章y(n)=2,1,0.5,2,1,4.5,2,1;n=2,1,0,1,2,3,4,5时域离散信号和时域离散系统第 1 章解法（二）采用解析法。按照题7图写出x(n)和h(n)的表达式分别为x(n)=(n+2)+(n1)+2(n3)h(n)=2(n)+(n1)+(n2)由于x(n)*(n)=x(n)x(n)*A(nk)=Ax(nk)故时域离散信号和时域离散系统第 1 章y(n)=x(n)*h(n)=x(n)*2(n)+(n1)+(n2)=2x(n)+x(n1)+x(n2)将x(n)的表示式代入上式，得到 y(n)=2(n+2)(n+1)0.5(n)+2(n1)+(n2)+4.5(n3)+2(n4)+(n5)时域离散信号和时域离散系统第 1 章8.设线性时不变系统的单位脉冲响应h(n)和输入x(n)分别有以下三种情况，分别求出输出y(n)。（1）h(n)=R4(n),x(n)=R5(n)（2）h(n)=2R4(n),x(n)=(n)(n2)（3）h(n)=0.5nu(n),xn=R5(n)解解：（1）y(n)=x(n)*h(n)=R4(m)R5(nm)先确定求和域。由R4(m)和R5(nm)确定y（n）对于m的非零区间如下:0m34mn时域离散信号和时域离散系统第 1 章根据非零区间，将n分成四种情况求解：n7时，y(n)=0时域离散信号和时域离散系统第 1 章最后结果为 0 n7 n+1 0n3 8n4n7y(n)的波形如题8解图（一）所示。(2)y(n)=2R4(n)*(n)(n2)=2R4(n)2R4(n2)=2(n)+(n1)(n+4)(n+5)y(n)的波形如题8解图（二）所示y(n)=时域离散信号和时域离散系统第 1 章题8解图（一）时域离散信号和时域离散系统第 1 章题8解图（二）时域离散信号和时域离散系统第 1 章(3)y(n)=x(n)*h(n)=R5(m)0.5nmu(nm)=0.5nR5(m)0.5mu(nm)y（n）对于m 的非零区间为 0m4,mn n0时，y(n)=0 0n4时，时域离散信号和时域离散系统第 1 章=(10.5n1)0.5n=20.5n n5时最后写成统一表达式：y(n)=(20.5n)R5(n)+310.5nu(n5)时域离散信号和时域离散系统第 1 章9 证明线性卷积服从交换律、结合律和分配律，即证明下面等式成立：（1）x(n)*h(n)=h(n)*x(n)（2）x(n)*(h1(n)*h2(n)=(x(n)*h1(n)*h2(n)（3）x(n)*(h1(n)+h2(n)=x(n)*h1(n)+x(n)*h2(n)证明：（1）因为令m=nm,则时域离散信号和时域离散系统第 1 章(2)利用上面已证明的结果，得到时域离散信号和时域离散系统第 1 章交换求和号的次序，得到时域离散信号和时域离散系统第 1 章10 设系统的单位脉冲响应h(n)=(3/8)0.5nu(n)，系统的输入x(n)是一些观测数据，设x(n)=x0,x1,x2,xk,，试利用递推法求系统的输出y(n)。递推时设系统初始状态为零状态。时域离散信号和时域离散系统第 1 章解解：n=0时，n0n=1时，时域离散信号和时域离散系统第 1 章n=2时，最后得到11 设系统由下面差分方程描述：设系统是因果的，利用递推法求系统的单位脉冲响应。时域离散信号和时域离散系统第 1 章解解：令x(n)=(n),则n=0时，n=1时，时域离散信号和时域离散系统第 1 章n=2时，n=3时，归纳起来，结果为时域离散信号和时域离散系统第 1 章12.设系统用一阶差分方程y(n)=ay(n1)+x(n)描述，初始条件y(-1)=0，试分析该系统是否是线性非时变系统。解解：分析的方法是让系统输入分别为(n)、(n1)、(n)+(n1)时，求它的输出，再检查是否满足线性叠加原理和非时变性。（1）令x(n)=(n),这时系统的输出用y1(n)表示。该情况在教材例1.4.1 中已求出，系统的输出为y1(n)=anu(n)时域离散信号和时域离散系统第 1 章(2)令x(n)=(n1)，这时系统的输出用y2(n)表示。n=0时，n=1时，n=2时，任意 n 时，时域离散信号和时域离散系统第 1 章最后得到(3)令x(n)=(n)+(n1),系统的输出用y3(n)表示。n=0时，n=1时，n=2时，时域离散信号和时域离散系统第 1 章n=3时，任意 n 时，最后得到时域离散信号和时域离散系统第 1 章由（1）和（2）得到y1(n)=T(n),y2(n)=T(n1)y1(n)=y2(n1)因此可断言这是一个时不变系统。情况（3）的输入信号是情况（1）和情况（2）输入信号的相加信号，因此y3(n)=T(n)+(n1)。观察y1(n)、y2(n)、y3(n)，得到y3(n)=y1(n)+y2(n)，因此该系统是线性系统。最后得到结论：用差分方程y(n)=ay(n1)+x(n),0a1描写的系统，当初始条件为零时，是一个线性时不变系统。时域离散信号和时域离散系统第 1 章13 有一连续信号xa(t)=cos(2ft+j)，式中，f=20 Hz,j=/2。（1）求出xa(t)的周期；（2）用采样间隔T=0.02 s对xa(t)进行采样，试写出采样信号 的表达式；（3）画出对应 的时域离散信号（序列）x(n)的波形，并求出x(n)的周期。解解：（1）xa(t)的周期为时域离散信号和时域离散系统第 1 章（2）（3）x(n)的数字频率=0.8，故，因而周期N=5，所以 x(n)=cos(0.8n+/2)画出其波形如题13解图所示。时域离散信号和时域离散系统第 1 章题13解图时域离散信号和时域离散系统第 1 章14.已知滑动平均滤波器的差分方程为（1）求出该滤波器的单位脉冲响应；（2）如果输入信号波形如前面例1.3.4的图1.3.1所示，试求出y(n)并画出它的波形。解：（1）将题中差分方程中的x(n)用(n)代替，得到该滤波器的单位脉冲响应，即时域离散信号和时域离散系统第 1 章（2）已知输入信号，用卷积法求输出。输出信号y(n)为表1.4.1表示了用列表法解卷积的过程。计算时，表中x(k)不动，h(k)反转后变成h(k)，h(nk)则随着n的加大向右滑动，每滑动一次，将h(nk)和x(k)对应相乘，再相加和平均，得到相应的y(n)。“滑动平均”清楚地表明了这种计算过程。最后得到的输出波形如前面图1.3.2所示。该图清楚地说明滑动平均滤波器可以消除信号中的快速变化，使波形变化缓慢。时域离散信号和时域离散系统第 1 章时域离散信号和时域离散系统第 1 章15*.已知系统的差分方程和输入信号分别为用递推法计算系统的零状态响应。解：求解程序ex115.m如下：%程序ex115.m%调用filter解差分方程y(n)+0.5y(n1)=x(n)+2x(n2)xn=1，2，3，4，2，1，zeros(1，10)；%x(n)=单位脉冲序列，长度N=31B=1，0，2；A=1，0.5；%差分方程系数时域离散信号和时域离散系统第 1 章yn=filter(B，A，xn)%调用filter解差分方程，求系统输出信号y(n)n=0：length(yn)1；subplot(3，2，1)；stem(n，yn，.)；axis(1，15，2，8)title(系统的零状态响应)；xlabel(n)；ylabel(y(n)程序运行结果：时域离散信号和时域离散系统第 1 章yn=1.0000 1.5000 4.2500 5.8750 5.0625 6.4688 0.7656 1.6172 -0.8086 0.4043-0.2021 0.1011 -0.0505 0.0253 -0.0126 0.0063 -0.0032 0.0016 -0.0008 0.0004 -0.0002 0.0001 -0.0000 0.0000 -0.0000 0.0000程序运行结果的y(n)波形图如题15*解图所示。时域离散信号和时域离散系统第 1 章题15*解图时域离散信号和时域离散系统第 1 章16*.已知两个系统的差分方程分别为 （1）y(n)=0.6y(n1)0.08y(n2)+x(n)（2）y(n)=0.7y(n1)0.1y(n2)+2x(n)x(n2)分别求出所描述的系统的单位脉冲响应和单位阶跃响应。解解：（1）系统差分方程的系数向量为B1=1，A1=1，0.6，0.08（2）系统差分方程的系数向量为B2=2，0，1,A2=1，0.7，0.1时域离散信号和时域离散系统第 1 章调用MATLAB函数filter计算两个系统的系统的单位脉冲响应和单位阶跃响应的程序ex116.m如下：%程序ex116.mB1=1；A1=1，0.6，0.08；%设差分方程（1）系数向量B2=2，0，1；A2=1，0.7，0.1；%设差分方程（2）系数向量%=%系统1xn=1，zeros(1，30)；%xn=单位脉冲序列，长度N=31xi=filtic(B1，A1，ys)；%由初始条件计算等效初始条件输入序列xi时域离散信号和时域离散系统第 1 章hn1=filter(B1，A1，xn，xi)；%调用filter解差分方程，求系统输出信号hn1n=0：length(hn1)-1；subplot(3，2，1)；stem(n，hn1，.)title(a)系统1的系统单位脉冲响应)；xlabel(n)；ylabel(h(n)xn=ones(1，30)；%xn=单位阶跃序列，长度N=31sn1=filter(B1，A1，xn，xi)；%调用filter解差分方程，求系统输出信号sn1n=0：length(sn1)1；subplot(3，2，2)；stem(n，sn1，.)时域离散信号和时域离散系统第 1 章itle(b)系统1的单位阶跃响应)；xlabel(n)；ylabel(s(n)%=%系统2xn=1，zeros(1，30)；%xn=单位脉冲序列，长度N=31xi=filtic(B2，A2，ys)；%由初始条件计算等效初始条件输入序列xihn2=filter(B2，A2，xn，xi)；%调用filter解差分方程，求系统输出信号hn2n=0：length(hn2)1；subplot(3，2，5)；stem(n，hn2，.)时域离散信号和时域离散系统第 1 章title(a)系统2的系统单位脉冲响应)；xlabel(n)；ylabel(h(n)xn=ones(1，30)；%xn=单位阶跃序列，长度N=31sn2=filter(B2，A2，xn，xi)；%调用filter解差分方程，求系统输出信号sn2n=0：length(sn2)1；subplot(3，2，6)；stem(n，sn2，.)title(b)系统2的单位阶跃响应)；xlabel(n)；ylabel(s(n)程序运行结果如题16*解图所示。时域离散信号和时域离散系统第 1 章题16*解图时域离散信号和时域离散系统第 1 章17*.已知系统的差分方程为y(n)=a1y(n1)a2y(n2)+bx(n)其中,a1=0.8,a2=0.64,b=0.866。（1）编写求解系统单位脉冲响应h(n)(0n49)的程序，并画出h(n)(0n49)；（2）编写求解系统零状态单位阶跃响应s(n)(0n100）的程序，并画出s(n)(0n100）。时域离散信号和时域离散系统第 1 章解解：调用MATLAB函数filter计算该系统的系统响应的程序ex117.m如下：%程序ex117.m%调用filter解差分方程，求系统单位脉冲响应和单位阶跃响应B=0.866；A=1，0.8，0.64；%差分方程系数向量%=%（1）求解系统单位脉冲响应，并画出h(n)xn=1，zeros(1，48)；%xn=单位脉冲序列，长度N=31时域离散信号和时域离散系统第 1 章hn=filter(B1，A1，xn)；%调用filter解差分方程，求系统输出信号hnn=0：length(hn)1；subplot(3，2，1)；stem(n，hn，.)title(a)系统的单位脉冲响应)；xlabel(n)；ylabel(h(n)%=%（2）求解系统单位阶跃响应，并画出h(n)xn=ones(1，100)；%xn=单位阶跃序列，长度N=100时域离散信号和时域离散系统第 1 章sn=filter(B，A，xn)；%调用filter解差分方程，求系统单位阶跃响应snn=0：length(sn)1；subplot(3，2，2)；stem(n，sn，.)；axis(0，30，0，2)title(b)系统的单位阶跃响应)；xlabel(n)；ylabel(s(n)%=程序运行结果如题17*解图所示。时域离散信号和时域离散系统第 1 章题17*解图时域离散信号和时域离散系统第 1 章18*.在题18*图中，有四个分系统T1、T2、T3和T4，四个分系统分别用下面的单位脉冲响应或者差分方程描述：时域离散信号和时域离散系统第 1 章题18*图编写程序计算整个系统的单位脉冲响应h(n),0n99。时域离散信号和时域离散系统第 1 章解解：由题18*图可知，可以采用以下步骤计算整个系统的单位脉冲响应h(n)。设x(n)=(n)，则v(n)=h1(n)*h2(n)+h3(n)该式调用conv函数计算。h(n)=T4v(n)该式调用filter函数计算。调用MATLAB函数conv和filter计算该系统的系统响应的程序ex118.m如下：时域离散信号和时域离散系统第 1 章%程序ex118.m%调用conv和filter求总系统单位脉冲响应序列h1n=1，1/2，1/4，1/8，1/16，1/32；%对h1n赋值 h2n=ones(1，6)；h3n=1/4，1/2，1/4，zeros(1，97)；%计算v(n)=h1(n)*h2(n)+h3(n)h12n=conv(h1n，h2n)；h12n=h12n，zeros(1，89)；vn=h12n+h3n；时域离散信号和时域离散系统第 1 章%调用filer计算hn等于T4对vn响应B4=1，1；A4=1，-0.9，0.81；hn=filter(B4，A4，vn)；%以下为绘图部分n=0：length(hn)-1；subplot(2，1，1)；stem(n，hn，.)xlabel(n)；ylabel(h(n)程序运行结果如题18*解图所示。时域离散信号和时域离散系统第 1 章题18*解图