数字信号处理上机实验滤波器

《数字信号处理上机实验滤波器》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、数字信号处理上机练习7.12(1)双线性变换法程序：%双线性变换法设计切比雪大1型带通滤波器%w计算式：w=(f/ADfs)*2*pi单位：radclc；clear；closeall；ADfs=2000;%ADC的采样频率，中•位：Hzwsl=(100/ADfs)*2*pi;S下附带截I卜.频率，单位：radwpl=(200/ADfs)*2*pi;%通带卜截止频率，单位:radwp2=(400/ADfs)*2*pi;%通带上截止频率,单位.•radws2=(600/ADfs)*2*pi;%上阻带截止频率，1丫1•位

2、：raddisplay(1数字滤波器的边界频率中•位：rad');wsl,wplzwp2,ws2Rp=2;%通祇最大袋减：革位：dBAs=20;%阻带最小衰减，牟位：dBT=l;%模拟滤波器数字化时的采样间隔，单位：sWsl=(2/T)*tan(wsl/2)；Wpl=(2/T)*tan(wpl/2)；Wp2=(2/T)*tan(wp2/2);Ws2=(2/T)*tan(ws2/2);%转换为模拟滤波器指标Ws=[Wsl,Ws2];Wp=[Wpl,Wp2];display(*转换后滤波器的边界频率单位：rad/s')

3、；WslzWpl,Wp2,Ws2[n,wn]=cheblord(Wp,Ws,Rp#As,'s');%计算切比雪夫1型滤波器阶数和截止频率display(1滤波器阶数*);ndisplay('2dB截止频率1);wn[b,a]=chebyl(n,Rp,wn,'s');%计算模拟带通滤波器的传输函数Ha(s>display('模拟原®系统函数Ha(s)分子多项成系数：1);bdisplay('模拟原型系统函数Ha(s)分W多项式系数')；a[Ha,Wa]=freqS(b,a,512);%计算并画出模拟原型滤波器的幅度频

4、谱图subplot(211);plot(Wa/pi,20*logl0(abs(Ha)),'LineWidth',2);axis([0,max(Wa/pi),min(20*logl0(abs(Ha))),5]);xlabel('州频率rad/s(Xji)')；ylabelC对数幅度dB1);titleC模拟原型带通滤波器的幅度谱(对数幅度)•)；grid;[bz'az]=bilinear(b,a,1/T);%川双线性变换法转换为数字滤波器display('数字滤波器的系统函数H(z)分•多项式系数：');bzdisp

5、lay(•数字滤波器的系统函数H(z)分母多项式系数•)；az[Hz,Wz]=freqz(bz,az,512,l/T);%计算并岡出数'f•滤波器的幅度频®图，Wz是以Hz为单位subplot(212);plot(Wz*2*T,20*logl0(abs(Hz)),1LineWidth1,2);axis([0,1,-100,5]);xlabel(1数字域频率(Xk)•)；ylabelU对数幅度dB')；title(

6、双线性变换法设计的数字滤波器的幅度谱(对数幅度)•)；grid；结果：111■■———————————

7、———*■1111幽1~~——~~————~~——~應111藝1•111111111•■■■■■■■■■■■■■■麵籌隹■■■■■■■■■■■■■«!■—1111111111111111r11111111111事1111麵1AM«■»MBMB«a«■»MB■*«■»«■MB«■«■■11111100.511.522.53角频率rad/s(xn}模拟原型带通滤波器的幅度谱(对数幅度)0•10-20■30-40-50-60-70-80双线性变换法没计的数字滤波器的幅度唞(对数幅度)0-20II111I■I1-40-60

8、-80-10000.10.20.30.40.50.60.70.80.91数字域频率(。川数字滤波器的边界频率单位：ndwsl=0.3142wpl=0.6283wp2=1.2566ws2=1.8850转换滤波器的边界频率单位：rad/sWsl=0.3168Wpl=0.6498Wp2=1.4531Ws2=2.7528滤波器阶数2dB截止频率wn=0.64981.4531模拟原型系统函数Ha(s)分子多项式系数:b=000.421800模拟原型系统函数Ha(s)分母多项式系数a=1.00000.64572.41960.6

9、0970.8916数字滤波器的系统阑数H(z)分子多项式系数：bz=0.05120.0000-0.10240.00000.0512数字滤波器的系统阑数H(z＞分母多项式系数az=1.0000-2.07332.4881-1.59440.6125(2)冲激响应不变法程序:%冲激响应不变法设计切比雪夫1型带通滤波器计兑式：w=(f/ADfs)*2*pi单位：rad