资源描述:
《武汉科技大学《数字信号处理》实验报告三》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、信息科学与工程学院实验报告三专业班级姓名学号_实验时间2013年12月5口指导教师陈华丽成绩验称实名计设的器波滤数RII器波滤拟模三验实验容实内减频减频为t袞幅袞幅别M*-⑽±-的-^ft雜浦雜灘解通據通城齡■2绘fe绘咏器HZ求)H求pli颇=5要=20要«z,繼齡30率15b球k$II频--OH要丄S,1一S11o截M献2为dB逍小通小带§跳欄>£:鶴通啸s阻g阻蹤删12HZ1■以-阻波1波1器B滤一-幻鼢波3d拟t«i贝拟J贝肺袋模趾分模扯分熟大幻设沃歸度•織W度趙带隨g满阻嗰—阻喇概通雌
2、腳f26/线个⑽线个OHMJ船一=2曲*=3曲一3度熟计=〃性计性计仏幅t设汉特设汉特设5HCf-L2.3.止的RP和端波截器。数低滤带波d/s阶应通阻滤rad的对带,:曰y器素和HZ绘妓波元通35H求g滤个低3要M为一得z>。S别第获OHdBW分的以J130和vnJ*柯为为PWSN数别fW和、函分拓,NS此率似禪数晌用德S/關餓J利。訨带:城SK数浦阻有腿回假解系關3,数聽返S«
3、邶的:3dba秈)ow.截A)足为tla蟀dB和IdB<满减Ma通<辉3试W',哀界减w'^l和项’•kJ^),
4、边衰S,娜数多器卯啟训带小樓地阶母波1?綠常S1)通最嫩識的分滤,嘯忡S/S是带瞰镧器妒阻,5研R,知阻識熇波>带4Z砂娥RP分和畸m滤(B拟0H拥器么“威矛对为娜3磯波,wW氣通素S*)别式斯z、(te诫WP和如碑ivn;分g沃H。拟(v^pft"T个WC多怜15线滿一⑽带对仁(N,aW抒巴为曲2?ut嫌趾r^oQter和W个别性概=b参为脾ributN姗H-勒跳Ewn'p射止频Al=中函f胭幌:在N.V,其t截界B,A其统i驴巾示[NRB边[B系4.提①和3d的②器-8O-9OO巴特沃斯低泊
5、滤波器的to频特性50100150200250猙率/HzoO20■-30-40-50-60-703004G巴特沃斯岛通滤波器的幅频特性50100150200频率/Hz2(K«8((K2(ap/M/w250[B,Al=butter(N,Wn,'high,;s,)可以获得高通滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。[B,A]=butter(N,Wn,'stop7s*)可以获得带阻滤波器系统函数的分子多项式(B)和分母多项式(A)的系数。利用freqs函数计算模拟滤波器的频率响应:H
6、=freqs(B,A,w)其中,B和A分别表示滤波器系统阑数的分子多项式和分母多项式的系数。该函数返回矢量w指定的那些频率点上的频率响成,w的单位是rad/s。不带输出变量的freqs函数,将绘制出幅频和相频曲线。1.Wp=2*pi*5;Ws=2*pi*12;Rp=2;Rs=30;实验[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s*);内容[B,A]=butter(N,Wn,’s’);W=0:300;h=freqs(B,A,W);H=20*logl0(abs(h));plot(W,
7、H);titleC巴特沃斯低通滤波器的幅频特性*〉;xlabeK^^^/Hz1);ylabeldS®:/ciBf);2.Wp=2*pi*20;Ws=2*pi*10;Rp=3;Rs=15;[N,Wn]=buttord(Wp,Ws,Rp,Rs,’s’);[B,A]=butter(N,Wnz1high1z▼s▼);W=0:250;h=freqs(B,A,W);H=20*logl0(abs(h));plot(W,H);title(»巴特沃斯高通滤波器的幅频特性▼);巴特沃斯帝通滤波器的嘛抽特性1002
8、00300400500600700800猢牟/Hz50D5000500050•-1-1-2-2§墓1-35050O300巴特沃飫带PU劝波戕的怊頻犄性50100150200250300频氺/HzO5000501W"藝•aw墓-200-250xlabel(▼频率/Hz!);:/1北61(1幅度/犯,〉;3.Wp=[2*pi*10,2*pi*25];Ws=[2*pi*5,2*pi*30];Rp=3;Rs=30;[N,Wn]=buttord(Wp,Ws,Rp,Rs,’s▼);[B,A]=butter(
9、N,Wn,*s*);W=0:800;h=freqs(BzA,W);H=20*logl0(abs(h));plot(W,H);titleC巴特沃斯带通滤波器的幅频特性M;xlabel(*频率/Hz’>;ylabel(M^fg/dBf);4.实验Wp=[2*pi*10z2*pi*35];内容Ws=[2*pi*15,2*pi*30];Rp=3;Rs=30;[N,Wn]=buttord(Wp,Ws,Rp,Rs,’s▼);[B,A]=butter(N,Wn,▼stop▼,▼s’};W=0:300;h=fr