第5章离散时间信号与系统的时域分析 前4章讨论了连续时间信号与系统的分析方法,这些方法在通信和自动控制等领域有着广泛的应用。随着计算机和微处理器的出现,数字信号处理技术得到了迅速发展。数字信号处理具有精度高、稳定性好、灵活性强和集成度高等特点,逐渐被图像处理、工业控制、军事以及人工智能等各领域所采用。数字信号是一种离散时间信号,其相关的处理技术离不开离散时间信号和系统的分析方法。 本章主要讨论离散时间信号与系统的基本概念和分析方法。首先介绍离散时间信号描述方法和运算规则,然后分析离散时间系统的特性与建模方法,最后讨论系统响应的时域求解方法。 5.1离散时间信号 离散时间信号与连续时间信号联系紧密,在研究内容和分析方法上具有一定的相似性。本节主要介绍离散时间信号的概念、常用的离散时间信号以及运算方法。 5.1.1离散时间信号及表示方法 若信号只在一系列离散的时间点上有定义,则该信号称为离散时间信号。在实际应用中,通常信号取值的时刻为某个时间间隔t0的整数倍,所以离散时间信号可以用f(nt0)来表示,通常简写为f(n)。其中自变量n体现了函数值的出现次序,称为序号,故离散时间信号又称为序列。 离散时间信号可以由连续时间信号采样得到,时域采样的具体过程在3.7节进行过讨论。实际生活中许多物理现象都可以利用采样信号来描述。例如在气象观测中,要绘制出一天内的气温变化,通常选取若干特定时刻对气温进行测量,因而此时描述气温变化的信号为离散时间信号,如图51所示。 图51气温变化图 离散时间信号也可以从应用中直接获得。例如在生物学、遗传学、概率、统计及社会学等领域,存在大量离散的数据,其自变量本身就具有离散性。此外,自变量也不一定始终为时间,例如温度随高度的变化等。但是为了表示方便,一般统一用n来表示自变量。 离散时间信号常用以下三种方式描述: (1) 数学表达式。当f(n)遵循某种变化规律,则可用确定性的函数表达式来描述。如 图52序列的图形 f(n)=12(n+1),n≥0(5.11) (2) 图形。用图形描述可直观看出离散时间信号的变化情况,例如式(5.11)中的信号f(n)的图形如图52所示。 (3) 数组。用数组表示离散时间信号,是将其信号值一一罗列出来,并标明位置。式(5.11)所描述的信号f(n)用数组可表示为f(n)=0.5↑,1,1.5,2,…,式中箭头表示n=0的位置。 5.1.2常用序列 1. 单位样值序列 图53单位样值序列 单位样值序列又称为单位脉冲序列,用δ(n)来表示,其定义为 δ(n)=1,n=0 0,n≠0(5.12) 需要注意,δ(n)在n=0时的函数值等于1,而不是无穷量,波形如图53所示。 2. 单位阶跃序列 单位阶跃序列u(n)定义为 u(n)=1,n≥0 0,n<0(5.13) 单位阶跃序列u(n)与单位阶跃信号u(t)类似,但两者的区别除了自变量取值是连续还 是离散外,还需要注意u(n)在0时刻给出了确定的函数值1。u(n)的波形如图54所示。 u(n)可以用来描述信号的存在范围,例如f(n)=12(n+1),n≥0又可以表示为 f(n)=12(n+1)u(n) u(n)和δ(n)存在如下关系: u(n)=∑+∞m=0δ(n-m)(5.14) δ(n)=u(n)-u(n-1)(5.15) 3. 单位矩形序列 单位矩形序列用RN(n)来表示,其定义为 RN(n)=1,0≤n≤N-1 0,其他(5.16) RN(n)中的“N”表示从n=0到n=N-1共有N个数值。单位矩形序列RN(n)的波形如图55所示,该序列又可通过单位样值序列和单位阶跃序列来描述,即 RN(n)=∑N-1m=0δ(n-m)=u(n)-u(n-N)(5.17) 图54单位阶跃序列 图55单位矩形序列 4. 单边实指数序列 f(n)=anu(n)(5.18) 指数序列的变化规律由底数a取决。当a>1时,序列是发散的; 00,则序列右移m个单位可表示为 y(n)=f(n-m)(5.117) 类似地,序列左移m个单位可表示为 y(n)=f(n+m)(5.118) 例53已知序列f(n)={2,1↑,3},求序列y(n)=f(n)+f(n-1)f(n-2)。 解: 已知f(n)={2,1↑,3},根据序列的移位运算,得 f(n-1)={2↑,1,3} f(n-2)={0↑,2,1,3} 将以上两个序列进行相乘运算,得 f(n-1)f(n-2)={0↑,2,3} 最后将上式与原序列f(n)相加后,得 y(n)={2,1↑,5,3} 另外还可以依次画出以上各序列的波形,根据波形进行运算,具体过程如图510(a)~(e)所示。 图510例53中各序列的波形 根据图510(e)的结果,可知 y(n)={2,1↑,5,3} 5. 尺度变换 序列的尺度变换是指将f(n)波形压缩或扩展为一个新的序列。若将自变量n乘以大于1的正整数a,构成序列f(an),则该序列是原序列f(n)的波形压缩,按每隔a-1个点取值并将取出的值重新排列。若将自变量n除以大于1的正整数a,构成序列f(n/a),则该序列是原序列f(n)的波形扩展,即在原序列相邻两点间插入a-1个零点值并重新排列。 例54若f(n)波形如图511(a)所示,画出f(2n)和f(n/2)的波形。 解: f(2n)是对原序列f(n)的波形压缩,具体方法是将原序列每隔1个点抽取1个样值,也就是抽取了原序列偶数点的值,然后再重新排序得到,如图511(b)所示。 f(n/2)是对原序列f(n)的波形扩展,即将原序列每相邻2个点之间插入1个零点值,后对每个序列值进行自左向右重新排序,使得原序列的点位于偶数点上,奇数点为插入的零点值,如图511(c)所示。 图511例54中各序列的波形 6. 差分 序列的差分运算是指原序列f(n)与其移位序列f(n-m)的相减运算,具体有前向差分和后向差分两种形式。 前向差分记为 Δx(n)=x(n+1)-x(n)(5.119) 后向差分记为 x(n)=x(n)-x(n-1)(5.120) 5.2离散时间系统 图512离散时间系统的框图 离散时间系统是指输入、输出信号均为离散时间信号的系统。离散时间系统的框图如图512所示。经过系统的作用后,输入f(n)转变为输出y(n)。本节针对单输入单输出系统,介绍其建模方法和特性分析。 5.2.1离散时间系统模型 系统模型是系统物理特性的数学抽象。要进行系统分析,首先要建立系统模型。离散时间系统模型主要包括数学模型和系统模拟图。 1. 数学模型——差分方程 下面通过两个例题来说明如何建立系统的数学模型,以及离散时间系统数学模型的具体形式。 例55某储户每月月初定期在银行存款。设第n个月的存款是e(n),银行支付月息为a,每月利息不取出。设第n个月初的本息和为r(n),试列写r(n)满足的方程。 解: 第n个月的本息和r(n)包括三个部分: 第n-1个月初的本息和r(n-1)、r(n-1)在第n-1个月所产生的利息ar(n-1)及第n个月初的存款e(n),因而可得 r(n)=r(n-1)+ar(n-1)+e(n) 整理得 r(n)-(1+a)r(n-1)=e(n)(5.21) 例56图513所示为梯形电阻网络。令每个结点相对地的电压为v(n),n=1,2,…,N。 试建立当n0时,由于δ(n)的函数值为0,因而将单位样值响应的数学模型简化为齐次差分方程,求解过程与零输入响应的求解类似,关键在于确定单位样值响应的初始条件h(0),h(1),…,h(N-1)。 例516已知描述某系统的差分方程r(n)-5r(n-1)=e(n),求该系统的单位样值响应。 解: 求解单位样值响应时,系统的数学模型为 h(n)-5h(n-1)=δ(n)(5.41) 激励δ(n)只在n=0时取值为1,因而当n>0时,数学模型转化为 h(n)-5h(n-1)=0 利用齐次差分方程的求解方法,可得 h(n)=c·5n 由单位样值响应的定义可知h(-1)=0,将其代入式(5.41),可得 h(0)-5h(-1)=δ(0) 解得h(0)=1。当n>0时,激励δ(n)=0,但由于δ(n)在0时刻的作用转化成初始状态h(0)=1,故使系统在n>0的区间产生响应。 将h(0)=1代入h(n)=c·5n,可得 c=1 系统的单位样值响应为 h(n)=5n,n≥0 或者写为 h(n)=5nu(n) 例517已知某系统的差分方程r(n)-5r(n-1)+6r(n-2)=e(n)+e(n-1),求该系统的单位样值响应。 解: 求解单位样值响应时,系统的数学模型为 h(n)-5h(n-1)+6h(n-2)=δ(n)+δ(n-1)(5.42) 激励仅作用于n=0及n=1两个时刻点,当n>1时,数学模型转化为 h(n)-5h(n-1)+6h(n-2)=0 由齐次差分方程的求解方法,可得 h(n)=c1·2n+c2·3n(5.43) 利用h(-2)=h(-1)=0代入式(5.42),可得 h(0)-5h(-1)+6h(-2)=δ(0)+δ(-1) h(1)-5h(0)+6h(-1)=δ(1)+δ(0) 解得h(0)=1,h(1)=6,将其代入式(5.43),得 c1+c2=1 2c1+3c2=6,c1=-3 c2=4 故单位样值响应为 h(n)=(-3·2n+4·3n)u(n) 另解: 假设激励δ(n)和δ(n-1)分别单独作用于系统,先求δ(n)作用时的单位样值响应h1(n)。 求解齐次解为 h1(n)=c1·2n+c2·3n 利用迭代法,有 h(0)-5h(-1)+6h(-2)=δ(0) h(1)-5h(0)+6h(-1)=δ(1) 解得h(0)=1,h(1)=5,将其代入齐次解的形式中,得 c1+c2=1 2c1+3c2=5,c1=-2 c2=3 由此可得 h1(n)=(3·3n-2·2n)u(n) 当δ(n-1)单独作用于系统时,由系统时不变特性可知 δ(n-1)→h1(n-1)=(3·3n-1-2·2n-1)u(n-1) =(3n-2n)u(n-1) 将以上结果叠加,利用系统叠加性,可得 δ(n)+δ(n-1)→h1(n)+h1(n-1)=(3·3n-2·2n)u(n)+(3n-2n)u(n-1) =(4·3n-3·2n)u(n-1)+δ(n) =(4·3n-3·2n)u(n) 2. 传输算子法 利用初始条件等效法求解高阶差分方程时需要求解齐次方程,并等效出系统N个初始条件,即h(0),h(1),…,h(N-1),有时比较烦琐。为了简化运算,可以使用移位算子求解单位样值响应。 利用移位算子描述差分方程,则式(5.31)的算子方程为 ∑Nk=0akE-kr(n)=∑Mm=0bmE-me(n)(5.44) 当激励为δ(n)时,系统模型为 ∑Nk=0akE-kh(n)=∑Mm=0bmE-mδ(n) 故单位样值响应为 h(n)=b0+b1E-1+…+bME-Ma0+a1E-1+…+aNE-Nδ(n)=H(E)δ(n)(5.45) 式(5.45)中,H(E)=b0+b1E-1+…+bME-Ma0+a1E-1+…+aNE-N称为传输算子。 假设系统特征根无重根,对式(5.45)进行部分分式展开,可得 h(n)=csE-s+cs-1E-s+1+…+c1E-1+c0+k11-λ1E-1+k21-λ2E-1+…+kN1-λNE-1δ(n)(5.46) 式(5.46)是M≥N且系统特征根为单根。若Mn0时刻的输入无关。满足这一个特性的系统称为因果系统,不满足则为非因果系统。 离散时间系统满足因果性的充分必要条件是 h(n)=0,n<0(5.49) 例如,r(n)=e(n+1)-e(n),求解该系统的单位样值响应,可得h(n)=δ(n+1)-δ(n)。当n=-1时,h(-1)=1,故可以判断出该系统为非因果系统。 2) 稳定性判断 系统的稳定性定义为有界的输入产生有界的输出。满足稳定性的系统称为稳定系统,不满足则为不稳定系统。 对于离散时间系统,稳定性的充分必要条件是系统的单位样值响应满足绝对可和,即 ∑+∞n=-∞h(n)≤M,M为有界正值(5.410) 如果某离散时间系统既满足因果性,也满足稳定性,则该系统称为因果稳定系统。 例519已知某系统的单位样值响应为h(n)=anu(n),试判断该系统的因果性和稳定性。 解: 由表达式可知,当n<0时,单位样值响应h(n)=0,则该系统具有因果性。 根据指数序列的求和公式可知,和是否存在与a的数值有关。若|a|<1,则∑+∞n=0|a|n=11-|a|,系统是稳定的; 若|a|≥1,则∑+∞n=0|a|n发散,系统是不稳定的。 5.5零状态响应的卷积法 已知序列e(n)=1,2↑,3,5,若用单位样值序列表示,则为 e(n)=δ(n+1)+2δ(n)+3δ(n-1)+5δ(n-2) =e(-1)δ(n+1)+e(0)δ(n)+e(1)δ(n-1)+e(2)δ(n-2) 同理,任意序列可以表示为 e(n)=∑+∞m=-∞e(m)δ(n-m)(5.51) 式(5.51)说明任意序列可分解为不同时刻不同幅度的单位样值序列之和。因而求任意序列通过系统所产生的零状态响应,可以先求出各单位样值序列产生的响应,然后利用系统的线性时不变特性,求出最终的零状态响应。具体分析过程如下。 已知单位样值序列δ(n)通过系统所产生的零状态响应为h(n),即 δ(n)→h(n) 根据系统的时不性,可得 δ(n-m)→h(n-m) 再根据系统的齐次性,有 e(m)δ(n-m)→e(m)h(n-m) 最后由系统的叠加性,可得 ∑+∞m=-∞e(m)δ(n-m)→∑+∞m=-∞e(m)h(n-m)(5.52) 由此可得 e(n)→∑+∞m=-∞e(m)h(n-m)(5.53) 即激励e(n)作用于系统时,产生的零状态响应为 rzs(n)=∑+∞m=-∞e(m)h(n-m)(5.54) 式(5.54)中激励与单位样值响应之间的运算称为卷积和。 5.5.1卷积和的定义及计算 任意序列f1(n)和f2(n)进行如下运算: y(n)=∑+∞m=-∞f1(m)f2(n-m)(5.55) 式(5.55)称为卷积和(简称为卷积),卷积和的结果为相同自变量的序列y(n)。卷积和运算与连续时间信号中的卷积积分具有类似的运算规律,区别在于此处是对离散量的叠加,因而将积分运算改成了求和运算,卷积和的运算符仍用“”来表示。根据卷积和的定义,式(5.54)可以写为 rzs(n)=e(n)h(n)(5.56) 式(5.56)表明系统零状态响应等于激励和系统单位样值响应的卷积,可用图524来描述。 图524信号经过系统框图 根据卷积的定义,卷积的计算过程可分为换元、反褶、移位、相乘及求和这几个步骤。卷积的计算方法包括图解法、解析式法以及对位相乘求和法等。 1. 图解法 根据式(5.55),从波形运算的角度可以看出,卷积包含以下步骤。 (1) 换元: 将两个信号的自变量由n变为m,得到f1(m)和f2(m)的波形; (2) 反褶: 将f2(m)的波形反褶得到f2(-m)的波形; (3) 移位: 对f2(-m)波形右移n个单位,得到f2(n-m)的波形; (4) 相乘: 将f1(m)和f2(n-m)进行对位相乘; (5) 求和: 将对位相乘后的各位乘积求和。 由于n的不同,相乘的两个函数会发生变化,要根据n的变化重复步骤(3)~(5),故两个序列卷积后还是一个关于n的序列。 下面通过例题来看一下卷积和的图解过程。 例520如图525所示的两个离散时间信号f1(n)和f2(n),用图解法计算卷积和y(n)=f1(n)f2(n)。 图525例520序列波形 解: 图解法的具体步骤如下。 (1) 换元: 将自变量n替换为m,得到f1(m)和f2(m)的波形,如图526(a)和(b)所示。 图526例520图解法示意图 (2) 反褶: 将f2(m)反褶,得到f2(-m)的波形,如图526(c)所示。 (3) 移位、相乘并求和: 将f2(-m)沿横轴平移,得出f2(n-m)的波形,再对f1(m)与f2(n-m)进行对位相乘及求和。f2(n-m)的位置取决于n的取值。 当n<0时,f2(-m)沿横轴向左移动,波形如图527(a)所示。此时f1(m)和f2(n-m)无重叠,f1(m)f2(n-m)=0,故y(n)=0; 图527例520图解法示意图 当n=0时,f2(n-m)即为f2(-m),波形如图527(b)所示。 y(0)=∑+∞m=-∞f1(m)f2(-m)=f1(0)f2(0)=6 当n=1时,f2(1-m)波形如图527(c)所示。 y(1)=∑+∞m=-∞f1(m)f2(1-m)=f1(0)f2(1)+f1(1)f2(0)=2×1+1×3=5 同理可得,y(2)=23,y(3)=12,y(4)=21,y(5)=5。 当n>5时,x1(m)和x2(n-m)无重叠,故y(n)=0。 综上,卷积和可表示为 y(n)={6↑,5,23,12,21,5} 由以上过程可以看出,当一个序列固定不动,另一个序列从左至右沿横轴进行移动过程中,需要将两个序列在重叠位置上的信号值进行相乘,然后再将所有乘积求和,即可以得到该移动值n所对应的卷积和。 2. 解析式法 解析式法是直接将序列表达式代入式(5.55)进行计算的方法。下面举两个例题具体 看一下解析式法的过程。 例521已知f1(n)=f2(n)=u(n),求y(n)=f1(n)f2(n)。 解: 根据卷积和定义 y(n)=∑+∞m=-∞f1(m)f2(n-m)=∑+∞m=-∞u(m)u(n-m) 当m≥0且n-m≥0,即0≤m≤n时,u(m)u(n-m)=1,可得 y(n)=∑nm=01u(n)=(n+1)u(n) 例522已知f(n)=2nu(n),h(n)=u(n-1),求y(n)=f(n)h(n)。 解: y(n)=f(n)h(n)=∑+∞m=-∞f(m)h(n-m)=∑+∞m=-∞2mu(m)u(n-m-1) 当m≥0且n-m-1≥0,即0≤m≤n-1时,u(m)u(n-m)=1,可得 y(n)=∑n-1m=02mu(n-1)=1-2n1-2u(n-1)=(2n-1)u(n-1) 在解析式法求卷积和的过程中,要注意利用阶跃序列的特性,对求和范围加以限制以便于简化运算。 3. 对位相乘法 对位相乘法巧妙地将图解过程转化为“对位排列”,大大简化了有限长序列的卷积过程。 例523利用对位相乘法求解例520。 解: 采用左对齐或右对齐的方式将两个序列的值进行排列,然后将序列值逐个对位相乘,最后把同一列上的乘积对位求和,和值不进位,即可得y(n)。下面以右对齐为例,介绍求解过程。 f1(n):2141 f2(n): 315 105205 2141 63123 652312215 f1(n)的起点为n=0,f2(n)的起点同样为n=0,卷积和y(n)的起点为两序列的起点序列号之和,则y(n)的起点也为n=0,故y(n)={6↑,5,23,12,21,5}。一般地,长度分别为N和M的两个有限长序列进行卷积后,卷积和的长度为N+M-1。 由对位相乘法可直观看出任意序列f(n)与δ(n)卷积还是任意序列本身。 即 f(n)*δ(n)=f(n) 图解法、解析式法以及对位相乘法各有优缺点。图解法和对位相乘法形象直观,而解析式法则直接利用表达式代入卷积和的定义进行计算。在计算卷积和的过程中要注意结合序列的特点灵活选择计算方法。 5.5.2卷积和的性质 卷积和存在一些运算规律可以简化运算过程。 设序列f1(n) 、f2(n) 和f3(n) 中任意两序列的卷积和均存在,则卷积和运算存在下列运算规律。 1. 移序 若有f(n)=f1(n)*f2(n) ,则有 f(n)=f1(n-m)f2(n+m)(5.57) f(n+m)=f1(n)f2(n+m)=f1(n+m)f2(n)(5.58) f(n-m)=f1(n)f2(n-m)=f1(n-m)f2(n)(5.59) 证明: 根据定义式,有 f1(n-m)f2(n+m)=∑+∞k=-∞f1(k-m)f2(n-k+m)(5.510) 令式(5.510)中k-m=λ,则n-k+m=n-λ,当k=-∞时,λ=-∞,当k=∞时,λ=∞,故上式改写为 f1(n-m)f2(n+m)=∑+∞λ=-∞f1(λ)f2(n-λ)=f1(n)f2(n)(5.511) 式(5.58)和式(5.59)的证明可以采用类似的方法,在此就不一一加以证明了。 2. 交换律 f1(n)f2(n)=f2(n)f1(n)(5.512) 该性质表明: 交换两个序列的位置,不改变卷积和的运算结果。在计算卷积时,可利用该性质,将形式相对复杂的序列自变量替换为m,形式相对简单的序列自变量替换为n-m,如此可以一定程度地减少运算量和计算复杂度。 证明: f1(n)f2(n)=∑+∞m=-∞f1(m)f2(n-m)(5.513) 令式(5.513)中n-m=λ,则m=n-λ,当m=-∞时,λ=∞,当m=∞时,λ=-∞,上式可改写为 f1(n)f2(n)=∑+∞λ=-∞f2(λ)f1(n-λ)=f2(n)f1(n)(5.514) 证毕。 3. 结合律 f1(n)f2(n)f3(n)=f1(n)[f2(n)f3(n)](5.515) 由式(5.515)可以看出,要计算多个序列的卷积和,可以改变卷积的运算次序,先计算其中任意两个序列的卷积和,再与第三个序列相卷积,不改变卷积和的最终计算结果。 将结合律运用到系统分析中,其物理意义在于,由若干LTI系统级联构成的复合系统,可从系统的角度将其等效成为一个系统,该系统的单位样值响应等于所有级联子系统单位样值响应的卷积,如图528所示。 图528级联系统方框图 4. 分配律 f1(n)[f2(n)+f3(n)]=f1(n)f2(n)+f1(n)f3(n)(5.516) 该性质说明,先求和再卷积的运算次序可以转变为先卷积再求和。分配律的物理意义在于,由若干LTI系统并联构成的复合系统,可以等价为一个系统,该系统的单位样值响应等于所有并联子系统单位样值响应之和,如图529所示。 图529并联系统方框图 例524已知f1(n)=2δ(n)+δ(n-3),f2(n)=3δ(n)+5δ(n-1),求y(n)=f1(n)f2(n)。 解: y(n)=f1(n)f2(n)=[2δ(n)+δ(n-3)][3δ(n)+5δ(n-1)] =2δ(n)3δ(n)+δ(n-3)3δ(n)+2δ(n)5δ(n-1)+δ(n-3)5δ(n-1) =6δ(n)+10δ(n-1)+3δ(n-3)+5δ(n-4) 例525某复合系统的结构如图530所示,已知各子系统的单位样值响应分别为h1(n)=δ(n-1),h2(n)=u(n)-u(n-3),求该复合系统的h(n)。 图530例525系统框图 解: 根据单位样值响应的定义,当激励为δ(n)时,所产生的零状态响应为h(n)。 设e(n)=δ(n),则产生的响应r(n)为单位样值响应h(n),故 r(n)=h(n)=[δ(n)+δ(n)h1(n)]h2(n) =h2(n)+h1(n)h2(n)=u(n)-u(n-3)+δ(n-1)[u(n)-u(n-3)] =[u(n)-u(n-3)]+[u(n-1)-u(n-4)] =δ(n)+2δ(n-1)+2δ(n-2)+δ(n-3) 利用卷积和的分配律和结合律可以将复杂系统等效为一个整体,利用总的单位样值响应来表征系统的特性,进而统一利用卷积和运算来描述输入、系统单位样值响应及输出三者之间的关系,在已知输入和系统单位样值响应的条件下就可以实现系统零状态响应的分析。 5.5.3零状态响应的求解 由前述可知,系统的零状态响应等于激励与系统单位样值响应的卷积和。因而若已知系统模型,要求任意激励所产生的零状态响应,需先求出系统的单位样值响应,然后利用激励与单位样值响应的卷积和运算即可求出零状态响应。 例526某LTI离散时间系统的系统框图如图531所示,激励e(n)=2nu(n),求当n≥0时系统的零状态响应r(n)。 图531例526系统框图 解: (1) 列写系统的差分方程。 根据加法器的输入输出关系,可得 r(n)=e(n)-3r(n-1)-2r(n-2) 整理得 r(n)+3r(n-1)+2r(n-2)=e(n) (2) 求系统的单位样值响应。 利用移序算子,写出单位样值响应所对应的算子方程 (1+3E-1+2E-2)h(n)=δ(n) h(n)=11+3E-1+2E-2δ(n)=21+2E-1δ(n)-11+E-1δ(n) 故有 h(n)=[2(-2)n-(-1)n]u(n) (3) 求系统的零状态响应r(n)。 r(n)=e(n)h(n)=∑+∞m=-∞[2(-2)m-(-1)m]u(m)2n-mu(n-m) =2n∑nm=02(-2)m12m-∑nm=0(-1)m12mu(n) =2n21-(-1)n+11-(-1)-1-(-1/2)n+11-(-1/2)u(n) =132n+(-2)n-13(-1)nu(n) 例527已知描述某系统的差分方程r(n)-0.9r(n-1)=0.05e(n),e(n)=u(n),其中r(-1)=1,求当n≥0时系统的响应r(n)。 解: (1) 求零输入响应rzi(n)。 系统的齐次差分方程为 r(n)-0.9r(n-1)=0 由齐次差分方程,列写特征方程得 λ-0.9=0 由特征根λ=0.9,写出零输入响应的形式为 rzi(n)=c·0.9n 将r(-1)=1代入系统的齐次差分方程,可得 r(0)-0.9r(-1)=0 解得r(0)=0.9,将其代入零输入响应的表达式,可得 rzi(0)=c·(0.9)0=0.9 解得c=0.9,故可得零输入响应为 rzi(n)=0.9n+1u(n) (2) 求系统的单位样值响应h(n)。 求单位样值响应时,系统模型为 h(n)-0.9h(n-1)=0.05δ(n) 由于h(n)具有与零输入响应相同的形式,故可表示为 h(n)=k·0.9nu(n) h(-1)=0 将h(-1)=0代入求单位样值响应时的系统模型,得 h(0)-0.9h(-1)=0.05 解得h(0)=0.05,将其代入单位样值响应的表达式,得 k·0.90=0.05 解得k=0.05,故单位样值响应为 h(n)=0.05·0.9nu(n) (3) 求零状态响应rzs(n)。 在e(n)作用下,系统的零状态响应为 rzs(n)=e(n)h(n)=∑+∞m=-∞0.05·0.9mu(m)u(n-m)=0.05∑nm=00.9mu(n) =0.05×1-0.9n+11-0.9u(n)=0.5(1-0.9n+1)u(n) (4) 求系统的全响应r(n)。 r(n)=rzi(n)+rzs(n)=0.9n+1u(n)+0.5(1-0.9n+1)u(n) =(0.45×0.9n+0.5)u(n) 求解离散时间系统响应时,既可以采用经典法求解差分方程,也可以将系统响应分为零输入响应和零状态响应,再分别求解。其中利用卷积和求零状态响应的方法,物理意义明确,同时也是变换域分析法求零状态响应的理论基础。后续将在时域分析的基础上变换分析问题的角度,通过序列及系统差分方程的z变换,将时域差分方程转换为z域代数方程,将卷积和运算转换为相乘运算,简化系统分析的过程,具体内容将在下一章进行讨论。 习题5 51试画出下列序列的图形: (1) f1(n)=n+1,-30 -12n,n≤0 ; (4) f4(n)={3,1↑,2,-5,4}。 52试画出下列序列的图形: (1) f1(n)=12n-1u(n); (2) f2(n)=12n[δ(n+1)-δ(n)+δ(n-1)]; (3) f3(n)=-12n-1[u(n)-u(n-3)]; (4) f4(n)=(-2)nu(-n); (5) f3(n)=2n-1u(n-1)。 53写出图532所示各序列的表达式。 图532题53图 54设序列f(n)={2,5↑,3,1,-1},请画出下列各序列的波形图: (1) f1(n)=f(n+2); (2) f2(n)=f(-n+1); (3) f3(n)=f(n+2)+f(n-2); (4) f4(n)=f(1-n)+f(n+1); (5) f5(n)=f(n)·f(1-n); (6) f6(n)=f(2n); (7) f7(n)=fn2; (8) f8(n)=f(2n)·f(1-n)。 55请绘出下列序列的图形: (1) f1(n)=cosnπ2[u(n)-u(n-8)]; (2) f2(n)=(-1)nu(2n); (3) f3(n)=(n-1)u(n); (4) f4(n)=nu(n-1); (5) f5(n)=n[δ(n+1)+δ(n)+δ(n-3)]; (6) f6(n)=R4(n+3); (7) f7(n)=u(-n+5)-u(-n-2); (8) f8(n)=nR4(n+1)。 56判断下列各序列是否是周期序列。若是,请求出其周期N。 (1) f1(n)=sinnπ2+π6; (2) f2(n)=cosn; (3) f3(n)=ejπ4n+π3; (4) f4(n)=cos43πn+cos17n; (5) f4(n)=cos43πn-sin45πn。 57试绘出下列离散时间系统的系统框图: (1) r(n)=e(n)+e(n-2); (2) r(n)+3r(n-1)+2r(n-2)=e(n)+e(n-1)。 58试写出图533所示各系统的差分方程,并指出系统的阶次。 图533题58图 59对于下列系统,试判断系统是否具有线性和时不变性。 (1) r(n)=e(n-n0); (2) r(n)=2e(n)+1; (3) r(n)=[e(n)]2; (4) r(n)=∑nm=-∞e(m); (5) r(n)=ee(n); (6) r(n)=e(n)sinπ3n+π6。 510某人向银行贷款20万,采取逐月计息偿还方式。从贷款后第1个月开始每月定时定额还款,每月固定还款金额0.4万元,贷款月息为β。设第n个月的欠款额为y(n),试写出y(n)满足的方程。 511同一平面上有k条直线,均两两相交,但没有三条或三条以上直线交于同一点。问: 满足这一条件的k条直线能把平面分成多少块? 512已知系统差分方程及边界条件,求n≥0时系统零输入响应。 (1) r(n)-13r(n-1)=0,r(0)=1; (2) r(n)-2r(n-1)=0,r(0)=12; (3) r(n)+3r(n-1)+2r(n-2)=0,r(-1)=1,r(-2)=1; (4) r(n)+2r(n-1)+r(n-2)=0,r(-1)=r(-2)=2; (5) r(n)+r(n-2)=0,r(0)=1,r(1)=3。 513已知r(n)-r(n-1)=n,r(-1)=1,求n≥0时系统输出r(n)。 514已知r(n)+3r(n-1)+2r(n-2)=u(n),r(-1)=r(-2)=1,求n≥0时系统输出r(n)。 515求下列各离散时间系统的单位样值响应h(n): (1) r(n)-2r(n-1)=2e(n); (2) r(n)-0.5r(n-1)+0.06r(n-2)=e(n-1); (3) r(n)+3r(n-1)+2r(n-2)=e(n)-e(n-1); (4) r(n)+14r(n-1)-18r(n-2)=e(n); (5) r(n)=e(n)-e(n-1)-2e(n-3)+e(n-4)。 516已知下列LTI离散时间系统的单位样值响应h(n),试判断系统的因果性和稳定性,并简要说明理由: (1) h(n)=0.2nu(n); (2) h(n)=δ(n-1); (3) h(n)=nu(n); (4) h(n)=2nR5(n); (5) h(n)=3nu(-n); (6) h(n)=-12n+110nu(n)。 517已知某LTI离散时间系统的阶跃响应为g(n)=[(-2)n+1]u(n),求该系统的单位样值响应h(n)。 518已知某LTI离散时间系统的单位样值响应和激励分别如图534(a)、(b)所示,画出系统响应r(n)的波形。 图534题518图 519已知某LTI离散时间系统的系统框图如图535所示,各子系统分别为h1(n)=u(n),h2(n)=δ(n-2),h3(n)=u(n-1),求该系统的单位样值响应。 图535题519图 520已知系统差分方程为r(n)=e(n)-2e(n-1)-e(n-2),激励e(n)=u(n)-u(n-2)。 求: (1) 系统单位样值响应h(n); (2) n≥0时零状态响应r(n)。 521求下列LTI离散时间系统的零输入响应rzi(n)、零状态响应rzs(n)及全响应r(n)。 (1) r(n)+13r(n-1)=e(n),r(-1)=-1,e(n)=12nu(n); (2) r(n)-56r(n-1)+16r(n-2)=e(n)-e(n-1),r(-1)=0,r(-2)=1,e(n)=u(n)。 522已知描述某LTI离散时间系统的差分方程为 r(n)+5r(n-1)+6r(n-2)=e(n)-e(n-1) (1) 求系统的单位样值响应; (2) 判断系统的稳定性; (3) 请画出该系统的模拟图。 523某离散时间系统的模拟图如图536所示。 (1) 请写出系统的差分方程; (2) 求h(n); (3) 若激励e(n)=u(n),r(-1)=0,r(-2)=1,求响应r(n)。 图536题523图 524已知二阶LTI离散时间系统的单位样值响应为 h(n)=(1+2·3n)u(n) (1) 请写出该系统的差分方程; (2) 判断系统的因果、稳定性; (3) 若激励为e(n)=u(n)-u(n-2),求系统的零状态响应r(n)。