實(shí)驗(yàn)報(bào)告基于MATLAB的數(shù)字濾波器設(shè)計(jì)驕陽教育

1、實(shí)驗(yàn)78 基于matlab的數(shù)字濾波器設(shè)計(jì)實(shí)驗(yàn)?zāi)康模杭由顚?shù)字濾波器的常用指標(biāo)和設(shè)計(jì)過程的理解。實(shí)驗(yàn)原理：低通濾波器的常用指標(biāo)：通帶邊緣頻率：，阻帶邊緣頻率： ，通帶起伏：，通帶峰值起伏： ，阻帶起伏：，最小阻帶衰減： 。數(shù)字濾波器有iir和fir兩種類型，它們的特點(diǎn)和設(shè)計(jì)方法不同。 在matlab中，可以用b，a=butter（n,wn）等函數(shù)輔助設(shè)計(jì)iir數(shù)字濾波器, 也可以用b=fir1(n,wn,ftype) 等函數(shù)輔助設(shè)計(jì)fir數(shù)字濾波器。實(shí)驗(yàn)內(nèi)容： 利用matlab編程設(shè)計(jì)一個(gè)數(shù)字帶通濾波器，指標(biāo)要求如下：通帶邊緣頻率：，通帶峰值起伏：。阻帶邊緣頻率：，最小阻帶衰減： 。 分別用i

2、ir和fir兩種數(shù)字濾波器類型進(jìn)行設(shè)計(jì)。實(shí)驗(yàn)要求：給出iir數(shù)字濾波器參數(shù)和fir數(shù)字濾波器的沖激響應(yīng)，繪出它們的幅度和相位頻響曲線，討論它們各自的實(shí)現(xiàn)形式和特點(diǎn)。實(shí)驗(yàn)內(nèi)容：irr代碼：wp=0.45*pi,0.65*pi;ws=0.3*pi,0.75*pi;ap=1;as=40;n,wc=buttord(wp/pi,ws/pi,ap,as) ;b,a=butter(n,wc)%b,a = butter(n,wn,ftype) h,w=freqz(b,a);subplot(211);mag=abs(h);plot(w/pi,mag); xlabel(角頻率(w);ylabel(幅度|ha(jw

3、)|);title(iir|ha(jw)|);grid on;subplot(212);phase=angle(h);plot(w/pi,phase);xlabel(角頻率(w);ylabel(相位);title(iir相位響應(yīng)曲線);grid on;圖像：fir代碼： wp1 = 0.4*pi; wp2 = 0.6*pi; ws1 = 0.3*pi; ws2 = 0.7*pi; tr_width = min(wp1-ws1),(ws2-wp2) tr_width =0.3142 m = ceil(6.2*pi/tr_width) + 1 m = 63 n=0:1:m-1; wc1 = (ws

4、1+wp1)/2; wc2 = (wp2+ws2)/2; wc=wc1/pi,wc2/pi; window= hanning(m); h1,w=freqz(window,1); figure(1); subplot(2,1,1) stem(window); axis(0 60 0 1.2); grid; xlabel(n); title(hanning窗函數(shù)); subplot(2,1,2) plot(w/pi,20*log(abs(h1)/abs(h1(1); axis(0 1 -350 0); grid; xlabel(w/pi); ylabel(幅度(db); title(hanning

5、窗函數(shù)的頻譜); hn = fir1(m-1,wc, hanning (m); h2,w=freqz(hn,1,512); figure(2); subplot(2,1,1) stem(n,hn); axis(0 60 -0.25 0.25); grid; xlabel(n); ylabel(h(n); title(hanning窗函數(shù)的單位脈沖響應(yīng)); subplot(2,1,2) plot(w/pi,20*log(abs(h2)/abs(h2(1); grid; xlabel(w/pi); ylabel(幅度(db); figure(3); phase=angle(h1); plot(ph