信号与系统实验10.pdfVIP

一、实验要求 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 R2016 三、实验内容及分析 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(tmpmaxlevel)=maxlevel*ones(1,length(tmp(tmpmaxlevel))); q(tmp0)=zeros(tmp(tmp0)); q=(q-levels)*M/levels; qc(k,:)=q; end （d） (e) [b,a]=ellip(4,.2,40