【文章內(nèi)容簡介】 設(shè)計濾波器特曲線及輸出波形 武漢理工大學(xué)《專業(yè)課程設(shè)計 3（數(shù)字信號處理）》課程設(shè)計說明書 14 （ 4）改變 N 值大小，矩形窗 設(shè)計濾波器的幅頻響應(yīng)的變化情況 圖 16 N=13 和 N=23 時矩形窗濾波器頻率響應(yīng) 圖 17 N=35 和 N=45 時矩形窗濾波器頻率響應(yīng) 武漢理工大學(xué)《專業(yè)課程設(shè)計 3（數(shù)字信號處理）》課程設(shè)計說明書 15 4 源程序代碼 clear。 %清除工作空間 clc。 %清除屏幕顯示 N=23。 %濾波去階數(shù) a=(N1)/2。 %濾波器相位延遲 wc=pi/4。 %截止頻率 n=0:N1。 m=na+eps。 %避免被零除， eps 為 MATLAB 的系統(tǒng)精度 hd=sin(wc*m)./(pi*m)。 %理想濾波器的單位取樣響應(yīng) %%%選擇 窗函數(shù) 并求其幅頻響應(yīng) ，求 FIR 數(shù)字濾波器的單位取樣響應(yīng) 和 頻率響應(yīng) wd1=(boxcar(N))39。 [H10,w10]=freqz(wd1)。h1=hd.*wd1。[H1,w1]=freqz(h1)。%矩形窗 wd2=(triang(N))39。 [H20,w20]=freqz(wd2)。h2=hd.*wd2。[H2,w2]=freqz(h2)。%三角窗 wd3=(hanning(N))39。 [H30,w30]=freqz(wd3)。h3=hd.*wd3。[H3,w3]=freqz(h3)。%漢 寧 窗 wd4=(hamming(N))39。 [H40,w40]=freqz(wd4)。h4=hd.*wd4。[H4,w4]=freqz(h4)。%漢明 窗 wd5=(blackman(N))39。 [H50,w50]=freqz(wd5)。h5=hd.*wd5。[H5,w5]=freqz(h5)。%布萊克曼窗 wd6=(kaiser(N,))39。 [H60,w60]=freqz(wd6)h6=hd.*wd6。[H6,w6]=freqz(h6)。%凱澤窗 %%%繪制窗函數(shù)的沖激響應(yīng)、幅頻響應(yīng) subplot(3,2,1)。stem(n,wd1,39。.39。)。 %矩形窗沖激響應(yīng) title(39。矩形窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。 subplot(3,2,3)。plot(w10/pi,20*log10(abs(H10)/max(abs(H10))))。axis([0 1 300 0])。 %幅頻響應(yīng) title(39。矩形窗幅頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。幅頻響應(yīng) /dB39。)。 subplot(3,2,5)。plot(w10/pi,180/pi*unwrap(angle(H10)))。axis([0 1 2020 200])。 %相頻響應(yīng) title(39。矩形窗相頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。相位 /176。39。)。 subplot(3,2,2)。stem(n,wd2,39。.39。)。 title(39。三角窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。 subplot(3,2,4)。plot(w20/pi,20*log10(abs(H20)/max(abs(H20))))。axis([0 1 300 0])。 title(39。三角窗幅頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。幅頻響應(yīng) /dB39。)。 subplot(3,2,6)。plot(w20/pi,180/pi*unwrap(angle(H20)))。axis([0 1 2020 200])。 武漢理工大學(xué)《專業(yè)課程設(shè)計 3（數(shù)字信號處理）》課程設(shè)計說明書 16 title(39。三角窗相頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。相位 /176。39。)。 subplot(3,2,1)。stem(n,wd3,39。.39。)。 title(39。漢寧窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。 subplot(3,2,3)。plot(w30/pi,20*log10(abs(H30)/max(abs(H30))))。axis([0 1 300 0])。 title(39。漢寧窗幅頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。幅頻響應(yīng) /dB39。)。 subplot(3,2,5)。plot(w30/pi,180/pi*unwrap(angle(H30)))。axis([0 1 2020 200])。 title(39。漢寧窗相頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。相位 /176。39。)。 subplot(3,2,2)。stem(n,wd4,39。.39。)。 title(39。漢明窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。 subplot(3,2,4)。plot(w40/pi,20*log10(abs(H40)/max(abs(H40))))。axis([0 1 300 0])。 title(39。漢明窗幅頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。幅頻響應(yīng) /dB39。)。 subplot(3,2,6)。plot(w40/pi,180/pi*unwrap(angle(H40)))。axis([0 1 2020 200])。 title(39。漢明窗相頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。相位 /176。39。)。 subplot(3,2,1)。stem(n,wd5,39。.39。)。 title(39。布萊克曼窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。 subplot(3,2,3)。plot(w50/pi,20*log10(abs(H50)/max(abs(H50))))。axis([0 1 300 0])。 title(39。布萊克曼窗幅頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。幅頻響應(yīng) /dB39。)。 subplot(3,2,5)。plot(w50/pi,180/pi*unwrap(angle(H50)))。axis([0 1 2020 200])。 title(39。布萊克曼窗相頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。相位 /176。39。)。 subplot(3,2,2)。stem(n,wd6,39。.39。)。 title(39。凱澤窗 ()波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。 subplot(3,2,4)。plot(w60/pi,20*log10(abs(H60)/max(abs(H60))))。axis([0 1 300 0])。 title(39。凱澤窗 ()幅頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。幅頻響應(yīng) /dB39。)。 subplot(3,2,6)。plot(w60/pi,180/pi*unwrap(angle(H60)))。axis([0 1 2020 200])。 title(39。凱澤窗 ()相頻響應(yīng) 39。)。xlabel(39。歸一化頻率 39。)。ylabel(39。相位 /176。39。)。 %%%繪制各種窗函數(shù)的幅頻響應(yīng) subplot(2,2,1)。stem(n,wd1,39。.39。)。axis([0 30 0 1])。 title(39。矩形窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。 subplot(2,2,2)。stem(n,wd2,39。.39。)。 武漢理工大學(xué)《專業(yè)課程設(shè)計 3（數(shù)字信號處理）》課程設(shè)計說明書 17 title(39。三角窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。axis([0 30 0 1])。 subplot(2,2,3)。stem(n,wd3,39。.39。)。 title(39。漢寧 窗時域波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。axis([0 30 0 1])。 subplot(2,2,4)。stem(n,wd4,39。.39。)。 title(39。漢明 窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。axis([0 30 0 1])。 subplot(2,2,1)。stem(n,wd5,39。.39。)。 title(39。布萊克曼窗波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。axis([0 30 0 1])。 subplot(2,2,2)。stem(n,wd6,39。.39。)。 title(39。凱澤窗 ()波形 39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。axis([0 30 0 1])。 %對輸入信號進(jìn)行采樣濾波 NN=1000。Fs=1000。%數(shù)據(jù)總數(shù)和采樣頻率 nn=[0:NN1]。t=nn/Fs。 f1=50。f2= x=sin(2*pi*f1*t)+cos(2*pi*f2*t)。%待濾波的信號 n1=81:241。t1=t(n1)。%選擇采樣點間隔 x1=x(n1)。%與采樣點對應(yīng)的輸入信號 plot(t1,x1)。 title(39。帶濾波信號 39。)。xlabel(39。時間 /t39。)。ylabel(39。幅度 39。)。 %%%把輸入信號經(jīng)濾波器處理 y1=fftfilt(h1,x,256)。 y2=fftfilt(h