实验离散系统Z域分析.docVIP

实验九 离散系统Z域分析 1实验目的 (1) 理解并掌握系统函数的概念。 (2) 掌握利用MATLAB绘制系统函数的零极点分布图。 (3) 掌握系统函数零极点的分布与系统时域、频域特性之间的关系。 (4) 利用MATLAB求系统函数零、极点的方法。 2实验原理及方法 2.1离散时间LTI系统的Z域描述 线性时不变离散系统可以用如下所示的线性常系数差分方程来描述。 9-1 其中y(k)为系统输出序列，f (k)为输入系列。 将式9-1两边进行Z变换得： 9-2 式9-2中A(z)和B(z)分别是由描述系统的差分方程的系数决定的关于z的多项式，将式9-2因式分解后有： 9-3 其中C为常数，qj (j＝1,2,..,M)为H(z)的M个零点，pi (i＝1,2,…,N)为H(z)的N个极点。 由以上分析可以看出，系统函数H(z)的零、极点的分布完全决定了系统的特性，若某离散系统的零点、极点已知，则系统函数便可确定下来。 因此，系统函数的零极点分布对我们进行离散系统特性的分析具有非常重要的意义。通过对系统函数零极点的分析，我们可以分析离散系统以下几个方面的特性； (1) 系统单位响应h(k)的时域特性 (2) 离散系统的稳定性 (3) 离散系统的频率特性（幅频响应和相频响应） 要通过系统函数零极点来分析系统特性，首先就要求出系统函数的零极点，然后绘制零点、极点图。MATLAB为我们快速、高效地分析离散系统特性提供了强有力的工具。下面就介绍如何利用MATLAB实现这一过程。 设离散系统的系统函数为： 则系统函数的零点和极点可以用MATLAB的多项式求根函数roots()来实现，调用函数roots()的命令格式为： p=roots(A) 其中A为待求根的多项式的系数构成的行向量，返回向量p则是包含该多项式所有根位置的列向量。例如多项式为： 则求该多项式根的MATLAB命令应为： A=[1 3/4 1/8]; P=roots(A) 运行结果为： p＝ -0.5000 -0.2500 需要注意的是，在求系统函数零极点时，离散系统的系统函数可能有两种形式，一种是分子和分母多项式均按Z的降幂次序排列，如式9-4所示：另一种是分子多项式和分母多项式均按Z的升幂次序排列，如式9-5所示。上述两种方式在构造多项式系数向量时稍有不同。 9-4 9-5 若H(z)是以z的降幂形式排列，则系数向量一定要由多项式的最高幂次开始，一直到常数项，缺项要用0补齐。例如对式9-4所示的系统函数，其分子多项式的系数向量应为：B=［1 0 2 0]，分母多项式的系数向量应为：A=[1 3 2 2 1]。 若H(z)是以的升幂形式排列，则分子和分母多项式系数向量的维数一定要相同，不足的要用0补齐，否则Z=0的零点或极点就可能被漏掉，例如，对式9-5所示的系统函数，其分子多项式系数向量应为：B=[1 1 0]，分母多项式系数向量应为：A=[1 1/2 1/41。 2.2系统函数的零极点分布图 下面给出一个用于绘制离散时间LTI系统的零极点图的扩展函数ljdt(A,B)： function ljdt(A,B) p=roots(A); %求系统极点 q=roots(B); %求系统零点 p=p; q=q; x=max(abs([p q 1])); %确定纵坐标范围 x=x+0.1; y=x; %确定横坐标范围 hold on axis([-x x -y y]) %确定坐标轴显示范围 w=0:pi/300:2*pi; t=exp(i*w); plot(t) %画单位圆 axis(square) plot([-x x],[0 0]) %画横坐标轴 plot([0 0],[-y y]) %画纵坐标轴 text(0.1,x,jIm[z]) text(y,1/10,Re[z]) plot(real(p),imag(p), rx,MarkerSize,10) %画极点 plot(real(q),imag(q),go,M