湖南大学信号与系统实验九.docVIP

实验九 Z变换 §9.1 作离散时间零极点图 目的 学习如何将一个离散时间有理系统函数的零极点展现一张零极点图。 下面的练习需要用到函数plotdpz，其具体实现如下。函数plotdpz(b,a)画出离散时间系统的零极点。 plotdpz的输入和filter具有相同的格式，并且如果分子和分母多项式不是同阶次的话，plotdpz将自动的对a或b补上适当个数的零。同时，plotdpz在图中还包括单位圆，以及如果在单位圆内有多于一个零点或极点的话，还会指出极点或零点的数目。 function plotdpz(b,a) la=length(a); lb=length(b); if (lalb) b=[b zeros(la-lb)]; elseif (lbla) a=[a zeros(lb-la)]; end ps=roots(a); zs=roots(b); mx=max(abs([pszs.95]))+.05; clg; axis([-mx mx -mx mx]); axis(equal); hold on; w=[0:.01:2*pi]; plot(cos(w),sin(w),.); plot([-mx mx],[0 0]); plot([0 0],[-mx mx]); text(0.1,1.1,Im,sc); text(1.1,.1,Re,sc); plot(real(ps),imag(ps),X); plot(real(zs),imag(zs),o); numz=sum(abs(zs)==0); nump=sum(abs(ps)==0); if numz1 text(-.1,-.1,num2str(numz)); elseif nump1 text(-.1,-.1,num2str(nump)); end hold off 基本题 1．用plotdpz画出(9.3)或的零极点图。 解：不用plotdpz,自定义函数： 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;%确定横坐标范围 clf 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),x)%画极点 plot(real(q),imag(q),o)%画零点 title(pole-zero diagram for discrete system)%标注标题 hold off a=[1 3 2]; b=[1 -1 0]; ljdt(a,b) p=roots(a) q=roots(b) pa=abs(p) 使用plotdpz(): b=[1 -1 0]; a=[1 3 2]; plotdpz(b,a); 2．用plotdpz画出满足差分方程的滤波器的零极点图。 解：代码如下： 使用自定义函数ljdt(): a=[1 1 0.5]; b=[1 0 0]; ljdt(a,b) 使用plotdpz(): b=[1 0 0]; a=[1 1 0.5]; plotdpz(b,a); 3．用plotdpz画出满足差分方程 所描述滤波器的零极点图。 解： 使用自定义函数ljdt(): a=[1 -1.25 0.75 -0.125] b=[1 0.5] ljdt(a,b) 使用plotdpz（）： b=[1 0.5 0 0]; a=[1 -1.25 0.75 -0.125]; plotdpz(b,a); §9.2离散时间频率响应的几何解释 基本题 对于这些题将考察一个简单的二阶离散时间系统，其系统函数是 1．定义b1和a1包含按filter要求格式的分子和分母多项式的系数，按10.1节所介绍的方法求出并画出的零极点图。 解：代码如下： b1=[1 0 0]; a1=[1 -0.9 0.81]; plotdpz(b1,a1); 2．定义omega=[0:511]*pi/256和unitcirc=exp(j*omega)得到在单位圆上512等分点，在这些点上将要对频率响应求值。定义polevectors1是一个2×