【正文】 grd = grpdelay(b,a,w)。附錄2：function hd = ideal_lp(wc,M)。)。)。)。axis([0,10000,0,1000])。Y_fil=Y_fil(1:N/2)。grid on。ylabel(39。xlabel(39。dw=df/fs*2*pi。ylabel(39。xlabel(39。title(39。 % 計算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。% 程序作者：余霏霖/% 最后修改日期：201338[x,fs,bits]=wavread(39。對各時域、頻譜圖的范圍沒有好的預(yù)計，導(dǎo)致出圖時波形效果不理想。 % 調(diào)用自編函數(shù)計算理想帶阻濾波器的脈沖響應(yīng)h_bs=w_ham39。fsu=2xlabel(39。title(39。 %加上標(biāo)題和橫坐標(biāo)名稱grid on。 Y=Y(1:N/2)。調(diào)用程序如下：[x,fs,bits]=wavread(39。（2）根據(jù)待求濾波器的理想頻率響應(yīng)求出理想單位脈沖響應(yīng)hd(n)。20世紀(jì)60年代以來，隨著計算機和信息技術(shù)的飛速發(fā)展，數(shù)字信號處理技術(shù)應(yīng)運而生并得到迅速的發(fā)展。（3）處理結(jié)果和分析結(jié)論應(yīng)該一致，而且應(yīng)符合理論。通過本次課設(shè)，我們能夠?qū)W會如何綜合運用課堂上學(xué)會的理論知識，增強自己的動能力與聯(lián)系實際的能力，為以后的工作奠定基礎(chǔ)。MATLAB在信號與系統(tǒng)中的另一主要應(yīng)用是數(shù)值計算與仿真分析，主要包括函數(shù)波形繪制、函數(shù)運算、沖激響應(yīng)與階躍響應(yīng)仿真分析、信號的時域分析、信號的頻譜分析、系統(tǒng)的S域分析、零極點圖繪制等內(nèi)容。在大多數(shù)實際應(yīng)用中，數(shù)字濾波器常常被用來實現(xiàn)選頻操作，所以指標(biāo)的形式一般為在頻域中以分貝值給出的相對幅度響應(yīng)和相位響應(yīng)。CE軟件操作界面如圖32所示。調(diào)用程序如下：X=abs(fft(x))。幅度39。grid on。title(39。再調(diào)用freqz(b,a,512,fs)即可得到該濾波器的頻率響應(yīng)。 % 計算漢寧窗設(shè)計該濾波器時需要的階數(shù)n=0:M1。聲音中刺耳的聲音沒有了，幾乎恢復(fù)成原始的聲音，但較原始的聲音更平滑一些。K % 對原始信號和加噪信號進行fft變換X=X(1:N/2)。 %加上標(biāo)題和橫坐標(biāo)名稱grid on。title(39。xlabel(39。df=min((fsdfpd),(fpufsu))。title(39。xlabel(39。相位（pha）39。grid on。grid on。axis([0,2,])。)。)。)。db = 20*log10((mag+eps)/max(mag))。% grd = [0 grd(1:1:500)。m = n alpha + eps。plot(f,Y_fil)。)。)。)。 axis([0,2,])。axis([0,800,1])。grid on。ylabel(39。M=ceil(*pi/dw)+1。axis([0,3000,0,3000])。ylabel(39。xlabel(39。 % 布局為2*2的四個小圖title(39。 % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個參數(shù)是每個樣本的值，fs是生成該波形文件時的采樣率，bits是波形文件每樣本的編碼位數(shù)。5 結(jié)束語在此次課程設(shè)計中，我的任務(wù)是利用漢寧窗函數(shù)設(shè)計FIR濾波器，對吹管音樂信號進行濾波去噪處理。 % 調(diào)用自編函數(shù)計算濾波器的頻率特性通過繪圖工具可得出濾波器的波形圖，如圖34所示。As=37。ylabel(39。xlabel(39。title(39。 % 計算頻譜頻率范圍figure(1)。 % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個參數(shù)是每個樣本的值，fs是生成該波形文件時的采樣率，bits是波形文件每樣本的編碼位數(shù)。 常見的窗函數(shù)性能表如下表21所示。[5] FIR濾波器有限長單位脈沖響應(yīng)數(shù)字濾波器（Finite Impulse Response Digital Filter,縮寫FIRDF）簡稱FIR濾波器，是數(shù)字信號處理系統(tǒng)中最基本的原件，其最大優(yōu)點是可以實現(xiàn)線性相位濾波，可以在保證任意幅頻特性的同時具有嚴(yán)格的線性相頻特性，滿足了在數(shù)字通信和圖像傳輸與處理等應(yīng)用場合對線性相位的要求。MATLAB 的名稱源自 Matrix Laboratory ，它是美國MathWorks公司生產(chǎn)的一個為科學(xué)和工程計算專門設(shè)計的交互式大型軟件，是一個可以完成各種精確計算和數(shù)據(jù)處理的、可視化的、強大的計算工具。關(guān)鍵詞 濾波去噪；FIR濾波器；漢寧窗；MATLAB1 引 言本課程設(shè)計主要針對一段吹管音樂信號，在進行加噪后，利用窗函數(shù)設(shè)計法選擇漢寧窗設(shè)計的FIR濾波器，對加噪后的吹管音樂信號進行濾波去噪處理，并對前后時域波形和頻域波形進行對比分析的程序設(shè)計。[2]MATLAB在信號與系統(tǒng)中的應(yīng)用主要包括符號運算和數(shù)值計算仿真分析。為了區(qū)別于幅頻響應(yīng)函數(shù)和相頻響應(yīng)函數(shù)，稱為幅頻特性函數(shù)，稱為相頻特性函數(shù)。程序的設(shè)計流程圖如下圖31所示。y=x+sin(fn*2*pi*t)。時間(t)39。幅度譜3