基于matlab下數(shù)字濾波器設(shè)計(jì)程序(編輯修改稿)
2025-02-14 15:01 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 阻帶截止頻率mtr=floor(Wp*(M+1)/(2*pi))+2。%向負(fù)方向入floor()=3。floor()=4Ad=[Wm=Wp]。Ad(mtr)=。Hd=Ad.*exp(j**M*Wm)。%構(gòu)造頻域采樣向量H(k)Hd=[Hd conj(fliplr(Hd(2:(M+1)/2)))]。%fliplr函數(shù)實(shí)現(xiàn)矩陣的左右翻轉(zhuǎn)conj是求復(fù)數(shù)的共軛h=real(ifft(Hd))。%h（n）=IDFT[H(k)]w=linspace(0,pi,1000)。%用于產(chǎn)生0,pi之間的1000點(diǎn)行矢量H=freqz(h,[1],w)。%濾波器的幅頻特性圖figure(1)plot(w/pi,20*log10(abs(H)))。%參數(shù)分別是歸一化頻率與幅值xlabel(39。歸一化角頻率39。)。ylabel(39。增益/分貝39。)。title(39。濾波器的增益響應(yīng)39。)。axis([0 1 50 ])。f1=100。f2=300。f3=700。fs=2000。%待濾波正弦信號(hào)頻率及采樣頻率figure(2)subplot(211)t=0:1/fs:。%定義時(shí)間范圍和步長(zhǎng)s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*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)。 AFs=abs(Fs)。%將信號(hào)變換到頻域及信號(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(h,1,s)。%使用filter函數(shù)對(duì)信號(hào)進(jìn)行濾波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)。 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。)。（2）利用頻率取樣法設(shè)計(jì)高通濾波器M=32。%所需頻率采樣點(diǎn)個(gè)數(shù)Wp=*pi。%通帶截止頻率m=0:M/2。%阻頻帶上的采樣點(diǎn)Wm=2*pi*m./(M+1)。%阻帶截止頻率mtr=ceil(Wp*(M+1)/(2*pi))。%向正方向舍入ceil()=4。ceil()=3。Ad=[Wm=Wp]。Ad(mtr)=。Hd=Ad.*exp(j**M*Wm)。%構(gòu)造頻域采樣向量H(k)Hd=[Hd conj(fliplr(Hd(2:M/2+1)))]。%fliplr函數(shù)實(shí)現(xiàn)矩陣的左右翻轉(zhuǎn)conj是求復(fù)數(shù)的共軛h=real(ifft(Hd))。%h（n）=IDFT[H(k)]w=linspace(0,pi,1000)。%用于產(chǎn)生0,pi之間的1000點(diǎn)行矢量H=freqz(h,[1],w)。%濾波器的幅頻特性圖figure(1)plot(w/pi,20*log10(abs(H)))。%參數(shù)分別是歸一化頻率與幅值xlabel(39。歸一化頻率39。)。ylabel(39。增益/分貝39。)。title(39。濾波器的增益響應(yīng)39。)。axis([0 1 50 0])。f1=200。f2=700。f3=800%待濾波正弦信號(hào)頻率fs=2000。%采樣頻率figure(2)subplot(211)t=0:1/fs:。%定義時(shí)間范圍和步長(zhǎng)s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t)。%濾波前信
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1