基于matlab下數(shù)字濾波器設(shè)計(jì)程序(已修改)
2025-01-30 15:01 本頁(yè)面
【正文】 附錄二(1) 利用fir1函數(shù)及hamming窗設(shè)計(jì)低通濾波器f1=100。f2=200。%待濾波正弦信號(hào)頻率fs=2000。%采樣頻率m=(*f1)/(fs/2)。%定義過(guò)度帶寬M=round(8/m)。%定義窗函數(shù)的長(zhǎng)度N=M1。%定義濾波器的階數(shù)b=fir1(N,*f2/(fs/2))。%使用fir1函數(shù)設(shè)計(jì)濾波器%輸入的參數(shù)分別是濾波器的階數(shù)和截止頻率figure(1)[h,f]=freqz(b,1,512)。%濾波器的幅頻特性圖%[H,W]=freqz(B,A,N)當(dāng)N是一個(gè)整數(shù)時(shí)函數(shù)返回N點(diǎn)的頻率向量和幅頻響應(yīng)向量plot(f*fs/(2*pi),20*log10(abs(h)))%參數(shù)分別是頻率與幅值xlabel(39。頻率/赫茲39。)。ylabel(39。增益/分貝39。)。title(39。濾波器的增益響應(yīng)39。)。figure(2)subplot(211)t=0:1/fs:。%定義時(shí)間范圍和步長(zhǎng)s=sin(2*pi*f1*t)+sin(2*pi*f2*t)。%濾波前信號(hào)plot(t,s)。%濾波前的信號(hào)圖像xlabel(39。時(shí)間/秒39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波前時(shí)域圖39。)。subplot(212)Fs=fft(s,512)。%將信號(hào)變換到頻域AFs=abs(Fs)。%信號(hào)頻域圖的幅值f=(0:255)*fs/512。%頻率采樣plot(f,AFs(1:256))。%濾波前的信號(hào)頻域圖xlabel(39。頻率/赫茲39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波前頻域圖39。)。figure(3)sf=filter(b,1,s)。%使用filter函數(shù)對(duì)信號(hào)進(jìn)行濾波%參數(shù)分別為濾波器系統(tǒng)函數(shù)的分子和分母多項(xiàng)式系數(shù)向量和待濾波信號(hào)輸入subplot(211)plot(t,sf)%濾波后的信號(hào)圖像xlabel(39。時(shí)間/秒39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波后時(shí)域圖39。)。axis([ 2 2])。%限定圖像坐標(biāo)范圍subplot(212)Fsf=fft(sf,512)。%濾波后的信號(hào)頻域圖AFsf=abs(Fsf)。%信號(hào)頻域圖的幅值f=(0:255)*fs/512。%頻率采樣plot(f,AFsf(1:256))%濾波后的信號(hào)頻域圖xlabel(39。頻率/赫茲39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波后頻域圖39。)。（2）利用fir1函數(shù)及Kaiser窗設(shè)計(jì)帶通濾波器Rs=。fs=8000。%采樣頻率fcuts=[1000 1300 2210 2410]。a=[0 1 0]。dev=Rs*ones(1,length(a))。[M,Wc,beta,ftype]=kaiserord(fcuts,a,dev,fs)。%M為能夠滿足設(shè)計(jì)要求的濾波器的最小階數(shù)，Wc為濾波器的截止頻率點(diǎn)%第一個(gè)元素f為待設(shè)計(jì)濾波器的過(guò)渡帶的起始點(diǎn)和結(jié)束點(diǎn)%第二個(gè)元素a指定第一個(gè)元素頻率段的理想幅度值%第三個(gè)元素dev中的元素為各通帶和阻帶內(nèi)允許的幅度最大誤差M=mod(M,2)+M。window=Kaiser(M+1,beta)。b=fir1(M,Wc,ftype,window)。%輸入的第一個(gè)參數(shù)是濾波器的階數(shù)%第二個(gè)參數(shù)是濾波器的截止頻率%第三個(gè)參數(shù)是濾波器的類(lèi)型，stop為帶阻濾波器%第四個(gè)參數(shù)是采用的窗函數(shù)[h,f]=freqz(b,1,512)。%濾波器的幅頻特性圖%[H,W]=freqz(B,A,N)當(dāng)N是一個(gè)整數(shù)時(shí)函數(shù)返回N點(diǎn)頻率向量和幅頻響應(yīng)向量figure(1)plot(f*fs/(2*pi),20*log10(abs(h)))%參數(shù)分別是頻率與幅值xlabel(39。頻率/赫茲39。)。ylabel(39。增益/分貝39。)。title(39。濾波器的增益響應(yīng)39。)。f1=500。f2=1500。f3=2000。f4=3000。%待濾波正弦信號(hào)頻率t=(0:200)/fs。%定義時(shí)間的步長(zhǎng)s=sin(2*f1*pi*t)+sin(2*f2*pi*t)+sin(2*f
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1