參與評論

評論

登錄后參與評論

全部評論(1條)

• 

  我昱云南的人民 2017-06-28 00:00:00

  這個(gè)我剛好做過一個(gè)濾波器，事實(shí)上對時(shí)域信號(hào)做FFT，截取一定點(diǎn)數(shù)再做逆FFT相當(dāng)于理想濾波。設(shè)計(jì)濾波器代碼如下： f1=100;f2=200;%待濾波正弦信號(hào)頻率 fs=2000;%采樣頻率 m=(0.3*f1)/(fs/2);%定義過度帶寬 M=round(8/m);%定義窗函數(shù)的長度 N=M-1;%定義濾波器的階數(shù) b=fir1(N，f2/fs);%使用fir1函數(shù)設(shè)計(jì)濾波器 %輸入的參數(shù)分別是濾波器的階數(shù)和截止頻率 figure(1) [h，f]=freqz(b，1，512);%濾波器的幅頻特性圖 %[H，W]=freqz(B，A，N)當(dāng)N是一個(gè)整數(shù)時(shí)函數(shù)返回N點(diǎn)的頻率向量和幅頻響應(yīng)向量 plot(f*fs/(2*pi)，20*log10(abs(h)))%參數(shù)分別是頻率與幅值 xlabel('頻率/赫茲');ylabel('增益/分貝');title('濾波器的增益響應(yīng)'); figure(2) subplot(211) t=0:1/fs:0.5;%定義時(shí)間范圍和步長 s=sin(2*pi*f1*t)+sin(2*pi*f2*t);%濾波前信號(hào) plot(t，s);%濾波前的信號(hào)圖像 xlabel('時(shí)間/秒');ylabel('幅度');title('信號(hào)濾波前時(shí)域圖'); subplot(212) Fs=fft(s，512);%將信號(hào)變換到頻域 AFs=abs(Fs);%信號(hào)頻域圖的幅值 f=(0:255)*fs/512;%頻率采樣 plot(f，AFs(1:256));%濾波前的信號(hào)頻域圖 xlabel('頻率/赫茲');ylabel('幅度');title('信號(hào)濾波前頻域圖'); figure(3) sf=filter(b，1，s);%使用filter函數(shù)對信號(hào)進(jìn)行濾波 %參數(shù)分別為濾波器系統(tǒng)函數(shù)的分子和分母多項(xiàng)式系數(shù)向量和待濾波信號(hào)輸入 subplot(211) plot(t，sf)%濾波后的信號(hào)圖像 xlabel('時(shí)間/秒');ylabel('幅度');title('信號(hào)濾波后時(shí)域圖'); axis([0.2 0.5 -2 2]);%限定圖像坐標(biāo)范圍 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('頻率/赫茲');ylabel('幅度');title('信號(hào)濾波后頻域圖');

  贊(12)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論