1、 -9 布萊克曼窗(N=33)四 . 實(shí)驗(yàn)報(bào)告要求簡(jiǎn)述實(shí)驗(yàn)原理及目的。按照實(shí)驗(yàn)步驟以及要求，比較各種情況下的濾波性能，說(shuō)明窗口長(zhǎng)度 N 和窗函數(shù)類(lèi)型對(duì)濾波特性的影響。總結(jié)用窗函數(shù)法設(shè)計(jì)FIR 濾波器的主要特點(diǎn)。簡(jiǎn)要回答思考題。五 . 思考題如果給定通帶截止頻率和阻帶截止頻率以及阻帶最小衰減，如何用窗函數(shù)法設(shè)計(jì)線(xiàn)性相位低通濾波器，寫(xiě)出設(shè)計(jì)步驟。如果要求用窗函數(shù)法設(shè)計(jì)帶通濾波器，而且給定上、下邊帶截止頻率為 1 和 2，試求理想帶通的單位脈沖響應(yīng)hd n%實(shí)驗(yàn)中的子函數(shù)：產(chǎn)生理想低通濾波器單位脈沖響應(yīng)hd(n)function hd=ideal(w,N);alpha=(N-1)/2;n=0:(N-

2、1);m=n-alpha+eps;hd=sin(w*m)./(pi*m);相關(guān)練習(xí)：例 1 用凱塞窗設(shè)計(jì)一FIR低通濾波器，低通邊界頻率， 阻帶邊界頻率，阻帶衰減不小于50dB。解 首先由過(guò)渡帶寬和阻帶衰減來(lái)決定凱塞窗的N 和圖 4.1 給出了以上設(shè)計(jì)的頻率特性，(a) 為 N=30直接截取的頻率特性 (b) 為凱塞窗設(shè)計(jì)的頻率特性。凱塞窗設(shè)計(jì)對(duì)應(yīng)的MATLAB程序?yàn)椋簑n=kaiser(30,4.55);nn=0:1:29;alfa=(30-1)/2;hd=sin(0.4*pi*(nn-alfa)./(pi*(nn-alfa );h=hd.*wn;h1,w1=freqz(h,1);plot(

3、w1/pi,20*log10(abs(h1);axis(0,1,-80,10);grid;xlabel( 歸一化頻率/ )ylabel( 幅度 /dB)例 2 利用雷米茲交替算法，設(shè)計(jì)一個(gè)線(xiàn)性相位低通FIR 數(shù)字濾波器，其指標(biāo)為：通帶邊界頻率f c=800Hz，阻帶邊界fr=1000Hz，通帶波動(dòng)阻帶最小衰減At=40dB，采樣頻率fs=4000Hz。在 MATLAB中可以用remezord 和 remez兩個(gè)函數(shù)設(shè)計(jì)，其結(jié)果如圖4.2，MATLAB程序如下：fedge=800 1000;mval=1 0;dev=0.0559 0.01;fs=4000;N,fpts,mag,wt=remezord(fedge,mval,dev,fs);b=remez(N,fpts,mag,wt);h,w=freqz(b,1,256);plot(w*2000/pi,20*log10(abs(h);grid;xlabel( 頻率 /Hz)ylabel( 幅度 /dB)函數(shù) remezord 中的數(shù)組f