DSP第二章Z变换.pptVIP

数字信号处理 II Z变换 机械学院 邹鲲 概述 z变换是离散时间傅立叶变换的推广形式 z变换是数字滤波器设计和分析的重要工具 序列的Z变换 Z变换的收敛域Region of convergence, ROC Z变换的零极点 例题 序列特性对收敛域的影响 Z变换的性质和定理 逆Z变换 用留数定理求逆Z变换 部分分式展开法 幂级数法(长除法) 用Z变换解差分方程 Matlab数值求解 例2.31 求差分方程 输入信号： 初始状态： filtic() filter()见例的matlab程序 系统函数 离散系统的因果性与稳定性 因果稳定系统 因果：Rx-|z|≤∞ 稳定： Rx-1 频率响应的几何确定法 Matlab求解频率响应函数 X(z)是z的有理分式，可分解成部分分式： 对各部分分式求z反变换： X(z)= B(z) A(z) =X1(z) + X2(z) + + Xk(z) … x(n)=IZT [X(z)] = IZT[X1(z)] + IZT[X2(z)] + + IZT[Xk(z)] … 进一步不讲 把X(z)展开成幂级数 级数的系数就是序列x(n) = + x(-1)z1 + x(0)z0 + x(1)z-1 + x(2)z-2 + … … 1、单边Z变换 序列x(n)的单边Z变换定义为： X+(z)= ∑ n=0 ∞ x ( n ) z-n 表示为：Z+[x(n)] 或 x(n) X+(z) z+ 因果序列的单边Z变换与双边Z变换相同 2、单边Z变换的时移性 如果有 x(n) X+(z) z+ 则有 x(n-k) z-k [ X+(z) + ∑ n= -k -1 x ( n ) z-n ] z+ 1) 右移 x(n+k) zk [ X+(z) - ∑ n= 0 k-1 x ( n ) z-n ] z+ 2) 左移 其中， k0 ∑ i= 0 N ai y ( n - i ) ∑ j= 0 M bj x ( n - j ) = 差分方程 求解步骤： 对其进行单边Z变换 变换后的对代数方程求解，得到Y(z) 对Y(z)反变换 零输入响应 零状态响应 激励为因果序列时，零状态响应 例：已知差分方程 y(n) =by(n-1)+x(n)，式中x(n)=anu(n), y(-1)=2，求y(n)。 解：将已知方程进行单边Z变换： Y(z) = bz-1Y(z) + by(-1) + X(z) Y(z) = 2b + X(z) 1 - bz-1 式中 X(z) = 1 1 - az-1 ,|z||a| Y(z) = 2b 1 - bz-1 + 1 (1 - bz-1)(1 - az-1) y(n) = [ 2bn+1 + a - b 1 (an+1 – bn+1) ] u(n) 解析求解 例2.19：已知差分方程 y(n) -5y(n-1)+6(n-2)=0，初始条件y(0)=2，y(1)=3，求零输入响应 例2.19：已知差分方程 y(n) -by(n-1)=x(n)，初始条件y(-1)=0，x(n)=anu(n)，求零状态响应 例题2.22 已知线性时不变系统的系统函数为 求系统的差分公式。 解：展开H(z) 例题2.23 已知线性时不变系统的差分方程 求系统函数及单位响应；求激励为 的零状态响应 解： 系统函数 单位响应 零状态响应 Matlab filter数值求解 离散系统的系统函数和频率响应 LTI系统的系统函数H(z)：单位抽样响应h(n)的z变换 其中：y(n)=x(n) * h(n) Y(z)=X(z)H(z) 系统的频率响应H(ejω ）：单位抽样响应h(n)的FT变换 单位圆上的系统函数 H(ejω) = H(z) = FT[h(n)] z= ejω H(z) = ZT[h(n)] = ∑ n= -∞ ∞ h(n) z-n X(z) Y(z) = 归一化频率概念： LTI系统是因果系统的充要条件：h(n)=0 n0 因果序列 右序列 Z变换H(z)的收敛域(ROC)：圆外域 序列h(n)绝对可和，即 而h(n)的z变换的Roc： ∑ n= -∞ ∞ |h(n) z-n| ∞ 稳定系统的系统函数H(z)的Roc须包含单位圆 ∑ n= -∞ ∞ |h(n)| ∞ 1）因果： Rx-|z|≤∞ 2）稳定： 稳定因果系统：系统函数的极点全部落在单位圆内 * * Z变换的定义 z是一个复变量， 它所在的复平面称为z平面。 单边Z变换 双边Z变换 收敛域(ROC)：对于任意给定序列x(n)，使其z变换X(z)收敛的所有z值的集合称为