MATLAB离散系统z域分析.docVIP

实验八离散系统的z域分析 一、 目的 掌握利用MATLAB绘制系统零极点图的方法 掌握离散时间系统的零极点分析方法 掌握逆Z变换概念及MATLAB实现方法 二、 离散系统零极点 线性时不变离散系统可用线性常系数差分方程描述，即 N M (8-1)(8-2)[aiy(n-i) = ^bjX(n-j) (8-1) (8-2) /=0 7=0 其中:va)为系统的输出序列，对幻为输入序列。 Y(z)X(z)将式(8-1)两边进行Z Y(z) X(z) 一則 z=o 将式(8-2)因式分解后有: M n(z-a) H(z) = C^ (8-3) n(z-A) /=l 其中C为常数，人(J = l，2,…，M)为//(z)的M个零点，A(/ = l，2,…，A0为//⑵的TV个 极点。 系统函数H(z)的零极点分布完全决定了系统的特性，若某系统函数的零极点已知， 则系统函数便可确定下来。 因此，系统函数的零极点分布对离散系统特性的分析具冇非常重要意义。通过对系 统函数零极点的分析，可以分析离散系统以下几个方面的特性： ?系统单位样值响应A(n)的时域特性； ?离散系统的稳定性； ?离散系统的频率特性； 三、离散系统零极点图及零极点分析 零极点图的绘制 设离散系统的系统函数为 H(z) =B(z) H(z) = B(z) A⑵ 则系统的零极点可用MATLAB的多项式求根函数rootsO来实现，调用格式为： p=roots(A) 艽中A为待求多项式的系数构成的行矩阵，返冋向量p则是包含多项式所有根的列向 Q I 量。如多项式为B(z) = z2+-z + -,则求该多项式根的MATLAB命令为: 4 8 A=[l 3/4 1/8]; P=roots(A) 运行结果为： P = -0.5000 -0.2500 需注意的是，在求系统函数零极点时，系统函数可能有两种形式：一种是分子、分 母多项式均按z的降幂次序排列；另一种是分子、分母多项式均按厂1的升幂次序排列。 这两种方式在构造多项式系数向量时稍有不同。 H(z)按z的降幂次序排列：系数向量一定要由多项式最高次幂开始，一直到常数 项，缺项要用0补齐；如 H(z) =-—— z4+3? + 2?+2z + 1 其分子、分母多项式系数向量分别为A=[l 0 2 01、B=fl 3 22 IK //(幻按厂1的升幂次序排列：分子和分母多项式系数向量的维数一定要相同，不 足的要用0补齐，否则z = 0的零点或极点就可能被漏掉。如 H(z) =l + 2z,1 + 2Z_1其分子、分母多项式系数向量分别为A=[l 2 0]、B=[l 1/2 1/4]o用roots() H(z) = l + 2z, 1 + 2Z_1 其分子、分母多项式系数向量分别为A=[l 2 0]、B=[l 1/2 1/4]o 用roots()求得//(z)的零极点后，就可以用plot()函数绘制出系统的零极点图。T面 是求系统零极点，并绘制其零极点阁的MATLAB实用函数ljdt(),同时还绘制了单位圆。 function ljdt(A，B) % The function to draw the pole-zero diagram for discrete system p 二 roots(A); q=roots(B); p=pr; q=q’； x=max(abs([p q 1])); x=x+0.1; y=x； elf %求系统极点 %求系统零点 %将极点列向量转置为行向量 %将零点列向量转置为行向量 %确定纵坐标范围 %确定横坐标范围 hold on axis([-x x -y y]) w=0:pi/300:2*pi; t=exp(i*w); plot ⑴ axisCsquare’) plot([-x x]，[0 0]) plot([0 0]，[-y y]) text(0.1,x,jIm[z]) text(y，l/10，Re[z]’) plot(real(p)，imag(p)，’x’) %确定坐标轴显示范围 %岡单位园 %画横坐标轴 %画纵坐标轴 plot(real(q)，imag(q)，’o’) title(’pole-zero diagram for discrete system’) hold off %画极点 %画零点 %标注标题 例h绘制如下系统函数的零极点 (1) H(z) =3z3-5z (1) H(z) = 3z3-5z2 + 10z z3-3z2 + 7z-5 (2) //(z) = l-0.5z, 解：MATLAB命令如下 A=[l -3 7-5]; B=[3-5 100]; ljdt(A，B) 绘制的零极点图如图8-1 (a)所示。 A=[l 3/4 1/8]; B=[l -0.5 0]; ljdt