matlab分析z函数频率,信号与系统实训指导(MATLAB)第6章离散时间信号和系统的Z变换分析.ppt...

本章介绍离散信号的Z变换和系统的Z变换分析以及如何用MATLAB进行信号和系统的Z变换分析。6.1节简要介绍Z变换分析用到的MATLAB函数。6.2节介绍用MATLAB求离散信号的Z变换。6.3节介绍用MATLAB求离散信号的Z反变换。6.4节介绍离散系统的系统函数和用系统函数进行系统分析。6.5节为系统的Z变换分析实例。 　　　　信号与系统的Z变换分析涉及到的MATLAB函数主要有：　　 ztrans　　 iztrans　　 residuez　　 freqz　　下面简要介绍这些函数。 　　1.ztrans和iztrans　　这两个函数是SymbolicMathToolbox中的函数。ztrans函数实现因果离散时间信号的Z变换,命令为　　 Z＝ztrans(x)　　其中,x为信号的时域符号表达式,可以用sym函数定义。iztrans函数实现Z反变换,命令为　　 x＝iztrans(Z)　　其中,Z为Z域符号表达式。这两个函数给出的结果也是符号表达式。 　　2.residuez　　已知离散系统的系统函数,可以用residuez函数进行部分分式分解,命令为　　 [r,p,k]＝residuez(num,den)　　其中,num是系统函数分子多项式的系数向量,den是分母多项式的系数向量,系数均是按z－1的升幂排列的。 　　　　离散时间信号x(n)的Z变换定义为(6.1)其中,z是一个复数。信号x(t)的单边Z变换也记做X(z),定义为(6.2) 如果x(n)是因果信号,其单边Z变换和双边Z变换是一样的。我们这里着重讨论单边Z变换。从公式(6.1)和(6.2)可知,Z变换是z－1的幂级数,使级数收敛的复数z的集合称为Z变换的收敛域。　　对于因果信号x(n),当X(z)的收敛域包括|z|＝1时,它的离散时间傅里叶变换X(ejω)等于z＝ejω时的Z变换：(6.3)即单位圆上的Z变换就是信号x(n)的离散时间傅里叶变换。 　　离散时间信号的Z变换可以通过查表来确定,结合Z变换的性质,可以得到许多常用离散信号的变换对。MATLAB的SymbolicMathToolbox中提供了用于计算Z变换和反变换的函数ztrans和iztrans。　　例6-1 求序列x(n)＝12nu(n)的Z变换。　　解：MATLAB程序如下：　　 x＝sym(′(1/2)^n′)；　　 X＝ztrans(x)　　运行结果如下：　　 X＝2*z/(2*z－1) 　　　　若离散时间信号x(n)的Z变换为X(z),则信号x(n)可以由X(z)的Z反变换求得。Z反变换的定义为(6.4)式中,C是包围X(z)zn－1所有极点的逆时针闭合积分路线,通常选择Z平面收敛域内以原点为中心的圆。和拉普拉斯反变换一样,公式(6.4)的计算也是比较复杂的。 当离散时间信号x(n)的Z变换X(z)可以表示为z或z－1的多项式之比,即有理分式的形式时,可以通过幂级数展开法或部分分式展开法获得x(n) 。　　当X(z)表示为有理分式 (6.5) 时,函数residuez可以求X(z)的部分分式展开式所需的留数和极点,进而求出Z反变换。式(6.5)中,M和N为正整数,系数和是实数。分子多项式B(z)和分母多项式A(z)都写成z－1的升幂形式。一般假设B(z)和A(z)没有公用因子,如果有,应该约去。 　　函数residuez求X(z)的部分分式展开的基本形式是：(6.6)在查表求反变换时要注意。 例6-2 求的Z反变换。　　解：MATLAB程序如下：　　 num＝[1－0.5]；　　 den＝[13/41/8]；　　 [r,p,k]＝residuez(num,den) 　　运行结果产生下面的列向量：　　 r＝ 4　　－3　　 p＝－0.5000　　－0.2500　　 k＝[] 　　所以部分分式展开式为查表得 在用手工计算Z反变换时,我们更习惯的是把X(z)表示为如下形式：　　(6.7)先把式(6.7)两边都除z,得(6.8) 然后对式(6.8)进行部分分式展开,最后把展开结果两边再乘以z,查表获得反变换结果。这里就可以用函数residue求式(6.8)中的部分分式展开式,得到下面的基本形式：(6.9)在查表时可能会更方便。我们可以通过下面的例题体会一下这两种方法的异同。 　　例6-3 求的Z反变换。　　解：用函数residuez求解时,