【正文】 _hanning=fir1(N1,wc/pi,39。,hanning(N))。high39。 hn_blackman=fir1(N1,wc/pi,39。,blackman(N))。 plot(n,hn_boxcar)。:39。linewidth39。 plot(n,hn_hamming,39。,39。,2)。o39。 hold off xlabel(39。)。h(n)39。 legend(39。,39。,39。,39。)。單位沖激響應 39。 M=1024。k=0:M/21。 figure()。linewidth39。 Hk=fft(hn_hanning,M) hold on plot(wk/pi,20*log10(abs(Hk(k+1))),39。,39。,3)。o39。 Hk=fft(hn_blackman,M) plot(wk/pi,20*log10(abs(Hk(k+1))),39。)。矩形窗 39。漢寧窗 39。哈明窗 39。布萊克曼窗 39。 axis([0,1,80,5])。w/\pi39。ylabel(39。)。損耗函數 39。 (2)四 種窗函數設計的單位脈沖響應的比較如圖 5所示： 圖 5 單位脈沖響應比較圖 (3)四種窗函數設計的損耗函數的比較如圖 6所示： 圖 6 損耗函數比較圖 (5)運行結果分析：由圖 6可見，當濾波器長度 N不變時，矩形窗設計的濾波器的過渡帶最窄，阻帶最小衰減最??；布萊克曼窗設計的濾波器的過渡帶最寬，同時阻帶最小衰減最大。 要求通帶截止頻率為 rad，阻帶截止頻率為 ，通帶最大衰減為 dB，阻帶最小衰減為 45 dB。 (1)matlab 代碼： wp=*pi。 Bt=wpws。 N=N0+mod(N0+1,2)。 hn=fir1(N1,wc,39。,hamming(N))。 Hk=fft(hn,M)。 stem(n,hn)。n39。ylabel(39。)。單位沖激響應 39。 k=0:M/21。 figure(2)。 axis([0,1,80,5])。w/\pi39。ylabel(39。)。損耗函數 39。grid on (2)高通 FIR濾波器的單位脈沖響應、損耗函數如圖 8所示： 圖 7 單位脈沖響應 圖 8 損耗函數 0 5 10 15 20 25 30 35 40 45 0 . 4 0 . 3 0 . 2 0 . 100 . 10 . 20 . 30 . 40 . 5nh(n)單位沖激響應0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 8 0 7 0 6 0 5 0 4 0 3 0 2 0 1 00w/ ?幅度/dB損耗函數 MATLAB工具箱函數 fir1設計線 性相位帶通 FIR濾波器。顯示所設計的單位脈沖響應 h(n)的數據，并畫出損耗函數曲線。wp2=*pi。ws2=*pi。 N=ceil(*pi/Bt)。 hn=fir1(N1,wc,hanning(N))。 Hk=fft(hn,M)。 stem(n,hn)。n39。ylabel(39。)。單位沖激響應 39。 k=0:M/21。 figure(2)。 axis([0,1,80,5])。w/\pi39。ylabel(39。)。損耗函數 39。grid on (2)帶通 FIR濾波器的單位脈沖響應、損耗函數如圖 10所示： 圖 9 單位脈沖響應 圖 10 損耗函數 0 5 10 15 20 25 30 35 40 45 0 . 2 5 0 . 2 0 . 1 5 0 . 1 0 . 0 500 . 0 50 . 10 . 1 5nh(n)單位沖激響應0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 8 0 7 0 6 0 5 0 4 0 3 0 2 0 1 00w/ ?幅度/dB損耗函數某信號 ()xt 為： 1 2 3( ) c os( 2 ) c os( 2 ) c os( 2 )x t f t f t f t? ? ? ?? ? ? ?，其中 1 2 1100 , 130 , 600 .f H z f H z f H z? ? ?設計最低階 FIR 數字濾波器，按下圖所示對 ()xt 進行數字濾波處理，實現： F I R 數 字 濾波 器()xn ()ynA / D ( 周 期T 采 樣 )D / A (周 期T )() ()ytD/A(周 期T) 1） 將 3f 頻率分量以高于 50dB 的衰減抑制，同時以低于 2dB的衰減通過 1f 和 2f 頻率分量； 一、基于 窗函數法 設計 FIR 數字濾波器： (1)matlab 代碼： Fs=3800。 fs=600。 wp=2*pi*fp/Fs。 Bt=wswp。 N=ceil((rs8)/)。 hn=fir1(N,wc,kaiser(N+1,alph))。 Hk=fft(hn,M)。 wk=(2*pi/M)*k。 plot(wk/pi,20*log10(abs(Hk(k+1))))。 xlabel(39。)。幅度 /dB39。 title(39。)。grid on xlabel(39。)。相位 /\pi39。 title(39。)。n=0:N1。T=1/f。 x=*cos(2*pi*100*t)+*cos(2*pi*130*t+*pi)+*cos(2*pi*600*t)。m2=0。m4=0。m6=0。m8=0。m10=0。m12=0。m14=0。m16=0。m18=0。m20=0。m22=0。m24=0。 m24=m23。m22=m21。m20=m19。m18=m17。m16=m15。m14=m13。m12=m11。m10=m9。m8=m7。m6=m5。m4=m3。m2=m1。 end 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 8 0 7 0 6 0 5 0 4 0 3 0 2 0 1 00w/ ?幅度/dB損耗函數0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 11 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81w/ ?相位/?相頻特性曲線plot(n,x)。信號 x(n)39。ylabel(39。)。n39。 S=fft(x,N)。