数字信号处理练习2.docVIP

下载本文档

2
0
约6.39千字
约 24页
2018-05-08 发布于河南
举报
版权申诉

数字信号处理练习2.doc

1、本文档共24页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数字信号处理练习2

T9.22 解： λ===0.203 所求 ()===num1 = 0.3404*[1 2 1]; den1 = [1 0.1842 0.1776]; [h1,w] = freqz(num1,den1,512); plot(w/pi,20*log10(abs(h1))); xlabel(\omega/\pi);ylabel(gain,db); hold on; num2 = [0.21623 0.43245 0.21623]; den2 = [0.96993 -0.2863 0.181416]; [h2,w2]= freqz(num2,den2,512); plot(w2/pi,20*log10(abs(h2)) ,g); xlabel(\omega2/\pi);ylabel(gain,db); 截图： M9.1 解： Fp = 10000hz ; Fs = 30000hz ; Ft = 100000hz ; ap = 0.4db; as = 50db; wp = 2*pi*Fp/Ft ; 得wp = 0.628 ws = 2*pi*Fs/Ft = 1.885 = tan(wp/2)= 1.376 ; = tan(ws/2) = 1.376; k =/=4.235; = (10^0.1*ap-1) = 0.096; = (10^0.1*as)^0.5= 100000; k1 = ((a*a-1)/b)^0.5 =0.0009798; N = log10(k1)/log10(k)=4.8 ; 取N =5 =0.435 [z,p,k] = buttap(6); [num,den]= zp2tf(z,p,k); [z,p,k]= tf2zp(num,den); [zd,pd,kd]= bilinear(z,p,k,0.5); [n2,d2]= zp2tf(zd,pd,kd); [h,w] = freqz(n2,d2); figure(1); plot(w/pi,20*log10(abs(h))); axis([0 1 -50 10]); title(gain response); figure(2); plot(w/pi,unwrap(angle(h))); title(phase respond); M9.5 解：同M9.1 有wp = 2*pi*Fp/Ft = 0.628 ws = 2*pi*Fs/Ft = 1.885 = tan(wp/2)= 1.376 ; = tan(ws/2) = 1.376; k =/=4.235; = (10^0.1*ap-1) = 0.096; = (10^0.1*as)^0.5= 100000; 1/k1 = ((-1)/ ) =1020.62; N = log10(k1)/log10(k)=4.8 ; 取N =5 =0.435 程序： [z, p, k] = cheb1ap(6, 0.4); [B, A] = zp2tf(z, p, k); [BT, AT] = lp2lp(B, A, 0.628); [num, den] = impinvar(BT, AT, 1); [h, w] = freqz(num, den, 512); figure(1); plot(w/pi, 20*log10(abs(h))); grid; axis([0 1 -60 5]); xlabel(\omega/\pi); ylabel(Gain, in dB); title(Impulse Invariance Method); figure(2); plot(w/pi, unwrap(angle(h))); grid axis([0 1 -8 1]); xlabel(\omega/\pi); ylabel(Phase, in radians); title(Impulse Invariance Method); % Bilinear transformation method [z, p, k] = cheb1ap(4, 0.4); [B, A] = zp2tf(z, p, k); [BT, AT] = lp2lp(B, A, 0.325); [num, den] = bilinear(BT, AT, 0.5); [h, w] = freqz(num, den, 512); figure(3); plot(w/pi, 20*log10(abs(h))); grid; axis([0 1 -60 5]); xlabel(\omega/\pi); ylabel(Gain, in dB); title(Bilinear Transformation Method); figure(4); plot(w/