实验4离散系统的频率响应分析和零极点分布-Read.DOCVIP

参考《数字信号处理》—基于计算机的方法（第三版） z变换 一、实验原理： 1.LTI离散时间系统的时域差分方程为： 传递函数： 对上面的差分方程两边求z变换，得： 我们定义LTI离散时间系统的输出的Z变换Y(z)与输入的Z变换X(z)的比值为该系统的传递函数，即为系统的传递函数。 分解因式 ，其中和称为零、极点。 利用系统的传递函数，我们可以分析系统的零极点，稳定性及实现结构等特点。 2.使用到的MATLAB命令： 分析系统传递函数的函数：tf2zp，zp2tf，zplane，zp2sos，residuez。 函数tf2zp的使用形式为[z,p,k]=tf2zp(num,den)，其中输入num,den分别表示有理分式Z变换的分子和分母多项式系数（按z的降幂排列），输出为系统的零点，极点和增益因子。函数zplane可以对输入的有理分式Z变换的分子和分母多项式或零极点坐标画系统的零极点图。zp2tf的使用形式为[num,den]= zp2tf(z,p,k)，其输入输出关系刚好与tf2zp相反。函数zp2sos的使用形式为sos=zp2sos(z,p,k),可以将一个用零极点及增益因子描述的系统分解为多个级联的二阶子系统。函数residuez可以用来创建一个有理分式Z变换的部分分式展开式。其使用形式有[r,p,k]=residuez(num,den)和[num,den]=residuez (r, p,k)，前者输入为有理分式Z变换的分子和分母多项式系数（按z的降幂排列），输出的矢量r表示留数和分子常数，矢量p表示极点，矢量k包含常数；而后者输入输出刚好相反。另外还有根据系统的传递函数求冲击响应的函数impz，filter等 二、实验目的： 1.能利用MATLAB确定有理z变换的因式分解形式，并确定零极点。 2.能利用MATLAB，根据零极点确定有理z变换的表达式 3.能利用MATLAB，用部分分式展开的方法求逆z变换 4.能利用MATLAB，用长除法求逆z变换 5.能利用MATLAB，完成系统的稳定性判断 三、实验内容： 1. Using Program6_1, determine the factored form and pole-zero plot of EXAMPLE6.9 2. Using Program6_2, determine the rational form of z-transform of EXAMPLE6.10 from its zero and pole locations 3.Using Program6_3，determine the inverse z-transform of EXAMPLE6.16 with partial-fraction expansion method 4. Using Program6_4, determine the rational form of z-transform of EXAMPLE6.17 from its partial-fraction expansion 5. Using Program6_5, determine the inverse z-transform of EXAMPLE6.20 with long division first method 6. Using Program6_6，determine the inverse z-transform of EXAMPLE6.20 with long division second method 7. Using Program6_7，test the stability of transfer function in EXAMPLE6.39 四、实验程序及结果分析： 五、参考程序： 1. % Program6_1 % Determination of the Factored Form % of a Rational z-Transform % num = input(Type in the numerator coefficients = ); den = input(Type in the denominator coefficients = ); K = num(1)/den(1); Numfactors = factorize(num) Denfactors = factorize(den) disp(Numerator factors);disp(Numfactors); disp(Denominator factors);disp(Denfactors); disp(Gain constant);disp(K); zplane(num,den) 2. % P