《第二章反z变换ppt课件.ppt》由会员分享，可在线阅读，更多相关《第二章反z变换ppt课件.ppt（88页珍藏版）》请在三一办公上搜索。
1、2.7.3 Z反变换,即由Z变换式X(z)求相应的序列x(n)， 常用Z-1x(z)表示，逆z变换是一个对X(z)zn-1进行的围线积分，积分路径C是一条在X(z)收敛环域（Rx-，Rx+）以内反时针方向绕原点一周的单围线。,求解反z变换的常用方法有常见的方法有幂级数法部分分式法留数法,1.幂级数法,由上已知，Z变换是一个幂级数表示式,那么，求X(z)的反变换只要将其展开为幂级数形式，再与上式相比较，其系数便是所求的序列x(n)。 幂级数的展开形式还必须依据收敛域,例2.5：已知X(z)=e1/z,|z|0,求其反Z变换,所以得,解：将其展开为幂级数形式,解：,再利用Z变换的线性和位移特性,于
2、是可知X(z)所对应的序列为,一般地，Z变换式为有理分式的情形，可以利用长除法来得到其幂级数展开式,a.先按降幂排列，同上。,b. 先按升幂排列,利用多项式除法得,2. 部分分式法,设X（z）可以分解成,其中,是简单的分式，可以通过Z变换表,查得其对应的反Z变换.根据Z 变换的线性,得到对应的序列,例如: 1. 如可将X(z)表示成,2. 如可将X(z)表示为,于是，,例2.7 已知,将其化为部分分式之和,解：,于是，,3. 留数法,由Z变换X(z)求其相应的序列x(n)，有下面的Z反变换关系式：,积分路径C是一条在X(z)收敛环域（Rx-，Rx+）内的一条包围原点的闭合曲线。,上式中的积分，
3、应用留数定理来求,若X(z)zn-1是z的有理函数，设z0是它的一个s阶极点，可以将X(z)zn-1表示为,则 在z=z0解析， X(z)zn-1在z=z0的留数为,若z0是一阶极点，即s=1，则此留数为,留数法思路,由于积分围线c在X(z)的收敛域内，所以首先要确定收敛域，如未给定要根据其极点确定。考虑被积函数X(z)zn-1在c内或c外的极点的情况确定用(a)式或(b)式来计算。原则是选择X(z)zn-1有有限个极点且极点阶次有限的区域来求留数，而尽量避免求z=的留数。如收敛域在一圆外，应计算n0时的x(n)，选择（a）式，因为此时在c内有有限个阶次有限的极点，zn-1在z=0处解析； 如
4、收敛域在一圆内，应计算n0的x(n)，而利用c外的极点求得n0的情况同样处理,例2.17 用留数法求 的 z 反变换。 解： ，显然， 并且 m=1。X(z)有两个极点：z1=1 及 z2 = 1/3，故有三种可能的收敛域。,(1) 收敛域|z|1: 此时收敛域在|z|=1的园外，围线c之内包含X0(z)的两个极点，所以有：当n1-m=0时，x(n)=0；而当n1-m=0时，有：,(2）收敛域|z|1/3: 此时收敛域在|z|=1/3的园内，围线c之外包含X0(z)的两个极点，所以有： 当n1=m=0时，x(n)=0；而当n1-m=0时，有：,(3) 收敛域1/3 z 1: 此时收敛域在一个环
5、内，X0(z)的极点 z1=1在围线c之外，而 z2=1/3在围线c之内, 于是有：当n1-m=0时， 当n1-m=0时， 因此，当收敛域在环内，,例,则被积函数可写为,n0时：X1在c内和c上解析，故Z反变换积分为0，即有x(n)=0;n0时： X1在c外只有一个一阶极点z=2。此时,故有,N=0时：,在c内只有一个一阶极点z=0，此时,所以,综上可知,求下列Z反变换,例：设：,解：|z|3，所以x(n)为右边序列，X(z)应按z的降幂排列。,z2-6z+9,3z,3z-1,3z-18+27z-1,18-27z-1,+18z-2,18-108z-1+162z-2,81z-1-162z-2,+
6、81z-3,81z-1-486z-2+729z-3,+,得：X(z)=3z-1+18z-2+81z-3+,=3z-1+232z-2+333z-3+, x(n)=n3nu(n-1),例：设：,解：,收敛域|z|2,例：求z的反变换，设：,解：, 当 n0 时，围线内有一个极点：z=1/2, 当 n=0 时，围线内有两个极点：z=0和z=1/2, 当 n0 时， X(z) zn-1 的分母多项式z的阶次比分子多项 式的阶次高二阶或二阶以上，故可用式来计算x(n)。由 于此时在围线外无极点，所以：x(n)=0。,综合、，得到x(n)：,2.6 单边Z变换,2.6.1单边Z变换的定义2.6.2单边Z反
7、变换 其计算方法与双边Z反变换相同，但是其解并不一定是唯一的，因为单边Z反变换算出的只是序列n0时的表达式，n0部分未定。,可以看出，与双边变换的不同就在于对于n的求和范围不同，其幂级数中只含z的负指数项，单边Z变换的收敛域在半径为某个值的圆周之外，包括z=在内。如果序列在n0的定义相同，则它们的单边Z变换就相同，尽管n0的定义可能不同。,2.6.3单边Z变换的性质,前节中的双边Z变换的所有性质除与序列位移有关的性质之外都适用于单边Z变换。单边Z变换的位移特性 设n0为正整数， 右移时：,-1-n0的序列移入正向区间，使Z变换的结果增加了n0项,左移时：,设n2=0,1,n0-1时，,左移时，
8、虽然表达式与双边变换相同，但原正向区间的前n0个值不能参加级数运算，需假设这n0个值都为零，才能与原Z变换序列建立联系,单边Z变换与双边Z变换的差异,单边Z变换适用与需要初值条件解决因果系统的响应问题，可以从某个需要的时刻开始，初始条件体现了在n0的情形对n0以后的反效作用。在无法提供初始条件的场合，如噪声激励，或只需了解稳定状态的场合如滤波器的设计，则使用双边Z变换。,用单边Z变换解线性差分方程,线性差分方程描述的因果系统，当给定适当的初始条件时，需用单边Z变换求解，以得到系统是输出响应。,2.7.4 Z变换与傅氏变换的关系,序列x(n)的Z变换为,令复变量,表示z在单位圆上取值，因此，傅氏
9、变换就是单位圆上的Z变换。,同理，Z反变换中,设单位圆在X（z)的收敛域内，因此可将单位圆作为积分围线c，即有,此即离散信号x(n)的傅氏变换的反变换。,Z变换和傅氏变换实际上都是级数求和，因此都存在收敛问题，这就对序列x(n)有一定的要求：,傅氏变换：,Z变换：,同理,可见，傅氏变换的收敛对x(n)的要求强于Z变换收敛对x(n)的要求，如果傅氏变换不收敛，还可以找到适当的模值r使得其Z变换收敛。,综上可知,模拟信号及其抽样信号在时域和频域中的相互关系, ,2.8系统的差分方程描述,1.非递归型 指输出对输入无反馈，输出值仅取决于系统，在n时刻的输出值可以表示为,对线性非移变系统,若此系统是因
10、果的,若又有iN,N阶线性差分方程,输出对输入有反馈，输出值不仅取决于输入值而且与输出值有关。对于线性、非移变、因果系统，有,2.递归型,可见，当bi=0时，递归型变为非递归系统。非递归是递归的特例,2.8.2系统函数 一个线性、非移变、因果系统的差分方程为,对上式两边作双边Z变换，得,定义该系统的传递函数，即系统函数为,而我们已经知道,由序列卷积的Z变换特性有,由（3）式和（5）式比较可知，系统函数H(z)实际上就是系统的单位取样响应h(n)的Z变换。,系统函数H(z) 表征系统在稳态时的特性，Y(z)经反变换后的输出也是系统在稳态的输出。,当需要求解在一定的初始条件下因果系统的响应问题时，
11、要采用单边Z变换。则求解（1）式两边的单边Z变换得,故可求得,由此可见，输出响应的单边Z变换，不仅与系统的传递函数有关，而且与输入输出的初始状态有关，只有当初始值为零时才与双边Z变换的结果相同。,如果令,此即系统的频率响应,一线性时不变系统y(n)=x(n)-1/2y(n-1)求单位抽样响应,2.8.3 系统函数的零、极点,前面得到,可见，H(z)是两个Z-1的多项式之比，于是可以对其分子分母因式分解，得,这里假设,可见，ci就是系统函数的零点，而dj就是其极点。也就是说系统函数可用其零点和极点来表示。,系统的频率响应，上式中令,将向量Ci和Dj用极坐标表示：,而且,因此,其中，,系统的幅频响
12、应,系统的相频响应,在系统函数的零、极点已知的情况下，可以利用几何作图法求系统的频率响应,可以看出零点、极点在z平面上的位置关系，,在单位圆上的位置移动而改变：在零点附近，幅频响应将出现谷值（最小值），特别地，当零点位于单位圆上时，谷值为零；在极点附近，幅值将出现峰值（最大值），若极点位于单位圆上，峰值为无穷大。,2.8.4 线性非移变因果系统的稳定性,这里来讨论系统函数的极点位置与系统的稳定性的关系。对于线性非移变系统，只需考察,就可以断定该系统是否稳定。,设一N阶线性非移变系统，系统函数为H(z)，由于其单位取样响应是因果的，H(z)的收敛域在一半径为R-的圆外。为方便讨论，假设H(z)只
13、有一阶极点，用pi表示(i=1,2,N),(1)设R-1，即所有极点都在单位圆内。,利用Z反变换的留数法求解可知,这说明，极点都在单位圆内是系统稳定的充分条件。,（2）假设有一个极点pk在单位圆外。,同理，可求得,由于,所以有,可见，单位圆外极点的存在使得系统不稳定,综上所见，一个线性时不变的因果系统稳定的充分必要条件是系统函数（传递函数）H(z)的所有极点都在z平面的单位圆内。显然，对于因果线性时不变系统，单位圆必定在其系统H(z)的收敛域内。,2.9 Matlab方法 2.9.1 常用序列及序列运算的Matlab实现1 单位抽样序列函数 zeros(1,N) 可以产生一个包含 N 个零的行
14、向量，在给定的区间上，可以用这个函数来产生。这个函数的输入参数应该满足条件。,2单位阶跃序列函数 ones(1,N) 产生一个由 N 个 1 组成的行向量,在给定的区间上，可以用它来产生。这个函数的输入参数应该满足条件 。3矩形序列其 Matlab 实现为： rect = zeros(1,N),ones(1,M),zeros(1,P),4实指数序列符号“.”用来实现一个实指数序列。 例2.18 用Matlab实现 。 n =0:10; x = (0.5).n; stem(n,x);,图2.26 例2.18的图形,5正弦序列函数sin（或cos）产生正（余）弦序列。 例2.19 用Matlab
15、实现 x(n)=2sin(0.6n) + 3cos(0.3n+ /3)， 0n10。 n = 0:0.1:10; x = 2*sin(0.6*pi*n) + 3*cos(0.3*pi*n+ pi/3); plot(n,x);,图2.27 例2.19的图形,6序列的翻褶y(n)=x(-n)的Matlab实现为： y=fliplr (x); n=-fliplr (n);,7信号的能量： 在Matlab中采用函数conj来求一个复数的共轭，而离散序列的能量的 Matlab 实现可以采用下述任一种方法。（1）E=sum(x.*conj(x);（2）E=sum(abs(x)，.2);,例2.20 用Ma
16、tlab实现下列序列，并画出相应图形。 解： n=0:10; x=n*Unitstepseq(0,0,10)+3*(0.5).(3*n); stem(n,x); xlabel(n); ylabel(x(n);,图2.28 例2.20 的图形,8序列的离散线性卷积计算Matlab 中计算两个有限长序列的线性卷积的函数是conv ，该函数假设两个序列都是从n=0开始的，其调用格式如下： y = conv(x,h),例2.21 求以下两个序列的线性卷积。x(n)=11,6,3,6,-9，-3n1h(n)=8,17,3,20,9,14，-1n4。 x=11,6,3,6,-9; h=8,17,3,20,
17、9,14; y=conv(x,h),于是用Matlab求得： y = 88 235 159 337 258 133 204 -84 3 -126y(n)定义的区间可以这样求出： 因为 ，其中x(k)的非零区间为,-3k1而h(n-k)的非零区间为 -1n-k4将这两个不等式相加就得到y(n)的非零区间： -4n5。,2.9.2 离散信号变换的Matlab 实现1 离散信号的DTFTDTFT就是2.6节讨论的离散信号的傅里叶变换。在Matlab中，可以利用freqz 函数计算序列的DTFT在给定的离散频率点上的抽样值。假设 可以表示为 ，则freqz函数有如下几种调用方式：,（1）H,w = f
18、reqz(b,a,N)其中，b和a 分别表示X(ej)的分子和分母多项式的系数向量。此函数在单位园上半部上等间隔地计算N个点处的频率响应，返回该系统的N点频率响应矢量 w 和 N 点复数频率响应矢量 H。如果 N 没有说明，则缺省值为 512。,（2）H = freqz(b,a,w)它返回矢量 w 指定的那些频率点上的频率响应，频率范围在0到之间。（3）H = freqz (b,a,F,Fs) 给定单位为 Hz 的抽样频率 Fs，返回矢量 F 指定的那些频率点上的复数频率响应，单位也是 Hz。,（4）H,w = freqz(b,a,N,whole) 在整个单位园上等间隔地计算N 点频率响应，即
19、频率的范围是02。（5）H,F = freqz (b,a,N,Fs) 和 H,F = freqz (b,a,N,whole,Fs),给定抽样频率 Fs，单位为 Hz；返回单位为 Hz 的频率矢量 F。也可以利用Matlab提供的函数abs、angle、real、image等来计算DTFT的幅度（ ）、相位 以及实部和虚部。,例 2.22 已知因果系统 ，试画出 的幅度响应 和相位响应 。,图2.29 例2.22的系统的频率响应,1 z变换与z反变换 （1） 函数tf2zp和zp2tf 函数tf2zp和zp2tf可以进行系统函数的不同表示形式之间的转换。假设 ，利用函数z,p,k=tf2zp(b
20、,a)，可以将H(z) 转换成零、极点的表示形式。,其中输入变量b、a分别是按z的降幂排列的分子、分母多项式的系数向量；输出变量z表示H(z)的零点，p表示H(z)的极点，k表示增益。 函数b,a=zp2tf(z,p,k)用来实现相反的过程。,（1） 函数zplane 函数zplane可以用来画出z变换的零、极点图，该函数有以下两种调用方式：zplane(zeros,poles)，其中zeros、poles分别为z变换的零点和极点；zplane(b,a)，其中b、a分别为z变换中分子和分母多项式的系数向量，注意这里的多项式按照z的降幂排列。,例2.23 已知离散系统的差分方程为 求其z变换，画
21、出零、极点示意图，并判断系统的稳定性。 解： 由差分方程可得,图2.30 例2.23 的系统的零极点示意图,由于系统的极点全部在单位圆内，所以系统是稳定的。,（3）函数residuez residuez 函数可以计算有理函数的留数和直接项（即多项式项），因此可以用来求z反变换。设多项式为： residuez函数的调用有以下两种方式：,用语句 R,p,C = residuez(b,a) 可以求得的留数、极点和直接项，其中输入数据b、a分别是分子多项式和分母多项式的系数向量（这些多项式都按z的降幂排列），输出数据 R 包含着留数，p 包含着极点，C 包含着直接项。语句 b,a = residuez
22、(R,p,C) 有三个输入变量和两个输出变量，它把部分分式变成多项式的系数行向量 b 和 a。,例2.24 求 （ ）的 z 反变换。 解： b = 0,1; a = 3,-4,1; R,p,C = residuez(b,a) % b,a = residuez(R,p,C),运行结果如下（留数、极点、直接项以及分子分母多项式的系数）：R = 0.5000 -0.5000p = 1.0000 0.3333C = a = -0.0000 0.3333b = 1.0000 -1.3333 0.3333,因此可以得到： 所以z反变换的结果为：,而由a、b 可以得到X(z)原来的形式：,3 求解差分方程
23、Matlab中用 filter 函数求解给定输入时差分方程的解，该函数调用形式为： y= filter(b,a,x)其中 b = b0,b1,.,bM；a = a0,a1,.,aN 是差分方程的系数，而 x 则是输入序列数组。输出 y(n) 和输入 x(n) 的长度一样，这里必须保证系数 a0 不为零。,例2.25 离散系统的差分方程为 y(n) y(n-1) + 0.5y(n-2) = x(n)（1）计算并画出冲激响应 h(n)， n=-10,50。（2）由此 h(n) 确定系统是否稳定。 解：,图2.31 例2.25的系统的冲激响应,b=1;a=1,-1,0.5;% 求冲激响应x=impseq(0,-10,50);n=-10:50;h=filter(b,a,x);stem(n,h)axis(-10,50,-1,1.5),title(Impulse Response);xlabel(n);ylabel(h(n)% 求出冲激响应的和，以判断系统是否稳定sum(abs(h)ans = 3.3333看出，的曲线逐渐趋于零，而 求得为3.3333，这意味着系统是稳定的,
