基于matlab下數(shù)字濾波器設計程序-wenkub
2023-02-02 15:01:58 本頁面
【正文】 39。%濾波前的信號頻域圖xlabel(39。subplot(212)Fs=fft(s,512)。)。時間/秒39。%定義時間的步長s=sin(2*f1*pi*t)+sin(2*f2*pi*t)+sin(2*f3*pi*t)+sin(2*f4*pi*t)。f2=1500。title(39。)。b=fir1(M,Wc,ftype,window)。dev=Rs*ones(1,length(a))。（2）利用fir1函數(shù)及Kaiser窗設計帶通濾波器Rs=。)。頻率/赫茲39。%限定圖像坐標范圍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:。)。頻率/赫茲39。%定義窗函數(shù)的長度N=M1。f2=200。%待濾波正弦信號頻率fs=2000。%定義濾波器的階數(shù)b=fir1(N,*f2/(fs/2))。)。title(39。%定義時間范圍和步長s=sin(2*pi*f1*t)+sin(2*pi*f2*t)。)。title(39。%將信號變換到頻域AFs=abs(Fs)。頻率/赫茲39。)。figure(3)sf=filter(b,1,s)。ylabel(39。信號濾波后時域圖39。%濾波后的信號頻域圖AFsf=abs(Fsf)。)。title(39。fs=8000。[M,Wc,beta,ftype]=kaiserord(fcuts,a,dev,fs)。%輸入的第一個參數(shù)是濾波器的階數(shù)%第二個參數(shù)是濾波器的截止頻率%第三個參數(shù)是濾波器的類型，stop為帶阻濾波器%第四個參數(shù)是采用的窗函數(shù)[h,f]=freqz(b,1,512)。ylabel(39。濾波器的增益響應39。f3=2000。sf=filter(b,1,s)。)。title(39。AFs=abs(Fs)。頻率/赫茲39。)。figure(3)subplot(211)plot(t,sf)%濾波后的信號圖像xlabel(39。幅度39。)。%信號頻域圖的幅值f=(0:255)*fs/512。ylabel(39。信號濾波后頻域圖39。%采樣頻率fcuts=[10 20 40 50 60 70 80 90]。M=mod(M,2)+M。%濾波器的幅頻特性圖figure(1)plot(f*fs/(2*pi),20*log10(abs(h)))%參數(shù)分別是頻率與幅值xlabel(39。增益/分貝39。)。f4=55。%定義時間的步長s1=sin(2*f1*pi*t)+sin(2*f2*pi*t)+sin(2*f3*pi*t)。%濾波前的信號圖像xlabel(39。幅度39。)。plot(f,AFs(1:256))。ylabel(39。信號濾波前頻域圖39。)。title(39。 AFsf=abs(Fsf)。)。title(39。%所需頻率采樣點個數(shù)及通帶截止頻率m=0:(M+1)/2。floor()=4Ad=[Wm=Wp]。%fliplr函數(shù)實現(xiàn)矩陣的左右翻轉conj是求復數(shù)的共軛h=real(ifft(Hd))。%參數(shù)分別是歸一化頻率與幅值xlabel(39。增益/分貝39。)。f3=700。%濾波前信號plot(t,s)。ylabel(39。信號濾波前時域圖39。%將信號變換到頻域及信號頻域圖的幅值f=(0:255)*fs/512。)。title(39。%使用filter函數(shù)對信號進行濾波subplot(211)plot(t,sf)%濾波后的信號圖像xlabel(39。幅度39。)。%濾波后的信號頻域圖及信號頻域圖的幅值f=(0:255)*fs/512。ylabel(39。信號濾波后頻域圖39。%通帶截止頻率m=0:M/2。ceil()=3。%構造頻域采樣向量H(k)Hd=[Hd conj(fliplr(Hd(2:M/2+1)))]。%濾波器的幅頻特性圖figure(1)plot(w/pi,20*log10(abs(H)))。ylabel(39。濾波器的增益響應39。f2=700。%濾
研究報告相關推薦
鄂ICP備17016276號-1