西南交通大学信息科学与技术学院数字信号处理课件 第六章.pptVIP

第六章 Z变换 z-transform Z变换（z-transform） DTFT提供了离散时间信号与线性时不变系统在频域的表示方法。 由于收敛条件，在许多情况下一个序列的DTFT可能不存在。 因此，在很多情况下它们的频域特征是不能使用的。 6.2 有理Z变换 因此 6.4.4使用MATLAB进行部分分式展开 [r,p,k]= residuez(num,den)用分子多项式的系数向量num和分母多项式的系数向量den 发展了有理z变换的部分分式展开法 向量 r表示留数 向量p表示极点 向量k 表示常数——假分式化成真分式所得多项式的系数（降幂排列）。 6.4.4使用MATLAB进行部分分式展开 语句[num,den]=residuez(r,p,k) 用来实现一个z变换的由部分分式形式到有理形式。 参看程序6_4 对于因果序列{g[n]} ，z变换G(z)可以展开成 z-1 的幂级数 在这个级数展开式中, 乘以 z-n 项的系数就是第n个样本 g[n]。 对于有理G(Z),常用的方法是用幂级数展开的方法来以z-1 的形式表示分子和分母，然后通过长除法获得幂级数展开。 用长除法计算逆z变换 例-考虑 6.4.6 使用MATLAB计算逆z变换 函数impz 能计算出有理z变换G(z)的逆变换。 该函数计算出G(z) 的幂级数展开式的系数向量。 系数的个数可以是用户决定也可以自动决定。 表6.2: z变换的性质 记x[n]=?n?[n] 以及 y[n]=-?n ?[-n-1] 使X(z) 和Y(z) 分别是x[n]和y[n]的z变换。 V(z)的收敛域是 |z||?| 和|z||?|区域的重叠部分。 若|?| |?|,则存在重叠部分，收敛域是环域 |?||z||?| 若|?||?| , 则不存在重叠，V(z) 不存在。 例- 确定因果序列的z变换和收敛域。 其传输函数由上式两边方程式进行z变换得到。 因此 或者, 还可写为 ?1, ?2,…, ?M 是H(z)的有限个零点 ， ?1, ? 2,…, ? N 是H(z)的有限个极点 。 若 N M, 在z = 0有 (N-M) 个零点 。 若 N M, 在z = 0有(M-N) 个极点。 对于因果 IIR数字滤波器，其冲激响应是因果序列 因果传输函数的收敛域 例- 考虑M点移动平均FIR 滤波器，该滤波器的冲激响应为 例 –因果 LTI IIR数字滤波器由一个常系数差分方程描述 ： y[n]=x[n-1]-1.2x[n-2]+x[n-3]+1.3y[n-1] -1.04y[n-2]+0.222y[n-3] 它的传输函数为： 变形为: 若传输函数 H(z)的收敛域包含单位圆，则该LTI数字滤波器的频率响应 H(ej?) 可以简单地在单位圆上计算它的传输函数H(Z)得到即: 对于稳定的有理传输函数其形式可写为 从频率响应的因式形式我们可方便地看出传输函数H(Z)的零点因式 (z-?k)和极点因式 (z-?k) 对整个频率响应的贡献。 幅度函数的表达式为 频率响应的几何形式： 频率响应相当于得到频率响应函数的缩略图。 频率响应的一个典型形式为 当ω 取一个特定值时，相位响应是通过把项p0/d0的相位、线性相位 ω(N-M) 的相位加上所有零点向量的角度，再减去所有极点向量的角度得到的。 因此，LTI数字滤波器传输函数的幅度和相位响应的近似图可以通过研究它的极点和零点位置得到。 当w = f是零点（极点）时，该零点向量有最小（最大)幅度。 为了在特定频率范围内很大程度地衰减信号，需要在该频率范围内将传输函数的零点放在距单位圆非常接近的位置或直接放在单位圆上。 类似地，若需要在某个频率范围内很大程度地强化信号，就需要在该频率范围内将传输函数的极点放到与单位圆非常接近的位置。 作业 阅读课本 p.242 到 275 作业 6.1, 6.5, 6.8, 6.10, 6.14, 6.16, 6.24, 6.37, 6.40, 6.44, 6.45, 6.58 M6.1, M6.3 6.7 传输函数 FIR数字滤波器 冲激响应h[n]定义区间 收敛域 在经过离原点最远的极点的圆的外部 即收敛域为 6.7 传输函数 它的传输函数为 6.7 传输函数 M = 8 6.7 传输函数 传输函数在单位圆上有M个零点，其位置为 Mz=ej2?k/M, 0 ? k ? M-1。 在 z = 0 处有一个(M-1)阶极点，在 z = 1处有一个单极点。 但在z = 1 处的极点恰好抵消 该处的一个零点。 收敛域是除了原点z=0