資源描述：
《Matlab實(shí)現(xiàn)數(shù)字FIR的高通 、帶通、低通、帶阻濾波器的程序.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、利用漢寧窗設(shè)計(jì)Ⅰ型數(shù)字高通濾波器clearall;Wp=0.6*pi;Ws=0.4*pi;tr_width=Wp-Ws;????????????????????????%過渡帶寬度N=ceil(6.2*pi/tr_width)???????????????????????%濾波器長度n=0:1:N-1;Wc=(Ws+Wp)/2;????????????????????????%理想低通濾波器的截止頻率hd=ideal_hp1(Wc,N);?????????????????????%理想低通濾波器的單位沖激響應(yīng)w_han=(hanning(N))';?????????
2、???????????%漢寧窗h=hd.*w_han;???????????????????????????%截取得到實(shí)際的單位脈沖響應(yīng)[db,mag,pha,w]=freqz_m2(h,[1]);?????????????%計(jì)算實(shí)際濾波器的幅度響應(yīng)delta_w=2*pi/1000;Ap=-(min(db(Wp/delta_w+1:1:501)))????????%實(shí)際通帶紋波As=-round(max(db(1:1:Ws/delta_w+1)))????????%實(shí)際阻帶紋波subplot(221)stem(n,hd)title('理想單位脈沖響應(yīng)hd(n)')
3、subplot(222)stem(n,w_han)title('漢寧窗w(n)')subplot(223)stem(n,h)title('實(shí)際單位脈沖響應(yīng)hd(n)')subplot(224)plot(w/pi,db)title('幅度響應(yīng)(dB)')axis([0,1,-100,10])clearall;Wp=0.6*pi;Ws=0.4*pi;tr_width=Wp-Ws;????????????????????????%過渡帶寬度N=ceil(6.2*pi/tr_width)???????????????????????%濾波器長度n=0:1:N-1;Wc=(W
4、s+Wp)/2;????????????????????????%理想低通濾波器的截止頻率hd=ideal_hp1(Wc,N);?????????????????????%理想低通濾波器的單位沖激響應(yīng)w_han=(hanning(N))';????????????????????%漢寧窗h=hd.*w_han;???????????????????????????%截取得到實(shí)際的單位脈沖響應(yīng)[db,mag,pha,w]=freqz_m2(h,[1]);?????????????%計(jì)算實(shí)際濾波器的幅度響應(yīng)delta_w=2*pi/1000;Ap=-(min(db(Wp
5、/delta_w+1:1:501)))????????%實(shí)際通帶紋波As=-round(max(db(1:1:Ws/delta_w+1)))????????%實(shí)際阻帶紋波subplot(221)stem(n,hd)title('理想單位脈沖響應(yīng)hd(n)')subplot(222)stem(n,w_han)title('漢寧窗w(n)')subplot(223)stem(n,h)title('實(shí)際單位脈沖響應(yīng)hd(n)')subplot(224)plot(w/pi,db)title('幅度響應(yīng)(dB)')axis([0,1,-100,10])基于切比雪夫一致逼近法設(shè)
6、計(jì)FIR數(shù)字低通濾波器clearall;??????????????f=[00.60.71];????????????????????%給定頻率軸分點(diǎn)A=[1100];????????????????????%給定在這些頻率分點(diǎn)上理想的幅頻響應(yīng)weigh=[110];?????????????????%給定在這些頻率分點(diǎn)上的加權(quán)b=remez(32,f,A,weigh);???????????%設(shè)計(jì)出切比雪夫最佳一致逼近濾波器[h,w]=freqz(b,1,256,1);h=abs(h);h=20*log10(h);subplot(211)stem(b,'.');g
7、rid;title('切比雪夫逼近濾波器的抽樣值')subplot(212)plot(w,h);grid;title('濾波器幅頻特性(dB)')利用漢寧窗設(shè)計(jì)Ⅰ型數(shù)字帶阻濾波器clearall;????????????????????Wpl=0.2*pi;Wph=0.8*pi;Wsl=0.4*pi;Wsh=0.6*pi;tr_width=min((Wsl-Wpl),(Wph-Wsh));???????????????????????????%過渡帶寬度N=ceil(6.2*pi/tr_width)????????????????????????????????
