一、实验规定10.1 作离散时间零极图（a）（b）（c）10.2 离散时间频率响应的几何解释基本题（a）（b）（c）（d）（e）10.3 离散时间滤波器构造的量化效应基本题（a）（b）（c）（d）（e）（f）（g）10.4 用欧拉近似法设计离散时间滤波器基本题（a）（b）（c）（d）二、实验设备（环境）Windows 10；Matlab R三、实验内容及分析10.1 作离散时间零极图（a）代码以下：a=[1 3 2];b=[1 -1 0];draw(b,a);title('H(z)µÄÁ㼫µãͼ');draw.m:function draw(num,den)p=roots(den);q=roots(num);p=p';q=q';x=max(abs([p q]));x=x+1;y=x;plot([-x x],[0 0],':');hold on;plot([0 0],[-y y],':');hold on;plot(real(p),imag(p),'*');hold on;plot(real(q),imag(q),'o');hold on;xlabel('ʵ²¿');ylabel('Ð鲿');axis([-x x -y y]);(b)代码以下：a=[1 1 0.5];b=[1];draw(b,a);title('Â˲¨Æ÷Á㼫µãͼ');（c）代码以下：a=[1 -1.25 0.75 -0.125];b=[1 0.5];draw(b,a);title('Â˲¨Æ÷µÄÁ㼫µãͼ');10.2 离散时间频率响应的几何解释（a）b=[1 0 0]; a=[1 -0.9 0.81]; zs=roots(b); ps=roots(a); plot(real(zs),imag(zs),'o'); hold on; plot(real(ps),imag(ps),'x');grid on;axis([-4 4 -4 4]);xlabel('ʵ²¿');ylabel('Ð鲿');title('H1(z)µÄÁ㼫µãͼ');10.3 离散时间滤波器构造的量化效应（a）[b,a]=ellip(4,.2,40,[.41 .47]);[H,w]=freqz(b,a,4096);plot(w/pi,20*log10(abs(H)));axis([0 1 -80 10]);title('HµÄ¶ÔÊý·ùֵͼ');（b）[b,a]=ellip(4,.2,40,[.41 .47]);[H,w]=freqz(b,a,4096);plot(w/pi,20*log10(abs(H)));axis([0 1 -160 20]);title('HµÄ¶ÔÊý·ùֵͼ');（c） [b,a]=ellip(4,.2,40,[.41 .47]);M=max(abs([b a]));a16=quant(a,16,M);b16=quant(b,16,M); [H,w]=freqz(b16,a16,4096);plot(w/pi,20*log10(abs(H)));axis([0 1 -80 10]);title('HµÄ¶ÔÊý·ùֵͼ');quant.m:function qc=quant(x,N,M) [mm,nn]=size(x);qc=zeros(mm,nn);levels=2^(N-1);maxlevel=2^N-1; for k=1:mm, tmp=fix((x(k, :)+M) ./ (M/levels)) q=zeros(1,nn); q(tmp<=maxlevel)=tmp(tmp<=maxlevel); q(tmp>maxlevel)=maxlevel*ones(1,length(tmp(tmp>maxlevel))); q(tmp<0)=zeros(tmp(tmp<0)); q=(q-levels)*M/levels; qc(k,:)=q;end（d）(e)[b,a]=ellip(4,.2,40,[.41 .47]);M=max(abs([b a]));a16=quant(a,16,M);b16=quant(b,16,M); draw(a16,b16)title('Â˲¨Æ÷µÄÁ㼫µãͼ');10.4 用欧拉近似法设计离散时间滤波器（a）a1p=[1/8];b1p=[1 1 0.5 1/8];w=[-10:0.25:10]; [H,w]=freqz(b1p,a1p,w);plot(w,abs(H));title('Â˲¨Æ÷µÄƵÂÊ·ùÖµÌØÐÔ');figure;draw(b1p,a1p);title('H1p(s)µÄÁ㼫µã');四、实验总结掌握了如何使用 fft、freqz、compose、dpzplot 简朴函数的使用，学会了使用Matlab 实现 Z 变换。教师评语：签名：日期：成绩：
