Z变换实验报告

实验九Z变换§9.1作离散时间零极点图目的学习如何将一个离散时间有理系统函数的零极点展现一张零极点图。下面的练习需要用到函数plotdpz，其具体实现如下。函数plotdpz(b,a)画出离散时间系统的零极点。plotdpz的输入和filter具有相同的格式，并且如果分子和分母多项式不是同阶次的话，plotdpz将自动的对a或b补上适当个数的零。同时，plotdpz在图中还包括单位圆，以及如果在单位圆内有多于一个零点或极点的话，还会指出极点或零点的数目。functionplotdpz(b,a)la=length(a);lb=length(b);if(la>lb)b=[bzeros(la-lb)];elseif(lb>la)a=[azeros(lb-la)];endps=roots(a);zs=roots(b);mx=max(abs([ps'zs'.95]))+.05;clg;axis([-mxmx-mxmx]);axis('equal');holdon;w=[0:.01:2*pi];plot(cos(w),sin(w),'.');plot([-mxmx],[00]);plot([00],[-mxmx]);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);ifnumz>1text(-.1,-.1,num2str(numz));elseifnump>1text(-.1,-.1,num2str(nump));endholdoff基本题1．用plotdpz画出(9.3)或的零极点图。解：不用plotdpz,自定义函数：functionljdt(A,B)p=roots(A);%求系统极点q=roots(B);%求系统零点p=p';%将极点列向量转置为行向量q=q';%将零点列向量转置为行向量x=max(abs([pq1]));%确定纵坐标范围x=x+0.1;y=x;%确定横坐标范围clfholdonaxis([-xx-yy])%确定坐标轴显示范围w=0:pi/300:2*pi;t=exp(i*w);plot(t)%画单位圆axis('square')plot([-xx],[00])%画横坐标轴plot([00],[-yy])%画纵坐标轴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-zerodiagramfordiscretesystem')%标注标题holdoffa=[132];b=[1-10];ljdt(a,b)p=roots(a)q=roots(b)pa=abs(p)使用plotdpz():b=[1-10];a=[132]; plotdpz(b,a);2．用plotdpz画出满足差分方程的滤波器的零极点图。解：代码如下：使用自定义函数ljdt():a=[110.5];b=[100];ljdt(a,b)使用plotdpz():b=[100];a=[110.5];plotdpz(b,a);3．用plotdpz画出满足差分方程所描述滤波器的零极点图。解：使用自定义函数ljdt():a=[1-1.250.75-0.125]b=[10.5]ljdt(a,b)使用plotdpz（）：b=[10.500];a=[1-1.250.75-0.125];plotdpz(b,a);§9.2离散时间频率响应的几何解释基本题对于这些题将考察一个简单的二阶离散时间系统，其系统函数是1．定义b1和a1包含按filter要求格式的分子和分母多项式的系数，按10.1节所介绍的方法求出并画出的零极点图。解：代码如下：b1=[100];a1=[1-0.90.81];plotdpz(b1,a1);2．定义omega=[0:511]*pi/256和unitcirc=exp(j*omega)得到在单位圆上512等分点，在这些点上将要对频率响应求值。定义polevectors1是一个2×512的矩阵，其中每一行包含一些复数，这些复数是由unitcirc的相应列减去一个极点位置得到的。如果ps1是一个列向量，它包含了极点的位置，那么就能用下列命令来完成>>polevectors1=ones(2,1)*unitcirc-ps1*ones(1,512);在MATLAB中总是愿意用矩阵和向量运算，因为MATLAB几乎总是计算矩阵运算比用for循环构成的运算要快。用abs和atan2定义polelength1和poleangle1作为polevectors1中每一元素的幅值和相位。解：b1=[100];a1=[1-0.90.81];ps1=roots(a1);omega=[0:511].*pi/256;unitcirc=exp(i.*omega);polevectors1=zeros(2,512);polevectors1=ones(2,1)*unitcirc-ps1*ones(1,512);polelength1=abs(polevectors1)poleangle1=atan2(imag(polevectors1),real(polevectors1))3．类似于polevecters1定义zerovectors1，使得它是2×512的矩阵，其中包含从零点位置到unitcirc元素的向量。定义zerolength1和zeroangle1分别是这些向量的幅值和相位。解：b1=[100];a1=[1-0.90.81];zs1=roots(a1);omega=[0:511].*pi/256;unitcirc=exp(i.*omega);zerovectors1=zeros(2,512);zerovectors1=ones(2,1)*unitcirc-ps1*ones(1,512);zerolength1=abs(zerovectors1)zeroangle1=atan2(imag(zerovectors1),real(zerovectors1))4．画出polelength1和zerolength1对omega的图。根据这些图，你预计在哪里有最大值和最小值？解：b1=[100];a1=[1-0.90.81];ps1=roots(a1);zs1=roots(b1);omega=[0:511]*pi/256;unitcirc=exp(i*omega);zerovectors1=zeros(2,512);zerovectors1=ones(2,1)*unitcirc-zs1*ones(1,512);zerolength1=abs(zerovectors1);subplot(2,1,1)plot(omega,zerolength1);polevectors1=zeros(2,512);polevectors1=ones(2,1)*unitcirc-ps1*ones(1,512);polelength1=abs(polevectors1);subplot(2,1,2)plot(omega,polelength1);5．用polelength1和zerolength1计算，并将这个结果存入geonH1mag中。用poleangle1和zeroangle1计算，并将这个结果存入gemoH1phase中。在定义geonH1mag和gemoH1phase时，你或许会发现prod和sum是很有用的。画出从几何上到处的幅值和相位，并将它和经由下面计算所得结果比较：>>H1=freqz(b1,a1,512,'whole');你在4中估计的峰值频率正确吗？解：b1=[100];a1=[1-0.90.81];ps1=roots(a1);zs1=roots(b1);omega=[0:511]*pi/256;unitcirc=exp(i*omega);zerovectors1=zeros(2,512);zerovectors1=ones(2,1)*unitcirc-zs1*ones(1,512);zerolength1=abs(zerovectors1);zeroangle1=atan2(imag(zerovectors1),real(zerovectors1));polevectors1=zeros(2,512);polevectors1=ones(2,1)*unitcirc-ps1*ones(1,512);polelength1=abs(polevectors1);poleangle1=atan2(imag(polevectors1),real(polevectors1));geonH1mag=zerolength1./polelength1;gemoH1phase=zeroangle1-poleangle1;subplot(2,1,1)plot(omega,geonH1mag);subplot(2,1,2)plot(omega,gemoH1phase);freqz(b1,a1,512,'whole');中等题下面这些题的第一组将考查移动一个零点的效果。具体说来，考虑系统函数为5．求出并画出的零极点图。你怎样预计这个系统的polevectors2或者zerovectors2与的有什么不同吗？解：b2=[1-0.50];a2=[1-0.90.81];plotdpz(b2,a2);结果分析：这个系统的polevectors2与相同，但zerovectors2与零点值不同，矩阵第二行的列值发生变化6．计算对polevectors2和zerovectors2以及所有这些向量的幅值和相位。画出幅值和相位对omega的图。你在5中的预计正确吗？解：b2=[1-0.50];a2=[1-0.90.81];ps2=roots(a2);zs2=roots(b2);omega=[0:511]*pi/256;unitcirc=exp(i*omega);zerovectors2=zeros(2,512);zerovectors2=ones(2,1)*unitcirc-zs2*ones(1,512);zerolength2=abs(zerovectors2);zeroangle2=atan2(imag(zerovectors2),real(zerovectors2));polevectors2=zeros(2,512);polevectors2=ones(2,1)*unitcirc-ps2*ones(1,512);polelength2=abs(polevectors2);poleangle2=atan2(imag(polevectors2),real(polevectors2));geonH1mag=zerolength2./polelength2;gemoH1phase=zeroangle2-poleangle2;subplot(2,1,1)plot(omega,geonH1mag);subplot(2,1,2)plot(omega,gemoH1phase);freqz(b2,a2,512,'whole');结果分析：在5中的预计正确7．根据这些零点的变化，预计与将会有怎样的不同。用freqz计算H2验证你的答案。解：b2=[1-0.50];a2=[1-0.90.81];H2=freqz(b2,a2,512,'whole')结果分析：与相比，相位和幅值都发生了变化以下的题将考虑系统函数为8．求出并画出的零极点图。这些零点和极点的位置是如何关联的？解：代码如下：b3=[0.25-1.732/21];a3=[1-1.732/20.25];plotdpz(b3,a3);结果分析：这些零点和极点的位置关于x轴对称9．类似于你在基本题中所做的方式定义polevectors3和zerovectors3。定义polelength3和zerolength3是这些复数的幅值。在同一组坐标轴上画出所有这些幅值，也就是polelength3和zerolength3每一行的幅值。这些幅值有什么关系？根据这一点，你如何预计频率响应的幅值岁频率变化？利用这些长度计算频率响应幅值，并将它存入gemoH3mag中。画出gemoH3mag对于omega的图。解：b3=[0.25-1.732/21];a3=[1-1.732/20.25];ps3=roots(a3);zs3=roots(b3);omega=[0:511]*pi/256;unitcirc=exp(i*omega);zerovectors3=zeros(2,512);zerovectors3=ones(2,1)*unitcirc-zs3*ones(1,512);zerolength3=abs(zerovectors3);zeroangle3=atan2