【正文】 。歸一化頻率w39。)。)。 figure(2) subplot(2,1,1) plot(w2,20*log10(abs(hh2))) xlabel(39。歸一化頻率w39。2）濾波器實現(xiàn)形式及特點： 由于一般的濾波器在利用窗函數(shù)是其通帶波紋和阻帶波紋不同（一般為第一個阻帶波紋最大）因此，在滿足第一個阻帶衰減旁瓣時，比其頻率高的旁瓣，它們的衰減都大大超出要求。 y2=sigsum(0,16,[0:16],(abs(delta2).^2))/17。用noscale不對濾波器歸一化。 n=:。三、實驗分析：1）總結(jié)：FIR濾波器實現(xiàn)一般采用窗函數(shù)法和等紋波設(shè)計法。 grid 圖4423）用等波紋法設(shè)計：調(diào)用函數(shù)[n,fpts,mag,wt]=remezord(f，a，dev)f=[ ]a=[0 1 0]dev=[ ] [n,fpts,mag,wt]=remezord([ ],[0 1 0],[ ])。 h2=fir1(n,wn,39。 h1=2*cos(pi*(n29)/2).*hd.*win。歸一化頻率w39。ylabel(39。 )。 解 在MATLAB中可以用remezord 和remez兩個函數(shù)設(shè)計，其結(jié)果如圖2，MATLAB程序如下： fedge=[800 1000]。幅度/dB39。 上圖給出了以上設(shè)計的頻率特性，(a) 為N=30直接截取的頻率特性(b)為凱塞窗設(shè)計的頻率特性。凱塞窗設(shè)計對應(yīng)的MATLAB程序為： wn=kaiser(30,)。) 。 mval=[1 0]。 幅度/db39。)。 [hh1,w1]=freqz(h1,1,256)。bandpass39。 h2=remez(n,fpts,mag,wt)。窗函數(shù)法還包含兩個分支，一種是用公式先手動算出N值和其他對應(yīng)得窗函數(shù)參數(shù)值，再代入窗函數(shù)和fir1實現(xiàn)，一種是用函數(shù)*rord估算出N和相應(yīng)參數(shù)再用fir1實現(xiàn)。 h=2*cos(*pi*(n15)).*sin(*pi*(n15))./(pi*(n15))。 y1=sigsum(0,30,[0:30],(abs(delta1).^2))/31。如果截止頻率在0或Nyquist頻率附近，或者設(shè)定的dev值較大，則得不到正確結(jié)果。 figure(2) subplot(2,1,1) plot(w2/pi,20*log10(abs(hh2))) grid xlabel(39。 [hh2,w2]=freqz(h2,1,256)。數(shù)字頻率w/rad39。幅度/db39。 subplot(2,1,2) plot(w1/pi,angle(hh1)) grid xlabel(39。 [N,fpts,mag,wt]=remezord(fedge,