基于matlab下數(shù)字濾波器設(shè)計(jì)程序-全文預(yù)覽
【正文】 bplot(211)t=0:1/fs:。axis([0 1 50 0])。)。歸一化頻率39。%h（n）=IDFT[H(k)]w=linspace(0,pi,1000)。Ad(mtr)=。%阻帶截止頻率mtr=ceil(Wp*(M+1)/(2*pi))。（2）利用頻率取樣法設(shè)計(jì)高通濾波器M=32。)。頻率/赫茲39。%限定圖像坐標(biāo)范圍subplot(212)Fsf=fft(sf,512)。title(39。)。)。幅度39。%濾波前的信號頻域圖xlabel(39。subplot(212)Fs=fft(s,512)。)。時間/秒39。%待濾波正弦信號頻率及采樣頻率figure(2)subplot(211)t=0:1/fs:。f1=100。title(39。)。%用于產(chǎn)生0,pi之間的1000點(diǎn)行矢量H=freqz(h,[1],w)。Hd=Ad.*exp(j**M*Wm)。%通頻帶上的采樣點(diǎn)及阻帶截止頻率mtr=floor(Wp*(M+1)/(2*pi))+2。)(1) 利用頻率取樣法設(shè)計(jì)低通濾波器M=63。幅度39。%頻率采樣plot(f,AFsf(1:256))%濾波后的信號頻域圖xlabel(39。)。幅度39。figure(3)subplot(211)plot(t,sf)%濾波后的信號圖像xlabel(39。)。頻率/赫茲39。AFs=abs(Fs)。title(39。)。%濾波前信號sf=filter(b,1,s)。f6=95。f2=20。title(39。)。b=fir1(M,Wc,ftype,window)。dev=Rs*ones(1,length(a))。（3）利用fir1函數(shù)及Kaiser窗設(shè)計(jì)多通帶濾波器Rs=。)。頻率/赫茲39。subplot(212)Fsf=fft(sf,512)。title(39。)。信號濾波前頻域圖39。ylabel(39。plot(f,AFs(1:256))。)。幅度39。%濾波前的信號圖像xlabel(39。%待濾波正弦信號頻率t=(0:200)/fs。f1=500。)。頻率/赫茲39。window=Kaiser(M+1,beta)。a=[0 1 0]。)。幅度39。%頻率采樣plot(f,AFsf(1:256))%濾波后的信號頻域圖xlabel(39。axis([ 2 2])。)。時間/秒39。信號濾波前頻域圖39。ylabel(39。%頻率采樣plot(f,AFs(1:256))。)。幅度39。%濾波前的信號圖像xlabel(39。)。增益/分貝39。%濾波器的幅頻特性圖%[H,W]=freqz(B,A,N)當(dāng)N是一個整數(shù)時函數(shù)返回N點(diǎn)的頻率向量和幅頻響應(yīng)向量plot(f*fs/(2*pi),20*log10(abs(h)))%參數(shù)分別是頻率與幅值xlabel(39。%定義過度帶寬M=round(8/m)。附錄二(1) 利用fir1函數(shù)及hamming窗設(shè)計(jì)低通濾波器f1=100。%采樣頻率m=(*f1)/(fs/2)。%使用fir1函數(shù)設(shè)計(jì)濾波器%輸入的參數(shù)分別是濾波器的階數(shù)和截止頻率figure(1)[h,f]=freqz(b,1,512)。ylabel(39。濾波器的增益響應(yīng)39。%濾波前信號plot(t,s)。ylabel(39。信號濾波前時域圖39。%信號頻域圖的幅值f=(0:255)*fs/512。)。title(39。%使用filter函數(shù)對信號進(jìn)行濾波%參數(shù)分別為濾波器系統(tǒng)函數(shù)的分子和分母多項(xiàng)式系數(shù)向量和待濾波信號輸入subplot(211)plot(t,sf)%濾波后的信號圖像xlabel(39。幅度39。)。%信號頻域圖的幅值f=(0:255)*fs/512。ylabel(39。信號濾波后頻域圖39。%采樣頻率fcuts=[1000 1300 2210 2410]。%M為能夠滿足設(shè)計(jì)要求的濾波器的最小階數(shù)，Wc為濾波器的截止頻率點(diǎn)%第一個元素f為待設(shè)計(jì)濾波器的過渡帶的起始點(diǎn)和結(jié)束點(diǎn)%第二個元素a指定第一個元素頻率段的理想幅度值%第三個元素dev中的元素為各通帶和阻帶內(nèi)允許的幅度最大誤差M=mod(M,2)+M。%濾波器的
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號-1