离散信号与系统的时域分析

1、引引 言言连续时间系统连续时间系统：这类系统用于传输和处理连续时间信号离散系统离散系统：用于传输和处理离散时间信号的系统称为离散时间系统，数字计算机是典型的离散系统例子，数据控制系统和数字通信系统的核心组成部分也都是离散系统。 混合系统混合系统：连续系统与离散系统组合起来使用。5.1 5.1 离散时间基本信号离散时间基本信号5.1.1 5.1.1 离散时间信号离散时间信号 连续时间信号，在数学上可以表示为连续时间变量t的函数。这类信号 的特点是：在时间定义域内，除有限个不连续点外， 对任一给定时刻都对应有确定的信号值。 离散时间信号，简称离散信号，它是离散时间变量tk(k=0,1, 2, )的

2、函数。信号仅在规定的离散时间点上有意义，而在其它时间则没有定义，如图 5.1-1(a)所 示。鉴于tk按一定顺序变化时，其相应的信号值组成一个数值序列，通常把离散时间信号定义为如下有序信号值的集合：fk=f(tk) k=0, 1, 2, (5.1-1)式中，k为整数，表示信号值在序列中出现的序号。图 5.1 1 离散时间信号f (tk)t-3t-2t-1t1ot2t3tk(a)f (kT) 3TokT(b) 2T TT 2T 3Tf (k) 3ok(c) 211 2 3 式(5.1-1)中tk和tk-1之间的间隔(tk-tk-1)可以是常数，也可以随k变化。在实际应用中，一般取为常数。例如，对

3、连续时间信号均匀取样后得到的离散时间信号便是如此。对于这类离散时间信号，若令tk-tk-1=T，则信号仅在均匀时刻t=kT(k=0,1,2,)上取值。此时，式(5.1 - 1)中的f(tk)可以改写为f(kT)，信号图形如图 5.1-1(b)所示。 为了简便，我们用序列值的通项f(kT)表示集合f(kT)，并将常数T省略，则式(5.1-1)可简写为fk=f(k) k=0, 1, 2, (5.1-2)工程应用中，常将定义在等间隔离散时刻点上的离散时间信号称为离散时间序列离散时间序列 ，简称序列序列。5.1.2 5.1.2 离散时间基本信号离散时间基本信号1. 1. 单位脉冲序列单位脉冲序列单位脉

4、冲序列定义为单位脉冲序列定义为01)(k00kk图 5.1 2 单位脉冲序列(k) 1ok 2121位移单位脉冲序列01)(0kk00kkkk或01)(0kk00kkkk(a)(kk0)okk011k0k01o k011 k0 k01(b)k(k k0)图5.1-3 移位单位脉冲序列2. 2. 正弦序列正弦序列正弦序列的一般形式为)cos()(0kAkf由于)(cos2cos)2cos()cos()(00000NkAmkAmkAkAkf式中，m、N均为整数。式(5.1-5)表明，只有当 为整数，或者02mNmN02为有理数时，正弦序列才是周期序列；否则为非周期序列。(5.1 - 6) 当正弦序

5、列是通过抽取连续时间正弦信号的样本获得时， 如果假设正弦信号 的周期为T0，取样间隔为Ts，那么，经过抽样得到的正弦序列可表示为)cos(0t)cos(2cos)cos()(000kkTTtkfskTts式中， ， 将它代入式(5.1 - 6)可 得002TTsmNTTs002 对于连续时间正弦信号 , 按几种不同间隔Ts抽样得到的正弦序列示于图 5.1-4 中。当 时，有此时， ， 是一个周期为 16 的周期性正弦序列，其 图形如图 5.1-4(a)所示。当 ，可得到如图 5.1 - 4(b)所示的序列，其 ，是一个周期为23 的周期性正弦序列。 当 ，序列图形如图5.1 - 4(c)所示，

6、其 ，由于 ,是一无理数，故f(k)是一非周期正弦序列，值得注意的是此时它的包络函数f(t)仍具有周期性。ttf2cos)(161sT162cos)2cos(2cos)(kkTttfskTts16206cos)(121kkfTs时，238cos)(234kkfTs时，238061012200sTT图 5.14 正弦序列k(a)(b)kk(c)3. 3. 指数序列指数序列指数序列的一般形式为kAekf)( (1)若A和 均为实数，且设 则 为实指数序列。 当a1时，f(k)随k单调指数增长。当0a 1时，f(k)随k单调指数衰减； 当a-1时，f(k)的绝对值随k按指数规律增长。 当-1a0时，

7、f(k)绝对值随k按指数 规律衰减。 且两者的序列值符号呈现正、负交替变化； 当a=1时，f(k)为常数序列。当a=-1时，f(k)符号也呈现正、 负交替变化。ekAekf)(图 5.1 5 实指数序列(b) 6 4 2o246f (k)1k(a) 6 4 2o246f (k)1kof (k)of (k)01k 6 4 224610(d)(c)(e) 4 20246f (k)1k 6 4 20246f (k)1k8(f) 6 4 2246k(2) 若A=1，=j0，则 kjekf0)(是虚指数序列。 我们已经知道，连续时间虚指数信号e j0t是周期信号。然而，离散 时间虚指数序列ej0k则只有

8、满足一定条件时才是周期的， 否则是非周 期的。根据欧拉公式，式(5.1 - 9)可写成 kjkekj00sincos0可见，e j0k的实部和虚部都是正弦序列，只有其实部和虚部同时为周 期序列时，才能保证ej0k是周期的。 (3) 若A和均为复数，则f(k)=Aek为一般形式的复指数序列 。设复数A=|A|ej, =+j0，并记e=r， 则有 )sin()cos()(00)()()(000kjkrAerAeeAeeAAekfkkjkjkkjjk可见，复指数序列f(k)的实部和虚部均为幅值按指数规律变化的正弦序列。 图 5.1 6 复指数序列 4. Z序列序列Z序列的一般形式为 kzkf)(式中

9、，z为复数。通常，称序列值为复值的序列为复序列复序列。显然， Z 序列是一复序列。若将z表示为极坐标形式 0jezzkjkjkezezzkf00)()()sin(cos)(00kjkzzkfkk根据欧拉公式， 还可写成 5.2 卷卷 积积 和和5.2.1 卷积和的定义卷积和的定义 定义两个连续时间信号f1(t)和f2(t)的卷积运算为 dtfftftf)()()()(2121同样地， 我们定义 iikfifkfkfkf)()()()()(2121为序列f1(k)和f2(k)的卷积和运算，简称卷积和 ( Convolution Sum)。 (5.2 - 2) 如果f1(k)为因果序列，由于k0时

10、，f1(k)=0，故式(5.2 - 2)中求和下限可 改写为零，即 02121)()()()(iikfifkfkf 如果f2(k)为因果序列，而f1(k)不受限制，那么式(5.2 - 2)中，当(k-i) 0，即ik时，f2(k-i)=0, 因而和式的上限可改写为k，也就是 kiikfifkfkf)()()()(2121如果f1(k)和f2(k)均为因果序列， 则有 (5.2 - 5)kiikfifkfkf02121)()()()(考虑到f1(k)、f2(k)均为因果序列，根据式(5.2 - 5)，可将上式表示为 iiikiekfkf)()()()(121例例 5.2 1 设f1(k)=e-k

11、( k)，f2(k)=(k), 求f1(k)*f2(k)。解解 由卷积和定义式(5.2 - 2)得 1)1(1100211111)()()(eeeeeeikekfkfkkkiiii显然，上式中k0，故应写为 )(11)()()()(1)1(21keekkekfkfkk 与卷积运算一样，用图解法求两序列的卷积和运算也包括信号的翻转、平移、相乘 、求和等四个基本步骤。 例例 5.2 2 已知离散信号 0231)(1kf其他210kkk 04)(2kkf其他3 , 2 , 1 , 0k求卷积和f1(k)*f2(k)。 解解 记卷积和运算结果为f(k)，由式(5.2 - 2)得 iikfifkfkfk

12、f)()()()()(2121 第一步，画出f1(i)、f2(i)图形，分别如图 5.2 - 1(a)、 (b)所示。 第二步，将f2(i)图形以纵坐标为轴线翻转 180，得到f2(-i)图形，如图 5.2 - 1(c)所示。 第三步，将f2(-i)图形沿i轴左移(k0)或右移(k0)|k|个时间单位，得到f2(k-i) 图形。例如，当k=-1和k=1时，f2(k-i)图形分别如图 5.2 - 1(d)、 (e)所示。 第四步第四步，对任一给定值k，按式(5.2 - 6)进行相乘、求和运算，得到序号为k的卷 积 和序列值f(k)。若令k由-至变化，f2(k-i)图形将从-处开始沿i轴自左向右移

13、动 ，并由式(5.2 - 6)计算求得卷积和序列f(k)。对于本例中给定的f1(k)和f2(k) ，具体计算过程如下： 1023if1(i)13214(a) 1023if2(i)3214(b)451 1 0f2( i)321(c)4i1 2 3 4 1 0f2( 1 i)321(d)4i1 2 3 4 5 1 0f2(1 i)321(e)4i1 2 3 42 1023kf(k)151914(f)51347267图5.2-1 卷积和计算于是，其卷积和为 002713191540)(kkf 对于两个有限长序列的卷积和计算， 可以采用下面介绍的更为简便实用的方法计算。 这种 方法不需要画出序列图形，

14、 只要把两个序列排成两行，按普通乘法运算进行相乘， 但中 间结果不进位，最后将位于同一列的中间结果相加得到卷积和序列。 例如，对于例5. 2 - 2 中给定的f1(k)和f2(k)，为了方便，将f2(k)写在第一行， f1( k)写在第二行， 经序列值相乘和中间结果相加运算后得到 例：一个系统，其单位抽样响应为： 其输入序列为： ，求输出响应000)(nnanhn)()()(Nnununx)(ny 解：讨论解析法求卷积和。首先要根据卷积和的变化情况，按转折点划段，然后对每段的卷积和确定上下限。确定上下限的原则是：若给定两序列 和 的非零值的下限分别为 ，上限分别为 ，则选 中大者作为卷积和的下

15、限，选 中小者作为卷积和的上限。 )(kx)(knh21LL 和21VV 和21LL 和21VV 和5.2.2 卷积和的性质卷积和的性质 性质性质1 离散信号的卷积和运算服从交换律、结合律和分配律，即)()()()()()()()()()()()()()()()()(31213213213211221kfkfkfkfkfkfkfkfkfkfkfkfkfkfkfkfkf 性质性质 2 任一序列f(k)与单位脉冲序列(k)的卷 积和等于序列f(k)本身， 即 )()()()()(kfkfkkkf性质性质 3 若f1(k)*f2(k)=f(k)，则 )()()()()(1211121kkfkfkkf

16、kkfkf)()()()()(2112212211kkkfkkfkkfkkfkkf式中k1 , k2均为整数。 例例 5.2-3 已知序列x(k)=(3)-k(k) ，y(k)=1, -k, 试验证x(k)和y(k)的卷积和运算满足交换律，即 )()()()(kxkykykx 解解 先计算x(k)*y(k)，考虑到x(k)是因果序列，根据式(5.2-3)，有 5 . 1233111) 3(1)() 3()()()()(0iiiiiiikyixkykx再计算y(k)*x(k)，同样考虑到x(k)是因果序列，可得 求解过程中对k没有限制，故上式可写为x(k)*y(k)=y(k)*x(k)=1.5

17、-k可见，x(k)*y(k)运算满足交换律。 5 . 1313) 3() 3(3) 3() 3()() 3(1)()()()()()(kkkikiikikiikiikikxiykxky所以 5 . 1)()()()(kxkykykx 例例 5.2-4 已知序列f1(k)=2-(k+1) (k+1)和f2(k)=(k-2)，试计算卷积和f1(k)*f2(k)。 解解 用下面两种方法计算。 方法一：方法一：图解法。将序列f1(k), f2(k)的自变量换为i，画出f 1(i)和f2(i)的图形如图 5.2-2(a), (b)所示。 将f2(i)图形翻转 180后，得f2(-i)，如图5.2-2(c

18、)所示。 当k1时，由图 5.2-2(d)可知，其乘积项f1(i)f2(k-i)为零，故f1(k)*f2(k)=0。 图 5.2-2 f1(i)0i(a)11 2 3 4 51f2(i)0i(b)1 2 3 4 516f2( i)0(c)21i46f2(k i)0(d)2i424(k1)f2(k i)0(e)2i24(k1)f1(i)f1(i)当k1时，按卷积和定义，参见图 5.2-2(e)，可得 )21 (2222212222222)2(2)2() 1(2)()(11)2(121121)1(1)1()1(21kkkkiikiiiiiiikikikfkf于是)21 (20)()(21kkfkf

19、11kk故有) 1()21 (2)()(21kkfkfk方法二方法二： 应用卷积和性质 3。 先计算 kkkiiiikikikkkf22212212)()()2()()()2()(110上式中k0， 故有 )()22()()(2)(kkkkfkk再应用卷积和性质 3，求得 ) 1()21 (2) 1(22) 1()21()2() 1(2)()()1()1(21kkkfkfkkkfkfkkk5.2.3 常用序列的卷积和公式常用序列的卷积和公式 表表 5.1 常用序列的卷积和公式常用序列的卷积和公式 5.3 离散系统的算子方程离散系统的算子方程 5.3.1 LTI离散时间系统离散时间系统 图图 5

20、.3-1 离散系统的输入输出模型离散系统的输入输出模型 )()()()(kykfkfTky离散系统f (k)y(k) 离散时间系统的状态和状态变量。离散时间系统在k0时刻的状态是指 满足如下条件的数目最少的一组数据x1(k0), x2(k0), , xn(k0)。 这组 数据连同k0k上的输入f(k)就可以惟一地确定k时刻的输出y(k)，而不需具体知道k 0以前的输入情况。n称为离散系统的阶数。 在实际工作过程中，系统的状态x1(k0), x2(k0), , xn(k0)随k0不同 而变化，我们把描述系统状态变化的变量称作状态变量， 它是一组序 列信号，记为x1(k), x2(k), , xn

21、(k)。 离散时间系统的零输入响应、零状态响应和完全响应。设k0为初始观察 时刻，则可将系统的输入区分为两部分，称k0以前的输入为历史输入信号，称k0及k0以后的输入为当前输入信号或简称输入信号。我们将仅由k0时刻的初始状态或历史输入信号引起的响应称作零输入响应，记为yx(k)；仅由当前输入信号引起的响应称作零状态响应，记为yf(k)。而将零输入响应、零状态响应之和 称作系统的完全响应，记为y(k)。 离散时间系统的齐次性、叠加性和线性特性。设离散系统的输入输出关系为f(k) y(k) 所谓齐次性是指对于任意常数a、 输入f(k)和输出y(k)，恒有af(k) ay(k) (5.3-3) 所谓

22、叠加性叠加性是指对于输入f1(k)、f2(k)和输出y(k)，若设f1(k)y1(k)，f2(k) y2( k)，则恒有f1(k), f2(k) y1(k)+y2(k) (5.3- 4)式中，f1(k), f2(k)表示f1(k)和f2(k)同时作为系统的输入。 齐次性和叠加性统称为线性特性。 对于任意常数a和b，输入f1(k)和 f2(k)共同作用时，系统的线性特性可表示为af1(k), bf2(k) ay1(k)+by2(k) (5.3 - 5)它同时体现了式(5.3-3)的齐次性和式(5.3-4)的叠加性。 线性离散时间系统和非线性离散时间系统。 若离散时间系统的响应可 分 解为零输入响

23、应和零状态响应两部分， 且零输入响应与初始状态或历史输入信号、 零状态 响应与当前输入信号之间分别满足齐次性和叠加性，则称该系统为线性离散时间 系统，否则称为非线性离散时间系统。 时不变离散时间系统和时变离散时间系统。 设离散时间系统的输入输出关系为 )()(kykfT若对于任意整数k0， 恒有 )()(00kkykkfT则称该系统为时不变离散时间系统时不变离散时间系统，否则称为时变离散时间时变离散时间系系 统。统。 因果离散时间系统和非因果离散时间系统。 如果系统始终不会在 输入加入之前产生响应， 这种系统称为因果系统因果系统， 否则称为非因果系统。非因果系统。 例如，有三个系统的输入输出关

24、系如下： 系统 1 y(k)=kf(k) 系统 2 y(k)=|f(k)| 系统 3 y(k)=2f(k)+3f(k-1) 根据定义容易验证： 系统 1 是线性时变离散时间系统， 系统 2 是非线性时不变离散时间 系统， 而系统 3 是线性时不变离散时间系统。 根据第 1 章讨论结果，一个n阶线性时不变离散时间系统，若其输入为f(k)，全响应为y (k)，那么，描述该系统输入输出关系的数学模型是n阶线性常系数差分方程，它可以表 示为 )() 1()()() 1()(0101mkfbkfbkfbnkyakyakymmn式中，ai(i=0, 1, , n-1)，bj(j=0, 1, , m)均为常

25、数。 (5.3-7)5.3.2 离散系统算子方程离散系统算子方程 在连续时间系统分析中，我们曾用微分算子p和积分算子p-1分别表示对函数的微分 和积分运算。与此类似，在离散系统分析中，我们引入E算子(超前算子)，表示将序列提前一个单位时间的运算；E-1算子(迟后算子 )，表示将序列延迟一个单位时间的运算，即：),1()(),1()(1kfkfEkfkEf)()()()(nkfkfEnkfkfEnn应用中，统称E算子和E-1算子为差分算子。 利用差分算子，可将差分方程式(5.3-7)写成下述形式： )()()()()()(011011kfEbkfEbkfbkyEakyEakymmmnn或写成)(

26、)()()1 (011011kfEbEbbkyEaEammmnn进一步写成)()()()(1)(011011kfEAEBkfEaEaEbEbbkynnmmm式中：nnmmmEaEaEAEbEbbEB0110111)()(若令 011011011011)(1)()()(aEaEbEbEbEEaEaEbEbbEAEBEHnnnmmmmmnnnmmm则式(5.3-9)可表示为 )()()(kfEHky此式称为离散时间系统的算子方程。式中的H(E)称为离散系统离散系统 的传输算子的传输算子。H(E)在离散系统分析中的作用与H(p)在连续系统分析中的作用相同 ，它完整地描述了离散系统的输入输出关， 或者

27、说集中反映了系统对输入序列的传输特 性。例如，设某离散系统的差分方程为 )()() 1(kfkayky以单位延迟算子E-1作用于方程两边后，得到 ）（1) 1()(kfkayky图 5.3-2 用H(E)表示离散系统 H(E)f (k)y(k)根据差分算子的定义，容易证明： )()(1)(1)()(1)(1mnkfkfEEkfEEkfkEfEkfEEnmmn可见，对于同一序列而言，超前算子与迟后算子的作用可以互相抵消， 或者说作用于同 一序列的差分算子公式中，分子分母中的算子公因子允许消去。 例例 5.3-1 设描述某离散时间系统的差分方程为 )()2() 1()(kfkbykayky求其传输

28、算子H(E), 并画出系统的模拟框图和信号流图表示。 解解 写出系统的算子方程为 )()()1 (21kfkybEaE所以，系统的传输算子 2111)(bEaEEH再将算子方程改写成 )()()()(21kybEkyaEkfky图 5.3-3 例 5.3-1图 f (k)E1 ay(k)f (k) ay(k)11E11(a)(b)E1 b bE11例例 5.3-2 某离散时间系统的输入输出算子 方程为 )()()()(kfEBkyEA式中： 201122011)(1)(EbEbbEBEaEaEA试画出系统的模拟框图和信号流图。 解解 如同连续系统那样，选择中间变量x(k)，并令 )(11)()

29、(1)(2011kfEaEakfEAkx则有 )()()()()(20112kxEbEbbkxEBky)()()()(2011kfkxEakxEakx)()()()(20112kxEbkxEbkxbky图 5.3-4 例 5.3-2图 f (k)E1 a1(a)E1 a0 x(k)E1x(k)E2x(k)f (k)E1 a1(b)E1 a0y(k)b0b1b2f (k)1E1E1b0b1b21y(k) a1 a0(c)5.4 离散系统的零输入响应离散系统的零输入响应 根据线性系统定义，系统的完全响应由零输入响应和零状态响应两部分组成。 在连续时间系统的时域分析中，我们从描述系统的微分方程或传输

30、算子H(p)出发，分别求 出系统的零输入响应和零状态响应， 然后把它们叠加起来得到系统的完全响应。 这种做法 同样适用于离散系统的时域分析。 只是在离散时间系统分析中， 我们讨论问题的出发点是 描述系统的差分方程或传输算子H(E)。 此外，求解系统零状态响应时， 与连续时间信号 的卷积积分相对应， 需要进行离散时间信号的卷积和计算。 如前所述，一个描述n阶线性时不变离散时间系统的差分方程，若应用差分算子E，则可 表示为 )()()()1 (011011kfEbEbbkyEaEammmnn或者写为 )()()()()()()(kfEHkykfEBkyEA式中： nnmmmEaEaEbEbbEAE

31、BEH0110111)()()( 根据系统零输入响应的定义，如果假定初始观察时刻为k0，那么，离散系统的零 输入响应就是k0及k0以后的输入为零时，仅由k0以前的输入或k0时刻的状态引起 的响应，常记为yx(k)。 由此可见，在系统差分方程式(5.4-1)中，只需 令输入信号f(k)为零，就可得到求解零输入响应yx(k)的方程，其一般形式为 0)()1 (011kyEaEaxnn0kk 或者简写为 0)()(kyEAx0kk 具体地说，离散系统的零输入响应就是上面齐次差分方程满足给定初始条件yx(0 )，yx(1)，yx(n-1)时的解。 5.4.1 简单系统的零输入响应简单系统的零输入响应

32、如果离散系统传输算子H(E)仅含有单个极点r， 这时式(5.4-6)可表示为 0)()(kyrEx0kk 这是一个一阶齐次差分方程，将上式改写为 0)() 1(krykyxx0kk 于是有 rkykyxx)() 1( 此式表明，序列yx(k)是一个以r为公比的几何级数，它具有以下形式： kxrcky1)(0kk 式中，c1是常数，由系统零输入响应的初始条件确定。上述结果与一阶齐次微分方程 解c1et的形式非常类似，因为当时间t按t=kT离散变化时，其解可改写成c1et=c1ekT=c1(eT)k，令eT=r时，就是差分方程式(5.4-7)的解。 因此，我们有如下结论： kxrckyEEBEH1

33、)()()(0kk 如果系统传输算子仅含有g个单极点r1, r2, , rg，则相应齐次差分方程可写成 0)()(21kyrEEExg）（0kk 显然，满足以下方程 0)()(kyExigi, 2 , 1的解，必定也满足式(5.4-10)。仿照微分方程解结构定理的证明，可导得式(5 .4-10)的解为 kngkkxrcrcrcky2211)(式中, 待定系数值c1, c2, , cg由系统零输入响应的初始条件确定。 于是，有结论 gikiixgrckyrErErEEBEH121)()()()()(0kk 为了考察H(E)含有重极点的情况，我们假定对于一极小值，其系统齐次差分方程为 0)()()

34、(kyrErEx且系统初始条件为 ) 1 ()0(xxyyyx(k)表示为 kkxrcrcky)()(21代入初始条件，有 )(2121rcrccc解得 rc1rc2)()()()()()(1 )(313212131321211rrcrkrarrrcrcrrrrrrrkykkkkkkkkkx现在，令0取极限，使得H(E)的两个极点相重合，于是有 kxrkrky)()(lim10或写成 kxrkccky)()(1110式中： )(11110rcc同样道理，如果传输算子H(E)仅含有r的d重极点，这时系统的齐次差分方程为 kddxrkckckccky)()(112210相应的零输入响应可表示为 0

35、)()(kyExd 式中，常数c0, c1, , cd-1由系统零输入响应的初始条件确定。 因此 kdjjjxdrkckyrEEBEH)()()()()(100kk 5.4.2 一般系统的零输入响应一般系统的零输入响应 设n阶离散时间系统的齐次差分方程为 0)()1 (02211kyEaEaEaxnnn其传输算子H(E)含有g个相异极点r1, r2, , rg，对应的重数分别是d1, d2, , dg。 这里， (d1+d2+dg)=n。显然，若di(i=1, 2, , g)为 1 时， 表示相应的极点ri是单极点。此时式可表示为 0)()()()(2121kyrErErExgddgdn阶LT

36、I离散系统的零输入响应为 10)(idjkijijxirkckygixixkyky1)()(0kk 式中： gi, 2 , 1式中，各待定系数由系统零输入响应yx(k)的初始条件确定。 综上所述，由LTI离散系统传输算子H(E)求零输入响应yx(k)的具体步骤可归纳如下： 第一步，第一步，求解方程A(E)=0，得到H(E)的相异极点r1, r2，, rg及相应的重数d1, d2, , dg。将系统齐次差分方程表示为 0)()(1kyrExgidii第二步第二步，求解方程0)()(kyrExidiigi, 2 , 1得到各极点相应输入响应分量10)(idjkijijxirkckygi, 2 ,

37、1第三步，第三步，写出系统的零输入响应gikijijdjgixixrkckykyi1101)()(0kk 第四步第四步，由零输入响应初始条件确定式（5.4-22）中的各个待定系数cij，并最后求出系统的零输入响应yx(k)。（5.4-22）例例 5.4-1 已知离散时间系统传输算子 2)5 . 0)(3 . 0)(2 . 0(2)(EEEEEH及初始条件yx(0)=12，yx(1)=4.9， yx(2)=2.47，yx(3) =1.371。 求该系统的零输入响应。 解解 因为传输算子H(E)极点为r1=0.2，r2=0.3，r3=0.5(二重极点)。所以，可得 kkxkkxkkxkccrkcc

38、kycrckycrcky)5 . 0()()() 3 . 0()()2 . 0()(3130331303202202101101kkkkxkccccky)5 . 0()5 . 0() 3 . 0()2 . 0()(31302010上式中令k=0, 1, 2, 3, 代入初始条件后得到 371. 1375. 0125. 0027. 0008. 0) 3(47. 25 . 025. 009. 004. 0)2(9 . 45 . 05 . 03 . 02 . 0) 1 (12)0(313020103130201031302010302010ccccyccccyccccycccyxxxx联立上述方程，

39、求解得c10=5, c20=3, c30=4, c31=2。于是，系统的零输入响应为kkkxkky)5 . 0)(24() 3 . 0( 3)2 . 0( 5)(0k 与连续系统中的H(p)一样，H(E)中若有复极点，则必定共轭成对。若设H(E)的共轭复极点为 jjerrer*121)sincos()()()()(2120102010kckcececececkykjkjkkkjkjx式中： )(2010220101ccjcccc例例 5.4-2 设描述离散时间系统的差分方程 为 )2(2) 1()2(25. 0)(kfkfkyky系统初始条件为yx(0)=2, yx(1)=3。试求k0时系统的

40、零输入响应。 解解 写出系统传输算子 25. 0225. 012)(2221EEEEEEH其极点是一对共轭复极点： 22215 . 05 . 0,5 . 05 . 0jjejrejr由式(5.4-22)或式(5.4-23)， 得 2sin2cos)5 . 0()(21kckckykx利用初始条件，得到35 . 02sin2cos5 . 0) 1 (2)0(2211cccycyxx即c1=2,c6=6于是得出系统的零输入响应2sin62cos25 . 0)(kkkykx0k5.5 离散系统的零状态响应离散系统的零状态响应 设系统的初始观察时间为k0，所谓离散时间系统的零状态响应，是指该系统在k0

41、时刻 的状态或者历史输入为零时，仅由kk0时加入的输入所引起的响应，通常记为yf(k)。 在连续系统的时域法分析中，我们根据信号的分解特性和LTI系统的线性时不变特性， 导出了系统零状态响应的计算公式。具体做法包括： (1) 将一般信号分解为众多基本信号单元的线性组合； (2) 求出基本信号激励下系统的零状态响应； (3) 导出一般信号激励下系统零状态响应的计算公式。5.5.1 离散信号的时域分解离散信号的时域分解 根据单位脉冲序列定义和序列位移的概念， 我们有 01)(mkmkmk于是可得 0)()()(mfmkkfmkmk因此，对于任意序列f(k)，可写成 )2()2() 1() 1 ()

42、()0() 1() 1()2()2()(kfkfkfkfkfkf即 )()()(mkmfkfm(5.5-1)图 5.5-1 离散信号的时域分解 1023kf (k)321451 2 3 1 2可以将图 5.5-1 所示的序列分解表示为 dtftf)()()()4()2(2)(3) 1()2(2)(kkkkkkf显然，式(5.5-1)是与连续时间信号f(t)的时域分解公式: 相对应的。在连续系统时域分析中，我们还给出了另一个分解公式 dtftf)()( )(容易得到相应的分解公式为 mmkmfmfkf)()1()(5.5.2 基本信号基本信号(k)激励下的零状态响应激励下的零状态响应 设系统初始

43、观察时刻k0=0，则离散系统对于单位脉冲序列(k)的零状态响应称为系统的 单位脉冲响应，或简称为单位响应， 记作h(k)。 LTI离散系统的单位响应可由系统的传输算子H(E)求出。例例 5.5-1 单极点情况。若系统传输算子 rEEEH)(具有单极点E=r，则相应的差分方程为 )()()(kEfkyrE令f(k)=(k)时，其yf(k)=h(k)， 故有)()()(kEkhrE即 ) 1()() 1(kkrhkh移项后有 ) 1()() 1(kkrhkh根据系统的因果性，当k-1时，有h(k)=0。以此为初始条件， 对式(5.5-6)进 行递推运算得出 krkkrhkhrrhhrrhhrhh)

44、() 1()()2() 1 ()2() 1 ()0() 1 (1)0() 1()0(2因此有 )()()(krkhrEEEHk例例 5.5-2 重极点情况。设系统传输 算子 2)(rEEEH在E=r处有二阶重极点。写出系统的差分方程 )()()(2kEfkyrE同样，令f(k)=(k)，得到单位响应h(k)的求解方程为 )()()(2kEkhrE将该方程改写为 )()()(kEkhrErE可将上式方括号中的(E-r)h(k)表示为 )()()(krkhrErEk或者写成 )()() 1(krkrhkhk采用例 5.5-1 类似求解方法，可求得系统的单位响应 ）（kkrkhk1)(于是有 )()

45、()()(12kkrkhrEEEHk同理，可得 )(! 2) 1()()()(23krkkkhrEEEHk以及d阶重极点相应的单位响应 )()2() 1()!1(1)()()(1krdkkkdkhrEEEHdkd设LTI离散系 统的传输算子为 0111101111011011)(1)(aEaEaEbEbEbEbEEaEaEbEbbEHnnnmmmmmnnnmmm求单位响应h(k)的具体步骤是：第一步， 将H(E)除以E得到 ;EEH）（ 第二步， 将 展开成部分分式和的形式； 第三步， 将上面得到的部分分式展开式两边乘以E， 得到H(E)的部分分式展开式 ;EEH）（gidiigiiirEEK

46、EHEEH11)()(）（第四步，由式(5.5-11)求得各Hi(E)对应的单位响应分量hi(k)；第五步， 求出系统的单位响应 giikhkh1)()(例例 5.5-3 求图 5.5-2 所示离 散系统的单位响应h(k)。 图 5.5-2 例 5.5-3图 f (k)D1D2y(k)y(k1 )y(k2 )解解 )()2(2) 1()(kfkykyky或写为 )()2(2) 1()(kfkykyky相应的传输算子为2211)(2221EEEEEEH将 进行部分分式展开，得EEH)(23/213/1)2)(1()(EEEEEEEH由于 )() 1(31)(1311kkhEEk)()2(32)(

47、2322kkhEEk所以，系统的单位响应 )()2(32) 1(31)()()(21kkhkhkhkk于是232131)(EEEEEH例例 5.5-4 如图 5.5-3 的离散 系统，求其单位响应h(k)。 图 5.5-3 例 5.5-4图 f (k)E11E12y(k)1x(k)E1x(k)E 2x(k)1解解(1) 列算子方程。 )()(2)()(21kfkxEkxEkx它可写为)(211)(21kfEEkx由右端加法器的输出端可列出方程)()1 ()()()(22kxEkxEkxky)(211)(212kfEEEky系统的输入输出算子方程 (2) 求单位响应。 )1)()1 (2211)

48、(20222212EEHEEEEEEEEH)()2(32) 1(31)()(00kkhEHkk)2()2(32) 1(31)()(222020kkhEEHkkkkkkkkkkkkkkh)2(21)2(32) 1(310)2()2(32) 1(31)()2(32) 1(31)(22将上面两个单位响应分量相减，即可得到系统的单位响应 21 , 00kkk例例 5.5-5 设描述离散时间系统的差分 方程为 ) 1(25. 0)2(3) 3(11)(05. 0)2(45. 0)2(2 . 1) 3(kfkfkfkykykyky求系统的单位响应。解解 由已知差分方程得系统传输算子05. 045. 02

49、. 125. 0311)(2323EEEEEEEHEEH)(222232)5 . 0(55 . 0102 . 01)5 . 0)(2 . 0(25. 031105. 045. 02 . 125. 0311)(EEEEEEEEEEEEEEH将 进行部分分式展开，得即2)5 . 0(55 . 0102 . 0)(EEEEEEEH由式（5.5-11）得)()5 . 0(5)5 . 0(5)()5 . 0(105 . 010)(2 . 02 . 012kkEEkEEkEEkkk因此，系统单位响应为)()5 . 0(5)5 . 0(102 . 0)(1kkkhkkk5.5.3 一般信号一般信号f(k)激

50、励下的零状态响应激励下的零状态响应 设离散时间系统的输入为f(k)，对应的零状态响应为yf(k)。由离散时间信号的时 域分解公式(5.5-1)知道，可将任一输入序列f(k)分解表示成众多移位脉冲序列的 线性组合，即 )()()(mkmfkfm 根据LTI离散系统的特性，应用单位响应h(k)可以分别求出每个移位脉冲序列f(m)(k-m)作用于系统的零状态响应。然后， 把它们叠加起来就可以得到系统对输 入f(k)的零状态响应yf(k)。 )()()()()()()()()()()()()()()(khkfkfmkhmfmkmfmkhmfmkmfmkhmkkhkmm单位响应定义 系统的时不变特 性

51、yf(k) 的齐次性 yf(k)的叠加性 信号的分解公式和卷积和运算 定义 于是，得到系统在一般信号f(k)激励下的零状态响应为 mfkhkfmkhmfky)()()()()(5.5 - 18)可将离散时间系统的完全响应表示为 gidjkijijfxikhkfrkckykyky110)()()()()( 这一结果表明：LTI离散时间系统的零状态响应等于输入序列f(k)和单位响应h(k)的卷 积和。 例例 5.5-6 已知离散系统的输入序列f(k)和 单位响应h(k)如下： )(21)()5()()(kkhkkkfk求系统的零状态响应yf(k)。 解解 根据式(5.5 - 18)，有 )()5(

52、)()()()(khkkkhkfkyf由卷积和的分配律，将上式写成 )()5()()()(khkkhkkyf查卷积和计算公式表 5.1，得 )(212211211)(21)()()()(11kkkkhkkykkkf由系统的时不变特性，得 )5(212)5()()5()(512kkykhkkykff于是，系统的零状态响应为 )5(212)(212)()()(521kkkykykykkfff例例 5.5-7 描述某离散系统的差分方程 为y(k)-0.7y(k-1)+0.12y(k-2)=2f(k)-f(k-1)若输入f(k)=(0.2)k(k)，零输入响应初始条件yx(0)=8, yx(1)=3。 试求系统的零输入响应、零状态响应和完全响应。 解解 写出系统的算子方程