1.MATLAB源碼

M=63; Wp=0.5*pi;%所需頻率采樣點(diǎn)個(gè)數(shù)及通帶截止頻率 m=0:(M+1)/2; Wm=2*pi*m./(M+1);%通頻帶上的采樣點(diǎn)及阻帶截止頻率 mtr=floor(Wp*(M+1)/(2*pi))+2;%向負(fù)方向入floor(3.5)=3;floor(-3.2)=-4 Ad=[Wm<=Wp]; Ad(mtr)=0.38; Hd=Ad.*exp(-1i*0.5*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('歸一化角頻率');ylabel('增益/分貝');title('濾波器的增益響應(yīng)'); axis([0 1 -50 0.5]); f1=100;f2=300;f3=700;fs=2000;%待濾波正弦信號(hào)頻率及采樣頻率 figure(2) subplot(211) t=0:1/fs:0.25;%定義時(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('時(shí)間/秒');ylabel('幅度');title('信號(hào)濾波前時(shí)域圖'); 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('頻率/赫茲');ylabel('幅度');title('信號(hào)濾波前頻域圖'); figure(3) sf=filter(h,1,s);%使用filter函數(shù)對(duì)信號(hào)進(jìn)行濾波 subplot(211) plot(t,sf)%濾波后的信號(hào)圖像 xlabel('時(shí)間/秒');ylabel('幅度');title('信號(hào)濾波后時(shí)域圖'); axis([0.2 0.25 -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('頻率/赫茲');ylabel('幅度');title('信號(hào)濾波后頻域圖');

2.結(jié)果

總結(jié)

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

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。