1.MATLAB源碼

M=63; Wp=0.5*pi;%所需頻率采樣點個數及通帶截止頻率 m=0:(M+1)/2; Wm=2*pi*m./(M+1);%通頻帶上的采樣點及阻帶截止頻率 mtr=floor(Wp*(M+1)/(2*pi))+2;%向負方向入floor(3.5)=3;floor(-3.2)=-4 Ad=[Wm<=Wp]; Ad(mtr)=0.38; Hd=Ad.*exp(-1i*0.5*M*Wm);%構造頻域采樣向量H(k) Hd=[Hd conj(fliplr(Hd(2:(M+1)/2)))]; %fliplr函數實現矩陣的左右翻轉conj是求復數的共軛 h=real(ifft(Hd));%h（n）=IDFT[H(k)] w=linspace(0,pi,1000);%用于產生0,pi之間的1000點行矢量 H=freqz(h,[1],w);%濾波器的幅頻特性圖 figure(1) plot(w/pi,20*log10(abs(H)));%參數分別是歸一化頻率與幅值 xlabel('歸一化角頻率');ylabel('增益/分貝');title('濾波器的增益響應'); axis([0 1 -50 0.5]); f1=100;f2=300;f3=700;fs=2000;%待濾波正弦信號頻率及采樣頻率 figure(2) subplot(211) t=0:1/fs:0.25;%定義時間范圍和步長 s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);%濾波前信號 plot(t,s);%濾波前的信號圖像 xlabel('時間/秒');ylabel('幅度');title('信號濾波前時域圖'); subplot(212) Fs=fft(s,512); AFs=abs(Fs);%將信號變換到頻域及信號頻域圖的幅值 f=(0:255)*fs/512;%頻率采樣 plot(f,AFs(1:256));%濾波前的信號頻域圖 xlabel('頻率/赫茲');ylabel('幅度');title('信號濾波前頻域圖'); figure(3) sf=filter(h,1,s);%使用filter函數對信號進行濾波 subplot(211) plot(t,sf)%濾波后的信號圖像 xlabel('時間/秒');ylabel('幅度');title('信號濾波后時域圖'); axis([0.2 0.25 -2 2]);%限定圖像坐標范圍 subplot(212) Fsf=fft(sf,512); AFsf=abs(Fsf);%濾波后的信號頻域圖及信號頻域圖的幅值 f=(0:255)*fs/512;%頻率采樣 plot(f,AFsf(1:256))%濾波后的信號頻域圖 xlabel('頻率/赫茲');ylabel('幅度');title('信號濾波后頻域圖');

2.結果

總結

以上是生活随笔為你收集整理的FIR数字滤波器设计频率抽样法MATLAB仿真的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。