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