信号与系统实验.doc

(a)(b)(c)(a)(b)(c)(d)(e)(a)(b)(c)(d)(e)(f)(g)(a)(b)(c)(d)实验设备(环境)Windows10;(a)代码如下:a=[132];b=[1-10];draw(b,a);title('H(z)µÄÁ㼫µãͼ');:functiondraw(num,den)p=roots(den);q=roots(num);p=p';q=q';x=max(abs([pq]));x=x+1;y=x;plot([-xx],[00],':');holdon;plot([00],[-yy],':');holdon;plot(real(p),imag(p),'*');holdon;plot(real(q),imag(q),'o');holdon;xlabel('ʵ²¿');ylabel('Ð鲿');axis([-xx-yy]);(b)代码如下:a=[];b=[1];draw(b,a);title('Â˲¨Æ÷Á㼫µãͼ');(c)代码如下:a=[1--];b=[];draw(b,a);title('Â˲¨Æ÷µÄÁ㼫µãͼ');(a)b=[100];a=[1-];zs=roots(b);ps=roots(a);plot(real(zs),imag(zs),'o');holdon;plot(real(ps),imag(ps),'x');gridon;axis([-44-44]);xlabel('ʵ²¿');ylabel('Ð鲿');title('H1(z)µÄÁ㼫µãͼ');(a)[b,a]=ellip(4,.2,40,[.]);[H,w]=freqz(b,a,4096);plot(w/pi,20*log10(abs(H)));axis([01-8010]);title('HµÄ¶ÔÊý·ùֵͼ');(b)[b,a]=ellip(4,.2,40,[.]);[H,w]=freqz(b,a,4096);plot(w/pi,20*log10(abs(H)));axis([01-16020]);title('HµÄ¶ÔÊý·ùֵͼ');(c)[b,a]=ellip(4,.2,40,[.]);M=max(abs([ba]));a16=quant(a,16,M);b16=quant(b,16,M);[H,w]=freqz(b16,a16,4096);plot(w/pi,20*log10(abs(H)));axis([01-8010]);title('HµÄ¶ÔÊý·ùֵͼ');:functionqc=quant(x,N,M)[mm,nn]=size(x);qc=zeros(mm,nn);levels=2^(N-1);maxlevel=2