基于matlab下數(shù)字濾波器設(shè)計(jì)程序-文庫(kù)吧
2025-01-03 15:01 本頁面
【正文】 3*pi*t)+sin(2*f4*pi*t)。sf=filter(b,1,s)。%使用filter函數(shù)對(duì)信號(hào)進(jìn)行濾波figure(2)subplot(211)plot(t,s)。%濾波前的信號(hào)圖像xlabel(39。時(shí)間/秒39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波前時(shí)域圖39。)。subplot(212)Fs=fft(s,512)。AFs=abs(Fs)。f=fs/512*(0:255)。plot(f,AFs(1:256))。%濾波前的信號(hào)頻域圖xlabel(39。頻率/赫茲39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波前頻域圖39。)。figure(3)subplot(211)plot(t,sf)%濾波后的信號(hào)圖像xlabel(39。時(shí)間/秒39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波后時(shí)域圖39。)。axis([ 4 4])。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。)。（3）利用fir1函數(shù)及Kaiser窗設(shè)計(jì)多通帶濾波器Rs=。fs=200。%采樣頻率fcuts=[10 20 40 50 60 70 80 90]。a=[0,1,0,1,0]。dev=Rs*ones(1,length(a))。[M,Wc,beta,ftype]=kaiserord(fcuts,a,dev,fs)。M=mod(M,2)+M。window=Kaiser(M+1,beta)。b=fir1(M,Wc,ftype,window)。[h,f]=freqz(b,1,512)。%濾波器的幅頻特性圖figure(1)plot(f*fs/(2*pi),20*log10(abs(h)))%參數(shù)分別是頻率與幅值xlabel(39。頻率/赫茲39。)。ylabel(39。增益/分貝39。)。title(39。濾波器的增益響應(yīng)39。)。f1=5。f2=20。f3=30。f4=55。f5=75。f6=95。%待濾波正弦信號(hào)頻率t=(0:200)/fs。%定義時(shí)間的步長(zhǎng)s1=sin(2*f1*pi*t)+sin(2*f2*pi*t)+sin(2*f3*pi*t)。s=s1+sin(2*f4*pi*t)+sin(2*f5*pi*t)+sin(2*f6*pi*t)。%濾波前信號(hào)sf=filter(b,1,s)。%使用filter函數(shù)對(duì)信號(hào)進(jìn)行濾波figure(2)subplot(211)plot(t,s)。%濾波前的信號(hào)圖像xlabel(39。時(shí)間/秒39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波前時(shí)域圖39。)。subplot(212)Fs=fft(s,512)。AFs=abs(Fs)。f=fs/512*(0:255)。plot(f,AFs(1:256))。%濾波前的信號(hào)頻域圖xlabel(39。頻率/赫茲39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波前頻域圖39。)。figure(3)subplot(211)plot(t,sf)%濾波后的信號(hào)圖像xlabel(39。時(shí)間/秒39。)。ylabel(39。幅度39。)。title(39。信號(hào)濾波后時(shí)域圖39。)。subplot(212)Fsf=fft(sf,512)。 AFsf=abs(Fsf)。%濾波后的信號(hào)頻域圖及信號(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。)(1) 利用頻率取樣法設(shè)計(jì)低通濾波器M=63。 Wp=*pi。%所需頻率采樣點(diǎn)個(gè)數(shù)及通帶截止頻率m=0:(M+1)/2。 Wm=2*pi*m./(M+1)。%通頻帶上的采樣點(diǎn)及
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1