1、2.2 离散系统的时域分析离散系统分析与连续系统分析在于多方面是相互平行的，有许多类似之处。连续系统可用微分方程描述，离散系统可用差分方程描述。差分方程和微分方程的许多求解方法在许多方面是相互对应的。在连续系统分析中，卷积积分具有重要的意义；在离散系统中，卷积和也具有同等重要的地位。连续系统分析与离散系统分析的相似性为读者学习本章节提供了有力条件，不过，读者应该十分注意他们之间存在着的重要差异。在离散系统分析中，激励（输入）用 表示，响应（输出）用 表示，其中()fn()yn为整数；初始状态用 表示，其中 为正常数，通常取 。n0()xn00下面，从离散系统的差分方程（或系统框图）及其求解开始
2、，研究 离散系统的时LTI域分析。2.2.1 LTI 离散系统的响应差分与差分方程与连续时间信号的微分与积分运算相对应，离散时间信号有差分及序列求和运算。设有序列 ，则称 ， ， ， 等为 的位移()fn(2)fn (1)f()fn (2)f ()fn序列。序列的差分可分为前向差分和后向差分。一阶前向差分定义为（2.2-1）()1)(deff一阶后向差分定义为（2.2-2）()()defnfn式中 和 称为差分算子。由式（2.2-1）和（2.2-2 ）可见，前向差分与后向差分的关系为（2.2-3）1ff二者仅位移不同，没有原则的差别，因而它们的性质也相同。本书主要采用后向差分，并简称其为差分。
3、由差分的定义，若有序列 、 和常数 、 ，则1()fn2f1a212 2()()()()()afnfafnf122ff（2.1-4）12()()ff这表明差分运算具有线性性质。二阶差分可定义为 2(n)()()1)()1)def fnfnf（2.2-5）()21)(2)fnffn类似的，可定义三阶、四阶、五阶差分。一般的， 阶差分k（2.2-6）k10()()(1)()kdefkjjffnj式中（2.2-7）!,2.()kjkjj为二项式系数。序列 的求和运算为()fn（2.2-8）()nif差分方程是包含关于变量 的未知序列 及各阶差分方程的方程式，它的一般形式()fn可写为（2.2-9a）
4、,(),(),()0r NFnffnyy 式中差分的最高阶为 阶，称为 阶差分方程。由式（2.2-6）可知，各阶差分均可写成及其各位移序列的线性组合，故上式常写为()y（2.2-9b ）,()1,(),(1),()0ffrn 通常所说的差分方程是指式（2.2-9b）形式的方程。差分方程是具有递推关系的代数方程，若已知初始条件和激励，利用迭代法可求得差分方程的数值解。例 2.2-1 若描述某离散系统的差分方程为 ()31)2()(ynynf已知初始条件 激励 ，求 。0,f)yn解：将差分方程中除 以外的各项都移到等号右端，得()312()(ynynf对于 将已知初始值 代入上式，得2,0),(
5、)(f类似地，一次迭代可得 3)213)0yy(4)(4f由上例可见，用迭代法求解差分方程思路清楚，便于编写计算机程序，能得到方程的数值解，但它常常不宜得出解析形式（或称闭式）的解。2.2.2 差分方程的经典解一般而言，如果单输入-单输出的 系统的激励为 ，其全响应为 ，那么，LTI()fn()yn描述该系统激励 与响应 之间的关系的数学模型是 阶常系数线性差分方程，它()fn()y可以写为 10()()()NyaanN（2.2-10a）MbfnfbfM式中 ， 都是常数。上式可缩写为(0,2,)i (,1)j(式中 ) （2.2-10b）N00i ji jayfnN1a与微分方程的经典解相类
6、似，上述差分方程的解由齐次节和特解两部分组成。齐次解用 表示，特解用 表示。即()hn()p（2.2-11）hyny齐次解当式（2.2-10）中的 及其各位移项均为零时，齐次方程()f（2.2-12）10()()nyaaynN（2.2-12）的解称为齐次解。首先分析最简单的一阶差分方程，若一阶差分方程的齐次方程为（2.2-13）()它可改写为 (1)yna与 之比等于 表明，序列 是一个公比为 的等比序列，因此()()(n)y()a应有如下形式（2.2-14 ）nyCa式中 为常数，由初始条件确定。C对于 阶齐次差分方程，它的齐次解由形式为 的序列组合而成，将 代入到时n nnC（2.2-12
7、） ，得 1110nnNnNaa由于 ，消去 ；且 ，以 除上式，得00（2.2-15）110Nnaa上式称为差分方程（2.1-10）和（ 2.1-12）的特征方程，它有 个根 ，称为n(1,2)in差分方程的特征根。显然，形式为 的序列都满足式（2.2-12） ，因而它们是式（2.1-Cni10）方程的解。依特征根的不同取值，差分方程齐次解的形式如表 2.2-1，其中 、 、CiDi、 等为待定常数。Aii表 2.2-1 不同特征根所对应的齐次解特征根 齐次解单实根 nC重实根r1210CCrnrnkk一对共轭复根1,2jajbe 或 ，其中cos()i()nDcos()Ajej特解特解的函
8、数形式与激励的函数形式有关，表 2.2-2 列出了几种典型的激励 所对()fn应的特解的 。选定特解后代入原差分方程，求出其待定系数 （或 、 ）等，(n)py iPA就得出方程的特解。表 2.2-2 不同激励所对应的齐次解激励 (n)f 齐次解m所有特征根均不等于 1 时；11mPnP当有 重等于 1 的特征根10nrm r时。na， 当 不等于特征根时；na当 是特征根时；10mnPa， 当 是 重特征根110rrnPa r时。或cos()ni()当所有的特征根均不等于Pcos()in()Qje或 ，其中A=P+jjAeQ全解式（2.2-10）的线性差分方程的完全解是其齐次解与特解之和。如
9、果方程的特征根均为单根，则差分方程的全解为（2.2-16 ）1()()()NnhpipynyCy如果特征根 为 重根，而其余 个特征根为单根时，差分方程的全解为1rr（2.2-17 ）11() ()Nrinnjpi jryny式中系数 、 由初始条件确定。iCj如果激励信号是在 时接入的，差分方程的解适用于 。对于 阶差分方程，0n0nN用给定的 初始条件 就可确定全部待定系数 和 。如果差分方n(),1,()yyN iCj程的特征根都是单根，则方程的全解为式（2.2-16 ） ，将给定的初始条件分别代入到式（2.2-16） ，可得(0),1.,yN（2.2-18）121112(0)(0)()
10、 ()npNNnpyCyC 由以上方程可求得全部待定系数 。(,)i例 2.2-2 若描述某系统的差分方程为（2.2-19）()41)2)nynyf已知初始条件 ， ；激励 。求方程的全解。0()2,0n解：首先求齐次解。上述差分方程的特征方程为 24可解得特征根 ，为二重根，由表 2.2-1 可知，其齐次解122()()()nnhyC其次求特解。由表 2.2-2，根据 的形式可知特解f(),0npyP将 ， 和 代入到式（2.2-19 ） ，得n1(2)py1224()nnnPf上式中消去 ，可解得 ，于是得特解(),04npy差分方程的全解 121,04nnnhpyC将已知的初始条件代入上
11、式，有 21(0)04y12.1由上式可求得 ， ，最后得方程的全解为C（2.2-20）()()(2n04nny自 由 响 应 强 迫 响 应 ，差分方程的齐次解也称为系统的自由响应，特解也称为强迫响应。本例中由于 ，故1其自由响应随着 的增大而增大。n例 2.2-3 若描述某离散系统的差分方程为（2.2-21 ）6()51)(2)(yynf已知初始条件 ， ，激励为有始的周期序列 ，0()10cos02nf,求其全解。解：首先求齐次解。差分方程的特征方程为 26510其特征根为 ， ，方程的齐次解231nnhyC其次求特解。由表 2.2-2 可知，特解cossin2pynPQ其移位序列 11
12、1cssicossin22py PQ(2)(2)(2)cossincossin2pnynPQPQ将 ， ， 代入到式（2.2-21）并稍加整理，得p1p()py n(65)cos65)sin()10cos222nQf由于上式对于任何 成立，因为等号两端的正、余弦序列的系数应相等，于是有01P由上式可解得 ，于是特解（2.2-cosin2cos,n024pyn22）方程的全解 12cosin,032nnphynyC将已知的初始条件代入上式，有 120013y由上式可解得 ， ，最后解得全解12C(2.2-23)cosin32204ny强 迫 响 应自 由 响 应 稳 态 响 应瞬 态 响 应 ，
13、由上式可见，由于本例中特征根 ，因而其自由响应是衰减的。一般而言，如果差分1,2方程所有特征根均满足 那么自由响应将随着 的增大而逐渐衰减趋近in n于零。这样系统称为稳定的系统，这时自由响应也称为瞬态响应。稳定系统在阶跃序列或有始周期序列作用下，其强迫响应也称为稳态响应。2.2.3 零输入响应和零状态响应系统的全响应 也可分为零输入响应和零状态响应。零输入响应是激励为零LTIny时仅由初始状态所引起的响应，用 表示；零状态响应是系统的初始状态为零时，仅x由输入信号 所引起的响应，用 表示。这样 系统的全响应将是零输入响应f fLTI与零状态响应之和，即（2.2-24）xfynyn在零输入条件
14、下，式（2.2-10 ）等号右端为零，化为齐次方程。若其特征根均为单根，则其零状态响应（2.2-25）1NnxxikC式中 为待定常数。i若系统的初始状态为零，这时方程（2.2-10 ）仍是非齐次方程，若其特征根均为单根，则其零状态响应（2.2-26 ）1Nnpiffiyy式中 为待定常数。fiC系统的全响应可分为自由响应和强迫响应，也可分为零输入响应和零状态响应，它们的关系是（2.2-27）111nNnnipxiipfiiyyCy自 由 响 应 零 输 入 响 应 零 状 态 响 应强 迫 响 应式中（2.2-28）11=NnnixiifC可见，两种分解方式有明显的区别。虽然自由响应与零输入
15、响应都是齐次解的形式，但它们的系数并不相同， 仅由系统的初始状态所决定，而 是由初始状态和激励共同xi iC决定。在用经典法分别求解系统的零输入响应和零状态响应时，需要已知各响应的初始值，用以分别确定常数 和 。由式（2.2-25 ）可知，在时刻xiCfi nj（2.2-29）0,12,fyjyjN这时各初始值 中不仅包含有零输入响应的初始值 ，也包含零状态响应的初始xyj值 ，这常常不便区分。为此，如果激励是在 时接入的，通常以fyj 0n描述系统的初始状态。因为在 时，激励尚未接入，显然1,2,yN 零状态响应在这些时刻的值为零，即（2.2-30）0fffy由式（3.1-25）可知，这时（
16、2.2-31）1,2,xxxyyN它们给出了该系统已往历史的全部信息，称其为系统的初始状态。这样，在求解时可分别根据 和 ，利用该系统所满足的差分方程，用迭代,xxyN 1,ff法分别求得零输入响应和零状态响应的初始值 和 ，并xyj 0,12xjN进一步确定系数 和 。xiCfi例 2.2-4 若描述某离散系统的差分方程为（2.2-32）312nynyf已知激励 ，初始状态 ，求系统的零输入响应、,0nf110,2y零状态响应和全响应。解：（1 ） 零输入响应根据定义，零输入响应满足方程（2.2-33）3120xxxynyn由式（2.2-32） ，其初始状态 。11,2xy首先求出初始值 ，
17、式（2.2-33 ）可写为0,xy312xxynn令 ，并将 代入，得0,x，1xxy1323xy式（2.2-34）的特征值为 ，由表 2.2-1，其齐次解12,（2.2-34）12nnxxxyC将初始值代入得 120xx3y可解得 ，于是得该系统的零输入响应12,xxC（2.2-35 ）,0nny实际上，式（2.2-34）满足齐次方程式（ 2.2-33） ，而初始值 ， 也是有该方程0xy1x递推出的，因而直接用 确定待定常数 将更加简便。即在式12xxy， 12C，（2.2-34）中令 ，有,n1204xxyC可解得 ，与前述结果相同。12,xx（2 ） 零状态响应根据定义，零状态响应满足方程（2.2-36 ）32fffynynf和初始状态 。10ff首先求出初始值 ，将式（2.2-37 ）写为,1ffy32fffynnf令 ，并代入 和 ，得0,10ffy,1f（2.2-37 ）3021fffy系统的零状态响应是非齐次差分方程式（2.2-36 ）的全解，分别求出方程的齐次解和特解，得 n121213nnnnfffpffyCyC将式（2.2-37）的初始值代入上式，有
