离散时间系统变换域分析难点.doc

第八章 离散时间系统的变换域分析8-1 引 言一、 变换域分析的目的：类似于连续时间系统的L.T.，离散时间系统通过Z变换（Z.T.），可以将原来求解差分方程的问题转变为求解代数方程的问题，其目的是通过变换域分析将原来的求解问题简化。 二、 Z变换的发展史 十八世纪，英国数学家棣莫弗（De Moivre）提出生成函数，并应用于概率论。实质上，生成函数与Z变换的形式相同。从十九世纪拉普拉斯（P.S.Laplace）到二十世纪沙尔(H.L.Seal)等人都对其进行了进一步深入研究。二十世纪六十年代起，由于计算机技术和控制技术的飞速发展，抽样控制理论的应用，离散信号处理和数字信号处理得到了广泛应用。作为离散时间系统分析的重要工具，Z.T.得到了很大的发展，其用途甚至超过了L.T.三、 离散时间系统的分析方法1 离散时间系统的Z域分析法，这在本课程进行研究。2 离散时间系统的频域分析法，即利用离散傅里叶变换（DFT）在离散时间系统分析中同样占用很重要的地位，而DFT的快速算法FFT的提出使得DFT在各种信号处理场合得到的广泛的应用。这在数字信号处理课程中进行。3 除了DFT以外，还有如沃尔什变换等分析方法，在离散信号处理中同样得到的很广泛的应用。这在数字信号处理课程中进行。8-2 Z变换定义及其收敛区域一、 Z变换的定义Z变换的定义可以从纯数学的角度进行，也可以通过信号分解的角度提出，后者更加容易理解。本课程中，通过连续时间系统的F.T.导出Z.T.。离散时间信号f(k)可以看成是连续时间信号通过抽样而得到的冲激序列：对其进行F.T.：根据Dirichlet条件，只有在信号满足绝对可积条件的情况下才成立，即满足绝对可和条件：时，FT才存在。如果不满足，可以利用LT中的方法，在信号上首先乘以一个衰减因子，然后再求FT。这样一来上式就可以变成为：令，代入上式，得：上式称为序列f(k)的Z变换。F(z)被称为序列f(k)的生成函数。l 上面的推导反映了抽样信号的FT与用其冲激序列的强度构成的信号序列的ZT之间的关系，即： 而抽样信号的LT与用其冲激序列的强度构成的信号序列的ZT之间的关系为：l 在某些情况下，Z变换的求和限可以简化：1、 如果f(k)是一个左边序列（其在k0时才有非零值），则： 2、 如果f(k)是一个右边序列，则： 3、 如果f(k)是一个有限长序列，则： 二、 单边Z变换与双边Z变换双边Z变换与单边Z变换的区别从应用上考虑，从实际的因果系统和非因果系统上考虑。本课程主要考虑单边Z变换： 三、 Z变换的收敛域 ZT是一个级数求和问题，ZT存在意味着级数收敛。Z变换的收敛域也就是使这个级数收敛的全部Z的集合，即对于任意序列f（k）的z变换，使存在且有限的z值的取值范围称为的收敛区。1、 级数收敛的判别方法：1） 比值法：2） 根值法：其中，为其级数的第k项。2、 几种常见序列的收敛域：1） 右边序列：利用根值法，有：所以，右边信号的收敛域为是半径为R、圆心在原点的圆以外的全部区域。例：单边指数序列的收敛域。解：用上面的结论（根值法）：思考：如果右边序列的起始点不在0,收敛区间应该怎样？提示：收敛域是否包含+？2） 左边序列同上可得左边序列的收敛域为：即左边信号的收敛域为是半径为R、圆心在原点的圆以内的全部区域。例：单边指数序列的收敛域。解：用上面的结论（根值法）：思考：如果左边序列的起始点不在-1,收敛区间应该怎样？提示：收敛域是否包含原点？3） 双边序列与连续时间系统一样，双边序列也可以看成右边序列和左边序列之和，收敛域为两个序列的公共收敛域。收敛域可能存在（当两个序列的收敛域有公共区间时），也可能不存在（当两个序列的收敛域没有公共区间）。如果存在，其收敛域为一个环行区域。 例：求序列的收敛区。解：它的收敛域为左边序列和右边序列的公共收敛区间。1、 当时，两者没有公共收敛区间，Z变换不存在。2、 当时，收敛域为4） 有限长序列：（实际上，收敛域是按照1）-3）的情况确定。）a、 当或,收敛域b、 当,收敛域c、 当或,收敛域d、 当，收敛域四、 常见右边序列的ZT1、 单位函数：，收敛域：全平面。2、 单位阶跃信号：收敛域：3、 单边指数序列： ，收敛域：4、 单边正弦和余弦序列：可以通过上面指数序列推导出，见P386-387其它常见ZT：见P387,表8-1五、 左边和双边序列的ZT计算方法：1、 左边序列ZT求法： 由此可以得到由右边序列计算左边序列ZT计算方法：1） 将序列f(k)反褶，称为右边序列f(-k)；2） 求f(-k)的右边ZT，假设为，收敛域为；3） 得到左边序列的ZT：，收敛域为2、 双边序列ZT求法：与双边信号的LT一样，可以将双边序列分解为左边序列和右边序列之和，分别求解。例：求的ZT解：其中：1），收敛域：2）为了求，a、 将信号反褶，成为新的右边序列：b、 求右边序列ZT：，收敛域：c、 得到原序列ZT：，收敛域：4） 综合得到双边序列的ZT：a、 如果，则f(k)的双边ZT不存在（两个收敛域没有公共部分）；b、 如果，则f(k)的双边ZT为：收敛域：8-3 Z变换的性质1、 线性：2、 移序特性：1） 单边ZT移序特性：a、 增序：可推导出： . b、 减序：推广：l 移序算子q的作用相当于乘z;l 移序计算不影响收敛域；l 移序特性与LT中的微分特性很相似：2） 双边序列移序： ， 3、 （z域）尺度变换特性：若，收敛区域，则：收敛区域。证：4、 （z域）微分特性： 证：即， 。例：求斜变函数的ZT，见P61。5、 卷积定理：单双边相同。6、 初值和终值定理：在f(0)存在的条件下，（如果和存在，且的LT也存在，则：）在f()存在并且有限的条件下，（如果和存在，的LT也存在，且F(s)的极点位于s平面的左半平面，在s=0上至多存在单极点，则：）习题：. 8.2(2); .8.3(1)、(6)；.8.5.*8-4 反Z变换反Z变换有三种方法：1） 级数展开法；2） 部分分式展开法；3） 留数法。一、 级数展开法：指导思想：将F（z）表示成Z变换的原始形式，将各个元素与f(k)对号入座。与z变换的定义式：比较，即得到：实现途径：长除。例：求的原函数。见P397l 用这种方法容易求得信号的前面的几个点上的值，但是无法得到解析表达式。l 用这种方法可能得到多个解。(在升幂和降幂问题的处理上，见P397-398)l 这种方法无法与收敛域相结合，得到正确的原函数。 二、部分分式展开法：同LT中的Heaviside分解法。1、 其用到的基本变换为：2、 对进行部分分式展开，对应于上面的基本的ZT公式，就可以得到原函数。3、 也可以采用另外一种基本函数：这时候只要对F（z）进行部分分式展开即可。4、 上面讨论的是单边ZT的反变换。与LT一样，在双边ZT中，F(z)的原函数与其收敛区间有关。可以是右边序列的象函数，也可以是左边序列的象函数，差别在于收敛域不同。所以，必须根据收敛域，决定部分分式展开式中各项的归属。例1：同上。P399例2：见P399,8.10解:,解得:A=2/3，B=1/3，C=-2，所以：当时，右边函数，有：，当时，左边函数，有：，当时，双边函数，有：三、留数法：1、 通过计算留数，可以得到原函数：很多教材上将其作为反Z变换的定义：（这样注意围线C的定义）证明从略。例：P400-401小结：l 留数法不仅可以用于计算单边反Z变换，而且可以用于计算双边反Z变换。l 用留数法进行计算，可能会遇到计算z=0点的各阶留数的计算，不很方便。（主要由于公式中（当k-1理想抽样序列离散序列 F(s) F(z)已知信号的F(s),通过可以得到：对f(t)理想抽样，其冲激幅度序列为：对序列求ZT：l 假设：，则：可见：F(s)在处有极点，而F(z)在处有极点。l 假设（假设没有重极点），则有： 在F(s)没有重根的情况下,可以通过部分分解的方法得到F(z)。l 从上面可以看出：F(s)的极点和F(z)的极点之间的关系为：，或 或：，可见，F(s) 和F(z)的极点的映射关系与上面的关系相同。这里同样有多点映射的问题。8-6 离散时间系统ZT分析法时域差分方程Z域代数方程系统响应时域解。将系统响应依然分成和两部分讨论。 一、的ZT求解法：在输入信号为零的条件下，差分方程变为了一个齐次差分方程。其一般形式为：对其求ZT，可以得到：（移位特性）所以，有了初始条件，就可以通过直接写出，再由反ZT就可以得到。但是，这种方法比时域解法复杂，因为：1、 形式复杂，难于记忆；2、 要进行反ZT计算。二、的ZT求解法：零状态响应有很多推导方法。教材上提出了直接用这个差分方程的求解方法。这里给出一个更加简单的方法。为此将差分方程改写为： 然后对方程两边求ZT（注意：1、系统初始状态为零;2、同时激励信号也是一个有始信号;3、对于因果系统，m=n）：定义：则：其中，H(z)称为系统的转移函数。它可以根据系统的差分方程的系数直接写出。可见，离散时间系统的零状态响应的求解方法与连续时间系统中的方法完全一样，只不过LT变成了ZT而已。这里，同样可以证明，H(z)是系统的单位函数响应的ZT，即：三、系统的全响应求解：1、通过和分别求。综合上面的结果，可以得到：在实际应用中一般不直接使用这个公式，而是利用其原理进行计算。2、直接求解法：对差分方程两边直接求ZT，并代入初始条件。原始差分方程：或：两边同时求ZT：注意这个公式与前面公式的差别：1、响应初始条件的含义；2、是否考虑激励的初始值。例1：P411例2：P430，8.18（1）解:做Z变换:，(*推导见后)，解得：A=1，B=0.9，得：，得到：。推导:四、 离散时间系统的稳定性l 与连续时间系统中的结论相似，离散时间系统稳定的充分必要条件是其单位函数响应绝对可和。l 这要求系统传输函数的极点都在单位圆的内部。l 如果系统在单位圆上有单极点，则系统是临界稳定的。l 如何根据H(z)的分母D(z)的系数判断系统是否稳定？直接根据D(z)的系数不能判断其根是否处于单位圆以内