【正文】 h(n)39。)。xlabel(39。布萊克曼窗波形39。)。stem(n,wd5,39。axis([0 30 0 1])。h(n)39。)。xlabel(39。漢明窗波形39。)。stem(n,wd4,39。axis([0 30 0 1])。h(n)39。)。xlabel(39。漢寧窗時域波形39。)。stem(n,wd3,39。axis([0 30 0 1])。h(n)39。)。xlabel(39。三角窗波形39。)。stem(n,wd2,39。)。ylabel(39。n39。)。title(39。)。stem(n,wd1,39。)。相位/176。)。xlabel(39。凱澤窗()相頻響應39。axis([0 1 2000 200])。subplot(3,2,6)。幅頻響應/dB39。)。xlabel(39。凱澤窗()幅頻響應39。axis([0 1 300 0])。subplot(3,2,4)。h(n)39。)。xlabel(39。凱澤窗()波形39。)。stem(n,wd6,39。)。相位/176。)。xlabel(39。布萊克曼窗相頻響應39。axis([0 1 2000 200])。subplot(3,2,5)。幅頻響應/dB39。)。xlabel(39。布萊克曼窗幅頻響應39。axis([0 1 300 0])。subplot(3,2,3)。h(n)39。)。xlabel(39。布萊克曼窗波形39。)。stem(n,wd5,39。)。相位/176。)。xlabel(39。漢明窗相頻響應39。axis([0 1 2000 200])。subplot(3,2,6)。幅頻響應/dB39。)。xlabel(39。漢明窗幅頻響應39。axis([0 1 300 0])。subplot(3,2,4)。h(n)39。)。xlabel(39。漢明窗波形39。)。stem(n,wd4,39。)。相位/176。)。xlabel(39。漢寧窗相頻響應39。axis([0 1 2000 200])。subplot(3,2,5)。幅頻響應/dB39。)。xlabel(39。漢寧窗幅頻響應39。axis([0 1 300 0])。subplot(3,2,3)。h(n)39。)。xlabel(39。漢寧窗波形39。)。stem(n,wd3,39。)。相位/176。)。xlabel(39。三角窗相頻響應39。axis([0 1 2000 200])。subplot(3,2,6)。幅頻響應/dB39。)。xlabel(39。三角窗幅頻響應39。axis([0 1 300 0])。subplot(3,2,4)。h(n)39。)。xlabel(39。三角窗波形39。)。stem(n,wd2,39。)。相位/176。)。xlabel(39。矩形窗相頻響應39。axis([0 1 2000 200])。subplot(3,2,5)。幅頻響應/dB39。)。xlabel(39。矩形窗幅頻響應39。axis([0 1 300 0])。subplot(3,2,3)。h(n)39。)。xlabel(39。矩形窗波形39。)。stem(n,wd1,39。[H6,w6]=freqz(h6)。%布萊克曼窗wd6=(kaiser(N,))39。h5=hd.*wd5。%漢明窗wd5=(blackman(N))39。h4=hd.*wd4。%漢寧窗wd4=(hamming(N))39。h3=hd.*wd3。%三角窗wd3=(hanning(N))39。h2=hd.*wd2。%矩形窗wd2=(triang(N))39。h1=hd.*wd1。 %理想濾波器的單位取樣響應%%%選擇窗函數(shù)并求其幅頻響應，求FIR數(shù)字濾波器的單位取樣響應和頻率響應wd1=(boxcar(N))39。m=na+eps。 %濾波器相位延遲wc=pi/4。 %清除屏幕顯示 N=23。開始 程序流程圖設置濾波器基本參數(shù)求理想低通濾波器沖激響應設定窗函數(shù)繪制窗函數(shù)特性曲線求FIR低通濾波器沖激響應和頻率響應設定采樣頻率、輸入信號并繪制輸入信號用上述濾波器對輸入信號進行濾波繪制各濾波器特性曲線及輸出信號結束圖5 程序設計流程圖 濾波器仿真結果（1）N=23，ωC=π/4時各種窗函數(shù)的時域波形及幅頻響應圖6 矩形窗和三角窗波形及幅頻響應圖7 漢寧窗和漢明窗波形及幅頻響應圖8 布萊克曼和凱澤（）窗波形及幅頻響應（2）輸入待濾波信號圖9 待濾波信號波形（3）N=23時，用窗函數(shù)設計濾波器的沖激響應、幅頻響應及輸出波形圖10 矩形窗設計濾波器的特性曲線及輸出波形圖11 三角窗設計濾波器特曲線及輸出波形圖12 漢寧窗設計濾波器特曲線及輸出波形圖13 漢寧窗設計濾波器特曲線及輸出波形圖14 布萊克曼窗設計濾波器特曲線及輸出波形圖15 凱澤窗（β=）設計濾波器特曲線及輸出波形（4）改變N值大小，矩形窗設計濾波器的幅頻響應的變化情況圖16 N=13和N=23時矩形窗濾波器頻率響應圖17 N=35和N=45時矩形窗濾波器頻率響應4 源程序代碼clear。注：漢寧窗通過調用函數(shù)w = hanning(n)來生成，不返回零點的窗函數(shù)的首尾兩個元素。 窗函數(shù)法編程實現(xiàn)濾波器功能 窗函數(shù)的實現(xiàn)表3 不同窗函數(shù)的調用函數(shù)窗函數(shù)調用函數(shù)矩形窗w=boxcar（n）三角窗w= triang(n)漢寧窗w = hanning(n) 或 w = hanning(n,‘periodic’)漢明窗w = hamming(n) 或 w = hamming(n,sflag)布萊克曼窗w = blackman (n) 或w = blackman (n,sflag)凱澤窗w = kaiser(n,beta)其中n為窗的長度，beta為凱澤窗參數(shù)β值，輸出參數(shù)w是由窗函數(shù)的值組成的N階向量。圖2 FIR數(shù)字濾波器系統(tǒng)圖3 濾波前信號 圖4 濾波后信號通過仿真結果可以看出該濾波器可以有效抑制干擾的高頻信號，而是有效低頻信號得以保留。濾波前后的時域波形分別經(jīng)各自的示波器模塊顯示。輸入的兩個信號經(jīng)過加法器進行疊加，其中sine wave輸出為有用低頻信號，sine wave1輸出為高頻信號。建立仿真系統(tǒng)時，注意選用數(shù)字模塊，并進行相應參數(shù)的設計。3 濾波器的實現(xiàn) Simulink模塊實現(xiàn)濾波器功能Simulink是MATLAB里面的重要工具箱之一．其主要功能是實現(xiàn)系統(tǒng)建模、仿真與分析，從而可以在實際系統(tǒng)制作出來之前，預先對系統(tǒng)進行仿真與分析，并可以對系統(tǒng)做適當?shù)膶崟r修正或者調試及系統(tǒng)的參數(shù)，以提高系統(tǒng)的性能，減少系統(tǒng)設計的修改的時間，高效率地開發(fā)系統(tǒng)。表2歸納了以上6種窗的主要性能，可供設計FIR濾波器時參考。當β=，相當于漢明窗，當凱澤窗旁瓣頻譜收斂得更快，能量更集中在主瓣中。90177。70177。50177。30177。凱澤窗在不同β值下的性能歸納在表1中。改變β值可對主瓣寬度與旁瓣衰減進行選擇。凱澤窗可以在主瓣寬度和旁瓣高度之間自由選擇它們的比重，適應性較強。為了更進一步抑制旁瓣，可再加上余弦的二次諧波分量，得到布萊克曼窗：其頻譜的幅度函數(shù)為： 此時主瓣寬度為矩形窗譜主瓣寬的的3倍，即12π/N。對余弦窗加以改進，可以使旁瓣更小，窗函數(shù)為：其頻率響應的幅度函數(shù)為：漢明窗函數(shù)的最大旁瓣