《数字信号处理教程——MATLAB释义与实现》第五章-精品.ppt

第五章 变换域中的离散时间系统 5.1 z变换 傅氏变换有两个缺点，其一，在实际中许多有用的信号,如μ(n)和n* μ(n)...等,它们的离散傅利叶变换不存在。其二，系统对初始条件的暂态响应，或由时变输入引起的系统响应，都无法用离散傅利叶变换方法来计算。 为了克服上述的两个缺点，要把离散傅利叶变换方法进行推广，推广后的方法称为z变换。 z变换的定义 定义：对于一个序列x(n)，其双边z变换定义为: 其中 是复变量，其向径为ρ，幅角为φ，构成一个复平面，称为z平面。 序列x(n)的单边z变换则定义为： . z变换的定义 例如，对于序列 它的双边z变换是： 其单边z变换是: 如果序数向量n只取正值，称为右序列，右序列的z变换X(z)中只出现z的负次幂；当采样序数n取负值时，X(z)中将出现z的正次幂。在工程中，人们感兴趣的主要是右序列。 z变换的收敛域 z变换的收敛性1.有限长序列：在除原点外的全z平面上收敛； 2。无限长右序列：在一个半径为r（称为收敛半径）的圆外的全z平面上收敛； 3。无限长左序列：在一个半径为r（也称为收敛半径）的圆内收敛； 4。双向无限序列：右序列和右序列收敛区的（环形）公共区，也可能没有； 分别见下图中子图(a),(b),(c),(d)。 z变换的收敛域 z反变换 单边或双边z变换的反变换定义为。 其中，积分路径是在复数平面中处于收敛域中的一条围线。 在数字信号处理中，不需要用围线积分来求z反变换。5.2节中，将专门讨论求z反变换的其它方法。 z反变换的非单值性 右序列： 和左序列： 两个不同的序列具有相同的z变换 因此必须规定z变换的收敛域，本书限定研究在∞邻域收敛的右序列。故其反变换唯一。 z变换的重要特性 z变换的特性与DTFT和DFT的特性有很多相似之处，其证明都可以类比或从定义直接导出,所以不再重复推证。这里着重讨论几个重要特性的意义和应用。因为讨论都限于右序列，所以也免除了对收敛域的说明。 1.线性特性： 设Z[g(n)]=G(z), Z[h(n)]=H(z), α,β为常数， 则 z变换的重要特性 2。样本的移位： 3。序列卷积： 4。初值定理 z变换的重要特性 5。终值定理 ：若x(n)是因果序列，且其z变换的极点均在单位圆内部，最多只有一个一阶极点在z=1上，则x(n)在n趋于无穷时的终值等于 在信号处理中有用的z变换特性，主要是上面几个，下面的几个特性可用于解析分析，对工程计算用处不大，供查考。 z变换的重要特性 6。乘以指数序列 7。序列乘以n(z变换微分) 8。时域折叠 z变换的重要特性 9。复序列共軛 10。序列乘积： 用z变换计算卷积实例5.1.3 设x1(n) = [2,3,4] ，x2(n) = [3,4,5,6]，求它们的z变换及两者的卷积输出y(n)。 解：由z变换特性3，可得： X1(z)=2+3z-1+4z-2 和 X2(z)=3+4z-1+5z-2+6z-3 为求x1(n)和x2(n)的卷积，先把X1(z)和X2(z)相乘 Y(z)= (2+3z-1+4z-2)(3+4z-1+5z-2+6z-3） = 6+17z-1+34z-2+43z-3+38z-4+24z-5 它的反变换就是卷积输出序列： y(n) ＝ 6 17 34 43 38 24 z变换的计算实例 其实多项式相乘和卷积计算相仿，在MATLAB中用的是同一个函数conv。在本例中要求z-3的系数，可以把第二组系数反过来排列,与第一组系数对齐如下。 2 3 4 6 5 4 3 12+ 15 + 16+ 0= 43 把对应项相乘（空项看作零）并逐项相加，得到上面的结果43。求z的其他幂次 的系数时只需把第二组系数向左或向右移位即可，所以其计算和卷积过程相同。 z变换的计算实例5.1.4 求长度为N的方波函数x(n)=RN(n)的z变换 解： RN(n)是一个因果的有限长序列。当N趋向于无穷时， 。这就是阶跃函数的z变换。这类常用序列的z变换示于表5.1.1中 z变换的计算实例5.1.5 用z变换性质和z变换表求下面序列的z变换。 解：这是一个解析推导题。依次运用z变换的移位特性、微分特性和查表，得到 5.2 z反变换和差分方程的解 由定义（5.1.13）可知，求z反变换需要在一复围线上求积分： 工程上计算z反变换并不需要求积分。解决这个问题有三种方法。（1）极点留数法；（2）部分分式法；（3）幂级数法（长除法）。MATLAB为